Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Diversified quantity, gene structure, and expression profile of OPR gene family of A. annua

Artemisia annua, the source of artemisinin production, is a traditional herb used for treating malaria for thousand years. The genetic background is of high heterozygosity and traits (plant height, biomass, artemisinin content, etc.) are diverse across different germplasms. Unraveling the key genes associated with growth and secondary metabolism is essential for the efficient production of artemisinin. The 12-oxo-phytodienoic acid reductase (OPR) genes, crucial for plant growth and development and stress resistance, remain unexplored in A. annua. In this study, nine OPR genes (named as AaOPR1 to AaOPR9) were identified in A. annua, including two pairs of genes formed from recent tandem duplications. The number of OPRs varied among different haplotype genomes, and each OPR gene exhibiting distinct expression pattern. Moreover, the OPR family displayed evolutionarily activity with significant variations in numbers and gene structures observed across different plant species. Widespread gene duplication of OPRs, observed in the majority of analyzed plant genomes, brought evolutionary potential. DBR2, a member of AaOPRs involved in artemisinin biosynthesis, had two copies (AaOPR1/DBR2.1 and AaOPR2/DBR2.2) with different expression patterns, one of which was a recently generated copy with a significant 7-amino acids truncation. Heterologous protein expression and functional characterization of the two copies of DBR2 yielded multiple isomers with identical molecular weights but different arrangements, indicating neofunctionalization of the newly generated copy. The polymorphism within the OPR gene family merely scratches the surface of the genetic diversity in A. annua, and further investigation of genetic features is needed for the screening of elite germplasm resources.

PnOPR6 from Antarctic moss mediates JA-ABA crosstalk and enhances abiotic stress tolerance in Arabidopsis

Jasmonates (JAs) and abscisic acid (ABA) are vital plant hormones that are integral to the plant's response mechanisms against various abiotic stresses. These hormones also function in an antagonistic manner to regulate seed germination and dormancy. However, little is known about the molecular mechanism underlying the interaction between ABA and JA signaling. Here, seven 12-oxo-phytodienoic acid reductase genes (PnOPR1-7), a key enzyme in the JA biosynthesis pathway, were identified in the Antarctic moss Pohlia nutans transcriptome, and their expressions in response to abiotic stress were examined. Among these, PnOPR6 expression levels rose most under cold and UV-B stresses. Transgenic Arabidopsis overexpressing PnOPR6 demonstrated increased tolerance to salt, cold, dehydration, glucose, and ABA, but also greater sensitivity to methyl jasmonate (MeJA) during seed germination or early root growth. Furthermore, in the transgenic Arabidopsis, PnOPR6 suppressed the expression of genes involved in the ABA pathway and ABI3/5-responsive JA receptor COI1. Additionally, phytohormone metabolomics investigations revealed a significant rise in JA precursor (OPDA, OPC-6, and OPC-4), JA, and its derivative 12-OH-JA in PnOPR6-overexpressing line. Moreover, the accumulation of flavonoid in Arabidopsis was enhanced by heterologous expression of PnOPR6. These findings imply that PnOPR6 functions as a signaling regulator, improving plant resistance to abiotic stress through flavonoid accumulation and JA-ABA antagonistic crosstalk, therefore aiding P. nutans in adjusting to polar climates.

Genome-wide analysis of OPR family genes in Vitis vinifera and the role of VvOPR1 in copper, zinc tolerance

12-oxo-phytodienoic acid reductase (OPR) is one of the key enzymes in the octadecanoid pathway, and it controls the last step of jasmonic acid (JA) biosynthesis. Although multiple isoforms and functions of OPRs have been identified in various plants, no OPR genes have been identified, and their possible roles in grapevine development and defense mechanisms remain unknown. In this study, nine VvOPR genes were identified from grapevine genome and classified into two subfamilies. Systematic analyses of the physical and chemical properties, the expression and structure of the VvOPR genes, promoter elements, and chromosome locations were performed via bioinformatics and molecular biology methods. In addition, we described the characterization of the OPRI gene VvOPR1, which was synthesized via a PCR-based two-step DNA synthesis quantification reverse-transcription (PTDS) method. VvOPR1 expression is tissue-specific and induced by various stresses. The overexpression of VvOPR1 in Arabidopsis and rice (OT) significantly increased tolerance to Cu, Zn stress, and Cu, Zn stress-induced restriction of the germination rate, root/shoot length and fresh weight was significantly alleviated in OT. In OT, VvOPR1 enhanced the photosynthetic capacity, promoted ABA synthesis and the ABA-dependent stress response pathway, improved the antioxidation capacity by increasing the activities of ROS scavengers and the expression level of the related genes, while enhancing the accumulation of proline, AsA, GSH and reducing MDA and H2O2 levels. Moreover, VvOPR1 reduced Cu2+, Zn2+ accumulation and translocation. Together, we first systematically characterized the grapevine OPR gene family and reported that VvOPR1 responded to Cu, Zn stress in an ABA-dependent manner, and was quite independent of JA synthesis and signaling. All of the above results provide an important research basis and theoretical basis for further revealing the functions of VvOPR in grapevines in the future.

Genome-wide identification, gene cloning, subcellular location and expression analysis of the OPR gene family under salt stress in sweetpotato

BACKGROUND

The 12-oxo-phytodienoic acid reductase (OPR) enzyme is crucial for the synthesis of jasmonates (JAs), and is involved in the plant stress response. However, the OPR gene family in sweetpotato, an important horticultural crop, remains unidentified.

RESULTS

In this study, we employed bioinformatics techniques to identify nine IbOPR genes. Phylogenetic analysis revealed that these genes could be divided into Group I and Group II. Synteny analysis indicated that IbOPR evolution was driven by tandem duplication, whole-genome duplication (WGD), and segmental duplication events. The promoter sequences of IbOPRs were found to be associated with stress and hormonal responses. Additionally, we successfully cloned four IbOPRs from "Haida HD7791" and "Haida HD7798" using homologous cloning technology. These sequences were 1203 bp, 1200 bp, 1134 bp, and 1137 bp in length and encoded 400, 399, 377, and 378 amino acids, respectively. The protein sequence similarity between the salt-tolerant variety "Haida HD7791" and the salt-sensitive variety "Haida HD7798" was determined to be 96.75% for IbOPR2, 99.75% for IbOPR3, 92.06% for IbOPR6, and 98.68% for IbOPR7. Phylogenetic analysis categorized IbOPR2 and IbOPR3 proteins into Group II, while IbOPR6 and IbOPR7 proteins belonged to Group I. Subcellular localization experiments showed IbOPR2 protein present in the peroxisome, while IbOPR3, IbOPR6, and IbOPR7 proteins were found in the cytoplasm and nucleus. Salt stress induction experiments demonstrated that IbOPR2, IbOPR3, and IbOPR7 were significantly upregulated only in 'Haida HD7791' after 6 h. In contrast, IbOPR6 was induced in 'Haida HD7798' at 6 h but inhibited in 'Haida HD7791' at later time points (12, 24, 48, and 72 h), highlighting functional differences in salt stress responses.

CONCLUSIONS

Our findings suggest that IbOPR2 may play a crucial role in sweetpotato's response to salt stress by participating in JAs synthesis. These results provide a foundation for future functional analyses of OPR genes in sweetpotato.

Genome-wide identification and expression analysis of ferric reductase oxidase (FRO) genes in Gossypium spp. reveal their crucial role in iron homeostasis under abiotic and biotic stress

