Natural variation in ozone sensitivity among Arabidopsis thaliana accessions and its relation to stomatal conductance

Genetic variation between naturally occurring populations provides a unique source to unravel the complex mechanisms of stress tolerance. Here, we have analysed O₃ sensitivity of 93 natural Arabidopsis thaliana accessions together with five O₃‐sensitive mutants to acute O₃ exposure. The variation in O₃ sensitivity among the natural accessions was much higher than among the O₃‐sensitive mutants and corresponding wild types. A subset of nine accessions with major variation in their O₃ responses was studied in more detail. Among the traits assayed, stomatal conductance (g_st) was an important factor determining O₃ sensitivity of the selected accessions. The most O₃‐sensitive accession, Cvi‐0, had constitutively high g_st, leading to high initial O₃ uptake rate and dose received during the first 30 min of exposure. Analyzing O₃‐induced changes in stress hormone concentrations indicated that jasmonate (JA) concentration was also positively correlated with leaf damage. Quantitative trait loci (QTL) mapping in a Col‐0 × Cvi‐0 recombinant inbred line (RIL) population identified three QTLs for O₃ sensitivity, and one for high water loss of Cvi‐0. The major O₃ QTL mapped to the same position as the water loss QTL further supporting the role of stomata in regulating O₃ entry and damage.     

Functional screening of willow alleles in Arabidopsis combined with QTL mapping in willow (Salix) identifies SxMAX4 as a coppicing response gene

Willows (Salix spp.) are important biomass crops due to their ability to grow rapidly with low fertilizer inputs and ease of cultivation in short‐rotation coppice cycles. They are relatively undomesticated and highly diverse, but functional testing to identify useful allelic variation is time‐consuming in trees and transformation is not yet possible in willow. Arabidopsis is heralded as a model plant from which knowledge can be transferred to advance the improvement of less tractable species. Here, knowledge and methodologies from Arabidopsis were successfully used to identify a gene influencing stem number in coppiced willows, a complex trait of key biological and industrial relevance. The strigolactone‐related More AXillary growth (MAX) genes were considered candidates due to their role in shoot branching. We previously demonstrated that willow and Arabidopsis show similar response to strigolactone and that transformation rescue of Arabidopsis max mutants with willow genes could be used to detect allelic differences. Here, this approach was used to screen 45 SxMAX1, SxMAX2, SxMAX3 and SxMAX4 alleles cloned from 15 parents of 11 mapping populations varying in shoot‐branching traits. Single‐nucleotide polymorphism (SNP) frequencies were locus dependent, ranging from 29.2 to 74.3 polymorphic sites per kb. SxMAX alleles were 98%–99% conserved at the amino acid level, but different protein products varying in their ability to rescue Arabidopsis max mutants were identified. One poor rescuing allele, SxMAX4D, segregated in a willow mapping population where its presence was associated with increased shoot resprouting after coppicing and colocated with a QTL for this trait.     

A naturally occurring promoter polymorphism of the Arabidopsis FUM2 gene causes expression variation,and is associated with metabolic and growth traits

Fumarate and malate are known intermediates of the TCA cycle, a mitochondrial metabolic pathway generating NADH for respiration. Arabidopsis thaliana and other Brassicaceae contain an additional cytosolic fumarase (FUM2) that functions in carbon assimilation and nitrogen use. Here, we report the identification of a hitherto unknown FUM2 promoter insertion/deletion (InDel) polymorphism found between the Col‐0 and C24 accessions, which also divides a large number of Arabidopsis accessions carrying either the Col‐0 or the C24 allele. The polymorphism consists of two stretches of 2.1 and 3.8 kb, which are both absent from the promotor region of Col‐0 FUM2. By analysing mutants as well as mapping and natural populations with contrasting FUM2 alleles, the promotor insertion was linked to reduced FUM2 mRNA expression, reduced fumarase activity and reduced fumarate/malate ratio in leaves. In a large population of 174 natural accessions, the polymorphism was also found to be associated with the fumarate/malate ratio, malate and fumarate levels, and with dry weight at 15 days after sowing (DAS). The association with biomass production was confirmed in an even larger (251) accession population for dry weight at 22 DAS. The dominant Col‐0 allele that results in increased fumarate/malate ratios and enhanced biomass production is predominantly found in central/eastern European accessions, whereas the C24 type allele is prevalent on the Iberian Peninsula, west of the Rhine and in the British Isles. Our findings support the role of FUM2 in diurnal carbon storage, and point to a growth advantage of accessions carrying the FUM2 Col‐0 allele.     

