Improving the Annotation of Arabidopsis lyrata Using RNA-Seq Data

Differential expression and evolutionary diversification of RNA helicases in Boechera sexual and apomictic reproduction

In higher plants, sexual reproduction is characterized by meiosis of the first cells of the germlines, and double fertilization of the egg and central cell after gametogenesis. In contrast, in apomicts of the genus Boechera, meiosis is omitted or altered and only the central cell requires fertilization, while the embryo forms parthenogenetically from the egg cell. To deepen the understanding of the transcriptional basis underlying these differences, we applied RNA-seq to compare expression in reproductive tissues of different Boechera accessions. This confirmed previous evidence of an enrichment of RNA helicases in plant germlines. Furthermore, few RNA helicases were differentially expressed in female reproductive ovule tissues harboring mature gametophytes from apomictic and sexual accessions. For some of these genes, we further found evidence for a complex recent evolutionary history. This included a homolog of Arabidopsis thaliana FASCIATED STEM4 (FAS4). In contrast to AtFAS4, which is a single-copy gene, FAS4 is represented by three homologs in Boechera, suggesting a potential for subfunctionalization to modulate reproductive development. To gain first insights into functional roles of FAS4, we studied Arabidopsis lines carrying mutant alleles. This identified the crucial importance of AtFAS4 for reproduction, as we observed developmental defects and arrest during male and female gametogenesis.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

Defective pollen tube tip growth induces neo-polyploid infertility

Genome duplication (generating polyploids) is an engine of novelty in eukaryotic evolution and a promising crop improvement tool. Yet newly formed polyploids often have low fertility. Here we report that a severe fertility-compromising defect in pollen tube tip growth arises in new polyploids of Arabidopsis arenosa. Pollen tubes of newly polyploid A. arenosa grow slowly, have aberrant anatomy and disrupted physiology, often burst prematurely, and have altered gene expression. These phenotypes recover in evolved polyploids. We also show that gametophytic (pollen tube) genotypes of two tip-growth genes under selection in natural tetraploid A. arenosa are strongly associated with pollen tube performance in the tetraploid. Our work establishes pollen tube tip growth as an important fertility challenge for neo-polyploid plants and provides insights into a naturally evolved multigenic solution.

Structural evidence for MADS-box type I family expansion seen in new assemblies of Arabidopsis arenosa and A. lyrata

Arabidopsis thaliana diverged from A. arenosa and A. lyrata at least 6 million years ago. The three species differ by genome-wide polymorphisms and morphological traits. The species are to a high degree reproductively isolated, but hybridization barriers are incomplete. A special type of hybridization barrier is based on the triploid endosperm of the seed, where embryo lethality is caused by endosperm failure to support the developing embryo. The MADS-box type I family of transcription factors is specifically expressed in the endosperm and has been proposed to play a role in endosperm-based hybridization barriers. The gene family is well known for its high evolutionary duplication rate, as well as being regulated by genomic imprinting. Here we address MADS-box type I gene family evolution and the role of type I genes in the context of hybridization. Using two de-novo assembled and annotated chromosome-level genomes of A. arenosa and A. lyrata ssp. petraea we analyzed the MADS-box type I gene family in Arabidopsis to predict orthologs, copy number, and structural genomic variation related to the type I loci. Our findings were compared to gene expression profiles sampled before and after the transition to endosperm cellularization in order to investigate the involvement of MADS-box type I loci in endosperm-based hybridization barriers. We observed substantial differences in type-I expression in the endosperm of A. arenosa and A. lyrata ssp. petraea, suggesting a genetic cause for the endosperm-based hybridization barrier between A. arenosa and A. lyrata ssp. petraea.

Role of phenotypic and transcriptomic plasticity in alpine adaptation of Arabidopsis arenosa

Plasticity is an important component of the response of organism to environmental changes, but whether plasticity facilitates adaptation is still largely debated. Using transcriptomic and phenotypic data, we explored the evolution of ancestral plasticity during alpine colonization in Arabidopsis arenosa. We leveraged naturally replicated adaptation in four distinct mountain regions in Central Europe. We sampled seeds from ancestral foothill and independently formed alpine populations in each region and raised them in growth chambers under conditions approximating their natural environments. We gathered RNA-seq and genetic data of 48 and 63 plants and scored vegetative and flowering traits in 203 and 272 plants respectively. Then, we compared gene expression and trait values over two treatments differing in temperature and irradiance and elevations of origin and quantified the extent of ancestral and derived plasticity. At the transcriptomic level, initial plastic changes tended to be more reinforced than reversed in adapted alpine populations. Genes showing reinforcement were involved in the stress response, developmental processes and morphogenesis and those undergoing reversion were related to the stress response (light and biotic stress). At the phenotypic level, initial plastic changes in all but one trait were also reinforced supporting a facilitating role of phenotypic plasticity during colonization of an alpine environment. Our results contrasted with previous studies that showed generally higher reversion than reinforcement and supported the idea that ancestral plasticity tends to be reinforced in the context of alpine adaptation. However, plasticity may also be the source of potential maladaptation, especially at the transcriptomic level.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

Multiomics Reveals the Regulatory Mechanisms of Arabidopsis Tissues under Heat Stress

Understanding the mechanisms of responses to high temperatures in Arabidopsis will provide insights into how plants may mitigate heat stress under global climate change. And exploring the interconnections of different modification levels in heat stress response could help us to understand the molecular mechanism of heat stress response in Arabidopsis more comprehensively and precisely. In this paper, we combined multiomics analyses to explore the common heat stress-responsive genes and specific heat-responsive metabolic pathways in Arabidopsis leaf, seedling, and seed tissues. We found that genes such as AT1G54050 play a role in promoting proper protein folding in response to HS (Heat stress). In addition, it was revealed that the binding profile of A1B is altered under elevated temperature conditions. Finally, we also show that two microRNAs, ath-mir156h and ath-mir166b-5p, may be core regulatory molecules in HS. Also elucidated that under HS, plants can regulate specific regulatory mechanisms, such as oxygen levels, by altering the degree of CHH methylation.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Shedding light on AT1G29480 of Arabidopsis thaliana-An enigmatic locus restricted to Brassicacean genomes

A key feature of C₄ Kranz anatomy is the presence of an enlarged, photosynthetically highly active bundle sheath whose cells contain large numbers of chloroplasts. With the aim to identify novel candidate regulators of C₄ bundle sheath development, we performed an activation tagging screen with Arabidopsis thaliana. The reporter gene used encoded a chloroplast-targeted GFP protein preferentially expressed in the bundle sheath, and the promoter of the C₄ phosphoenolpyruvate carboxylase gene from Flaveria trinervia served as activation tag because of its activity in all chlorenchymatous tissues of A. thaliana. Primary mutants were selected based on their GFP signal intensity, and one stable mutant named kb-1 with a significant increase in GFP fluorescence intensity was obtained. Despite the increased GFP signal, kb-1 showed no alterations to bundle sheath anatomy. The causal locus, AT1G29480, is specific to the Brassicaceae with its second exon being conserved. Overexpression and reconstitution studies confirmed that AT1G29480, and specifically its second exon, were sufficient for the enhanced GFP phenotype, which was not dependent on translation of the locus or its parts into protein. We conclude, therefore, that the AT1G29480 locus enhances the GFP reporter gene activity via an RNA-based mechanism.

