Transcriptomics of C4 photosynthesis in rice paddy

This is the author accepted manuscript. The final version is available from the American Society of Plant Biologists via http://dx.doi.org/​10.​1104/​pp.​15.​00889The C₄ pathway is a highly complex trait that increases photosynthetic efficiency in over sixty plant lineages. Although the majority of C₄ plants occupy disturbed, arid and nutrient-poor habitats, some grow in high-nutrient, waterlogged conditions. One such example is Echinochloa glabrescens, which is an aggressive weed of rice paddies. We generated comprehensive transcriptome datasets for C₄ E. glabrescens and C₃ rice to identify genes associated with adaption to waterlogged, nutrient-replete conditions, but also used the data to better understand how C₄ photosynthesis operates in these conditions. Leaves of E. glabrescens exhibited classical Kranz anatomy with lightly lobed mesophyll cells having low chloroplast coverage. As with rice and other hygrophytic C₃ species, leaves of E. glabrescens accumulated a chloroplastic phosphoenolpyruvate carboxylase protein, albeit at reduced amounts relative to rice. The arid-grown species Setaria italica (C₄) and Brachypodium distachyon (C₃) were also found to accumulate chloroplastic PEPC. We identified a molecular signature associated with C₄ photosynthesis in nutrient-replete, waterlogged conditions that is highly similar to those previously reported from C₄ plants that grow in more arid conditions. We also identified a cohort of genes that have been subjected to a selective sweep associated with growth in paddy conditions. Overall, this approach highlights the value of using wild species such as weeds to identify adaptions to specific conditions associated with high-yielding crops in agriculture

Targeted knockdown of ribulose-1, 5-bisphosphate carboxylase-oxygenase in rice mesophyll cells.

We generated antisense constructs targeting two of the five Rubisco small subunit genes (OsRBCS2 and 4) which account for between 30-40 % of the RBCS transcript abundance in leaf blades. The constructs were driven by a maize phosphoenolpyruvate carboxylase (PEPC) promoter known to have enriched expression in mesophyll cells (MCs). In the resulting lines leaf, Rubisco protein content was reduced by between 30-50 % and CO2 assimilation rate was limited under photorespiratory and non-photorespiratory conditions. A relationship between Rubisco protein content and CO2 assimilation rate was found. This was associated with a significant reduction in dry biomass accumulation and grain yield of between 37-70%. In addition to serving as a resource for reducing Rubisco accumulation in a cell-preferential manner, these lines allow us to characterize gene function and isoform specific suppression on photosynthesis and growth. Our results suggest that the knockdown of multiple genes is required to completely reduce Rubisco accumulation in MCs

A Partial C4 Photosynthetic Biochemical Pathway in Rice.

Introduction of a C4 photosynthetic pathway into C3 rice (Oryza sativa) requires installation of a biochemical pump that concentrates CO2 at the site of carboxylation in modified bundle sheath cells. To investigate the feasibility of this, we generated a quadruple line that simultaneously accumulates four of the core C4 photosynthetic enzymes from the NADP-malic enzyme subtype, phosphoenolpyruvate carboxylase (ZmPEPC), NADP-malate dehydrogenase (ZmNADP-MDH), NADP-malic enzyme (ZmNADP-ME), and pyruvate phosphate dikinase (ZmPPDK). This led to enhanced enzyme activity and mild phenotypic perturbations but was largely neutral in its effects on photosynthetic rate. Measurements of the flux of 13CO2 through photosynthetic metabolism revealed a significant increase in the incorporation of 13C into malate, consistent with increased fixation of 13CO2 via PEP carboxylase in lines expressing the maize PEPC enzyme. However, there was no significant differences in labeling of 3-phosphoglycerate (3PGA) indicating that there was no carbon flux through NADP-ME into the Calvin-Benson cycle. There was also no significant difference in labeling of phosphoenolpyruvate (PEP) indicating that there was no carbon flux through PPDK. Crossing the quadruple line with a line with reduced glycine decarboxylase H-protein (OsGDCH) abundance led to a photosynthetic phenotype characteristic of the reduced OsGDCH line and higher labeling of malate, aspartate and citrate than in the quintuple line. There was evidence of 13C labeling of aspartate indicating 13CO2 fixation into oxaloacetate by PEPC and conversion to aspartate by the endogenous aspartate aminotransferase activity. While Kranz anatomy or other anatomical modifications have not yet been installed in these plants to enable a fully functional C4 cycle, these results demonstrate for the first-time a partial flux through the carboxylation phase of NADP-ME C4 metabolism in transgenic rice containing two of the key metabolic steps in the C4 pathway

Knockdown of glycine decarboxylase complex alters photorespiratory carbon isotope fractionation in Oryza sativa leaves