Ferric Reductase Oxidase (FRO) genes are pivotal in iron uptake and homeostasis in plants, yet they are not studied in cotton. Here, we identify and analyze 65 FRO homologs (21 GhFRO, 21 GbFRO, 11 GaFRO, 12 GrFRO) across four Gossypium species (G. hirsutum, G. barbadense, G. arboreum, G. raimondii). FRO exhibit conserved ferric reductase activity and conserved domain structures; Ferric_reduct (PF01794), FAD_binding_8 (PF08022), and NAD_binding_6 (PF08030) across species. Physicochemical properties and subcellular localization analysis provided insights into FRO proteins' functional characteristics, mainly localized to the plasma membrane. Phylogenetic analysis delineates 11 groups, indicating both conserved and divergent evolutionary patterns. Gene structure analysis unveils varying exon-intron compositions. Chromosomal localization shows distribution across A and D genomes, suggesting evolutionary dynamics. Synteny analysis reveals paralogous and orthologous gene pairs subjected to purifying selection. The cis-regulatory elements analysis implicates diverse regulatory mechanisms. Expression profiling highlights dynamic regulation across developmental stages, abiotic and biotic stress conditions. GhFRO interacts with Ca++-dependent protein kinases-10/28-like (CDPKs10/28-like) and metal transporter Natural resistance-associated macrophage protein 6 (Nramp6) to regulate metal ion transport and iron homeostasis. The three-dimensional protein structure prediction suggests potential ligand-binding sites in FRO proteins. Moreover, qRT-PCR analysis of selected eight GhFROs in leaves treated with stress elicitors, MeJA, SA, NaCl, and PEG for 1h, 2h, 4h, and 6h revealed significant downregulation. Overall, this comprehensive study provides insights into FRO gene diversity, evolution, structure, regulation, and function in cotton, with implications for understanding plant iron homeostasis and stress responses.

Multiomics Combined with Expression Pattern Analysis Reveals the Regulatory Response of Key Genes in Potato Jasmonic Acid Signaling Pathways to Cadmium Stress

Jasmonic acid (JA) is an endogenous phytohormone that regulates plant physiological metabolism and stress response processes, either independently or through hormone crosstalk. Our phytohormone assay and transcriptome-metabolome analysis revealed the key genes and metabolites involved in the JA pathway in response to 0-250 μM cadmium (Cd) in potato seedlings. Transcriptome gene set enrichment and gene ontology analysis indicated that JA-related genes were significantly enriched. Specifically, members from the StOPR and StJAZ gene families showed pronounced responses to Cd stress and methyl jasmonate treatment. As a negative regulatory transcription factor of the JA signaling pathway, StJAZ14 exhibited a decreasing trend under Cd stress. Yeast two-hybrid assay identified an interaction between StJAZ14 and StBZR1, which is located on the brassinolide pathway. In addition to unveiling the critical role of the JA pathway in regulating potato response to Cd stress, the functional mechanism was preliminarily explored.

Inactivation of a Wheat Ribosomal Silencing Factor Gene TaRsfS Confers Resistance to Both Powdery Mildew and Stripe Rust

Powdery mildew and stripe rust are major diseases on wheat worldwide that cause significant reductions in wheat production. The ribosomal silencing factor (RsfS) has been proven to regulate protein biosynthesis by inhibiting the translation process in bacterial response to stress. However, the role of RsfS in plant resistance to biotic stresses remains unclear. In this study, the RsfS homolog, TaRsfS was isolated from wheat. Overexpression of TaRsfS (TaRsfS-OE) reduces wheat resistance to powdery mildew and stripe rust and TaRsfS knockout (TaRsfS-KO) increases wheat resistance to both diseases without affecting key agronomic traits. The interaction protein of TaRsfS, 12-oxo-phytodienoic acid reductase 1 (TaOPR1), a key enzyme in the biosynthesis of jasmonic acid (JA), was screened and identified. Knocking-down and overexpression of TaOPR1 indicated that TaOPR1 positively regulates wheat resistance to powdery mildew and stripe rust. TaRsfS may regulate TaOPR1 at upstream, bind to the enzyme activity pocket of TaOPR1 and affect TaOPR1 enzyme activity, resulting in a reduced JA biosynthesis and wheat susceptible to powdery mildew and stripe rust. Collectively, TaRsfS is a susceptibility gene and negatively regulates wheat resistance to powdery mildew and stripe rust, and it has good potential for improving wheat resistance by genetic modifications.

Evolutionarily conserved 12-oxophytodienoate reductase trans-lncRNA pair affects disease resistance in tea (Camellia sinensis) via the jasmonic acid signaling pathway

Long non-coding RNAs (lncRNAs) have gathered significant attention due to their pivotal role in plant growth, development, and biotic and abiotic stress resistance. Despite this, there is still little understanding regarding the functions of lncRNA in these domains in the tea plant (Camellia sinensis), mainly attributable to the insufficiencies in gene manipulation techniques for tea plants. In this study, we designed a novel strategy to identify evolutionarily conserved trans-lncRNA (ECT-lncRNA) pairs in plants. We used highly consistent base sequences in the exon-overlapping region between trans-lncRNAs and their target gene transcripts. Based on this method, we successfully screened 24 ECT-lncRNA pairs from at least two or more plant species. In tea, as observed in model plants such as Arabidopsis, alfalfa, potatoes, and rice, there exists a trans-lncRNA capable of forming an ECT-lncRNA pair with transcripts of the 12-oxophytodienoate reductase (OPR) family, denoted as the OPRL/OPR pair. Considering evolutionary perspectives, the OPRL gene cluster in each species likely originates from a replication event of the OPR gene cluster. Gene manipulation and gene expression analysis revealed that CsOPRL influences disease resistance by regulating CsOPR expression in tea plants. Furthermore, the knockout of StOPRL1 in Solanum tuberosum led to aberrant growth characteristics and strong resistance to fungal infection. This study provides insights into a strategy for the screening and functional verification of ECT-lncRNA pairs.

A peroxisomal cinnamate:CoA ligase-dependent phytohormone metabolic cascade in submerged rice germination

The mechanism underlying the ability of rice to germinate underwater is a largely enigmatic but key research question highly relevant to rice cultivation. Moreover, although rice is known to accumulate salicylic acid (SA), SA biosynthesis is poorly defined, and its role in underwater germination is unknown. It is also unclear whether peroxisomes, organelles essential to oilseed germination and rice SA accumulation, play a role in rice germination. Here, we show that submerged imbibition of rice seeds induces SA accumulation to promote germination in submergence. Two submergence-induced peroxisomal Oryza sativa cinnamate:CoA ligases (OsCNLs) are required for this SA accumulation. SA exerts this germination-promoting function by inducing indole-acetic acid (IAA) catabolism through the IAA-amino acid conjugating enzyme GH3. The metabolic cascade we identified may potentially be adopted in agriculture to improve the underwater germination of submergence-intolerant rice varieties. SA pretreatment is also a promising strategy to improve submerged rice germination in the field.

Functional analysis of a rice 12-oxo-phytodienoic acid reductase gene (OsOPR1) involved in Cd stress tolerance

BACKGROUND

The accumulation of cadmium (Cd) in plants may compromise the growth and development of plants, thereby endangering human health through the food chain. Understanding how plants respond to Cd is important for breeding low-Cd rice cultivars.

METHODS

In this study, the functions of 12-oxo-phytodienoic acid reductase 1 (OsOPR1) were predicted through bioinformatics analysis. The expression levels of OsOPR1 under Cd stress were analyzed by using qRT-PCR. Then, the role that OsOPR1 gene plays in Cd tolerance was studied in Cd-sensitive yeast strain (ycf1), and the Cd concentration of transgenic yeast was analyzed using inductively coupled plasma mass spectrometry (ICP-MS).

RESULTS

Bioinformatics analysis revealed that OsOPR1 was a protein with an Old yellow enzyme-like FMN (OYE_like_FMN) domain, and the cis-acting elements which regulate hormone synthesis or responding abiotic stress were abundant in the promoter region, which suggested that OsOPR1 may exhibit multifaceted biological functions. The expression pattern analysis showed that the expression levels of OsOPR1 were induced by Cd stress both in roots and roots of rice plants. However, the induced expression of OsOPR1 by Cd was more significant in the roots compared to that in roots. In addition, the overexpression of OsOPR1 improved the Cd tolerance of yeast cells by affecting the expression of antioxidant enzyme related genes and reducing Cd content in yeast cells.

