Genome structures and transcriptomes signify niche adaptation for the multiple-ion-tolerant extremophyte Schrenkiella parvula

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

Balancing growth amidst salt stress - lifestyle perspectives from the extremophyte model Schrenkiella parvula

Schrenkiella parvula, a leading extremophyte model in Brassicaceae, can grow and complete its lifecycle under multiple environmental stresses, including high salinity. Yet, the key physiological and structural traits underlying its stress-adapted lifestyle are unknown along with trade-offs when surviving salt stress at the expense of growth and reproduction. We aimed to identify the influential adaptive trait responses that lead to stress-resilient and uncompromised growth across developmental stages when treated with salt at levels known to inhibit growth in Arabidopsis and most crops. Its resilient growth was promoted by traits that synergistically allowed primary root growth in seedlings, the expansion of xylem vessels across the root-shoot continuum, and a high capacity to maintain tissue water levels by developing thicker succulent leaves while enabling photosynthesis during salt stress. A successful transition from vegetative to reproductive phase was initiated by salt-induced early flowering, resulting in viable seeds. Self-fertilization in salt-induced early flowering was dependent upon filament elongation in flowers otherwise aborted in the absence of salt during comparable plant ages. The maintenance of leaf water status promoting growth, and early flowering to ensure reproductive success in a changing environment, were among the most influential traits that contributed to the extremophytic lifestyle of S. parvula.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Effective root responses to salinity stress include maintained cell expansion and carbon allocation

Acclimation of root growth is vital for plants to survive salt stress. Halophytes are great examples of plants that thrive even under severe salinity, but their salt tolerance mechanisms, especially those mediated by root responses, are still largely unknown. We compared root growth responses of the halophyte Schrenkiella parvula with its glycophytic relative species Arabidopsis thaliana under salt stress and performed transcriptomic analysis of S. parvula roots to identify possible gene regulatory networks underlying their physiological responses. Schrenkiella parvula roots do not avoid salt and experience less growth inhibition under salt stress. Salt-induced abscisic acid levels were higher in S. parvula roots compared with Arabidopsis. Root transcriptomic analysis of S. parvula revealed the induction of sugar transporters and genes regulating cell expansion and suberization under salt stress. ¹⁴ C-labeled carbon partitioning analyses showed that S. parvula continued allocating carbon to roots from shoots under salt stress while carbon barely allocated to Arabidopsis roots. Further physiological investigation revealed that S. parvula roots maintained root cell expansion and enhanced suberization under severe salt stress. In summary, roots of S. parvula deploy multiple physiological and developmental adjustments under salt stress to maintain growth, providing new avenues to improve salt tolerance of plants using root-specific strategies.

Living with high potassium: Balance between nutrient acquisition and K-induced salt stress signaling

High potassium (K) in the growth medium induces salinity stress in plants. However, the molecular mechanisms underlying plant responses to K-induced salt stress are virtually unknown. We examined Arabidopsis (Arabidopsis thaliana) and its extremophyte relative Schrenkiella parvula using a comparative multiomics approach to identify cellular processes affected by excess K and understand which deterministic regulatory pathways are active to avoid tissue damages while sustaining growth. Arabidopsis showed limited capacity to curb excess K accumulation and prevent nutrient depletion, contrasting to S. parvula which could limit excess K accumulation without restricting nutrient uptake. A targeted transcriptomic response in S. parvula promoted nitrogen uptake along with other key nutrients followed by uninterrupted N assimilation into primary metabolites during excess K-stress. This resulted in larger antioxidant and osmolyte pools and corresponded with sustained growth in S. parvula. Antithetically, Arabidopsis showed increased reactive oxygen species levels, reduced photosynthesis, and transcriptional responses indicative of a poor balance between stress signaling, subsequently leading to growth limitations. Our results indicate that the ability to regulate independent nutrient uptake and a coordinated transcriptomic response to avoid nonspecific stress signaling are two main deterministic steps toward building stress resilience to excess K+-induced salt stress.

Genome-wide analysis of annexin gene family in Schrenkiella parvula and Eutrema salsugineum suggests their roles in salt stress response

Annexins (Anns) play an important role in plant development, growth and responses to various stresses. Although Ann genes have been characterized in some plants, their role in adaptation mechanisms and tolerance to environmental stresses have not been studied in extremophile plants. In this study, Ann genes in Schrenkiella parvula and Eutrema salsugineum were identified using a genome-wide method and phylogenetic relationships, subcellular distribution, gene structures, conserved residues and motifs and also promoter prediction have been studied through bioinformatics analysis. We identified ten and eight encoding putative Ann genes in S. parvula and E. salsugineum genome respectively, which were divided into six subfamilies according to phylogenetic relationships. By observing conservation in gene structures and protein motifs we found that the majority of Ann members in two extremophile plants are similar. Furthermore, promoter analysis revealed a greater number of GATA, Dof, bHLH and NAC transcription factor binding sites, as well as ABRE, ABRE3a, ABRE4, MYB and Myc cis-acting elements in compare to Arabidopsis thaliana. To gain additional insight into the putative roles of candidate Ann genes, the expression of SpAnn1, SpAnn2 and SpAnn6 in S. parvula was studied in response to salt stress, which indicated that their expression level in shoot increased. Similarly, salt stress induced expression of EsAnn1, 5 and 7, in roots and EsAnn1, 2 and 5 in leaves of E. salsugineum. Our comparative analysis implies that both halophytes have different regulatory mechanisms compared to A. thaliana and suggest SpAnn2 gene play important roles in mediating salt stress.

Where whole-genome duplication is most beneficial: Adaptation of mangroves to a wide salinity range between land and sea

Whole-genome duplication (WGD) is believed to increase the chance of adaptation to a new environment. This conjecture may apply particularly well to new environments that are not only different but also more variable than ancestral habitats. One such prominent environment is the interface between land and sea, which has been invaded by woody plants, collectively referred as mangroves, multiple times. Here, we use two distantly related mangrove species (Avicennia marina and Rhizophora apiculata) to explore the effects of WGD on the adaptive process. We found that a high proportion of duplicated genes retained after WGD have acquired derived differential expression in response to salt gradient treatment. The WGD duplicates differentially expressed in at least one copy usually (>90%) diverge from their paralogues' expression profiles. Furthermore, both species evolved in parallel to have one paralogue expressed at a high level in both fresh water and hypersaline conditions but at a lower level at medium salinity. The pattern contrasts with the conventional view of monotone increase/decrease as salinity increases. Differentially expressed copies have thus probably acquired a new role in salinity tolerance. Our results indicate that the WGD duplicates may have evolved to function collaboratively in coping with different salinity levels, rather than specializing in the intermediate salinity optimal for mangrove plants. In conclusion, WGD and the retained duplicates appear to be an effective solution for adaptation to new and unstable environments.

Positive selection and heat-response transcriptomes reveal adaptive features of the Brassicaceae desert model, Anastatica hierochuntica

Plant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts. We generated an A. hierochuntica reference transcriptome and identified extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat, and detecting positively selected genes in A. hierochuntica. The two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat 'stress-ready' state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype. We suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

Divergence in the ABA gene regulatory network underlies differential growth control

