Large-scale plant protein subcellular location prediction

Genome-Wide Identification of the Alfin-like Gene Family in Cotton (Gossypium hirsutum) and the GhAL19 Gene Negatively Regulated Drought and Salt Tolerance

Alfin-like (AL) is a small plant-specific gene family characterized by a PHD-finger-like structural domain at the C-terminus and a DUF3594 structural domain at the N-terminus, and these genes play prominent roles in plant development and abiotic stress response. In this study, we conducted genome-wide identification and analyzed the AL protein family in Gossypium hirsutum cv. NDM8 to assess their response to various abiotic stresses for the first time. A total of 26 AL genes were identified in NDM8 and classified into four groups based on a phylogenetic tree. Moreover, cis-acting element analysis revealed that multiple phytohormone response and abiotic stress response elements were highly prevalent in AL gene promoters. Further, we discovered that the GhAL19 gene could negatively regulate drought and salt stresses via physiological and biochemical changes, gene expression, and the VIGS assay. The study found there was a significant increase in POD and SOD activity, as well as a significant change in MDA in VIGS-NaCl and VIGS-PEG plants. Transcriptome analysis demonstrated that the expression levels of the ABA biosynthesis gene (GhNCED1), signaling genes (GhABI1, GhABI2, and GhABI5), responsive genes (GhCOR47, GhRD22, and GhERFs), and the stress-related marker gene GhLEA14 were regulated in VIGS lines under drought and NaCl treatment. In summary, GhAL19 as an AL TF may negatively regulate tolerance to drought and salt by regulating the antioxidant capacity and ABA-mediated pathway.

Genomic survey and expression analysis of LcARFs reveal multiple functions to somatic embryogenesis in Liriodendron

BACKGROUND

Auxin response factors (ARFs) are critical transcription factors that mediate the auxin signaling pathway and are essential for regulating plant growth. However, there is a lack of understanding regarding the ARF gene family in Liriodendron chinense, a vital species in landscaping and economics. Thus, further research is needed to explore the roles of ARFs in L. chinense and their potential applications in plant development.

RESULT

In this study, we have identified 20 LcARF genes that belong to three subfamilies in the genome of L. chinense. The analysis of their conserved domains, gene structure, and phylogeny suggests that LcARFs may be evolutionarily conserved and functionally similar to other plant ARFs. The expression of LcARFs varies in different tissues. Additionally, they are also involved in different developmental stages of somatic embryogenesis. Overexpression of LcARF1, LcARF2a, and LcARF5 led to increased activity within callus. Additionally, our promoter-GFP fusion study indicated that LcARF1 may play a role in embryogenesis. Overall, this study provides insights into the functions of LcARFs in plant development and embryogenesis, which could facilitate the improvement of somatic embryogenesis in L. chinense.

CONCLUSION

The research findings presented in this study shed light on the regulatory roles of LcARFs in somatic embryogenesis in L. chinense and may aid in accelerating the breeding process of this tree species. By identifying the specific LcARFs involved in different stages of somatic embryogenesis, this study provides a basis for developing targeted breeding strategies aimed at optimizing somatic embryogenesis in L. chinense, which holds great potential for improving the growth and productivity of this economically important species.

Structure and Expression Analysis of PtrSUS, PtrINV, PtrHXK, PtrPGM, and PtrUGP Gene Families in Populus trichocarpa Torr. and Gray

Exogenous nitrogen and carbon can affect plant cell walls, which are composed of structural carbon. Sucrose synthase (SUS), invertase (INV), hexokinase (HXK), phosphoglucomutase (PGM), and UDP-glucose pyrophosphorylase (UGP) are the key enzymes of sucrose metabolism involved in cell wall synthesis. To understand whether these genes are regulated by carbon and nitrogen to participate in structural carbon biosynthesis, we performed genome-wide identification, analyzed their expression patterns under different carbon and nitrogen treatments, and conducted preliminary functional verification. Different concentrations of nitrogen and carbon were applied to poplar (Populus trichocarpa Torr. and Gray), which caused changes in cellulose, lignin, and hemicellulose contents. In poplar, 6 SUSs, 20 INVs, 6 HXKs, 4 PGMs, and 2 UGPs were identified. Moreover, the physicochemical properties, collinearity, and tissue specificity were analyzed. The correlation analysis showed that the expression levels of PtrSUS3/5, PtrNINV1/2/3/5/12, PtrCWINV3, PtrVINV2, PtrHXK5/6, PtrPGM1/2, and PtrUGP1 were positively correlated with the cellulose content. Meanwhile, the knockout of PtrNINV12 significantly reduced the cellulose content. This study could lay the foundation for revealing the functions of SUSs, INVs, HXKs, PGMs, and UGPs, which affected structural carbon synthesis regulated by nitrogen and carbon, proving that PtrNINV12 is involved in cell wall synthesis.

Identification of WUSCHEL-related homeobox gene and truncated small peptides in transformation efficiency improvement in Eucalyptus

BACKGROUND

The WUSCHEL-related Homeobox (WOX) genes, which encode plant-specific homeobox (HB) transcription factors, play crucial roles in regulating plant growth and development. However, the functions of WOX genes are little known in Eucalyptus, one of the fastest-growing tree resources with considerable widespread cultivation worldwide.

RESULTS

A total of nine WOX genes named EgWOX1-EgWOX9 were retrieved and designated from Eucalyptus grandis. From the three divided clades marked as Modern/WUS, Intermediate and Ancient, the largest group Modern/WUS (6 EgWOXs) contains a specific domain with 8 amino acids: TLQLFPLR. The collinearity, cis-regulatory elements, protein-protein interaction network and gene expression analysis reveal that the WUS proteins in E. grandis involve in regulating meristems development and regeneration. Furthermore, by externally adding of truncated peptides isolated from WUS specific domain, the transformation efficiency in E. urophylla × E. grandis DH32-29 was significant enhanced. The transcriptomics data further reveals that the use of small peptides activates metabolism pathways such as starch and sucrose metabolism, phenylpropanoid biosynthesis and flavonoid biosynthesis.

CONCLUSIONS

Peptides isolated from WUS protein can be utilized to enhance the transformation efficiency in Eucalyptus, thereby contributing to the high-efficiency breeding of Eucalyptus.

Functional Study of Amorpha fruticosa WRKY20 Gene in Response to Drought Stress

The WRKY gene family in plants regulates the plant's response to drought through regulatory networks and hormone signaling. AfWRKY20 (MT859405) was cloned from Amorpha fruticosa (A. fruticosa) seedlings using RT-PCR. The binding properties of the AfWRKY20 protein and the W-box (a DNA cis-acting element) were verified both in vivo and in vitro using EMSA and Dual-Luciferase activity assays. RT-qPCR detected that the total expression level of AfWRKY20 in leaves and roots was 22 times higher in the 30% PEG6000 simulated drought treatment compared to the untreated group. Under the simulated drought stress treatments of sorbitol and abscisic acid (ABA), the transgenic tobacco with the AfWRKY20 gene showed enhanced drought resistance at the germination stage, with significantly increased germination rate, green leaf rate, fresh weight, and root length compared to the wild-type (WT) tobacco. In addition, the superoxide dismutase (SOD) activity, chlorophyll content, and Fv/Fm ratio of AfWRKY20 transgenic tobacco were significantly higher than those of the WT tobacco under natural drought stress, while the malondialdehyde (MDA) content and 3,3'-diaminobenzidine (DAB) and nitroblue tetrazolium (NBT) staining levels were lower. The expression levels of oxidation kinase genes (NbSOD, NbPOD, and NbCAT) in transgenic tobacco under drought stress were significantly higher than those in WT tobacco. This enhancement in gene expression improved the ability of transgenic tobacco to detoxify reactive oxygen species (ROS). The survival rate of transgenic tobacco after natural drought rehydration was four times higher than that of WT tobacco. In summary, this study revealed the regulatory mechanism of AfWRKY20 in response to drought stress-induced ABA signaling, particularly in relation to ROS. This finding provides a theoretical basis for understanding the pathways of WRKY20 involved in drought stress, and offers genetic resources for molecular plant breeding aimed at enhancing drought resistance.

The Moonlighting Function of Soybean Disordered Methyl-CpG-Binding Domain 10c Protein

