A. Mead

High-Resolution Temporal Profiling of Transcripts during Arabidopsis Leaf Senescence Reveals a Distinct Chronology of Processes and Regulation

This work presents a high-resolution time-course analysis of gene expression during development of a leaf from expansion through senescence. Enrichment in ontologies, sequence motifs, and transcription factor families within genes showing altered expression over time identified both metabolic pathways and potential regulators active at different stages of leaf development and senescence. Leaf senescence is an essential developmental process that impacts dramatically on crop yields and involves altered regulation of thousands of genes and many metabolic and signaling pathways, resulting in major changes in the leaf. The regulation of senescence is complex, and although senescence regulatory genes have been characterized, there is little information on how these function in the global control of the process. We used microarray analysis to obtain a high-resolution time-course profile of gene expression during development of a single leaf over a 3-week period to senescence. A complex experimental design approach and a combination of methods were used to extract high-quality replicated data and to identify differentially expressed genes. The multiple time points enable the use of highly informative clustering to reveal distinct time points at which signaling and metabolic pathways change. Analysis of motif enrichment, as well as comparison of transcription factor (TF) families showing altered expression over the time course, identify clear groups of TFs active at different stages of leaf development and senescence. These data enable connection of metabolic processes, signaling pathways, and specific TF activity, which will underpin the development of network models to elucidate the process of senescence.

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.

Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis

The authors generated a high-resolution time series of Arabidopsis thaliana gene expression following infection with the fungal pathogen Botrytis cinerea. Computational analysis of this large data set identified the timing of specific processes and regulatory events in the host plant and showed a role for the transcription factor TGA3 in the defense response against the fungal pathogen. Transcriptional reprogramming forms a major part of a plant’s response to pathogen infection. Many individual components and pathways operating during plant defense have been identified, but our knowledge of how these different components interact is still rudimentary. We generated a high-resolution time series of gene expression profiles from a single Arabidopsis thaliana leaf during infection by the necrotrophic fungal pathogen Botrytis cinerea. Approximately one-third of the Arabidopsis genome is differentially expressed during the first 48 h after infection, with the majority of changes in gene expression occurring before significant lesion development. We used computational tools to obtain a detailed chronology of the defense response against B. cinerea, highlighting the times at which signaling and metabolic processes change, and identify transcription factor families operating at different times after infection. Motif enrichment and network inference predicted regulatory interactions, and testing of one such prediction identified a role for TGA3 in defense against necrotrophic pathogens. These data provide an unprecedented level of detail about transcriptional changes during a defense response and are suited to systems biology analyses to generate predictive models of the gene regulatory networks mediating the Arabidopsis response to B. cinerea.

Cesium toxicity in Arabidopsis

Cesium (Cs) is chemically similar to potassium (K). However, although K is an essential element, Cs is toxic to plants. Two contrasting hypotheses to explain Cs toxicity have been proposed: (1) extracellular Cs+ prevents K+ uptake and, thereby, induces K starvation; and (2) intracellular Cs+ interacts with vital K+-binding sites in proteins, either competitively or noncompetitively, impairing their activities. We tested these hypotheses with Arabidopsis (Arabidopsis thaliana). Increasing the Cs concentration in the agar ([Cs]agar) on which Arabidopsis were grown reduced shoot growth. Increasing the K concentration in the agar ([K]agar) increased the [Cs]agar at which Cs toxicity was observed. However, although increasing [Cs]agar reduced shoot K concentration ([K]shoot), the decrease in shoot growth appeared unrelated to [K]shoot per se. Furthermore, the changes in gene expression in Cs-intoxicated plants differed from those of K-starved plants, suggesting that Cs intoxication was not perceived genetically solely as K starvation. In addition to reducing [K]shoot, increasing [Cs]agar also increased shoot Cs concentration ([Cs]shoot), but shoot growth appeared unrelated to [Cs]shoot per se. The relationship between shoot growth and [Cs]shoot/[K]shoot suggested that, at a nontoxic [Cs]shoot, growth was determined by [K]shoot but that the growth of Cs-intoxicated plants was related to the [Cs]shoot/[K]shoot quotient. This is consistent with Cs intoxication resulting from competition between K+ and Cs+ for K+-binding sites on essential proteins.

A Virulent Strain of Deformed Wing Virus (DWV) of Honeybees (Apis mellifera) Prevails after Varroa destructor-Mediated, or In Vitro, Transmission

The globally distributed ectoparasite Varroa destructor is a vector for viral pathogens of the Western honeybee (Apis mellifera), in particular the Iflavirus Deformed Wing Virus (DWV). In the absence of Varroa low levels DWV occur, generally causing asymptomatic infections. Conversely, Varroa-infested colonies show markedly elevated virus levels, increased overwintering colony losses, with impairment of pupal development and symptomatic workers. To determine whether changes in the virus population were due Varroa amplifying and introducing virulent virus strains and/or suppressing the host immune responses, we exposed Varroa-naïve larvae to oral and Varroa-transmitted DWV. We monitored virus levels and diversity in developing pupae and associated Varroa, the resulting RNAi response and transcriptome changes in the host. Exposed pupae were stratified by Varroa association (presence/absence) and virus levels (low/high) into three groups. Varroa-free pupae all exhibited low levels of a highly diverse DWV population, with those exposed per os (group NV) exhibiting changes in the population composition. Varroa-associated pupae exhibited either low levels of a diverse DWV population (group VL) or high levels of a near-clonal virulent variant of DWV (group VH). These groups and unexposed controls (C) could be also discriminated by principal component analysis of the transcriptome changes observed, which included several genes involved in development and the immune response. All Varroa tested contained a diverse replicating DWV population implying the virulent variant present in group VH, and predominating in RNA-seq analysis of temporally and geographically separate Varroa-infested colonies, was selected upon transmission from Varroa, a conclusion supported by direct injection of pupae in vitro with mixed virus populations. Identification of a virulent variant of DWV, the role of Varroa in its transmission and the resulting host transcriptome changes furthers our understanding of this important viral pathogen of honeybees.

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 F1 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 Ca2+ and Mg2+ uptake and homeostasis at the molecular level.