The phytohormone abscisic acid (ABA) is a central regulator of acclimation to environmental stress; however, its contribution to differences in stress tolerance between species is unclear. To establish a comparative framework for understanding how stress hormone signalling pathways diverge across species, we studied the growth response of four Brassicaceae species to ABA treatment and generated transcriptomic and DNA affinity purification and sequencing datasets to construct a cross-species gene regulatory network (GRN) for ABA. Comparison of genes bound directly by ABA-responsive element binding factors suggests that cis-factors are most important for determining the target loci represented in the ABA GRN of a particular species. Using this GRN, we reveal how rewiring of growth hormone subnetworks contributes to stark differences in the response to ABA in the extremophyte Schrenkiella parvula. Our study provides a model for understanding how divergence in gene regulation can lead to species-specific physiological outcomes in response to hormonal cues.

Identification and molecular characterization of the alternative spliced variants of beta carbonic anhydrase 1 (βCA1) from Arabidopsis thaliana

Carbonic anhydrases (CAs) are ubiquitous zinc metalloenzymes that catalyze the interconversion of carbon dioxide and bicarbonate. Higher plants mainly contain the three evolutionarily distinct CA families αCA, βCA, and γCA, with each represented by multiple isoforms. Alternative splicing (AS) of the CA transcripts is common. However, there is little information on the spliced variants of individual CA isoforms. In this study, we focused on the characterization of spliced variants of βCA1 from Arabidopsis. The expression patterns and subcellular localization of the individual spliced variants of βCA1 were examined. The results showed that the spliced variants of βCA1 possessed different subcellular and tissue distributions and responded differently to environmental stimuli. Additionally, we addressed the physiological role of βCA1 in heat stress response and its protein-protein interaction (PPI) network. Our results showed that βCA1 was regulated by heat stresses, and βca1 mutant was hypersensitive to heat stress, indicating a role for βCA1 in heat stress response. Furthermore, PPI network analysis revealed that βCA1 interacts with multiple proteins involved in several processes, including photosynthesis, metabolism, and the stress response, and these will provide new avenues for future investigations of βCA1.

Disruption of AtHAK/KT/KUP9 enhances plant cesium accumulation under low potassium supply

Understanding the molecular mechanisms that underlie cesium (Cs⁺ ) transport in plants is important to limit the entry of its radioisotopes from contaminated areas into the food chain. The potentially toxic element Cs⁺ , which is not involved in any biological process, is chemically closed to the macronutrient potassium (K⁺ ). Among the multiple K⁺ carriers, the high-affinity K⁺ transporters family HAK/KT/KUP is thought to be relevant in mediating opportunistic Cs⁺ transport. Of the 13 KUP identified in A. thaliana, only HAK5, the major contributor to root K⁺ acquisition under low K⁺ supply, has been functionally demonstrated to be involved in Cs⁺ uptake in planta. In the present study, we showed that accumulation of Cs⁺ increased by up to 30% in two A. thaliana mutant lines lacking KUP9 and grown under low K⁺ supply. Since further experiments revealed that Cs⁺ release from contaminated plants to the external medium is proportionally lower in the two kup9 mutant alleles, we proposed that KUP9 disruption could impair Cs⁺ efflux. By contrast, K⁺ status in kup9 mutants is not affected, suggesting that KUP9 disruption does not alter substantially K⁺ transport in experimental conditions used. The putative primary role of KUP9 in plants is further discussed.

Cross species multi-omics reveals cell wall sequestration and elevated global transcript abundance as mechanisms of boron tolerance in plants

Boron toxicity is a world-wide problem for crops, yet we have a limited understanding of the genetic responses and adaptive mechanisms to this stress in plants. We employed a cross-species comparison between boron stress-sensitive Arabidopsis thaliana and its boron stress-tolerant extremophyte relative Schrenkiella parvula, and a multi-omics approach integrating genomics, transcriptomics, metabolomics and ionomics to assess plant responses and adaptations to boron stress. Schrenkiella parvula maintains lower concentrations of total boron and free boric acid than Arabidopsis when grown with excess boron. Schrenkiella parvula excludes excess boron more efficiently than Arabidopsis, which we propose is partly driven by SpBOR5, a boron transporter that we functionally characterize in this study. Both species use cell walls as a partial sink for excess boron. When accumulated in the cytoplasm, excess boron appears to interrupt RNA metabolism. The extremophyte S. parvula facilitates critical cellular processes while maintaining the pool of ribose-containing compounds that can bind with boric acid. The S. parvula transcriptome is pre-adapted to boron toxicity. It exhibits substantial overlaps with the Arabidopsis boron-stress responsive transcriptome. Cell wall sequestration and increases in global transcript levels under excess boron conditions emerge as key mechanisms for sustaining plant growth under boron toxicity.

A comparison of shared patterns of differential gene expression and gene ontologies in response to water-stress in roots and leaves of four diverse genotypes of Lolium and Festuca spp. temperate pasture grasses

Ryegrasses (Lolium spp.) and fescues (Festuca spp.) are closely related and widely cultivated perennial forage grasses. As such, resilience in the face of abiotic stresses is an important component of their traits. We have compared patterns of differentially expressed genes (DEGs) in roots and leaves of two perennial ryegrass genotypes and a single genotype of each of a festulolium (predominantly Italian ryegrass) and meadow fescue with the onset of water stress, focussing on overall patterns of DEGs and gene ontology terms (GOs) shared by all four genotypes. Plants were established in a growing medium of vermiculite watered with nutrient solution. Leaf and root material were sampled at 35% (saturation) and, as the medium dried, at 15%, 5% and 1% estimated water contents (EWCs) and RNA extracted. Differential gene expression was evaluated comparing the EWC sampling points from RNAseq data using a combination of analysis methods. For all genotypes, the greatest numbers of DEGs were identified in the 35/1 and 5/1 comparisons in both leaves and roots. In total, 566 leaf and 643 root DEGs were common to all 4 genotypes, though a third of these leaf DEGs were not regulated in the same up/down direction in all 4 genotypes. For roots, the equivalent figure was 1% of the DEGs. GO terms shared by all four genotypes were often enriched by both up- and down-regulated DEGs in the leaf, whereas generally, only by either up- or down-regulated DEGs in the root. Overall, up-regulated leaf DEGs tended to be more genotype-specific than down-regulated leaf DEGs or root DEGs and were also associated with fewer GOs. On average, only 5–15% of the DEGs enriching common GO terms were shared by all 4 genotypes, suggesting considerable variation in DEGs between related genotypes in enacting similar biological processes.

HKT sodium and potassium transporters in Arabidopsis thaliana and related halophyte species

High salinity induces osmotic stress and often leads to sodium ion-specific toxicity, with inhibitory effects on physiological, biochemical and developmental pathways. To cope with increased Na⁺ in soil water, plants restrict influx, compartmentalize ions into vacuoles, export excess Na⁺ from the cell, and distribute ions between the aerial and root organs. In this review, we discuss our current understanding of how high-affinity K⁺ transporters (HKT) contribute to salinity tolerance, focusing on HKT1-like family members primarily involved in long-distance transport, and in the recent research in the model plant Arabidopsis and its halophytic counterparts of the Eutrema genus. Functional characterization of the salt overly sensitive (SOS) pathway and HKT1-type transporters in these species indicate that they utilize similar approaches to deal with salinity, regardless of their tolerance.

