Genome assembly and annotation ofArabidopsis halleri, a model for heavy metal hyperaccumulation and evolutionary ecology

Convergent evolution in Arabidopsis halleri and Arabidopsis arenosa on calamine metalliferous soils

It is a plausible hypothesis that parallel adaptation events to the same environmental challenge should result in genetic changes of similar or identical effects, depending on the underlying fitness landscapes. However, systematic testing of this is scarce. Here we examine this hypothesis in two closely related plant species, Arabidopsis halleri and Arabidopsis arenosa, which co-occur at two calamine metalliferous (M) sites harbouring toxic levels of the heavy metals zinc and cadmium. We conduct individual genome resequencing alongside soil elemental analysis for 64 plants from eight populations on M and non-metalliferous (NM) soils, and identify genomic footprints of selection and local adaptation. Selective sweep and environmental association analyses indicate a modest degree of gene as well as functional network convergence, whereby the proximal molecular factors mediating this convergence mostly differ between site pairs and species. Notably, we observe repeated selection on identical single nucleotide polymorphisms in several A. halleri genes at two independently colonized M sites. Our data suggest that species-specific metal handling and other biological features could explain a low degree of convergence between species. The parallel establishment of plant populations on calamine M soils involves convergent evolution, which will probably be more pervasive across sites purposely chosen for maximal similarity in soil composition. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Convergent evolution in Arabidopsis halleri and Arabidopsis arenosa on calamine metalliferous soils

It is a plausible hypothesis that parallel adaptation events to the same environmental challenge should result in genetic changes of similar or identical effects, depending on the underlying fitness landscapes. However, systematic testing of this is scarce. Here we examine this hypothesis in two closely related plant species, Arabidopsis halleri and Arabidopsis arenosa, which co-occur at two calamine metalliferous (M) sites harbouring toxic levels of the heavy metals zinc and cadmium. We conduct individual genome resequencing alongside soil elemental analysis for 64 plants from eight populations on M and non-metalliferous (NM) soils, and identify genomic footprints of selection and local adaptation. Selective sweep and environmental association analyses indicate a modest degree of gene as well as functional network convergence, whereby the proximal molecular factors mediating this convergence mostly differ between site pairs and species. Notably, we observe repeated selection on identical single nucleotide polymorphisms in several A. halleri genes at two independently colonized M sites. Our data suggest that species-specific metal handling and other biological features could explain a low degree of convergence between species. The parallel establishment of plant populations on calamine M soils involves convergent evolution, which will probably be more pervasive across sites purposely chosen for maximal similarity in soil composition. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Differential retention of transposable element-derived sequences in outcrossing Arabidopsis genomes

Background

Transposable elements (TEs) are genomic parasites with major impacts on host genome architecture and host adaptation. A proper evaluation of their evolutionary significance has been hampered by the paucity of short scale phylogenetic comparisons between closely related species. Here, we characterized the dynamics of TE accumulation at the micro-evolutionary scale by comparing two closely related plant species, Arabidopsis lyrata and A. halleri.

Results

Joint genome annotation in these two outcrossing species confirmed that both contain two distinct populations of TEs with either 'recent' or 'old' insertion histories. Identification of rare segregating insertions suggests that diverse TE families contribute to the ongoing dynamics of TE accumulation in the two species. Orthologous TE fragments (i.e. those that have been maintained in both species), tend to be located closer to genes than those that are retained in one species only. Compared to non-orthologous TE insertions, those that are orthologous tend to produce fewer short interfering RNAs, are less heavily methylated when found within or adjacent to genes and these tend to have lower expression levels. These findings suggest that long-term retention of TE insertions reflects their frequent acquisition of adaptive roles and/or the deleterious effects of removing nearly neutral TE insertions when they are close to genes.

Conclusion

Our results indicate a rapid evolutionary dynamics of the TE landscape in these two outcrossing species, with an important input of a diverse set of new insertions with variable propensity to resist deletion.

Seasonal Stability and Dynamics of DNA Methylation in Plants in a Natural Environment

