Wheat MYOSIN-RESEMBLING CHLOROPLAST PROTEIN controls B-type starch granule initiation timing during endosperm development.

Molecular factors that contribute to the diverse spatial and temporal patterns of starch granule initiation between species and organs are poorly understood. Wheat (Triticum sp.) endosperm contains both large A-type granules initiated during early grain development and small B-type granules that initiate about 10 to 15 days later. Here, we identify that the MYOSIN-RESEMBLING CHLOROPLAST PROTEIN (MRC) is required for the correct timing of B-type granule initiation in wheat endosperm during grain development. MRC is expressed in the endosperm exclusively in early grain development, before B-type granule initiation. We isolated three independent TILLING mutants of tetraploid wheat (Triticum turgidum cv. 'Kronos') with premature stop or missense mutations in the A-genome homeolog, which we showed to be the only active homeolog in tetraploid wheat due to a disruption of the B-genome homeolog. The mrc mutants had significantly smaller A-type granules and a higher relative volume of B-type granules in the endosperm than the wild type. Whereas B-type granules initiated 15 to 20 days post-anthesis (dpa) in the wild type, they appeared as early as 10 dpa in the mrc-1 mutant. These results suggest a temporal role for MRC in repressing B-type granule initiation, providing insight into how the distinct biochemical mechanisms that control A- and B-type granule initiation are regulated. This role of MRC in the wheat endosperm is distinct from the previously described role of Arabidopsis (Arabidopsis thaliana) MRC in promoting granule initiation in leaves, providing an example of functional diversification among granule initiation proteins.

Gene expression profile of the developing endosperm in durum wheat provides insight into starch biosynthesis.

BACKGROUND: Durum wheat (Triticum turgidum subsp. durum) is widely grown for pasta production, and more recently, is gaining additional interest due to its resilience to warm, dry climates and its use as an experimental model for wheat research. Like in bread wheat, the starch and protein accumulated in the endosperm during grain development are the primary contributors to the calorific value of durum grains. RESULTS: To enable further research into endosperm development and storage reserve synthesis, we generated a high-quality transcriptomics dataset from developing endosperms of variety Kronos, to complement the extensive mutant resources available for this variety. Endosperms were dissected from grains harvested at eight timepoints during grain development (6 to 30 days post anthesis (dpa)), then RNA sequencing was used to profile the transcriptome at each stage. The largest changes in gene expression profile were observed between the earlier timepoints, prior to 15 dpa. We detected a total of 29,925 genes that were significantly differentially expressed between at least two timepoints, and clustering analysis revealed nine distinct expression patterns. We demonstrate the potential of our dataset to provide new insights into key processes that occur during endosperm development, using starch metabolism as an example. CONCLUSION: We provide a valuable resource for studying endosperm development in this increasingly important crop species.

Spatial expression patterns of genes encoding sugar sensors in leaves of C4 and C3 grasses.

BACKGROUND AND AIMS: The mechanisms of sugar sensing in grasses remain elusive, especially those using C4 photosynthesis even though a large proportion of the world's agricultural crops utilize this pathway. We addressed this gap by comparing the expression of genes encoding components of sugar sensors in C3 and C4 grasses, with a focus on source tissues of C4 grasses. Given C4 plants evolved into a two-cell carbon fixation system, it was hypothesized this may have also changed how sugars were sensed. METHODS: For six C3 and eight C4 grasses, putative sugar sensor genes were identified for target of rapamycin (TOR), SNF1-related kinase 1 (SnRK1), hexokinase (HXK) and those involved in the metabolism of the sugar sensing metabolite trehalose-6-phosphate (T6P) using publicly available RNA deep sequencing data. For several of these grasses, expression was compared in three ways: source (leaf) versus sink (seed), along the gradient of the leaf, and bundle sheath versus mesophyll cells. KEY RESULTS: No positive selection of codons associated with the evolution of C4 photosynthesis was identified in sugar sensor proteins here. Expressions of genes encoding sugar sensors were relatively ubiquitous between source and sink tissues as well as along the leaf gradient of both C4 and C3 grasses. Across C4 grasses, SnRK1ß1 and TPS1 were preferentially expressed in the mesophyll and bundle sheath cells, respectively. Species-specific differences of gene expression between the two cell types were also apparent. CONCLUSIONS: This comprehensive transcriptomic study provides an initial foundation for elucidating sugar-sensing genes within major C4 and C3 crops. This study provides some evidence that C4 and C3 grasses do not differ in how sugars are sensed. While sugar sensor gene expression has a degree of stability along the leaf, there are some contrasts between the mesophyll and bundle sheath cells.