Gene duplication and stress genomics in Brassicas: Current understanding and future prospects

Polyploidy or whole genome duplication (WGD) is an evolutionary phenomenon that happened in all angiosperms multiple times over millions of years. Extensive studies on the model plant Arabidopsis thaliana genome have revealed that it has undergone five rounds of WGDs followed, in the Brassicaceae tribe, by a characteristic whole genome triplication (WGT). In addition, small-scale events such as tandem or segmental duplications and retrotransposition also enable plants to reshape their genomes. Over the decades, extensive research efforts have been undertaken to understand the evolutionary significance of polyploidy. On the other hand, much less attention has been paid to understanding the impact of gene duplication on the diversification of important stress response genes. The main objective of this review is to discuss key aspects of gene and genome duplications with a focus on genes primarily regulated by osmotic stresses. The focal family is the Brassicaceae, since it (i) underwent multiple rounds of WGDs plus WGTs, (ii) hosts many economically important crops and wild relatives that are tolerant to a range of stresses, and (iii) comprises many species that have already been sequenced. Diverse molecular mechanisms that lead to structural and regulatory alterations of duplicated genes are discussed. Examples are drawn from recent literature to elucidate expanded, stress responsive gene families identified from different Brassica crops. A combined bioinformatic and transcriptomic method has been proposed and tested on a known stress-responsive gene pair to prove that stress-responsive duplicated allelic variants can be identified by this method. Finally, future prospects for engineering these genes into crops to enhance stress tolerance are discussed, and important resources for Brassica genome research are provided.

Genomic Convergence in the Adaptation to Extreme Environments

Convergent evolution is especially common in plants that have independently adapted to the same extreme environments (i.e., extremophile plants). The recent burst of omics data has alleviated many limitations that have hampered molecular convergence studies of non-model extremophile plants. In this review, we summarize cases of genomic convergence in these taxa to examine the extent and type of genomic convergence during the process of adaptation to extreme environments. Despite being well studied by candidate gene approaches, convergent evolution at individual sites is rare and often has a high false-positive rate when assessed in whole genomes. By contrast, genomic convergence at higher genetic levels has been detected during adaptation to the same extreme environments. Examples include the convergence of biological pathways and changes in gene expression, gene copy number, amino acid usage, and GC content. Higher convergence levels play important roles in the adaptive evolution of extremophiles and may be more frequent and involve more genes. In several cases, multiple types of convergence events have been found to co-occur. However, empirical and theoretical studies of this higher level convergent evolution are still limited. In conclusion, both the development of powerful approaches and the detection of convergence at various genetic levels are needed to further reveal the genetic mechanisms of plant adaptation to extreme environments.

Duplication and functional divergence of a calcium sensor in the Brassicaceae

The presence of varied numbers of CALCINEURIN B-LIKE10 (CBL10) calcium sensor genes in species across the Brassicaceae and the demonstrated role of CBL10 in salt tolerance in Arabidopsis thaliana and Eutrema salsugineum provided a unique opportunity to determine if CBL10 function is modified in different species and linked to salt tolerance. Salinity effects on species growth and cross-species complementation were used to determine the extent of conservation and divergence of CBL10 function in four species representing major lineages within the core Brassicaceae (A. thaliana, E. salsugineum, Schrenkiella parvula, and Sisymbrium irio) as well as the first diverging lineage (Aethionema arabicum). Evolutionary and functional analyses indicate that CBL10 duplicated within expanded lineage II of the Brassicaceae and that, while portions of CBL10 function are conserved across the family, there are species-specific variations in CBL10 function. Paralogous CBL10 genes within a species diverged in expression and function probably contributing to the maintenance of the duplicated gene pairs. Orthologous CBL10 genes diverged in function in a species-specific manner, suggesting that functions arose post-speciation. Multiple CBL10 genes and their functional divergence may have expanded calcium-mediated signaling responses and contributed to the ability of certain members of the Brassicaceae to maintain growth in salt-affected soils.

Salt Tolerance Mechanisms of Plants

Crop loss due to soil salinization is an increasing threat to agriculture worldwide. This review provides an overview of cellular and physiological mechanisms in plant responses to salt. We place cellular responses in a time- and tissue-dependent context in order to link them to observed phases in growth rate that occur in response to stress. Recent advances in phenotyping can now functionally or genetically link cellular signaling responses, ion transport, water management, and gene expression to growth, development, and survival. Halophytes, which are naturally salt-tolerant plants, are highlighted as success stories to learn from. We emphasize that (a) filling the major knowledge gaps in salt-induced signaling pathways, (b) increasing the spatial and temporal resolution of our knowledge of salt stress responses, (c) discovering and considering crop-specific responses, and (d) including halophytes in our comparative studies are all essential in order to take our approaches to increasing crop yields in saline soils to the next level.

Occurrence and cycling of trace elements in ultramafic soils and their impacts on human health: A critical review

The transformation of trace metals (TMs) in natural environmental systems has created significant concerns in recent decades. Ultramafic environments lead to potential risks to the agricultural products and, subsequently, to human health. This unique review presents geochemistry of ultramafic soils, TM fractionation (i.e. sequential and single extraction techniques), TM uptake and accumulation mechanisms of ultramafic flora, and ultramafic-associated health risks to human and agricultural crops. Ultramafic soils contain high levels of TMs (i.e. Cr, Ni, Mn, and Co) and have a low Ca:Mg ratio together with deficiencies in essential macronutrients required for the growth of crops. Even though a higher portion of TMs bind with the residual fraction of ultramafic soils, environmental changes (i.e. natural or anthropogenic) may increase the levels of TMs in the bioavailable or extractable fractions of ultramafic soils. Extremophile plants that have evolved to thrive in ultramafic soils present clear examples of evolutionary adaptations to TM resistance. The release of TMs into water sources and accumulation in food crops in and around ultramafic localities increases health risks for humans. Therefore, more focused investigations need to be implemented to understand the mechanisms related to the mobility and bioavailability of TMs in different ultramafic environments. Research gaps and directions for future studies are also discussed in this review. Lastly, we consider the importance of characterizing terrestrial ultramafic soil and its effect on crop plants in the context of multi-decadal plans by NASA and other space agencies to establish human colonies on Mars.

Landscape of gene transposition-duplication within the Brassicaceae family

We developed the CLfinder-OrthNet pipeline that detects co-linearity among multiple closely related genomes, finds orthologous gene groups, and encodes the evolutionary history of each orthologue group into a representative network (OrthNet). Using a search based on network topology, we identified 1,394 OrthNets that included gene transposition-duplication (tr-d) events, out of 17,432 identified in six Brassicaceae genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. These were derived from older tr-d events and had a higher chance of being expressed. We also found an enrichment of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3 that relate to their adaptive evolution. CLfinder-OrthNet provides a flexible toolkit to compare gene order, visualize evolutionary paths among orthologues as networks, and identify gene loci that share an evolutionary history.

Role and Functional Differences of HKT1-Type Transporters in Plants under Salt Stress

