MicroRNA Superfamilies Descended from miR390 and Their Roles in Secondary Small Interfering RNA Biogenesis in Eudicots

Trans-acting small interfering RNAs (tasiRNAs) are a major class of small RNAs performing essential biological functions in plants. The first reported tasiRNA pathway, that of miR173-TAS1/2, produces tasiRNAs regulating a set of pentatricopeptide repeat (PPR) genes and has been characterized only in Arabidopsis thaliana to date. Here, we demonstrate that the microRNA (miRNA)-trans-acting small interfering RNA gene (TAS)-pentatricopeptide repeat-containing gene (PPR)-small interfering RNA pathway is a highly dynamic and widespread feature of eudicots. Nine eudicot plants, representing six different plant families, have evolved similar tasiRNA pathways to initiate phased small interfering RNA (phasiRNA) production from PPR genes. The PPR phasiRNA production is triggered by different 22-nucleotide miRNAs, including miR7122, miR1509, and fve-PPRtri1/2, and through distinct mechanistic strategies exploiting miRNA direct targeting or indirect targeting through TAS-like genes (TASL), one-hit or two-hit, or even two layers of tasiRNA–TASL interactions. Intriguingly, although those miRNA triggers display high sequence divergence caused by the occurrence of frequent point mutations and splicing shifts, their corresponding MIRNA genes show pronounced identity to the Arabidopsis MIR173, implying a common origin of this group of miRNAs (super-miR7122). Further analyses reveal that super-miR7122 may have evolved from a newly defined miR4376 superfamily, which probably originated from the widely conserved miR390. The elucidation of this evolutionary path expands our understanding of the course of miRNA evolution, especially for relatively conserved miRNA families.     

Evolutionary Dynamics of the Leucine-Rich Repeat Receptor-Like Kinase (LRR-RLK) Subfamily in Angiosperms

