Phylogenetic variation in the shoot mineral concentration of angiosperms

The calcium (Ca) concentration of plant shoot tissues varies systematically between angiosperm orders. The phylogenetic variation in the shoot concentration of other mineral nutrients has not yet been described at an ordinal level. The aims of this study were (1) to quantify the shoot mineral concentration of different angiosperm orders, (2) to partition the phylogenetic variation in shoot mineral concentration between and within orders, (3) to determine if the shoot concentration of different minerals are correlated across angiosperm species, and (4) to compare experimental data with published ecological survey data on 81 species sampled from their natural habitats. Species, selected pro rata from different angiosperm orders, were grown in a hydroponic system under a constant external nutrient regime. Shoots of 117 species were sampled during vegetative growth. Significant variation in shoot carbon (C), calcium (Ca), potassium (K), and magnesium (Mg) concentration occurred between angiosperm orders. There was no evidence for systematic differences in shoot phosphorus (P) or organic-nitrogen (N) concentration between orders. At a species level, there were strong positive correlations between shoot Ca and Mg concentration, between shoot P and organic-N concentration, and between shoot K concentration and shoot fresh weight:dry weight ratio. Shoot C and cation concentration correlated negatively at a species level. Species within the Poales and the Caryophyllales had distinct shoot mineralogies in both the designed experiment and in the ecological survey.

Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea

Calcium (Ca) and magnesium (Mg) are the most abundant group II elements in both plants and animals. Genetic variation in shoot Ca and shoot Mg concentration (shoot Ca and Mg) in plants can be exploited to biofortify food crops and thereby increase dietary Ca and Mg intake for humans and livestock. We present a comprehensive analysis of within-species genetic variation for shoot Ca and Mg, demonstrating that shoot mineral concentration differs significantly between subtaxa (varietas). We established a structured diversity foundation set of 376 accessions to capture a high proportion of species-wide allelic diversity within domesticated Brassica oleracea, including representation of wild relatives (C genome, 1n = 9) from natural populations. These accessions and 74 modern F-1 hybrid cultivars were grown in glasshouse and field environments. Shoot Ca and Mg varied 2- and 2.3-fold, respectively, and was typically not inversely correlated with shoot biomass, within most subtaxa. The closely related capitata (cabbage) and sabauda (Savoy cabbage) subtaxa consistently had the highest mean shoot Ca and Mg. Shoot Ca and Mg in glasshouse-grown plants was highly correlated with data from the field. To understand and dissect the genetic basis of variation in shoot Ca and Mg, we studied homozygous lines from a segregating B. oleracea mapping population. Shoot Ca and Mg was highly heritable (up to 40). Quantitative trait loci (QTL) for shoot Ca and Mg were detected on chromosomes C2, C6, C7, C8, and, in particular, C9, where QTL accounted for 14 to 55 of the total genetic variance. The presence of QTL on C9 was substantiated by scoring recurrent backcross substitution lines, derived from the same parents. This also greatly increased the map resolution, with strong evidence that a 4-cM region on C9 influences shoot Ca. This region corresponds to a 0.41-Mb region on Arabidopsis (Arabidopsis thaliana) chromosome 5 that includes 106 genes. There is also evidence that pleiotropic loci on C8 and C9 affect shoot Ca and Mg. Map-based cloning of these loci will reveal how shoot-level phenotypes relate to Ca 21 and Mg 21 uptake and homeostasis at the molecular level.

Shoot zinc (Zn) concentration varies widely within Brassica oleracea L. and is affected by soil Zn and phosphorus (P) levels