DNA methylation has been considered a stable epigenetic mark but may respond to fluctuating environments. However, it is unclear how they behave in natural environments. Here, we analyzed seasonal patterns of genome-wide DNA methylation in a single clone from a natural population of the perennial Arabidopsis halleri. The genome-wide pattern of DNA methylation was primarily stable, and most of the repetitive regions were methylated across the year. Although the proportion was small, we detected seasonally methylated cytosines (SeMCs) in the genome. SeMCs in the CHH context were detected predominantly at repetitive sequences in intergenic regions. In contrast, gene-body CG methylation (gbM) itself was generally stable across seasons, but the levels of gbM were positively associated with seasonal stability of RNA expression of the genes. These results suggest the existence of two distinct aspects of DNA methylation in natural environments: sources of epigenetic variation and epigenetic marks for stable gene expression.

Family quarrels in seeds and rapid adaptive evolution in Arabidopsis

Evolutionary conflict can drive rapid adaptive evolution, sometimes called an arms race, because each party needs to respond continually to the adaptations of the other. Evidence for such arms races can sometimes be seen in morphology, in behavior, or in the genes underlying sexual interactions of host-pathogen interactions, but is rarely predicted a priori. Kin selection theory predicts that conflicts of interest should usually be reduced but not eliminated among genetic relatives, but there is little evidence as to whether conflict within families can drive rapid adaptation. Here we test multiple predictions about how conflict over the amount of resources an offspring receives from its parent would drive rapid molecular evolution in seed tissues of the flowering plant Arabidopsis As predicted, there is more adaptive evolution in genes expressed in Arabidopsis seeds than in other specialized organs, more in endosperms and maternal tissues than in embryos, and more in the specific subtissues involved in nutrient transfer. In the absence of credible alternative hypotheses, these results suggest that kin selection and conflict are important in plants, that the conflict includes not just the mother and offspring but also the triploid endosperm, and that, despite the conflict-reducing role of kinship, family members can engage in slow but steady tortoise-like arms races.

Genetic architecture of a plant adaptive trait: QTL mapping of intraspecific variation for tolerance to metal pollution in Arabidopsis halleri

Anthropogenic activities are among the main drivers of global change and result in drastic habitat modifications, which represent strong evolutionary challenges for biological species that can either migrate, adapt, or disappear. In this context, understanding the genetics of adaptive traits is a prerequisite to enable long-term maintenance of populations under strong environmental constraints. To examine these processes, a QTL approach was developed here using the pseudometallophyte Arabidopsis halleri, which displays among-population adaptive divergence for tolerance to metallic pollution in soils. An F2 progeny was obtained by crossing individuals from metallicolous and non-metallicolous populations from Italian Alps, where intense metallurgic activities have created strong landscape heterogeneity. Then, we combined genome de novo assembly and genome resequencing of parental genotypes to obtain single-nucleotide polymorphism markers and achieve high-throughput genotyping of the progeny. QTL analysis was performed using growth parameters and photosynthetic yield to assess zinc tolerance levels. One major QTL was identified for photosynthetic yield. It explained about 27% of the phenotypic variance. Functional annotation of the QTL and gene expression analyses highlighted putative candidate genes. Our study represents a successful approach combining evolutionary genetics and advanced molecular tools, helping to better understand how a species can face new selective pressures of anthropogenic origin.

Heavy metal rich stone-processing wastewater inhibits the growth and development of plants

Large amounts of wastewater are generated from stone processing, which are toxic and cause serious environmental and health risks. To quantify the content of stone processing wastewater and estimate its effects on plant growth, we collected water samples from sewage outfall of four stone processing factories and nearby water bodies. The concentration of potential toxic metals were much higher in the wastewater than background controls. Wastewater inhibited plant primary root elongation, lateral root formation, and growth of aerial part. Seedlings treated with the effluents were unhealthy with deep purple leaves and usually died before flowering. Chlorophyll a/b contents and chloroplast number were reduced in those abnormal mesophyll cells. Transcriptional levels were decreased for chloroplast formation genes, but increased for those participated in chloroplast degradation and catabolism. Six out of nine tested senescence-associated genes were up-regulated. Furthermore, our results show that endogenous toxic metal levels indeed increased after wastewater treatment. Altogether, these results indicated that the potential toxic metals rich wastewater had significant inhibition on plant growth and led to senescence-associated program cell death, which could be helpful for the government and enterprises to understand the environmental risks and formulate reasonable wastewater emission standards for the stone processing industry.