Gene duplications are an important factor in plant evolution, and lineage-specific expanded (LSE) genes are of particular interest. Receptor-like kinases expanded massively in land plants, and leucine-rich repeat receptor-like kinases (LRR-RLK) constitute the largest receptor-like kinases family. Based on the phylogeny of 7,554 LRR-RLK genes from 31 fully sequenced flowering plant genomes, the complex evolutionary dynamics of this family was characterized in depth. We studied the involvement of selection during the expansion of this family among angiosperms. LRR-RLK subgroups harbor extremely contrasting rates of duplication, retention, or loss, and LSE copies are predominantly found in subgroups involved in environmental interactions. Expansion rates also differ significantly depending on the time when rounds of expansion or loss occurred on the angiosperm phylogenetic tree. Finally, using a d_N/d_S-based test in a phylogenetic framework, we searched for selection footprints on LSE and single-copy LRR-RLK genes. Selective constraint appeared to be globally relaxed at LSE genes, and codons under positive selection were detected in 50% of them. Moreover, the leucine-rich repeat domains, and specifically four amino acids in them, were found to be the main targets of positive selection. Here, we provide an extensive overview of the expansion and evolution of this very large gene family.Receptor-like kinases (RLKs) constitute one of the largest gene families in plants and expanded massively in land plants (Embryophyta; Lehti-Shiu et al., 2009, 2012). For plant RLK gene families, the functions of most members are often not known (especially in recently expanded families), but some described functions include innate immunity (Albert et al., 2010), pathogen response (Dodds and Rathjen, 2010), abiotic stress (Yang et al., 2010), development (De Smet et al., 2009), and sometimes multiple functions (Lehti-Shiu et al., 2012). The RLKs usually consist of three domains: an N-terminal extracellular domain, a transmembrane domain, and a C-terminal kinase domain (KD). In plants, the KD usually has a Ser/Thr specificity (Shiu and Bleecker, 2001), but Tyr-specific RLKs were also described (e.g. BRASSINOSTEROID INSENSITIVE1; Oh et al., 2009). Interestingly, it was estimated that approximately 20% of RLKs contain a catalytically inactive KD (e.g. STRUBBELIG and CORYNE; Chevalier et al., 2005; Castells and Casacuberta, 2007; Gish and Clark, 2011). In Arabidopsis (Arabidopsis thaliana), 44 RLK subgroups (SGs) were defined by inferring the phylogenetic relationships between the KDs (Shiu and Bleecker, 2001). Interestingly, different SGs show different duplication/retention rates (Lehti-Shiu et al., 2009). Specifically, RLKs involved in stress responses show a high number of tandemly duplicated genes whereas those involved in development do not (Shiu et al., 2004), which suggests that some RLK genes are important for the responses of land plants to a changing environment (Lehti-Shiu et al., 2012). There seem to be relatively few RLK pseudogenes compared with other large gene families, and copy retention was argued to be driven by both drift and selection (Zou et al., 2009; Lehti-Shiu et al., 2012). As most SGs are relatively old and RLK subfamilies expanded independently in several plant lineages, duplicate retention cannot be explained by drift alone, and natural selection is expected to be an important driving factor in RLK gene family retention (Lehti-Shiu et al., 2009).Leucine-rich repeat-receptor-like kinases (LRR-RLKs), which contain up to 30 leucine-rich repeat (LRRs) in their extracellular domain, constitute the largest RLK family (Shiu and Bleecker, 2001). Based on the KD, 15 LRR-RLK SGs have been established in Arabidopsis (Shiu et al., 2004; Lehti-Shiu et al., 2009). So far, two major functions have been attributed to them: defense against pathogens and development (Tang et al., 2010b). LRR-RLKs involved in defense are predominantly found in lineage-specific expanded (LSE) gene clusters, whereas LRR-RLKs involved in development are mostly found in nonexpanded groups (Tang et al., 2010b). It was also discovered that the LRR domains are significantly less conserved than the remaining domains of the LRR-RLK genes (Tang et al., 2010b). In addition, a study of four plant genomes (Arabidopsis, grape [Vitis vinifera], poplar [Populus trichocarpa], and rice [Oryza sativa]) showed that LRR-RLK genes from LSE gene clusters show significantly more indications of positive selection or relaxed constraint than LRR-RLKs from nonexpanded groups (Tang et al., 2010b).The genomes of flowering plants (angiosperms) have been shown to be highly dynamic compared with most other groups of land plants (Leitch and Leitch, 2012). This dynamic is mostly caused by the frequent multiplication of genetic material, followed by a complex pattern of differential losses (i.e. the fragmentation process) and chromosomal rearrangements (Langham et al., 2004; Leitch and Leitch, 2012). Most angiosperm genomes sequenced so far show evidence for at least one whole-genome multiplication event during their evolution (Jaillon et al., 2007; D’Hont et al., 2012; Tomato Genome Consortium, 2012). At a smaller scale, tandem and segmental duplications are also very common in angiosperms (Arabidopsis Genome Initiative, 2000; International Rice Genome Sequencing Project, 2005; Rizzon et al., 2006). Although the most common fate of duplicated genes is to be progressively lost, in some cases they can be retained in the genome, and adaptive as well as nonadaptive scenarios have been discussed to play a role in this preservation process (for review, see Moore and Purugganan, 2005; Hahn, 2009; Innan, 2009; Innan and Kondrashov, 2010). Whole-genome sequences also revealed that the same gene may undergo several rounds of duplication and retention. These LSE genes were shown to evolve under positive selection more frequently than single-copy genes in angiosperms (Fischer et al., 2014). That study analyzed general trends over whole genomes. Here, we ask if, and to what extent, this trend is observable at LRR-RLK genes. As this gene family is very dynamic and large, and in accordance with the results of Tang et al. (2010b), we expect the effect of positive selection to be even more pronounced than in the whole-genome average.We analyzed 33 Embryophyta genomes to investigate the evolutionary history of the LRR-RLK gene family in a phylogenetic framework. Twenty LRR-RLK SGs were identified, and from this data set, we deciphered the evolutionary dynamics of this family within angiosperms. The expansion/reduction rates were contrasted between SGs and species as well as in ancestral branches of the angiosperm phylogeny. We then focused on genes whose number increased dramatically in an SG- and/or species-specific manner (i.e. LSE genes). Those genes are likely to be involved in species-specific cellular processes or adaptive interactions and were used as a template to infer the potential occurrence of positive selection. This led to the identification of sites at which positive selection likely acted. We discuss our results in the light of angiosperm genome evolution and current knowledge of LRR-RLK functions. Positive selection footprints identified in LSE genes highlight the importance of combining evolutionary analysis and functional knowledge to guide further investigations.     