Intrinsically disordered proteins (IDPs) are multifunctional due to their ability to adopt different structures depending on the local conditions. The intrinsically disordered regions of methyl-CpG-binding domain (MBD) proteins play important roles in regulating growth and development by interpreting DNA methylation patterns. However, whether MBDs have a stress-protective function is far from clear. In this paper, soybean GmMBD10c protein, which contains an MBD and is conserved in Leguminosae, was predicted to be located in the nucleus. It was found to be partially disordered by bioinformatic prediction, circular dichroism and a nuclear magnetic resonance spectral analysis. The enzyme activity assay and SDS-PAGE results showed that GmMBD10c can protect lactate dehydrogenase and a broad range of other proteins from misfolding and aggregation induced by the freeze-thaw process and heat stress, respectively. Furthermore, overexpression of GmMBD10c enhanced the salt tolerance of Escherichia coli. These data validate that GmMBD10c is a moonlighting protein with multiple functions.

Computational models for prediction of protein-protein interaction in rice and Magnaporthe grisea

INTRODUCTION

Plant-microbe interactions play a vital role in the development of strategies to manage pathogen-induced destructive diseases that cause enormous crop losses every year. Rice blast is one of the severe diseases to rice Oryza sativa (O. sativa) due to Magnaporthe grisea (M. grisea) fungus. Protein-protein interaction (PPI) between rice and fungus plays a key role in causing rice blast disease.

METHODS

In this paper, four genomic information-based models such as (i) the interolog, (ii) the domain, (iii) the gene ontology, and (iv) the phylogenetic-based model are developed for predicting the interaction between O. sativa and M. grisea in a whole-genome scale.

RESULTS AND DISCUSSION

A total of 59,430 interacting pairs between 1,801 rice proteins and 135 blast fungus proteins are obtained from the four models. Furthermore, a machine learning model is developed to assess the predicted interactions. Using composition-based amino acid composition (AAC) and conjoint triad (CT) features, an accuracy of 88% and 89% is achieved, respectively. When tested on the experimental dataset, the CT feature provides the highest accuracy of 95%. Furthermore, the specificity of the model is verified with other pathogen-host datasets where less accuracy is obtained, which confirmed that the model is specific to O. sativa and M. grisea. Understanding the molecular processes behind rice resistance to blast fungus begins with the identification of PPIs, and these predicted PPIs will be useful for drug design in the plant science community.

NF-YB family transcription factors in Arabidopsis: Structure, phylogeny, and expression analysis in biotic and abiotic stresses

Nuclear factor-Y (NF-Y) transcription factors (TFs) are conserved heterotrimeric complexes present and widespread across eukaryotes. Three main subunits make up the structural and functional aspect of the NF-Y TFs: NF-YA, NF-YB and NF-YC, which bind to the conserved CCAAT- box of the promoter region of specific genes, while also interacting with each other, thereby forming myriad combinations. The NF-YBs are expressed differentially in various tissues and plant development stages, likely impacting many of the cellular processes constitutively and under stress conditions. In this study, ten members of NF-YB family from Arabidopsis thaliana were identified and expression profiles were mined from microarray data under different biotic and abiotic conditions, revealing key insights into the involvement of this class of proteins in the cellular and biological processes in Arabidopsis. Analysis of cis-acting regulatory elements (CAREs) indicated the presence of abiotic and biotic stress-related transcription factor binding sites (TFBs), shedding light on the multifaceted roles of these TFs. Microarray data analysis inferred distinct patterns of expression in various tissues under differing treatments such as drought, cold and heat stress as well as bacterial, fungal, and viral stress, indicating their likelihood of having an expansive range of regulatory functions under native and stressed conditions; while quantitative real-time PCR (qRT-PCR) based expression analysis revealed that these TFs get real-time-modulated in a stress dependent manner. This study, overall, provides an understanding of the AtNF-YB family of TFs in their regulation and participation in various morphogenetic and defense- related pathways and can provide insights for development of transgenic plants for trait dependent studies.

Genome-wide characterization of the SHORT INTER-NODES/STYLISH and Shi-Related Sequence family in Gossypium hirsutum and functional identification of GhSRS21 under salt stress

Saline stress is a significant factor that caused crop growth inhibition and yield decline. SHORT INTERNODES/STYLISH (SHI/STY) and SHI-RELATED SEQUENCE (SRS) transcription factors are specific to plants and share a conserved RING-like zinc-finger domain (CX₂CX₇CX₄CX₂C₂X₆C). However, the functions of SHI/STY and SRS genes in cotton responses to salt stress remain unclear. In this study, 26 GhSRSs were identified in Gossypium hirsutum, which further divided into three subgroups. Phylogenetic analysis of 88 SRSs from8 plant species revealed independent evolutionary pattern in some of SRSs derived from monocots. Conserved domain and subcellular location predication of GhSRSs suggested all of them only contained the conserved RING-like zinc-finger domain (DUF702) domain and belonged to nucleus-localized transcription factors except for the GhSRS22. Furthermore, synteny analysis showed structural variation on chromosomes during the process of cotton polyploidization. Subsequently, expression patterns of GhSRS family members in response to salt and drought stress were analyzed in G. hirsutum and identified a salt stress-inducible gene GhSRS21. The GhSRS21 was proved to localize in the nuclear and silencing it in G. hirsutum increased the cotton resistance to salt using the virus-induced gene silencing (VIGS) system. Finally, our transcriptomic data revealed that GhSRS21 negatively controlled cotton salt tolerance by regulating the balance between ROS production and scavenging. These results will increase our understanding of the SRS gene family in cotton and provide the candidate resistant gene for cotton breeding.

Phylogenomics-Based Reconstruction and Molecular Evolutionary Histories of Brassica Photoreceptor Gene Families

Photosensory proteins known as photoreceptors (PHRs) are crucial for delineating light environments in synchronization with other environmental cues and regulating their physiological variables in plants. However, this has not been well studied in the Brassica genus, which includes several important agricultural and horticultural crops. Herein, we identified five major PHR gene families—phytochrome (PHY), cryptochrome (CRY), phototropin (PHOT), F-box containing flavin binding proteins (ZTL/FKF1/LKP2), and UV RESISTANCE LOCUS 8 (UVR8)—genomic scales and classified them into subfamilies based on their phylogenetic clustering with Arabidopsis homologues. The molecular evolution characteristics of Brassica PHR members indicated indirect expansion and lost one to six gene copies at subfamily levels. The segmental duplication was possibly the driving force of the evolution and amplification of Brassica PHRs. Gene replication retention and gene loss events of CRY, PHY, and PHOT members found in diploid progenitors were highly conserved in their tetraploid hybrids. However, hybridization events were attributed to quantitative changes in UVR8 and ZTL/FKF1/LKP2 members. All PHR members underwent purifying selection. In addition, the transcript expression profiles of PHR genes in different tissue and in response to exogenous ABA, and abiotic stress conditions suggested their multiple biological significance. This study is helpful in understanding the molecular evolution characteristics of Brassica PHRs and lays the foundation for their functional characterization.

Genome-wide identification and expression analysis of the ftsH protein family and its response to abiotic stress in Nicotiana tabacum L

BACKGROUND

The filamentous temperature-sensitive H protease (ftsH) gene family plays an important role in plant growth and development. FtsH proteins belong to the AAA protease family. Studies have shown that it is a key gene for plant chloroplast development and photosynthesis regulation. In addition, the ftsH gene is also involved in plant response to stress. At present, the research and analysis of the ftsH gene family are conducted in microorganisms such as Escherichia coli and Oenococcus and various plants such as Arabidopsis, pear, rice, and corn. However, analysis reports on ftsH genes from tobacco (Nicotiana tabacum L.), an important model plant, are still lacking. Since ftsH genes regulate plant growth and development, it has become necessary to systematically study this gene in an economically important plant like tobacco.

RESULTS

This is the first study to analyze the ftsH gene from Nicotiana tabacum L. K326 (NtftsH). We identified 20 ftsH genes from the whole genome sequence, renamed them according to their chromosomal locations, and divided them into eight subfamilies. These 20 NtftsH genes were unevenly distributed across the 24 chromosomes. We found four pairs of fragment duplications. We further investigated the collinearity between these genes and related genes in five other species. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis identified differential expression patterns of NtftsH in different tissues and under various abiotic stress conditions.

CONCLUSIONS

