Quantitative trait loci mapping in hybrids between Dent and Flint maize multiparental populations reveals group-specific QTL for silage quality traits with variable pleiotropic effects on yield.

KEY MESSAGE: Silage quality traits of maize hybrids between the Dent and Flint heterotic groups mostly involved QTL specific of each parental group, some of them showing unfavorable pleiotropic effects on yield. Maize (Zea mays L.) is commonly used as silage for cattle feeding in Northern Europe. In addition to biomass production, improving whole-plant digestibility is a major breeding objective. To identify loci involved in the general (GCA, parental values) and specific combining ability (SCA, cross-specific value) components of hybrid value, we analyzed an incomplete factorial design of 951 hybrids obtained by crossing inbred lines issued from two multiparental connected populations, each specific to one of the heterotic groups used for silage in Europe ("Dent" and "Flint"). Inbred lines were genotyped for approximately 20K single nucleotide polymorphisms, and hybrids were phenotyped in eight environments for seven silage quality traits measured by near-infrared spectroscopy, biomass yield and precocity (partly analyzed in a previous study). We estimated variance components for GCA and SCA and their interaction with environment. We performed QTL detection using different models adapted to this hybrid population. Strong family effects and a predominance of GCA components compared to SCA were found for all traits. In total, 230 QTL were detected, with only two showing SCA effects significant at the whole-genome level. More than 80% of GCA QTL were specific of one heterotic group. QTL explained individually less than 5% of the phenotypic variance. QTL co-localizations and correlation between QTL effects of quality and productivity traits suggest at least partial pleiotropic effects. This work opens new prospects for improving maize hybrid performances for both biomass productivity and quality accounting for complementarities between heterotic groups.

A glycosyl transferase family 43 protein involved in xylan biosynthesis is associated with straw digestibility in Brachypodium distachyon.

The recalcitrance of secondary plant cell walls to digestion constrains biomass use for the production of sustainable bioproducts and for animal feed. We screened a population of Brachypodium recombinant inbred lines (RILs) for cell wall digestibility using commercial cellulases and detected a quantitative trait locus (QTL) associated with this trait. Examination of the chromosomal region associated with this QTL revealed a candidate gene that encodes a putative glycosyl transferase family (GT) 43 protein, orthologue of IRX14 in Arabidopsis, and hence predicted to be involved in the biosynthesis of xylan. Arabinoxylans form the major matrix polysaccharides in cell walls of grasses, such as Brachypodium. The parental lines of the RIL population carry alternative nonsynonymous polymorphisms in the BdGT43A gene, which were inherited in the RIL progeny in a manner compatible with a causative role in the variation in straw digestibility. In order to validate the implied role of our candidate gene in affecting straw digestibility, we used RNA interference to lower the expression levels of the BdGT43A gene in Brachypodium. The biomass of the silenced lines showed higher digestibility supporting a causative role of the BdGT43A gene, suggesting that it might form a good target for improving straw digestibility in crops.

Natural diversity in daily rhythms of gene expression contributes to phenotypic variation.

Daily rhythms of gene expression provide a benefit to most organisms by ensuring that biological processes are activated at the optimal time of day. Although temporal patterns of expression control plant traits of agricultural importance, how natural genetic variation modifies these patterns during the day and how precisely these patterns influence phenotypes is poorly understood. The circadian clock regulates the timing of gene expression, and natural variation in circadian rhythms has been described, but circadian rhythms are measured in artificial continuous conditions that do not reflect the complexity of biologically relevant day/night cycles. By studying transcriptional rhythms of the evening-expressed gene gigantea (GI) at high temporal resolution and during day/night cycles, we show that natural variation in the timing of GI expression occurs mostly under long days in 77 Arabidopsis accessions. This variation is explained by natural alleles that alter light sensitivity of GI, specifically in the evening, and that act at least partly independent of circadian rhythms. Natural alleles induce precise changes in the temporal waveform of GI expression, and these changes have detectable effects on phytochrome interacting factor 4 expression and growth. Our findings provide a paradigm for how natural alleles act within day/night cycles to precisely modify temporal gene expression waveforms and cause phenotypic diversity. Such alleles could confer an advantage by adjusting the activity of temporally regulated processes without severely disrupting the circadian system.