The Coronatine Toxin of Pseudomonas syringae Is a Multifunctional Suppressor of Arabidopsis Defense

The phytotoxin coronatine (COR) promotes various aspects of Pseudomonas syringae virulence, including invasion through stomata, growth in the apoplast, and induction of disease symptoms. COR is a structural mimic of active jasmonic acid (JA) conjugates. Known activities of COR are mediated through its binding to the F-box–containing JA coreceptor CORONATINE INSENSITIVE1. By analyzing the interaction of P. syringae mutants with Arabidopsis thaliana mutants, we demonstrate that, in the apoplastic space of Arabidopsis, COR is a multifunctional defense suppressor. COR and the critical P. syringae type III effector HopM1 target distinct signaling steps to suppress callose deposition. In addition to its well-documented ability to suppress salicylic acid (SA) signaling, COR suppresses an SA-independent pathway contributing to callose deposition by reducing accumulation of an indole glucosinolate upstream of the activity of the PEN2 myrosinase. COR also suppresses callose deposition and promotes bacterial growth in coi1 mutant plants, indicating that COR may have multiple targets inside plant cells.     

One Divinyl Reductase Reduces the 8-Vinyl Groups in Various Intermediates of Chlorophyll Biosynthesis in a Given Higher Plant Species,But the Isozyme Differs between Species