Transmembrane transport and stress response genes play an important role in adaptation of Arabidopsis halleri to metalliferous soils

When plants adapt to local environments, strong signatures of selection are expected in the genome, particularly in high-stress environments such as trace metal element enriched (metalliferous) soils. Using Arabidopsis halleri, a model species for metal homeostasis and adaptation to extreme environments, we identifid genes, gene variants, and pathways that are associated with soil properties and may thus contribute to adaptation to high concentrations of trace metal elements. We analysed whole-genome Pool-seq data from two metallicolous (from metalliferous soils) and two non-metallicolous populations (in total 119 individuals) and associated allele frequencies of the identified single-nucleotide polymorphisms (SNPs) with soil variables measured on site. Additionally, we accounted for polygenic adaptation by searching for gene pathways showing enrichment of signatures of selection. Out of >2.5 million SNPs, we identified 57 SNPs in 19 genes that were significantly associated with soil variables and are members of three enriched pathways. At least three of these candidate genes and pathways are involved in transmembrane transport and/or associated with responses to various stresses such as oxidative stress. We conclude that both allocation and detoxification processes play a crucial role in A. halleri for coping with these unfavourable conditions.

Higher Rates of Protein Evolution in the Self-Fertilizing Plant Arabidopsis thaliana than in the Out-Crossers Arabidopsis lyrata and Arabidopsis halleri

The common transition from out-crossing to self-fertilization in plants decreases effective population size. This is expected to result in a reduced efficacy of natural selection and in increased rates of protein evolution in selfing plants compared with their outcrossing congeners. Prior analyses, based on a very limited number of genes, detected no differences between the rates of protein evolution in the selfing Arabidopsis thaliana compared with the out-crosser Arabidopsis lyrata. Here, we reevaluate this trend using the complete genomes of A. thaliana, A. lyrata, Arabidopsis halleri, and the outgroups Capsella rubella and Thellungiella parvula. Our analyses indicate slightly but measurably higher nonsynonymous divergences (dN), synonymous divergences (dS) and dN/dS ratios in A. thaliana compared with the other Arabidopsis species, indicating that purifying selection is indeed less efficacious in A. thaliana.

CAX1 suppresses Cd-induced generation of reactive oxygen species in Arabidopsis halleri

The molecular analysis of metal hyperaccumulation in species such as Arabidopsis halleri offers the chance to gain insights into metal homeostasis and into the evolution of adaptation to extreme habitats. A prerequisite of metal hyperaccumulation is metal hypertolerance. Genetic analysis of a backcross population derived from Arabidopsis lyrata × A. halleri crosses revealed three quantitative trait loci for Cd hypertolerance. A candidate gene for Cdtol2 is AhCAX1, encoding a vacuolar Ca²⁺ /H⁺ antiporter. We developed a method for the transformation of vegetatively propagated A. halleri plants and generated AhCAX1-silenced lines. Upon Cd²⁺ exposure, several-fold higher accumulation of reactive oxygen species (ROS) was detectable in roots of AhCAX1-silenced plants. In accordance with the dependence of Cdtol2 on external Ca²⁺ concentration, this phenotype was exclusively observed in low Ca²⁺ conditions. The effects of external Ca²⁺ on Cd accumulation cannot explain the phenotype as they were not influenced by the genotype. Our data strongly support the hypothesis that higher expression of CAX1 in A. halleri relative to other Arabidopsis species represents a Cd hypertolerance factor. We propose a function of AhCAX1 in preventing a positive feedback loop of Cd-elicited ROS production triggering further Ca²⁺ -dependent ROS accumulation.

Patterns of polymorphism and selection in the subgenomes of the allopolyploid Arabidopsis kamchatica

Genome duplication is widespread in wild and crop plants. However, little is known about genome-wide selection in polyploids due to the complexity of duplicated genomes. In polyploids, the patterns of purifying selection and adaptive substitutions may be affected by masking owing to duplicated genes or homeologs as well as effective population size. Here, we resequence 25 accessions of the allotetraploid Arabidopsis kamchatica, which is derived from the diploid species A. halleri and A. lyrata. We observe a reduction in purifying selection compared with the parental species. Interestingly, proportions of adaptive non-synonymous substitutions are significantly positive in contrast to most plant species. A recurrent pattern observed in both frequency and divergence–diversity neutrality tests is that the genome-wide distributions of both subgenomes are similar, but the correlation between homeologous pairs is low. This may increase the opportunity of different evolutionary trajectories such as in the HMA4 gene involved in heavy metal hyperaccumulation.