The low availability of zinc (Zn) in soils and crops affects dietary Zn intake worldwide. This study sought to determine if the natural genetic variation in shoot Zn concentrations (Zn(shoot)) is sufficient to pursue a crop improvement breeding strategy in a leafy vegetable crop. The gene-pool of Brassica oleracea L. was sampled using a large (n = 376) diversity foundation set (DFS), representing almost all species-wide common allelic variation, and 74 commercial varieties (mostly F(1)). The DFS genotypes were grown at low and high soil phosphorus (P) levels under glasshouse and field conditions, and also in a Zn-deficient soil, with or without Zn-fertilisation, in a glasshouse. Despite the large variation in Zn(shoot) among genotypes, environment had a profound effect on Zn(shoot) The heritability of Zn(shoot) was significant, but relatively low, among 90 doubled-haploid (DH) lines from a mapping population. While several quantitative trait loci (QTL) associated with Zn(shoot) occurred on chromosomes C2, C3, C5, C7, and C9, these were generally weak and conditional upon growth conditions. Breeding for Zn(shoot) in B. oleracea is therefore likely to be challenging. Shoot P concentrations increased substantially in all genotypes under low soil Zn conditions. Conversely, only some genotypes had increased Zn(shoot) at low soil P levels. Sufficient natural genetic variation may therefore exist to study some of the interactions between Zn and P nutrition.

A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes

Whole-genome transcriptome profiling is revealing how biological systems are regulated at the transcriptional level. This study reports the development of a robust method to profile and compare the transcriptomes of two nonmodel plant species, Thlaspi caerulescens, a zinc (Zn) hyperaccumulator, and Thlaspi arvense, a nonhyperaccumulator, using Affymetrix Arabidopsis thaliana ATH1-121501 GeneChip (R) arrays (Affymetrix, Santa Clara, CA, USA). Transcript abundance was quantified in the shoots of agar- and compost-grown plants of both species. Analyses were optimized using a genomic DNA (gDNA)-based probe-selection strategy based on the hybridization efficiency of Thlaspi gDNA with corresponding A. thaliana probes. In silico alignments of GeneChip (R) probes with Thlaspi gene sequences, and quantitative real-time PCR, confirmed the validity of this approach. Approximately 5000 genes were differentially expressed in the shoots of T. caerulescens compared with T. arvense, including genes involved in Zn transport and compartmentalization. Future functional analyses of genes identified as differentially expressed in the shoots of these closely related species will improve our understanding of the molecular mechanisms of Zn hyperaccumulation.

Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits

The environmental and financial costs of using inorganic phosphate fertilizers to maintain crop yield and quality are high. Breeding crops that acquire and use phosphorus (P) more efficiently could reduce these costs. The variation in shoot P concentration (shoot-P) and various measures of P use efficiency (PUE) were quantified among 355 Brassica oleracea L. accessions, 74 current commercial cultivars, and 90 doubled haploid (DH) mapping lines from a reference genetic mapping population. Accessions were grown at two or more external P concentrations in glasshouse experiments; commercial and DH accessions were also grown in replicated field experiments. Within the substantial species-wide diversity observed for shoot-P and various measures of PUE in B. oleracea, current commercial cultivars have greater PUE than would be expected by chance. This may be a consequence of breeding for increased yield, which is a significant component of most measures of PUE, or early establishment. Root development and architecture correlate with PUE; in particular, lateral root number, length, and growth rate. Significant quantitative trait loci associated with shoot-P and PUE occur on chromosomes C3 and C7. These data provide information to initiate breeding programmes to improve PUE in B. oleracea.

A baseline study of vicine–convicine levels in faba bean (Vicia faba L.) germplasm

Vicine and convicine are anti-nutritional compounds that accumulate in the cotyledons of faba beans. When humans consume beans with high levels of these compounds, it can cause a condition called favism in individuals harbouring a deficiency in the activity of their glucose-6-phosphate dehydrogenase. When faba beans are used in animal feeds, there can be effects on performance. These concerns have resulted in increasing interest within plant breeding in developing low vicine and convicine faba bean germplasm. In order to facilitate this objective, we developed a rapid and robust screening method for vicine and convicine, capable of distinguishing between faba beans that are either high (wild type) or low in vicine and convicine. In the absence of reliable commercial reference materials, we report an adaptation of a previously published method where a biochemical assay and spectral data were used to confirm the identity of our analytes, vicine and convicine. This method could be readily adopted in other facilities and open the way to the efficient exploitation of diverse germplasm in regions where faba beans play a significant role in human nutrition. We screened a collection of germplasm of interest to a collaborative plant breeding programme developing between the National Institute for Agricultural Botany in the UK and L'Institut Nationale d'Agronomie de Tunisie in Tunisia. We report the results obtained and discuss the prospects for developing molecular markers for the low vicine and convicine trait.

Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.)