Divinyl reductase (DVR) converts 8-vinyl groups on various chlorophyll intermediates to ethyl groups, which is indispensable for chlorophyll biosynthesis. To date, five DVR activities have been detected, but adequate evidence of enzymatic assays using purified or recombinant DVR proteins has not been demonstrated, and it is unclear whether one or multiple enzymes catalyze these activities. In this study, we systematically carried out enzymatic assays using four recombinant DVR proteins and five divinyl substrates and then investigated the in vivo accumulation of various chlorophyll intermediates in rice (Oryza sativa), maize (Zea mays), and cucumber (Cucumis sativus). The results demonstrated that both rice and maize DVR proteins can convert all of the five divinyl substrates to corresponding monovinyl compounds, while both cucumber and Arabidopsis (Arabidopsis thaliana) DVR proteins can convert three of them. Meanwhile, the OsDVR (Os03g22780)-inactivated 824ys mutant of rice exclusively accumulated divinyl chlorophylls in its various organs during different developmental stages. Collectively, we conclude that a single DVR with broad substrate specificity is responsible for reducing the 8-vinyl groups of various chlorophyll intermediates in higher plants, but DVR proteins from different species have diverse and differing substrate preferences, although they are homologous.Chlorophyll (Chl) molecules universally exist in photosynthetic organisms. As the main component of the photosynthetic pigments, Chl molecules perform essential processes of absorbing light and transferring the light energy in the reaction center of the photosystems (Fromme et al., 2003). Based on the number of vinyl side chains, Chls are classified into two groups, 3,8-divinyl (DV)-Chl and 3-monovinyl (MV)-Chl. The DV-Chl molecule contains two vinyl groups at positions 3 and 8 of the tetrapyrrole macrocycle, whereas the MV-Chl molecule contains a vinyl group at position 3 and an ethyl group at position 8 of the macrocycle. Almost all of the oxygenic photosynthetic organisms contain MV-Chls, with the exceptions of some marine picophytoplankton species that contain only DV-Chls as their primary photosynthetic pigments (Chisholm et al., 1992; Goericke and Repeta, 1992; Porra, 1997).The classical single-branched Chl biosynthetic pathway proposed by Granick (1950) and modified by Jones (1963) assumed the rapid reduction of the 8-vinyl group of DV-protochlorophyllide (Pchlide) catalyzed by a putative 8-vinyl reductase. Ellsworth and Aronoff (1969) found evidence for both MV and DV forms of several Chl biosynthetic intermediates between magnesium-protoporphyrin IX monomethyl ester (MPE) and Pchlide in Chlorella spp. mutants. Belanger and Rebeiz (1979, 1980) reported that the Pchlide pool of etiolated higher plants contains both MV- and DV-Pchlide. Afterward, following the further detection of MV- and DV-tetrapyrrole intermediates and their biosynthetic interconversion in tissues and extracts of different plants (Belanger and Rebeiz, 1982; Duggan and Rebeiz, 1982; Tripathy and Rebeiz, 1986, 1988; Parham and Rebeiz, 1992, 1995; Kim and Rebeiz, 1996), a multibranched Chl biosynthetic heterogeneity was proposed (Rebeiz et al., 1983, 1986, 1999; Whyte and Griffiths, 1993; Kolossov and Rebeiz, 2010).Biosynthetic heterogeneity refers to the biosynthesis of a particular metabolite by an organelle, tissue, or organism via multiple biosynthetic routes. Varieties of reports lead to the assumption that Chl biosynthetic heterogeneity originates mainly in parallel DV- and MV-Chl biosynthetic routes. These routes are interconnected by 8-vinyl reductases that convert DV-tetrapyrroles to MV-tetrapyrroles by conversion of the vinyl group at position 8 of ring B to the ethyl group (Parham and Rebeiz, 1995; Rebeiz et al., 2003). DV-MPE could be converted to MV-MPE in crude homogenates from etiolated wheat (Triticum aestivum) seedlings (Ellsworth and Hsing, 1974). Exogenous DV-Pchlide could be partially converted to MV-Pchlide in barley (Hordeum vulgare) plastids (Tripathy and Rebeiz, 1988). 8-Vinyl chlorophyllide (Chlide) a reductases in etioplast membranes isolated from etiolated cucumber (Cucumis sativus) cotyledons and barley and maize (Zea mays) leaves were found to be very active in the conversion of exogenous DV-Chlide a to MV-Chlide a (Parham and Rebeiz, 1992, 1995). Kim and Rebeiz (1996) suggested that Chl biosynthetic heterogeneity in higher plants may originate at the level of DV magnesium-protoporphyrin IX (Mg-Proto) and would be mediated by the activity of a putative 8-vinyl Mg-Proto reductase in barley etiochloroplasts and plastid membranes. However, since these reports did not use purified or recombinant enzyme, it is not clear whether the reductions of the 8-vinyl groups of various Chl intermediates are catalyzed by one enzyme of broad specificity or by multiple enzymes of narrow specificity, which actually has become one of the focus issues in Chl biosynthesis.Nagata et al. (2005) and Nakanishi et al. (2005) independently identified the AT5G18660 gene of Arabidopsis (Arabidopsis thaliana) as an 8-vinyl reductase, namely, divinyl reductase (DVR). Chew and Bryant (2007) identified the DVR BciA (CT1063) gene of the green sulfur bacterium Chlorobium tepidum, which is homologous to AT5G18660. An enzymatic assay using a recombinant Arabidopsis DVR (AtDVR) on five DV substrates revealed that the major substrate of AtDVR is DV-Chlide a, while the other four DV substrates could not be converted to corresponding MV compounds (Nagata et al., 2007). Nevertheless, a recombinant BciA is able to reduce the 8-vinyl group of DV-Pchlide to generate MV-Pchlide (Chew and Bryant, 2007). Recently, we identified the rice (Oryza sativa) DVR encoded by Os03g22780 that has sequence similarity with the Arabidopsis DVR gene AT5G18660. We also confirmed that the recombinant rice DVR (OsDVR) is able to not only convert DV-Chlide a to MV-Chlide a but also to convert DV-Chl a to MV-Chl a (Wang et al., 2010). Thus, it is possible that the reductions of the 8-vinyl groups of various Chl biosynthetic intermediates are catalyzed by one enzyme of broad specificity.In this report, we extended our studies to four DVR proteins and five DV substrates. First, ZmDVR and CsDVR genes were isolated from maize and cucumber genomes, respectively, using a homology-based cloning approach. Second, enzymatic assays were systematically carried out using recombinant OsDVR, ZmDVR, CsDVR, and AtDVR as representative DVR proteins and using DV-Chl a, DV-Chlide a, DV-Pchlide a, DV-MPE, and DV-Mg-Proto as DV substrates. Third, we examined the in vivo accumulations of various Chl intermediates in rice, maize, and cucumber. Finally, we systematically investigated the in vivo accumulations of Chl and its various intermediates in the OsDVR (Os03g22780)-inactivated 824ys mutant of rice (Wang et al., 2010). The results strongly suggested that a single DVR protein with broad substrate specificity is responsible for reducing the 8-vinyl groups of various intermediate molecules of Chl biosynthesis in higher plants, but DVR proteins from different species could have diverse and differing substrate preferences even though they are homologous.     