Genetic Variation for Adaptive Traits in Bottlebrush Squirreltail in the Northern Intermountain West,United States

Adaptive‐trait correlations in plant ecology are often calculated among species, but in order to develop and characterize plant materials of target species for restoration, intraspecific comparisons are of greatest relevance. Elymus elymoides (Raf.) Swezey (bottlebrush squirreltail) is an important component of sagebrush‐steppe communities in the northern Intermountain West, United States. We evaluated 32 accessions of E. elymoides subspecies C, a newly recognized unnamed taxon, in the field and greenhouse. Our objectives were to assess genetic diversity for putatively adaptive traits, to elucidate biological relationships among biomass, morphological, and phenological traits through correlation analysis, and to gather evidence suggesting whether these traits might be truly adaptive, that is, related to collection‐site variables. We observed a positive correlation (r = 0.73;p < 0.01) between greenhouse shoot and root biomass among accessions, suggesting that shoot and root biomass are not in an inherent trade‐off relationship across accessions. In addition, accessions with higher greenhouse shoot biomass possessed lower specific leaf area (r = ?0.43;p < 0.05) and lower specific root length (r = ?0.47,p < 0.05). Correlations between greenhouse and field‐measured productivity traits were not significant (p > 0.05), indicating seedling performance is not predictive of mature‐plant performance. Elevation was the collection‐site variable most closely correlated with plant‐measured traits, particularly phenological dates, whereas average annual precipitation was the least significant variable. Therefore, elevation may be used as an easily applied metric to match subspecies C plant material to restoration site in the northern Intermountain West.     

Patterns of univariate and multivariate plasticity to elevated carbon dioxide in six European populations of Arabidopsis thaliana

The impact of elevated carbon dioxide on plants is a growing concern in evolutionary ecology and global change biology. Characterizing patterns of phenotypic integration and multivariate plasticity to elevated carbon dioxide can provide insights into ecological and evolutionary dynamics in future human‐altered environments. Here, we examined univariate and multivariate responses to carbon enrichment in six functional traits among six European accessions of Arabidopsis thaliana. We detected phenotypic plasticity in both univariate and multivariate phenotypes, but did not find significant variation in plasticity (genotype by environment interactions) within or among accessions. Eigenvector, eigenvalue variance, and common principal components analyses showed that elevated carbon dioxide altered patterns of trait covariance, reduced the strength of phenotypic integration, and decreased population‐level differentiation in the multivariate phenotype. Our data suggest that future carbon dioxide conditions may influence evolutionary dynamics in natural populations of A. thaliana.     

The wheat TaGBF1 gene is involved in the blue‐light response and salt tolerance

Light and abiotic stress both strongly modulate plant growth and development. However, the effect of light‐responsive factors on growth and abiotic stress responses in wheat (Triticum aestivum) is unknown. G–box binding factors (GBFs) are blue light‐specific components, but their function in abiotic stress responses has not been studied. Here we identified a wheat GBF1 gene that mediated both the blue light‐ and abiotic stress‐responsive signaling pathways. TaGBF1 was inducible by blue light, salt and exposure to abscisic acid (ABA). TaGBF1 interacted with a G–box light‐responsive element in vitro and promoted a blue‐light response in wheat and Aradidopsis thaliana. Both TaGBF1 over‐expression in wheat and its heterologous expression in A. thaliana heighten sensitivity to salinity and ABA, but its knockdown in wheat conferred resistance to high salinity and ABA. The expression of AtABI5, a key component of the ABA signaling pathway in A. thaliana, and its homolog Wabi5 in wheat was increased by transgenic expression of TaGBF1. The hypersensitivity to salt and ABA caused by TaGBF1 was not observed in the abi5 mutant background, showing that ABI5 is the mediator in TaGBF1‐induced abiotic stress responses. However, the hypersensitivity to salt conferred by TaGBF1 is not dependent on light. This suggests that TaGBF1 is a common component of blue light‐ and abiotic stress‐responsive signaling pathways.     

Imputation of 3 million SNPs in the Arabidopsis regional mapping population