CONCLUSION

Overall, these results suggested that OsOPR1 is a Cd-responsive gene and may has a potential for breeding low-Cd or Cd-tolerant rice cultivars and for phytoremediation of Cd-contaminated in farmland.

Identification of the MYB gene family in Sorghum bicolor and functional analysis of SbMYBAS1 in response to salt stress

Salt stress adversely affects plant growth and development. It is necessary to understand the underlying salt response mechanism to improve salt tolerance in plants. MYB transcription factors can regulate plant responses to salt stress. However, only a few studies have explored the role of MYB TFs in Sorghum bicolor (L.) Moench. So we decided to make a systematic analysis and research on the sorghum MYB family. A total of 210 MYB genes in sorghum were identified in this study. Furthermore, 210 MYB genes were distributed across ten chromosomes, named SbMYB1-SbMYB210. To study the phylogeny of the identified TFs, 210 MYB genes were divided into six subfamilies. We further demonstrated that SbMYB genes have evolved under strong purifying selection. SbMYBAS1 (SbMYB119) was chosen as the study object, which the expression decreased under salt stress conditions. Further study of the SbMYBAS1 showed that SbMYBAS1 is located in the nucleus. Under salt stress conditions, Arabidopsis plants overexpressed SbMYBAS1 showed significantly lower dry/fresh weight and chlorophyll content but significantly higher membrane permeability, MDA content, and Na+/K+ ratio than the wild-type Arabidopsis plants. Yeast two-hybrid screening result showed that SbMYBAS1 might interact with proteins encoded by SORBI_302G184600, SORBI_3009G247900 and SORBI_3004G59600. Results also showed that SbMYBAS1 could regulate the expression of AtGSTU17, AtGSTU16, AtP5CS2, AtUGT88A1, AtUGT85A2, AtOPR2 and AtPCR2 under salt stress conditions. This work laid a foundation for the study of the response mechanism of sorghum MYB gene family to salt stress.

Genomic insight into domestication of rubber tree

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

Axenic in vitro cultivation and genome diploidization of the moss Vesicularia montagnei for horticulture utilization

Mosses are widely used in the establishment of greenery. However, little research has been conducted to choose a suitable species or improve their performance for this application. In our study, we examined Vesicularia montagnei (V. montagnei), a robust moss that is widely distributed in temperate, subtropical, and tropical Asia with varying environmental conditions. Axenic cultivation system of V. montagnei was developed on modified BCD medium, which enabled its propagation and multiplication in vitro. In this axenic cultivation environment, several diploid V. montagnei lines with enhancement of rhizoid system were generated through artificial induction of diploidization. Transcriptomic analysis revealed that several genes responsible for jasmonic acid (JA) biosynthesis and signaling showed significant higher expression levels in the diploid lines compared to the wild type. These results are consistent with the increasement of JA content in the diploid lines. Our establishment of the axenic cultivation method may provide useful information for further study of other Vesicularia species. The diploid V. montagnei lines with improved rhizoid system may hold promising potential for greenery applications. Additionally, our study sheds light on the biosynthesis and functions of JA in the early landed plants.

Dosage differences in 12-OXOPHYTODIENOATE REDUCTASE genes modulate wheat root growth

Wheat, an essential crop for global food security, is well adapted to a wide variety of soils. However, the gene networks shaping different root architectures remain poorly understood. We report here that dosage differences in a cluster of monocot-specific 12-OXOPHYTODIENOATE REDUCTASE genes from subfamily III (OPRIII) modulate key differences in wheat root architecture, which are associated with grain yield under water-limited conditions. Wheat plants with loss-of-function mutations in OPRIII show longer seminal roots, whereas increased OPRIII dosage or transgenic over-expression result in reduced seminal root growth, precocious development of lateral roots and increased jasmonic acid (JA and JA-Ile). Pharmacological inhibition of JA-biosynthesis abolishes root length differences, consistent with a JA-mediated mechanism. Transcriptome analyses of transgenic and wild-type lines show significant enriched JA-biosynthetic and reactive oxygen species (ROS) pathways, which parallel changes in ROS distribution. OPRIII genes provide a useful entry point to engineer root architecture in wheat and other cereals.

The genome wide analysis of Tryptophan Aminotransferase Related gene family, and their relationship with related agronomic traits in Brassica napus

Tryptophan Aminotransferase of Arabidopsis1/Tryptophan Aminotransferase-Related (TAA1/TAR) proteins are the enzymes that involved in auxin biosynthesis pathway. The TAA1/TAR gene family has been systematically characterized in several plants but has not been well reported in Brassica napus. In the present study, a total of 102 BnTAR genes with different number of introns were identified. It was revealed that these genes are distributed unevenly and occurred as clusters on different chromosomes except for A4, A5, A10 and C4 in B. napus. Most of the these BnTAR genes are conserved despite of existing of gene loss and gene gain. In addition, the segmental replication and whole-genome replication events were both play an important role in the BnTAR gene family formation. Expression profiles analysis indicated that the expression of BnTAR gene showed two patterns, part of them were mainly expressed in roots, stems and leaves of vegetative organs, and the others were mainly expressed in flowers and seeds of reproductive organs. Further analysis showed that many of BnTAR genes were located in QTL intervals of oil content or seed weight, for example BnAMI10 was located in cqOC-C5-4 and cqSW-A2-2, it indicated that some of the BnTAR genes might have relationship with these two characteristics. This study provides a multidimensional analysis of the TAA1/TAR gene family and a new insight into its biological function in B. napus.

Identification of the 12-oxo-phytoeienoic acid reductase (OPR) gene family in pepper (Capsicum annuum L.) and functional characterization of CaOPR6 in pepper fruit development and stress response

The 12-oxo-phytoeienoic acid reductase (OPR) is a kind of enzyme in octadecanoid biosynthesis pathway, which determines the biosynthesis of jasmonic acid. Although the roles of OPRs have been extensively studied in several crop plants, little is known about the biological functions of OPR encoding genes in Capsicum annuum plants. In this study, seven OPR family genes (CaOPR1-7) were identified from the C. annuum genome. The physical and chemical properties of CaOPR1-7 were further analyzed, including gene expression patterns, promoter elements and chromosomal locations. The results showed that the seven CaOPR homologous could be divided into two subgroups, and CaOPR6 was highly similar to AtOPR3 in Arabidopsis. The expression of CaOPR6 was significantly induced by various stresses such as cold, salt and pathogen infection, indicating that CaOPR6 plays important roles in response to abiotic and biotic stresses. Overall, these findings improve the understanding of the biological functions of CaOPR6 in the development of pepper fruit and stress response of pepper plants, and facilitate further studies on the molecular biology of OPR proteins in Solanaceae vegetables.

Plant protein-coding gene families: Their origin and evolution

Steady advances in genome sequencing methods have provided valuable insights into the evolutionary processes of several gene families in plants. At the core of plant biodiversity is an extensive genetic diversity with functional divergence and expansion of genes across gene families, representing unique phenomena. The evolution of gene families underpins the evolutionary history and development of plants and is the subject of this review. We discuss the implications of the molecular evolution of gene families in plants, as well as the potential contributions, challenges, and strategies associated with investigating phenotypic alterations to explain the origin of plants and their tolerance to environmental stresses.

Phytohormone biosynthesis and signaling pathways of mosses