Abiotic stresses generally cause a series of morphological, biochemical and molecular changes that unfavorably affect plant growth and productivity. Among these stresses, soil salinity is a major threat that can seriously impair crop yield. To cope with the effects of high salinity on plants, it is important to understand the mechanisms that plants use to deal with it, including those activated in response to disturbed Na⁺ and K⁺ homeostasis at cellular and molecular levels. HKT1-type transporters are key determinants of Na⁺ and K⁺ homeostasis under salt stress and they contribute to reduce Na⁺-specific toxicity in plants. In this review, we provide a brief overview of the function of HKT1-type transporters and their importance in different plant species under salt stress. Comparison between HKT1 homologs in different plant species will shed light on different approaches plants may use to cope with salinity.

Salt tolerance mechanisms in three Irano-Turanian Brassicaceae halophytes relatives of Arabidopsis thaliana

Salt tolerance mechanisms were studied in three Irano-Turanian halophytic species from the Brassicaceae ‎‎(Lepidium latifolium, L. perfoliatum and Schrenkiella parvula) and compared with the glycophyte Arabidopsis thaliana. According to seed germination under salt stress, L. perfoliatum was the most tolerant species, while L. latifolium and S. parvula were rather susceptible. Contrastingly, based on biomass production L. perfoliatum was more salt sensitive than the other two species. In S. parvula biomass was increased up to 2.8-fold by 100 mM NaCl; no significant growth reduction was observed even when exposed to 400 mM NaCl. Stable activities of antioxidative defense enzymes, nil or negligible accumulation of superoxide anion and hydrogen peroxide, as well as stable membrane integrity in the three halophytes revealed that no oxidative stress occurred in these tolerant species under salt stress. Proline levels increased in response to salt treatment. However, it contributed only by 0.3‒2.0% to the total osmolyte concentration in the three halophytes (at 400 mM NaCl) and even less (0.04%) in the glycophyte, A. thaliana (at 100 mM NaCl). Soluble sugars in all three halophytes and free amino acids pool in S. parvula decreased under salt treatment in contrast to the glycophyte, A. thaliana. The contribution of organic osmolytes to the total osmolyte pool increased by salt treatment in the roots, while decreased in halophyte and glycophyte, A. thaliana leaves. Interestingly, this reduction was compensated by a higher relative contribution of K in the leaves of the halophytes, but of Na in A. thaliana. Taken together, biomass data and biochemical indicators show that S. parvula is more salt tolerant than the two Lepidium species. Our data indicate that L. latifolium, as a perennial halophyte with a large biomass, is highly suitable for both restoration of saline habitats and saline agriculture.

Halophytism: What Have We Learnt From Arabidopsis thaliana Relative Model Systems?

Halophytes are able to thrive in salt concentrations that would kill 99% of other plant species, and identifying their salt-adaptive mechanisms has great potential for improving the tolerance of crop plants to salinized soils. Much research has focused on the physiological basis of halophyte salt tolerance, whereas the elucidation of molecular mechanisms has traditionally lagged behind due to the absence of a model halophyte system. However, over the last decade and a half, two Arabidopsis (Arabidopsis thaliana) relatives, Eutrema salsugineum and Schrenkiella parvula, have been established as transformation-competent models with various genetic resources including high-quality genome assemblies. These models have facilitated powerful comparative analyses with salt-sensitive Arabidopsis to unravel the genetic adaptations that enable a halophytic lifestyle. The aim of this review is to explore what has been learned about halophytism using E. salsugineum and S. parvula We consider evidence from physiological and molecular studies suggesting that differences in salt tolerance between related halophytes and salt-sensitive plants are associated with alterations in the regulation of basic physiological, biochemical, and molecular processes. Furthermore, we discuss how salt tolerance mechanisms of the halophytic models are reflected at the level of their genomes, where evolutionary processes such as subfunctionalization and/or neofunctionalization have altered the expression and/or functions of genes to facilitate adaptation to saline conditions. Lastly, we summarize the many areas of research still to be addressed with E. salsugineum and S. parvula as well as obstacles hindering further progress in understanding halophytism.

Transcriptomic view of survival during early seedling growth of the extremophyte Haloxylon ammodendron

Seedling establishment in an extreme environment requires an integrated genomic and physiological response to survive multiple abiotic stresses. The extremophyte, Haloxylon ammodendron is a pioneer species capable of colonizing temperate desert sand dunes. We investigated the induced and basal transcriptomes in H. ammodendron under water-deficit stress during early seedling establishment. We find that not only drought-responsive genes, but multiple genes in pathways associated with salt, osmotic, cold, UV, and high-light stresses were induced, suggesting an altered regulatory stress response system. Additionally, H. ammodendron exhibited enhanced biotic stress tolerance by down-regulation of genes that were generally up-regulated during pathogen entry in susceptible plants. By comparing the H. ammodendron basal transcriptome to six closely related transcriptomes in Amaranthaceae, we detected enriched basal level transcripts in H. ammodendron that shows preadaptation to abiotic stress and pathogens. We found transcripts that were generally maintained at low levels and some induced only under abiotic stress in the stress-sensitive model, Arabidopsis thaliana to be highly expressed under basal conditions in the Amaranthaceae transcriptomes including H. ammodendron. H. ammodendron shows coordinated expression of genes that regulate stress tolerance and seedling development resource allocation to support survival against multiple stresses in a sand dune dominated temperate desert environment.

The High-Affinity Potassium Transporter EpHKT1;2 From the Extremophile Eutrema parvula Mediates Salt Tolerance

To survive salt stress, plants must maintain a balance between sodium and potassium ions. High-affinity potassium transporters (HKTs) play a key role in reducing Na+ toxicity through K+ uptake. Eutrema parvula (formerly known as Thellungiella parvula), a halophyte closely related to Arabidopsis, has two HKT1 genes that encode EpHKT1;1 and EpHKT1;2. In response to high salinity, the EpHKT1;2 transcript level increased rapidly; by contrast, the EpHKT1;1 transcript increased more slowly in response to salt treatment. Yeast cells expressing EpHKT1;2 were able to tolerate high concentrations of NaCl, whereas EpHKT1;1-expressing yeast cells remained sensitive to NaCl. Amino acid sequence alignment with other plant HKTs showed that EpHKT1;1 contains an asparagine residue (Asn-213) in the second pore-loop domain, but EpHKT1;2 contains an aspartic acid residue (Asp-205) at the same position. Yeast cells expressing EpHKT1;1, in which Asn-213 was substituted with Asp, were able to tolerate high concentrations of NaCl. In contrast, substitution of Asp-205 by Asn in EpHKT1;2 did not enhance salt tolerance and rather resulted in a similar function to that of AtHKT1 (Na+ influx but no K+ influx), indicating that the presence of Asn or Asp determines the mode of cation selectivity of the HKT1-type transporters. Moreover, Arabidopsis plants (Col-gl) overexpressing EpHKT1;2 showed significantly higher tolerance to salt stress and accumulated less Na+ and more K+ compared to those overexpressing EpHKT1;1 or AtHKT1. Taken together, these results suggest that EpHKT1;2 mediates tolerance to Na+ ion toxicity in E. parvula and is a major contributor to its halophytic nature.