Root tip contact with low-phosphate media reprograms plant root architecture.

Plant roots are able to sense soil nutrient availability. In order to acquire heterogeneously distributed water and minerals, they optimize their root architecture. One poorly understood plant response to soil phosphate (P(i)) deficiency is a reduction in primary root growth with an increase in the number and length of lateral roots. Here we show that physical contact of the Arabidopsis thaliana primary root tip with low-P(i) medium is necessary and sufficient to arrest root growth. We further show that loss-of-function mutations in Low Phosphate Root1 (LPR1) and its close paralog LPR2 strongly reduce this inhibition. LPR1 was previously mapped as a major quantitative trait locus (QTL); the molecular origin of this QTL is explained by the differential allelic expression of LPR1 in the root cap. These results provide strong evidence for the involvement of the root cap in sensing nutrient deficiency, responding to it, or both. LPR1 and LPR2 encode multicopper oxidases (MCOs), highlighting the essential role of MCOs for plant development.

Natural variation for anthocyanin accumulation under high-light and low-temperature stress is attributable to the ENHANCER OF AG-4 2 (HUA2) locus in combination with PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1) and PAP2.

Growing conditions combining high light intensities and low temperatures lead to anthocyanin accumulation in plants. This response was contrasted between two Arabidopsis thaliana accessions, which were used to decipher the genetic and molecular bases underlying the variation of this response. Quantitative trait loci (QTLs) for flowering time (FT) and anthocyanin accumulation under a high-light and low-temperature scenario versus a control environment were mapped. Major QTLs were confirmed using near-isogenic lines. Candidate genes were examined using mutants and gene expression studies as well as transgenic complementation. Several QTLs were found for FT and for anthocyanin content, of which one QTL co-located at the ENHANCER OF AG-4 2 (HUA2) locus. That HUA2 is a regulator of both pathways was confirmed by the analysis of loss-of-function mutants. For a strong expression of anthocyanin, additional allelic variation was detected for the PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1) and PAP2 genes which control the anthocyanin pathway. The genetic control of variation for anthocyanin content was dissected in A. thaliana and shown to be affected by a common regulator of flowering and anthocyanin biosynthesis together with anthocyanin-specific regulators.

Conserved histidine of metal transporter AtNRAMP1 is crucial for optimal plant growth under manganese deficiency at chilling temperatures.

Manganese (Mn) is an essential nutrient required for plant growth, in particular in the process of photosynthesis. Plant performance is influenced by various environmental stresses including contrasting temperatures, light or nutrient deficiencies. The molecular responses of plants exposed to such stress factors in combination are largely unknown. Screening of 108 Arabidopsis thaliana (Arabidopsis) accessions for reduced photosynthetic performance at chilling temperatures was performed and one accession (Hog) was isolated. Using genetic and molecular approaches, the molecular basis of this particular response to temperature (G × E interaction) was identified. Hog showed an induction of a severe leaf chlorosis and impaired growth after transfer to lower temperatures. We demonstrated that this response was dependent on the nutrient content of the soil. Genetic mapping and complementation identified NRAMP1 as the causal gene. Chlorotic phenotype was associated with a histidine to tyrosine (H239Y) substitution in the allele of Hog NRAMP1. This led to lethality when Hog seedlings were directly grown at 4°C. Chemical complementation and hydroponic culture experiments showed that Mn deficiency was the major cause of this G × E interaction. For the first time, the NRAMP-specific highly conserved histidine was shown to be crucial for plant performance.

Flagellin perception varies quantitatively in Arabidopsis thaliana and its relatives.