A local regulatory network around three NAC transcription factors in stress responses and senescence in Arabidopsis leaves

Summary A model is presented describing the gene regulatory network surrounding three similar NAC transcription factors that have roles in Arabidopsis leaf senescence and stress responses. ANAC019, ANAC055 and ANAC072 belong to the same clade of NAC domain genes and have overlapping expression patterns. A combination of promoter DNA/protein interactions identified using yeast 1-hybrid analysis and modelling using gene expression time course data has been applied to predict the regulatory network upstream of these genes. Similarities and divergence in regulation during a variety of stress responses are predicted by different combinations of upstream transcription factors binding and also by the modelling. Mutant analysis with potential upstream genes was used to test and confirm some of the predicted interactions. Gene expression analysis in mutants of ANAC019 and ANAC055 at different times during leaf senescence has revealed a distinctly different role for each of these genes. Yeast 1-hybrid analysis is shown to be a valuable tool that can distinguish clades of binding proteins and be used to test and quantify protein binding to predicted promoter motifs.

Contrasting arbuscular mycorrhizal communities colonizing different host plants show a similar response to a soil phosphorus concentration gradient

High soil phosphorus (P) concentration is frequently shown to reduce root colonization by arbuscular mycorrhizal (AM) fungi, but the influence of P on the diversity of colonizing AM fungi is uncertain. We used terminal restriction fragment length polymorphism (T-RFLP) of 18S rDNA and cloning to assess diversity of AM fungi colonizing maize (Zea mays), soybean (Glycene max) and field violet (Viola arvensis) at three time points in one season along a P gradient of 10–280 mg l−1 in the field. Percentage AM colonization changed between sampling time points but was not reduced by high soil P except in maize. There was no significant difference in AM diversity between sampling time points. Diversity was reduced at concentrations of P > 25 mg l−1, particularly in maize and soybean. Both cloning and T-RFLP indicated differences between AM communities in the different host species. Host species was more important than soil P in determining the AM community, except at the highest P concentration. Our results show that the impact of soil P on the diversity of AM fungi colonizing plants was broadly similar, despite the fact that different plants contained different communities. However, subtle differences in the response of the AM community in each host were evident.

A method to assess taxonomic variation in shoot caesium concentration among flowering plants

A method was developed to obtain relative shoot caesium (Cs) concentration data from the literature and assess the influence of plant taxonomy on these values. A residual maximum likelihood (REML) analysis was performed on data from 14 published studies, after these data were log-transformed to adjust for between-study differences in means and variances. There were two orders of magnitude difference between the lowest and highest relative shoot Cs concentration of the 136 taxa. Hierarchical nested analysis of variance revealed more than 40% of the variation in relative shoot Cs concentration was at the level of family or above. Dicotyledons (Magnoliopsida) had three-fold higher mean relative shoot Cs concentrations than monocotyledons (Liliopsida) whilst differences were also observed at lower taxonomic levels. The Caryophyllidae had the highest mean relative shoot Cs concentration among superorders; this group of plants contains many halophyte taxa and crop derivatives (e.g. beets, quinoas, buckwheats and amaranths). This method could inform soil-to-plant Cs transfer models and identify taxa with high Cs accumulation patterns that may have phytoremediation potential. The method reported could be used to study the accumulation of other elements in plants.

Temperature, light and nitrate sensing coordinate Arabidopsis seed dormancy cycling, resulting in winter and summer annual phenotypes

Seeds use environmental cues to sense the seasons and their surroundings to initiate the life cycle of the plant. The dormancy cycling underlying this process is extensively described, but the molecular mechanism is largely unknown. To address this we selected a range of representative genes from published array experiments in the laboratory, and investigated their expression patterns in seeds of Arabidopsis ecotypes with contrasting life cycles over an annual dormancy cycle in the field. We show how mechanisms identified in the laboratory are coordinated in response to the soil environment to determine the dormancy cycles that result in winter and summer annual phenotypes. Our results are consistent with a seed-specific response to seasonal temperature patterns (temporal sensing) involving the gene DELAY OF GERMINATION 1 (DOG1) that indicates the correct season, and concurrent temporally driven co-opted mechanisms that sense spatial signals, i.e. nitrate, via CBL-INTERACTING PROTEIN KINASE 23 (CIPK23) phosphorylation of the NITRATE TRANSPORTER 1 (NRT1.1), and light, via PHYTOCHROME A (PHYA). In both ecotypes studied, when all three genes have low expression there is enhanced GIBBERELLIN 3 BETA-HYDROXYLASE 1 (GA3ox1) expression, exhumed seeds have the potential to germinate in the laboratory, and the initiation of seedling emergence occurs following soil disturbance (exposure to light) in the field. Unlike DOG1, the expression of MOTHER of FLOWERING TIME (MFT) has an opposite thermal response in seeds of the two ecotypes, indicating a role in determining their different dormancy cycling phenotypes.