Most known phytohormones regulate moss development. We present a comprehensive view of the synthesis and signaling pathways for the most investigated of these compounds in mosses, focusing on the model Physcomitrium patens. The last 50 years of research have shown that most of the known phytohormones are synthesized by the model moss Physcomitrium patens (formerly Physcomitrella patens) and regulate its development, in interaction with responses to biotic and abiotic stresses. Biosynthesis and signaling pathways are best described in P. patens for the three classical hormones auxins, cytokinins and abscisic acid. Furthermore, their roles in almost all steps of development, from early filament growth to gametophore development and sexual reproduction, have been the focus of much research effort over the years. Evidence of hormonal roles exist for ethylene and for CLE signaling peptides, as well as for salicylic acid, although their possible effects on development remain unclear. Production of brassinosteroids by P. patens is still debated, and modes of action for these compounds are even less known. Gibberellin biosynthesis and signaling may have been lost in P. patens, while gibberellin precursors such as ent-kaurene derivatives could be used as signals in a yet to discover pathway. As for jasmonic acid, it is not used per se as a hormone in P. patens, but its precursor OPDA appears to play a corresponding role in defense against abiotic stress. We have tried to gather a comprehensive view of the biosynthesis and signaling pathways for all these compounds in mosses, without forgetting strigolactones, the last class of plant hormones to be reported. Study of the strigolactone response in P. patens points to a novel signaling compound, the KAI2-ligand, which was likely employed as a hormone prior to land plant emergence.

Genome-Wide Identification, Genomic Organization, and Characterization of Potassium Transport-Related Genes in Cajanus cajan and Their Role in Abiotic Stress

Potassium is the most important and abundant inorganic cation in plants and it can comprise up to 10% of a plant's dry weight. Plants possess complex systems of transporters and channels for the transport of K+ from soil to numerous parts of plants. Cajanus cajan is cultivated in different regions of the world as an economical source of carbohydrates, fiber, proteins, and fodder for animals. In the current study, 39 K+ transport genes were identified in C. cajan, including 25 K+ transporters (17 carrier-like K+ transporters (KUP/HAK/KTs), 2 high-affinity potassium transporters (HKTs), and 6 K+ efflux transporters (KEAs) and 14 K+ channels (9 shakers and 5 tandem-pore K+ channels (TPKs). Chromosomal mapping indicated that these genes were randomly distributed among 10 chromosomes. A comparative phylogenetic analysis including protein sequences from Glycine max, Arabidopsis thaliana, Oryza sativa, Medicago truncatula Cicer arietinum, and C. cajan suggested vital conservation of K+ transport genes. Gene structure analysis showed that the intron/exon organization of K+ transporter and channel genes is highly conserved in a family-specific manner. In the promoter region, many cis-regulatory elements were identified related to abiotic stress, suggesting their role in abiotic stress response. Abiotic stresses (salt, heat, and drought) adversely affect chlorophyll, carotenoids contents, and total soluble proteins. Furthermore, the activities of catalase, superoxide, and peroxidase were altered in C. cajan leaves under applied stresses. Expression analysis (RNA-seq data and quantitative real-time PCR) revealed that several K+ transport genes were expressed in abiotic stress-responsive manners. The present study provides an in-depth understanding of K+ transport system genes in C. cajan and serves as a basis for further characterization of these genes.

The OPR gene family in watermelon: Genome-wide identification and expression profiling under hormone treatments and root-knot nematode infection

The enzyme 12-oxo-phytodienoic acid reductase (OPR) is important in the jasmonic acid (JA) biosynthesis pathway and thus plays a vital role in plant defence. However, systematic and comprehensive analyses of OPR genes in watermelon and their roles in defence responses are extremely limited. The physicochemical properties, phylogenetic tree, gene structure and cis-acting elements of watermelon OPR genes were analysed using bioinformatics, and qRT-PCR and RNA-Seq were applied to assay expression of OPR genes in watermelon. A total of five OPR family genes were identified in watermelon, which were unevenly distributed across the four chromosomes. Phylogenetic analysis assigned OPR members from different plant species to five subfamilies (OPRI-OPRV). The motif compositions of OPR members were relatively conserved. Expression analysis using qRT-PCR revealed that ClOPR genes, except for ClOPR5, were highly expressed in the flower and fruit. RNA-seq analysis showed that the ClOPR genes had different expression patterns during flesh and rind development. Furthermore, the ClOPR genes, particularly ClOPR2 and ClOPR4, were significantly upregulated by exogenous JA, salicylic acid (SA) and ethylene (ET) treatments. In addition, red light induced expression of ClOPR2 and ClOPR4 in leaves and roots of root-knot nematode (RKN)-infected watermelon plants, suggesting their involvement in red light-induced defence against RKN. These results provide a theoretical basis for elucidating the diverse functions of OPR family genes in watermelon.

Genome-wide identification and characterization of the OFP gene family in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Ovate family proteins (OFPs) are a class of proteins with a conserved OVATE domain that contains approximately 70 amino acid residues. OFP proteins are plant-specific transcription factors that participate in regulating plant growth and development and are widely distributed in many plants. Little is known about OFPs in Brassica rapa to date. We identified 29 OFP genes in Brassica rapa and found that they were unevenly distributed on 10 chromosomes. Intron gain events may have occurred during the structural evolution of BraOFP paralogues. Syntenic analysis verified Brassica genome triplication, and whole genome duplication likely contributed to the expansion of the OFP gene family. All BraOFP genes had light responsive- and phytohormone-related cis-acting elements. Expression analysis from RNA-Seq data indicated that there were obvious changes in the expression levels of six OFP genes in the Brassica rapa hybrid, which may contribute to the formation of heterosis. Finally, we found that the paralogous genes had different expression patterns among the hybrid and its parents. These results provide the theoretical basis for the further analysis of the biological functions of OFP genes across the Brassica species.

Comprehensive analysis of polygalacturonase genes offers new insights into their origin and functional evolution in land plants

Polygalacturonase (PG) is a hydrolase that participates in pectin degradation, pod shattering and fruit softening. Here, we identified 2786 PG genes across 54 plants, which could be divided into three groups. Evolutionary analysis suggested that PG family originated from the charophyte green algae, and Subgroups A2-A4 evolved from the Subgroup A1 after the tracheophyte-angiosperm split. Whole-genome duplication was the major force leading to PG gene expansion. Interestingly, the PG genes continuously expanded in eudicots, whereas it contracted in monocots after the eudicot-monocot split. PG genes in Group A are expressed at high levels in floral organs, whereas genes in Groups B and C are expressed at high levels in various tissues. Moreover, three BnaPG15 members were found for their potential possibility in pod shattering in Brassica napus. Our results provide new insight into the evolutionary history of PG family, and their potentially functional role in plants.

Genome-Wide Analysis of OPR Family Genes in Cotton Identified a Role for GhOPR9 in Verticillium dahliae Resistance

The 12-oxo-phytodienoic acid reductases (OPRs) have been proven to play a major role in plant development and growth. Although the classification and functions of OPRs have been well understood in Arabidopsis, tomato, rice, maize, and wheat, the information of OPR genes in cotton genome and their responses to biotic and abiotic stresses have not been reported. In this study, we found 10 and 9 OPR genes in Gossypium hirsutum and Gossypium barbadense, respectively. They were classified into three groups, based on the similar gene structure and conserved protein motifs. These OPR genes just located on chromosome 01, chromosome 05, and chromosome 06. In addition, the whole genome duplication (WGD) or segmental duplication events contributed to the evolution of the OPR gene family. The analyses of cis-acting regulatory elements of GhOPRs showed that the functions of OPR genes in cotton might be related to growth, development, hormone, and stresses. Expression patterns showed that GhOPRs were upregulated under salt treatment and repressed by polyethylene glycol 6000 (PEG6000). The expression patterns of GhOPRs were different in leaf, root, and stem under V. dahliae infection. GhOPR9 showed a higher expression level than other OPR genes in cotton root. The virus-induced gene silencing (VIGS) analysis suggested that knockdown of GhOPR9 could increase the susceptibility of cotton to V. dahliae infection. Furthermore, GhOPR9 also modulated the expressions of jasmonic acid (JA) pathway-regulated genes under the V. dahliae infection. Overall, our results provided the evolution and potential functions of the OPR genes in cotton. These findings suggested that GhOPR9 might play an important role in cotton resistance to V. dahliae.