The Tropical Invasive Seagrass, Halophila stipulacea, Has a Superior Ability to Tolerate Dynamic Changes in Salinity Levels Compared to Its Freshwater Relative, Vallisneria americana

The tropical seagrass species, Halophila stipulacea, originated from the Indian Ocean and the Red Sea, subsequently invading the Mediterranean and has recently established itself in the Caribbean Sea. Due to its invasive nature, there is growing interest in understanding this species' capacity to adapt to new conditions. One approach to understanding the natural tolerance of a plant is to compare the tolerant species with a closely related non-tolerant species. We compared the physiological responses of H. stipulacea exposed to different salinities, with that of its nearest freshwater relative, Vallisneria americana. To achieve this goal, H. stipulacea and V. americana plants were grown in dedicated microcosms, and exposed to the following salt regimes: (i) H. stipulacea: control (40 PSU, practical salinity units), hyposalinity (25 PSU) and hypersalinity (60 PSU) for 3 weeks followed by a 4-week recovery phase (back to 40 PSU); (ii) V. americana: control (1 PSU), and hypersalinity (12 PSU) for 3 weeks, followed by a 4-week recovery phase (back to 1 PSU). In H. stipulacea, leaf number and chlorophyll content showed no significant differences between control plants and plants under hypo and hypersalinities, but a significant decrease in leaf area under hypersalinity was observed. In addition, compared with control plants, H. stipulacea plants exposed to hypo and hypersalinity were found to have reduced below-ground biomass and C/N ratios, suggesting changes in the allocation of resources in response to both stresses. There was no significant effect of hypo/hypersalinity on dark-adapted quantum yield of photosystem II (Fv/Fm) suggesting that H. stipulacea photochemistry is resilient to hypo/hypersalinity stress. In contrast to the seagrass, V. americana exposed to hypersalinity displayed significant decreases in above-ground biomass, shoot number, leaf number, blade length and Fv/Fm, followed by significant recoveries of all these parameters upon return of the plants to non-saline control conditions. These data suggest that H. stipulacea shows remarkable tolerance to both hypo and hypersalinity. Resilience to a relatively wide range of salinities may be one of the traits explaining the invasive nature of this species in the Mediterranean and Caribbean Seas.

Q&A: How do gene regulatory networks control environmental responses in plants?

A gene regulatory network (GRN) describes the hierarchical relationship between transcription factors, associated proteins, and their target genes. Studying GRNs allows us to understand how a plant's genotype and environment are integrated to regulate downstream physiological responses. Current efforts in plants have focused on defining the GRNs that regulate functions such as development and stress response and have been performed primarily in genetically tractable model plant species such as Arabidopsis thaliana. Future studies will likely focus on how GRNs function in non-model plants and change over evolutionary time to allow for adaptation to extreme environments. This broader understanding will inform efforts to engineer GRNs to create tailored crop traits.

Conceptualizing plant systems evolution

Organisms inhabiting extreme environments are emerging models in systems evolution, enabling us to identify the molecular alterations effecting major phenotypic divergence through comparative approaches. Here I discuss possible physiological mechanisms underlying evolutionary adaptations to extreme environments both theoretically and in relation to experimental observations. Reasoning leads me to the 'conserved steady-state' hypothesis of evolutionary adaptation: Between closely related plants adapted to differently composed soils, the homeostatically controlled steady-state set point cytosolic (buffered) concentrations of mineral ions are conserved. Subsequently, I compare molecular alterations expected to contribute to physiological adaptations with our present knowledge. Key roles of enhanced gene product dosage in plant evolutionary adaptations question the widespread stimulus response-centric paradigm. As a broader implication, co-regulation networks can lack decisive functional network elements. With this article, I hope to stimulate a discussion across research fields and provide an initial conceptual framework for future experimental testing and for quantitative modelling.

Transcriptome Signatures of Selection, Drift, Introgression, and Gene Duplication in the Evolution of an Extremophile Endemic Plant

Plants on serpentine soils provide extreme examples of adaptation to environment, and thus offer excellent models for the study of evolution at the molecular and genomic level. Serpentine outcrops are derived from ultramafic rock and have extremely low levels of essential plant nutrients (e.g., N, P, K, and Ca), as well as toxic levels of heavy metals (e.g., Ni, Cr, and Co) and low moisture availability. These outcrops provide habitat to a number of endemic plant species, including the annual mustard Caulanthus amplexicaulis var. barbarae (Cab) (Brassicaceae). Its sister taxon, C. amplexicaulis var. amplexicaulis (Caa), is intolerant to serpentine soils. Here, we assembled and annotated comprehensive reference transcriptomes of both Caa and Cab for use in protein coding sequence comparisons. A set of 29,443 reciprocal best Blast hit (RBH) orthologs between Caa and Cab was compared with identify coding sequence variants, revealing a high genome-wide dN/dS ratio between the two taxa (mean = 0.346). We show that elevated dN/dS likely results from the composite effects of genetic drift, positive selection, and the relaxation of negative selection. Further, analysis of paralogs within each taxon revealed the signature of a period of elevated gene duplication (∼10 Ma) that is shared with other species of the tribe Thelypodieae, and may have played a role in the striking morphological and ecological diversity of this tribe. In addition, distribution of the synonymous substitution rate, dS, is strongly bimodal, indicating a history of reticulate evolution that may have contributed to serpentine adaptation.

A comparative transcriptomic analysis reveals the core genetic components of salt and osmotic stress responses in Braya humilis

Braya humilis is a member of the Euclidieae tribe within the family Brassicaceae. This species exhibits a broad range of adaptations to different climatic zones and latitudes as it has a distribution that ranges from northern Asia to the arctic-alpine regions of northern North America. In China, B. humilis is mainly found on the Qinghai-Tibetan Plateau (QTP) and in adjacent arid regions. In this study, we sequenced a sample from an arid region adjacent to the QTP using the Illumina platform generating a total of 46,485 highly accurate unigenes, of which 78.41% were annotated by BLASTing versus public protein databases. The B. humilis transcriptome is characterized by a high level of sequence conservation compared with its close relative, Arabidopsis thaliana. We also used reciprocal blast to identify shared orthologous genes between B. humilis and four other sequenced Brassicaceae species (i.e. A. thaliana, A. lyrata, Capsella rubella, and Thellungiella parvula). To enable precise characterization of orthologous genes, the early-diverging basal angiosperm Amborella trichopoda was also included. A total of 6,689 orthologous genes were identified before stricter criteria for the determination of e-values, amino acid hit lengths, and identity values was applied to further reduce this list. This led to a final list of 381 core orthologous genes for B. humilis; 39 out of these genes are involved in salt and osmotic stress responses and estimations of nonsynonymous/synonymous substitution ratios for this species and A. thaliana orthologs show that these genes are under purifying selection in B. humilis. Expression of six genes was detected in B. humilis seedlings under salt and osmotic stress treatments. Comparable expression patterns to their counterparts in Arabidopsis suggest that these orthologous genes are both sequence and functional conservation. The results of this study demonstrate that the environmental adaptations of B. humilis are mainly the results of preexisting genetic components. Future work will be required to characterize the expression patterns of these orthologous genes in natural populations and will provide further insights into the adaptive mechanisms underlying the wide range of B. humilis adaptations.

