Model‐assisted identification of metabolic engineering strategies for Jatropha curcas lipid pathways

Efficient approaches to increase plant lipid production are necessary to meet current industrial demands for this important resource. While Jatropha curcas cell culture can be used for in vitro lipid production, scaling up the system for industrial applications requires an understanding of how growth conditions affect lipid metabolism and yield. Here we present a bottom‐up metabolic reconstruction of J. curcas supported with labeling experiments and biomass characterization under three growth conditions. We show that the metabolic model can accurately predict growth and distribution of fluxes in cell cultures and use these findings to pinpoint energy expenditures that affect lipid biosynthesis and metabolism. In addition, by using constraint‐based modeling approaches we identify network reactions whose joint manipulation optimizes lipid production. The proposed model and computational analyses provide a stepping stone for future rational optimization of other agronomically relevant traits in J. curcas.     

Expansion of MIR482/2118 by a class‐II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide‐binding site leucine‐rich repeat (NBS‐LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target‐gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class‐II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The ‘GosTE’ involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3‐bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co‐evolution of miR482/2118 and its NBS‐LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.     

Leucine‐rich repeat receptor‐like kinase II phylogenetics reveals five main clades throughout the plant kingdom

Receptor‐like kinases (RLKs) represent the largest group of cell surface receptors in plants. The monophyletic leucine‐rich repeat (LRR)‐RLK subfamily II is considered to contain the somatic embryogenesis receptor kinases (SERKs) and NSP‐interacting kinases known to be involved in developmental processes and cellular immunity in plants. There are only a few published studies on the phylogenetics of LRR‐RLKII; unfortunately these suffer from poor taxon/gene sampling. Hence, it is not clear how many and what main clades this family contains, let alone what structure–function relationships exist. We used 1342 protein sequences annotated as ‘SERK’ and ‘SERK‐like’ plus related sequences in order to estimate phylogeny within the LRR‐RLKII clade, using the nematode protein kinase Pelle as an outgroup. We reconstruct five main clades (LRR‐RLKII 1–5), in each of which the main pattern of land plant relationships re‐occurs, confirming previous hypotheses that duplication events happened in this gene subfamily prior to divergence among land plant lineages. We show that domain structures and intron–exon boundaries within the five clades are well conserved in evolution. Furthermore, phylogenetic patterns based on the separate LRR and kinase parts of LRR‐RLKs are incongruent: whereas the LRR part supports a LRR‐RLKII 2/3 sister group relationship, the kinase part supports clades 1/2. We infer that the kinase part includes few ‘radical’ amino acid changes compared with the LRR part. Finally, our results confirm that amino acids involved in each LRR‐RLKII–receptor complex interaction are located at N‐capping residues, and that the short amino acid motifs of this interaction domain are highly conserved throughout evolution within the five LRR‐RLKII clades.     

Chitosan‐triggered immunity to Fusarium in chickpea is associated with changes in the plant extracellular matrix architecture,stomatal closure and remodeling of the plant metabolome and proteome

Pathogen‐/microbe‐associated molecular patterns (PAMPs/MAMPs) initiate complex defense responses by reorganizing the biomolecular dynamics of the host cellular machinery. The extracellular matrix (ECM) acts as a physical scaffold that prevents recognition and entry of phytopathogens, while guard cells perceive and integrate signals metabolically. Although chitosan is a known MAMP implicated in plant defense, the precise mechanism of chitosan‐triggered immunity (CTI) remains unknown. Here, we show how chitosan imparts immunity against fungal disease. Morpho‐histological examination revealed stomatal closure accompanied by reductions in stomatal conductance and transpiration rate as early responses in chitosan‐treated seedlings upon vascular fusariosis. Electron microscopy and Raman spectroscopy showed ECM fortification leading to oligosaccharide signaling, as documented by increased galactose, pectin and associated secondary metabolites. Multiomics approach using quantitative ECM proteomics and metabolomics identified 325 chitosan‐triggered immune‐responsive proteins (CTIRPs), notably novel ECM structural proteins, LYM2 and receptor‐like kinases, and 65 chitosan‐triggered immune‐responsive metabolites (CTIRMs), including sugars, sugar alcohols, fatty alcohols, organic and amino acids. Identified proteins and metabolites are linked to reactive oxygen species (ROS) production, stomatal movement, root nodule development and root architecture coupled with oligosaccharide signaling that leads to Fusarium resistance. The cumulative data demonstrate that ROS, NO and eATP govern CTI, in addition to induction of PR proteins, CAZymes and PAL activities, besides accumulation of phenolic compounds downstream of CTI. The immune‐related correlation network identified functional hubs in the CTI pathway. Altogether, these shifts led to the discovery of chitosan‐responsive networks that cause significant ECM and guard cell remodeling, and translate ECM cues into cell fate decisions during fusariosis.     

Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes

If two related plant species hybridize, their genomes may be combined and duplicated within a single nucleus, thereby forming an allotetraploid. How the emerging plant balances two co‐evolved genomes is still a matter of ongoing research. Here, we focus on satellite DNA (satDNA), the fastest turn‐over sequence class in eukaryotes, aiming to trace its emergence, amplification, and loss during plant speciation and allopolyploidization. As a model, we used Chenopodium quinoa Willd. (quinoa), an allopolyploid crop with 2n = 4x = 36 chromosomes. Quinoa originated by hybridization of an unknown female American Chenopodium diploid (AA genome) with an unknown male Old World diploid species (BB genome), dating back 3.3–6.3 million years. Applying short read clustering to quinoa (AABB), C. pallidicaule (AA), and C. suecicum (BB) whole genome shotgun sequences, we classified their repetitive fractions, and identified and characterized seven satDNA families, together with the 5S rDNA model repeat. We show unequal satDNA amplification (two families) and exclusive occurrence (four families) in the AA and BB diploids by read mapping as well as Southern, genomic, and fluorescent in situ hybridization. Whereas the satDNA distributions support C. suecicum as possible parental species, we were able to exclude C. pallidicaule as progenitor due to unique repeat profiles. Using quinoa long reads and scaffolds, we detected only limited evidence of intergenomic homogenization of satDNA after allopolyploidization, but were able to exclude dispersal of 5S rRNA genes between subgenomes. Our results exemplify the complex route of tandem repeat evolution through Chenopodium speciation and allopolyploidization, and may provide sequence targets for the identification of quinoa's progenitors.     

Ptr1 evolved convergently with RPS2 and Mr5 to mediate recognition of AvrRpt2 in diverse solanaceous species

The Ptr1 (Pseudomonas tomato race 1) locus in Solanum lycopersicoides confers resistance to strains of Pseudomonas syringae pv. tomato expressing AvrRpt2 and Ralstonia pseudosolanacearum expressing RipBN. Here we describe the identification and phylogenetic analysis of the Ptr1 gene. A single recombinant among 585 F2 plants segregating for the Ptr1 locus was discovered that narrowed the Ptr1 candidates to eight nucleotide‐binding leucine‐rich repeat protein (NLR)‐encoding genes. From analysis of the gene models in the S. lycopersicoides genome sequence and RNA‐Seq data, two of the eight genes emerged as the strongest candidates for Ptr1. One of these two candidates was found to encode Ptr1 based on its ability to mediate recognition of AvrRpt2 and RipBN when it was transiently expressed with these effectors in leaves of Nicotiana glutinosa. The ortholog of Ptr1 in tomato and in Solanum pennellii is a pseudogene. However, a functional Ptr1 ortholog exists in Nicotiana benthamiana and potato, and both mediate recognition of AvrRpt2 and RipBN. In apple and Arabidopsis, recognition of AvrRpt2 is mediated by the Mr5 and RPS2 proteins, respectively. Phylogenetic analysis places Ptr1 in a distinct clade compared with Mr5 and RPS2, and it therefore appears to have arisen by convergent evolution for recognition of AvrRpt2.     

Overexpression of the persimmon abscisic acid β‐glucosidase gene (DkBG1) alters fruit ripening in transgenic tomato

β‐Glucosidases (BG) are present in many plant tissues. Among these, abscisic acid (ABA) β‐glucosidases are thought to take part in the adjustment of cellular ABA levels, however the role of ABA‐BG in fruits is still unclear. In this study, through RNA‐seq analysis of persimmon fruit, 10 full‐length DkBG genes were isolated and were all found to be expressed. In particular, DkBG1 was highly expressed in persimmon fruits with a maximum expression 95 days after full bloom (DAFD). We verified that, in vitro, DkBG1 protein can hydrolyze ABA‐glucose ester (ABA‐GE) to release free ABA. Compared with wild‐type, tomato plants that overexpressed DkBG1 significantly upregulated the expression of ABA receptor PYL3/7 genes and showed typical symptoms of ABA hypersensitivity in fruits. DkBG1 overexpression (DkBG1‐OE) accelerated fruit ripening onset by 3–4 days by increasing ABA levels at the pre‐breaker stage and induced early ethylene release compared with wild‐type fruits. DkBG1‐OE altered the expression of ripening regulator NON‐RIPENING (NOR) and its target genes; this in turn altered fruit quality traits such as coloration. Our results demonstrated that DkBG1 plays an important role in fruit ripening and quality by adjusting ABA levels via hydrolysis of ABA‐GE.