Natural variation has become a prime resource to identify genetic variants that contribute to phenotypic variation. The regional mapping (RegMap) population is one of the most important populations for studying natural variation in Arabidopsis thaliana, and has been used in a large number of association studies and in studies on climatic adaptation. However, only 413 RegMap accessions have been completely sequenced, as part of the 1001 Genomes (1001G) Project, while the remaining 894 accessions have only been genotyped with the Affymetrix 250k chip. As a consequence, most association studies involving the RegMap are either restricted to the sequenced accessions, reducing power, or rely on a limited set of SNPs. Here we impute millions of SNPs to the 894 accessions that are exclusive to the RegMap, using the 1135 accessions of the 1001G Project as the reference panel. We assess imputation accuracy using a novel cross‐validation scheme, which we show provides a more reliable measure of accuracy than existing methods. After filtering out low accuracy SNPs, we obtain high‐quality genotypic information for 2029 accessions and 3 million markers. To illustrate the benefits of these imputed data, we reconducted genome‐wide association studies on five stress‐related traits and could identify novel candidate genes.     

Genetic interactions reveal the antagonistic roles of FT/TSF and TFL1 in the determination of inflorescence meristem identity in Arabidopsis

During the transition to the reproductive phase, the shoot apical meristem switches from the developmental program that generates vegetative organs to instead produce flowers. In this study, we examined the genetic interactions of FLOWERING LOCUS T (FT)/TWIN SISTER OF FT (TSF) and TERMINAL FLOWER 1 (TFL1) in the determination of inflorescence meristem identity in Arabidopsis thaliana. The ft‐10 tsf‐1 mutants produced a compact inflorescence surrounded by serrated leaves (hyper‐vegetative shoot) at the early bolting stage, as did plants overexpressing TFL1. Plants overexpressing FT or TSF (or both FT and TFL1) generated a terminal flower, as did tfl1‐20 mutants. The terminal flower formed in tfl1‐20 mutants converted to a hyper‐vegetative shoot in ft‐10 tsf‐1 mutants. Grafting ft‐10 tsf‐1 or ft‐10 tsf‐1 tfl1‐20 mutant scions to 35S::FT rootstock plants produced a normal inflorescence and a terminal flower in the scion plants, respectively, although both scions showed similar early flowering. Misexpression of FT in the vasculature and in the shoot apex in wild‐type plants generated a normal inflorescence and a terminal flower, respectively. By contrast, in ft‐10 tsf‐1 mutants the vasculature‐specific misexpression of FT converted the hyper‐vegetative shoot to a normal inflorescence, and in the ft‐10 tsf‐1 tfl1‐20 mutants converted the shoot to a terminal flower. TFL1 levels did not affect the inflorescence morphology caused by FT/TSF overexpression at the early bolting stage. Taking these results together, we proposed that FT/TSF and TFL1 play antagonistic roles in the determination of inflorescence meristem identity, and that FT/TSF are more important than TFL1 in this process.     

A VIGS screen identifies immunity in the Arabidopsis Pla‐1 accession to viruses in two different genera of the Geminiviridae

Geminiviruses are DNA viruses that cause severe crop losses in different parts of the world, and there is a need for genetic sources of resistance to help combat them. Arabidopsis has been used as a source for virus‐resistant genes that derive from alterations in essential host factors. We used a virus‐induced gene silencing (VIGS) vector derived from the geminivirus Cabbage leaf curl virus (CaLCuV) to assess natural variation in virus–host interactions in 190 Arabidopsis accessions. Silencing of CH‐42, encoding a protein needed to make chlorophyll, was used as a visible marker to discriminate asymptomatic accessions from those showing resistance. There was a wide range in symptom severity and extent of silencing in different accessions, but two correlations could be made. Lines with severe symptoms uniformly lacked extensive VIGS, and lines that showed attenuated symptoms over time (recovery) showed a concomitant increase in the extent of VIGS. One accession, Pla‐1, lacked both symptoms and silencing, and was immune to wild‐type infectious clones corresponding to CaLCuV or Beet curly top virus (BCTV), which are classified in different genera in the Geminiviridae. It also showed resistance to the agronomically important Tomato yellow leaf curl virus (TYLCV). Quantitative trait locus mapping of a Pla‐1 X Col‐0 F₂ population was used to detect a major peak on chromosome 1, which is designated gip‐1 (geminivirus immunity Pla‐1‐1). The recessive nature of resistance to CaLCuV and the lack of obvious candidate genes near the gip‐1 locus suggest that a novel resistance gene(s) confers immunity.     