An expanded role for the transcription factor WRINKLED1 in the biosynthesis of triacylglycerols during seed development

The transcription factor WRINKLED1 (WRI1) is known as a master regulator of fatty acid synthesis in developing oilseeds of Arabidopsis thaliana and other species. WRI1 is known to directly stimulate the expression of many fatty acid biosynthetic enzymes and a few targets in the lower part of the glycolytic pathway. However, it remains unclear to what extent and how the conversion of sugars into fatty acid biosynthetic precursors is controlled by WRI1. To shortlist possible gene targets for future in-planta experimental validation, here we present a strategy that combines phylogenetic foot printing of cis-regulatory elements with additional layers of evidence. Upstream regions of protein-encoding genes in A. thaliana were searched for the previously described DNA-binding consensus for WRI1, the ASML1/WRI1 (AW)-box. For about 900 genes, AW-box sites were found to be conserved across orthologous upstream regions in 11 related species of the crucifer family. For 145 select potential target genes identified this way, affinity of upstream AW-box sequences to WRI1 was assayed by Microscale Thermophoresis. This allowed definition of a refined WRI1 DNA-binding consensus. We find that known WRI1 gene targets are predictable with good confidence when upstream AW-sites are phylogenetically conserved, specifically binding WRI1 in the in vitro assay, positioned in proximity to the transcriptional start site, and if the gene is co-expressed with WRI1 during seed development. When targets predicted in this way are mapped to central metabolism, a conserved regulatory blueprint emerges that infers concerted control of contiguous pathway sections in glycolysis and fatty acid biosynthesis by WRI1. Several of the newly predicted targets are in the upper glycolysis pathway and the pentose phosphate pathway. Of these, plastidic isoforms of fructokinase (FRK3) and of phosphoglucose isomerase (PGI1) are particularly corroborated by previously reported seed phenotypes of respective null mutations.

High transcriptome plasticity drives phosphate starvation responses in tomato

Tomato is an important vegetable crop and fluctuating available soil phosphate (Pi) level elicits several morpho-physiological responses driven by underlying molecular responses. Therefore, understanding these molecular responses at the gene and isoform levels has become critical in the quest for developing crops with improved Pi use efficiency. A quantitative time-series RNA-seq analysis was performed to decipher the global transcriptomic changes that accompany Pi starvation in tomato. Apart from changes in the expression levels of genes, there were also alterations in the expression of alternatively-spliced transcripts. Physiological responses such as anthocyanin accumulation, reactive oxygen species generation and cell death are obvious 7 days after Pi deprivation accompanied with the maximum amount of transcriptional change in the genome making it an important stage for in-depth study while studying Pi stress responses (PSR). Our study demonstrates that transcriptomic changes under Pi deficiency are dynamic and complex in tomato. Overall, our study dwells on the dynamism of the transcriptome in eliciting a response to adapt to low Pi stress and lays it bare. Findings from this study will prove to be an invaluable resource for researchers using tomato as a model for understanding nutrient deficiency.

Male and female recombination landscapes of diploid Arabidopsis arenosa

The number and placement of meiotic crossover events during meiosis have important implications for the fidelity of chromosome segregation as well as patterns of inheritance. Despite the functional importance of recombination, recombination landscapes vary widely among and within species, and this can have a strong impact on evolutionary processes. A good knowledge of recombination landscapes is important for model systems in evolutionary and ecological genetics, since it can improve interpretation of genomic patterns of differentiation and genome evolution, and provides an important starting point for understanding the causes and consequences of recombination rate variation. Arabidopsis arenosa is a powerful evolutionary genetic model for studying the molecular basis of adaptation and recombination rate evolution. Here, we generate genetic maps for 2 diploid A. arenosa individuals from distinct genetic lineages where we have prior knowledge that meiotic genes show evidence of selection. We complement the genetic maps with cytological approaches to map and quantify recombination rates, and test the idea that these populations might have distinct patterns of recombination. We explore how recombination differs at the level of populations, individuals, sexes and genomic regions. We show that the positioning of crossovers along a chromosome correlates with their number, presumably a consequence of crossover interference, and discuss how this effect can cause differences in recombination landscape among sexes or species. We identify several instances of female segregation distortion. We found that averaged genome-wide recombination rate is lower and sex differences subtler in A. arenosa than in Arabidopsis thaliana.

Parallel adaptation in autopolyploid Arabidopsis arenosa is dominated by repeated recruitment of shared alleles

Relative contributions of pre-existing vs de novo genomic variation to adaptation are poorly understood, especially in polyploid organisms. We assess this in high resolution using autotetraploid Arabidopsis arenosa, which repeatedly adapted to toxic serpentine soils that exhibit skewed elemental profiles. Leveraging a fivefold replicated serpentine invasion, we assess selection on SNPs and structural variants (TEs) in 78 resequenced individuals and discover significant parallelism in candidate genes involved in ion homeostasis. We further model parallel selection and infer repeated sweeps on a shared pool of variants in nearly all these loci, supporting theoretical expectations. A single striking exception is represented by TWO PORE CHANNEL 1, which exhibits convergent evolution from independent de novo mutations at an identical, otherwise conserved site at the calcium channel selectivity gate. Taken together, this suggests that polyploid populations can rapidly adapt to environmental extremes, calling on both pre-existing variation and novel polymorphisms.

Relative contributions of pre-existing versus de novo genomic variation to adaptation remain unclear. Here, the authors address this problem by examining the adaptation of autotetraploid Arabidopsis arenosa to serpentine soils and find that both types of variations contribute to rapid adaptation.

A high-quality draft genome for Melaleuca alternifolia (tea tree): a new platform for evolutionary genomics of myrtaceous terpene-rich species

The economically important Melaleuca alternifolia (tea tree) is the source of a terpene-rich essential oil with therapeutic and cosmetic uses around the world. Tea tree has been cultivated and bred in Australia since the 1990s. It has been extensively studied for the genetics and biochemistry of terpene biosynthesis. Here, we report a high quality de novo genome assembly using Pacific Biosciences and Illumina sequencing. The genome was assembled into 3128 scaffolds with a total length of 362 Mb (N50  = 1.9 Mb), with significantly higher contiguity than a previous assembly (N50  = 8.7 Kb). Using a homology-based, RNA-seq evidence-based and ab initio prediction approach, 37,226 protein-coding genes were predicted. Genome assembly and annotation exhibited high completeness scores of 98.1% and 89.4%, respectively. Sequence contiguity was sufficient to reveal extensive gene order conservation and chromosomal rearrangements in alignments with Eucalyptus grandis and Corymbia citriodora genomes. This new genome advances currently available resources to investigate the genome structure and gene family evolution of M. alternifolia. It will enable further comparative genomic studies in Myrtaceae to elucidate the genetic foundations of economically valuable traits in this crop.