Genome wide survey, evolution and expression analysis of PHD finger genes reveal their diverse roles during the development and abiotic stress responses in Brassica rapa L

BACKGROUND

Plant homeodomain (PHD) finger proteins are widely present in all eukaryotes and play important roles in chromatin remodeling and transcriptional regulation. The PHD finger can specifically bind a number of histone modifications as an "epigenome reader", and mediate the activation or repression of underlying genes. Many PHD finger genes have been characterized in animals, but only few studies were conducted on plant PHD finger genes to this day. Brassica rapa (AA, 2n = 20) is an economically important vegetal, oilseed and fodder crop, and also a good model crop for functional and evolutionary studies of important gene families among Brassica species due to its close relationship to Arabidopsis thaliana.

RESULTS

We identified a total of 145 putative PHD finger proteins containing 233 PHD domains from the current version of B. rapa genome database. Gene ontology analysis showed that 67.7% of them were predicted to be located in nucleus, and 91.3% were predicted to be involved in protein binding activity. Phylogenetic, gene structure, and additional domain analyses clustered them into different groups and subgroups, reflecting their diverse functional roles during plant growth and development. Chromosomal location analysis showed that they were unevenly distributed on the 10 B. rapa chromosomes. Expression analysis from RNA-Seq data showed that 55.7% of them were constitutively expressed in all the tested tissues or organs with relatively higher expression levels reflecting their important housekeeping roles in plant growth and development, while several other members were identified as preferentially expressed in specific tissues or organs. Expression analysis of a subset of 18 B. rapa PHD finger genes under drought and salt stresses showed that all these tested members were responsive to the two abiotic stress treatments.

CONCLUSIONS

Our results reveal that the PHD finger genes play diverse roles in plant growth and development, and can serve as a source of candidate genes for genetic engineering and improvement of Brassica crops against abiotic stresses. This study provides valuable information and lays the foundation for further functional determination of PHD finger genes across the Brassica species.

Genome-Wide Identification and Characterization of the Vacuolar H+-ATPase Subunit H Gene Family in Crop Plants

The vacuolar H⁺-ATPase (V-ATPase) plays many important roles in cell growth and in response to stresses in plants. The V-ATPase subunit H (VHA-H) is required to form a stable and active V-ATPase. Genome-wide analyses of VHA-H genes in crops contribute significantly to a systematic understanding of their functions. A total of 22 VHA-H genes were identified from 11 plants representing major crops including cotton, rice, millet, sorghum, rapeseed, maize, wheat, soybean, barley, potato, and beet. All of these VHA-H genes shared exon-intron structures similar to those of Arabidopsis thaliana. The C-terminal domain of VHA-H was shorter and more conserved than the N-terminal domain. The VHA-H gene was effectively used as a genetic marker to infer the phylogenetic relationships among plants, which were congruent with currently accepted taxonomic groupings. The VHA-H genes from six species of crops (Gossypium raimondii, Brassica napus, Glycine max, Solanum tuberosum, Triticum aestivum, and Zea mays) showed high gene structural diversity. This resulted from the gains and losses of introns. Seven VHA-H genes in six species of crops (Gossypium raimondii, Hordeum vulgare, Solanum tuberosum, Setaria italica, Triticum aestivum, and Zea mays) contained multiple transcript isoforms arising from alternative splicing. The study of cis-acting elements of gene promoters and RNA-seq gene expression patterns confirms the role of VHA-H genes as eco-enzymes. The gene structural diversity and proteomic diversity of VHA-H genes in our crop sampling facilitate understanding of their functional diversity, including stress responses and traits important for crop improvement.

Expansion and Functional Divergence of Inositol Polyphosphate 5-Phosphatases in Angiosperms