The roles of divergent and parallel molecular evolution contributing to thermal adaptive strategies in trees.

Local adaptation is a driver of biological diversity, and species may develop analogous (parallel evolution) or alternative (divergent evolution) solutions to similar ecological challenges. We expect these adaptive solutions would culminate in both phenotypic and genotypic signals. Using two Eucalyptus species (Eucalyptus grandis and Eucalyptus tereticornis) with overlapping distributions grown under contrasting 'local' temperature conditions to investigate the independent contribution of adaptation and plasticity at molecular, physiological and morphological levels. The link between gene expression and traits markedly differed between species. Divergent evolution was the dominant pattern driving adaptation (91% of all significant genes); but overlapping gene (homologous) responses were dependent on the determining factor (plastic, adaptive or genotype by environment interaction). Ninety-eight percent of the plastic homologs were similarly regulated, while 50% of the adaptive homologs and 100% of the interaction homologs were antagonistical. Parallel evolution for the adaptive effect in homologous genes was greater than expected but not in favour of divergent evolution. Heat shock proteins for E. grandis were almost entirely driven by adaptation, and plasticity in E. tereticornis. These results suggest divergent molecular evolutionary solutions dominated the adaptive mechanisms among species, even in similar ecological circumstances. Suggesting that tree species with overlapping distributions are unlikely to equally persist in the future.

High intrinsic water use efficiency is underpinned by high stomatal aperture and guard cell potassium flux in C3 and C4 grasses grown at glacial CO2 and low light.

We compared how stomatal morphology and physiology control intrinsic leaf water use efficiency (iWUE) in two C3 and six C4 grasses grown at ambient (400 µmol mol-1) or glacial CO2 (180 µmol mol-1) and high (1000 µmol m-2 s-1) or low light intensity (200 µmol m-2 s-1). C4 grasses tended to have higher iWUE and CO2 assimilation rates, and lower stomatal conductance (gs), operational stomatal aperture (aop), and guard cell K+ influx rate relative to C3 grasses, while stomatal size (SS) and stomatal density (SD) did not vary according to the photosynthetic type. Overall, iWUE and gs depended most on aop and density of open stomata. In turn, aop correlated with K+ influx, stomatal opening speed on transition to high light, and SS. Species with higher SD had smaller and faster-opening stomata. Although C4 grasses operated with lower gs and aop at ambient CO2, they showed a greater potential to open stomata relative to maximal stomatal conductance (gmax), indicating heightened stomatal sensitivity and control. We uncovered promising links between aop, gs, iWUE, and K+ influx among C4 grasses, and differential K+ influx responses of C4 guard cells to low light, revealing molecular targets for improving iWUE in C4 crops.

Loss of PROTEIN TARGETING TO STARCH 2 has variable effects on starch synthesis across organs and species.

Recent work has identified several proteins involved in starch granule initiation, the first step of starch synthesis. However, the degree of conservation in the granule initiation process remains poorly understood, especially among grass species differing in patterns of carbohydrate turnover in leaves, and granule morphology in the endosperm. We therefore compared mutant phenotypes of Hordeum vulgare (barley), Triticum turgidum (durum wheat), and Brachypodium distachyon defective in PROTEIN TARGETING TO STARCH 2 (PTST2), a key granule initiation protein. We report striking differences across species and organs. Loss of PTST2 from leaves resulted in fewer, larger starch granules per chloroplast and normal starch content in wheat, fewer granules per chloroplast and lower starch content in barley, and almost complete loss of starch in Brachypodium. The loss of starch in Brachypodium leaves was accompanied by high levels of ADP-glucose and detrimental effects on growth and physiology. Additionally, we found that loss of PTST2 increased granule initiation in Brachypodium amyloplasts, resulting in abnormal compound granule formation throughout the seed. These findings suggest that the importance of PTST2 varies greatly with the genetic and developmental background and inform the extent to which the gene can be targeted to improve starch in crops.

The transcriptomic responses of C4 grasses to subambient CO2 and low light are largely species specific and only refined by photosynthetic subtype.