This study provides a comprehensive analysis of the NtftsH gene family. The exon-intron structure and motif composition are highly similar in NtftsH genes that belong to the same evolutionary tree branch. Homology analysis and phylogenetic comparison of ftsH genes from several different plants provide valuable clues for studying the evolutionary characteristics of NtftsH genes. The NtftsH genes play important roles in plant growth and development, revealed by their expression levels in different tissues as well as under different stress conditions. Gene expression and phylogenetic analyses will provide the basis for the functional analysis of NtftsH genes. These results provide a valuable resource for a better understanding of the biological role of the ftsH genes in the tobacco plant.

The LONELY GUY gene family: from mosses to wheat, the key to the formation of active cytokinins in plants

LONELY GUY (LOG) was first identified in a screen of rice mutants with defects in meristem maintenance. In plants, LOG codes for cytokinin riboside 5'-monophosphate phosphoribohydrolase, which converts inactive cytokinin nucleotides directly to the active free bases. Many enzymes with the PGGxGTxxE motif have been misannotated as lysine decarboxylases; conversely not all enzymes containing this motif are cytokinin-specific LOGs. As LOG mutants clearly impact yield in rice, we investigated the LOG gene family in bread wheat. By interrogating the wheat (Triticum aestivum) genome database, we show that wheat has multiple LOGs. The close alignment of TaLOG1, TaLOG2 and TaLOG6 with the X-ray structures of two functional Arabidopsis thaliana LOGs allows us to infer that the wheat LOGs 1-11 are functional LOGs. Using RNA-seq data sets, we assessed TaLOG expression across 70 tissue types, their responses to various stressors, the pattern of cis-regulatory elements (CREs) and intron/exon patterns. TaLOG gene family members are expressed variously across tissue types. When the TaLOG CREs are compared with those of the cytokinin dehydrogenases (CKX) and glucosyltransferases (CGT), there is close alignment of CREs between TaLOGs and TaCKXs reflecting the key role of CKX in maintaining cytokinin homeostasis. However, we suggest that the main homeostatic mechanism controlling cytokinin levels in response to biotic and abiotic challenge resides in the CGTs, rather than LOG or CKX. However, LOG transgenics and identified mutants in rice variously impact yield, providing interesting avenues for investigation in wheat.

Metal Tolerance Protein Encoding Gene Family in Fagopyrum tartaricum: Genome-Wide Identification, Characterization and Expression under Multiple Metal Stresses

Metal tolerance proteins (MTP) as divalent cation transporters are essential for plant metal tolerance and homeostasis. However, the characterization and the definitive phylogeny of the MTP gene family in Fagopyrum tartaricum, and their roles in response to metal stress are still unknown. In the present study, MTP genes in Fagopyrum tartaricum were identified, and their phylogenetic relationships, structural characteristics, physicochemical parameters, as well as expression profiles under five metal stresses including Fe, Mn, Cu, Zn, and Cd were also investigated. Phylogenetic relationship analysis showed that 12 Fagopyrum tartaricum MTP genes were classified into three major clusters and seven groups. All FtMTPs had typical structural features of the MTP gene family and were predicted to be located in the cell vacuole. The upstream region of FtMTPs contained abundant cis-acting elements, implying their functions in development progress and stress response. Tissue-specific expression analysis results indicated the regulation of FtMTPs in the growth and development of Fagopyrum tataricum. Besides, the expression of most FtMTP genes could be induced by multiple metals and showed different expression patterns under at least two metal stresses. These findings provide useful information for the research of the metal tolerance mechanism and genetic improvement of Fagopyrum tataricum.

Genome-Wide Identification, Evolution and Expression Analysis of the Glutathione S-Transferase Supergene Family in Euphorbiaceae

Euphorbiaceae, a family of plants mainly grown in the tropics and subtropics, is also widely distributed all over the world and is well known for being rich in rubber, oil, medicinal materials, starch, wood and other economically important plant products. Glutathione S-transferases (GSTs) constitute a family of proteins encoded by a large supergene family and are widely expressed in animals, bacteria, fungi and plants, but with few reports of them in Euphorbiaceae plants. These proteins participate in and regulate the detoxification and oxidative stress response of heterogeneous organisms, resistance to stress, growth and development, signal transduction and other related processes. In this study, we identified and analyzed the whole genomes of four species of Euphorbiaceae, namely Ricinus communis, Jatropha curcas, Hevea brasiliensis, and Manihot esculenta, which have high economic and practical value. A total of 244 GST genes were identified. Based on their sequence characteristics and conserved domain types, the GST supergene family in Euphorbiaceae was classified into 10 subfamilies. The GST supergene families of Euphorbiaceae and Arabidopsis have been found to be highly conserved in evolution, and tandem repeats and translocations in these genes have made the greatest contributions to gene amplification here and have experienced strong purification selection. An evolutionary analysis showed that Euphorbiaceae GST genes have also evolved into new subtribes (GSTO, EF1BG, MAPEG), which may play a specific role in Euphorbiaceae. An analysis of expression patterns of the GST supergene family in Euphorbiaceae revealed the functions of these GSTs in different tissues, including resistance to stress and participation in herbicide detoxification. In addition, an interaction analysis was performed to determine the GST gene regulatory mechanism. The results of this study have laid a foundation for further analysis of the functions of the GST supergene family in Euphorbiaceae, especially in stress and herbicide detoxification. The results have also provided new ideas for the study of the regulatory mechanism of the GST supergene family, and have provided a reference for follow-up genetics and breeding work.

Genome-Wide Identification and Expression Analysis of Pseudouridine Synthase Family in Arabidopsis and Maize

Pseudouridine (Ψ), the isomer of uridine (U), is the most abundant type of RNA modification, which is crucial for gene regulation in various cellular processes. Pseudouridine synthases (PUSs) are the key enzymes for the U-to-Ψ conversion. However, little is known about the genome-wide features and biological function of plant PUSs. In this study, we identified 20 AtPUSs and 22 ZmPUSs from Arabidopsis and maize (Zea mays), respectively. Our phylogenetic analysis indicated that both AtPUSs and ZmPUSs could be clustered into six known subfamilies: RluA, RsuA, TruA, TruB, PUS10, and TruD. RluA subfamily is the largest subfamily in both Arabidopsis and maize. It's noteworthy that except the canonical XXHRLD-type RluAs, another three conserved RluA variants, including XXNRLD-, XXHQID-, and XXHRLG-type were also identified in those key nodes of vascular plants. Subcellular localization analysis of representative AtPUSs and ZmPUSs in each subfamily revealed that PUS proteins were localized in different organelles including nucleus, cytoplasm and chloroplasts. Transcriptional expression analysis indicated that AtPUSs and ZmPUSs were differentially expressed in various tissues and diversely responsive to abiotic stresses, especially suggesting their potential roles in response to heat and salt stresses. All these results would facilitate the functional identification of these pseudouridylation in the future.

Genome-wide identification of the Liriodendron chinense WRKY gene family and its diverse roles in response to multiple abiotic stress

BACKGROUND

Liriodendron chinense (Lchi) is a tree species within the Magnoliaceae family and is considered a basal angiosperm. The too low or high temperature or soil drought will restrict its growth as the adverse environmental conditions, thus improving L. chinense abiotic tolerance was the key issues to study. WRKYs are a major family of plant transcription factors known to often be involved in biotic and abiotic stress responses. So far, it is still largely unknown if and how the LchiWRKY gene family is tied to regulating L. chinense stress responses. Therefore, studying the involvement of the WRKY gene family in abiotic stress regulation in L. chinense could be very informative in showing how this tree deals with such stressful conditions.

RESULTS