BSKs are partially redundant positive regulators of brassinosteroid signaling in Arabidopsis

Arabidopsis thaliana brassinosteroid signaling kinases (BSKs) constitute a receptor‐like cytoplasmic kinase sub‐family (RLCK‐XII) with 12 members. Previous analysis demonstrated a positive role for BSK1 and BSK3 in the initial steps of brassinosteroid (BR) signal transduction. To investigate the function of BSKs in plant growth and BR signaling, we characterized T‐DNA insertion lines for eight BSK genes (BSK1–BSK8) and multiple mutant combinations. Simultaneous elimination of three BSK genes caused alterations in growth and the BR response, and the most severe phenotypes were observed in the bsk3,4,7,8 quadruple and bsk3,4,6,7,8 pentuple mutants, which displayed reduced rosette size, leaf curling and enhanced leaf inclination. In addition, upon treatment with 24‐epibrassinolide, these mutants showed reduced hypocotyl elongation, enhanced root growth and alteration in the expression of BR‐responsive genes. Some mutant combinations also showed antagonistic interactions. In support of a redundant function in BR signaling, multiple BSKs interacted in vivo with the BR receptor BRI1, and served as its phosphorylation substrates in vitro. The BIN2 and BIL2 GSK3‐like kinases, which are negative regulators of BR signaling, interacted in vivo with BSKs and phosphorylated them in vitro, probably at different sites to BRI1. This study demonstrates redundant biological functions for BSKs, and suggests the existence of a regulatory link between BSKs and GSK3‐like kinases.     

Responses of antioxidant enzymes to cold and high light are not correlated to freezing tolerance in natural accessions of Arabidopsis thaliana

Low temperatures and high light cause imbalances in primary and secondary reactions of photosynthesis, and thus can result in oxidative stress. Plants employ a range of low‐molecular weight antioxidants and antioxidant enzymes to prevent oxidative damage, and antioxidant defence is considered an important component of stress tolerance. To figure out whether oxidative stress and antioxidant defence are key factors defining the different cold acclimation capacities of natural accessions of the model plant Arabidopsis thaliana, we investigated hydrogen peroxide (H₂O₂) production, antioxidant enzyme activity and lipid peroxidation during a time course of cold treatment and exposure to high light in four differentially cold‐tolerant natural accessions of Arabidopsis (C24, Nd, Rsch, Te) that span the European distribution range of the species. All accessions except Rsch (from Russia) had elevated H₂O₂ in the cold, indicating that production of reactive oxygen species is part of the cold response in Arabidopsis. Glutathione reductase activity increased in all but Rsch, while ascorbate peroxidase and superoxide dismutase were unchanged and catalase decreased in all but Rsch. Under high light, the Scandinavian accession Te had elevated levels of H₂O₂. Te appeared most sensitive to oxidative stress, having higher malondialdehyde (MDA) levels in the cold and under high light, while only high light caused elevated MDA in the other accessions. Although the most freezing‐tolerant, Te had the highest sensitivity to oxidative stress. No correlation was found between freezing tolerance and activity of antioxidant enzymes in the four accessions investigated, arguing against a key role for antioxidant defence in the differential cold acclimation capacities of Arabidopsis accessions.     

Genetic and evolution analysis of extrafloral nectary in cotton

Extrafloral nectaries are a defence trait that plays important roles in plant–animal interactions. Gossypium species are characterized by cellular grooves in leaf midribs that secret large amounts of nectar. Here, with a panel of 215 G. arboreum accessions, we compared extrafloral nectaries to nectariless accessions to identify a region of Chr12 that showed strong differentiation and overlapped with signals from GWAS of nectaries. Fine mapping of an F₂ population identified GaNEC1, encoding a PB1 domain‐containing protein, as a positive regulator of nectary formation. An InDel, encoding a five amino acid deletion, together with a nonsynonymous substitution, was predicted to cause 3D structural changes in GaNEC1 protein that could confer the nectariless phenotype. mRNA‐Seq analysis showed that JA‐related genes are up‐regulated and cell wall‐related genes are down‐regulated in the nectary. Silencing of GaNEC1 led to a smaller size of foliar nectary phenotype. Metabolomics analysis identified more than 400 metabolites in nectar, including expected saccharides and amino acids. The identification of GaNEC1 helps establish the network regulating nectary formation and nectar secretion, and has implications for understanding the production of secondary metabolites in nectar. Our results will deepen our understanding of plant–mutualism co‐evolution and interactions, and will enable utilization of a plant defence trait in cotton breeding efforts.     