De Novo Mutation and Rapid Protein (Co-)evolution during Meiotic Adaptation in Arabidopsis arenosa

A sudden shift in environment or cellular context necessitates rapid adaptation. A dramatic example is genome duplication, which leads to polyploidy. In such situations, the waiting time for new mutations might be prohibitive; theoretical and empirical studies suggest that rapid adaptation will largely rely on standing variation already present in source populations. Here, we investigate the evolution of meiosis proteins in Arabidopsis arenosa, some of which were previously implicated in adaptation to polyploidy, and in a diploid, habitat. A striking and unexplained feature of prior results was the large number of amino acid changes in multiple interacting proteins, especially in the relatively young tetraploid. Here, we investigate whether selection on meiosis genes is found in other lineages, how the polyploid may have accumulated so many differences, and whether derived variants were selected from standing variation. We use a range-wide sample of 145 resequenced genomes of diploid and tetraploid A. arenosa, with new genome assemblies. We confirmed signals of positive selection in the polyploid and diploid lineages they were previously reported in and find additional meiosis genes with evidence of selection. We show that the polyploid lineage stands out both qualitatively and quantitatively. Compared with diploids, meiosis proteins in the polyploid have more amino acid changes and a higher proportion affecting more strongly conserved sites. We find evidence that in tetraploids, positive selection may have commonly acted on de novo mutations. Several tests provide hints that coevolution, and in some cases, multinucleotide mutations, might contribute to rapid accumulation of changes in meiotic proteins.

Drivers of linkage disequilibrium across a species' geographic range

While linkage disequilibrium (LD) is an important parameter in genetics and evolutionary biology, the drivers of LD remain elusive. Using whole-genome sequences from across a species’ range, we assessed the impact of demographic history and mating system on LD. Both range expansion and a shift from outcrossing to selfing in North American Arabidopsis lyrata were associated with increased average genome-wide LD. Our results indicate that range expansion increases short-distance LD at the farthest range edges by about the same amount as a shift to selfing. However, the extent over which LD in genic regions unfolds was shorter for range expansion compared to selfing. Linkage among putatively neutral variants and between neutral and deleterious variants increased to a similar degree with range expansion, providing support that genome-wide LD was positively associated with mutational load. As a consequence, LD combined with mutational load may decelerate range expansions and set range limits. Finally, a small number of genes were identified as LD outliers, suggesting that they experience selection by either of the two demographic processes. These included genes involved in flowering and photoperiod for range expansion, and the self-incompatibility locus for mating system.

Nearby genomic variants are often co-inherited because of limited recombination. The extent of non-random association of alleles at different loci is called linkage disequilibrium (LD) and is commonly used in genomic analyses, for example to detect regions under selection or to determine effective population size. Here we reversed testing and addressed how demographic history may affect LD within a species. Using genomic data from more than a thousand individuals of North American Arabidopsis lyrata from across the entire species’ range, we quantified the effect of postglacial range expansion and a shift in mating system from outcrossing to selfing on LD. We show that both factors lead to increased LD, and that the maximal effect of range expansion is comparable with a shift in mating system to selfing. Heightened LD involves deleterious mutations, and therefore, LD can also serve as an indicator of mutation accumulation. Furthermore, we provide evidence that some genes experienced stronger increases in LD possibly due to selection associated with the two demographic changes. Our results provide a novel and broad view on the evolutionary factors shaping LD that may also apply to the very many species that underwent postglacial range expansion.

Parallelism in gene expression between foothill and alpine ecotypes in Arabidopsis arenosa

Parallel adaptation results from the independent evolution of similar traits between closely related lineages and allows us to test to what extent evolution is repeatable. Similar gene expression changes are often detected but the identity of genes shaped by parallel selection and the causes of expression parallelism remain largely unknown. By comparing genomes and transcriptomes of four distinct foothill-alpine population pairs across four treatments, we addressed the genetic underpinnings, plasticity and functional consequences of gene expression parallelism in alpine adaptation. Seeds of eight populations of Arabidopsis arenosa were raised under four treatments that differed in temperature and irradiance, factors varying strongly with elevation. Parallelism in differential gene expression between the foothill and alpine ecotypes was quantified by RNA-seq in leaves of young plants. By manipulating temperature and irradiance, we also tested for parallelism in plasticity (i.e., gene-environment interaction, GEI). In spite of global non-parallel patterns transcriptome wide, we found significant parallelism in gene expression at the level of individual loci with an over-representation of genes involved in biotic stress response. In addition, we demonstrated significant parallelism in GEI, indicating a shared differential response of the originally foothill versus alpine populations to environmental variation across mountain regions. A fraction of genes showing expression parallelism also encompassed parallel outliers for genomic differentiation, with greater enrichment of such variants in cis-regulatory elements in some mountain regions. In summary, our results suggest frequent evolutionary repeatability in gene expression changes associated with the colonization of a challenging environment that combines constitutive expression differences and plastic interaction with the surrounding environment.

New model species for arctic-alpine plant molecular ecology

Arctic and alpine, high latitude and high elevation environments are one of the most stressful environments for species to inhabit. This harshness manifests itself in lower species richness in comparison to more southern vegetation zones (Francis & Currie, 2003). Furthermore, the climatic oscillations-past and predicted-have the most dramatic effect on these ecosystems. For example, in regions of continental ice sheets-the northernmost part of Western Europe and North America-the Arctic species assemblages are no older than a few thousands of years, which is a relatively short period from an evolutionary perspective. Although similar environments may have existed further south during the Ice Age, allowing some preadaptation for the Arctic species, the current habitat is a unique combination of environmental factors such as the climate, soil, bedrock, and photoperiod. Hence, understanding the evolutionary forces shaping Arctic-alpine species will be important for predicting these vulnerable environments' population viability and adaptive potential in the future. In this issue of Molecular Ecology Resources, Nowak et al. (Molecular Ecology Resources) present extensive genome-wide resources for an Arctic-alpine plant Draba nivalis. This adds a valuable new member into the cabbage family models for evolutionary genetics and adaptation studies, to accompany e.g., Arabidopsis (Nature Genetics, 43, 476; Nature, 408, 796), Arabis (Nature Plants, 1, 14023) and Capsella (Nature Genetics, 45, 831). A whole new avenue will open up for molecular ecological studies not only for D. nivalis, but the whole large Draba genus with its diverse ecological and evolutionary characteristics.

Gradual evolution of allopolyploidy in Arabidopsis suecica