In this research, we performed a genome-wide analysis of the Liriodendron chinense (Lchi) WRKY gene family, studying their classification relationships, gene structure, chromosomal locations, gene duplication, cis-element, and response to abiotic stress. The 44 members of the LchiWRKY gene family contain a significant amount of sequence diversity, with their lengths ranging from 525 bp to 40,981 bp. Using classification analysis, we divided the 44 LchiWRKY genes into three phylogenetic groups (I, II, II), with group II then being further divided into five subgroups (IIa, IIb, IIc, IId, IIe). Comparative phylogenetic analysis including the WRKY families from 17 plant species suggested that LchiWRKYs are closely related to the Magnolia Cinnamomum kanehirae WRKY family, and has fewer family members than higher plants. We found the LchiWRKYs to be evenly distributed across 15 chromosomes, with their duplication events suggesting that tandem duplication may have played a major role in LchiWRKY gene expansion model. A Ka/Ks analysis indicated that they mainly underwent purifying selection and distributed in the group IId. Motif analysis showed that LchiWRKYs contained 20 motifs, and different phylogenetic groups contained conserved motif. Gene ontology (GO) analysis showed that LchiWRKYs were mainly enriched in two categories, i.e., biological process and molecular function. Two group IIc members (LchiWRKY10 and LchiWRKY37) contain unique WRKY element sequence variants (WRKYGKK and WRKYGKS). Gene structure analysis showed that most LchiWRKYs possess 3 exons and two different types of introns: the R- and V-type which are both contained within the WRKY domain (WD). Additional promoter cis-element analysis indicated that 12 cis-elements that play different functions in environmental adaptability occur across all LchiWRKY groups. Heat, cold, and drought stress mainly induced the expression of group II and I LchiWRKYs, some of which had undergone gene duplication during evolution, and more than half of which had three exons. LchiWRKY33 mainly responded to cold stress and LchiWRKY25 mainly responded to heat stress, and LchiWRKY18 mainly responded to drought stress, which was almost 4-fold highly expressed, while 5 LchiWRKYs (LchiWRKY5, LchiWRKY23, LchiWRKY14, LchiWRKY27, and LchiWRKY36) responded equally three stresses with more than 6-fold expression. Subcellular localization analysis showed that all LchiWRKYs were localized in the nucleus, and subcellular localization experiments of LchiWRKY18 and 36 also showed that these two transcription factors were expressed in the nucleus.

CONCLUSIONS

This study shows that in Liriodendron chinense, several WRKY genes like LchiWRKY33, LchiWRKY25, and LchiWRKY18, respond to cold or heat or drought stress, suggesting that they may indeed play a role in regulating the tree's response to such conditions. This information will prove a pivotal role in directing further studies on the function of the LchiWRKY gene family in abiotic stress response and provides a theoretical basis for popularizing afforestation in different regions of China.

Ralstonia solanacearum Effectors Localize to Diverse Organelles in Solanum Hosts

Phytopathogenic bacteria secrete type III effector (T3E) proteins directly into host plant cells. T3Es can interact with plant proteins and frequently manipulate plant host physiological or developmental processes. The proper subcellular localization of T3Es is critical for their ability to interact with plant targets, and knowledge of T3E localization can be informative for studies of effector function. Here we investigated the subcellular localization of 19 T3Es from the phytopathogenic bacteria Ralstonia pseudosolanacearum and Ralstonia solanacearum. Approximately 45% of effectors in our library localize to both the plant cell periphery and the nucleus, 15% exclusively to the cell periphery, 15% exclusively to the nucleus, and 25% to other organelles, including tonoplasts and peroxisomes. Using tomato hairy roots, we show that T3E localization is similar in both leaves and roots and is not impacted by Solanum species. We find that in silico prediction programs are frequently inaccurate, highlighting the value of in planta localization experiments. Our data suggest that Ralstonia targets a wide diversity of cellular organelles and provides a foundation for developing testable hypotheses about Ralstonia effector function.

Genome-Wide Analysis of SRNF Genes in Gossypium hirsutum Reveals the Role of GhSRNF18 in Primary Root Growth

Root systems are instrumental for water and nutrient uptake and the anchorage of plants in the soil. Root regulating GL2-interacting repressors (GIRs) contain a Short RING-like Zinc-Finger (SRNF) domain, but there has been no comprehensive characterization about this gene family in any plant species. Here, we renamed the GIR-like proteins as SRNF proteins due to their conserved domain and identified 140 SRNF genes from 16 plant species including 24 GhSRNF genes in Gossypium hirsutum. Phylogenetic analysis of the SRNFs revealed both similarities and divergences between five subfamilies. Notably, synteny analysis revealed that polyploidization and whole-genome duplication contribute to the expansion of the GhSRNF gene family. Various cis-acting regulatory elements were shown to be pertinent to light, phytohormone, defense responsive, and meristem regulation. Furthermore, GhSRNF2/15 were predominantly expressed in root, whereas the expression of GhSRNF18 is positively correlated with the primary root (PR) length in G. hirsutum, quantified by quantitative real-time PCR (qRT-PCR). Over-expression of GhSRNF18 in Arabidopsis and virus-induced gene silencing (VIGS) of GhSRNF18 in G. hirsutum has revealed the role of GhSRNF18 in PR growth. The over-expression of GhSRNF18 in Arabidopsis resulted in an increase of meristematic activities and auxin accumulations in PRs, which were consistent with the transcriptomic data. Our results suggested that GhSRNF18 positively regulates PR growth. This study increased our understanding of the SRNF gene family in plants and provided a novel rationale for the further investigation of cotton root morphogenesis regulated by the GhSRNFs.

Cytokinin glucosyl transferases, key regulators of cytokinin homeostasis, have potential value for wheat improvement

The cytokinins, which are N⁶ -substituted adenine derivatives, control key aspects of crop productivity. Cytokinin levels are controlled via biosynthesis by isopentenyl transferase (IPT), destruction by cytokinin oxidase/dehydrogenase (CKX), and inactivation via glucosylation by cytokinin glucosyl transferases (CGTs). While both yield components and tolerance to drought and related abiotic stressors have been positively addressed via manipulation of IPT and/or CKX expression, much less attention has been paid to the CGTs. As naming of the CGTs has been unclear, we suggest COGT, CNGT, CONGT and CNOGT to describe the O-, N- and dual function CGTs. As specific CGT mutants of both rice and arabidopsis showed impacts on yield components, we interrogated the wheat genome database, IWGSC RefSeq v1.0 & v2.0, to investigate wheat CGTs. Besides providing unambiguous names for the 53 wheat CGTs, we show their expression patterns in 70 developmental tissues and their response characteristics to various stress conditions by reviewing more than 1000 RNA-seq data sets. These revealed various patterns of responses and showed expression generally being more limited in reproductive tissues than in vegetative tissues. Multiple cis-regulatory elements are present in the 3 kb upstream of the start codons of the 53 CGTs. Elements associated with abscisic acid, light and methyl jasmonate are particularly over-represented, indicative of the responsiveness of CGTs to the environment. These data sets indicate that CGTs have potential value for wheat improvement and that these could be targeted in TILLING or gene editing wheat breeding programmes.

Structure and Expression Analysis of Sucrose Phosphate Synthase, Sucrose Synthase and Invertase Gene Families in Solanum lycopersicum

Sucrose phosphate synthase (SPS), sucrose synthase (SUS) and invertase (INV) are all encoded by multigene families. In tomato (Solanum lycopersicum), a comprehensive analysis of structure characteristics of these family genes is still lacking, and the functions of individual isoforms of these families are mostly unclear under stress. Here, the structure characteristics of the three families in tomato were analyzed; moreover, as a first step toward understanding the functions of isoforms of these proteins under stress, the tissue expression pattern and stress response of these genes were also investigated. The results showed that four SPS genes, six SUS genes and nineteen INV genes were identified in tomato. The subfamily differentiation of SlSPS and SlSUS might have completed before the split of monocotyledons and dicotyledons. The conserved motifs were mostly consistent within each protein family/subfamily. These genes demonstrated differential expressions among family members and tissues, and in response to polyethylene glycerol, NaCl, H₂O₂, abscisic acid or salicylic acid treatment. Our results suggest that each isoform of these families may have different functions in different tissues and under environmental stimuli. SlSPS1, SlSPS3, SlSUS1, SlSUS3, SlSUS4, SlINVAN5 and SlINVAN7 demonstrated consistent expression responses and may be the major genes responding to exogenous stimuli.

Recent Advances in the Prediction of Protein Structural Classes: Feature Descriptors and Machine Learning Algorithms