Much is known about the evolution of plant immunity components directed against specific pathogen strains: They show pervasive functional variation and have the potential to coevolve with pathogen populations. However, plants are effectively protected against most microbes by generalist immunity components that detect conserved pathogen-associated molecular patterns (PAMPs) and control the onset of PAMP-triggered immunity. In Arabidopsis thaliana, the receptor kinase flagellin sensing 2 (FLS2) confers recognition of bacterial flagellin (flg22) and activates a manifold defense response. To decipher the evolution of this system, we performed functional assays across a large set of A. thaliana genotypes and Brassicaceae relatives. We reveal extensive variation in flg22 perception, most of which results from changes in protein abundance. The observed variation correlates with both the severity of elicited defense responses and bacterial proliferation. We analyzed nucleotide variation segregating at FLS2 in A. thaliana and detected a pattern of variation suggestive of the rapid fixation of a novel adaptive allele. However, our study also shows that evolution at the receptor locus alone does not explain the evolution of flagellin perception; instead, components common to pathways downstream of PAMP perception likely contribute to the observed quantitative variation. Within and among close relatives, PAMP perception evolves quantitatively, which contrasts with the changes in recognition typically associated with the evolution of R genes.

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana.

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment-specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2-5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.

Targeted linkage map densification to improve cell wall related QTL detection and interpretation in maize.

Several QTLs for cell wall degradability and lignin content were previously detected in the F288 × F271 maize RIL progeny, including a set of major QTLs located in bin 6.06. Unexpectedly, allelic sequencing of genes located around the bin 6.06 QTL positions revealed a monomorphous region, suggesting that these QTLs were likely "ghost" QTLs. Refining the positions of all QTLs detected in this population was thus considered, based on a linkage map densification in most important QTL regions, and in several large still unmarked regions. Re-analysis of data with an improved genetic map (173 markers instead of 108) showed that ghost QTLs located in bin 6.06 were then fractionated over two QTL positions located upstream and downstream of the monomorphic region. The area located upstream of bin 6.06 position carried the major QTLs, which explained from 37 to 59 % of the phenotypic variation for per se values and extended on only 6 cM, corresponding to a physical distance of 2.2 Mbp. Among the 92 genes present in the corresponding area of the B73 maize reference genome, nine could putatively be considered as involved in the formation of the secondary cell wall [bHLH, FKBP, laccase, fasciclin, zinc finger C2H2-type and C3HC4-type (two genes), NF-YB, and WRKY]. In addition, based on the currently improved genetic map, eight QTLs were detected in bin 4.09, while only one QTL was highlighted in the initial investigation. Moreover, significant epistatic interaction effects were shown for all traits between these QTLs located in bin 4.09 and the major QTLs located in bin 6.05. Three genes related to secondary cell wall assembly (ZmMYB42, COV1-like, PAL-like) underlay QTL support intervals in this newly identified bin 4.09 region. The current investigations, even if they were based only on one RIL progeny, illustrated the interest of a targeted marker mapping on a genetic map to improve QTL position.

Natural variation for seed dormancy in Arabidopsis is regulated by additive genetic and molecular pathways.

Timing of germination is presumably under strong natural selection as it determines the environmental conditions in which a plant germinates and initiates its postembryonic life cycle. To investigate how seed dormancy is controlled, quantitative trait loci (QTL) analyses has been performed in six Arabidopsis thaliana recombinant inbred line populations by analyzing them simultaneously using a mixed model QTL approach. The recombinant inbred line populations were derived from crosses between the reference accession Landsberg erecta (Ler) and accessions from different world regions. In total, 11 delay of germination (DOG) QTL have been identified, and nine of them have been confirmed by near isogenic lines (NILs). The absence of strong epistatic interactions between the different DOG loci suggests that they affect dormancy mainly by distinct genetic pathways. This was confirmed by analyzing the transcriptome of freshly harvested dry seeds of five different DOG NILs. All five DOG NILs showed discernible and different expression patterns compared with the expression of their genetic background Ler. The genes identified in the different DOG NILs represent largely different gene ontology profiles. It is proposed that natural variation for seed dormancy in Arabidopsis is mainly controlled by different additive genetic and molecular pathways rather than epistatic interactions, indicating the involvement of several independent pathways.

Precise control of water stress in the field reveals different response thresholds for forage yield and digestibility of maize hybrids.