Despite the prevalence of genome duplication in plants, little is known about the evolutionary patterns of entire subgenomes. Here the authors resequence allopolyploid Arabidopsis kamchatica genome to estimate diversity, linkage disequilibrium and strengths of both positive and purifying selection.

Evolutionary dynamics of quantitative variation in an adaptive trait at the regional scale: The case of zinc hyperaccumulation in Arabidopsis halleri

Metal hyperaccumulation in plants is an ecological trait whose biological significance remains debated, in particular because the selective pressures that govern its evolutionary dynamics are complex. One of the possible causes of quantitative variation in hyperaccumulation may be local adaptation to metalliferous soils. Here, we explored the population genetic structure of Arabidopsis halleri at fourteen metalliferous and nonmetalliferous sampling sites in southern Poland. The results were integrated with a quantitative assessment of variation in zinc hyperaccumulation to trace local adaptation. We identified a clear hierarchical structure with two distinct genetic groups at the upper level of clustering. Interestingly, these groups corresponded to different geographic subregions, rather than to ecological types (i.e., metallicolous vs. nonmetallicolous). Also, approximate Bayesian computation analyses suggested that the current distribution of A. halleri in southern Poland could be relictual as a result of habitat fragmentation caused by climatic shifts during the Holocene, rather than due to recent colonization of industrially polluted sites. In addition, we find evidence that some nonmetallicolous lowland populations may have actually derived from metallicolous populations. Meanwhile, the distribution of quantitative variation in zinc hyperaccumulation did separate metallicolous and nonmetallicolous accessions, indicating more recent adaptive evolution and diversifying selection between metalliferous and nonmetalliferous habitats. This suggests that zinc hyperaccumulation evolves both ways-towards higher levels at nonmetalliferous sites and lower levels at metalliferous sites. Our results open a new perspective on possible evolutionary relationships between A. halleri edaphic types that may inspire future genetic studies of quantitative variation in metal hyperaccumulation.

Genome Editing Weds CRISPR: What Is in It for Phytoremediation?

The arrival of sequence-specific endonucleases that allow genome editing has shaken the pillars of basic and applied plant biology. Clustered regularly interspaced palindromic repeats (CRISPR) is a revolutionary genome-engineering tool that enables the enhancement of targeted traits in plants. Numerous plants, including energy crops, known for their potential to tolerate, immobilize, and stabilize inorganic and organic pollutants, have already been edited using different CRISPR systems. Moreover, a large array of genes responsible for increased metal tolerance, metal uptake and hyperaccumulation have already been identified. Thus, the CRISPR-mediated genome reprogramming of plants, including its use in gene expression regulation through transcriptional repression or activation (CRISPRi and CRISPRa), could be of paramount importance for phytoremediation. The simplicity, inexpensiveness, and capabilities of this gene editing technique could soon be used to enhance plants and bacteria involved in phytotechnologies, such as phystabilization, phytoextraction, phytomining, phytovolatilization, and bio-energy generation. In this brief viewpoint piece, we posit some of the potential benefits of CRISPR for phytoremediation.

Agrobacterium-mediated floral dip transformation of the model polyploid species Arabidopsis kamchatica

Polyploidization has played an important role in the speciation and diversification of plant species. However, genetic analyses of polyploids are challenging because the vast majority of the model species are diploids. The allotetraploid Arabidopsis kamchatica, which originated through the hybridization of the diploid Arabidopsis halleri and Arabidopsis lyrata, is an emerging model system for studying various aspects of polyploidy. However, a transgenic method that allows the insertion of a gene of interest into A. kamchatica is still lacking. In this study, we investigated the early development of pistils in A. kamchatica and confirmed the formation of open pistils in young flower buds (stages 8-9), which is important for allowing Agrobacterium to access female reproductive tissues. We established a simple Agrobacterium-mediated floral dip transformation method to transform a gene of interest into A. kamchatica by dipping A. kamchatica inflorescences bearing many young flower buds into a 5% sucrose solution containing 0.05% Silwet L-77 and Agrobacterium harboring the gene of interest. We showed that a screenable marker comprising fluorescence-accumulating seed technology with green fluorescent protein was useful for screening the transgenic seeds of two accessions of A. kamchatica subsp. kamchatica and an accession of A. kamchatica subsp. kawasakiana.