Leaf functional trait variation associated with salt tolerance in perennial ryegrass

Salinity is one of major environmental stresses that dramatically threaten plant growth, and variations in genetic structure and functional traits have important effects on the salt tolerance of perennial ryegrass (Lolium perenne L.). The objectives of this study were to: (i) assess the inter‐clonal variation of functional traits of accessions among geographic groups or between wild and commercial groups in response to salt stress; (ii) develop a mathematical model to effectively assess salt tolerance of perennial ryegrass accessions originating from different geographic populations; and (iii) determine the relation between spatial genetic structure and salt tolerance in perennial ryegrass. Wide variations were found among the accessions for seven functional traits. One regression model (F = 0.49 × F₁ + 0.303 × F₂ + 0.207 × F₃) was established to ascertain salt tolerance of each accession. The highest variation of the traits and salt tolerance were obtained for accessions from the European group. Wild accessions exhibited more variation in functional traits and salt tolerance than commercial cultivars. Both molecular marker techniques and functional traits were used to conduct phylogenetic analysis, and the majority of accessions from the same or adjacent regions were clustered into the same group or subgroup. The perennial ryegrass accessions with similar salt tolerance had a close phylogenetic background. The patterns in functional trait variations associated with salt tolerance might allow acceleration of the process for improving salt stress resistance in perennial ryegrass.     

Expression of Arabidopsis Bax Inhibitor‐1 in transgenic sugarcane confers drought tolerance

The sustainability of global crop production is critically dependent on improving tolerance of crop plants to various types of environmental stress. Thus, identification of genes that confer stress tolerance in crops has become a top priority especially in view of expected changes in global climatic patterns. Drought stress is one of the abiotic stresses that can result in dramatic loss of crop productivity. In this work, we show that transgenic expression of a highly conserved cell death suppressor, Bax Inhibitor‐1 from Arabidopsis thaliana (AtBI‐1), can confer increased tolerance of sugarcane plants to long‐term (>20 days) water stress conditions. This robust trait is correlated with an increased tolerance of the transgenic sugarcane plants, especially in the roots, to induction of endoplasmic reticulum (ER) stress by the protein glycosylation inhibitor tunicamycin. Our findings suggest that suppression of ER stress in C₄ grasses, which include important crops such as sorghum and maize, can be an effective means of conferring improved tolerance to long‐term water deficit. This result could potentially lead to improved resilience and yield of major crops in the world.     

Two combinatorial patterns of telomere histone marks in plants with canonical and non‐canonical telomere repeats

Telomeres, nucleoprotein structures at the ends of linear eukaryotic chromosomes, are crucial for the maintenance of genome integrity. In most plants, telomeres consist of conserved tandem repeat units comprising the TTTAGGG motif. Recently, non‐canonical telomeres were described in several plants and plant taxons, including the carnivorous plant Genlisea hispidula (TTCAGG/TTTCAGG), the genus Cestrum (Solanaceae; TTTTTTAGGG), and plants from the Asparagales order with either a vertebrate‐type telomere repeat TTAGGG or Allium genus‐specific CTCGGTTATGGG repeat. We analyzed epigenetic modifications of telomeric histones in plants with canonical and non‐canonical telomeres, and further in telomeric chromatin captured from leaves of Nicotiana benthamiana transiently transformed by telomere CRISPR‐dCas9‐eGFP, and of Arabidopsis thaliana stably transformed with TALE_telo C‐3×GFP. Two combinatorial patterns of telomeric histone modifications were identified: (i) an Arabidopsis‐like pattern (A. thaliana, G. hispidula, Genlisea nigrocaulis, Allium cepa, Narcissus pseudonarcissus, Petunia hybrida, Solanum tuberosum, Solanum lycopersicum) with telomeric histones decorated predominantly by H3K9me2; (ii) a tobacco‐like pattern (Nicotiana tabacum, N. benthamiana, C. elegans) with a strong H3K27me3 signal. Our data suggest that epigenetic modifications of plant telomere‐associated histones are related neither to the sequence of the telomere motif nor to the lengths of the telomeres. Nor the phylogenetic position of the species plays the role; representatives of the Solanaceae family are included in both groups. As both patterns of histone marks are compatible with fully functional telomeres in respective plants, we conclude that the described specific differences in histone marks are not critical for telomere functions.