The influence of reduced glycine decarboxylase complex (GDC) activity on leaf atmosphere CO2 and 13CO2 exchange was tested in transgenic Oryza sativa with the GDC H-subunit knocked down in leaf mesophyll cells. Leaf measurements on transgenic gdch knockdown and wild-type plants were carried out in the light under photorespiratory and low photorespiratory conditions (i.e. 18.4 kPa and 1.84 kPa atmospheric O2 partial pressure, respectively), and in the dark. Under approximately current ambient O2 partial pressure (18.4 kPa pO2), the gdch knockdown plants showed an expected photorespiratory-deficient phenotype, with lower leaf net CO2 assimilation rates (A) than the wild-type. Additionally, under these conditions, the gdch knockdown plants had greater leaf net discrimination against 13CO2 (Δo) than the wild-type. This difference in Δo was in part due to lower 13C photorespiratory fractionation (f) ascribed to alternative decarboxylation of photorespiratory intermediates. Furthermore, the leaf dark respiration rate (Rd) was enhanced and the 13CO2 composition of respired CO2 (δ13CRd) showed a tendency to be more depleted in the gdch knockdown plants. These changes in Rd and δ13CRd were due to the amount and carbon isotopic composition of substrates available for dark respiration. These results demonstrate that impairment of the photorespiratory pathway affects leaf 13CO2 exchange, particularly the 13C decarboxylation fractionation associated with photorespiration.Research was funded by a C4 Rice Project grant from The Bill and Melinda Gates Foundation to IRRI (2012–2015) and to the University of Oxford (2015–2019); by the National Science Foundation, grant MCB-1146928; by the National Science Foundation, grant MRI0923562; and by the Russian Science Foundation, grant 16-16-00089

Evolutionary convergence of cell-specific gene expression in independent lineages of C-4 grasses

Leaves of almost all C(4) lineages separate the reactions of photosynthesis into the mesophyll (M) and bundle sheath (BS). The extent to which messenger RNA profiles of M and BS cells from independent C(4) lineages resemble each other is not known. To address this, we conducted deep sequencing of RNA isolated from the M and BS of Setaria viridis and compared these data with publicly available information from maize (Zea mays). This revealed a high correlation (r = 0.89) between the relative abundance of transcripts encoding proteins of the core C(4) pathway in M and BS cells in these species, indicating significant convergence in transcript accumulation in these evolutionarily independent C(4) lineages. We also found that the vast majority of genes encoding proteins of the C(4) cycle in S. viridis are syntenic to homologs used by maize. In both lineages, 122 and 212 homologous transcription factors were preferentially expressed in the M and BS, respectively. Sixteen shared regulators of chloroplast biogenesis were identified, 14 of which were syntenic homologs in maize and S. viridis. In sorghum (Sorghum bicolor), a third C(4) grass, we found that 82% of these trans-factors were also differentially expressed in either M or BS cells. Taken together, these data provide, to our knowledge, the first quantification of convergence in transcript abundance in the M and BS cells from independent lineages of C(4) grasses. Furthermore, the repeated recruitment of syntenic homologs from large gene families strongly implies that parallel evolution of both structural genes and trans-factors underpins the polyphyletic evolution of this highly complex trait in the monocotyledons

RNA-Seq based phylogeny recapitulates previous phylogeny of the genus Flaveria (Asteraceae) with some modifications

Abstract Background The genus Flaveria has been extensively used as a model to study the evolution of C4 photosynthesis as it contains C3 and C4 species as well as a number of species that exhibit intermediate types of photosynthesis. The current phylogenetic tree of the genus Flaveria contains 21 of the 23 known Flaveria species and has been previously constructed using a combination of morphological data and three non-coding DNA sequences (nuclear encoded ETS, ITS and chloroplast encoded trnL-F). Results Here we developed a new strategy to update the phylogenetic tree of 16 Flaveria species based on RNA-Seq data. The updated phylogeny is largely congruent with the previously published tree but with some modifications. We propose that the data collection method provided in this study can be used as a generic method for phylogenetic tree reconstruction if the target species has no genomic information. We also showed that a “F. pringlei” genotype recently used in a number of labs may be a hybrid between F. pringlei (C3) and F. angustifolia (C3-C4). Conclusions We propose that the new strategy of obtaining phylogenetic sequences outlined in this study can be used to construct robust trees in a larger number of taxa. The updated Flaveria phylogenetic tree also supports a hypothesis of stepwise and parallel evolution of C4 photosynthesis in the Flavaria clade

RNA-Seq based phylogeny recapitulates previous phylogeny of the genus Flaveria (Asteraceae) with some modifications.