Functional divergence of duplicate genes several million years after gene duplication in Arabidopsis

Lineage-specific duplicated genes likely contribute to the phenotypic divergence in closely related species. However, neither the frequency of duplication events nor the degree of selection pressures immediately after gene duplication is clear in the speciation process. Here, using Illumina DNA-sequencing reads from Arabidopsis halleri, which has multiple closely related species with high-quality genome assemblies (A. thaliana and A. lyrata), we succeeded in generating orthologous gene groups in Brassicaceae. The duplication frequency of retained genes in the Arabidopsis lineage was ∼10 times higher than the duplication frequency inferred by comparative genomics of Arabidopsis, poplar, rice and moss (Physcomitrella patens). The difference of duplication frequencies can be explained by a rapid decay of anciently duplicated genes. To examine the degree of selection pressure on genes duplicated in either the A. halleri-lyrata or the A. halleri lineage, we examined positive and purifying selection in the A. halleri-lyrata and A. halleri lineages throughout the ratios of nonsynonymous to synonymous substitution rates (KA/KS). Duplicate genes tended to have a higher proportion of positive selection compared with non-duplicated genes. Interestingly, we found that functional divergence of duplicated genes was accelerated several million years after gene duplication compared with immediately after gene duplication.

Reference-guided de novo assembly approach improves genome reconstruction for related species

Background

The development of next-generation sequencing has made it possible to sequence whole genomes at a relatively low cost. However, de novo genome assemblies remain challenging due to short read length, missing data, repetitive regions, polymorphisms and sequencing errors. As more and more genomes are sequenced, reference-guided assembly approaches can be used to assist the assembly process. However, previous methods mostly focused on the assembly of other genotypes within the same species. We adapted and extended a reference-guided de novo assembly approach, which enables the usage of a related reference sequence to guide the genome assembly. In order to compare and evaluate de novo and our reference-guided de novo assembly approaches, we used a simulated data set of a repetitive and heterozygotic plant genome.

Results

The extended reference-guided de novo assembly approach almost always outperforms the corresponding de novo assembly program even when a reference of a different species is used. Similar improvements can be observed in high and low coverage situations. In addition, we show that a single evaluation metric, like the widely used N50 length, is not enough to properly rate assemblies as it not always points to the best assembly evaluated with other criteria. Therefore, we used the summed z-scores of 36 different statistics to evaluate the assemblies.

Conclusions

The combination of reference mapping and de novo assembly provides a powerful tool to improve genome reconstruction by integrating information of a related genome. Our extension of the reference-guided de novo assembly approach enables the application of this strategy not only within but also between related species. Finally, the evaluation of genome assemblies is often not straight forward, as the truth is not known. Thus one should always use a combination of evaluation metrics, which not only try to assess the continuity but also the accuracy of an assembly.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1911-6) contains supplementary material, which is available to authorized users.

Chloroplast genomes of Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea: Structures and comparative analysis

We investigated the complete chloroplast (cp) genomes of non-model Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea using Illumina paired-end sequencing to understand their genetic organization and structure. Detailed bioinformatics analysis revealed genome sizes of both subspecies ranging between 154.4~154.5 kbp, with a large single-copy region (84,197~84,158 bp), a small single-copy region (17,738~17,813 bp) and pair of inverted repeats (IRa/IRb; 26,264~26,259 bp). Both cp genomes encode 130 genes, including 85 protein-coding genes, eight ribosomal RNA genes and 37 transfer RNA genes. Whole cp genome comparison of A. halleri ssp. gemmifera and A. lyrata ssp. petraea, along with ten other Arabidopsis species, showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. The location and distribution of repeat sequences were determined, and sequence divergences of shared genes were calculated among related species. Comparative phylogenetic analysis of the entire genomic data set and 70 shared genes between both cp genomes confirmed the previous phylogeny and generated phylogenetic trees with the same topologies. The sister species of A. halleri ssp. gemmifera is A. umezawana, whereas the closest relative of A. lyrata spp. petraea is A. arenicola.