Anastatica hierochuntica, an Arabidopsis Desert Relative, Is Tolerant to Multiple Abiotic Stresses and Exhibits Species-Specific and Common Stress Tolerance Strategies with Its Halophytic Relative, Eutrema (Thellungiella) salsugineum

The search for novel stress tolerance determinants has led to increasing interest in plants native to extreme environments - so called "extremophytes." One successful strategy has been comparative studies between Arabidopsis thaliana and extremophyte Brassicaceae relatives such as the halophyte Eutrema salsugineum located in areas including cold, salty coastal regions of China. Here, we investigate stress tolerance in the desert species, Anastatica hierochuntica (True Rose of Jericho), a member of the poorly investigated lineage III Brassicaceae. We show that A. hierochuntica has a genome approximately 4.5-fold larger than Arabidopsis, divided into 22 diploid chromosomes, and demonstrate that A. hierochuntica exhibits tolerance to heat, low N and salt stresses that are characteristic of its habitat. Taking salt tolerance as a case study, we show that A. hierochuntica shares common salt tolerance mechanisms with E. salsugineum such as tight control of shoot Na+ accumulation and resilient photochemistry features. Furthermore, metabolic profiling of E. salsugineum and A. hierochuntica shoots demonstrates that the extremophytes exhibit both species-specific and common metabolic strategies to cope with salt stress including constitutive up-regulation (under control and salt stress conditions) of ascorbate and dehydroascorbate, two metabolites involved in ROS scavenging. Accordingly, A. hierochuntica displays tolerance to methyl viologen-induced oxidative stress suggesting that a highly active antioxidant system is essential to cope with multiple abiotic stresses. We suggest that A. hierochuntica presents an excellent extremophyte Arabidopsis relative model system for understanding plant survival in harsh desert conditions.

Plant Carbonic Anhydrases: Structures, Locations, Evolution, and Physiological Roles

Carbonic anhydrases (CAs) are zinc metalloenzymes that catalyze the interconversion of CO2 and HCO3- and are ubiquitous in nature. Higher plants contain three evolutionarily distinct CA families, αCAs, βCAs, and γCAs, where each family is represented by multiple isoforms in all species. Alternative splicing of CA transcripts appears common; consequently, the number of functional CA isoforms in a species may exceed the number of genes. CAs are expressed in numerous plant tissues and in different cellular locations. The most prevalent CAs are those in the chloroplast, cytosol, and mitochondria. This diversity in location is paralleled in the many physiological and biochemical roles that CAs play in plants. In this review, the number and types of CAs in C3, C4, and crassulacean acid metabolism (CAM) plants are considered, and the roles of the α and γCAs are briefly discussed. The remainder of the review focuses on plant βCAs and includes the identification of homologs between species using phylogenetic approaches, a consideration of the inter- and intracellular localization of the proteins, along with the evidence for alternative splice forms. Current understanding of βCA tissue-specific expression patterns and what controls them are reviewed, and the physiological roles for which βCAs have been implicated are presented.

Gene mining in halophytes: functional identification of stress tolerance genes in Lepidium crassifolium

Extremophile plants are valuable sources of genes conferring tolerance traits, which can be explored to improve stress tolerance of crops. Lepidium crassifolium is a halophytic relative of the model plant Arabidopsis thaliana, and displays tolerance to salt, osmotic and oxidative stresses. We have employed the modified Conditional cDNA Overexpression System to transfer a cDNA library from L. crassifolium to the glycophyte A. thaliana. By screening for salt, osmotic and oxidative stress tolerance through in vitro growth assays and non-destructive chlorophyll fluorescence imaging, 20 Arabidopsis lines were identified with superior performance under restrictive conditions. Several cDNA inserts were cloned and confirmed to be responsible for the enhanced tolerance by analysing independent transgenic lines. Examples include full-length cDNAs encoding proteins with high homologies to GDSL-lipase/esterase or acyl CoA-binding protein or proteins without known function, which could confer tolerance to one or several stress conditions. Our results confirm that random gene transfer from stress tolerant to sensitive plant species is a valuable tool to discover novel genes with potential for biotechnological applications.

A New Insight of Salt Stress Signaling in Plant

Many studies have been conducted to understand plant stress responses to salinity because irrigation-dependent salt accumulation compromises crop productivity and also to understand the mechanism through which some plants thrive under saline conditions. As mechanistic understanding has increased during the last decades, discovery-oriented approaches have begun to identify genetic determinants of salt tolerance. In addition to osmolytes, osmoprotectants, radical detoxification, ion transport systems, and changes in hormone levels and hormone-guided communications, the Salt Overly Sensitive (SOS) pathway has emerged to be a major defense mechanism. However, the mechanism by which the components of the SOS pathway are integrated to ultimately orchestrate plant-wide tolerance to salinity stress remains unclear. A higher-level control mechanism has recently emerged as a result of recognizing the involvement of GIGANTEA (GI), a protein involved in maintaining the plant circadian clock and control switch in flowering. The loss of GI function confers high tolerance to salt stress via its interaction with the components of the SOS pathway. The mechanism underlying this observation indicates the association between GI and the SOS pathway and thus, given the key influence of the circadian clock and the pathway on photoperiodic flowering, the association between GI and SOS can regulate growth and stress tolerance. In this review, we will analyze the components of the SOS pathways, with emphasis on the integration of components recognized as hallmarks of a halophytic lifestyle.

The Cytoplasmic Carbonic Anhydrases βCA2 and βCA4 Are Required for Optimal Plant Growth at Low CO2

Carbonic anhydrases (CAs) are zinc metalloenzymes that interconvert CO2 and HCO3 (-) In plants, both α- and β-type CAs are present. We hypothesize that cytoplasmic βCAs are required to modulate inorganic carbon forms needed in leaf cells for carbon-requiring reactions such as photosynthesis and amino acid biosynthesis. In this report, we present evidence that βCA2 and βCA4 are the two most abundant cytoplasmic CAs in Arabidopsis (Arabidopsis thaliana) leaves. Previously, βCA4 was reported to be localized to the plasma membrane, but here, we show that two forms of βCA4 are expressed in a tissue-specific manner and that the two proteins encoded by βCA4 localize to two different regions of the cell. Comparing transfer DNA knockout lines with wild-type plants, there was no reduction in the growth rates of the single mutants, βca2 and βca4 However, the growth rate of the double mutant, βca2βca4, was reduced significantly when grown at 200 μL L(-1) CO2 The reduction in growth of the double mutant was not linked to a reduction in photosynthetic rate. The amino acid content of leaves from the double mutant showed marked reduction in aspartate when compared with the wild type and the single mutants. This suggests the cytoplasmic CAs play an important but not previously appreciated role in amino acid biosynthesis.

Heritable gene expression differences between apomictic clone members in Taraxacum officinale: Insights into early stages of evolutionary divergence in asexual plants

BACKGROUND

Asexual reproduction has the potential to enhance deleterious mutation accumulation and to constrain adaptive evolution. One source of mutations that can be especially relevant in recent asexuals is activity of transposable elements (TEs), which may have experienced selection for high transposition rates in sexual ancestor populations. Predictions of genomic divergence under asexual reproduction therefore likely include a large contribution of transposable elements but limited adaptive divergence. For plants empirical insight into genome divergence under asexual reproduction remains limited. Here, we characterize expression divergence between clone members of a single apomictic lineage of the common dandelion (Taraxacum officinale) to contribute to our knowledge of genome evolution under asexuality.