Key message We have identified QTLs for stomatal characteristics on chromosome II of faba bean by applying SNPs derived from M. truncatula , and have identified candidate genes within these QTLs using synteny between the two species. Abstract Faba bean (Vicia faba L.) is a valuable food and feed crop worldwide, but drought often limits its production, and its genome is large and poorly mapped. No information is available on the effects of genomic regions and genes on drought adaptation characters such as stomatal characteristics in this species, but the synteny between the sequenced model legume, Medicago truncatula, and faba bean can be used to identify candidate genes. A mapping population of 211 F5 recombinant inbred lines (Mélodie/2 × ILB 938/2) were phenotyped to identify quantitative trait loci (QTL) affecting stomatal morphology and function, along with seed weight, under well-watered conditions in a climate-controlled glasshouse in 2013 and 2014. Canopy temperature (CT) was evaluated in 2013 under water-deficit (CTd). In total, 188 polymorphic single nucleotide polymorphisms (SNPs), developed from M. truncatula genome data, were assigned to nine linkage groups that covered ~928 cM of the faba bean genome with an average inter-marker distance of 5.8 cM. 15 putative QTLs were detected, of which eight (affecting stomatal density, length and conductance and CT) co-located on chromosome II, in the vicinity of a possible candidate gene—a receptor-like protein kinase found in the syntenic interval of M. truncatula chromosome IV. A ribose-phosphate pyrophosphokinase from M. truncatula chromosome V, postulated as a possible candidate gene for the QTL for CTd, was found some distance away in the same chromosome. These results demonstrate that genomic information from M. truncatula can successfully be translated to the faba bean genome.

A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.)

Faba bean (Vicia faba L.) is a globally important nitrogen-fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence-based consensus map was used to explore synteny with the most closely-related crop species, lentil, and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2 cM interval of Vf chromosome 2 which was highly collinear with Mt3. The obvious candidate gene from 77 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene (Mt3g092830, Mt3g092840) controlling anthocyanin production in Medicago and re-sequencing of the Vf orthologue showed a putative causative deletion of the entire 5’ end of the gene.

Flanking SNP markers for vicine–convicine concentration in faba bean (Vicia faba L).

The pyrimidine glycosides, vicine and convicine, limit the use of faba bean (Vicia faba L.) as food and feed. A single recessive gene, vc-, is responsible for a lowered vicine–convicine concentration. The biosynthetic pathway of these closely related compounds is not known, and the nearest available markers are several cM away from vc-. Improved markers would assist breeding and help to identify candidate genes. A segregating population of 210 F5 recombinant inbred lines was developed from the cross of Mélodie/2 (low vicine–convicine) × ILB 938/2 (normal vicine–convicine), and vicine–convicine concentrations were determined twice on each line. The population was genotyped with a set of 188 SNPs. A strong, single QTL for vicine–convicine concentration was identified on chromosome I, flanked by markers 1.0 cM away on one side and 2.6 cM on the other. The interval defined by these markers in the model species Medicago truncatula includes about 340 genes, but no candidate genes were identified. Further fine mapping should lead to the identification of tightly linked markers as well as narrowing down the search for candidate regulatory or biosynthetic genes which could underlie the vc- locus.

High-resolution melt and morphological analyses of mealybugs (Hemiptera: Pseudococcidae) from cacao: tools for the control of Cacao Swollen Shoot Virus spread