In the postgenomic age, rapid growth in the number of sequence-known proteins has been accompanied by much slower growth in the number of structure-known proteins (as a result of experimental limitations), and a widening gap between the two is evident. Because protein function is linked to protein structure, successful prediction of protein structure is of significant importance in protein function identification. Foreknowledge of protein structural class can help improve protein structure prediction with significant medical and pharmaceutical implications. Thus, a fast, suitable, reliable, and reasonable computational method for protein structural class prediction has become pivotal in bioinformatics. Here, we review recent efforts in protein structural class prediction from protein sequence, with particular attention paid to new feature descriptors, which extract information from protein sequence, and the use of machine learning algorithms in both feature selection and the construction of new classification models. These new feature descriptors include amino acid composition, sequence order, physicochemical properties, multiprofile Bayes, and secondary structure-based features. Machine learning methods, such as artificial neural networks (ANNs), support vector machine (SVM), K-nearest neighbor (KNN), random forest, deep learning, and examples of their application are discussed in detail. We also present our view on possible future directions, challenges, and opportunities for the applications of machine learning algorithms for prediction of protein structural classes.

Identify RNA-associated subcellular localizations based on multi-label learning using Chou's 5-steps rule

BACKGROUND

Biological functions of biomolecules rely on the cellular compartments where they are located in cells. Importantly, RNAs are assigned in specific locations of a cell, enabling the cell to implement diverse biochemical processes in the way of concurrency. However, lots of existing RNA subcellular localization classifiers only solve the problem of single-label classification. It is of great practical significance to expand RNA subcellular localization into multi-label classification problem.

RESULTS

In this study, we extract multi-label classification datasets about RNA-associated subcellular localizations on various types of RNAs, and then construct subcellular localization datasets on four RNA categories. In order to study Homo sapiens, we further establish human RNA subcellular localization datasets. Furthermore, we utilize different nucleotide property composition models to extract effective features to adequately represent the important information of nucleotide sequences. In the most critical part, we achieve a major challenge that is to fuse the multivariate information through multiple kernel learning based on Hilbert-Schmidt independence criterion. The optimal combined kernel can be put into an integration support vector machine model for identifying multi-label RNA subcellular localizations. Our method obtained excellent results of 0.703, 0.757, 0.787, and 0.800, respectively on four RNA data sets on average precision.

CONCLUSION

To be specific, our novel method performs outstanding rather than other prediction tools on novel benchmark datasets. Moreover, we establish user-friendly web server with the implementation of our method.

Using the Chou’s Pseudo Component to Predict the ncRNA Locations Based on the Improved K-Nearest Neighbor (iKNN) Classifier

Background:

The non-coding RNA identification at the organelle genome level is a challenging task. In our previous work, an ncRNA dataset with less than 80% sequence identity was built, and a method incorporating an increment of diversity combining with support vector machine method was proposed.

Objective:

Based on the ncRNA_361 dataset, a novel decision-making method-an improved KNN (iKNN) classifier was proposed.

Methods:

In this paper, based on the iKNN algorithm, the physicochemical features of nucleotides, the degeneracy of genetic codons, and topological secondary structure were selected to represent the effective ncRNA characters. Then, the incremental feature selection method was utilized to optimize the feature set.

Results:

The results of iKNN indicated that the decision-making method of mean value is distinctly superior to the traditional decision-making method of majority vote the Increment of Diversity Combining Support Vector Machine (ID-SVM). The iKNN algorithm achieved an overall accuracy of 97.368% in the jackknife test, when k=3.

Conclusion:

It should be noted that the triplets of the structure-sequence mode under reading frames not only contains the entire sequence information but also reflects whether the base was paired or not, and the secondary structural topological parameters further describe the ncRNA secondary structure on the spatial level. The ncRNA dataset and the iKNN classifier are freely available at http://202.207.14.87:8032/fuwu/iKNN/index.asp.

Genome-wide identification and expression analysis of the SR gene family in longan (Dimocarpus longan Lour.)

Longan (Dimocarpus longan Lour.) is an important commercial fruit tree in southern China. The embryogenesis of longan affects the quality and yield of fruit. A large number of alternative splicing events occurs during somatic embryogenesis (SE), which is regulated by serine/arginine-rich (SR) proteins. However, the functions of SR proteins in longan are poorly understood. In this study, 21 Dlo-SR gene family members belonging to six subfamilies were identified, among which Dlo-RSZ20a, Dlo-SR30, Dlo-SR17, Dlo-SR53 and Dlo-SR32 were localized in the nucleus, Dlo-RSZ20b, Dlo-RSZ20c, Dlo-RSZ20d, Dlo-SC18, Dlo-RS2Z29, Dlo-SCL41, and Dlo-SR33 were localized in chloroplasts, and Dlo-RS43, Dlo-SC33, Dlo-SC37, Dlo-RS2Z33, Dlo-RS2Z16, Dlo-RS2Z24, Dlo-SCL43, Dlo-SR112, and Dlo-SR59 were localized in the nucleus and chloroplasts. The Dlo-SR genes exhibited differential expression patterns in different tissues of longan. The transcript levels of Dlo-RSZ20a, Dlo-SC18, Dlo-RS2Z29, DLo-SR59, Dlo-SR53, and Dlo-SR17 were low in all analyzed tissues, whereas Dlo-RS43, Dlo-RS2Z16, Dlo-RS2Z24, and Dlo-SR30 were highly expressed in all tissues. To clarify their function during SE, the transcript levels of Dlo-SR genes were analyzed at different four stages of SE, comprising non-embryonic callus (NEC), friable-embryogenic callus (EC), incomplete compact pro-embryogenic culture (ICpEC) and globular embryo (GE). Interestingly, the transcript levels of Dlo-RS2Z29 and Dlo-SR112 were increased in embryogenic cells compared with the NEC stage, whereas transcript levels of Dlo-RSZ20a, Dlo-RS43, Dlo-SC37, and Dlo-RS2Z16 were especially increased at the GE stage compared with the other stages. Alternative splicing events of Dlo-SR mRNA precursors (pre-mRNAs) was detected during SE, with totals of 41, 29, 35, and 44 events detected during NEC, EC, ICpEC, and GE respectively. Protein–protein interaction analysis showed that SR proteins were capable of interaction with each other. The results indicate that the alternative splicing of Dlo-SR pre-mRNAs occurs during SE and that Dlo-SR proteins may interact to regulate embryogenesis of longan.

Isolation, characterization and expression analysis of stress responsive plant nuclear transcriptional factor subunit (NF-YB2) from commercial Saccharum hybrid and wild relative Erianthus arundinaceus

Plant nuclear factor (NF-Y) is a transcription activating factor, consisting of three subunits, and plays a key regulatory role in many stress-responsive mechanisms including drought and salinity stresses. NF-Ys function both as complex and individual subunits. Considering the importance of sugarcane as a commercial crop with high socio-economic importance and the crop being affected mostly by water deficit stress and salinity stress causing significant yield loss, nuclear transcriptional factor NF-YB2 was focused in this study. Plant nuclear factor subunit B2 from Erianthus arundinaceus (EaNF-YB2), a wild relative of sugarcane which is known for its drought and salinity stress tolerance, and commercial Saccharum hybrid Co 86032 (ShNF-YB2) was isolated and characterized. Both EaNF-YB2 and ShNF-YB2 genes are 543 bp long that encodes for a polypeptide of 180 amino acid residues. Comparison of EaNF-YB2 and ShNF-YB2 gene sequences revealed nucleotide substitutions at nine positions corresponding to three synonymous and six nonsynonymous amino acid substitutions that resulted in variations in physiochemical properties. However, multiple sequence alignment (MSA) of NF-YB2 proteins showed conservation of functionally important amino acid residues. In silico analysis revealed NF-YB2 to be a hydrophilic and intracellular protein, and EaNF-YB2 is thermally more stable than that of ShNF-YB2. Phylogenetic analysis suggested the lower rate of evolution of NF-YB2. Subcellular localization in sugarcane callus revealed NF-YB2 localization at nucleus that further evidenced it to be a transcription activation factor. Comparative RT-qPCR experiments showed a significantly higher level of NF-YB2 expression in E. arundinaceus when compared to that in the commercial Saccharum hybrid Co 86032 under drought and salinity stresses. Hence, EaNF-YB2 could be an ideal candidate gene, and its overexpression in sugarcane through genetic engineering approach might enhance tolerance to drought and salinity stresses.

Using the Chou's 5-steps rule, transient overexpression technique, subcellular location, and bioinformatic analysis to verify the function of Vitis vinifera O-methyltranferase 3 (VvOMT3) protein