RESULTS

Using RNA-Seq, we show that about one third of heritable divergence within the apomictic lineage is driven by TEs and TE-related gene activity. In addition, we identify non-random transcriptional differences in pathways related to acyl-lipid and abscisic acid metabolisms which might reflect functional divergence within the apomictic lineage. We analyze SNPs in the transcriptome to assess genetic divergence between the apomictic clone members and reveal that heritable expression differences between the accessions are not explained simply by genome-wide genetic divergence.

CONCLUSION

The present study depicts a first effort towards a more complete understanding of apomictic plant genome evolution. We identify abundant TE activity and ecologically relevant functional genes and pathways affecting heritable within-lineage expression divergence. These findings offer valuable resources for future work looking at epigenetic silencing and Cis-regulation of gene expression with particular emphasis on the effects of TE activity on asexual species' genome.

Evolutionary Dynamics of the Leucine-Rich Repeat Receptor-Like Kinase (LRR-RLK) Subfamily in Angiosperms

Gene duplications are an important factor in plant evolution, and lineage-specific expanded (LSE) genes are of particular interest. Receptor-like kinases expanded massively in land plants, and leucine-rich repeat receptor-like kinases (LRR-RLK) constitute the largest receptor-like kinases family. Based on the phylogeny of 7,554 LRR-RLK genes from 31 fully sequenced flowering plant genomes, the complex evolutionary dynamics of this family was characterized in depth. We studied the involvement of selection during the expansion of this family among angiosperms. LRR-RLK subgroups harbor extremely contrasting rates of duplication, retention, or loss, and LSE copies are predominantly found in subgroups involved in environmental interactions. Expansion rates also differ significantly depending on the time when rounds of expansion or loss occurred on the angiosperm phylogenetic tree. Finally, using a dN/dS-based test in a phylogenetic framework, we searched for selection footprints on LSE and single-copy LRR-RLK genes. Selective constraint appeared to be globally relaxed at LSE genes, and codons under positive selection were detected in 50% of them. Moreover, the leucine-rich repeat domains, and specifically four amino acids in them, were found to be the main targets of positive selection. Here, we provide an extensive overview of the expansion and evolution of this very large gene family.

Functional Conservation in the SIAMESE-RELATED Family of Cyclin-Dependent Kinase Inhibitors in Land Plants

The best-characterized members of the plant-specific SIAMESE-RELATED (SMR) family of cyclin-dependent kinase inhibitors regulate the transition from the mitotic cell cycle to endoreplication, also known as endoreduplication, an altered version of the cell cycle in which DNA is replicated without cell division. Some other family members are implicated in cell cycle responses to biotic and abiotic stresses. However, the functions of most SMRs remain unknown, and the specific cyclin-dependent kinase complexes inhibited by SMRs are unclear. Here, we demonstrate that a diverse group of SMRs, including an SMR from the bryophyte Physcomitrella patens, can complement an Arabidopsis thaliana siamese (sim) mutant and that both Arabidopsis SIM and P. patens SMR can inhibit CDK activity in vitro. Furthermore, we show that Arabidopsis SIM can bind to and inhibit both CDKA;1 and CDKB1;1. Finally, we show that SMR2 acts to restrict cell proliferation during leaf growth in Arabidopsis and that SIM, SMR1/LGO, and SMR2 play overlapping roles in controlling the transition from cell division to endoreplication during leaf development. These results indicate that differences in SMR function in plant growth and development are primarily due to differences in transcriptional and posttranscriptional regulation, rather than to differences in fundamental biochemical function.

Improving the Annotation of Arabidopsis lyrata Using RNA-Seq Data

Gene model annotations are important community resources that ensure comparability and reproducibility of analyses and are typically the first step for functional annotation of genomic regions. Without up-to-date genome annotations, genome sequences cannot be used to maximum advantage. It is therefore essential to regularly update gene annotations by integrating the latest information to guarantee that reference annotations can remain a common basis for various types of analyses. Here, we report an improvement of the Arabidopsis lyrata gene annotation using extensive RNA-seq data. This new annotation consists of 31,132 protein coding gene models in addition to 2,089 genes with high similarity to transposable elements. Overall, ~87% of the gene models are corroborated by evidence of expression and 2,235 of these models feature multiple transcripts. Our updated gene annotation corrects hundreds of incorrectly split or merged gene models in the original annotation, and as a result the identification of alternative splicing events and differential isoform usage are vastly improved.

A novel workflow correlating RNA-seq data to Phythophthora infestans resistance levels in wild Solanum species and potato clones

Comparative transcriptomics between species can provide valuable understanding of plant-pathogen interactions. Here, we focus on wild Solanum species and potato clones with varying degree of resistance against Phytophthora infestans, which causes the devastating late blight disease in potato. The transcriptomes of three wild Solanum species native to Southern Sweden, Solanum dulcamara, Solanum nigrum, and Solanum physalifolium were compared to three potato clones, Desiree (cv.), SW93-1015 and Sarpo Mira. Desiree and S. physalifolium are susceptible to P. infestans whereas the other four have different degrees of resistance. By building transcript families based on de novo assembled RNA-seq across species and clones and correlating these to resistance phenotypes, we created a novel workflow to identify families with expanded or depleted number of transcripts in relation to the P. infestans resistance level. Analysis was facilitated by inferring functional annotations based on the family structure and semantic clustering. More transcript families were expanded in the resistant clones and species and the enriched functions of these were associated to expected gene ontology (GO) terms for resistance mechanisms such as hypersensitive response, host programmed cell death and endopeptidase activity. However, a number of unexpected functions and transcripts were also identified, for example transmembrane transport and protein acylation expanded in the susceptible group and a cluster of Zinc knuckle family proteins expanded in the resistant group. Over 400 expressed putative resistance (R-)genes were identified and resistant clones Sarpo Mira and SW93-1015 had ca 25% more expressed putative R-genes than susceptible cultivar Desiree. However, no differences in numbers of susceptibility (S-)gene homologs were seen between species and clones. In addition, we identified P. infestans transcripts including effectors in the early stages of P. infestans-Solanum interactions.

The complete chloroplast genome of Schrenkiella parvula (Brassicaceae)

Schrenkiella parvula is an Arabidopsis-related model species used here for studying plant stress tolerance. In this study, the complete chloroplast genome sequence of S. parvula has been reported for the first time. The total length of the chloroplast genome was 153 979 bp, which had a typical quadripartite structure. The annotated plastid genome includes 87 protein-coding genes, 39 tRNA genes and 8 ribosomal RNA genes. The evolutionary relationships revealed by our phylogenetic analysis indicated that S. parvula is closer to the Brassiceae species when compared with Eutrema salsugineum.

Cell type-specific responses to salinity - the epidermal bladder cell transcriptome of Mesembryanthemum crystallinum