Inositol polyphosphate 5-phosphatase (5PTase), a key enzyme that hydrolyzes the 5` position of the inositol ring, has essential functions in growth, development, and stress responses in plants, yeasts, and animals. However, the evolutionary history and patterns of 5PTases have not been examined systematically. Here, we report a comprehensive molecular evolutionary analysis of the 5PTase gene family and define four groups. These four groups are different from former classifications, which were based on in vitro substrate specificity. Most orthologous groups appear to be conserved as single or low-copy genes in all lineages in Groups II-IV, whereas 5PTase genes in Group I underwent several duplication events in angiosperm, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for 5PTase duplications in angiosperm. Plant 5PTases have more members than that of animals, and most plant 5PTase genes appear to have evolved under strong purifying selection. The paralogs have diverged in substrate specificity and expression pattern, showing evidence of selection pressure. Meanwhile, the increase in 5PTases and divergences in sequence, expression, and substrate might have contributed to the divergent functions of 5PTase genes, allowing the angiosperms to successfully adapt to a great number of ecological niches.

Genome-Wide Identification and Characterization of the OPR Gene Family in Wheat (Triticum aestivum L.)

The 12-oxo-phytodienoic acid reductases (OPRs), which belong to the old yellow enzyme (OYE) family, are flavin mononucleotide (FMN)-dependent oxidoreductases with critical functions in plants. Despite the clear characteristics of growth and development, as well as the defense responses in Arabidopsis, tomato, rice, and maize, the potential roles of OPRs in wheat are not fully understood. Here, forty-eight putative OPR genes were found and classified into five subfamilies, with 6 in sub. I, 4 in sub. II, 33 in sub. III, 3 in sub. IV, and 2 in sub. V. Similar gene structures and conserved protein motifs of TaOPRs in wheat were identified in the same subfamilies. An analysis of cis-acting elements in promoters revealed that the functions of OPRs in wheat were mostly related to growth, development, hormones, biotic, and abiotic stresses. A total of 14 wheat OPR genes were identified as tandem duplicated genes, while 37 OPR genes were segmentally duplicated genes. The expression patterns of TaOPRs were tissue- and stress-specific, and the expression of TaOPRs could be regulated or induced by phytohormones and various stresses. Therefore, there were multiple wheat OPR genes, classified into five subfamilies, with functional diversification and specific expression patterns, and to our knowledge, this was the first study to systematically investigate the wheat OPR gene family. The findings not only provide a scientific foundation for the comprehensive understanding of the wheat OPR gene family, but could also be helpful for screening more candidate genes and breeding new varieties of wheat, with a high yield and stress resistance.

Metabolomic and transcriptomic changes underlying cold and anaerobic stresses after storage of table grapes

The currently accepted paradigm is that fruits and vegetables should be consumed fresh and that their quality deteriorates during storage; however, there are indications that some metabolic properties can, in fact, be improved. We examined the effects of low temperature and high-CO₂ conditions on table grapes, Vitis vinifera L. cv. 'Superior Seedless'. Berries were sampled at harvest (T0) and after low-temperature storage for 6 weeks under either normal atmosphere conditions (TC) or under an O₂ level of 5 kPa and elevated CO₂ levels of 5, 10 or 15 kPa (T5, T10, T15). Accumulation of 10 stilbenes, including E-ε-viniferin, E-miyabenol C and piceatannol, significantly increased under TC treatment as compared to T0 or T15. Sensory analysis demonstrated that elevated CO₂ elicited dose-dependent off-flavor accumulation. These changes were accompanied by an accumulation of 12 volatile metabolites, e.g., ethyl acetate and diacetyl, that imparted disagreeable flavors to fresh fruit. Transcriptome analysis revealed enrichment of genes involved in pyruvate metabolism and the phenylpropanoid pathway. One of the transcription factors induced at low temperature but not under high CO₂ was VvMYB14, which regulates stilbene biosynthesis. Our findings reveal the potential to alter the levels of targeted metabolites in stored produce through understanding the effects of postharvest treatments.

Transcriptome sequencing of Antarctic moss under salt stress emphasizes the important roles of the ROS-scavenging system

Mosses are predominant terrestrial vegetation in Antarctica. Their distributions appear to be controlled more by water and salinity. The Antarctic moss Pohlia nutans can tolerate high salt stress. Here, high-throughput sequencing was employed to study the transcriptional characteristics of P. nutans under salt stress. Differentially expressed genes (DEGs) analysis showed that 1340 genes were significantly upregulated and 831 genes were markedly downregulated. The expression of representative DEGs including abscisic acid (ABA) and Jasmonates (JAs) pathway-related genes, antioxidant enzyme genes, and flavonoid biosynthesis-related genes were analyzed by real-time PCR and most were upregulated after salt stress. Furthermore, malondialdehyde (MDA) content was significantly increased after salt treatment. The levels of hydroxyl free radical (∙OH) first rose then quickly decreased. In addition, the activities of antioxidant enzymes, such as catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD), and the flavonoid content were enhanced after salt stress. Exogenous application of ABA, Methyl jasmonate (MeJA) or proanthocyanidins (PA) improved the performance of P. nutans in response to high salt stress. Furthermore, real-time PCR showed that ABA or MeJA treatment upregulated the gene expression of antioxidant and flavonoid biosynthesis-related enzymes. These results suggest that the responses of P. nutans under salt stress are involved in activating phytohormone signaling pathways which trigger two main antioxidant defense systems (i.e., antioxidant enzymes and flavonoids) for protecting cell and scavenging reactive oxygen species.

RNA-Seq in the discovery of a sparsely expressed scent-determining monoterpene synthase in lavender (Lavandula)

Using RNA-Seq, we cloned and characterized a unique monoterpene synthase responsible for the formation of a scent-determining S-linalool constituent of lavender oils from Lavandula × intermedia. Several species of Lavandula produce essential oils (EOs) consisting mainly of monoterpenes including linalool, one of the most abundant and scent-determining oil constituents. Although R-linalool dominates the EOs of lavenders, varying amounts (depending on the species) of the S-linalool enantiomer can also be found in these plants. Despite its relatively low abundance, S-linalool contributes a sweet, pleasant scent and is an important constituent of lavender EOs. While several terpene synthase genes including R-linalool synthase have been cloned from lavenders many important terpene synthases including S-linalool synthase have not been described from these plants. In this study, we employed RNA-Seq and other complementary sequencing data to clone and functionally characterize the sparsely expressed S-linalool synthase cDNA (LiS-LINS) from Lavandula × intermedia. Recombinant LiS-LINS catalyzed the conversion of the universal monoterpene precursor geranyl diphosphate to S-linalool as the sole product. Intriguingly, LiS-LINS exhibited very low (~ 30%) sequence similarity to other Lavandula terpene synthases, including R-linalool synthase. However, the predicted 3D structure of this protein, including the composition and arrangement of amino acids at the active site, is highly homologous to known terpene synthase proteins. LiS-LINS transcripts were detected in flowers, but were much less abundant than those corresponding to LiR-LINS, paralleling enantiomeric composition of linalool in L. × intermedia oils. These data indicate that production of S-linalool is at least partially controlled at the level of transcription from LiS-LINS. The cloned LiS-LINS cDNA may be used to enhance oil composition in lavenders and other plants through metabolic engineering.

Proteomic analysis dissects the impact of nodulation and biological nitrogen fixation on Vicia faba root nodule physiology

Symbiotic nitrogen fixation in root nodules of legumes is a highly important biological process which is only poorly understood. Root nodule metabolism differs from that of roots. Differences in root and nodule metabolism are expressed by altered protein abundances and amenable to quantitative proteome analyses. Differences in the proteomes may either be tissue specific and related to the presence of temporary endosymbionts (the bacteroids) or related to nitrogen fixation activity. An experimental setup including WT bacterial strains and strains not able to conduct symbiotic nitrogen fixation as well as root controls enables identification of tissue and nitrogen fixation specific proteins. Root nodules are specialized plant organs housing and regulating the mutual symbiosis of legumes with nitrogen fixing rhizobia. As such, these organs fulfill unique functions in plant metabolism. Identifying the proteins required for the metabolic reactions of nitrogen fixation and those merely involved in sustaining the rhizobia:plant symbiosis, is a challenging task and requires an experimental setup which allows to differentiate between these two physiological processes. Here, quantitative proteome analyses of nitrogen fixing and non-nitrogen fixing nodules as well as fertilized and non-fertilized roots were performed using Vicia faba and Rhizobium leguminosarum. Pairwise comparisons revealed altered enzyme abundance between active and inactive nodules. Similarly, general differences between nodules and root tissue were observed. Together, these results allow distinguishing the proteins directly involved in nitrogen fixation from those related to nodulation. Further observations relate to the control of nodulation by hormones and provide supportive evidence for the previously reported correlation of nitrogen and sulfur fixation in these plant organs. Additionally, data on altered protein abundance relating to alanine metabolism imply that this amino acid may be exported from the symbiosomes of V. faba root nodules in addition to ammonia. Data are available via ProteomeXchange with identifier PXD008548.

Comparative in Silico Analysis of Ferric Reduction Oxidase (FRO) Genes Expression Patterns in Response to Abiotic Stresses, Metal and Hormone Applications

The ferric reduction oxidase (FRO) gene family is involved in various biological processes widely found in plants and may play an essential role in metal homeostasis, tolerance and intricate signaling networks in response to a number of abiotic stresses. Our study describes the identification, characterization and evolutionary relationships of FRO genes families. Here, total 50 FRO genes in Plantae and 15 ‘FRO like’ genes in non-Plantae were retrieved from 16 different species. The entire FRO genes have been divided into seven clades according to close similarity in biological and functional behavior. Three conserved domains were common in FRO genes while in two FROs sub genome have an extra NADPH-Ox domain, separating the function of plant FROs. OsFRO1 and OsFRO7 genes were expressed constitutively in rice plant. Real-time RT-PCR analysis demonstrated that the expression of OsFRO1 was high in flag leaf, and OsFRO7 gene expression was maximum in leaf blade and flag leaf. Both genes showed vigorous expressions level in response to different abiotic and hormones treatments. Moreover, the expression of both genes was also substantial under heavy metal stresses. OsFRO1 gene expression was triggered following 6 h under Zn, Pb, Co and Ni treatments, whereas OsFRO7 gene expression under Fe, Pb and Ni after 12 h, Zn and Cr after 6 h, and Mn and Co after 3 h treatments. These findings suggest the possible involvement of both the genes under abiotic and metal stress and the regulation of phytohormones. Therefore, our current work may provide the foundation for further functional characterization of rice FRO genes family.

Ligand-receptor co-evolution shaped the jasmonate pathway in land plants

The phytohormone jasmonoyl-isoleucine (JA-Ile) regulates defense, growth and developmental responses in vascular plants. Bryophytes have conserved sequences for all JA-Ile signaling pathway components but lack JA-Ile. We show that, in spite of 450 million years of independent evolution, the JA-Ile receptor COI1 is functionally conserved between the bryophyte Marchantia polymorpha and the eudicot Arabidopsis thaliana but COI1 responds to different ligands in each species. We identified the ligand of Marchantia MpCOI1 as two isomeric forms of the JA-Ile precursor dinor-OPDA (dinor-cis-OPDA and dinor-iso-OPDA). We demonstrate that AtCOI1 functionally complements Mpcoi1 mutation and confers JA-Ile responsiveness and that a single-residue substitution in MpCOI1 is responsible for the evolutionary switch in ligand specificity. Our results identify the ancestral bioactive jasmonate and clarify its biosynthetic pathway, demonstrate the functional conservation of its signaling pathway, and show that JA-Ile and COI1 emergence in vascular plants required co-evolution of hormone biosynthetic complexity and receptor specificity.

OPEN GLUME1: a key enzyme reducing the precursor of JA, participates in carbohydrate transport of lodicules during anthesis in rice

OG1 is involved in JA-regulated anthesis by modulating carbohydrate transport of lodicules in rice. Flowering plants have evolved a sophisticated regulatory network to coordinate anthesis and maximize reproductive success. In addition to various environmental conditions, the plant hormone jasmonic acid and its derivatives (JAs) are involved in anthesis. However, the underlying mechanism remains largely unexplored. Here, we report a JA-defective mutant in rice (Oryza sativa), namely open glume 1, which has dysfunctional lodicules that lead to open glumes following anthesis. Map-based cloning and subsequent complementation tests confirmed that OG1 encodes a peroxisome-localized 12-oxo-phytodienoic acid reductase-a key enzyme that reduces the precursor of JA. Loss-of-function of OG1 resulted in almost no JA accumulation. Exogenous JA treatment completely rescued the defects caused by the og1 mutation. Further studies revealed that intracellular metabolism was disrupted in the lodicules of og1 mutant. At the mature plant stage, most seeds of the mutant were malformed with significantly reduced starch content. We speculate that JA or JA signaling mediates the carbohydrate transport of lodicules during anthesis, and signal the onset of cell degradation in lodicules after anthesis. We conclude that the OPEN GLUME 1 gene that produces a key enzyme involved in reducing the precursor of JA in JA biosynthesis and is involved in carbohydrate transport underlying normal lodicule function during anthesis in rice.

An OPR3-independent pathway uses 4,5-didehydrojasmonate for jasmonate synthesis

Biosynthesis of the phytohormone jasmonoyl-isoleucine (JA-Ile) requires reduction of the JA precursor 12-oxo-phytodienoic acid (OPDA) by OPDA reductase 3 (OPR3). Previous analyses of the opr3-1 Arabidopsis mutant suggested an OPDA signaling role independent of JA-Ile and its receptor COI1; however, this hypothesis has been challenged because opr3-1 is a conditional allele not completely impaired in JA-Ile biosynthesis. To clarify the role of OPR3 and OPDA in JA-independent defenses, we isolated and characterized a loss-of-function opr3-3 allele. Strikingly, opr3-3 plants remained resistant to necrotrophic pathogens and insect feeding, and activated COI1-dependent JA-mediated gene expression. Analysis of OPDA derivatives identified 4,5-didehydro-JA in wounded wild-type and opr3-3 plants. OPR2 was found to reduce 4,5-didehydro-JA to JA, explaining the accumulation of JA-Ile and activation of JA-Ile-responses in opr3-3 mutants. Our results demonstrate that in the absence of OPR3, OPDA enters the β-oxidation pathway to produce 4,5-ddh-JA as a direct precursor of JA and JA-Ile, thus identifying an OPR3-independent pathway for JA biosynthesis.

Classification of Protein Structure Classes on Flexible Neutral Tree

Accurate classification on protein structural is playing an important role in Bioinformatics. An increase in evidence demonstrates that a variety of classification methods have been employed in such a field. In this research, the features of amino acids composition, secondary structure's feature, and correlation coefficient of amino acid dimers and amino acid triplets have been used. Flexible neutral tree (FNT), a particular tree structure neutral network, has been employed as the classification model in the protein structures' classification framework. Considering different feature groups owing diverse roles in the model, impact factors of different groups have been put forward in this research. In order to evaluate different impact factors, Impact Factors Scaling (IFS) algorithm, which aim at reducing redundant information of the selected features in some degree, have been put forward. To examine the performance of such framework, the 640, 1189, and ASTRAL datasets are employed as the low-homology protein structure benchmark datasets. Experimental results demonstrate that the performance of the proposed method is better than the other methods in the low-homology protein tertiary structures.

Whole-Genome Identification and Expression Pattern of the Vicinal Oxygen Chelate Family in Rapeseed (Brassica napus L.)

Vicinal oxygen chelate proteins (VOC) are members of the metalloenzyme superfamily, which plays roles in many biological reactions. Some members of the VOC superfamily have been systematically characterized but not in Brassica napus. In this study, 38 VOC genes were identified based on their conserved domains. The present results revealed that most of the BnaVOC genes have few introns, and all contained the typical VOC structure of βαβββ modules. The BnaVOC genes are distributed unevenly across 15 chromosomes in B. napus and occur as gene clusters on chromosomes C5 and A6. The synteny and phylogenetic analyses revealed that the VOC gene family is a consequence of mesopolyploidy events that occurred in Brassica evolution, and whole-genome duplication and segmental duplication played a major role in the expansion of the BnaVOC gene family. The expression profile analysis indicated that the expression of most BnaVOCs was increased in the leaves and late stage seeds. Further results indicated that seeds of B. napus with a high oil content show higher expression levels under drought stress conditions, suggesting that BnaVOCs not only respond to abiotic stress but may also affect lipid metabolism in drought stress. This present study provides a comprehensive overview of the VOC gene family and provides new insights into their biological function in B. napus evolution.

Identification of Jasmonic Acid and Jasmonoyl-Isoleucine, and Characterization of AOS, AOC, OPR and JAR1 in the Model Lycophyte Selaginella moellendorffii

Jasmonic acid (JA) is involved in a variety of physiological responses in seed plants. However, the detection and role of JA in lycophytes, a group of seedless vascular plants, have remained elusive until recently. This study provides the first evidence of 12-oxo-phytodienoic acid (OPDA), JA and jasmonoyl-isoleucine (JA-Ile) in the model lycophyte Selaginella moellendorffii. Mechanical wounding stimulated the accumulation of OPDA, JA and JA-Ile. These data were corroborated by the detection of enzymatically active allene oxide synthase (AOS), allene oxide cyclase (AOC), 12-oxo-phytodienoic acid reductase 3 (OPR3) and JA-Ile synthase (JAR1) in S. moellendorffii. SmAOS2 is involved in the first committed step of JA biosynthesis. SmAOC1 is a crucial enzyme for generating the basic structure of jasmonates and is actively involved in the formation of OPDA. SmOPR5, a functionally active OPR3-like enzyme, is also vital for the reduction of (+)-cis-OPDA, the only isomer of the JA precursor. The conjugation of JA to Ile by SmJAR1 demonstrates that S. moellendorffii produces JA-Ile. Thus, the four active enzymes have characteristics similar to those in seed plants. Wounding and JA treatment induced the expression of SmAOC1 and SmOPR5. Furthermore, JA inhibited the growth of shoots in S. moellendorffii, which suggests that JA functions as a signaling molecule in S. moellendorffii. This study proposes that JA evolved as a plant hormone for stress adaptation, beginning with the emergence of vascular plants.

A putative 12-oxophytodienoate reductase gene CsOPR3 from Camellia sinensis, is involved in wound and herbivore infestation responses

12-Oxophytodienoate reductase (OPR) is a key enzyme in the biosynthesis of jasmonic acid (JA), which plays an important role in plant defense responses. Although multiple isoforms of OPRs have been identified in various annual herbaceous plants, genes encoding these enzymes in perennial woody plants have yet to be fully investigated. In the tea plant, Camellia sinensis (L.), no OPR genes have been isolated, and their possible roles in tea plant development and defense mechanism remain unknown. In this study, a putative OPR gene, designated as CsOPR3, was isolated from tea plants for the first time through the rapid amplification of cDNA ends. The open reading frame of CsOPR3 is 1197bp in length, and encodes a protein of 398 amino acids. Real-time qPCR analysis revealed that CsOPR3 was expressed in different organs. In particular, CsOPR3 was highly expressed in flowers, leaves and stems but was weakly expressed in roots and seeds. CsOPR3 expression could be rapidly induced by mechanical wounding, and increased JA levels were correlated with the wound-induced CsOPR3 expression. The infestation of the tea geometrid (TG) Ectropis obliqua Prout, regurgitant derived from TG and exogenous JA application could enhance the CsOPR3 expression. Our study is the first to report that CsOPR3 plays an important role in JA biosynthesis and tea plant defense against herbivorous insects.

Evolution of jasmonate biosynthesis and signaling mechanisms

Jasmonates are phytohormones that modulate a wide spectrum of plant physiological processes, especially defense against herbivores and necrotrophs. The molecular mechanisms of jasmonate biosynthesis and signaling have been well characterized in model plants. In this review, we provide an in-depth analysis and overview of the origin and evolution of the jasmonate biosynthesis and signaling pathways. Furthermore, we discuss the striking parallels between jasmonate and auxin signaling mechanisms, which reveals a common ancestry of these signaling mechanisms. Finally, we highlight the importance of studying jasmonate biosynthesis and signaling in lower plants.

Plant Genome Duplication Database

Genome duplication, widespread in flowering plants, is a driving force in evolution. Genome alignments between/within genomes facilitate identification of homologous regions and individual genes to investigate evolutionary consequences of genome duplication. PGDD (the Plant Genome Duplication Database), a public web service database, provides intra- or interplant genome alignment information. At present, PGDD contains information for 47 plants whose genome sequences have been released. Here, we describe methods for identification and estimation of dates of genome duplication and speciation by functions of PGDD.The database is freely available at http://chibba.agtec.uga.edu/duplication/.

Characterization of 12-Oxophytodienoic Acid Reductases from Rose-scented Geranium (Pelargonium graveolens)

Pelargonium graveolens L'Hér, also referred to as rose geranium, is a popular herbal plant with typical rosy fragrance largely based on the blend of monoterpenoid constituents. Among them, citronellol, which is biosynthesized from geraniol via double bond reduction, is the most abundant scent compound. In this study, three 12-oxophytodienoic acid reductases (PgOPR1–3) have been cloned from P. graveolens, as possible candidates for the double-bond reductase involved in citronellol biosynthesis. The bacterially expressed recombinant PgOPRs did not reduce geraniol to citronellol, but stereoselectively converted citral into ( S)-citronellal in the presence of NADPH. Thus, the α,β-unsaturated carbonyl moiety in the substrate is essential for the catalytic activity of PgOPRs, as reported for OPRs from other plants and structurally related yeast old yellow enzymes. PgOPRs promiscuously accepted linear and cyclic α,β-unsaturated carbonyl substrates, including methacrolein, a typical reactive carbonyl compound. The possible biotechnological applications for PgOPRs in plant metabolic engineering, based on their catalytic properties, are discussed herein.

Learning the Languages of the Chloroplast: Retrograde Signaling and Beyond

The chloroplast can act as an environmental sensor, communicating with the cell during biogenesis and operation to change the expression of thousands of proteins. This process, termed retrograde signaling, regulates expression in response to developmental cues and stresses that affect photosynthesis and yield. Recent advances have identified many signals and pathways-including carotenoid derivatives, isoprenes, phosphoadenosines, tetrapyrroles, and heme, together with reactive oxygen species and proteins-that build a communication network to regulate gene expression, RNA turnover, and splicing. However, retrograde signaling pathways have been viewed largely as a means of bilateral communication between organelles and nuclei, ignoring their potential to interact with hormone signaling and the cell as a whole to regulate plant form and function. Here, we discuss new findings on the processes by which organelle communication is initiated, transmitted, and perceived, not only to regulate chloroplastic processes but also to intersect with cellular signaling and alter physiological responses.

Comprehensive Genomic Analysis and Expression Profiling of the NOX Gene Families under Abiotic Stresses and Hormones in Plants

Plasma membrane NADPH oxidases (NOXs) are key producers of reactive oxygen species under both normal and stress conditions in plants and they form functional subfamilies. Studies of these subfamilies indicated that they show considerable evolutionary selection. We performed a comparative genomic analysis that identified 50 ferric reduction oxidases (FRO) and 77 NOX gene homologs from 20 species representing the eight major plant lineages within the supergroup Plantae: glaucophytes, rhodophytes, chlorophytes, bryophytes, lycophytes, gymnosperms, monocots, and eudicots. Phylogenetic and structural analysis classified these FRO and NOX genes into four well-conserved groups represented as NOX, FRO I, FRO II, and FRO III. Further analysis of NOXs of phylogenetic and exon/intron structures showed that single intron loss and gain had occurred, yielding the diversified gene structures during the evolution of NOXs family genes and which were classified into four conserved subfamilies which are represented as Sub.I, Sub.II, Sub.III, and Sub.IV. Additionally, both available global microarray data analysis and quantitative real-time PCR experiments revealed that the NOX genes in Arabidopsis and rice (Oryza sativa) have different expression patterns in different developmental stages, various abiotic stresses and hormone treatments. Finally, coexpression network analysis of NOX genes in Arabidopsis and rice revealed that NOXs have significantly correlated expression profiles with genes which are involved in plants metabolic and resistance progresses. All these results suggest that NOX family underscores the functional diversity and divergence in plants. This finding will facilitate further studies of the NOX family and provide valuable information for functional validation of this family in plants.

TaOPR2 encodes a 12-oxo-phytodienoic acid reductase involved in the biosynthesis of jasmonic acid in wheat (Triticum aestivum L.)

The 12-oxo-phytodienoic acid reductases (OPRs) are involved in the various processes of growth and development in plants, and classified into the OPRⅠ and OPRⅡ subgroups. In higher plants, only OPRⅡ subgroup genes take part in the biosynthesis of endogenous jasmonic acid. In this study, we isolated a novel OPRⅡ subgroup gene named TaOPR2 (GeneBank accession: KM216389) from the thermo-sensitive genic male sterile (TGMS) wheat cultivar BS366. TaOPR2 was predicted to encode a protein with 390 amino acids. The encoded protein contained the typical oxidored_FMN domain, the C-terminus peroxisomal-targeting signal peptide, and conserved FMN-binding sites. TaOPR2 was mapped to wheat chromosome 7B and located on peroxisome. Protein evolution analysis revealed that TaOPR2 belongs to the OPRⅡ subgroup and shares a high degree of identity with other higher plant OPR proteins. The quantitative real-time PCR results indicated that the expression of TaOPR2 is inhibited by abscisic acid (ABA), salicylic acid (SA), gibberellic acid (GA3), low temperatures and high salinity. In contrast, the expression of TaOPR2 can be induced by wounding, drought and methyl jasmonate (MeJA). Furthermore, the transcription level of TaOPR2 increased after infection with Puccinia striiformis f. sp. tritici and Puccinia recondite f. sp. tritici. TaOPR2 has NADPH-dependent oxidoreductase activity. In addition, the constitutive expression of TaOPR2 can rescue the male sterility phenotype of Arabidopsis mutant opr3. These results suggest that TaOPR2 is involved in the biosynthesis of jasmonic acid (JA) in wheat.

Identification and characterization of TIFY family genes in Brachypodium distachyon

The TIFY family is a plant-specific gene family encoding proteins characterized by a conserved TIFY domain. This family encodes four subfamilies of proteins, including ZIM-like (ZML), TIFY, PPD and JASMONATE ZIM-Domain (JAZ) proteins. TIFY proteins play important roles in plant development and stress responses. In this study, 21 BdTIFYs were identified in Brachypodium distachyon through genome-wide analysis, including 15 JAZ and 6 ZML genes. Analysis of the distribution of conserved domains showed that there are three additional domains (CCT domain, GATA domain and Jas domain) in the BdTIFY proteins besides the TIFY domain. Phylogenetic analysis indicated that these 21 proteins were classified into two major groups. Expression profile of BdTIFY genes in response to abiotic stresses and phytohormones was analyzed using quantitative real-time RT-PCR. Among 21 BdTIFY genes, 12 of them were induced by JA treatment, and 4 of them were induced by ABA treatment. Most of BdTIFY genes were responsive to one or more abiotic stresses including drought, salinity, low temperature and heat. Especially, BdTIFY5, 9a, 9b, 10c and 11a were significantly up-regulated by multiple abiotic stresses. These results provided important clues for functional analysis of TIFY family genes in B. distachyon.