Introduction: With dwindling global freshwater supplies and increasing water stress, agriculture is coming under increasing pressure to reduce water use. Plant breeding requires high analytical capabilities. For this reason, near-infrared spectroscopy (NIRS) has been used to develop prediction equations for whole-plant samples, particularly for predicting dry matter digestibility, which has a major impact on the energy value of forage maize hybrids and is required for inclusion in the official French catalogue. Although the historical NIRS equations have long been used routinely in seed company breeding programmes, they do not predict all variables with the same accuracy. In addition, little is known about how accurate their predictions are under different water stress-environments. Methods: Here, we examined the effects of water stress and stress intensity on agronomic, biochemical, and NIRS predictive values in a set of 13 modern S0-S1 forage maize hybrids under four different environmental conditions resulting from the combination of a northern and southern location and two monitored water stress levels in the south. Results: First, we compared the reliability of NIRS predictions for basic forage quality traits obtained using the historical NIRS predictive equations and the new equations we recently developed. We found that NIRS predicted values were affected to varying degrees by environmental conditions. We also showed that forage yield gradually decreased as a function of water stress, whereas both dry matter and cell wall digestibilities increased regardless of the intensity of water stress, with variability among the tested varieties decreasing under the most stressed conditions. Discussion: By combining forage yield and dry matter digestibility, we were able to quantify digestible yield and identify varieties with different strategies for coping with water stress, raising the exciting possibility that important potential selection targets still exist. Finally, from a farmer's perspective, we were able to show that late silage harvest has no effect on dry matter digestibility and that moderate water stress does not necessarily result in a loss of digestible yield.

Incremental steps toward incompatibility revealed by Arabidopsis epistatic interactions modulating salicylic acid pathway activation.

Plant growth is influenced by genetic factors and environmental cues. Genotype-by-environment interactions are governed by complex genetic epistatic networks that are subject to natural selection. Here we describe a novel epistatic interaction modulating growth in response to temperature common to 2 Arabidopsis recombinant inbred line (RIL) populations (Ler x Kas-2 and Ler x Kond). At 14 degrees C, lines with specific allele combinations at interacting loci (incompatible interactions) have severe growth defects. These lines exhibit deregulated cell death programs and enhanced disease resistance. At 20 degrees C, growth defects are suppressed, but a positive trait of enhanced resistance is retained. Mapping of 1 interacting QTL to a cluster of RPP1-like TIR-NB-LRR genes on chromosome 3 is consistent with our finding that environmentally conditioned epistasis depends on activation of the salicylic acid (SA) stress signaling pathway. The nature of the epistatic interaction conforms to the Dobzhansky-Muller model of genetic incompatibility with incomplete penetrance for reproductive isolation. Variation in fitness of different incompatible lines reveals the presence of additional modifiers in the genetic background. We propose that certain interacting loci lead to an optimal balance between growth and resistance to pathogens by modulating SA signaling under specific environments. This could allow the accumulation of additional incompatibilities before reaching complete reproductive isolation.

Relationships between growth, growth response to nutrient supply, and ion content using a recombinant inbred line population in Arabidopsis.

Growth is an integrative trait that responds to environmental factors and is crucial for plant fitness. A major environmental factor influencing plant growth is nutrient supply. In order to explore this relationship further, we quantified growth-related traits, ion content, and other biochemical traits (protein, hexose, and chlorophyll contents) of a recombinant inbred line population of Arabidopsis (Arabidopsis thaliana) grown on different levels of potassium and phosphate. Performing an all subsets multiple regression analyses revealed a link between growth-related traits and mineral nutrient content. Based on our results, up to 85% of growth variation can be explained by variation in ion content, highlighting the importance of ionomics for a broader understanding of plant growth. In addition, quantitative trait loci (QTLs) were detected for growth-related traits, ion content, further biochemical traits, and their responses to reduced supplies of potassium or phosphate. Colocalization of these QTLs is explored, and candidate genes are discussed. A QTL for rosette weight response to reduced potassium supply was identified on the bottom of chromosome 5, and its effects were validated using selected near isogenic lines. These lines retained over 20% more rosette weight in reduced potassium supply, accompanied by an increase in potassium content in their leaves.