BACKGROUND: Mealybugs (Hemiptera: Coccoidea: Pseudococcidae) are key vectors of badnaviruses, including Cacao Swollen Shoot Virus (CSSV) the most damaging virus affecting cacao (Theobroma cacao L.). The effectiveness of mealybugs as virus vectors is species dependent and it is therefore vital that CSSV resistance breeding programmes in cacao incorporate accurate mealybug identification. In this work the efficacy of a CO1-based DNA barcoding approach to species identification was evaluated by screening a range of mealybugs collected from cacao in seven countries. RESULTS: Morphologically similar adult females were characterised by scanning electron microscopy and then, following DNA extraction, were screened with CO1 barcoding markers. A high degree of CO1 sequence homology was observed for all 11 individual haplotypes including those accessions from distinct geographical regions. This has allowed for the design of a High Resolution Melt (HRM) assay capable of rapid identification of the commonly encountered mealybug pests of cacao. CONCLUSIONS: HRM Analysis (HRMA) readily differentiated between mealybug pests of cacao that can not necessarily be identified by conventional morphological analysis. This new approach, therefore, has potential to facilitate breeding for resistance to CSSV and other mealybug transmitted diseases.

Field phenotyping of potato to assess root and shoot characteristics associated with drought tolerance

Aims Potatoes are a globally important source of food whose production requires large inputs of fertiliser and water. Recent research has highlighted the importance of the root system in acquiring resources. Here measurements, previously generated by field phenotyping, tested the effect of root size on maintenance of yield under drought (drought tolerance). Methods Twelve potato genotypes, including genotypes with extremes of root size, were grown to maturity in the field under a rain shelter and either irrigated or subjected to drought. Soil moisture, canopy growth, carbon isotope discrimination and final yields were measured. Destructively harvested field phenotype data were used as explanatory variables in a general linear model (GLM) to investigate yield under conditions of drought or irrigation. Results Drought severely affected the small rooted genotype Pentland Dell but not the large rooted genotype Cara. More plantlets, longer and more numerous stolons and stolon roots were associated with drought tolerance. Previously measured carbon isotope discrimination did not correlate with the effect of drought. Conclusions These data suggest that in-field phenotyping can be used to identify useful characteristics when known genotypes are subjected to an environmental stress. Stolon root traits were associated with drought tolerance in potato and could be used to select genotypes with resilience to drought.

Insect pollination reduces yield loss following heat stress in faba bean (Vicia faba L.)

Global food security, particularly crop fertilization and yield production, is threatened by heat waves that are projected to increase in frequency and magnitude with climate change. Effects of heat stress on the fertilization of insect-pollinated plants are not well understood, but experiments conducted primarily in self-pollinated crops, such as wheat, show that transfer of fertile pollen may recover yield following stress. We hypothesized that in the partially pollinator-dependent crop, faba bean (Vicia faba L.), insect pollination would elicit similar yield recovery following heat stress. We exposed potted faba bean plants to heat stress for 5 days during floral development and anthesis. Temperature treatments were representative of heat waves projected in the UK for the period 2021-2050 and onwards. Following temperature treatments, plants were distributed in flight cages and either pollinated by domesticated Bombus terrestris colonies or received no insect pollination. Yield loss due to heat stress at 30°C was greater in plants excluded from pollinators (15%) compared to those with bumblebee pollination (2.5%). Thus, the pollinator dependency of faba bean yield was 16% at control temperatures (18 to 26°C) and extreme stress (34°C), but was 53% following intermediate heat stress at 30°C. These findings provide the first evidence that the pollinator dependency of crops can be modified by heat stress, and suggest that insect pollination may become more important in crop production as the probability of heat waves increases.

Cacao swollen shoot virus in Nigeria: analysis of a pathogen and its vectors