3-Isobutyl-2-methoxypyrazine (IBMP) is an important odor compound that revives unripe grapes or poor-quality wine. The biosynthesis of IBMP in grape berries is under the catalysis of Vitis vinifera O-methyltranferase 3 (VvOMT3). The homologous verification in this paper was carried out with the transient overexpression technique. The results showed that both the expression levels of the VvOMT3 gene and the IBMP concentration in 'Red globe' grapes increased significantly, which suggested that VvOMT3 could function in the biosynthesis of IBMP. Based on β-glucuronidase (GUS) staining results, blue color was only observed in grape pulp, not in grape skin, which indicated that VvOMT3 was expressed in grape pulp. The outcomes of the subcellular location examination performed on the protoplasts of Arabidopsis thaliana showed that the VvOMT3 protein was located on the inner surface of the cytoplasmic membrane. In summary, the VvOMT3 enzyme may function at the inner surface of the cytoplasmic membrane of pulp cells during grape development. These results will provide a background for future research on the catalytic mechanisms of VvOMT3.

ECMPride: prediction of human extracellular matrix proteins based on the ideal dataset using hybrid features with domain evidence

Extracellular matrix (ECM) proteins play an essential role in various biological processes in multicellular organisms, and their abnormal regulation can lead to many diseases. For large-scale ECM protein identification, especially through proteomic-based techniques, a theoretical reference database of ECM proteins is required. In this study, based on the experimentally verified ECM datasets and by the integration of protein domain features and a machine learning model, we developed ECMPride, a flexible and scalable tool for predicting ECM proteins. ECMPride achieved excellent performance in predicting ECM proteins, with appropriate balanced accuracy and sensitivity, and the performance of ECMPride was shown to be superior to the previously developed tool. A new theoretical dataset of human ECM components was also established by applying ECMPride to all human entries in the SwissProt database, containing a significant number of putative ECM proteins as well as the abundant biological annotations. This dataset might serve as a valuable reference resource for ECM protein identification.

Blast resistance gene Pi54 over-expressed in rice to understand its cellular and sub-cellular localization and response to different pathogens

Rice blast resistance gene, Pi54 provides broad-spectrum resistance against different strains of Magnaporthe oryzae. Understanding the cellular localization of Pi54 protein is an essential step towards deciphering its place of interaction with the cognate Avr-gene. In this study, we investigated the sub-cellular localization of Pi54 with Green Fluorescent Protein (GFP) as a molecular tag through transient and stable expression in onion epidermal cells (Allium cepa) and susceptible japonica cultivar rice Taipei 309 (TP309), respectively. Confocal microscopy based observations of the onion epidermal cells revealed nucleus and cytoplasm specific GFP signals. In the stable transformed rice plants, GFP signal was recorded in the stomata, upper epidermal cells, mesophyll cells, vascular bundle, and walls of bundle sheath and bulliform cells of leaf tissues. These observations were further confirmed by Immunocytochemical studies. Using GFP specific antibodies, it was found that there was sufficient aggregation of GFP::Pi54protein in the cytoplasm of the leaf mesophyll cells and periphery of the epidermal cells. Interestingly, the transgenic lines developed in this study could show a moderate level of resistance to Xanthomonas oryzae and Rhizoctonia solani, the causal agents of the rice bacterial blight and sheath blight diseases, respectively. This study is a first detailed report, which emphasizes the cellular and subcellular distribution of the broad spectrum blast resistance gene Pi54 in rice and the impact of its constitutive expression towards resistance against other fungal and bacterial pathogens of rice.

The Genome-Wide Analysis of RALF-Like Genes in Strawberry (Wild and Cultivated) and Five Other Plant Species (Rosaceae)

The rapid alkalinization factor (RALF) gene family is essential for the plant growth and development. However, there is little known about these genes among Rosaceae species. Here, we identify 124 RALF-like genes from seven Rosaceae species, and 39 genes from Arabidopsis, totally 163 genes, divided into four clades according to the phylogenetic analysis, which includes 45 mature RALF genes from Rosaceae species. The YISY motif and RRXL cleavage site are typical features of true RALF genes, but some variants were detected in our study, such as YISP, YIST, NISY, YINY, YIGY, YVGY, FIGY, YIAY, and RRVM. Motif1 is widely distributed among all the clades. According to screening of cis-regulatory elements, GO annotation, expression sequence tags (EST), RNA-seq, and RT-qPCR, we reported that 24 RALF genes coding mature proteins related to tissue development, fungal infection, and hormone response. Purifying selection may play an important role in the evolutionary process of RALF-like genes among Rosaceae species according to the result from ka/ks. The tandem duplication event just occurs in four gene pairs (Fv-RALF9 and Fv-RALF10, Md-RALF7 and Md-RALF8, Pm-RALF2 and Pm-RALF8, and Pp-RALF11 and Pp-RALF14) from four Rosaceae species. Our research provides a wide overview of RALF-like genes in seven Rosaceae species involved in identification, classification, structure, expression, and evolution analysis.

The effector GpRbp-1 of Globodera pallida targets a nuclear HECT E3 ubiquitin ligase to modulate gene expression in the host

Plant-parasitic nematodes secrete effectors that manipulate plant cell morphology and physiology to achieve host invasion and establish permanent feeding sites. Effectors from the highly expanded SPRYSEC (SPRY domain with a signal peptide for secretion) family in potato cyst nematodes have been implicated in activation and suppression of plant immunity, but the mechanisms underlying these activities remain largely unexplored. To study the host mechanisms used by SPRYSEC effectors, we identified plant targets of GpRbp-1 from the potato cyst nematode Globodera pallida. Here, we show that GpRbp-1 interacts in yeast and in planta with a functional potato homologue of the Homology to E6-AP C-Terminus (HECT)-type ubiquitin E3 ligase UPL3, which is located in the nucleus. Potato lines lacking StUPL3 are not available, but the Arabidopsis mutant upl3-5 displaying a reduced UPL3 expression showed a consistently small but not significant decrease in susceptibility to cyst nematodes. We observed a major impact on the root transcriptome by the lower levels of AtUPL3 in the upl3-5 mutant, but surprisingly only in association with infections by cyst nematodes. To our knowledge, this is the first example that a HECT-type ubiquitin E3 ligase is targeted by a pathogen effector and that a member of this class of proteins specifically regulates gene expression under biotic stress conditions. Together, our data suggest that GpRbp-1 targets a specific component of the plant ubiquitination machinery to manipulate the stress response in host cells.

Higher expression of the strawberry xyloglucan endotransglucosylase/hydrolase genes FvXTH9 and FvXTH6 accelerates fruit ripening

Fruit softening in Fragaria (strawberry) is proposed to be associated with the modification of cell wall components such as xyloglucan by the action of cell wall-modifying enzymes. This study focuses on the in vitro and in vivo characterization of two recombinant xyloglucan endotransglucosylase/hydrolases (XTHs) from Fragaria vesca, FvXTH9 and FvXTH6. Mining of the publicly available F. vesca genome sequence yielded 28 putative XTH genes. FvXTH9 showed the highest expression level of all FvXTHs in a fruit transcriptome data set and was selected with the closely related FvXTH6 for further analysis. To investigate their role in fruit ripening in more detail, the coding sequences of FvXTH9 and FvXTH6 were cloned into the vector pYES2 and expressed in Saccharomyces cerevisiae. FvXTH9 and FvXTH6 displayed xyloglucan endotransglucosylase (XET) activity towards various acceptor substrates using xyloglucan as the donor substrate. Interestingly, FvXTH9 showed activity of mixed-linkage glucan:xyloglucan endotransglucosylase (MXE) and cellulose:xyloglucan endotransglucosylase (CXE). The optimum pH of both FvXTH9 and FvXTH6 was 6.5. The prediction of subcellular localization suggested localization to the secretory pathway, which was confirmed by localization studies in Nicotiana tabacum. Overexpression showed that Fragaria × ananassa fruits infiltrated with FvXTH9 and FvXTH6 ripened faster and showed decreased firmness compared with the empty vector control pBI121. Thus FvXTH9 and also FvXTH6 might promote strawberry fruit ripening by the modification of cell wall components.

Chemodiversity of the Glucosinolate-Myrosinase System at the Single Cell Type Resolution