Three subtypes of C4 photosynthesis exist (NADP-ME, NAD-ME and PEPCK), each known to be beneficial under specific environmental conditions. However, the influence of photosynthetic subtype on transcriptomic plasticity, as well as the genes underpinning this variability, remain largely unknown. Here, we comprehensively investigate the responses of six C4 grass species, spanning all three C4 subtypes, to two controlled environmental stresses: low light (200 µmol m-2  sec-1 ) and glacial CO2 (subambient; 180 ppm). We identify a susceptibility within NADP-ME species to glacial CO2 . Notably, although glacial CO2 phenotypes could be tied to C4 subtype, biochemical and transcriptomic responses to glacial CO2 were largely species specific. Nevertheless, we were able to identify subtype specific subsets of significantly differentially expressed transcripts which link resource acquisition and allocation to NADP-ME species susceptibility to glacial CO2 . Here, low light phenotypes were comparable across species with no clear subtype response, while again, transcriptomic responses to low light were largely species specific. However, numerous functional similarities were noted within the transcriptomic responses to low light, suggesting these responses are functionally relatively conserved. Additionally, PEPCK species exhibited heightened regulation of transcripts related to metabolism in response to both stresses, likely tied to their C4 metabolic pathway. These results highlight the influence that both species and subtype can have on plant responses to abiotic stress, building on our mechanistic understanding of acclimation within C4 grasses and highlighting avenues for future crop improvements.

STARCH SYNTHASE 4 is required for normal starch granule initiation in amyloplasts of wheat endosperm.

Starch granule initiation is poorly understood at the molecular level. The glucosyltransferase, STARCH SYNTHASE 4 (SS4), plays a central role in granule initiation in Arabidopsis leaves, but its function in cereal endosperms is unknown. We investigated the role of SS4 in wheat, which has a distinct spatiotemporal pattern of granule initiation during grain development. We generated TILLING mutants in tetraploid wheat (Triticum turgidum) that are defective in both SS4 homoeologs. The morphology of endosperm starch was examined in developing and mature grains. SS4 deficiency led to severe alterations in endosperm starch granule morphology. During early grain development, while the wild-type initiated single 'A-type' granules per amyloplast, most amyloplasts in the mutant formed compound granules due to multiple initiations. This phenotype was similar to mutants deficient in B-GRANULE CONTENT 1 (BGC1). SS4 deficiency also reduced starch content in leaves and pollen grains. We propose that SS4 and BGC1 are required for the proper control of granule initiation during early grain development that leads to a single A-type granule per amyloplast. The absence of either protein results in a variable number of initiations per amyloplast and compound granule formation.

Genome-wide identification and characterisation of Aquaporins in Nicotiana tabacum and their relationships with other Solanaceae species.

BACKGROUND: Cellular membranes are dynamic structures, continuously adjusting their composition, allowing plants to respond to developmental signals, stresses, and changing environments. To facilitate transmembrane transport of substrates, plant membranes are embedded with both active and passive transporters. Aquaporins (AQPs) constitute a major family of membrane spanning channel proteins that selectively facilitate the passive bidirectional passage of substrates across biological membranes at an astonishing 108 molecules per second. AQPs are the most diversified in the plant kingdom, comprising of five major subfamilies that differ in temporal and spatial gene expression, subcellular protein localisation, substrate specificity, and post-translational regulatory mechanisms; collectively providing a dynamic transportation network spanning the entire plant. Plant AQPs can transport a range of solutes essential for numerous plant processes including, water relations, growth and development, stress responses, root nutrient uptake, and photosynthesis. The ability to manipulate AQPs towards improving plant productivity, is reliant on expanding our insight into the diversity and functional roles of AQPs. RESULTS: We characterised the AQP family from Nicotiana tabacum (NtAQPs; tobacco), a popular model system capable of scaling from the laboratory to the field. Tobacco is closely related to major economic crops (e.g. tomato, potato, eggplant and peppers) and itself has new commercial applications. Tobacco harbours 76 AQPs making it the second largest characterised AQP family. These fall into five distinct subfamilies, for which we characterised phylogenetic relationships, gene structures, protein sequences, selectivity filter compositions, sub-cellular localisation, and tissue-specific expression. We also identified the AQPs from tobacco's parental genomes (N. sylvestris and N. tomentosiformis), allowing us to characterise the evolutionary history of the NtAQP family. Assigning orthology to tomato and potato AQPs allowed for cross-species comparisons of conservation in protein structures, gene expression, and potential physiological roles. CONCLUSIONS: This study provides a comprehensive characterisation of the tobacco AQP family, and strengthens the current knowledge of AQP biology. The refined gene/protein models, tissue-specific expression analysis, and cross-species comparisons, provide valuable insight into the evolutionary history and likely physiological roles of NtAQPs and their Solanaceae orthologs. Collectively, these results will support future functional studies and help transfer basic research to applied agriculture.

Sugar sensing responses to low and high light in leaves of the C4 model grass Setaria viridis.