Most diploid organisms have polyploid ancestors. The evolutionary process of polyploidization is poorly understood but has frequently been conjectured to involve some form of 'genome shock', such as genome reorganization and subgenome expression dominance. Here we study polyploidization in Arabidopsis suecica, a post-glacial allopolyploid species formed via hybridization of Arabidopsis thaliana and Arabidopsis arenosa. We generated a chromosome-level genome assembly of A. suecica and complemented it with polymorphism and transcriptome data from all species. Despite a divergence around 6 million years ago (Ma) between the ancestral species and differences in their genome composition, we see no evidence of a genome shock: the A. suecica genome is colinear with the ancestral genomes; there is no subgenome dominance in expression; and transposon dynamics appear stable. However, we find changes suggesting gradual adaptation to polyploidy. In particular, the A. thaliana subgenome shows upregulation of meiosis-related genes, possibly to prevent aneuploidy and undesirable homeologous exchanges that are observed in synthetic A. suecica, and the A. arenosa subgenome shows upregulation of cyto-nuclear processes, possibly in response to the new cytoplasmic environment of A. suecica, with plastids maternally inherited from A. thaliana. These changes are not seen in synthetic hybrids, and thus are likely to represent subsequent evolution.

Demographic Processes Linked to Genetic Diversity and Positive Selection across a Species' Range

Demography determines the strength of genetic drift, which generally reduces genetic variation and the efficacy of selection. Here, we disentangled the importance of demographic processes at a local scale (census size and mating system) and at a species-range scale (old split between population clusters, recolonization after the last glaciation cycle, and admixture) in determining within-population genomic diversity and genomic signatures of positive selection. Analyses were based on re-sequence data from 52 populations of North American Arabidopsis lyrata collected across its entire distribution. The mating system and range dynamics since the last glaciation cycle explained around 60% of the variation in genomic diversity among populations and 52% of the variation in the signature of positive selection. Diversity was lowest in selfing compared with outcrossing populations and in areas further away from glacial refugia. In parallel, reduced positive selection was found in selfing populations and in populations with a longer route of postglacial range expansion. The signature of positive selection was also reduced in populations without admixture. We conclude that recent range expansion can have a profound influence on diversity in coding and non-coding DNA, similar in magnitude to the shift toward selfing. Distribution limits may in fact be caused by reduced effective population size and compromised positive selection in recently colonized parts of the range.

Resistance Gene Analogs in the Brassicaceae: Identification, Characterization, Distribution, and Evolution