Glucosinolates (GLSs) are a well-defined group of specialized metabolites, and like any other plant specialized metabolites, their presence does not directly affect the plant survival in terms of growth and development. However, specialized metabolites are essential to combat environmental stresses, such as pathogens and herbivores. GLSs naturally occur in many pungent plants in the order of Brassicales. To date, more than 200 different GLS structures have been characterized and their distribution differs from species to species. GLSs co-exist with classical and atypical myrosinases, which can hydrolyze GLS into an unstable aglycone thiohydroximate-O-sulfonate, which rearranges to produce different degradation products. GLSs, myrosinases, myrosinase interacting proteins, and GLS degradation products constitute the GLS-myrosinase (GM) system ("mustard oil bomb"). This review discusses the cellular and subcellular organization of the GM system, its chemodiversity, and functions in different cell types. Although there are many studies on the functions of GLSs and/or myrosinases at the tissue and whole plant levels, very few studies have focused on different single cell types. Single cell type studies will help to reveal specific functions that are missed at the tissue and organismal level. This review aims to highlight (1) recent progress in cellular and subcellular compartmentation of GLSs, myrosinases, and myrosinase interacting proteins; (2) molecular and biochemical diversity of GLSs and myrosinases; and (3) myrosinase interaction with its interacting proteins, and how it regulates the degradation of GLSs and thus the biological functions (e.g., plant defense against pathogens). Future prospects may include targeted approaches for engineering/breeding of plants and crops in the cell type-specific manner toward enhanced plant defense and nutrition.

PMLPR: A novel method for predicting subcellular localization based on recommender systems

The importance of protein subcellular localization problem is due to the importance of protein's functions in different cell parts. Moreover, prediction of subcellular locations helps to identify the potential molecular targets for drugs and has an important role in genome annotation. Most of the existing prediction methods assign only one location for each protein. But, since some proteins move between different subcellular locations, they can have multiple locations. In recent years, some multiple location predictors have been introduced. However, their performances are not accurate enough and there is much room for improvement. In this paper, we introduced a method, PMLPR, to predict locations for a protein. PMLPR predicts a list of locations for each protein based on recommender systems and it can properly overcome the multiple location prediction problem. For evaluating the performance of PMLPR, we considered six datasets RAT, FLY, HUMAN, Du et al., DBMLoc and Höglund. The performance of this algorithm is compared with six state-of-the-art algorithms, YLoc, WOLF-PSORT, prediction channel, MDLoc, Du et al. and MultiLoc2-HighRes. The results indicate that our proposed method is significantly superior on RAT and Fly proteins, and decent on HUMAN proteins. Moreover, on the datasets introduced by Du et al., DBMLoc and Höglund, PMLPR has comparable results. For the case study, we applied the algorithms on 8 proteins which are important in cancer research. The results of comparison with other methods indicate the efficiency of PMLPR.

Identification and expression analysis of EDR1-like genes in tobacco (Nicotiana tabacum) in response to Golovinomyces orontii

ENHANCED DISEASE RESISTANCE1 (EDR1) encodes a Raf-like mitogen-activated protein kinase, and it acts as a negative regulator of disease resistance and ethylene-induced senescence. Mutations in the EDR1 gene can enhance resistance to powdery mildew both in monocotyledonous and dicotyledonous plants. However, little is known about EDR1-like gene members from a genome-wide perspective in plants. In this study, the tobacco (Nicotiana tabacum)EDR1-like gene family was first systematically analyzed. We identified 19 EDR1-like genes in tobacco, and compared them to those from Arabidopsis, tomato and rice. Phylogenetic analyses divided the EDR1-like gene family into six clades, among them monocot and dicot plants were respectively divided into two sub-clades. NtEDR1-1A and NtEDR1-1B were classified into clade I in which the other members have been reported to negatively regulate plant resistance to powdery mildew. The expression patterns of tobacco EDR1-like genes were analyzed after plants were challenged by Golovinomyces orontii, and showed that several other EDR1-like genes were induced after infection, as well as NtEDR1-1A and NtEDR1-1B. Expression analysis showed that NtEDR1-13 and NtEDR1-16 had exclusively abundant expression patterns in roots and leaves, respectively, and the remaining NtEDR1-like members were actively expressed in most of the tissue/organ samples investigated. Our findings will contribute to further study of the physiological functions of EDR1-like genes in tobacco.

HPSLPred: An Ensemble Multi-Label Classifier for Human Protein Subcellular Location Prediction with Imbalanced Source

Predicting the subcellular localization of proteins is an important and challenging problem. Traditional experimental approaches are often expensive and time-consuming. Consequently, a growing number of research efforts employ a series of machine learning approaches to predict the subcellular location of proteins. There are two main challenges among the state-of-the-art prediction methods. First, most of the existing techniques are designed to deal with multi-class rather than multi-label classification, which ignores connections between multiple labels. In reality, multiple locations of particular proteins imply that there are vital and unique biological significances that deserve special focus and cannot be ignored. Second, techniques for handling imbalanced data in multi-label classification problems are necessary, but never employed. For solving these two issues, we have developed an ensemble multi-label classifier called HPSLPred, which can be applied for multi-label classification with an imbalanced protein source. For convenience, a user-friendly webserver has been established at http://server.malab.cn/HPSLPred.

2L-piRNA: A Two-Layer Ensemble Classifier for Identifying Piwi-Interacting RNAs and Their Function

Involved with important cellular or gene functions and implicated with many kinds of cancers, piRNAs, or piwi-interacting RNAs, are of small non-coding RNA with around 19–33 nt in length. Given a small non-coding RNA molecule, can we predict whether it is of piRNA according to its sequence information alone? Furthermore, there are two types of piRNA: one has the function of instructing target mRNA deadenylation, and the other does not. Can we discriminate one from the other? With the avalanche of RNA sequences emerging in the postgenomic age, it is urgent to address the two problems for both basic research and drug development. Unfortunately, to the best of our knowledge, so far no computational methods whatsoever could be used to deal with the second problem, let alone deal with the two problems together. Here, by incorporating the physicochemical properties of nucleotides into the pseudo K-tuple nucleotide composition (PseKNC), we proposed a powerful predictor called 2L-piRNA. It is a two-layer ensemble classifier, in which the first layer is for identifying whether a query RNA molecule is piRNA or non-piRNA, and the second layer for identifying whether a piRNA is with or without the function of instructing target mRNA deadenylation. Rigorous cross-validations have indicated that the success rates achieved by the proposed predictor are quite high. For the convenience of most biologists and drug development scientists, the web server for 2L-piRNA has been established at http://bioinformatics.hitsz.edu.cn/2L-piRNA/, by which users can easily get their desired results without the need to go through the mathematical details.

iATC-mHyb: a hybrid multi-label classifier for predicting the classification of anatomical therapeutic chemicals

Recommended by the World Health Organization (WHO), drug compounds have been classified into 14 main ATC (Anatomical Therapeutic Chemical) classes according to their therapeutic and chemical characteristics. Given an uncharacterized compound, can we develop a computational method to fast identify which ATC class or classes it belongs to? The information thus obtained will timely help adjusting our focus and selection, significantly speeding up the drug development process. But this problem is by no means an easy one since some drug compounds may belong to two or more than two ATC classes. To address this problem, using the DO (Drug Ontology) approach based on the ChEBI (Chemical Entities of Biological Interest) database, we developed a predictor called iATC-mDO. Subsequently, hybridizing it with an existing drug ATC classifier, we constructed a predictor called iATC-mHyb. It has been demonstrated by the rigorous cross-validation and from five different measuring angles that iATC-mHyb is remarkably superior to the best existing predictor in identifying the ATC classes for drug compounds. To convenience most experimental scientists, a user-friendly web-server for iATC-mHyd has been established at http://www.jci-bioinfo.cn/iATC-mHyb, by which users can easily get their desired results without the need to go through the complicated mathematical equations involved.

Proteomic Profiling of the Microsomal Root Fraction: Discrimination of Pisum sativum L. Cultivars and Identification of Putative Root Growth Markers

