- Research article
- Open Access
Root morphology and seed and leaf ionomic traits in a Brassica napus L. diversity panel show wide phenotypic variation and are characteristic of crop habit
BMC Plant Biology volume 16, Article number: 214 (2016)
Mineral nutrient uptake and utilisation by plants are controlled by many traits relating to root morphology, ion transport, sequestration and translocation. The aims of this study were to determine the phenotypic diversity in root morphology and leaf and seed mineral composition of a polyploid crop species, Brassica napus L., and how these traits relate to crop habit. Traits were quantified in a diversity panel of up to 387 genotypes: 163 winter, 127 spring, and seven semiwinter oilseed rape (OSR) habits, 35 swede, 15 winter fodder, and 40 exotic/unspecified habits. Root traits of 14 d old seedlings were measured in a ‘pouch and wick’ system (n = ~24 replicates per genotype). The mineral composition of 3–6 rosette-stage leaves, and mature seeds, was determined on compost-grown plants from a designed experiment (n = 5) by inductively coupled plasma-mass spectrometry (ICP-MS).
Seed size explained a large proportion of the variation in root length. Winter OSR and fodder habits had longer primary and lateral roots than spring OSR habits, with generally lower mineral concentrations. A comparison of the ratios of elements in leaf and seed parts revealed differences in translocation processes between crop habits, including those likely to be associated with crop-selection for OSR seeds with lower sulphur-containing glucosinolates. Combining root, leaf and seed traits in a discriminant analysis provided the most accurate characterisation of crop habit, illustrating the interdependence of plant tissues.
High-throughput morphological and composition phenotyping reveals complex interrelationships between mineral acquisition and accumulation linked to genetic control within and between crop types (habits) in B. napus. Despite its recent genetic ancestry (<10 ky), root morphology, and leaf and seed composition traits could potentially be used in crop improvement, if suitable markers can be identified and if these correspond with suitable agronomy and quality traits.
Plants require at least 14 essential mineral elements to complete their life-cycles . These include nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg) and sulphur (S), which are macronutrients required in large amounts (typically 1000– > 10,000 mg kg−1 leaf dry weight, DW). The micronutrients chlorine (Cl), boron (B), iron (Fe), manganese (Mn), zinc (Zn), copper (Cu), nickel (Ni) and molybdenum (Mo) are required in smaller amounts (typically 0.1–100 mg kg−1 leaf DW) . Plants also accumulate non-essential elements, some of which have little or no effect on plant growth and development at the concentrations they occur in nature, and others of which may have beneficial and/or detrimental effects depending upon their concentrations in plant tissues. These include arsenic (As), cadmium (Cd), selenium (Se), silicon (Si) and sodium (Na).
Most mineral elements are taken up in ionic form from the soil solution by plant roots. Traits/phenes affecting root morphology and anatomy play a key role in the acquisition of mineral nutrients by plants and impact on crop yields [3–5]. For example, increased root hairs and shallower basal root growth angles can increase P uptake [6, 7]. Reduced allocation of carbon to root structures via increased aerenchyma and reduced cortical cell file formations  and smaller root diameter  may allow some plants more efficient access to larger soil volumes, and thereby water and nutrients. The subsequent uptake and utilisation of mineral elements by plants is controlled by traits affecting ion transport, translocation and sequestration . Mineral elements in both chelated and free-ionic forms move across the root via apoplastic (extracellular) and symplastic (intracellular) pathways to the stele. Following xylem loading and subsequent transport to transpiring leaf tissues, elements are taken up from the leaf apoplast by specific cell types. Translocation of mineral elements in the plants to non-transpiring or xylem-deficient tissues occurs via the phloem [10, 11]. Some elements are highly mobile in phloem tissues (K, Na, Mg, Cd, N, P, S, Se and Cl), some are relatively immobile in the phloem (Ca and Mn), and some elements have intermediate mobility (B, Fe, Zn, Cu, Mo and I) [10–12].
The term ‘ionome’ defines the complement of mineral elements in all of their chemical forms within an organism or tissue, irrespective of whether they are essential or non-essential . The ionome is thus the inorganic subset of the metabolome at a given moment in space and time, which varies at all scales. Within an individual plant, an ionome is specific to tissue type and developmental stage; e.g. seed, fruit and tuber concentrations of Ca are lower than leaf concentrations of Ca due to its limited phloem mobility [14, 15]. Between individuals, the ionome of a specific tissue type varies due to environmental and genetic factors at all scales and this can be observed as differences between populations, species, and plant families [13, 14, 16–18].
Variation in the ionomes of edible crop tissues has enabled identification of quantitative trait loci (QTL) linked to mineral composition and important to human and animal nutrition in several crop species [3, 19, 20]. For example, genetic loci affecting the mineral composition of leaves of Brassica oleracea , Brassica rapa , Brassica napus  and Lotus japonicus  have been identified. In the study of Bus et al. , there were strong pair-wise positive correlations in the shoot concentrations of many of the 11 mineral elements in 30 d old B. napus (>500 genotypes). Furthermore, there were many pair-wise negative correlations between the shoot concentrations of several elements, notably Ca and K, and numerous leaf and seedling size related traits. Plant ionomes are also amenable to genetic dissection using natural and induced genetic variation via mutagenesis, using association mapping and reverse genetic approaches. Several genes underlying variation in mineral nutrient acquisition and translocation have recently been identified. For example, in Arabidopsis thaliana, a deletion mutant with a reduced leaf Ca concentration led subsequently to the identification of ESB1 (Enhanced Suberin Biosynthesis 1) which affects Casparian Band formation [25, 26]. A mutant with reduced leaf Mg, Ca, Fe, and Mo and increased leaf Na and K concentration was similarly associated with reduced sphingolipid biosynthesis . A variety of other Arabidopsis genes are associated with phenotypic variation in leaf As , Cd , K , S and Se .
Brassica napus is an important crop in global terms, with crop types including oilseed rape (OSR), vegetable swede, and fodder crops. Currently, oilseed types of OSR are the third largest source of vegetable oil globally after soybean and oil palm. Worldwide production of OSR was 72.8 Mt in 2013 . Other uses for OSR oils include biodiesel and rape meal for animal feeds, and co-products, including vitamin E (tocopherol) and cholesterol lowering compounds (phytosterols) from the oil, and waxes from pod walls with medical/cosmetic properties. Further industrial oils are currently underexploited but could increase economic margins for farmers. There is considerable scope for improvement of yield of seeds and co-products if suitable traits can be identified and introduced into well-adapted varieties, for example, through improvements in yield and resource-use efficiency [33, 34]. Worldwide average yields for OSR have increased from 1.5 to 2 t ha−1 from 2000 to 2013. Yields are higher in Western Europe, with 2013 average yields of 3.5 t ha−1. The long term average yield of UK OSR is 3.1 t ha−1 , which is much less than UK wheat (8.1 t ha−1) and UK barley (6.4 t ha−1) yet it is similarly nutrient-intensive . The yields of UK OSR are also far less than their estimated potential of >6.5 t ha−1 .
The aim of this study was to determine the phenotypic diversity in root morphology, shoot ionomic (leaf and seed) and seed size/yield traits within a broad genetic diversity panel of B. napus (encompassing all crop types) and to identify their relationship to crop habit. Determining the phenotypic diversity in these traits, and their interrelationships, in this population could inform subsequent studies to dissect the genetic bases and identify markers in traits relevant for crop improvement . An increased understanding of these traits could also help in breeding strategies via more conventional means. To our knowledge, no previous studies have simultaneously characterised the phenotypic variation in root morphology, ionomes and seed size from such a large diversity panel, which is likely to capture most of the species-wide variation in these traits in B. napus.
Plant material for all experiments
Inbred lines of Brassica napus L. genotypes were used in this study. These were from the ERANET-ASSYST consortium diversity population [23, 38–40]. A core panel of 387 genotypes were selected, comprising 163 winter, 127 spring, and seven semiwinter oilseed rape (OSR), 35 swede, 15 winter fodder, and 40 exotic/unspecified habits (Additional file 1: Table S1). Two cultivation systems were deployed. Seedling root traits were determined in a ‘pouch and wick’ hydroponic system in a controlled environment (CE) room. Leaf and seed mineral composition traits were measured on compost-grown plants grown in a designed experiment in a polytunnel.
Root phenotyping in a pouch and wick system
The ‘pouch and wick’ high-throughput phenotyping (HTP) system was reported previously [5, 41]. This system comprised growth pouches assembled from blue germination paper (SD7640; Anchor Paper Company, St Paul, MN, USA), re-cut to 24 × 30 cm and overlain with black polythene (Cransford Polythene Ltd, Woodbridge, UK). Along one of the shorter edges, the paper and polythene were clipped together to an acrylic rod (Acrylic Online, Hull, UK) using ‘bulldog’-type fold-back clips. The growth pouches were suspended above plastic drip trays, supported within lightweight aluminium/polycarbonate frames (KJN Aluminium Profiles, Leicester, UK). Each drip tray contained 2 L of 25 % strength Hoagland’s solution (No. 2 Basal Salt Mixture, Sigma Aldrich, Dorset, UK) made with deionised water. Drip trays were replenished with 500 mL of deionised water every 3 d. Prior to sowing, the pouches were suspended above the nutrient solution for a minimum of 4 h to become fully saturated. Within each aluminium frame, nine drip trays were used, arranged in three columns and three rows. Pouches were allocated randomly to drip trays, 10 or 11 pouches per drip tray, thus 96 pouches and 192 plants per frame (i.e. a single plant on each side of the paper). A total of four frames were used in each experimental run, giving a potential sample size of 768 plants per run within the CE room. The CE room was 2.2 m width, 3.3 m length, 3.0 m height, set to a 12 h photoperiod 18/15 °C day/night temperatures and relative humidity of 60–80 %. Photosynthetically Active Radiation (PAR; measured at plant height with a 190 SB quantum sensor; LI-COR Inc., Lincoln, NE, USA) was approximately 207 μmol m−2 s−1, generated by 400 W white fluorescent lamps (HIT 400w/u/Euro/4 K, Venture Lighting, Rickmansworth, UK).
A single seed was sown on each germination paper, in the middle of the upper edge of the paper, by pressing the seed into the paper. The potential effect of seed size on root traits was controlled for by selecting individual seeds which spanned a range of sizes for each genotype, therefore giving a mean seed diameter of ~1.8 mm for each genotype. Seeds of each genotype were sieved using mesh with a diameter (Ø) of 1.4, 1.7, 2.0 and 2.36 mm (Scientific Laboratory Supplies Ltd, Hessle, UK). Seed retained within the mesh of each faction were selected such that 25 % of seed represented each Ø-category for each genotype. Where insufficient seeds were available for a given Ø-category, the next smallest Ø-category was used instead.
Fourteen days after sowing (DAS), the polythene sheets were removed from all pouches and images were taken of the germination paper and root system for downstream image analysis. Images were taken using a Digital Single Lens Reflex (DSLR) camera (Canon EOS 1100D, Canon Inc., Tokyo, Japan) with a focal length of 35 mm at a fixed height of 75 cm. The root images from the HTP system were renamed with each sample’s unique experimental design information using Bulk Rename Utility (Version 220.127.116.11, TGRMN Software, www.bulkrenameutility.co.uk). Images were cropped by reducing extraneous pixels on bulked images, using XnConvert (Version 1.66, www.xnconvert.com). Cropped images were analysed using RootReader2D (RR2D) . First, a ‘batch process’ was carried out which automatically ‘thresholds’, ‘skeletonises’ and ‘builds segments’ of all images in bulk. The root system was then measured on individual images by placing a marker at the base and tip of the primary root. From these markers, RR2D automatically calculates primary root length (PRL), lateral root length (LRL) of all laterals, and lateral root number (LRN). Further traits calculated from these data included total root length (TRL = PRL + LRL), mean lateral root length (MLRL = LRL/LRN) and lateral root density (LRD = LRN/PRL). A database was developed which integrated the experimental design information from the image name, with the RR2D measurements for each sample, using a programming script (2.7.10; Python Software Foundation, www.python.org).
Of the core panel of 387 genotypes, 354 genotypes comprising 156 winter, 124 spring and seven semiwinter OSR, 14 winter fodder, 33 swede and 20 exotic/unspecified types were screened. Two additional reference winter OSR lines were screened in each experimental run. Each experimental run comprised 32 genotypes, of 24 individuals per genotype. There were 16 experimental runs in total. This equates to a total of 11,176 potential images. An image was removed from analysis if the seed had failed to germinate, or if the seed had rolled down the paper and therefore the shoot failed to emerge above the pouch, or if the seedling was stunted with a radicle < 3 cm, or the radicle appeared deformed such as being twisted around the seed. Overall, 29 % of samples were removed from analysis; excluded data are noted in Additional file 1: Table S2 and all images are available on request.
The relative contribution of genotypic and non-genotypic variance components underlying variation in root traits were calculated using a REML (REsidual Maximum Likelihood) procedure according to the model [(run/frame/column/tray/paper-side) + habit + seed size + genotype]. Genotype was subsequently added as a fixed factor to estimate genotype-means of root traits.
Leaf and seed mineral composition traits in soil-grown plants
Growth of plant material
Seed of all genotypes were sown directly into fine-grade (<3 mm particle size) compost-based growing media (Levington Seed & Modular + Sand -F2S; Everris Ltd., Ipswich, UK) in modular propagation trays (650 plants m−2; internal Ø 2.5 cm, module volume 55 cm3; Type ‘104’, Desch Plantpak, Essex, UK). Sowing took place from 22 to 29 October 2013. The compost was covered with perlite and transferred to a glasshouse vented at 15 °C (controlled by TomTech μClimate, Spalding, Lincs). Supplementary, artificial lighting (Philips Master GreenPower SON-T 400 W bulbs controlled by Grasslin Uni 45 timer) was used to maintain day lengths of 12 h light d−1. Watering was once daily by hand as required until transplantation. From 16 to 29 January 2014, five plants of each genotype were transplanted into individual 5 L pots (internal Ø 22.5 cm; height 18 cm) containing Levington C2 compost (Scotts Professional, Ipswich, UK). Pots were arranged within two single-skinned polytunnels (with a Visqueen Luminance Skin, Northern Polytunnels, Colne, UK) with no additional lighting or heating, at the Sutton Bonington Campus of the University of Nottingham (52°49'58.9" N, 1°14'59.2" W).
Pots were arranged in a randomised block design of five replicate blocks using an R script (personal communications, Edmondson RA, superseded ). Three replicates were allocated to one polytunnel, two to the other. Each replicate comprised 432 units, including one of each of the 387 core genotypes, plus 16 reference genotypes added to enable more accurate normalisation. A further 29 genotypes were included to fill gaps. Each replicate block was split into 12 sub-blocks of 36 genotypes, which were allocated at random. Where a lack of germination meant that insufficient plants were available at the transplanting stage, empty, compost-filled pots were used in their place.
Automatic irrigation was controlled in each polytunnel by a Hunter Irrigation Controller (Hunter Industries, San Marcos, CA, USA, provided by Hortech Systems Ltd., Holbeach, UK). Water from a header tank was distributed by a pump (DAB Active JI112M; DAB Pumps Ltd. Bishop's Stortford, UK) such that each pot received 133 mL of water at 08:00, 12:00 and 16:00 each day, via a low density polyethylene (LDPE) pipe based irrigation system fitted with compensated, non-leaking (CNL) drippers at 4 L h−1 capacity. Each CNL dripper supplied four pots using an attached, four-tipped manifold (Netafim, Tel Aviv, Israel, provided by Hortech systems Ltd.). Each system was also fitted with a Dosatron D3GL-2 feed injector (Tresses, France) used to provide plants with Kristalon Red NPK fertiliser (Yara, Grimsby, UK) between 24 March and 22 May 2014. This was set to mix fertiliser from a stock solution (made up at 100 g fertiliser per litre water) into water at a ratio of one part stock solution to 100 parts water before being sent to the pots. Plants were covered by 380 × 900 mm micro-perforated pollination bags (Focus Packaging & Design Ltd, Brigg, UK) once inflorescences began to show to prevent cross pollination. Any side shoots that emerged after bagging were cut. Watering was reduced to 50 % from 2 June 2014 and switched off completely from 1 July 2014 to encourage senescence. All plants were sprayed with 0.1 % (v/v) azoxystrobin (Amistar, Syngenta, Cambridge, UK) to control first signs of Phoma and some Botrytis on 20 November 2013 and were sprayed again on 29 January 2014 and 17 February 2014. Tebuconazole (Folicur, Bayer, Cambridge, UK) and Amistar were applied at a rate of 0.06 % (v/v) on 28 April 2014. Aphid control was by 0.05 % (w/v) Pirimicarb (Aphox, Syngenta) on 20 May 2014 and 0.07 % (v/v) Deltamethrin (Decis, Bayer) on 2 June 2014.
The total quantity of experimental units was 2160. All plants were harvested from the polytunnels from 14 to 17 July 2014. Stems were cut just above the bottom of the micro-perforated bag containing the top of plants. Each bag was then tied up such that no plant material could escape. Labelled bags were placed into 1 m3 ventilated crates for storage prior to threshing. Crates containing plant material were transported to Elsoms Seeds (Spalding, Lincolnshire), where they were threshed for seed with an SRC single plant thresher (Nickerson Brothers Limited, Lincoln) and cleaned using a Selecta seed cleaner (Selecta Machinefabriek B.V., Enkhuizen, Netherlands). Thousand seed weight for each plant was measured using a Contador seed counter (Pfeuffer GmbH, Kitzingen, Germany). Total seed yield per plant are indicative data, since side-stems were removed where these grew outside of the bags.
Sampling, digestion and analysis of leaf samples
Leaves were sampled at the rosette stage (typically 6–8 true leaves showing) from 5 to 11 March 2014. A minimum of three fully expanded leaves were cut from the plant, weighed and photographed while fresh. Leaves from each plant were stored in separate labelled paper bags at -20 °C. All samples were freeze dried (CHRIST Alpha 2-4 LD freeze dryer; Martin Christ Gefriertrocknungsanlagen GmbH, Osterode, Germany) for 48–60 h, and re-weighed. Leaves were homogenised in liquid N2 using a pestle and mortar and kept frozen prior to analyses.
Subsamples (~0.20 g DW) of leaf were digested using a microwave system comprising a Multiwave 3000 platform with a 48-vessel MF50 rotor (Anton Paar GmbH, Graz, Austria); digestion vessels were perfluoroalkoxy (PFA) liner material and polyethylethylketone (PEEK) pressure jackets (Anton Paar GmbH). Leaf material was digested in 2 mL 70 % Trace Analysis Grade HNO3, 1 mL Milli-Q water (18.2 MΩ cm; Fisher Scientific UK Ltd, Loughborough, UK), and 1 mL H2O2 with microwave settings as follows: power = 1400 W, temp = 140 °C, pressure = 2 MPa, time = 45 min. Two operational blanks were included in each digestion run. Duplicate samples of certified reference material (CRM) of leaf (Tomato SRM 1573a, NIST, Gaithersburg, MD, USA) were included approximately every fourth digestion run; laboratory reference material (LRM) from pooled / freeze-dried Brassica napus leaves was also used for later digests. Following digestion, each tube was made up to a final volume of 15 mL by adding 11 mL Milli-Q water and transferred to a 25 mL universal tube (Sarstedt Ltd., Nümbrecht, Germany) and stored at room temperature.
Leaf digestates were diluted 1-in-5 using Milli-Q water prior to elemental analysis. The concentrations of 28 elements were obtained using inductively coupled plasma-mass spectrometry (ICP-MS; Thermo Fisher Scientific iCAPQ, Thermo Fisher Scientific, Bremen, Germany); Ag, Al, As, B, Ba, Ca, Cd, Cr, Co, Cs, Cu, Fe, K, Mg, Mn, Mo, Na, Ni, P, Pb, Rb, S, Se, Sr, Ti, U, V, Zn. Operational modes included: (i) a helium collision-cell (He-cell) with kinetic energy discrimination to remove polyatomic interferences, (ii) standard mode (STD) in which the collision cell was evacuated, and (iii) a hydrogen collision-cell (H2-cell). Samples were introduced from an autosampler incorporating an ASXpress™ rapid uptake module (Cetac ASX-520, Teledyne Technologies Inc., Omaha, NE, USA) through a PEEK nebulizer (Burgener Mira Mist, Mississauga, Burgener Research Inc., Canada). Internal standards were introduced to the sample stream on a separate line via the ASXpress unit and included Sc (20 μg L−1), Rh (10 μg L−1), Ge (10 μg L−1) and Ir (5 μg L−1) in 2 % trace analysis grade HNO3 (Fisher Scientific UK Ltd). External multi-element calibration standards (Claritas-PPT grade CLMS-2; SPEX Certiprep Inc., Metuchen, NJ, USA) included Ag, Al, As, B, Ba, Cd, Ca, Co, Cr, Cs, Cu, Fe, K, Mg, Mn, Mo, Na, Ni, P, Pb, Rb, S, Se, Sr, Ti (semi-quant), U, V and Zn, in the range 0–100 μg L−1 (0, 20, 40, 100 μg L−1). A bespoke external multi-element calibration solution (PlasmaCAL, SCP Science, Courtaboeuf, France) was used to create Ca, K, Mg and Na standards in the range 0–30 mg L−1. Boron, P and S calibration utilized in-house standard solutions (KH2PO4, K2SO4 and H3BO3). In-sample switching was used to measure B and P in STD mode, Se in H2-cell mode and all other elements in He-cell mode. Sample processing was undertaken using Qtegra™ software (Thermo Fisher Scientific) with external cross-calibration between pulse-counting and analogue detector modes when required. In total, 2096 samples were analysed in 14 ICP-MS runs.
Digestion and analysis of seed samples
Dry seeds (three seeds per tube and occasionally four for very small seeds) were transferred into Pyrex test tubes (16 × 100 mm). After weighing an appropriate number of samples the masses of the remaining samples were calculated using method of Danku et al. . The Seed samples were left overnight to pre-digest in 1.16 mL trace metal grade HNO3 (J. T. Baker Instra-Analyzed; Avantor Performance Materials; Scientific & Chemical Supplies Ltd, Aberdeen, UK) spiked with indium internal standard and 1.2 mL H2O2 (Primar-Trace analysis grade, 30 %; Fisher Scientific, Loughborough, UK) was also added. They were then digested in dry block heaters (DigiPREP MS, SCP Science; QMX Laboratories, Essex, UK) at 115 °C for 4 h.
Seed digestates were diluted to 11.5 mL with Milli-Q water (18.2 MΩ cm, Merck Millipore, Watford, UK) and aliquots transferred to 96-well deep well plates using adjustable multichannel pipette (Rainin; Anachem Ltd, Luton, UK) for analysis. Elemental analysis was performed with an ICP-MS (PerkinElmer NexION 300D equipped with Elemental Scientific Inc. autosampler and Apex HF sample introduction system; PerkinElmer LAS Ltd, Seer Green, UK and Elemental Scientific Inc., Omaha, NE, USA, respectively) in the standard mode. Twenty elements (Li, B, Na, Mg, P, S, K, Ca, Mn, Fe, Co, Ni, Cu, Zn, As, Se, Rb, Sr, Mo, and Cd) were monitored. Liquid reference material composed of pooled samples of the digested seed materials was prepared before the first sample run and was used throughout the remaining sample runs. The liquid reference material was included after every ninth sample in all ICP-MS sample sets to correct for variation between and within ICP-MS analysis runs . Sample concentrations were calculated using external calibration method within the instrument software. The calibration standards (with indium internal standard and blanks) were prepared from single element standards (Inorganic Ventures; Essex Scientific Laboratory Supplies Ltd, Essex, UK) solutions. In total, 1986 samples were analysed across four ICP-MS runs.
Data processing of leaf and seed mineral composition traits
For each data-point, an element-specific operational blank concentration (mean of each ICP-MS run) was subtracted. Data were then multiplied by initial sample volume, divided by the initial dry mass of plant material, and converted to mg element kg−1 of dry leaf or seed material. Element-specific limits of detection (LODs) were reported as three times the standard deviation (SD) of the operational blank concentrations, assuming a notional starting dry weight of 0.200 g for leaf and 0.015 g for seed data (Additional file 1: Table S3). For leaves, element-specific recoveries from CRMs ranged from 68 to 134 %, for 18 elements with certified CRM values (Additional file 1: Table S4). From leaf mineral concentration data, seven elements (Ag, Co, Cr, Ni, Pb, U, V) were removed from further analysis due to having mean mineral concentrations which were less than or close to the LOD (Additional file 1: Table S5). Likewise, seven elements (As, Co, Cr, Fe, Ni, Pb, Se) were removed from seed mineral concentration data (Additional file 1: Table S6). For those elements retained for analysis, data for individual leaf and seed element concentrations which were below element-specific LODs were replaced with half LOD values. Leaf and seed element concentrations >5 standard deviation (SDs) above the global arithmetic mean for each element were also removed as a precaution against using contaminated samples (125 out of 58,688 values for leaves; 107 out of 42,504 values for seed).
The relative contribution of genotypic and non-genotypic variance components underlying variation in leaf and seed composition traits was calculated using a REML procedure in GenStat. Genotype, habit and experimental sources of variation were classed as random factors according to the model [habit + genotype + polytunnel + polytunnel/replicate + polytunnel/replicate/sub-block]. For leaf composition traits, a further model was used [replicate + (replicate/sub-block) + genotype + (replicate/genotype)] in which genotype was subsequently added as a fixed factor to estimate genotype-means. For seeds, the arithmetic mean data were used for each genotype.
Multivariate analysis of root morphology and mineral composition traits
Correlation analysis was conducted on all 945 possible pairwise combinations of the 44 root, leaf and seed trait variate sets (genotype means). Five stepwise discriminant analyses were conducted in GenStat, one each for the root morphology-, leaf- and seed mineral- and seed weight variate sets, which contained 6, 21,15 and two variates, respectively, and one for the variate set of all traits combined. Genotypes were grouped according to ‘crop habit’. The Wilks’ Lambda ‘forward selection’ stepwise algorithm option was selected, which, at each step, adds the trait-variate which explains the most between-group variation from all of the remaining trait-variate sets. Specificity plots were drawn, to view the proportion of genotypes of each ‘crop habit’ correctly allocated to each group, at each step. Discrimination plots were drawn to represent the separation of variation in the crop habits in two dimensions. All statistical analyses were conducted using GenStat 15th Edition (VSN International Ltd, Hemel Hempstead, UK).
Results and discussion
Root growth was influenced strongly by seed size
Seed diameter accounted for a large proportion of the variation in total root length (TRL; 44 %), primary root length (PRL; 35 %), lateral root length (LRL; 41 %) and lateral root number (LRN; 41 %), but not for mean lateral root length (MLRL; 6 %) or lateral root density (LRD; 3 %) in 14 d old seedlings (Table 1; Fig. 1a). Genotype/habit factors accounted for between 7 % (MLRL) and 17 % (PRL) of the total variation in the six root traits. Residual (plant-to-plant) variation accounted for the largest single source of variation in the study, up to 75 and 81 % for LRD and MLRL, respectively, indicating that lateral roots traits are particularly responsive to the environment. This large residual source of variation is consistent with previous studies of Brassica seedling root traits, which show that large numbers of individuals are required to detect subtle differences in root traits between genotypes with confidence [5, 45]. Thousand seed weight (TSW) in the 2013 seed, from which all plants were grown, varied significantly within and between crop habits, from largest to smallest in: semiwinter OSR, winter OSR, spring OSR, winter fodder and swede types (P < 0.001, Fig. 1b). However, whilst seed diameter had a significant positive correlation with root length, based on the data for individual seedlings, potential correlations between TSW and root length could not be tested in this study because seeds were selected for uniformity between genotypes based on diameter classification and not by individual seed weights.
Positive relationships have been reported previously between seed size and the seminal root length and total root weight of barley (Hordeum vulgare) , and between seed size and total root size and lateral root number in tomato (Solanum lycopersicum) . Larger seeds have also been shown to improve seedling establishment, shoot weight, biomass and final yield in some, but not all field studies of OSR in Canada [48–51]. However, larger-sized seeds typically had more vigorous early growth . Thousand seed weight was also shown to correlate positively with absolute growth rate 21 days after germination . Improved seed size-related root growth of B. napus seedlings might also increase tolerance to shoot pests (e.g. flea beetle, Phyllotreta spp.)  and root diseases such as Rhizoctonia solani which can damage the primary roots of B. napus . Seed weight has previously been associated with pre-emergence growth in a bi-parental mapping population of Brassica oleracea, but under separate genetic control to germination [54, 55]. Additionally, the present study found that the thousand seed weight (TSW) in the winter OSR varieties from different release periods has steadily increased over time, suggesting that larger seeds may have been bred for (Additional file 2: Figure S6). This present study shows there is scope to exploit the genetic control of seed size-related root growth as a potential route to improve early vigour in the small-seeded B. napus.
Winter OSR and fodder types had larger root systems than other crop habits
Winter OSR and winter fodder types had a greater mean TRL, PRL, TLRL, and LNR at 14 d than the other crop habits (P < 0.001, Fig. 2; Additional file 1: Table S7). Semiwinter OSR had a shorter mean PRL than all other habits (P < 0.001, Fig. 2b), and a greater mean LRD (P < 0.001, Fig. 2f). It is important to note that these differences in root system size between OSR crop types were observed when seeds of uniform diameter were sown for each genotype. Increased root length in seedlings is likely to indicate increased early vigour. Velicka et al.  observed that early sowing afforded a greater root collar thickness and leaf number, and these earlier sown plants had greater over-winter survival and more rapid accumulation of matter in the apical bud in spring. Furthermore, Scott et al.  observed that earlier sowing significantly increased seed yield because of increased leaf and root growth. Seedling root-length traits, measured in this same ‘pouch and wick’ system, correlated with early plant growth and final seed yield in 30 commercial winter OSR B. napus genotypes . Finch-Savage et al.  suggested that vigorous early root growth is essential for small seeded crops such as Brassica to acquire resources before desiccation occurs. A fast-growing, thick root collar contains large amounts of soluble carbohydrates which will enable the plant to withstand frost and afford a rapid re-growth in spring . Thus, sufficient early root growth is necessary for winter crop survival and may have been selected for based on yield in previous breeding programs, whereas spring sown crops have less need for a rapid development to ensure hardiness.
Spring varieties typically had higher leaf concentrations of macronutrients and some micronutrients than winter varieties
The mean leaf concentration of 21 elements varied by more than six orders of magnitude across genotypes, from 0.01 mg kg−1 (As) to >50,000 mg kg−1 (K) (Fig. 3; Additional file 1: Table S7). Genotypic variation in leaf mineral concentration ranged from 1.8-fold (Fe) to >40-fold (Se). Among the macronutrients, leaf mineral concentrations varied 2.0-fold for P, 2.1-fold for K, 3.0-fold for Ca, 2.6-fold for Mg, and 2.5-fold variation in S. In comparison, among a panel of ~450 B. oleracea, also grown in compost and sampled during early vegetative growth, shoot concentrations of: Ca and Mg varied 2.0- and 2.3-fold , respectively; P and K varied 4.9-  and 3.4-fold , respectively. Among a panel of soil-grown 509 inbred lines of B. napus, the shoot mineral concentrations of 30 d old seedlings varied (approximately) 2.0-fold for Ca, 1.6-fold for Mg, 6.7-fold for P and 2.0-fold for K .
Winter OSR, winter fodder and swede had lower mean leaf macronutrient (Ca, Mg, P, K, S) concentrations than spring and semiwinter OSR (Ca, Mg; P < 0.001, Fig. 3). Semiwinter OSR also had higher mean leaf Ca and Mg concentrations compared to other habits (P < 0.001). Among the micronutrients, leaf Cu was greater in spring OSR than other habits (P < 0.001). Leaf Fe concentrations were greater in winter and spring OSR (P < 0.001). The mean leaf Mo concentrations were greatest in Spring OSR and swede (P < 0.001), although there was substantial variation within crop type. The mean leaf concentrations of beneficial and non-essential elements (As, Cd, Na, Se) were consistently higher in the semiwinter OSR leaves used in this study, typically followed by spring OSR (Fig. 3). Likewise, in the study of Bus et al. , winter OSR also had lower mean shoot Ca, K and S concentrations than spring and semiwinter OSR, but similar P concentrations and semiwinter OSR had the highest shoot S and Zn of the crop types. Despite these overall trends in the data, there is wide variation in shoot mineral composition within all crop types of B. napus, which will be influenced strongly by the nutritional environment in which the plant is grown as well as genotypic factors.
Variance components analysis (Table 1; Additional file 1: Table S10) shows that genotype had the largest influence on leaf S concentration (40 %) and the smallest influence on leaf Se concentrations (0 %). Habit accounted for the least variation in all traits (generally less than 10 %) but had the greatest effect on leaf Na, Mg and S concentration. The trends of heritabilities (i.e. genotype effect) for leaf composition traits follow a similar pattern to those observed previously in soil grown leaves of Arabidopsis , whereby leaf Mg was the most heritable macro nutrient in their study, and the second most heritable in this study. Leaf Ca, K and Mo concentration were ranked among the most heritable leaf composition traits in both studies. Leaf Fe, Mn and Cu concentration were among the least heritable traits in both studies. The variance components analysis indicates that the effect of experimental variance is generally higher for micronutrients than macronutrients.
Seed mineral concentrations were consistent across habits for many nutrients, but S concentrations were lower and Mo concentrations were higher in OSR types
The mean seed concentration of 15 elements varied by more than six orders of magnitude across genotypes, from 0.01 mg kg−1 (Cd) to >13,000 mg kg−1 (K) (Fig. 4; Additional file 1: Table S7). Genotypic variation in seed mineral concentration varied 1.7-fold (P) to 14-fold (Na). Among the macronutrients, seed mineral concentrations varied 3.1-fold for Ca, 1.9-fold for Mg, 2.0-fold for K, and 7.5-fold for S. We are not aware of previous reports of species-wide variation in seed mineral composition traits in a Brassica species. White and Broadley  reviewed variation in the mineral composition of edible cereal grains and dicot seeds for several species, typically core germplasm collections, which had been grown under comparative conditions. Among the dicots, seed Ca concentration varied 3.7-, 2.0-, 9.1-, 1.5- and 1.9-fold, and seed Mg concentration varied 2.4-, 1.4-, 2.3-, 1.3-, and 1.6-fold, for chickpea, peanut, pea, bean and soybean, respectively. Therefore, the seed macronutrient composition of this B. napus panel appears to be a similar range as other dicot species.
Winter and spring OSR had similar seed macronutrient concentrations, except for P, in which spring and semiwinter OSR had higher seed P concentrations than winter types (P < 0.001, Fig. 4) and Mg in which the semiwinter had higher concentrations than other types (P < 0.001). Winter fodder and swede had higher seed S concentrations than OSR crop types (P < 0.001), presumably due to the smaller proportion of “double-low” (low glucosinolate, low erucic acid) compared to “double-high” varieties, because winter fodder and swede have not been selected for low seed glucosinolate concentration. Glucosinolates are sulphur and nitrogen-containing secondary metabolites common in the Brassicaceae family [62, 63], and some genes of the sulphate assimilation pathway are members of the glucosinolate biosynthetic network . Seed Mo concentrations showed the opposite pattern, with higher concentrations in OSR than other types (P < 0.001). An antagonistic effect of sulphate on Mo concentration in the shoot, root and seeds of OSR has previously been observed [65, 66]. Given that S and Mo are known to share some assimilation and transport pathways [67, 68], this could imply a relationship between seed glucosinolates, S and Mo. It has previously been suggested that tissue specific demand can regulate the expression of sulphate transporters . It could be surmised that plants with low seed S/glucosinolates display a slightly different sulphate transporter profile leading to differences in Mo accumulation. Nevertheless, how S and Mo are interacting in relation to glucosinolates is yet to be determined. Consistent with S/Mo transport being under strong genetic control, variance components analysis (Table 1; Additional file 1: Table S10) showed that genotype had the largest influence on seed Mo concentration (41 %) out of all the elements analysed in this study. As observed previously in Arabidopsis seeds , Mo, Cu, S and K concentrations were influenced by genotype to a greater extent than seed Mg and P concentrations. Of the micronutrients, seed B, Cu and Mn concentrations were generally consistent across habits; B was highest in spring OSR (P < 0.001), Cu was highest in winter fodder and winter OSR (P < 0.001), Mn was highest in semiwinter OSR (P < 0.001). Zn concentrations were the highest in winter fodder (P < 0.001), Cd concentrations were highest in semiwinter OSR (P < 0.001), and seed Na concentrations were highest in spring OSR.
Differences in nutrient translocation between crop habits can be detected by comparing leaf and seed concentration ratios of elements with potentially similar transport or assimilation pathways
Four pairs of elements were selected, for which both leaf and seed concentration data were available, and for which there are reports of shared transport/assimilation pathways. These element pairs are S and Mo [67, 68], Ca and Sr [70, 71], K and Rb [72–74], and Zn and Cd . The second element of each of these pairs is either not a nutrient, or, for Mo, is only required in very small amounts compared to the first element in each pair. The first hypothesis tested was that the ratio of first:second element concentration in leaves is ≫ 1:1, as expected from external nutrient supply and plant requirements. The second hypothesis was that the seed ratio/leaf ratio of elements is >1:1. For example, in the case of S and Mo, this would be expressed as ([S]seed/[Mo]seed) / ([S]leaf/[Mo]leaf; Fig. 5a). A seed ratio/leaf ratio of elements >1:1 indicates that the net seed accumulation of the essential, or more abundant, nutrient element, is greater than the non-essential, or less abundant, nutrient. This could be due to selective processes in the source (e.g. increased mobilisation from the leaf) and/or sink (e.g. decreased remobilisation from the pod) tissues. As expected, the ratio of first:second element is ≫ 1:1 in leaves. For S:Mo, the mean of 387 genotypes is 15902:1 (range 4808–37661:1); for Ca:Sr 519:1 (433–604:1); for K:Rb 3106:1 (2652–3760:1); for Zn:Cd 200:1 (85–449:1). However, there was considerable variation in the seed/leaf ratio of elements (Fig. 5; Additional file 1: Table S7). For S:Mo, the differences between habits was most apparent. For the OSR habits, most OSR genotypes had seed ratio/leaf ratios <1:1, indicating net accumulation of Mo in seeds is greater than S in relative terms (Fig. 5c). In contrast, most swede genotypes and almost half of the winter fodder genotypes had a greater net accumulation of S than Mo. This relationship is indicative of a potential link between glucosinolates, S and Mo content of seeds but the nature of the interactions between these components has yet to be elucidated and is the focus of current studies on this population. For Ca:Sr, most genotypes showed slight increase in net accumulation of Ca compared to Sr in seeds, with swede habits having the highest mean (Fig. 5b). Interestingly, all genotypes showed increased net accumulation of Rb in seeds compared to K, with swede habits showing the lowest relative net accumulation of K compared to Rb (Fig. 5c). The reasons for this observation are not clear. For Zn:Cd, all genotypes showed increased net Zn accumulation in seeds compared to Cd. This was again greater in swede habits, and also winter fodder; both had greater net accumulation of Zn in seeds, compared to Cd, than did the OSR habits.
Root traits and leaf and seed mineral composition traits correlate within, but not between tissues
There were strong correlations between root traits. The strongest positive correlations were between traits relating directly to the total length of the root system (i.e. comprising PRL and LRL components, P < 0.001) and also between these traits and LRN (P < 0.001, Fig. 6; Additional file 1: Table S8). There was a weak negative relationship between PRL and LRD (r = -0.28; P < 0.001), since the latter is derived from LRN/PRL. In addition, there was also a weak negative relationship between MLRL and LRD (r = -0.38; P < 0.001), suggesting a trade-off between lateral root length and number. These correlations are consistent with previous observations on a much smaller panel, of 32 UK-field adapted OSR genotypes, grown in the same ‘pouch and wick’ system . There were many strong positive correlations between the leaf concentrations of pairs of elements (Fig. 6; Additional file 1: Table S8). The strongest positive relationships were between leaf concentrations of Ca and Sr (r = 0.97), Sr and Ba (r = 0.93), K and Rb (r = 0.92) and Ca and Mg (r = 0.87); all P < 0.001. Positive correlations between the leaf concentrations of Group II elements reflect the relative lack of selectivity between these elements during transport within the plant . Such relationships between Ca and Mg have been observed previously within-species, including among panels of diverse B. oleracea , OSR , Arabidopsis  and multi-species datasets [17, 73]. The strong positive relationship in the leaf concentration of Group I elements K and Rb is as expected from previous observation across many species . However, the transport of other Group I elements is typically much more selective than Group II elements . There were few negative relationships between the leaf concentrations of pairs of elements (Fig. 6), but weak negative relationships were observed between Rb and the Group II elements, Ba (r = -0.24), Sr (r = -0.18), and Ca (r = -0.16); all P < 0.001.
There were few strong positive correlations between the seed concentrations of pairs of elements (Fig. 6; Additional file 1: Table S8). The strongest positive correlation was between Ca and Sr (r = 0.51; P < 0.001), followed by between K and Rb (r = 0.40; P < 0.001). Both of these correlations were much weaker than those observed in leaves. The strongest negative correlations between the seed concentrations of pairs of elements were between S and Mo (r = -0.47), S and B (r = -0.46), and Ca and K (r = -0.40); all P < 0.001 (Fig. 6).
In general, leaf and seed mineral composition traits correlated very weakly (Fig. 6; Additional file 1: Table S8). The strongest positive correlation in a compositional trait between plant parts was a positive correlation between leaf Cd and seed Cd (r = 0.41; P < 0.001). All other correlations between leaf and seed mineral composition traits, and between root traits and leaf or seed composition traits were weaker, with correlation coefficients ranging from -0.26 to +0.33. Likewise, with the exception of P, no correlations were observed between elemental concentrations of leaf, root and seed tissues in Arabidopsis .
The strongest positive correlations between a seedling root trait and a leaf composition trait were between LRD and leaf Ca (r = 0.16; P = 0.006), Sr (r = 0.15; P = 0.008) and Ba (r = 0.14; P = 0.01) concentrations. A similar weak positive correlation between seedling LRD and leaf Ca (and Zn) concentrations was also seen previously in some field experiments . Interestingly, PRL had weak, but significant, negative correlations with most leaf composition traits e.g. Mo (r = 0.31, P = 0.001) and Na (r = 0.30, P = 0.001). These data provide some evidence that LRD might be a beneficial trait for nutrient resource and above-ground biomass acquisition in B. napus . Some leaf and seed mineral concentrations e.g. Ca and Mg and some beneficial elements (Fig. 3) of semiwinter OSR genotypes were greater than in other habits (Fig. 2f). Semiwinter OSR genotypes used as starting material for this study had the greatest LRD and greatest mean TSW of the five crop habits, which may have led to improved overall root size and function, and mineral acquisition in this study (Additional file 2: Figure S5). Semiwinter types are likely to have a distinct breeding history from winter and spring OSR habits due to having more introgressions from B. rapa and a longer period of domestication [76, 77]. These differences in pedigree might also explain some of the variation in root and shoot traits assigned in this present study to crop habit, and warrant further study in a range of environments.
Combining root, leaf and seed traits in a discriminant analysis is characteristic of crop habit
Combining all traits within the discriminant analysis provided the most accurate characterisation of crop habit (Fig. 7). Using the variate set of all traits combined, and after the addition of the 20 most informative trait variate sets, genotypes with the winter (86 %), spring (85 %) and semiwinter (81 %) OSR habits were correctly allocated to the correct group, as were 76 % of swede and 68 % of the winter fodder habits. The trait which contributed the most in terms of allocation to crop habit was leaf S concentration, followed by PRL, and seed Ca and Mo concentration, followed by the TSW of the 2013 seed from which plants were grown. The next most important root trait after PRL was MLRL, which was ranked 12th. The relative contributions of all 44 traits in this variate set to the discriminant analysis are presented in Additional file 1: Table S9. The discriminant analyses for each of the root-, leaf- and seed-trait variate sets was less accurate (Additional file 2: Figures S1–S4).
Brassica napus has been shown to be amenable to rapid marker identification linked to useful agronomic traits in wide diversity populations. For example, using associative transcriptomics (AT) with a panel of 84 B. napus genotypes from within the same panel used here, Harper et al.  showed that seed contents of both erucic acid and glucosinolates (GS) were associated with specific genes known to be involved in their biosynthetic pathways, whilst identifying additional new target loci of potential use in breeding. Similarly, Koprivova et al.  used a subset of this study’s population to identify novel loci linked to shoot anion accumulation. The AT technique is based on transcriptome-sequencing, combined with association mapping. It uses transcribed sequences (mRNA-seq) which allows variation of gene sequences to be detected (through single nucleotide polymorphisms; SNPs) whilst reducing the complexity of the analysis compared to typical genome-wide association analysis (GWAS). In addition, transcript abundance (gene expression markers; GEMs) can be quantified simultaneously. Transcript abundance is likely to be of particular relevance in the control of traits in complex polyploid species in which gene duplication may lead to unequal expression of gene paralogues .
A substantial proportion of the same inbred population used in this study have also been used to study other traits, including those relating to nutritional composition and seedling growth, in other environments [23, 38–40, 80] and those based on leaf ionome traits [25–31]. The most accurate characterisation of crop habit is when multiple plant part traits are combined in analysis. Thus, combining multiple datasets using this panel has the potential for more accurate candidate trait identification, and to dissect the transport pathways which lead to altered elemental accumulation for crop improvement. The volume and throughput of data obtained from root phenotyping and ionomics platforms is considerable and the challenge now is to combine datasets for a better understanding of the ionome and the key traits involved in elemental accumulation, and to mine for the underlying molecular mechanisms of these useful traits using associative transcriptomics analysis.
Lateral root density
Lateral root length
Lateral root number
Mean lateral root length
Primary root length
Total root length
Thousand seed weight
White PJ, Brown PH. Plant nutrition for sustainable development and global health. Ann Bot. 2010;105:1073–80.
Kirkby E. Introduction, Definition and Classification of Nutrients. In: Marschner P, editor. Marschner’s Mineral Nutrition of Higher Plants. 3rd ed. London: Academic; 2012. p. 3–5.
White PJ, Broadley MR, Gregory PJ. Managing the nutrition of plants and people. Appl Environ Soil Sci 2012. doi:10.1155/2012/104826.
White PJ, George TS, Gregory PJ, Bengough AG, Hallett PD, McKenzie BM. Matching roots to their environment. Ann Bot. 2013;112:207–22.
Thomas CL, Graham NS, Hayden R, Meacham MC, Neugebauer K, Nightingale M, et al. High throughput phenotyping (HTP) identifies seedling root traits linked to variation in seed yield and nutrient capture in field-grown oilseed rape (Brassica napus L.). Ann Bot. 2016;in press.
Bonser AM, Lynch J, Snapp S. Effect of phosphorus deficiency on growth angle of basal roots in Phaseolus vulgaris. New Phytol. 1996;132:281–8.
Miguel MA, Widrig A, Vieira RF, Brown KM, Lynch JP. Basal root whorl number: a modulator of phosphorus acquisition in common bean (Phaseolus vulgaris). Ann Bot. 2013;112:973–82.
Lynch JP. Root phenes for enhanced soil exploration and phosphorus acquisition: tools for future crops. Plant Physiol. 2011;156:1041–9.
Hufnagel B, Sousa SM, Assis L, Guimaraes CT, Leiser W, Azevedo GC, et al. Duplicate and conquer: Multiple homologs of PHOSPHORUS-STARVATION TOLERANCE1 enhance phosphorus acquisition and sorghum performance on low-phosphorus soils. Plant Physiol. 2014;166:659–77.
Karley AJ, White PJ. Moving cationic minerals to edible tissues: potassium, magnesium, calcium. Curr Opin Plant Biol. 2009;12:291–8.
White PJ. Long-distance transport in the xylem and phloem. In: Marschner P, editor. Marschner’s Mineral Nutrition of Higher Plants. 3rd ed. London: Academic; 2012. p. 49–70.
Brown PH, Bellaloui N, Wimmer MA, Bassil ES, Ruiz J, Hu H, et al. Boron in plant biology. Plant Biol. 2002;4:205–23.
Salt DE, Baxter I, Lahner B. Ionomics and the study of the plant ionome. Annu Rev Plant Biol. 2008;59:709–33.
White PJ, Broadley MR. Biofortification of crops with seven mineral elements often lacking in human diets – iron, zinc, copper, calcium, magnesium, selenium and iodine. New Phytol. 2009;182:49–84.
Broadley MR, White PJ. Eats roots and leaves. Can edible horticultural crops address dietary calcium (Ca), magnesium (Mg) and potassium (K) deficiencies in humans? Proc Nutr Soc. 2010;69:601–12.
Lahner B, Gong J, Mahmoudian M, Smith EL, Abid KB, Rogers EE, et al. Genomic scale profiling of nutrient and trace elements in Arabidopsis thaliana. Nature Biotechnol. 2003;21:1215–21.
Broadley MR, Bowen HC, Cotterill HL, Hammond JP, Meacham MC, Mead A, et al. Phylogenetic variation in the shoot mineral concentration of angiosperms. J Exp Bot. 2004;55:321–36.
Baxter I. Should we treat the ionome as a combination of individual elements, or should we be deriving novel combined traits? J Exp Bot. 2015;66:2127–31.
Baxter IR, Ziegler G, Lahner B, Mickelbart MV, Foley R, Danku J, et al. Single-kernel ionomic profiles are highly heritable indicators of genetic and environmental influences on elemental accumulation in maize grain (Zea mays). PLOS One. 2014. doi:10.1371/journal.pone.0087628.
Zhang M, Pinson SR, Tarpley L, Huang XY, Lahner B, Yakubova E, et al. Mapping and validation of quantitative trait loci associated with concentrations of 16 elements in unmilled rice grain. Theor Appl Genet. 2014;127:137–65.
Broadley MR, Hammond JP, King GJ, Astley D, Bowen HC, Meacham MC, et al. Shoot calcium (Ca) and magnesium (Mg) concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea. Plant Physiol. 2008;146:1707–20.
Graham NS, Hammond JP, Lysenko A, Mayes S, Ó Lochlainn S, Blasco B, et al. Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs. Plant Cell. 2014;26:2818–30.
Bus A, Körber N, Parkin IAP, Samans B, Snowdon RJ, Li J, et al. Species- and genome-wide dissection of the shoot ionome in Brassica napus and its relationship to seedling development. Front Plant Sci. 2014;5:485.
Chen Z, Watanabe T, Shinano T, Osaki M. Rapid characterization of plant mutants with an altered ion-profile: a case study using Lotus japonicus. New Phytol. 2009;181:795–801.
Baxter I, Hosmani PS, Rus A, Lahner B, Borevitz JO, Muthukumar B, et al. Root suberin forms an extracellular barrier that affects water relations and mineral nutrition in Arabidopsis. PLOS Genet. 2009. doi:10.1371/journal.pgen.1000492.
Hosmani PS, Kamiya T, Danku J, Naseer S, Geldner N, Guerinot ML, et al. Dirigent domain-containing protein is part of the machinery required for formation of the lignin-based Casparian strip in the root. Proc Natl Acad Sci U S A. 2013;110:14498–503.
Chao D-Y, Gable K, Chen M, Baxter I, Dietrich CR, Cahoon EB, et al. Sphingolipids in the root play an important role in regulating the leaf ionome in Arabidopsis thaliana. Plant Cell. 2011;23:1061–81.
Chao D-Y, Chen Y, Chen J, Shi S, Chen Z, Wang C, et al. 2013a. Genome-wide association mapping identifies a new arsenate reductase enzyme critical for limiting arsenic accumulation in plants. PLOS Biol. 2013. doi:10.1371/journal.pbio.1002009.
Chao D-Y, Silva A, Baxter I, Huang YS, Nordborg M, Danku J, et al. Genome-wide association studies identify Heavy Metal ATPase3 as the primary determinant of natural variation in leaf cadmium in Arabidopsis thaliana. PLOS Genet. 2012. doi:10.1371/journal.pgen.1002923.
Chao D-Y, Dilkes B, Luo H, Douglas A, Yakubova E, Lahner B, et al. Polyploids exhibit higher potassium uptake and salinity tolerance in Arabidopsis. Science. 2013;341:658–9.
Chao D-Y, Baraniecka P, Danku J, Koprivova A, Lahner B, Luo H, et al. Variation in sulfur and selenium accumulation is controlled by naturally occurring isoforms of the key sulfur assimilation enzyme ADENOSINE 5’-PHOSPHOSULFATE REDUCTASE2 across the Arabidopsis species range. Plant Physiol. 2014;166:1593–608.
Faostat. http://faostat.fao.org/. Accessed 04 Apr 2016.
Berry PM, Spink J, Foulkes MJ, White PJ. The physiological basis of genotypic differences in nitrogen use efficiency in oilseed rape (Brassica napus L.). Field Crops Res. 2010;119:365–73.
Koeslin-Findeklee F, Becker MA, van der Graaff E, Roitsch T, Horst WJ. Differences between winter oilseed rape (Brassica napus L.) cultivars in nitrogen starvation-induced leaf senescence are governed by leaf-inherent rather than root-derived signals. J Exp Bot. 2015;66:3669–81.
HGCA, Berry P. In: Cook S, Ellis S, Gladders P, Roques S, Boys E, editors. Oilseed rape guide. 2nd ed. Warwickshire: HGCA Publications; 2014.
Defra. Fertiliser Manual (RB209). 8th ed. Norwich: TSO (The Stationery Office); 2010.
Harper AL, Trick M, Higgins J, Fraser F, Clissold L, Wells R, et al. Associative transcriptomics of traits in the polyploid crop species Brassica napus. Nature Biotechnol. 2012;30:798–802.
Bus A, Körber N, Snowdon RJ, Stich B. Patterns of molecular variation in a species-wide germplasm set of Brassica napus. Theor Appl Genet. 2011;123:1413–23.
Körber N, Wittkop B, Bus A, Friedt W, Snowdon RJ, Stich B. Seedling development in a Brassica napus diversity set and its relationship to agronomic performance. Theor Appl Genet. 2012;125:1275–87.
Körber N, Bus A, Li J, Higgins J, Bancroft I, Higgins EE, et al. Seedling development traits in Brassica napus examined by gene expression analysis and association mapping. BMC Plant Biol. 2015;15:136.
Atkinson JA, Wingen LU, Griffiths M, Pound MP, Gaju O, Foulkes MJ, et al. Phenotyping pipeline reveals major seedling root growth QTL in hexaploid wheat. J Exp Bot. 2015;66:2283–92.
Clark RT, Famoso AN, Zhao K, Shaff JE, Craft EJ, Bustamante CD, et al. High-throughput two-dimensional root system phenotyping platform facilitates genetic analysis of root growth and development. Plant Cell Environ. 2013;36:454–66.
Edmondson RN. blocksdesign: Nested Block Designs for Unstructured Treatments. R package version 1.4. http://CRAN.R-project.org/package=blocksdesign. 2015.
Danku JMC, Lahner B, Yakubova E, Salt DE. Large-scale plant ionomics. In: Maathuis FJM, editor. Plant Mineral Nutrients: Methods and Protocols. New York: Humana Press; 2013. p. 255–76.
Adu MO, Chatot A, Wiesel L, Bennett MJ, Broadley MR, White PJ, et al. A scanner system for high-resolution quantification of variation in root growth dynamics of Brassica rapa genotypes. J Exp Bot. 2014;65:2039–48.
Bertholdsson N-O, Brantestam AK. A century of Nordic barley breeding- effects on early vigour root and shoot growth, straw length, harvest index and grain weight. Eur J Agron. 2009;30:266–74.
Khan N, Kazmi RH, Willems LAJ, van Heusden AW, Ligterink W, Hilhorst HWM. Exploring the natural variation for seedling traits and their link with seed dimensions in tomato. PLOS One. 2012. doi:10.1371/journal.pone.0043991.
Elliott RH, Franke C, Rakow GFW. Effects of seed size and seed weight on seedling establishment, vigour and tolerance of Argentine canola (Brassica napus) to flea beetles, Phyllotreta spp. Can J Plant Sci. 2008;88:207–17.
Hwang SF, Ahmed HU, Turnbull GD, Gossen BD, Strelkov SE. The effect of seed size, seed treatment, seeding date and depth on Rhizoctonia seedling blight of canola. Canadian J Plant Sci. 2014;94:311–21.
Hwang SF, Ahmed HU, Turnbull GD, Gossen BD, Strelkov SE. Effect of seeding date and depth, seed size and fungicide treatment on Fusarium and Pythium seedling blight of canola. Can J Plant Sci. 2015;95:293–301.
Harker KN, O'Donovan JT, Smith EG, Johnson EN, Peng G, Willenborg CJ, et al. Seed size and seeding rate effects on canola emergence, development, yield and seed weight. Can J Plant Sci. 2015;95:1–8.
Franzaring J, Holz I, Fangmeier A. Effects of increased autumn temperatures and sub-acute levels of ammonia on post-winter development of four cultivars of winter oilseed rape (Brassica napus L.). J Appl Bot Food Qual. 2012;85:134–43.
Sturrock CJ, Woodhall J, Brown M, Walker C, Mooney SJ, Ray RV. Effects of damping-off caused by Rhizoctonia solani anastomosis group 2-1 on roots of wheat and oil seed rape quantified using X-ray Computed Tomography and real-time PCR. Front Plant Sci. 2015;6:461.
Bettey M, Finch-Savage WE, King GJ, Lynn JR. Quantitative genetic analysis of seed vigour and pre-emergence seedling growth traits in Brassica oleracea. New Phytol. 2000;148:277–86.
Finch-Savage WE, Clay HA, Lynn JR, Morris K. Towards a genetic understanding of seed vigour in small-seeded crops using natural variation in Brassica oleracea. Plant Sci. 2010;179:582–9.
Velicka R, Pupaliene R, Butkeviciene LM, Kriauciuniene Z. Peculiarities of overwintering of hybrid and conventional cultivars of winter oilseed rape depending on the sowing date. Acta Sci Pol. 2012;11:53–66.
Scott RK, Ogunremi EA, Ivins JD, Mendham NJ. The effect of sowing date and season on growth and yield of oilseed rape (Brassica napus). J Agr Sci. 1973;81:277–85.
Lääniste P, Jõudu J, Eremeev V, Mäeorg E. Sowing date influence on winter oilseed rape overwintering in Estonia. Acta Agric Scand Sect B. 2007;57:342–8.
Hammond JP, Broadley MR, White PJ, King GJ, Bowen HC, Hayden R, et al. Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits. J Exp Bot. 2009;60:1953–68.
White PJ, Hammond JP, King GJ, Bowen HC, Hayden RM, Meacham MC, et al. Genetic analysis of potassium use efficiency in Brassica oleracea L. Ann Bot. 2010;105:1199–210.
Baxter I, Hermans C, Lahner B, Yakubova E, Tikhonova M, Verbruggen N, et al. Biodiversity of mineral nutrient and trace element accumulation in Arabidopsis thaliana. PLoS One. 2012. doi:10.1371/journal.pone.0035121.
Bloem E, Haneklaus S, Schnug E. Relation between total sulphur analysed by ICP-AES and glucosinolates in oilseed rape and Indian mustard seeds. Landbauforsch Volk. 2005;55:205–10.
Redovnikovic IR, Glivetic T, Delonga K, Vorkapic-Furac J. Glucosinolates and their potential role in plant. Period Biol. 2008;110:297–309.
Yatusevich R, Mugford SG, Matthewman C, Gigolashvili T, Frerigmann H, Delaney S, et al. Genes of primary sulfate assimilation are part of the glucosinolate biosynthetic network in Arabidopsis thaliana. Plant J. 2010;62:1–11.
Philips CJC, Chiy PC. Effects of applying sodium and sulphur fertilisers on the concentrations of elements in water leached from permanent pasture. J Sci Food Agr. 2002;82:806–15.
Balik J, Pavlikova D, Tlustos P, Sykora K, Cerny J. The fluctuation of molybdenum content in oilseed rape plants after the application of nitrogen and sulphur fertilisers. Plant Soil Environ. 2006;52:301–7.
Malagoli M, Schiavon M, Wirtz M, Hell R, Pilon-Smits EAH, Trevisan S, et al. Molybdate affects sulfate acquisition in Brassica juncea plants. In: DeKok LJ, Tausz M, Hawkesford MJ, Hoefgen R, McManus MT, Norton RM, Rennenberg H, Saito K, Schnug E, Tabe L, editors. Sulfur metabolism in plants: mechanisms and applications to food security and responses to climate change. Dordrecht: Springer; 2012. p. 169–73.
DeTar RA, Alford ÉR, Pilon-Smits EAH. Molybdenum accumulation, tolerance and molybdenum-selenium-sulfur interactions in Astragalus selenium hyperaccumulator and nonaccumulator species. J Plant Physiol. 2015;183:32–40.
Buchner P, Takahashi H, Hawkesford MJ. Plant sulphate transporters: co-ordination of uptake, intracellular and long-distance transport. J Exp B. 2004;55:1765–73.
White PJ. The pathways of calcium movement to the xylem. J Exp Bot. 2001;52:891–9.
Broadley MR, White PJ. Some elements are more equal than others: soil-to-plant transfer of radiocaesium and radiostrontium, revisited. Plant Soil. 2012;355:23–7.
White PJ, Clarkson DT, Earnshaw MJ. Acclimation of potassium influx in rye (Secale cereale) to low root temperatures. Planta. 1987;171:377–85.
Watanabe T, Broadley MR, Jansen S, White PJ, Takada J, Satake K, et al. Evolutionary control of leaf element composition in plants. New Phytol. 2007;174:516–23.
McGonigle TP, Grant CA. Variation in potassium and calcium uptake with time and root depth. Can J Plant Sci. 2015;95:771–7.
Martinka M, Vaculik M, Lux A. Plant cell responses to cadmium and zinc. In: Nick P, Opatrny Z, editors. Applied Plant Cell Biology: Cellular Tools And Approaches For Plant Biotechnology. Berlin: Springer; 2014. p. 209–46.
Qian W, Meng J, Li M, Frauen M, Sass O, Noack J, et al. Introgression of genomic components from Chinese Brassica rapa contributes to widening the genetic diversity in rapeseed (B. napus L.), with emphasis on the evolution of Chinese rapeseed. Theor Appl Genet. 2006;113:49–54.
Chalhoub B, Denoeud F, Liu S, Parkin IAP, Tang H, Wang X, et al. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014;345:950–3.
Koprivova A, Harper AL, Trick M, Bancroft I, Kopriva S. Dissection of the control of anion homeostasis by associative transcriptomics in Brassica napus. Plant Physiol. 2014;166:442–50.
Adams KL, Cronn R, Percifield R, Wendel JF. Genes duplicated by polyploidy show unequal contributions to the transcriptome and organ-specific reciprocal silencing. Proc Natl Acad Sci U S A. 2003;100:4649–54.
Hatzig SV, Frisch M, Breuer F, Nesi N, Ducournau S, Wagner M-H, et al. Genome-wide association mapping unravels the genetic control of seed germination and vigor in Brassica napus. Front Plant Sci. 2015;6:221.
We thank Jonathan Atkinson, Darren Wells, Andrew French and Michael Pound at the University of Nottingham for advice on root measurement methods.
The work was supported by the Biotechnology and Biological Sciences Research Council (BBSRC) Crop Improvement Research Club (CIRC) Grant BB/J019631/1 to MRB, BBSRC grant BB/L000113/1 to DES, the BBSRC Renewable Industrial Products from Rapeseed (RIPR) Programme Grant BB/L002124/1 to IB, and the Rural and Environment Science and Analytical Services Division (RESAS) of the Scottish Government through Work Package 3.3 (2011–2016) to PJW.
Availability of data and materials
The data sets supporting the results of this article are included within the article and its additional files. Data are also available at the Brassica Information Portal (BIP; https://bip.earlham.ac.uk/; The Earlham Institute, Norwich, UK), doi:10.5281/zenodo.59927; doi:10.5281/zenodo.59937; doi:10.5281/zenodo.59936.
CLT collected and analysed root and seed morphology trait data. TDA analysed leaf and seed mineral and discrimination data. RH and SM oversaw growth of plant material in polytunnels. LW and SDY oversaw ICP-MS analysis of leaf tissue. DES and JMCD oversaw ICP-MS analysis of seed tissue. AS advised on effects of S and Mo in plants. MRB and NSG oversaw project management and data analysis. All other authors contributed to project conception and advised on data analysis. All authors have read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Supplementary tables. A collection of tables containing a variety of extra data including raw data gathered in the experiments, limits of detection calculations, data used for analyses of traits, and data used for generation of some of the figures. Also included is a detailed ouput from variance components analysis and a list of the abbreviations used across the tables. (XLSX 1324 kb)
Supplementary figures. A collection of extra figures which may be of interest to readers but that aren’t in the main scope of the submission. Figures S1–S4 show step-wise discriminant analyses plots using different subsets of the traits measured (root morphology traits, leaf mineral composition traits, seed mineral composition traits, & seed yield traits respectively). Plots from analyses using the full set of traits are included as a main figure in the submission; see Fig. 7. Figures S5 and S6 are box and whisker plots of seed yield data by crop habit and thousand seed weight by genotype release date respectively. (PPTX 429 kb)
About this article
Cite this article
Thomas, C.L., Alcock, T.D., Graham, N.S. et al. Root morphology and seed and leaf ionomic traits in a Brassica napus L. diversity panel show wide phenotypic variation and are characteristic of crop habit. BMC Plant Biol 16, 214 (2016). https://doi.org/10.1186/s12870-016-0902-5
- Mineral concentration
- High-throughput phenotyping
- Root morphology
- Seed size
- Leaf/seed elemental ratios