The Brassicaceae consists of a wide range of species, including important Brassica crop species and the model plant Arabidopsis (Arabidopsis thaliana). Brassica spp. crop diseases impose significant yield losses annually. A major way to reduce susceptibility to disease is the selection in breeding for resistance gene analogs (RGAs). Nucleotide binding site-leucine rich repeats (NLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are the main types of RGAs; they contain conserved domains and motifs and play specific roles in resistance to pathogens. Here, all classes of RGAs have been identified using annotation and assembly-based pipelines in all available genome annotations from the Brassicaceae, including multiple genome assemblies of the same species where available (total of 32 genomes). The number of RGAs, based on genome annotations, varies within and between species. In total 34,065 RGAs were identified, with the majority being RLKs (21,691), then NLRs (8,588) and RLPs (3,786). Analysis of the RGA protein sequences revealed a high level of sequence identity, whereby 99.43% of RGAs fell into several orthogroups. This study establishes a resource for the identification and characterization of RGAs in the Brassicaceae and provides a framework for further studies of RGAs for an ultimate goal of assisting breeders in improving resistance to plant disease.

The Genome of Microthlaspi erraticum (Brassicaceae) Provides Insights Into the Adaptation to Highly Calcareous Soils

Microthlaspi erraticum is widely distributed in temperate Eurasia, but restricted to Ca²⁺-rich habitats, predominantly on white Jurassic limestone, which is made up by calcium carbonate, with little other minerals. Thus, naturally occurring Microthlaspi erraticum individuals are confronted with a high concentration of Ca²⁺ ions while Mg²⁺ ion concentration is relatively low. As there is a competitive uptake between these two ions, adaptation to the soil condition can be expected. In this study, it was the aim to explore the genomic consequences of this adaptation by sequencing and analysing the genome of Microthlaspi erraticum. Its genome size is comparable with other diploid Brassicaceae, while more genes were predicted. Two Mg²⁺ transporters known to be expressed in roots were duplicated and one showed a significant degree of positive selection. It is speculated that this evolved due to the pressure to take up Mg²⁺ ions efficiently in the presence of an overwhelming amount of Ca²⁺ ions. Future studies on plants specialized on similar soils and affinity tests of the transporters are needed to provide unequivocal evidence for this hypothesis. If verified, the transporters found in this study might be useful for breeding Brassicaceae crops for higher yield on Ca²⁺-rich and Mg²⁺ -poor soils.

mRNA Extraction from Gill Tissue for RNA-sequencing

Adaptation is thought to proceed in part through spatial and temporal changes in gene expression. Fish species such as the threespine stickleback are powerful vertebrate models to study the genetic architecture of adaptive changes in gene expression since divergent adaptation to different environments is common, they are abundant and easy to study in the wild and lab, and have well-established genetic and genomic resources. Fish gills, due to their respiratory and osmoregulatory roles, show many physiological adaptations to local water chemistry, including differences in gene expression. However, obtaining high-quality RNA using popular column-based extraction methods can be challenging from small tissue samples high in cartilage and bone such as fish gills. Here, we describe a bead-based mRNA extraction and transcriptome RNA-seq protocol that does not use purification columns. The protocol can be readily scaled according to sample size for the purposes of diverse gene expression experiments using animal or plant tissue.

Protein lysine methylation contributes to modulating the response of sensitive and tolerant Arabidopsis species to cadmium stress

The mechanisms underlying the response and adaptation of plants to excess of trace elements are not fully described. Here, we analysed the importance of protein lysine methylation for plants to cope with cadmium. We analysed the effect of cadmium on lysine-methylated proteins and protein lysine methyltransferases (KMTs) in two cadmium-sensitive species, Arabidopsis thaliana and A. lyrata, and in three populations of A. halleri with contrasting cadmium accumulation and tolerance traits. We showed that some proteins are differentially methylated at lysine residues in response to Cd and that a few genes coding KMTs are regulated by cadmium. Also, we showed that 9 out of 23 A. thaliana mutants disrupted in KMT genes have a tolerance to cadmium that is significantly different from that of wild-type seedlings. We further characterized two of these mutants, one was knocked out in the calmodulin lysine methyltransferase gene and displayed increased tolerance to cadmium, and the other was interrupted in a KMT gene of unknown function and showed a decreased capacity to cope with cadmium. Together, our results showed that lysine methylation of non-histone proteins is impacted by cadmium and that several methylation events are important for modulating the response of Arabidopsis plants to cadmium stress.

Biomolecular approaches to understanding metal tolerance and hyperaccumulation in plants

Trace metal elements are essential for plant growth but become toxic at high concentrations, while some non-essential elements, such as Cd and As, show toxicity even in traces.

What Can Long Terminal Repeats Tell Us About the Age of LTR Retrotransposons, Gene Conversion and Ectopic Recombination?

LTR retrotransposons constitute a significant part of plant genomes and their evolutionary dynamics play an important role in genome size changes. Current methods of LTR retrotransposon age estimation are based only on LTR (long terminal repeat) divergence. This has prompted us to analyze sequence similarity of LTRs in 25,144 LTR retrotransposons from fifteen plant species as well as formation of solo LTRs. We found that approximately one fourth of nested retrotransposons showed a higher LTR divergence than the pre-existing retrotransposons into which they had been inserted. Moreover, LTR similarity was correlated with LTR length. We propose that gene conversion can contribute to this phenomenon. Gene conversion prediction in LTRs showed potential converted regions in 25% of LTR pairs. Gene conversion was higher in species with smaller genomes while the proportion of solo LTRs did not change with genome size in analyzed species. The negative correlation between the extent of gene conversion and the abundance of solo LTRs suggests interference between gene conversion and ectopic recombination. Since such phenomena limit the traditional methods of LTR retrotransposon age estimation, we recommend an improved approach based on the exclusion of regions affected by gene conversion.

Relaxed purifying selection in autopolyploids drives transposable element over-accumulation which provides variants for local adaptation

Polyploidization is frequently associated with increased transposable element (TE) content. However, what drives TE dynamics following whole genome duplication (WGD) and the evolutionary implications remain unclear. Here, we leverage whole-genome resequencing data available for ~300 individuals of Arabidopsis arenosa, a well characterized natural diploid-autotetraploid plant species, to address these questions. Based on 43,176 TE insertions we detect in these genomes, we demonstrate that relaxed purifying selection rather than transposition bursts is the main driver of TE over-accumulation after WGD. Furthermore, the increased pool of TE insertions in tetraploids is especially enriched within or near environmentally responsive genes. Notably, we show that the major flowering-time repressor gene FLC is disrupted by a TE insertion specifically in the rapid-cycling tetraploid lineage that colonized mainland railways. Together, our findings indicate that tetrasomy leads to an enhanced accumulation of genic TE insertions, some of which likely contribute to local adaptation.

Why transposable elements (TEs) accumulate in polyploids and the evolutionary implications remain unclear. Here, the authors show that following whole genome duplication, relaxed purifying selection is the main driver of TE over-accumulation, which provides variants for rapid local adaptation.

Nested plant LTR retrotransposons target specific regions of other elements, while all LTR retrotransposons often target palindromes and nucleosome-occupied regions: in silico study

Background

Nesting is common in LTR retrotransposons, especially in large genomes containing a high number of elements.

Results

We analyzed 12 plant genomes and obtained 1491 pairs of nested and original (pre-existing) LTR retrotransposons. We systematically analyzed mutual nesting of individual LTR retrotransposons and found that certain families, more often belonging to the Ty3/gypsy than Ty1/copia superfamilies, showed a higher nesting frequency as well as a higher preference for older copies of the same family ("autoinsertions"). Nested LTR retrotransposons were preferentially located in the 3'UTR of other LTR retrotransposons, while coding and regulatory regions (LTRs) are not commonly targeted. Insertions displayed a weak preference for palindromes and were associated with a strong positional pattern of higher predicted nucleosome occupancy. Deviation from randomness in target site choice was also found in 13,983 non-nested plant LTR retrotransposons.

Conclusions

We reveal that nesting of LTR retrotransposons is not random. Integration is correlated with sequence composition, secondary structure and the chromatin environment. Insertion into retrotransposon positions with a low negative impact on family fitness supports the concept of the genome being viewed as an ecosystem of various elements.

Postglacial ecotype formation under outcrossing and self-fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate-related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole-genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate-related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate-related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome-wide linkage disequilibrium.

Contribution of Functional Divergence Through Copy Number Variations to the Inter-Species and Intra-Species Diversity in Specialized Metabolites

There is considerable diversity in the specialized metabolites within a single plant species (intra-species) and among different plant species (inter-species). The functional divergence associated with gene duplications largely contributes to the inter-species diversity in the specialized metabolites, whereas the intra-species diversity is due to gene dosage changes via gene duplications [i.e., copy number variants (CNVs)] at the intra-species level of evolution. This is because CNVs are thought to undergo associated with less functional divergence at the intra-species level of evolution. However, functional divergence caused by CNVs may induce specialized metabolite diversity at the intra-species and inter-species levels of evolution. We herein discuss the functional divergence of CNVs in metabolic quantitative trait genes (mQTGs). We focused on 5,654 previously identified mQTGs in 270 Arabidopsis thaliana accessions. The ratio of nonsynonymous to synonymous variations tends to be higher for mQTGs with CNVs than for mQTGs without CNVs within A. thaliana accessions, suggesting that CNVs are responsible for the functional divergence among mQTGs at the intra-species level of evolution. To evaluate the contribution of CNVs to inter-species diversity, we calculated the ratio of nonsynonymous to synonymous substitutions in the Arabidopsis lineage. The ratio tends to be higher for the mQTGs with CNVs than for the mQTGs without CNVs. Additionally, we determined that mQTGs with CNVs are subject to positive selection in the Arabidopsis lineage. Our data suggest that CNVs are closely related to functional divergence contributing to adaptations via the production of diverse specialized metabolites at the intra-species and inter-species levels of evolution.

Interspecific introgression mediates adaptation to whole genome duplication

Adaptive gene flow is a consequential phenomenon across all kingdoms. Although recognition is increasing, there is no study showing that bidirectional gene flow mediates adaptation at loci that manage core processes. We previously discovered concerted molecular changes among interacting members of the meiotic machinery controlling crossover number upon adaptation to whole-genome duplication (WGD) in Arabidopsis arenosa. Here we conduct a population genomic study to test the hypothesis that adaptation to WGD has been mediated by adaptive gene flow between A. arenosa and A. lyrata. We find that A. lyrata underwent WGD more recently than A. arenosa, suggesting that pre-adapted alleles have rescued nascent A. lyrata, but we also detect gene flow in the opposite direction at functionally interacting loci under the most extreme levels of selection. These data indicate that bidirectional gene flow allowed for survival after WGD, and that the merger of these species is greater than the sum of their parts.

Whole genome duplication (WGD) presents new challenges to the establishment of optimal allelic combinations and to the meiotic machinery. Here, the authors show that adaptive gene flow from Arabidopsis arenosa could rescue the nascent A. lyrata from extinction following WGD.

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'.

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'.

A stigmatic gene confers interspecies incompatibility in the Brassicaceae

Pre-zygotic interspecies incompatibility in angiosperms is a male-female relationship that inhibits the formation of hybrids between two species. Here, we report on the identification of STIGMATIC PRIVACY 1 (SPRI1), an interspecies barrier gene in Arabidopsis thaliana. We show that the rejection activity of this stigma-specific plasma membrane protein is effective against distantly related Brassicaceae pollen tubes and is independent of self-incompatibility. Point-mutation experiments and functional tests of synthesized hypothetical ancestral forms of SPRI1 suggest evolutionary decay of SPRI1-controlled interspecies incompatibility in self-compatible A. thaliana. Hetero-pollination experiments indicate that SPRI1 ensures intraspecific fertilization in the pistil when pollen from other species are present. Our study supports the idea that SPRI1 functions as a barrier mechanism that permits entrance of pollen with an intrinsic signal from self species.

Adaptive and degenerative evolution of the S-Phase Kinase-Associated Protein 1-Like family in Arabidopsis thaliana

Genome sequencing has uncovered tremendous sequence variation within and between species. In plants, in addition to large variations in genome size, a great deal of sequence polymorphism is also evident in several large multi-gene families, including those involved in the ubiquitin-26S proteasome protein degradation system. However, the biological function of this sequence variation is yet not clear. In this work, we explicitly demonstrated a single origin of retroposed Arabidopsis Skp1-Like (ASK) genes using an improved phylogenetic analysis. Taking advantage of the 1,001 genomes project, we here provide several lines of polymorphism evidence showing both adaptive and degenerative evolutionary processes in ASK genes. Yeast two-hybrid quantitative interaction assays further suggested that recent neutral changes in the ASK2 coding sequence weakened its interactions with some F-box proteins. The trend that highly polymorphic upstream regions of ASK1 yield high levels of expression implied negative expression regulation of ASK1 by an as-yet-unknown transcriptional suppression mechanism, which may contribute to the polymorphic roles of Skp1-CUL1-F-box complexes. Taken together, this study provides new evolutionary evidence to guide future functional genomic studies of SCF-mediated protein ubiquitylation.

Contribution of time of day and the circadian clock to the heat stress responsive transcriptome in Arabidopsis

In Arabidopsis, a large subset of heat responsive genes exhibits diurnal or circadian oscillations. However, to what extent the dimension of time and/or the circadian clock contribute to heat stress responses remains largely unknown. To determine the direct contribution of time of day and/or the clock to differential heat stress responses, we probed wild-type and mutants of the circadian clock genes CCA1, LHY, PRR7, and PRR9 following exposure to heat (37 °C) and moderate cold (10 °C) in the early morning (ZT1) and afternoon (ZT6). Thousands of genes were differentially expressed in response to temperature, time of day, and/or the clock mutation. Approximately 30% more genes were differentially expressed in the afternoon compared to the morning, and heat stress significantly perturbed the transcriptome. Of the DEGs (~3000) specifically responsive to heat stress, ~70% showed time of day (ZT1 or ZT6) occurrence of the transcriptional response. For the DEGs (~1400) that are shared between ZT1 and ZT6, we observed changes to the magnitude of the transcriptional response. In addition, ~2% of all DEGs showed differential responses to temperature stress in the clock mutants. The findings in this study highlight a significant role for time of day in the heat stress responsive transcriptome, and the clock through CCA1 and LHY, appears to have a more profound role than PRR7 and PRR9 in modulating heat stress responses during the day. Our results emphasize the importance of considering the dimension of time in studies on abiotic stress responses in Arabidopsis.

Accumulation of Mutational Load at the Edges of a Species Range

Why species have geographically restricted distributions is an unresolved question in ecology and evolutionary biology. Here, we test a new explanation that mutation accumulation due to small population size or a history of range expansion can contribute to restricting distributions by reducing population growth rate at the edge. We examined genomic diversity and mutational load across the entire geographic range of the North American plant Arabidopsis lyrata, including old, isolated populations predominantly at the southern edge and regions of postglacial range expansion at the northern and southern edges. Genomic diversity in intergenic regions declined toward distribution edges and signatures of mutational load in exon regions increased. Genomic signatures of mutational load were highly linked to phenotypically expressed load, measured as reduced performance of individual plants and lower estimated rate of population growth. The geographic pattern of load and the connection between load and population growth demonstrate that mutation accumulation reduces fitness at the edge and helps restrict species' distributions.

Divergence of RNA editing among Arabidopsis species

RNA editing altered the RNA sequence by replacing the C nucleotide to U in the organellar genomes of plants. RNA editing status sometimes differed among distant species. The pattern of conservation and variation of RNA editing status made it possible to evaluate evolutionary mechanisms impacting functional aspects of RNA editing. In this study, divergence of RNA editing in the chloroplast genome among Arabidopsis species was analyzed to determine 9 losses and 1 gain in RNA editing. All changes in A. thaliana lineage resulted from changes to the chloroplast genome sequence, whereas changes in the A. lyrata / halleri lineage were possibly due to exclusive changes in the nuclear editing factors. One loss of RNA editing in A. lyrata was caused by a deficiency in the PPR gene OTP80. The changes in RNA editing occurred approximately every two million years and were not observed at functionally important sites. These results highlight the conserved nature of RNA editing status suggesting the importance of RNA editing during evolution.

The genomic basis of adaptation to calcareous and siliceous soils in Arabidopsis lyrata

Edaphic conditions are important determinants of plant fitness. While much has been learnt in recent years about plant adaptation to heavy metal contaminated soils, the genomic basis underlying adaptation to calcareous and siliceous substrates remains largely unknown. We performed a reciprocal germination experiment and whole-genome resequencing in natural calcareous and siliceous populations of diploid Arabidopsis lyrata to test for edaphic adaptation and detect signatures of selection at loci associated with soil-mediated divergence. In parallel, genome scans on respective diploid ecotypes from the Arabidopsis arenosa species complex were undertaken, to search for shared patterns of adaptive genetic divergence. Soil ecotypes of A. lyrata display significant genotype-by-treatment responses for seed germination. Sequence (SNPs) and copy-number variants (CNVs) point towards loci involved in ion transport as the main targets of adaptive genetic divergence. Two genes exhibiting high differentiation among soil types in A. lyrata further share trans-specific single nucleotide polymorphisms with A. arenosa. This work applies experimental and genomic approaches to study edaphic adaptation in A. lyrata and suggests that physiological response to elemental toxicity and deficiency underlies the evolution of calcareous and siliceous ecotypes. The discovery of shared adaptive variation between sister species indicates that ancient polymorphisms contribute to adaptive evolution.

Genetic differentiation in life history traits and thermal stress performance across a heterogeneous dune landscape in Arabidopsis lyrata

Background and Aims

Over very short spatial scales, the habitat of a species can differ in multiple abiotic and biotic factors. These factors may impose natural selection on several traits and can cause genetic differentiation within a population. We studied multivariate genetic differentiation in a plant species of a sand dune landscape by linking environmental variation with differences in genotypic trait values and gene expression levels to find traits and candidate genes of microgeographical adaptation.

Methods

Maternal seed families of Arabidopsis lyrata were collected in Saugatuck Dunes State Park, Michigan, USA, and environmental parameters were recorded at each collection site. Offspring plants were raised in climate chambers and exposed to one of three temperature treatments: regular occurrence of frost, heat, or constant control conditions. Several traits were assessed: plant growth, time to flowering, and frost and heat resistance.

Key Results

The strongest trait-environment association was between a fast switch to sexual reproduction and weaker growth under frost, and growing in the open, away from trees. The second strongest association was between the trait combination of small plant size and early flowering under control conditions combined with large size under frost, and the combination of environmental conditions of growing close to trees, at low vegetation cover, on dune bottoms. Gene expression analysis by RNA-seq revealed candidate genes involved in multivariate trait differentiation.

Conclusions

The results support the hypothesis that in natural populations, many environmental factors impose selection, and that they affect multiple traits, with the relative direction of trait change being complex. The results highlight that heterogeneity in the selection environment over small spatial scales is a main driver of the maintenance of adaptive genetic variation within populations.

Combining RNA-seq data and homology-based gene prediction for plants, animals and fungi

Background

Genome annotation is of key importance in many research questions. The identification of protein-coding genes is often based on transcriptome sequencing data, ab-initio or homology-based prediction. Recently, it was demonstrated that intron position conservation improves homology-based gene prediction, and that experimental data improves ab-initio gene prediction.

Results

Here, we present an extension of the gene prediction program GeMoMa that utilizes amino acid sequence conservation, intron position conservation and optionally RNA-seq data for homology-based gene prediction. We show on published benchmark data for plants, animals and fungi that GeMoMa performs better than the gene prediction programs BRAKER1, MAKER2, and CodingQuarry, and purely RNA-seq-based pipelines for transcript identification. In addition, we demonstrate that using multiple reference organisms may help to further improve the performance of GeMoMa. Finally, we apply GeMoMa to four nematode species and to the recently published barley reference genome indicating that current annotations of protein-coding genes may be refined using GeMoMa predictions.

Conclusions

GeMoMa might be of great utility for annotating newly sequenced genomes but also for finding homologs of a specific gene or gene family. GeMoMa has been published under GNU GPL3 and is freely available at http://www.jstacs.de/index.php/GeMoMa.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2203-5) contains supplementary material, which is available to authorized users.

Multi-Omics Driven Assembly and Annotation of the Sandalwood (Santalum album) Genome

Indian sandalwood (Santalum album) is an important tropical evergreen tree known for its fragrant heartwood-derived essential oil and its valuable carving wood. Here, we applied an integrated genomic, transcriptomic, and proteomic approach to assemble and annotate the Indian sandalwood genome. Our genome sequencing resulted in the establishment of a draft map of the smallest genome for any woody tree species to date (221 Mb). The genome annotation predicted 38,119 protein-coding genes and 27.42% repetitive DNA elements. In-depth proteome analysis revealed the identities of 72,325 unique peptides, which confirmed 10,076 of the predicted genes. The addition of transcriptomic and proteogenomic approaches resulted in the identification of 53 novel proteins and 34 gene-correction events that were missed by genomic approaches. Proteogenomic analysis also helped in reassigning 1,348 potential noncoding RNAs as bona fide protein-coding messenger RNAs. Gene expression patterns at the RNA and protein levels indicated that peptide sequencing was useful in capturing proteins encoded by nuclear and organellar genomes alike. Mass spectrometry-based proteomic evidence provided an unbiased approach toward the identification of proteins encoded by organellar genomes. Such proteins are often missed in transcriptome data sets due to the enrichment of only messenger RNAs that contain poly(A) tails. Overall, the use of integrated omic approaches enhanced the quality of the assembly and annotation of this nonmodel plant genome. The availability of genomic, transcriptomic, and proteomic data will enhance genomics-assisted breeding, germplasm characterization, and conservation of sandalwood trees.

Microbiome and ecotypic adaption of Holcus lanatus (L.) to extremes of its soil pH range, investigated through transcriptome sequencing

BACKGROUND

Plants can adapt to edaphic stress, such as nutrient deficiency, toxicity and biotic challenges, by controlled transcriptomic responses, including microbiome interactions. Traditionally studied in model plant species with controlled microbiota inoculation treatments, molecular plant-microbiome interactions can be functionally investigated via RNA-Seq. Complex, natural plant-microbiome studies are limited, typically focusing on microbial rRNA and omitting functional microbiome investigations, presenting a fundamental knowledge gap. Here, root and shoot meta-transcriptome analyses, in tandem with shoot elemental content and root staining, were employed to investigate transcriptome responses in the wild grass Holcus lanatus and its associated natural multi-species eukaryotic microbiome. A full factorial reciprocal soil transplant experiment was employed, using plant ecotypes from two widely contrasting natural habitats, acid bog and limestone quarry soil, to investigate naturally occurring, and ecologically meaningful, edaphically driven molecular plant-microbiome interactions.

RESULTS

Arbuscular mycorrhizal (AM) and non-AM fungal colonization was detected in roots in both soils. Staining showed greater levels of non-AM fungi, and transcriptomics indicated a predominance of Ascomycota-annotated genes. Roots in acid bog soil were dominated by Phialocephala-annotated transcripts, a putative growth-promoting endophyte, potentially involved in N nutrition and ion homeostasis. Limestone roots in acid bog soil had greater expression of other Ascomycete genera and Oomycetes and lower expression of Phialocephala-annotated transcripts compared to acid ecotype roots, which corresponded with reduced induction of pathogen defense processes, particularly lignin biosynthesis in limestone ecotypes. Ascomycota dominated in shoots and limestone soil roots, but Phialocephala-annotated transcripts were insignificant, and no single Ascomycete genus dominated. Fusarium-annotated transcripts were the most common genus in shoots, with Colletotrichum and Rhizophagus (AM fungi) most numerous in limestone soil roots. The latter coincided with upregulation of plant genes involved in AM symbiosis initiation and AM-based P acquisition in an environment where P availability is low.

CONCLUSIONS

Meta-transcriptome analyses provided novel insights into H. lanatus transcriptome responses, associated eukaryotic microbiota functions and taxonomic community composition. Significant edaphic and plant ecotype effects were identified, demonstrating that meta-transcriptome-based functional analysis is a powerful tool for the study of natural plant-microbiome interactions.

Diversification of defensins and NLRs in Arabidopsis species by different evolutionary mechanisms

BACKGROUND

Genes encoding proteins underlying host-pathogen co-evolution and which are selected for new resistance specificities frequently are under positive selection, a process that maintains diversity. Here, we tested the contribution of natural selection, recombination and transcriptional divergence to the evolutionary diversification of the plant defensins superfamily in three Arabidopsis species. The intracellular NOD-like receptor (NLR) family was used for comparison because positive selection has been well documented in its members. Similar to defensins, NLRs are encoded by a large and polymorphic gene family and many of their members are involved in the immune response.

RESULTS

Gene trees of Arabidopsis defensins (DEFLs) show a high prevalence of clades containing orthologs. This indicates that their diversity dates back to a common ancestor and species-specific duplications did not significantly contribute to gene family expansion. DEFLs are characterized by a pervasive pattern of neutral evolution with infrequent positive and negative selection as well as recombination. In comparison, most NLR alignment groups are characterized by frequent occurrence of positive selection and recombination in their leucine-rich repeat (LRR) domain as well negative selection in their nucleotide-binding (NB-ARC) domain. While major NLR subgroups are expressed in pistils and leaves both in presence or absence of pathogen infection, the members of DEFL alignment groups are predominantly transcribed in pistils. Furthermore, conserved groups of NLRs and DEFLs are differentially expressed in response to Fusarium graminearum regardless of whether these genes are under positive selection or not.

CONCLUSIONS

The present analyses of NLRs expands previous studies in Arabidopsis thaliana and highlights contrasting patterns of purifying and diversifying selection affecting different gene regions. DEFL genes show a different evolutionary trend, with fewer recombination events and significantly fewer instances of natural selection. Their heterogeneous expression pattern suggests that transcriptional divergence probably made the major contribution to functional diversification. In comparison to smaller families encoding pathogenesis-related (PR) proteins under positive selection, DEFLs are involved in a wide variety of processes that altogether might pose structural and functional trade-offs to their family-wide pattern of evolution.

Thermal acclimation in Arabidopsis lyrata: genotypic costs and transcriptional changes

Frost and heat events can be challenging for sessile organisms that cannot escape thermal extremes. However, adverse effects of thermal stress on fitness may be reduced by pre-exposure to cold or heat, a process known as acclimation. To understand the ecological and evolutionary implications of acclimation, we investigated (1) the reduction in performance due to stress pre-exposure, (2) the magnitude of increased leaf resistance to subsequent stress, (3) the costs of acclimation and (4) the genes differing in expression due to stress pre-exposure. Plants of Arabidopsis lyrata were raised under three treatments of pre-exposure: bouts of frost, bouts of heat or constant temperature. Resistance of leaves to subsequent frost and heat stress was then measured by electrolyte leakage. RNA-seq analysis was performed to examine the genes differentially expressed between stress-pre-exposed and control plants. Pre-exposure to stress during growth decreased plant size and increased leaf resistance to subsequent stress independent of whether pre-exposure was to frost or heat. But the highest increase in leaf resistance to frost was found after pre-exposure to frost (as a trend) and in leaf resistance to heat after pre-exposure to heat. No evidence for costs of acclimation was detected. RNA-sequencing suggested that acclimation by frost and heat pre-exposure was caused by distinct mechanisms: modification of the chloroplast membrane and modification of the cell wall and membrane, respectively. Our results suggest that thermal resistance is a labile complex of traits, strongly affected by the previously experienced stress environment, with undetermined costs.

Pinpointing genes underlying annual/perennial transitions with comparative genomics

BACKGROUND

Transitions between perennial and an annual life history occur often in plant lineages, but the genes that control whether a plant is an annual or perennial are largely unknown. To identify genes that confer differences between annuals and perennials we compared the gene content of four pairs of sister lineages (Arabidopsis thaliana/Arabidopsis lyrata, Arabis montbretiana/Arabis alpina, Arabis verna/Aubrieta parviflora and Draba nemorosa/Draba hispanica) in the Brassicaceae in which each pair contains one annual and one perennial, plus one extra annual species (Capsella rubella).

RESULTS

After sorting all genes in all nine species into gene families, we identified five families in which well-annotated genes are present in the perennials A. lyrata and A. alpina, but are not present in any of the annual species. For the eleven genes in perennials in these families, an orthologous pseudogene or otherwise highly diverged gene was found in the syntenic region of the annual species in six cases. The five candidate families identified encode: a kinase, an oxidoreductase, a lactoylglutathione lyase, a F-box protein and a zinc finger protein. By comparing the active gene in the perennial to the pseudogene or heavily altered gene in the annual, dN and dS were calculated. The low dN/dS values in one kinase suggest that it became pseudogenized more recently, while the other kinase, F-box, oxidoreductase and zinc-finger became pseudogenized closer to the divergence between the annual-perennial pair.

CONCLUSIONS

We identified five gene families that may be involved in the life history switch from perennial to annual. Considering the dN and dS data and whether syntenic pseudogenes were found and the potential functions of the genes, the F-box family is considered the most promising candidate for future functional studies to determine if it affects life history.

Recurrent evolution of heat-responsiveness in Brassicaceae COPIA elements

BACKGROUND

The mobilization of transposable elements (TEs) is suppressed by host genome defense mechanisms. Recent studies showed that the cis-regulatory region of Arabidopsis thaliana COPIA78/ONSEN retrotransposons contains heat-responsive elements (HREs), which cause their activation during heat stress. However, it remains unknown whether this is a common and potentially conserved trait and how it has evolved.

RESULTS

We show that ONSEN, COPIA37, TERESTRA, and ROMANIAT5 are the major families of heat-responsive TEs in A. lyrata and A. thaliana. Heat-responsiveness of COPIA families is correlated with the presence of putative high affinity heat shock factor binding HREs within their long terminal repeats in seven Brassicaceae species. The strong HRE of ONSEN is conserved over millions of years and has evolved by duplication of a proto-HRE sequence, which was already present early in the evolution of the Brassicaceae. However, HREs of most families are species-specific, and in Boechera stricta, the ONSEN HRE accumulated mutations and lost heat-responsiveness.

CONCLUSIONS

Gain of HREs does not always provide an ultimate selective advantage for TEs, but may increase the probability of their long-term survival during the co-evolution of hosts and genomic parasites.

Borrowed alleles and convergence in serpentine adaptation

Serpentine barrens represent extreme hazards for plant colonists. These sites are characterized by high porosity leading to drought, lack of essential mineral nutrients, and phytotoxic levels of metals. Nevertheless, nature forged populations adapted to these challenges. Here, we use a population-based evolutionary genomic approach coupled with elemental profiling to assess how autotetraploid Arabidopsis arenosa adapted to a multichallenge serpentine habitat in the Austrian Alps. We first demonstrate that serpentine-adapted plants exhibit dramatically altered elemental accumulation levels in common conditions, and then resequence 24 autotetraploid individuals from three populations to perform a genome scan. We find evidence for highly localized selective sweeps that point to a polygenic, multitrait basis for serpentine adaptation. Comparing our results to a previous study of independent serpentine colonizations in the closely related diploid Arabidopsis lyrata in the United Kingdom and United States, we find the highest levels of differentiation in 11 of the same loci, providing candidate alleles for mediating convergent evolution. This overlap between independent colonizations in different species suggests that a limited number of evolutionary strategies are suited to overcome the multiple challenges of serpentine adaptation. Interestingly, we detect footprints of selection in A. arenosa in the context of substantial gene flow from nearby off-serpentine populations of A. arenosa, as well as from A. lyrata In several cases, quantitative tests of introgression indicate that some alleles exhibiting strong selective sweep signatures appear to have been introgressed from A. lyrata This finding suggests that migrant alleles may have facilitated adaptation of A. arenosa to this multihazard environment.