Legumes are a large and economically important family, containing a variety of crop plants. Alongside different cereals, some fruits, and tropical roots, a number of leguminosae evolved for millennia as crops with human society. One of these legumes is Pisum sativum L., the common garden pea. In the past, breeding has been largely selective on improved above-ground organs. However, parameters, such as root-growth, which determines acquisition of nutrients and water, have largely been underestimated. Although the genome of P. sativum is still not fully sequenced, multiple proteomic studies have been published on a variety of physiological aspects in the last years. The presented work focused on the connection between root length and the influence of the microsomal root proteome of four different pea cultivars after five days of germination (cultivar Vroege, Girl from the Rhineland, Kelvedon Wonder, and Blauwschokker). In total, 60 proteins were identified to have significantly differential abundances in the four cultivars. Root growth of five-days old seedlings and their microsomal proteome revealed a similar separation pattern, suggesting that cultivar-specific root growth performance is explained by differential membrane and ribosomal protein levels. Hence, we reveal and discuss several putative root growth protein markers possibly playing a key role for improved primary root growth breeding strategies.

Classify vertebrate hemoglobin proteins by incorporating the evolutionary information into the general PseAAC with the hybrid approach

Hemoglobin is an oxygen-binding protein widely present in all kingdoms of life from prokaryotic to eukaryotic, but well established in the vertebrate system. An attempt was made to determine the Vertebrate hemoglobin (VerHb) protein on their animal classifications, based on general pseudo amino acid composition (PseAAC)'s evolutionary profiles and hybrid approach. The support vector machine (SVM) has been applied to develop all models, the prediction results further compared according to their animal classification. The performance of the approaches estimated using five-fold cross-validation techniques. The prediction performance was further investigated by receiver operating characteristic (ROC) and prediction score graphs. The prediction accuracy (ACC), sensitivity (SN) and specificity (SP) were examined to find the accurate predictions on the threshold level. Based on the approach, a web-tool has been developed for identifying the VerHb proteins.

Nuclear Proteomics Reveals the Role of Protein Synthesis and Chromatin Structure in Root Tip of Soybean during the Initial Stage of Flooding Stress

To identify the upstream events controlling the regulation of flooding-responsive proteins in soybean, proteomic analysis of nuclear proteins in root tip was performed. By using nuclear fractions, which were highly enriched, a total of 365 nuclear proteins were changed in soybean root tip at initial stage of flooding stress. Four exon-junction complex-related proteins and NOP1/NOP56, which function in upstream of 60S preribosome biogenesis, were decreased in flooded soybean. Furthermore, proteomic analysis of crude protein extract revealed that the protein translation was suppressed by continuous flooding stress. Seventeen chromatin structure-related nuclear proteins were decreased in response to flooding stress. Out of them, histone H3 was clearly decreased with protein abundance and mRNA expression levels at the initial flooding stress. Additionally, a number of protein synthesis-, RNA-, and DNA-related nuclear proteins were decreased in a time-dependent manner. mRNA expressions of genes encoding the significantly changed flooding-responsive nuclear proteins were inhibited by the transcriptional inhibitor, actinomycin D. These results suggest that protein translation is suppressed through inhibition of preribosome biogenesis- and mRNA processing-related proteins in nuclei of soybean root tip at initial flooding stress. In addition, flooding stress may regulate histone variants with gene expression in root tip.

Human Protein Subcellular Localization with Integrated Source and Multi-label Ensemble Classifier

Predicting protein subcellular location is necessary for understanding cell function. Several machine learning methods have been developed for computational prediction of primary protein sequences because wet experiments are costly and time consuming. However, two problems still exist in state-of-the-art methods. First, several proteins appear in different subcellular structures simultaneously, whereas current methods only predict one protein sequence in one subcellular structure. Second, most software tools are trained with obsolete data and the latest new databases are missed. We proposed a novel multi-label classification algorithm to solve the first problem and integrated several latest databases to improve prediction performance. Experiments proved the effectiveness of the proposed method. The present study would facilitate research on cellular proteomics.

Over-Expression of the Pikh Gene with a CaMV 35S Promoter Leads to Improved Blast Disease (Magnaporthe oryzae) Tolerance in Rice

Magnaporthe oryzae is a rice blast fungus and plant pathogen that causes a serious rice disease and, therefore, poses a threat to the world's second most important food security crop. Plant transformation technology has become an adaptable system for cultivar improvement and to functionally analyze genes in plants. The objective of this study was to determine the effects (through over-expressing and using the CaMV 35S promoter) of Pikh on MR219 resistance because it is a rice variety that is susceptible to the blast fungus pathotype P7.2. Thus, a full DNA and coding DNA sequence (CDS) of the Pikh gene, 3172 bp, and 1206 bp in length, were obtained through amplifying the gDNA and cDNA template from a PH9-resistant rice variety using a specific primer. Agrobacterium-mediated transformation technology was also used to introduce the Pikh gene into the MR219 callus. Subsequently, transgenic plants were evaluated from the DNA to protein stages using polymerase chain reaction (PCR), semi-quantitative RT-PCR, real-time quantitative PCR and high performance liquid chromatography (HPLC). Transgenic plants were also compared with a control using a real-time quantification technique (to quantify the pathogen population), and transgenic and control plants were challenged with the local most virulent M. oryzae pathotype, P7.2. Based on the results, the Pikh gene encodes a hydrophilic protein with 18 sheets, 4 helixes, and 21 coils. This protein contains 401 amino acids, among which the amino acid sequence from 1 to 376 is a non-cytoplasmic region, that from 377 to 397 is a transmembrane region, and that from 398 to 401 is a cytoplasmic region with no identified disordered regions. The Pikh gene was up-regulated in the transgenic plants compared with the control plants. The quantity of the amino acid leucine in the transgenic rice plants increased significantly from 17.131 in the wild-type to 47.865 mg g(-1) in transgenic plants. The M. oryzae population was constant at 31, 48, and 72 h after inoculation in transgenic plants, while it was increased in the inoculated control plants. This study successfully clarified that over-expression of the Pikh gene in transgenic plants can improve their blast resistance against the M. oryzae pathotype P7.2.

C-terminal extension of calmodulin-like 3 protein from Oryza sativa L.: interaction with a high mobility group target protein

A large number of calmodulin-like (CML) proteins are present in plants, but there is little detailed information on the functions of these proteins in rice (Oryza sativa L.). Here, the CML3 protein from rice (OsCML3) and its truncated form lacking the C-terminal extension (OsCML3m) were found to exhibit a Ca2+-binding property and subsequent conformational change, but the ability to bind the CaM kinase II peptide was only observed for OsCML3m. Changes in their secondary structure upon Ca2+-binding measured by circular dichroism revealed that OsCML3m had a higher helical content than OsCML3. Moreover, OsCML3 was mainly localized in the plasma membrane, whereas OsCML3m was found in the nucleus. The rice high mobility group B1 (OsHMGB1) protein was identified as one of the putative OsCML3 target proteins. Bimolecular fluorescence complementation analysis revealed that OsHMGB1 bound OsCML3, OsCML3m or OsCML3s (cysteine to serine mutation at the prenylation site) in the nucleus presumably through the methionine and phenylalanine-rich hydrophobic patches, confirming that OsHMGB1 is a target protein in planta. The effect of OsCML3 or OsCML3m on the DNA-binding ability of OsHMGB1 was measured using an electrophoretic mobility shift assay. OsCML3m decreased the level of OsHMGB1 binding to pUC19 double-stranded DNA whereas OsCML3 did not. Taken together, OsCML3 probably provides a mechanism for manipulating the DNA-binding ability of OsHMGB1 in the nucleus and its C-terminal extension provides an intracellular Ca2+ regulatory switch.

The DnaJ OsDjA7/8 is essential for chloroplast development in rice (Oryza sativa)

DnaJ proteins belong to chaperones of Hsp40 family that ubiquitously participate in various cellular processes. Previous studies have shown chloroplast-targeted DnaJs are involved in the development of chloroplast in some plant species. However, little is known about the function of DnaJs in rice, one of the main staple crops. In this study, we characterized a type I DnaJ protein OsDjA7/8. We found that the gene OsDjA7/8 was expressed in all collected tissues, with a priority in the vigorous growth leaf. Subcellular localization revealed that the protein OsDjA7/8 was mainly distributed in chloroplast. Reduced expression of OsDjA7/8 in rice led to albino lethal at the seedling stage. Transmission electron microscopy observation showed that the chloroplast structures were abnormally developed in the plants silenced for OsDjA7/8. In addition, the transcriptional expression of the genes tightly associated with the development of chloroplast was deeply reduced in the plants silenced for OsDjA7/8. Collectively, our study reveals that OsDjA7/8 encodes a chloroplast-localized protein and is essential for chloroplast development and differentiation in rice.