Distinct patterns of genetic variation alter flowering responses of Arabidopsis accessions to different daylengths.

Many plants flower in response to seasonal changes in daylength. This response often varies between accessions of a single species. We studied the variation in photoperiod response found in the model species Arabidopsis (Arabidopsis thaliana). Seventy-two accessions were grown under six daylengths varying in 2-h intervals from 6 to 16 h. The typical response was sigmoidal, so that plants flowered early under days longer than 14 h, late under days shorter than 10 h, and at intermediate times under 12-h days. However, many accessions diverged from this pattern and were clustered into groups showing related phenotypes. Thirty-one mutants and transgenic lines were also scored under the same conditions. Statistical comparisons demonstrated that some accessions show stronger responses to different daylengths than are found among the mutants. Genetic analysis of two such accessions demonstrated that different quantitative trait loci conferred an enhanced response to shortening the daylength from 16 to 14 h. Our data illustrate the spectrum of daylength response phenotypes present in accessions of Arabidopsis and demonstrate that similar phenotypic variation in photoperiodic response can be conferred by different combinations of loci.

Responses of Maize Internode to Water Deficit Are Different at the Biochemical and Histological Levels.

Maize feeding value is strongly linked to plant digestibility. Cell wall composition and structure can partly explain cell wall digestibility variations, and we recently showed that tissue lignification and lignin spatial distribution also contribute to cell wall digestibility variations. Although the genetic determinism of digestibility and cell wall composition has been studied for more than 20 years, little is available concerning that of tissue lignification. Moreover, maize yield is negatively impacted by water deficit, and we newly highlighted the impact of water deficit on cell wall digestibility and composition together with tissue lignification. Consequently, the aim of this study was to explore the genetic mechanisms of lignin distribution in link with cell wall composition and digestibility under contrasted water regimes. Maize internodes from a recombinant inbred line (RIL) population grown in field trials with contrasting irrigation scenarios were biochemically and histologically quantified. Results obtained showed that biochemical and histological traits have different response thresholds to water deficit. Histological profiles were therefore only modified under pronounced water deficit, while most of the biochemical traits responded whatever the strength of the water deficit. Three main clusters of quantitative trait locus (QTL) for histological traits were detected. Interestingly, overlap between the biochemical and histological clusters is rare, and one noted especially colocalizations between histological QTL/clusters and QTL for p-coumaric acid content. These findings reinforce the suspected role of tissue p-coumaroylation for both the agronomic properties of plants as well as their digestibility.

A locus conferring resistance to Colletotrichum higginsianum is shared by four geographically distinct Arabidopsis accessions.

Colletotrichum higginsianum is a hemibiotrophic fungal pathogen that causes anthracnose disease on Arabidopsis and other crucifer hosts. By exploiting natural variation in Arabidopsis we identified a resistance locus that is shared by four geographically distinct accessions (Ws-0, Kondara, Gifu-2 and Can-0). A combination of quantitative trait loci (QTL) and Mendelian mapping positioned this locus within the major recognition gene complex MRC-J on chromosome 5 containing the Toll-interleukin-1 receptor/nucleotide-binding site/leucine-rich repeat (TIR-NB-LRR) genes RPS4 and RRS1 that confer dual resistance to C. higginsianum in Ws-0 (Narusaka et al., 2009). We find that the resistance shared by these diverse Arabidopsis accessions is expressed at an early stage of fungal invasion, at the level of appressorial penetration and establishment of intracellular biotrophic hyphae, and that this determines disease progression. Resistance is not associated with host hypersensitive cell death, an oxidative burst or callose deposition in epidermal cells but requires the defense regulator EDS1, highlighting new functions of TIR-NB-LRR genes and EDS1 in limiting early establishment of fungal biotrophy. While the Arabidopsis accession Ler-0 is fully susceptible to C. higginsianum infection, Col-0 displays intermediate resistance that also maps to MRC-J. By analysis of null mutants of RPS4 and RRS1 in Col-0 we show that these genes, individually, do not contribute strongly to C. higginsianum resistance but are both required for resistance to Pseudomonas syringae bacteria expressing the Type III effector, AvrRps4. We conclude that distinct allelic forms of RPS4 and RRS1 probably cooperate to confer resistance to different pathogens.