Metabolic Interactions between the Lands Cycle and the Kennedy Pathway of Glycerolipid Synthesis in Arabidopsis Developing Seeds

It has been widely accepted that the primary function of the Lands cycle is to provide a route for acyl remodeling to modify fatty acid (FA) composition of phospholipids derived from the Kennedy pathway. Lysophosphatidylcholine acyltransferase (LPCAT) is an evolutionarily conserved key enzyme in the Lands cycle. In this study, we provide direct evidence that the Arabidopsis thaliana LPCATs, LPCAT1 and LPCAT2, participate in the Lands cycle in developing seeds. In spite of a substantially reduced initial rate of nascent FA incorporation into phosphatidylcholine (PC), the PC level in the double mutant lpcat1 lpcat2-2 remained unchanged. LPCAT deficiency triggered a compensatory response of de novo PC synthesis and a concomitant acceleration of PC turnover that were attributable at least in part to PC deacylation. Acyl-CoA profile analysis revealed complicated metabolic alterations rather than merely reduced acyl group shuffling from PC in the mutant. Shifts in FA stereo-specific distribution in triacylglycerol of the mutant seed suggested a preferential retention of saturated acyl chains at the stereospecific numbering (sn)-1 position from PC and likely a channeling of lysophosphatidic acid, derived from PC, into the Kennedy pathway. Our study thus illustrates an intricate relationship between the Lands cycle and the Kennedy pathway.     

The Arabidopsis Ethylene/Jasmonic Acid-NRT Signaling Module Coordinates Nitrate Reallocation and the Trade-Off between Growth and Environmental Adaptation

Stresses decouple nitrate assimilation and photosynthesis through stress-initiated nitrate allocation to roots (SINAR), which is mediated by the nitrate transporters NRT1.8 and NRT1.5 and functions to promote stress tolerance. However, how SINAR communicates with the environment remains unknown. Here, we present biochemical and genetic evidence demonstrating that in Arabidopsis thaliana, ethylene (ET) and jasmonic acid (JA) affect the crosstalk between SINAR and the environment. Electrophoretic mobility shift assays and chromatin immunoprecipitation assays showed that ethylene response factors (ERFs), including OCTADECANOID-RESPONSIVE ARABIDOPSIS AP2/ERF59, bind to the GCC boxes in the NRT1.8 promoter region, while ETHYLENE INSENSITIVE3 (EIN3) binds to the EIN3 binding site motifs in the NRT1.5 promoter. Genetic assays showed that cadmium and sodium stresses initiated ET/JA signaling, which converged at EIN3/EIN3-Like1 (EIL1) to modulate ERF expression and hence to upregulate NRT1.8. By contrast, ET and JA signaling mediated the downregulation of NRT1.5 via EIN3/EIL1 and other, unknown component(s). SINAR enhanced stress tolerance and decreased plant growth under nonstressed conditions through the ET/JA-NRT1.5/NRT1.8 signaling module. Interestingly, when nitrate reductase was impaired, SINAR failed to affect either stress tolerance or plant growth. These data suggest that SINAR responds to environmental conditions through the ET/JA-NRT signaling module, which further modulates stress tolerance and plant growth in a nitrate reductase-dependent manner.