Mesembryanthemum crystallinum (ice plant) exhibits extreme tolerance to salt. Epidermal bladder cells (EBCs), developing on the surface of aerial tissues and specialized in sodium sequestration and other protective functions, are critical for the plant's stress adaptation. We present the first transcriptome analysis of EBCs isolated from intact plants, to investigate cell type-specific responses during plant salt adaptation. We developed a de novo assembled, nonredundant EBC reference transcriptome. Using RNAseq, we compared the expression patterns of the EBC-specific transcriptome between control and salt-treated plants. The EBC reference transcriptome consists of 37 341 transcript-contigs, of which 7% showed significantly different expression between salt-treated and control samples. We identified significant changes in ion transport, metabolism related to energy generation and osmolyte accumulation, stress signalling, and organelle functions, as well as a number of lineage-specific genes of unknown function, in response to salt treatment. The salinity-induced EBC transcriptome includes active transcript clusters, refuting the view of EBCs as passive storage compartments in the whole-plant stress response. EBC transcriptomes, differing from those of whole plants or leaf tissue, exemplify the importance of cell type-specific resolution in understanding stress adaptive mechanisms.

The evolution of halophytes, glycophytes and crops, and its implications for food security under saline conditions

The effective development of salt tolerant crops requires an understanding that the evolution of halophytes, glycophytes and our major grain crops has involved significantly different processes. Halophytes (and other edaphic endemics) generally arose through colonization of habitats in severe disequilibrium by pre-adapted individuals, rather than by gradual adaptation from populations of 'glycophytes'. Glycophytes, by contrast, occur in low sodium ecosystems, where sodium was and is the major limiting nutrient in herbivore diets, suggesting that their evolution reflects the fact that low sodium individuals experienced lower herbivory and had higher fitness. For domestication/evolution of crop plants, the selective pressure was human imposed and involved humans co-opting functions of defense and reproductive security. Unintended consequences of this included loss of tolerance to various stresses and loss of the genetic variability needed to correct that. Understanding, combining and manipulating all three modes of evolution are now critical to the development of salt tolerant crops, particularly those that will offer food security in countries with few economic resources and limited infrastructure. Such efforts will require exploiting the genetic structures of recently evolved halophytes, the genetic variability of model plants, and endemic halophytes and 'minor' crops that already exist.

Genetic mechanisms of abiotic stress tolerance that translate to crop yield stability

Crop yield reduction as a consequence of increasingly severe climatic events threatens global food security. Genetic loci that ensure productivity in challenging environments exist within the germplasm of crops, their wild relatives and species that are adapted to extreme environments. Selective breeding for the combination of beneficial loci in germplasm has improved yields in diverse environments throughout the history of agriculture. An effective new paradigm is the targeted identification of specific genetic determinants of stress adaptation that have evolved in nature and their precise introgression into elite varieties. These loci are often associated with distinct regulation or function, duplication and/or neofunctionalization of genes that maintain plant homeostasis.

Extensive copy number variations in admixed Indian population of African ancestry: potential involvement in adaptation

Admixture mapping has been enormously resourceful in identifying genetic variations linked to phenotypes, adaptation, and diseases. In this study through analysis of copy number variable regions (CNVRs), we report extensive restructuring in the genomes of the recently admixed African-Indian population (OG-W-IP) that inhabits a highly saline environment in Western India. The study included subjects from OG-W-IP (OG), five different Indian and three HapMap populations that were genotyped using Affymetrix version 6.0 arrays. Copy number variations (CNVs) detected using Birdsuite were used to define CNVRs. Population structure with respect to CNVRs was delineated using random forest approach. OG genomes have a surprising excess of CNVs in comparison to other studied populations. Individual ancestry proportions computed using STRUCTURE also reveals a unique genetic component in OGs. Population structure analysis with CNV genotypes indicates OG to be distant from both the African and Indian ancestral populations. Interestingly, it shows genetic proximity with respect to CNVs to only one Indian population IE-W-LP4, which also happens to reside in the same geographical region. We also observe a significant enrichment of molecular processes related to ion binding and receptor activity in genes encompassing OG-specific CNVRs. Our results suggest that retention of CNVRs from ancestral natives and de novo acquisition of CNVRs could accelerate the process of adaptation especially in an extreme environment. Additionally, this population would be enormously useful for dissecting genes and delineating the involvement of CNVs in salt adaptation.

Transcriptomic analysis of a psammophyte food crop, sand rice (Agriophyllum squarrosum) and identification of candidate genes essential for sand dune adaptation

Background

Sand rice (Agriophyllum squarrosum) is an annual desert plant adapted to mobile sand dunes in arid and semi-arid regions of Central Asia. The sand rice seeds have excellent nutrition value and have been historically consumed by local populations in the desert regions of northwest China. Sand rice is a potential food crop resilient to ongoing climate change; however, partly due to the scarcity of genetic information, this species has undergone only little agronomic modifications through classical breeding during recent years.

Results

We generated a deep transcriptomic sequencing of sand rice, which uncovers 67,741 unigenes. Phylogenetic analysis based on 221 single-copy genes showed close relationship between sand rice and the recently domesticated crop sugar beet. Transcriptomic comparisons also showed a high level of global sequence conservation between these two species. Conservation of sand rice and sugar beet orthologs assigned to response to salt stress gene ontology term suggests that sand rice is also a potential salt tolerant plant. Furthermore, sand rice is far more tolerant to high temperature. A set of genes likely relevant for resistance to heat stress, was functionally annotated according to expression levels, sequence annotation, and comparisons corresponding transcriptome profiling results in Arabidopsis.

Conclusions

The present work provides abundant genomic information for functional dissection of the important traits in sand rice. Future screening the genetic variation among different ecotypes and constructing a draft genome sequence will further facilitate agronomic trait improvement and final domestication of sand rice.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-872) contains supplementary material, which is available to authorized users.

Transcriptome sequencing and analysis of leaf tissue of Avicennia marina using the Illumina platform

Avicennia marina is a widely distributed mangrove species that thrives in high-salinity habitats. It plays a significant role in supporting coastal ecosystem and holds unique potential for studying molecular mechanisms underlying ecological adaptation. Despite and sometimes because of its numerous merits, this species is facing increasing pressure of exploitation and deforestation. Both study on adaptation mechanisms and conservation efforts necessitate more genomic resources for A. marina. In this study, we used Illumina sequencing of an A. marina foliar cDNA library to generate a transcriptome dataset for gene and marker discovery. We obtained 40 million high-quality reads and assembled them into 91,125 unigenes with a mean length of 463 bp. These unigenes covered most of the publicly available A. marina Sanger ESTs and greatly extended the repertoire of transcripts for this species. A total of 54,497 and 32,637 unigenes were annotated based on homology to sequences in the NCBI non-redundant and the Swiss-prot protein databases, respectively. Both Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed some transcriptomic signatures of stress adaptation for this halophytic species. We also detected an extraordinary amount of transcripts derived from fungal endophytes and demonstrated the utility of transcriptome sequencing in surveying endophyte diversity without isolating them out of plant tissues. Additionally, we identified 3,423 candidate simple sequence repeats (SSRs) from 3,141 unigenes with a density of one SSR locus every 8.25 kb sequence. Our transcriptomic data will provide valuable resources for ecological, genetic and evolutionary studies in A. marina.