BACKGROUND: The genus Flaveria has been extensively used as a model to study the evolution of C4 photosynthesis as it contains C3 and C4 species as well as a number of species that exhibit intermediate types of photosynthesis. The current phylogenetic tree of the genus Flaveria contains 21 of the 23 known Flaveria species and has been previously constructed using a combination of morphological data and three non-coding DNA sequences (nuclear encoded ETS, ITS and chloroplast encoded trnL-F). RESULTS: Here we developed a new strategy to update the phylogenetic tree of 16 Flaveria species based on RNA-Seq data. The updated phylogeny is largely congruent with the previously published tree but with some modifications. We propose that the data collection method provided in this study can be used as a generic method for phylogenetic tree reconstruction if the target species has no genomic information. We also showed that a "F. pringlei" genotype recently used in a number of labs may be a hybrid between F. pringlei (C3) and F. angustifolia (C3-C4). CONCLUSIONS: We propose that the new strategy of obtaining phylogenetic sequences outlined in this study can be used to construct robust trees in a larger number of taxa. The updated Flaveria phylogenetic tree also supports a hypothesis of stepwise and parallel evolution of C4 photosynthesis in the Flavaria clade

Transgenic maize phosphoenolpyruvate carboxylase alters leaf-atmosphere CO2 and 13CO2 exchanges in Oryza sativa.

The engineering process of C4 photosynthesis into C3 plants requires an increased activity of phosphoenolpyruvate carboxylase (PEPC) in the cytosol of leaf mesophyll cells. The literature varies on the physiological effect of transgenic maize (Zea mays) PEPC (ZmPEPC) leaf expression in Oryza sativa (rice). Therefore, to address this issue, leaf-atmosphere CO2 and 13CO2 exchanges were measured, both in the light (at atmospheric O2 partial pressure of 1.84 kPa and at different CO2 levels) and in the dark, in transgenic rice expressing ZmPEPC and wild-type (WT) plants. The in vitro PEPC activity was 25 times higher in the PEPC overexpressing (PEPC-OE) plants (~20% of maize) compared to the negligible activity in WT. In the PEPC-OE plants, the estimated fraction of carboxylation by PEPC (β) was ~6% and leaf net biochemical discrimination against 13CO2[Formula: see text] was ~ 2‰ lower than in WT. However, there were no differences in leaf net CO2 assimilation rates (A) between genotypes, while the leaf dark respiration rates (Rd) over three hours after light-dark transition were enhanced (~ 30%) and with a higher 13C composition [Formula: see text] in the PEPC-OE plants compared to WT. These data indicate that ZmPEPC in the PEPC-OE rice plants contributes to leaf carbon metabolism in both the light and in the dark. However, there are some factors, potentially posttranslational regulation and PEP availability, which reduce ZmPEPC activity in vivo

Lineage-specific gene radiations underlie the evolution of novel betalain pigmentation in Caryophyllales.

Betalain pigments are unique to the Caryophyllales and structurally and biosynthetically distinct from anthocyanins. Two key enzymes within the betalain synthesis pathway have been identified: 4,5-dioxygenase (DODA) that catalyzes the formation of betalamic acid and CYP76AD1, a cytochrome P450 gene that catalyzes the formation of cyclo-DOPA. We performed phylogenetic analyses to reveal the evolutionary history of the DODA and CYP76AD1 lineages and in the context of an ancestral reconstruction of pigment states we explored the evolution of these genes in relation to the complex evolution of pigments in Caryophylalles. Duplications within the CYP76AD1 and DODA lineages arose just before the origin of betalain pigmentation in the core Caryophyllales. The duplications gave rise to DODA-α and CYP76AD1-α isoforms that appear specific to betalain synthesis. Both betalain-specific isoforms were then lost or downregulated in the anthocyanic Molluginaceae and Caryophyllaceae. Our findings suggest a single origin of the betalain synthesis pathway, with neofunctionalization following gene duplications in the CYP76AD1 and DODA lineages. Loss of DODA-α and CYP76AD1-α in anthocyanic taxa suggests that betalain pigmentation has been lost twice in Caryophyllales, and exclusion of betalain pigments from anthocyanic taxa is mediated through gene loss or downregulation. [Correction added after online publication 13 May 2015: in the last two paragraphs of the Summary the gene name CYP761A was changed to CYP76AD1.].S.C. was supported by a grant to IRRI from the Bill and Melinda Gates Foundation and UKAID. This work was supported by a National Science Foundation award (grant numbers DEB 1354048 and DEB 1352907) to S.F.B., M.J.M. and S.A.S., and a NERC Independent Research Fellowship to S.F.B. The 1000 Plants (1KP) initiative, led by G.K.S.W., is funded by the Alberta Ministry of Enterprise and Advanced Education, Alberta Innovates Technology Futures (AITF), Innovates Centre of Research Excellence (iCORE), Musea Ventures and BGI-Shenzhen.This is the final version of the article. It first appeared from Wiley via http://dx.doi.org/10.1111/nph.1344