Although sugar regulates photosynthesis, the signalling pathways underlying this process remain elusive, especially for C4 crops. To address this knowledge gap and identify potential candidate genes, we treated Setaria viridis (C4 model) plants acclimated to medium light intensity (ML, 500 µmol m-2 s-1) with low (LL, 50 µmol m-2 s-1) or high (HL, 1000 µmol m-2 s-1) light for 4 d and observed the consequences on carbon metabolism and the transcriptome of source leaves. LL impaired photosynthesis and reduced leaf content of signalling sugars (glucose, sucrose, and trehalose-6-phosphate). In contrast, HL strongly induced sugar accumulation without repressing photosynthesis. LL more profoundly impacted the leaf transcriptome, including photosynthetic genes. LL and HL contrastingly altered the expression of hexokinase (HXK) and sucrose-non-fermenting 1 (Snf1)-related protein kinase 1 (SnRK1) sugar sensors and trehalose pathway genes. The expression of key target genes of HXK and SnRK1 were affected by LL and sugar depletion, while surprisingly HL and strong sugar accumulation only slightly repressed the SnRK1 signalling pathway. In conclusion, we demonstrate that LL profoundly impacted photosynthesis and the transcriptome of S. viridis source leaves, while HL altered sugar levels more than LL. We also present the first evidence that sugar signalling pathways in C4 source leaves may respond to light intensity and sugar accumulation differently from C3 source leaves.

The methylome is altered for plants in a high CO2 world: Insights into the response of a wild plant population to multigenerational exposure to elevated atmospheric [CO2 ].

Unravelling plant responses to rising atmospheric CO2 concentration ([CO2 ]) has largely focussed on plastic functional attributes to single generation [CO2 ] exposure. Quantifying the consequences of long-term, decadal multigenerational exposure to elevated [CO2 ] and the genetic changes that may underpin evolutionary mechanisms with [CO2 ] as a driver remain largely unexplored. Here, we investigated both plastic and evolutionary plant responses to elevated [CO2 ] by applying multi-omic technologies using populations of Plantago lanceolata L., grown in naturally high [CO2 ] for many generations in a CO2 spring. Seed from populations at the CO2 spring and an adjacent control site (ambient [CO2 ]) were grown in a common environment for one generation, and then offspring were grown in ambient or elevated [CO2 ] growth chambers. Low overall genetic differentiation between the CO2 spring and control site populations was found, with evidence of weak selection in exons. We identified evolutionary divergence in the DNA methylation profiles of populations derived from the spring relative to the control population, providing the first evidence that plant methylomes may respond to elevated [CO2 ] over multiple generations. In contrast, growth at elevated [CO2 ] for a single generation induced limited methylome remodelling (an order of magnitude fewer differential methylation events than observed between populations), although some of this appeared to be stably transgenerationally inherited. In all, 59 regions of the genome were identified where transcripts exhibiting differential expression (associated with single generation or long-term natural exposure to elevated [CO2 ]) co-located with sites of differential methylation or with single nucleotide polymorphisms exhibiting significant inter-population divergence. This included genes in pathways known to respond to elevated [CO2 ], such as nitrogen use efficiency and stomatal patterning. This study provides the first indication that DNA methylation may contribute to plant adaptation to future atmospheric [CO2 ] and identifies several areas of the genome that are targets for future study.

Comparative analysis of thylakoid protein complexes in the mesophyll and bundle sheath cells from C3 , C4 and C3 -C4 Paniceae grasses.

To better understand the coordination between dark and light reactions during the transition from C3 to C4 photosynthesis, we optimized a method for separating thylakoids from mesophyll (MC) and bundle sheath cells (BSCs) across different plant species. We grew six Paniceae grasses including representatives from the C3 , C3 -C4 and C4 photosynthetic types and all three C4 biochemical subtypes [nicotinamide adenine dinucleotide phosphate-dependent malic enzyme (NADP-ME), nicotinamide adenine dinucleotide-dependent malic enzyme (NAD-ME) and phosphoenolpyruvate carboxykinase (PEPCK)] in addition to Zea mays under control conditions (1000 µmol quanta m-2  s-1 and 400 ppm of CO2 ). Proteomics analysis of thylakoids under native conditions, using blue native polyacrylamide gel electrophoresis followed by liquid chromatography-mass spectrometry (LC-MS), demonstrated the presence of subunits of all light-reaction-related complexes in all species and cell types. C4 NADP-ME species showed a higher photosystems I/II ratio and a clear accumulation of the NADH dehydrogenase-like complexes in BSCs, while Cytb6 f was more abundant in BSCs of C4 NAD-ME species. The C4 PEPCK species showed no clear differences between cell types. Our study presents, for the first time, a good separation between BSC and MC for a C3 -C4 intermediate grass which did not show noticeable differences in the distribution of the thylakoid complexes. For the NADP-ME species Panicum antidotale, growth at glacial CO2 (180 ppm of CO2 ) had no effect on the distribution of the light-reaction complexes, while growth at low light (200 µmol quanta m-2  s-1 ) promoted the accumulation of light-harvesting proteins in both cell types. These results add to our understanding of thylakoid distribution across photosynthetic types and subtypes, and introduce thylakoid distribution between the MC and BSC of a C3 -C4 intermediate species.