Polyamines: molecules with regulatory functions in plant abiotic stress tolerance.

Early studies on plant polyamine research pointed to their involvement in responses to different environmental stresses. During the last few years, genetic, transcriptomic and metabolomic approaches have unravelled key functions of different polyamines in the regulation of abiotic stress tolerance. Nevertheless, the precise molecular mechanism(s) by which polyamines control plant responses to stress stimuli are largely unknown. Recent studies indicate that polyamine signalling is involved in direct interactions with different metabolic routes and intricate hormonal cross-talks. Here we discuss the integration of polyamines with other metabolic pathways by focusing on molecular mechanisms of their action in abiotic stress tolerance. Recent advances in the cross talk between polyamines and abscisic acid are discussed and integrated with processes of reactive oxygen species (ROS) signalling, generation of nitric oxide, modulation of ion channel activities and Ca(2+) homeostasis, amongst others.

Keep on growing under drought: genetic and developmental bases of the response of rosette area using a recombinant inbred line population.

Variation in leaf development caused by water deficit was analysed in 120 recombinant inbred lines derived from two Arabidopsis thaliana accessions, Ler and An-1. Main effect quantitative trait loci (QTLs) and QTLs in epistatic interactions were mapped for the responses of rosette area, leaf number and leaf 6 area to water deficit. An epistatic interaction between two QTLs affected the response of whole rosette area and individual leaf area but only with effects in well-watered condition. A second epistatic interaction between two QTLs controlled the response of rosette area and leaf number with specific effects in the water deficit condition. These effects were validated by generating and phenotyping new appropriate lines. Accordingly, a low reduction of rosette area was observed for lines with a specific allelic combination at the two interacting QTLs. This low reduction was accompanied by an increase in leaf number with a lengthening of the vegetative phase and a low reduction in individual leaf area with low reductions in epidermal cell area and number. Statistical analyses suggested that responses of epidermal cell area and number to water deficit in individual leaves were partly caused by delay in flowering time and reduction in leaf emergence rate, respectively.

Relation among plant growth, carbohydrates and flowering time in the Arabidopsis Landsberg erecta x Kondara recombinant inbred line population.

Arabidopsis thaliana natural variation was used to study plant performance viewed as the accumulation of photo-assimilates, their allocation and storage, in relation to other growth-related features and flowering-related traits. Quantitative trait locus (QTL) analysis using recombinant inbred lines derived from the cross between Landsberg erecta (originating from Poland) and Kondara (originating from Tajikistan) grown on hydroponics, revealed QTLs for the different aspects of plant growth-related traits, sugar and starch contents and flowering-related traits. Co-locations of QTLs for these different aspects were detected at different regions, mainly at the ER locus; the top of chromosomes 3, 4 and 5; and the bottom of chromosome 5. Increased plant growth was associated with early flowering and leaf transitory starch, and correlated negatively with the levels of soluble sugar at early phases of development. From the significant correlations and the co-locations of the QTLs for these aspects, we conclude that there is a complex relationship between plant growth-related traits, carbohydrate content and flowering-related traits.

Dataset of organic sample near infrared spectra acquired on different spectrometers.

This dataset presents 127 raw near infrared spectra of different organic samples acquired on three different spectrometers in three different labs. An example of data processing is shown to create six spectra transfer models between the three spectrometers (two by two). In order to build and validate these transfer models, the dataset was split into two sets of spectra: a first set was used to compute six spectra transfer models thanks to the Piecewise Direct standardisation function (PDS). A second set of spectra, independent of the first one was used to validate transfer models. Spectrum treatments and models were created on ChemFlow (https://vm-chemflow-francegrille.eu/), a free online chemometric software that includes all the necessary functions.