Cacao swollen shoot virus (CSSV) causes the Cacao swollen shoot virus disease (CSSVD) and significantly reduces production in West African cacao. This study characterised the current status of the disease in the major cacao growing States in Nigeria and attempted a clarification on the manner of CSSV transmission. Two separate field surveys and sample collections were conducted in Nigeria in summer 2012 and spring 2013. PCR-based screening of cacao leaf samples and subsequent DNA sequencing showed that the disease continues to persist in Ondo and Oyo States and in new cacao sites in Abia, Akwa Ibom, Cross River and Edo States. Mealybug samples collected were identified using a robust approach involving environmental scanning electron microscopy, histology and DNA barcoding, which highlighted the importance of integrative taxonomy in the study. The results show that the genus Planococcus (Planococcus citri (Risso) and/or Planococcus minor (Maskell)) was the most abundant vector (73.5%) at the sites examined followed by Formicococcus njalensis (Laing) (19.0 %). In a laboratory study, the feeding behaviour of Pl. citri, Pseudococcus longispinus (Targioni-Tozzetti) and Pseudococcus viburni (Signoret) on cacao were investigated using electrical penetration graph (EPG) analysis. EPG waveforms reflecting intercellular stylet penetration (C), extracellular salivation (E1e), salivation in sieve elements (E1), phloem ingestion (E2), derailed stylet mechanics (F), xylem ingestion (G) and non-probing phase (Np) were analysed. Individual mealybugs exhibited marked variation within species and significantly differed (p ≤ .05) between species for E1e and E1. PCR-based assessments of the retention time for CSSV in viruliferous Pl. citri, Ps. longispinus and Ps. viburni fed on a non-cacao diet showed that CSSV was still detectable after 144 hours. These unusually long durations for a pathogen currently classified as a semi-persistent virus have implications for the design of non-malvaceous barrier crops currently being considered for the protection of new cacao plantings.

Susceptibility of faba bean (Vicia faba L.) to heat stress during floral development and anthesis

Experiments were conducted over two years to quantify the response of faba bean (Vicia faba L.) to heat stress. Potted winter faba bean plants (cv. Wizard) were exposed to temperature treatments (18/10; 22/14; 26/18; 30/22; 34/26°C day/night) for five days during floral development and anthesis. Developmental stages of all flowers were scored prior to stress, plants were grown in exclusion from insect pollinators to prevent pollen movement between flowers, and yield was harvested at an individual pod scale, enabling effects of heat stress to be investigated at a high resolution. Susceptibility to stress differed between floral stages, flowers were most affected during initial green-bud stages. Yield and pollen germination of flowers present before stress showed threshold relationships to stress, with lethal temperatures (t50) ~28°C and ~32°C, while whole plant yield showed a linear negative relationship to stress with high plasticity in yield allocation, such that yield lost at lower nodes was partially compensated at higher nodal positions. Faba bean has many beneficial attributes for sustainable modern cropping systems but these results suggest that yield will be limited by projected climate change, necessitating the development of heat tolerant cultivars, or improved resilience by other mechanisms such as earlier flowering times.

A tillering inhibition gene influences root-shoot carbon partitioning and pattern of water use to improve wheat productivity in rainfed environments

Genetic modification of shoot and root morphology has potential to improve water and nutrient 19 uptake of wheat crops in rainfed environments. Near-isogenic lines (NILs) varying for a tillering 20 inhibition (tin) gene and representing multiple genetic backgrounds were investigated in contrasting 21 controlled environments for shoot and root growth. Leaf area, shoot and root biomass were similar 22 until tillering whereupon reduced tillering in tin-containing NILs produced reductions of up to 60% in 23 total leaf area and biomass, and increases in total root length of up to 120% and root biomass to 24 145%. Together, root-to-shoot ratio increased two-fold with the tin gene. The influence of tin on shoot 25 and root growth was greatest in the cv. Banks genetic background, particularly in the biculm-selected 26 NIL, and was typically strongest in cooler environments. A separate de-tillering study confirmed 27 greater root-to-shoot ratios with regular tiller removal in non-tin containing genotypes. In validating 28 these observations in a rainfed field study, the tin allele had a negligible effect on seedling growth but 29 was associated with significantly (P<0.05) reduced tiller number (-37%), leaf area index (-26%) and 30 spike number (-35%) to reduce plant biomass (-19%) at anthesis. Root biomass, root-to-shoot ratio at 31 early stem elongation and root depth at maturity were increased in tin-containing NILs. Soil water use 32 was slowed in tin-containing NILs resulting in greater water availability, greater stomatal 33 conductance, cooler canopy temperatures and maintenance of green leaf area during grain-filling. 34 Together these effects contributed to increases in harvest index and grain yield. In both the controlled 35 and field environments, the tin gene was commonly associated with increased root length and biomass 36 but the significant influence of genetic background and environment suggests careful assessment of 37 tin-containing progeny in selection for genotypic increases in root growth.