Investigating the NAD-ME biochemical pathway within C4 grasses using transcript and amino acid variation in C4 photosynthetic genes.

Expanding knowledge of the C4 photosynthetic pathway can provide key information to aid biological improvements to crop photosynthesis and yield. While the C4 NADP-ME pathway is well characterised, there is increasing agricultural and bioengineering interest in the comparably understudied NAD-ME and PEPCK pathways. Within this study, a systematic identification of key differences across species has allowed us to investigate the evolution of C4-recruited genes in one C3 and eleven C4 grasses (Poaceae) spanning two independent origins of C4 photosynthesis. We present evidence for C4-specific paralogs of NAD-malic enzyme 2, MPC1 and MPC2 (mitochondrial pyruvate carriers) via increased transcript abundance and associated rates of evolution, implicating them as genes recruited to perform C4 photosynthesis within NAD-ME and PEPCK subtypes. We then investigate the localisation of AspAT across subtypes, using novel and published evidence to place AspAT3 in both the cytosol and peroxisome. Finally, these findings are integrated with transcript abundance of previously identified C4 genes to provide an updated model for C4 grass NAD-ME and PEPCK photosynthesis. This updated model allows us to develop on the current understanding of NAD-ME and PEPCK photosynthesis in grasses, bolstering our efforts to understand the evolutionary 'path to C4' and improve C4 photosynthesis.

Plant adaptation or acclimation to rising CO2 ? Insight from first multigenerational RNA-Seq transcriptome.

Atmospheric carbon dioxide (CO2 ) directly determines the rate of plant photosynthesis and indirectly effects plant productivity and fitness and may therefore act as a selective pressure driving evolution, but evidence to support this contention is sparse. Using Plantago lanceolata L. seed collected from a naturally high CO2 spring and adjacent ambient CO2 control site, we investigated multigenerational response to future, elevated atmospheric CO2 . Plants were grown in either ambient or elevated CO2 (700 µmol mol-1 ), enabling for the first time, characterization of the functional and population genomics of plant acclimation and adaptation to elevated CO2 . This revealed that spring and control plants differed significantly in phenotypic plasticity for traits underpinning fitness including above-ground biomass, leaf size, epidermal cell size and number and stomatal density and index. Gene expression responses to elevated CO2 (acclimation) were modest [33-131 genes differentially expressed (DE)], whilst those between control and spring plants (adaptation) were considerably larger (689-853 DE genes). In contrast, population genomic analysis showed that genetic differentiation between spring and control plants was close to zero, with no fixed differences, suggesting that plants are adapted to their native CO2 environment at the level of gene expression. An unusual phenotype of increased stomatal index in spring but not control plants in elevated CO2 correlated with altered expression of stomatal patterning genes between spring and control plants for three loci (YODA, CDKB1;1 and SCRM2) and between ambient and elevated CO2 for four loci (ER, YODA, MYB88 and BCA1). We propose that the two positive regulators of stomatal number (SCRM2) and CDKB1;1 when upregulated act as key controllers of stomatal adaptation to elevated CO2 . Combined with significant transcriptome reprogramming of photosynthetic and dark respiration and enhanced growth in spring plants, we have identified the potential basis of plant adaptation to high CO2 likely to occur over coming decades.

Effect of reduced irradiance on 13C uptake, gene expression and protein activity of the seagrass Zostera muelleri.

Photosynthesis in the seagrass Zostera muelleri remains poorly understood. We investigated the effect of reduced irradiance on the incorporation of 13C, gene expression of photosynthetic, photorespiratory and intermediates recycling genes as well as the enzymatic content and activity of Rubisco and PEPC within Z. muelleri. Following 48 h of reduced irradiance, we found that i) there was a ∼7 fold reduction in 13C incorporation in above ground tissue, ii) a significant down regulation of photosynthetic, photorespiratory and intermediates recycling genes and iii) no significant difference in enzyme activity and content. We propose that Z. muelleri is able to alter its physiology in order to reduce the amount of C lost through photorespiration to compensate for the reduced carbon assimilation as a result of reduced irradiance. In addition, the first estimated rate constant (Kcat) and maximum rates of carboxylation (Vcmax) of Rubisco is reported for the first time for Z. muelleri.