Prediction of subcellular localization using sequence-biased recurrent networks

Multiplex CRISPR/Cas9-mediated genome editing to address drought tolerance in wheat

Genome editing tools have rapidly been adopted by plant scientists for crop improvement. Genome editing using a multiplex sgRNA-CRISPR/Cas9 genome editing system is a useful technique for crop improvement in monocot species. In this study, we utilized precise gene editing techniques to generate wheat 3'(2'), 5'-bisphosphate nucleotidase (TaSal1) mutants using a multiplex sgRNA-CRISPR/Cas9 genome editing system. Five active TaSal1 homologous genes were found in the genome of Giza168 in addition to another apparently inactive gene on chromosome 4A. Three gRNAs were designed and used to target exons 4, 5 and 7 of the five wheat TaSal1 genes. Among the 120 Giza168 transgenic plants, 41 lines exhibited mutations and produced heritable TaSal1 mutations in the M1 progeny and 5 lines were full 5 gene knock-outs. These mutant plants exhibit a rolled-leaf phenotype in young leaves and bended stems, but there were no significant changes in the internode length and width, leaf morphology, and stem shape. Anatomical and scanning electron microscope studies of the young leaves of mutated TaSal1 lines showed closed stomata, increased stomata width and increase in the size of the bulliform cells. Sal1 mutant seedlings germinated and grew better on media containing polyethylene glycol than wildtype seedlings. Our results indicate that the application of the multiplex sgRNA-CRISPR/Cas9 genome editing is efficient tool for mutating more multiple TaSal1 loci in hexaploid wheat.

Genome-wide identification of nitrate transporter 1/peptide transporter family (NPF) induced by arbuscular mycorrhiza in the maize genome

The Transporter 1/Peptide Transporter Family (NPF) is essential for the uptake and transport of nitrate nitrogen. Significant increases in nitrogen have been increasingly reported for many mycorrhizal plants, but there are few reports on maize. Here, we have identified the maize NPF family and screened for arbuscular mycorrhiza fungi (AMF) induced NPFs. In this study, a systematic analysis of the maize NPF gene family was performed. A total of 82 NPF genes were identified in maize. ZmNPF4.5 was strongly induced by AMF in both low and high nitrogen. Lotus japonicus hairy root-induced transformation experiments showed that ZmNPF4.5 promoter-driven GUS activity was restricted to cells containing tufts. Yeast backfill experiments indicate that ZmNPF4.5 functions in nitrate uptake. Therefore, we speculate that ZmNPF4.5 is a key gene for nitrate-nitrogen uptake in maize through the mycorrhizal pathway. This is a reference value for further exploring the acquisition of nitrate-nitrogen by maize through AMF pathway.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01464-3.

A review from biological mapping to computation-based subcellular localization

Subcellular localization is crucial to the study of virus and diseases. Specifically, research on protein subcellular localization can help identify clues between virus and host cells that can aid in the design of targeted drugs. Research on RNA subcellular localization is significant for human diseases (such as Alzheimer's disease, colon cancer, etc.). To date, only reviews addressing subcellular localization of proteins have been published, which are outdated for reference, and reviews of RNA subcellular localization are not comprehensive. Therefore, we collated (the most up-to-date) literature on protein and RNA subcellular localization to help researchers understand changes in the field of protein and RNA subcellular localization. Extensive and complete methods for constructing subcellular localization models have also been summarized, which can help readers understand the changes in application of biotechnology and computer science in subcellular localization research and explore how to use biological data to construct improved subcellular localization models. This paper is the first review to cover both protein subcellular localization and RNA subcellular localization. We urge researchers from biology and computational biology to jointly pay attention to transformation patterns, interrelationships, differences, and causality of protein subcellular localization and RNA subcellular localization.

Phylogenetic analysis and structural prediction reveal the potential functional diversity between green algae SWEET transporters

Sugar-Will-Eventually-be-Exported-Transporters (SWEETs) are an important family of sugar transporters that appear to be ubiquitous in all organisms. Recent research has determined the structure of SWEETs in higher plants, identified specific residues required for monosaccharide or disaccharide transport, and begun to understand the specific functions of individual plant SWEET proteins. However, in green algae (Chlorophyta) these transporters are poorly characterised. This study identified SWEET proteins from across representative Chlorophyta with the aim to characterise their phylogenetic relationships and perform protein structure modelling in order to inform functional prediction. The algal genomes analysed encoded between one and six SWEET proteins, which is much less than a typical higher plant. Phylogenetic analysis identified distinct clusters of over 70 SWEET protein sequences, taken from almost 30 algal genomes. These clusters remain separate from representative higher or non-vascular plant SWEETs, but are close to fungi SWEETs. Subcellular localisation predictions and analysis of conserved amino acid residues revealed variation between SWEET proteins of different clusters, suggesting different functionality. These findings also showed conservation of key residues at the substrate-binding site, indicating a similar mechanism of substrate selectivity and transport to previously characterised higher plant monosaccharide-transporting SWEET proteins. Future work is now required to confirm the predicted sugar transport specificity and determine the functional role of these algal SWEET proteins.

In Silico Characterisation of the Late Embryogenesis Abundant (LEA) Protein Families and Their Role in Desiccation Tolerance in Ramonda serbica Panc

Ramonda serbica Panc. is an ancient resurrection plant able to survive a long desiccation period and recover metabolic functions upon watering. The accumulation of protective late embryogenesis abundant proteins (LEAPs) is a desiccation tolerance hallmark. To propose their role in R. serbica desiccation tolerance, we structurally characterised LEAPs and evaluated LEA gene expression levels in hydrated and desiccated leaves. By integrating de novo transcriptomics and homologues LEAP domains, 318 R. serbica LEAPs were identified and classified according to their conserved motifs and phylogeny. The in silico analysis revealed that hydrophilic LEA4 proteins exhibited an exceptionally high tendency to form amphipathic α-helices. The most abundant, atypical LEA2 group contained more hydrophobic proteins predicted to fold into the defined globular domains. Within the desiccation-upregulated LEA genes, the majority encoded highly disordered DEH1, LEA1, LEA4.2, and LEA4.3 proteins, while the greatest portion of downregulated genes encoded LEA2.3 and LEA2.5 proteins. While dehydrins might chelate metals and bind DNA under water deficit, other intrinsically disordered LEAPs might participate in forming intracellular proteinaceous condensates or adopt amphipathic α-helical conformation, enabling them to stabilise desiccation-sensitive proteins and membranes. This comprehensive LEAPs structural characterisation is essential to understanding their function and regulation during desiccation aiming at crop drought tolerance improvement.

Auxin Metabolome Profiling in the Arabidopsis Endoplasmic Reticulum Using an Optimised Organelle Isolation Protocol

The endoplasmic reticulum (ER) is an extensive network of intracellular membranes. Its major functions include proteosynthesis, protein folding, post-transcriptional modification and sorting of proteins within the cell, and lipid anabolism. Moreover, several studies have suggested that it may be involved in regulating intracellular auxin homeostasis in plants by modulating its metabolism. Therefore, to study auxin metabolome in the ER, it is necessary to obtain a highly enriched (ideally, pure) ER fraction. Isolation of the ER is challenging because its biochemical properties are very similar to those of other cellular endomembranes. Most published protocols for ER isolation use density gradient ultracentrifugation, despite its suboptimal resolving power. Here we present an optimised protocol for ER isolation from Arabidopsis thaliana seedlings for the subsequent mass spectrometric determination of ER-specific auxin metabolite profiles. Auxin metabolite analysis revealed highly elevated levels of active auxin form (IAA) within the ER compared to whole plants. Moreover, samples prepared using our optimised isolation ER protocol are amenable to analysis using various “omics” technologies including analyses of both macromolecular and low molecular weight compounds from the same sample.

Plasma membrane phylloquinone biosynthesis in nonphotosynthetic parasitic plants

Nonphotosynthetic holoparasites exploit flexible targeting of phylloquinone biosynthesis to facilitate plasma membrane redox signaling. Phylloquinone is a lipophilic naphthoquinone found predominantly in chloroplasts and best known for its function in photosystem I electron transport and disulfide bridge formation of photosystem II subunits. Phylloquinone has also been detected in plasma membrane (PM) preparations of heterotrophic tissues with potential transmembrane redox function, but the molecular basis for this noncanonical pathway is unknown. Here, we provide evidence of PM phylloquinone biosynthesis in a nonphotosynthetic holoparasite Phelipanche aegyptiaca. A nonphotosynthetic and nonplastidial role for phylloquinone is supported by transcription of phylloquinone biosynthetic genes during seed germination and haustorium development, by PM-localization of alternative terminal enzymes, and by detection of phylloquinone in germinated seeds. Comparative gene network analysis with photosynthetically competent parasites revealed a bias of P. aegyptiaca phylloquinone genes toward coexpression with oxidoreductases involved in PM electron transport. Genes encoding the PM phylloquinone pathway are also present in several photoautotrophic taxa of Asterids, suggesting an ancient origin of multifunctionality. Our findings suggest that nonphotosynthetic holoparasites exploit alternative targeting of phylloquinone for transmembrane redox signaling associated with parasitism.

Genome-wide characterization, expression analyses, and functional prediction of the NPF family in Brassica napus

Background

NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER (NRT1/PTR) family (NPF) members are essential transporters for many substrates in plants, including nitrate, hormones, peptides, and secondary metabolites. Here, we report the global characterization of NPF in the important oil crop Brassica napus, including that for phylogeny, gene/protein structures, duplications, and expression patterns.

Results

A total of 199 B. napus (BnaNPFs) NPF-coding genes were identified. Phylogenetic analyses categorized these genes into 11 subfamilies, including three new ones. Sequence feature analysis revealed that members of each subfamily contain conserved gene and protein structures. Many hormone−/abiotic stress-responsive cis-acting elements and transcription factor binding sites were identified in BnaNPF promoter regions. Chromosome distribution analysis indicated that BnaNPFs within a subfamily tend to cluster on one chromosome. Syntenic relationship analysis showed that allotetraploid creation by its ancestors (Brassica rapa and Brassica oleracea) (57.89%) and small-scale duplication events (39.85%) contributed to rapid BnaNPF expansion in B. napus. A genome-wide spatiotemporal expression survey showed that NPF genes of each Arabidopsis and B. napus subfamily have preferential expression patterns across developmental stages, most of them are expressed in a few organs. RNA-seq analysis showed that many BnaNPFs (32.66%) have wide exogenous hormone-inductive profiles, suggesting important hormone-mediated patterns in diverse bioprocesses. Homologs in a clade or branch within a given subfamily have conserved organ/spatiotemporal and hormone-inductive profiles, indicating functional conservation during evolution. qRT-PCR-based comparative expression analysis of the 12 BnaNPFs in the NPF2–1 subfamily between high- and low-glucosinolate (GLS) content B. napus varieties revealed that homologs of AtNPF2.9 (BnaNPF2.12, BnaNPF2.13, and BnaNPF2.14), AtNPF2.10 (BnaNPF2.19 and BnaNPF2.20), and AtNPF2.11 (BnaNPF2.26 and BnaNPF2.28) might be involved in GLS transport. qRT-PCR further confirmed the hormone-responsive expression profiles of these putative GLS transporter genes.

Conclusion

We identified 199 B. napus BnaNPFs; these were divided into 11 subfamilies. Allopolyploidy and small-scale duplication events contributed to the immense expansion of BnaNPFs in B. napus. The BnaNPFs had preferential expression patterns in different tissues/organs and wide hormone-induced expression profiles. Four BnaNPFs in the NPF2–1 subfamily may be involved in GLS transport. Our results provide an abundant gene resource for further functional analysis of BnaNPFs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07274-7.

Genome-wide identification of the restorer-of-fertility-like (RFL) gene family in Brassica napus and expression analysis in Shaan2A cytoplasmic male sterility

Background

Cytoplasmic male sterility (CMS) is very important in hybrid breeding. The restorer-of-fertility (Rf) nuclear genes rescue the sterile phenotype. Most of the Rf genes encode pentatricopeptide repeat (PPR) proteins.

Results

We investigated the restorer-of-fertility-like (RFL) gene family in Brassica napus. A total of 53 BnRFL genes were identified. While most of the BnRFL genes were distributed on 10 of the 19 chromosomes, gene clusters were identified on chromosomes A9 and C8. The number of PPR motifs in the BnRFL proteins varied from 2 to 19, and the majority of BnRFL proteins harbored more than 10 PPR motifs. An interaction network analysis was performed to predict the interacting partners of RFL proteins. Tissue-specific expression and RNA-seq analyses between the restorer line KC01 and the sterile line Shaan2A indicated that BnRFL1, BnRFL5, BnRFL6, BnRFL8, BnRFL11, BnRFL13 and BnRFL42 located in gene clusters on chromosomes A9 and C8 were highly expressed in KC01.

Conclusions

In the present study, identification and gene expression analysis of RFL gene family in the CMS system were conducted, and seven BnRFL genes were identified as candidates for the restorer genes in Shaan2A CMS. Taken together, this method might provide new insight into the study of Rf genes in other CMS systems.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07163-z.

A Genome-Wide Survey of MATE Transporters in Brassicaceae and Unveiling Their Expression Profiles under Abiotic Stress in Rapeseed

The multidrug and toxic compound extrusion (MATE) protein family is important in the export of toxins and other substrates, but detailed information on this family in the Brassicaceae has not yet been reported compared to Arabidopsis thaliana. In this study, we identified 57, 124, 81, 85, 130, and 79 MATE genes in A. thaliana, Brassica napus, Brassica oleracea, Brassica rapa, Brassica juncea, and Brassica nigra, respectively, which were unevenly distributed on chromosomes owing to both tandem and segmental duplication events. Phylogenetic analysis showed that these genes could be classified into four subgroups, shared high similarity and conservation within each group, and have evolved mainly through purifying selection. Furthermore, numerous B. napusMATE genes showed differential expression between tissues and developmental stages and between plants treated with heavy metals or hormones and untreated control plants. This differential expression was especially pronounced for the Group 2 and 3 BnaMATE genes, indicating that they may play important roles in stress tolerance and hormone induction. Our results provide a valuable foundation for the functional dissection of the different BnaMATE homologs in B. napus and its parental lines, as well as for the breeding of more stress-tolerant B. napus genotypes.

Detecting sequence signals in targeting peptides using deep learning

In bioinformatics, machine learning methods have been used to predict features embedded in the sequences. In contrast to what is generally assumed, machine learning approaches can also provide new insights into the underlying biology. Here, we demonstrate this by presenting TargetP 2.0, a novel state-of-the-art method to identify N-terminal sorting signals, which direct proteins to the secretory pathway, mitochondria, and chloroplasts or other plastids. By examining the strongest signals from the attention layer in the network, we find that the second residue in the protein, that is, the one following the initial methionine, has a strong influence on the classification. We observe that two-thirds of chloroplast and thylakoid transit peptides have an alanine in position 2, compared with 20% in other plant proteins. We also note that in fungi and single-celled eukaryotes, less than 30% of the targeting peptides have an amino acid that allows the removal of the N-terminal methionine compared with 60% for the proteins without targeting peptide. The importance of this feature for predictions has not been highlighted before.

Detecting sequence signals in targeting peptides using deep learning

In bioinformatics, machine learning methods have been used to predict features embedded in the sequences. In contrast to what is generally assumed, machine learning approaches can also provide new insights into the underlying biology. Here, we demonstrate this by presenting TargetP 2.0, a novel state-of-the-art method to identify N-terminal sorting signals, which direct proteins to the secretory pathway, mitochondria, and chloroplasts or other plastids. By examining the strongest signals from the attention layer in the network, we find that the second residue in the protein, that is, the one following the initial methionine, has a strong influence on the classification. We observe that two-thirds of chloroplast and thylakoid transit peptides have an alanine in position 2, compared with 20% in other plant proteins. We also note that in fungi and single-celled eukaryotes, less than 30% of the targeting peptides have an amino acid that allows the removal of the N-terminal methionine compared with 60% for the proteins without targeting peptide. The importance of this feature for predictions has not been highlighted before.

Molecular Characterization of ZosmaNRT2, the Putative Sodium Dependent High-Affinity Nitrate Transporter of Zostera marina L

One of the most important adaptations of seagrasses during sea colonization was the capacity to grow at the low micromolar nitrate concentrations present in the sea. In contrast to terrestrial plants that use H⁺ symporters for high-affinity NO₃⁻ uptake, seagrasses such as Zostera marina L. use a Na⁺-dependent high-affinity nitrate transporter. Interestingly, in the Z. marina genome, only one gene (Zosma70g00300.1; NRT2.1) is annotated to this function. Analysis of this sequence predicts the presence of 12 transmembrane domains, including the MFS domains of the NNP transporter family and the “nitrate signature” that appears in all members of the NNP family. Phylogenetic analysis shows that this sequence is more related to NRT2.5 than to NRT2.1, sharing a common ancestor with both monocot and dicot plants. Heterologous expression of ZosmaNRT2-GFP together with the high-affinity nitrate transporter accessory protein ZosmaNAR2 (Zosma63g00220.1) in Nicotiana benthamiana leaves displayed four-fold higher fluorescence intensity than single expression of ZosmaNRT2-GFP suggesting the stabilization of NRT2 by NAR2. ZosmaNRT2-GFP signal was present on the Hechtian-strands in the plasmolyzed cells, pointing that ZosmaNRT2 is localized on the plasma membrane and that would be stabilized by ZosmaNAR2. Taken together, these results suggest that Zosma70g00300.1 would encode a high-affinity nitrate transporter located at the plasma membrane, equivalent to NRT2.5 transporters. These molecular data, together with our previous electrophysiological results support that ZosmaNRT2 would have evolved to use Na⁺ as a driving ion, which might be an essential adaptation of seagrasses to colonize marine environments.

Functional characterization of Gh_A08G1120 (GH3.5) gene reveal their significant role in enhancing drought and salt stress tolerance in cotton

BACKGROUND

Auxins play an important role in plant growth and development; the auxins responsive gene; auxin/indole-3-acetic acid (Aux/IAA), small auxin-up RNAs (SAUR) and Gretchen Hagen3 (GH3) control their mechanisms. The GH3 genes function in homeostasis by the catalytic activities in auxin conjugation and bounding free indole-3-acetic acid (IAA) to amino acids.

RESULTS

In our study, we identified the GH3 genes in three cotton species; Gossypium hirsutum, Gossypium arboreum and Gossypium raimondii, analyzed their chromosomal distribution, phylogenetic relationships, cis-regulatory element function and performed virus induced gene silencing of the novel Gh_A08G1120 (GH3.5) gene. The phylogenetic tree showed four clusters of genes with clade 1, 3 and 4 having mainly members of the GH3 of the cotton species while clade 2 was mainly members belonging to Arabidopsis. There were no paralogous genes, and few orthologous genes were observed between Gossypium and other species. All the GO terms were detected, but only 14 genes were found to have described GO terms in upland cotton, more biological functions were detected, as compared to the other functions. The GH3.17 subfamily harbored the highest number of the cis-regulatory elements, most having promoters towards dehydration-responsiveness. The RNA expression analysis revealed that 10 and 8 genes in drought and salinity stress conditions respectively were upregulated in G. hirsutum. All the genes that were upregulated in plants under salt stress conditions were also upregulated in drought stress; moreover, Gh_A08G1120 (GH3.5) exhibited a significant upregulation across the two stress factors. Functional characterization of Gh_A08G1120 (GH3.5) through virus-induced gene silencing (VIGS) revealed that the VIGS plants ability to tolerate drought and salt stresses was significantly reduced compared to the wild types. The chlorophyll content, relative leaf water content (RLWC), and superoxide dismutase (SOD) concentration level were reduced significantly while malondialdehyde concentration and ion leakage as a measure of cell membrane stability (CMS) increased in VIGS plants under drought and salt stress conditions.

CONCLUSION

This study revealed the significance of the GH3 genes in enabling the plant's adaptation to drought and salt stress conditions as evidenced by the VIGS results and RT-qPCR analysis.

Genome sequencing of Prototheca zopfii genotypes 1 and 2 provides evidence of a severe reduction in organellar genomes

Prototheca zopfii (P. zopfii, class Trebouxiophyceae, order Chlorellales, family Chlorellaceae), a non-photosynthetic predominantly free-living unicellular alga, is one of the few pathogens belonging to the plant kingdom. This alga can affect many vertebrate hosts, sustaining systemic infections and diseases such as mastitis in cows. The aim of our work was to sequence and assemble the P. zopfii genotype 1 and genotype 2 mitochondrial and plastid genomes. Remarkably, the P. zopfii mitochondrial (38 Kb) and plastid (28 Kb) genomes are models of compaction and the smallest known in the Trebouxiophyceae. As expected, the P. zopfii genotype 1 and 2 plastid genomes lack all the genes involved in photosynthesis, but, surprisingly, they also lack those coding for RNA polymerases. Our results showed that plastid genes are actively transcribed in P. zopfii, which suggests that the missing RNA polymerases are substituted by nuclear-encoded paralogs. The simplified architecture and highly-reduced gene complement of the P. zopfii mitochondrial and plastid genomes are closer to those of P. stagnora and the achlorophyllous obligate parasite Helicosporidium than to those of P. wickerhamii or P. cutis. This similarity is also supported by maximum likelihood phylogenetic analyses inferences. Overall, the P. zopfii sequences reported here, which include nuclear genome drafts for both genotypes, will help provide both a deeper understanding of the evolution of Prototheca spp. and insights into the corresponding host/pathogen interactions.

Alcohol Acetyltransferase Eat1 Is Located in Yeast Mitochondria

Eat1 is a recently discovered alcohol acetyltransferase responsible for bulk ethyl acetate production in yeasts such as Wickerhamomyces anomalus and Kluyveromyces lactis These yeasts have the potential to become efficient bio-based ethyl acetate producers. However, some fundamental features of Eat1 are still not understood, which hampers the rational engineering of efficient production strains. The cellular location of Eat1 in yeast is one of these features. To reveal its location, Eat1 was fused with yeast-enhanced green fluorescent protein (yEGFP) to allow intracellular tracking. Despite the current assumption that bulk ethyl acetate production occurs in the yeast cytosol, most of Eat1 localized to the mitochondria of Kluyveromyces lactis CBS 2359 Δku80 We then compared five bulk ethyl acetate-producing yeasts in iron-limited chemostats with glucose as the carbon source. All yeasts produced ethyl acetate under these conditions. This strongly suggests that the mechanism and location of bulk ethyl acetate synthesis are similar in these yeast strains. Furthermore, an in silico analysis showed that Eat1 proteins from various yeasts were mostly predicted as mitochondrial. Altogether, it is concluded that Eat1-catalyzed ethyl acetate production occurs in yeast mitochondria. This study has added new insights into bulk ethyl acetate synthesis in yeast, which is relevant for developing efficient production strains.IMPORTANCE Ethyl acetate is a common bulk chemical that is currently produced from petrochemical sources. Several Eat1-containing yeast strains naturally produce large amounts of ethyl acetate and are potential cell factories for the production of bio-based ethyl acetate. Rational design of the underlying metabolic pathways may result in improved production strains, but it requires fundamental knowledge on the function of Eat1. A key feature is the location of Eat1 in the yeast cell. The precursors for ethyl acetate synthesis can be produced in multiple cellular compartments through different metabolic pathways. The location of Eat1 determines the relevance of each pathway, which will provide future targets for the metabolic engineering of bulk ethyl acetate production in yeast.

Cotton Late Embryogenesis Abundant (LEA2) Genes Promote Root Growth and Confer Drought Stress Tolerance in Transgenic Arabidopsis thaliana

Late embryogenesis abundant (LEA) proteins play key roles in plant drought tolerance. In this study, 157, 85 and 89 candidate LEA2 proteins were identified in G. hirsutum, G. arboreum and G. raimondii respectively. LEA2 genes were classified into 6 groups, designated as group 1 to 6. Phylogenetic tree analysis revealed orthologous gene pairs within the cotton genome. The cotton specific LEA2 motifs identified were E, R and D in addition to Y, K and S motifs. The genes were distributed on all chromosomes. LEA2s were found to be highly enriched in non-polar, aliphatic amino acid residues, with leucine being the highest, 9.1% in proportion. The miRNA, ghr-miR827a/b/c/d and ghr-miR164 targeted many genes are known to be drought stress responsive. Various stress-responsive regulatory elements, ABA-responsive element (ABRE), Drought-responsive Element (DRE/CRT), MYBS and low-temperature-responsive element (LTRE) were detected. Most genes were highly expressed in leaves and roots, being the primary organs greatly affected by water deficit. The expression levels were much higher in G. tomentosum as opposed to G. hirsutum. The tolerant genotype had higher capacity to induce more of LEA2 genes. Over expression of the transformed gene Cot_AD24498 showed that the LEA2 genes are involved in promoting root growth and in turn confers drought stress tolerance. We therefore infer that Cot_AD24498, CotAD_20020, CotAD_21924 and CotAD_59405 could be the candidate genes with profound functions under drought stress in upland cotton among the LEA2 genes. The transformed Arabidopsis plants showed higher tolerance levels to drought stress compared to the wild types. There was significant increase in antioxidants, catalase (CAT), peroxidase (POD) and superoxide dismutase (SOD) accumulation, increased root length and significant reduction in oxidants, Hydrogen peroxide (H₂O₂) and malondialdehyde (MDA) concentrations in the leaves of transformed lines under drought stress condition. This study provides comprehensive analysis of LEA2 proteins in cotton thus forms primary foundation for breeders to utilize these genes in developing drought tolerant genotypes.

Genome-Wide Analysis of Multidrug and Toxic Compound Extrusion (MATE) Family in Gossypium raimondii and Gossypium arboreum and Its Expression Analysis Under Salt, Cadmium, and Drought Stress

The extrusion of toxins and substances at a cellular level is a vital life process in plants under abiotic stress. The multidrug and toxic compound extrusion (MATE) gene family plays a large role in the exportation of toxins and other substrates. We carried out a genome-wide analysis of MATE gene families in Gossypium raimondii and Gossypium arboreum and assessed their expression levels under salt, cadmium and drought stresses. We identified 70 and 68 MATE genes in G. raimondii and G. arboreum, respectively. The majority of the genes were predicted to be localized within the plasma membrane, with some distributed in other cell parts. Based on phylogenetic analysis, the genes were subdivided into three subfamilies, designated as M1, M2 and M3. Closely related members shared similar gene structures, and thus were highly conserved in nature and have mainly evolved through purifying selection. The genes were distributed in all chromosomes. Twenty-nine gene duplication events were detected, with segmental being the dominant type. GO annotation revealed a link to salt, drought and cadmium stresses. The genes exhibited differential expression, with GrMATE18, GrMATE34, GaMATE41 and GaMATE51 significantly upregulated under drought, salt and cadmium stress, and these could possibly be the candidate genes. Our results provide the first data on the genome-wide and functional characterization of MATE genes in diploid cotton, and are important for breeders of more stress-tolerant cotton genotypes.

Analysis of tail-anchored protein translocation pathway in plants

Tail-anchored (TA) proteins are a special class of membrane proteins that carry out vital functions in all living cells. Targeting mechanisms of TA proteins are investigated as the best example for post-translational protein targeting in yeast. Of the several mechanisms, Guided Entry of Tail-anchored protein (GET) pathway plays a major role in TA protein targeting. Many in silico and in vivo analyses are geared to identify TA proteins and their targeting mechanisms in different systems including Arabidopsis thaliana. Yet, crop plants that grow in specific and/or different conditions are not investigated for the presence of TA proteins and GET pathway. This study majorly investigates GET pathway in two crop plants, Oryza sativa subsp. Indica and Solanum tuberosum, through detailed in silico analysis. 508 and 912 TA proteins are identified in Oryza sativa subsp. Indica and Solanum tuberosum respectively and their localization with respect to endoplasmic reticulum (ER), mitochondria, and chloroplast has been delineated. Similarly, the associated GET proteins are identified (Get1, Get3 and Get4) and their structural inferences are elucidated using homology modelling. Get3 models are based on yeast Get3. The cytoplasmic Get3 from O. sativa is identified to be very similar to yeast Get3 with conserved P-loop and TA binding groove. Three cytoplasmic Get3s are identified for S. tuberosum. Taken together, this is the first study to identify TA proteins and GET components in Oryza sativa subsp. Indica and Solanum tuberosum, forming the basis for any further experimental characterization of TA targeting and GET pathway mechanisms in crop plants.

Decoding the Divergent Subcellular Location of Two Highly Similar Paralogous LEA Proteins

Many mitochondrial proteins are synthesized as precursors in the cytosol with an N-terminal mitochondrial targeting sequence (MTS) which is cleaved off upon import. Although much is known about import mechanisms and MTS structural features, the variability of MTS still hampers robust sub-cellular software predictions. Here, we took advantage of two paralogous late embryogenesis abundant proteins (LEA) from Arabidopsis with different subcellular locations to investigate structural determinants of mitochondrial import and gain insight into the evolution of the LEA genes. LEA38 and LEA2 are short proteins of the LEA_3 family, which are very similar along their whole sequence, but LEA38 is targeted to mitochondria while LEA2 is cytosolic. Differences in the N-terminal protein sequences were used to generate a series of mutated LEA2 which were expressed as GFP-fusion proteins in leaf protoplasts. By combining three types of mutation (substitution, charge inversion, and segment replacement), we were able to redirect the mutated LEA2 to mitochondria. Analysis of the effect of the mutations and determination of the LEA38 MTS cleavage site highlighted important structural features within and beyond the MTS. Overall, these results provide an explanation for the likely loss of mitochondrial location after duplication of the ancestral gene.

A Novel G-Protein-Coupled Receptors Gene from Upland Cotton Enhances Salt Stress Tolerance in Transgenic Arabidopsis

Plants have developed a number of survival strategies which are significant for enhancing their adaptation to various biotic and abiotic stress factors. At the transcriptome level, G-protein-coupled receptors (GPCRs) are of great significance, enabling the plants to detect a wide range of endogenous and exogenous signals which are employed by the plants in regulating various responses in development and adaptation. In this research work, we carried out genome-wide analysis of target of Myb1 (TOM1), a member of the GPCR gene family. The functional role of TOM1 in salt stress tolerance was studied using a transgenic Arabidopsis plants over-expressing the gene. By the use of the functional domain PF06454, we obtained 16 TOM genes members in Gossypium hirsutum, 9 in Gossypium arboreum, and 11 in Gossypium raimondii. The genes had varying physiochemical properties, and it is significant to note that all the grand average of hydropathy (GRAVY) values were less than one, indicating that all are hydrophobic in nature. In all the genes analysed here, both the exonic and intronic regions were found. The expression level of Gh_A07G0747 (GhTOM) was significantly high in the transgenic lines as compared to the wild type; a similar trend in expression was observed in all the salt-related genes tested in this study. The study in epidermal cells confirmed the localization of the protein coded by the gene TOM1 in the plasma membrane. Analysis of anti-oxidant enzymes showed higher concentrations of antioxidants in transgenic lines and relatively lower levels of oxidant substances such as H₂O₂. The low malondialdehyde (MDA) level in transgenic lines indicated that the transgenic lines had relatively low level of oxidative damage compared to the wild types. The results obtained indicate that Gh_A07G0747 (GhTOM) can be a putative target gene for enhancing salt stress tolerance in plants and could be exploited in the future for the development of salt stress-tolerant cotton cultivars.

A Guide to Computational Methods for Predicting Mitochondrial Localization

Predicting mitochondrial localization of proteins remains challenging for two main reasons: (1) Not only one but several mitochondrial localization signals exist, which primarily dictate the final destination of a protein in this organelle. However, most localization prediction algorithms rely on the presence of a so-called presequence (or N-terminal mitochondrial targeting peptide, mTP), which occurs in only ~70% of mitochondrial proteins. (2) The presequence is highly divergent on sequence level and therefore difficult to identify on the computer.In this chapter, we review a number of protein localization prediction programs and propose a strategy to predict mitochondrial localization. Finally, we give some helpful suggestions for bench scientists when working with mitochondrial protein candidates in silico.

The cinnamyl alcohol dehydrogenase family in flax: Differentiation during plant growth and under stress conditions

Cinnamyl alcohol dehydrogenase (CAD), which catalyzes the reduction of cinnamaldehydes to their alcohol derivatives, is represented by a large family of proteins. The aim of the study was to identify the CAD isoforms in flax (Linum usitatissimum L.) - LuCADs - and to determine their specificity to enhance knowledge of the mechanisms controlling cell wall lignification in flax under environmental stresses. On the basis of genome-wide analysis, we identified 15 isoforms (one in two copies) belonging to three major classes of the CAD protein family. Their specificity was determined at the transcriptomic level in different tissues/organs, under Fusarium infection and abiotic stresses. Considering the function of particular LuCADs, it was established that LuCAD1 and 2 belong to Class I and they take part in the lignification of maturing stem and in the response to cold and drought stress. The Class II members LuCAD3, LuCAD4, LuCAD5 and LuCAD6 play various roles in flax being putatively responsible for lignin synthesis in different organs or under certain conditions. The obtained results indicate that within Class II, LuCAD6 was the most abundant in seedlings and maturing stems, LuCAD3 in leaves, and LuCAD4 in stems. Comparative analysis showed that expression of LuCAD genes in roots after F. oxysporum infection had the greatest contribution to differentiation of LuCAD expression patterns. Surprisingly, most of the analyzed LuCAD isoforms had reduced expression after pathogen infection. The decrease in mRNA level was primarily observed for LuCAD6 and LuCAD4, but also LuCAD1 and 8. However, the induction of LuCAD expression was mostly characteristic for Class I LuCAD1 and 2 in leaves. For cold stress, a clear correlation with phylogenic class membership was observed. Low temperatures caused induction of CAD isoforms belonging to Class I and repression of LuCADs from Class III.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

Background

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

Results

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

Conclusion

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0596-1) contains supplementary material, which is available to authorized users.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

BACKGROUND

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

RESULTS

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

CONCLUSION

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

The Obvious and the Hidden: Prediction and Function of Fungal Peroxisomal Matrix Proteins

Fungal peroxisomes are characterized by a number of specific biological functions. To understand the physiology and biochemistry of these organelles knowledge of the proteome content is crucial. Here, we address different strategies to predict peroxisomal proteins by bioinformatics approaches. These tools range from simple text searches to network based learning strategies. A complication of this analysis is the existence of cryptic peroxisomal proteins, which are overlooked in conventional bioinformatics queries. These include proteins where targeting information results from transcriptional and posttranscriptional alterations. But also proteins with low efficiency targeting motifs that are predominantly localized in the cytosol, and proteins lacking any canonical targeting information, can play important roles within peroxisomes. Many of these proteins are so far unpredictable. Detection and characterization of these cryptic peroxisomal proteins revealed the presence of novel peroxisomal enzymatic reaction networks in fungi.

Prediction of Peroxisomal Matrix Proteins in Plants

Our knowledge of the proteome of plant peroxisomes is far from being complete, and the functional complexity and plasticity of this cell organelle are amazingly high particularly in plants, as exemplified by the model species Arabidopsis thaliana. Plant-specific peroxisome functions that have been uncovered only recently include, for instance, the participation of peroxisomes in phylloquinone and biotin biosynthesis. Experimental proteome studies have been proved very successful in defining the proteome of Arabidopsis peroxisomes but this approach also faces significant challenges and limitations. Complementary to experimental approaches, computational methods have emerged as important powerful tools to define the proteome of soluble matrix proteins of plant peroxisomes. Compared to other cell organelles such as mitochondria, plastids and the ER, the simultaneous operation of two major import pathways for soluble proteins in peroxisomes is rather atypical. Novel machine learning prediction approaches have been developed for peroxisome targeting signals type 1 (PTS1) and revealed high sensitivity and specificity, as validated by in vivo subcellular targeting analyses in diverse transient plant expression systems. Accordingly, the algorithms allow the correct prediction of many novel peroxisome-targeted proteins from plant genome sequences and the discovery of additional organelle functions. In contrast, the prediction of PTS2 proteins largely remains restricted to genome searches by conserved patterns contrary to more advanced machine learning methods. Here, we summarize and discuss the capabilities and accuracies of available prediction algorithms for PTS1 and PTS2 carrying proteins.

Predicting Peroxisomal Targeting Signals to Elucidate the Peroxisomal Proteome of Mammals

Peroxisomes harbor a plethora of proteins, but the peroxisomal proteome as the entirety of all peroxisomal proteins is still unknown for mammalian species. Computational algorithms can be used to predict the subcellular localization of proteins based on their amino acid sequence and this method has been amply used to forecast the intracellular fate of individual proteins. However, when applying such algorithms systematically to all proteins of an organism the prediction of its peroxisomal proteome in silico should be possible. Therefore, a reliable detection of peroxisomal targeting signals (PTS ) acting as postal codes for the intracellular distribution of the encoding protein is crucial. Peroxisomal proteins can utilize different routes to reach their destination depending on the type of PTS. Accordingly, independent prediction algorithms have been developed for each type of PTS, but only those for type-1 motifs (PTS1) have so far reached a satisfying predictive performance. This is partially due to the low number of peroxisomal proteins limiting the power of statistical analyses and partially due to specific properties of peroxisomal protein import, which render functional PTS motifs inactive in specific contexts. Moreover, the prediction of the peroxisomal proteome is limited by the high number of proteins encoded in mammalian genomes, which causes numerous false positive predictions even when using reliable algorithms and buries the few yet unidentified peroxisomal proteins. Thus, the application of prediction algorithms to identify all peroxisomal proteins is currently ineffective as stand-alone method, but can display its full potential when combined with other methods.

Multiple Localization by Functional Translational Readthrough

In a compartmentalized cell, correct protein localization is crucial for function of virtually all cellular processes. From the cytoplasm as a starting point, proteins are imported into organelles by specific targeting signals. Many proteins, however, act in more than one cellular compartment. In this chapter, we discuss mechanisms by which proteins can be targeted to multiple organelles with a focus on a novel gene regulatory mechanism, functional translational readthrough, that permits multiple targeting of proteins to the peroxisome and other organelles. In mammals, lactate and malate dehydrogenase are the best-characterized enzymes whose targeting is controlled by functional translational readthrough.

Subcellular Localization Screening of Colletotrichum higginsianum Effector Candidates Identifies Fungal Proteins Targeted to Plant Peroxisomes, Golgi Bodies, and Microtubules

The genome of the hemibiotrophic anthracnose fungus, Colletotrichum higginsianum, encodes a large inventory of putative secreted effector proteins that are sequentially expressed at different stages of plant infection, namely appressorium-mediated penetration, biotrophy and necrotrophy. However, the destinations to which these proteins are addressed inside plant cells are unknown. In the present study, we selected 61 putative effector genes that are highly induced in appressoria and/or biotrophic hyphae. We then used Agrobacterium-mediated transformation to transiently express them as N-terminal fusions with fluorescent proteins in cells of Nicotiana benthamiana for imaging by confocal microscopy. Plant compartments labeled by the fusion proteins in N. benthamiana were validated by co-localization with specific organelle markers, by transient expression of the proteins in the true host plant, Arabidopsis thaliana, and by transmission electron microscopy-immunogold labeling. Among those proteins for which specific subcellular localizations could be verified, nine were imported into plant nuclei, three were imported into the matrix of peroxisomes, three decorated cortical microtubule arrays and one labeled Golgi stacks. Two peroxisome-targeted proteins harbored canonical C-terminal tripeptide signals for peroxisome import via the PTS1 (peroxisomal targeting signal 1) pathway, and we showed that these signals are essential for their peroxisome localization. Our findings provide valuable information about which host processes are potentially manipulated by this pathogen, and also reveal plant peroxisomes, microtubules, and Golgi as novel targets for fungal effectors.

The evolutionary fate of the chloroplast and nuclear rps16 genes as revealed through the sequencing and comparative analyses of four novel legume chloroplast genomes from Lupinus

The Fabaceae family is considered as a model system for understanding chloroplast genome evolution due to the presence of extensive structural rearrangements, gene losses and localized hypermutable regions. Here, we provide sequences of four chloroplast genomes from the Lupinus genus, belonging to the underinvestigated Genistoid clade. Notably, we found in Lupinus species the functional loss of the essential rps16 gene, which was most likely replaced by the nuclear rps16 gene that encodes chloroplast and mitochondrion targeted RPS16 proteins. To study the evolutionary fate of the rps16 gene, we explored all available plant chloroplast, mitochondrial and nuclear genomes. Whereas no plant mitochondrial genomes carry an rps16 gene, many plants still have a functional nuclear and chloroplast rps16 gene. Ka/Ks ratios revealed that both chloroplast and nuclear rps16 copies were under purifying selection. However, due to the dual targeting of the nuclear rps16 gene product and the absence of a mitochondrial copy, the chloroplast gene may be lost. We also performed comparative analyses of lupine plastomes (SNPs, indels and repeat elements), identified the most variable regions and examined their phylogenetic utility. The markers identified here will help to reveal the evolutionary history of lupines, Genistoids and closely related clades.

Ancient signal for nitrogen status sensing in the green lineage: Functional evidence of CDPK repertoire in Ostreococcus tauri

Calcium-dependent protein kinases (CDPKs) regulate plant development and many stress signalling pathways through the complex cytosolic [Ca²⁺] signalling. The genome of Ostreococcus tauri (Ot), a model prasinophyte organism that is on the base of the green lineage, harbours three sequences homologous to those encoding plant CDPKs with the three characteristic conserved domains (protein kinase, autoregulatory/autoinhibitory, and regulatory domain). Phylogenetic and structural analyses revealed that putative OtCDPK proteins are closely related to CDPKs from other Chlorophytes. We functionally characterised the first marine picophytoeukaryote CDPK gene (OtCDPK1) and showed that the expression of the three OtCDPK genes is up-regulated by nitrogen depletion. We conclude that CDPK signalling pathway might have originated early in the green lineage and may play a key role in prasinophytes by sensing macronutrient changes in the marine environment.

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.

Protein Function Prediction

Protein function is a concept that can have different interpretations in different biological contexts, and the number and diversity of novel proteins identified by large-scale "omics" technologies poses increasingly new challenges. In this review we explore current strategies used to predict protein function focused on high-throughput sequence analysis, as for example, inference based on sequence similarity, sequence composition, structure, and protein-protein interaction. Various prediction strategies are discussed together with illustrative workflows highlighting the use of some benchmark tools and knowledge bases in the field.

Genome-wide characterisation and expression profile of the grapevine ATL ubiquitin ligase family reveal biotic and abiotic stress-responsive and development-related members

The Arabidopsis Tóxicos en Levadura (ATL) protein family is a class of E3 ubiquitin ligases with a characteristic RING-H2 Zn-finger structure that mediates diverse physiological processes and stress responses in plants. We carried out a genome-wide survey of grapevine (Vitis vinifera L.) ATL genes and retrieved 96 sequences containing the canonical ATL RING-H2 domain. We analysed their genomic organisation, gene structure and evolution, protein domains and phylogenetic relationships. Clustering revealed several clades, as already reported in Arabidopsis thaliana and rice (Oryza sativa), with an expanded subgroup of grapevine-specific genes. Most of the grapevine ATL genes lacked introns and were scattered among the 19 chromosomes, with a high level of duplication retention. Expression profiling revealed that some ATL genes are expressed specifically during early or late development and may participate in the juvenile to mature plant transition, whereas others may play a role in pathogen and/or abiotic stress responses, making them key candidates for further functional analysis. Our data offer the first genome-wide overview and annotation of the grapevine ATL family, and provide a basis for investigating the roles of specific family members in grapevine physiology and stress responses, as well as potential biotechnological applications.

Identification of New Fungal Peroxisomal Matrix Proteins and Revision of the PTS1 Consensus

The peroxisomal targeting signal type 1 (PTS1) is a seemingly simple peptide sequence at the C-terminal end of most peroxisomal matrix proteins. PTS1 can be described as a tripeptide with the consensus motif [S/A/C] [K/R/H] L. However, this description is neither necessary nor sufficient. It does not cover all cases of PTS1 proteins, and some proteins in accordance with this consensus do not target to the peroxisome. In order to find new PTS proteins in yeast and to arrive at a more complete description of the PTS1 consensus motif, we developed a machine learning approach that involves orthologue expansion of the set of known peroxisomal proteins. We performed a genome-wide in silico screen, characterised several PTS1-containing peptides and identified two new peroxisomal matrix proteins, which we named Pxp1 (Yel020c) and Pxp2 (Yjr111c). Based on these in silico and in vivo analyses, we revised the yeast PTS1 consensus which now includes all known PTS1 proteins.

Effects of lng Mutations on LngA Expression, Processing, and CS21 Assembly in Enterotoxigenic Escherichia coli E9034A

Enterotoxigenic Escherichia coli (ETEC) is a major cause of morbidity in children under 5 years of age in low- and middle-income countries and a leading cause of traveler's diarrhea worldwide. The ability of ETEC to colonize the intestinal epithelium is mediated by fimbrial adhesins, such as CS21 (Longus). This adhesin is a type IVb pilus involved in adherence to intestinal cells in vitro and bacterial self-aggregation. Fourteen open reading frames have been proposed to be involved in CS21 assembly, hitherto only the lngA and lngB genes, coding for the major (LngA) and minor (LngB) structural subunit, have been characterized. In this study, we investigated the role of the LngA, LngB, LngC, LngD, LngH, and LngP proteins in the assembly of CS21 in ETEC strain E9034A. The deletion of the lngA, lngB, lngC, lngD, lngH, or lngP genes, abolished CS21 assembly in ETEC strain E9034A and the adherence to HT-29 cells was reduced 90%, compared to wild-type strain. Subcellular localization prediction of CS21 proteins was similar to other well-known type IV pili homologs. We showed that LngP is the prepilin peptidase of LngA, and that ETEC strain E9034A has another peptidase capable of processing LngA, although with less efficiency. Additionally, we present immuno-electron microscopy images to show that the LngB protein could be localized at the tip of CS21. In conclusion, our results demonstrate that the LngA, LngB, LngC, LngD, LngH, and LngP proteins are essential for CS21 assembly, as well as for bacterial aggregation and adherence to HT-29 cells.

Phosphorylated and Nonphosphorylated PfMAP2 Are Localized in the Nucleus, Dependent on the Stage of Plasmodium falciparum Asexual Maturation

Plasmodium falciparum mitogen-activated protein (MAP) kinases, a family of enzymes central to signal transduction processes including inflammatory responses, are a promising target for antimalarial drug development. Our study shows for the first time that the P. falciparum specific MAP kinase 2 (PfMAP2) is colocalized in the nucleus of all of the asexual erythrocytic stages of P. falciparum and is particularly elevated in its phosphorylated form. It was also discovered that PfMAP2 is expressed in its highest quantity during the early trophozoite (ring form) stage and significantly reduced in the mature trophozoite and schizont stages. Although the phosphorylated form of the kinase is always more prevalent, its ratio relative to the nonphosphorylated form remained constant irrespective of the parasites' developmental stage. We have also shown that the TSH motif specifically renders PfMAP2 genetically divergent from the other plasmodial MAP kinase activation sites using Neighbour Joining analysis. Furthermore, TSH motif-specific designed antibody is crucial in determining the location of the expression of the PfMAP2 protein. However, by using immunoelectron microscopy, PPfMAP2 were detected ubiquitously in the parasitized erythrocytes. In summary, PfMAP2 may play a far more important role than previously thought and is a worthy candidate for research as an antimalarial.

Characterization, prediction and evolution of plant peroxisomal targeting signals type 1 (PTS1s)

Our knowledge of the proteome of plant peroxisomes and their functional plasticity is far from being complete, primarily due to major technical challenges in experimental proteome research of the fragile cell organelle. Several unexpected novel plant peroxisome functions, for instance in biotin and phylloquinone biosynthesis, have been uncovered recently. Nevertheless, very few regulatory and membrane proteins of plant peroxisomes have been identified and functionally described up to now. To define the matrix proteome of plant peroxisomes, computational methods have emerged as important powerful tools. Novel prediction approaches of high sensitivity and specificity have been developed for peroxisome targeting signals type 1 (PTS1) and have been validated by in vivo subcellular targeting analyses and thermodynamic binding studies with the cytosolic receptor, PEX5. Accordingly, the algorithms allow the correct prediction of many novel peroxisome-targeted proteins from plant genome sequences and the discovery of additional organelle functions. In this review, we provide an overview of methodologies, capabilities and accuracies of available prediction algorithms for PTS1 carrying proteins. We also summarize and discuss recent quantitative, structural and mechanistic information of the interaction of PEX5 with PTS1 carrying proteins in relation to in vivo import efficiency. With this knowledge, we develop a model of how proteins likely evolved peroxisomal targeting signals in the past and still nowadays, in which order the two import pathways might have evolved in the ancient eukaryotic cell, and how the secondary loss of the PTS2 pathway probably happened in specific organismal groups.

Genome-wide identification, structural analysis and new insights into late embryogenesis abundant (LEA) gene family formation pattern in Brassica napus

Late embryogenesis abundant (LEA) proteins are a diverse and large group of polypeptides that play important roles in desiccation and freezing tolerance in plants. The LEA family has been systematically characterized in some plants but not Brassica napus. In this study, 108 BnLEA genes were identified in the B. napus genome and classified into eight families based on their conserved domains. Protein sequence alignments revealed an abundance of alanine, lysine and glutamic acid residues in BnLEA proteins. The BnLEA gene structure has few introns (<3), and they are distributed unevenly across all 19 chromosomes in B. napus, occurring as gene clusters in chromosomes A9, C2, C4 and C5. More than two-thirds of the BnLEA genes are associated with segmental duplication. Synteny analysis revealed that most LEA genes are conserved, although gene losses or gains were also identified. These results suggest that segmental duplication and whole-genome duplication played a major role in the expansion of the BnLEA gene family. Expression profiles analysis indicated that expression of most BnLEAs was increased in leaves and late stage seeds. This study presents a comprehensive overview of the LEA gene family in B. napus and provides new insights into the formation of this family.

Quantification of expression of dehydrin isoforms in the desiccation tolerant plant Craterostigma plantagineum using specifically designed reference genes

Craterostigma plantagineum is a desiccation tolerant resurrection plant. Many genes are induced during desiccation. Dehydrins are a group of dehydration-induced genes present in all higher plants. The current study aims at classifying the most abundantly expressed dehydrin genes from vegetative tissues of C. plantagineum and quantifying their expression. To identify variations between dehydrin isoforms at different stages of desiccation and rehydration by RT-qPCR, the target mRNA requires an accurate and reliable normalization. Previously we reported that RNAs from leaves and roots of C. plantagineum are not degraded during desiccation and subsequent rehydration thus allowing the use of RT-qPCR to test the stability of reference genes. The expression stability of eight candidate reference genes was tested in leaves, roots and callus. These genes were ranked according to their stability of gene expression using GeNorm(PLUS) and RefFinder. The most consistently expressed reference genes in each tissue were identified and used to normalize gene expression data. Dehydrin isoforms were divided in three groups based on the expression level during the desiccation process in three different tissues (leaves, roots and callus).

Mutations in the plant-conserved MTERF9 alter chloroplast gene expression, development and tolerance to abiotic stress in Arabidopsis thaliana

The control of organelle gene expression in plants is far from fully understood. The characterization of mutants in Arabidopsis thaliana is assigning an increasingly prominent role to the mitochondrial transcription termination factors (mTERFs) in this process. To gain insight into the function of mTERF genes in plants, we took a reverse genetics approach to identify and characterize A. thaliana mTERF-defective mutants. Here we report the characterization of the mterf9 mutant, affected in an mTERF protein functionally conserved in plants and targeted to chloroplasts. Loss of MTERF9 results in defective chloroplast development, which is likely to cause paleness, stunted growth and reduced mesophyll cell numbers. Expression analysis of different plastid genes revealed reduced levels of plastid-encoded polymerase (PEP)-dependent transcripts and increased levels of transcripts dependent of nucleus-encoded polymerase. mterf9 plants exhibited altered responses to sugars, abscisic acid (ABA), salt and osmotic stresses, and the microarray data analysis showed modifications in MTERF9 expression after salt or mannitol treatments. Our genetic interactions results indicate a functional relationship between MTERF9 and the previously characterized MDA1 gene, and between MDA1 and some plastid ribosomal genes. MDA1 and MTERF9 were upregulated in the mterf9 and mda1 mutants, respectively. Moreover, 21 of 50 genes were commonly co-expressed with MDA1 and MTERF9. The analysis of the MDA1 and MTERF9 promoters showed that both were rich in stress-related cis-regulatory elements. Our results highlight the role of the MTERF9 gene in plant biology and deepens the understanding of the functional relationship of plant mTERF genes.

MPIC: a mitochondrial protein import components database for plant and non-plant species

In the 2 billion years since the endosymbiotic event that gave rise to mitochondria, variations in mitochondrial protein import have evolved across different species. With the genomes of an increasing number of plant species sequenced, it is possible to gain novel insights into mitochondrial protein import pathways. We have generated the Mitochondrial Protein Import Components (MPIC) Database (DB; http://www.plantenergy.uwa.edu.au/applications/mpic) providing searchable information on the protein import apparatus of plant and non-plant mitochondria. An in silico analysis was carried out, comparing the mitochondrial protein import apparatus from 24 species representing various lineages from Saccharomyces cerevisiae (yeast) and algae to Homo sapiens (human) and higher plants, including Arabidopsis thaliana (Arabidopsis), Oryza sativa (rice) and other more recently sequenced plant species. Each of these species was extensively searched and manually assembled for analysis in the MPIC DB. The database presents an interactive diagram in a user-friendly manner, allowing users to select their import component of interest. The MPIC DB presents an extensive resource facilitating detailed investigation of the mitochondrial protein import machinery and allowing patterns of conservation and divergence to be recognized that would otherwise have been missed. To demonstrate the usefulness of the MPIC DB, we present a comparative analysis of the mitochondrial protein import machinery in plants and non-plant species, revealing plant-specific features that have evolved.

Bioinformatic indications that COPI- and clathrin-based transport systems are not present in chloroplasts: an Arabidopsis model

Coated vesicle transport occurs in the cytosol of yeast, mammals and plants. It consists of three different transport systems, the COPI, COPII and clathrin coated vesicles (CCV), all of which participate in the transfer of proteins and lipids between different cytosolic compartments. There are also indications that chloroplasts have a vesicle transport system. Several putative chloroplast-localized proteins, including CPSAR1 and CPRabA5e with similarities to cytosolic COPII transport-related proteins, were detected in previous experimental and bioinformatics studies. These indications raised the hypothesis that a COPI- and/or CCV-related system may be present in chloroplasts, in addition to a COPII-related system. To test this hypothesis we bioinformatically searched for chloroplast proteins that may have similar functions to known cytosolic COPI and CCV components in the model plants Arabidopsis thaliana and Oryza sativa (subsp. japonica) (rice). We found 29 such proteins, based on domain similarity, in Arabidopsis, and 14 in rice. However, many components could not be identified and among the identified most have assigned roles that are not related to either COPI or CCV transport. We conclude that COPII is probably the only active vesicle system in chloroplasts, at least in the model plants. The evolutionary implications of the findings are discussed.

An efficient method for stable protein targeting in grasses (Poaceae): a case study in Puccinellia tenuiflora

BACKGROUND

An efficient transformation method is lacking for most non-model plant species to test gene function. Therefore, subcellular localization of proteins of interest from non-model plants is mainly carried out through transient transformation in homologous cells or in heterologous cells from model species such as Arabidopsis. Although analysis of expression patterns in model organisms like yeast and Arabidopsis can provide important clues about protein localization, these heterologous systems may not always faithfully reflect the native subcellular distribution in other species. On the other hand, transient expression in protoplasts from species of interest has limited ability for detailed sub-cellular localization analysis (e.g., those involving subcellular fractionation or sectioning and immunodetection), as it results in heterogeneous populations comprised of both transformed and untransformed cells.

RESULTS

We have developed a simple and reliable method for stable transformation of plant cell suspensions that are suitable for protein subcellular localization analyses in the non-model monocotyledonous plant Puccinellia tenuiflora. Optimization of protocols for obtaining suspension-cultured cells followed by Agrobacterium-mediated genetic transformation allowed us to establish stably transformed cell lines, which could be maintained indefinitely in axenic culture supplied with the proper antibiotic. As a case study, protoplasts of transgenic cell lines stably transformed with an ammonium transporter-green fluorescent protein (PutAMT1;1-GFP) fusion were successfully used for subcellular localization analyses in P. tenuiflora.

CONCLUSIONS

We present a reliable method for the generation of stably transformed P. tenuiflora cell lines, which, being available in virtually unlimited amounts, can be conveniently used for any type of protein subcellular localization analysis required. Given its simplicity, the method can be used as reference for other non-model plant species lacking efficient regeneration protocols.

Evolution of bacterial-like phosphoprotein phosphatases in photosynthetic eukaryotes features ancestral mitochondrial or archaeal origin and possible lateral gene transfer

Protein phosphorylation is a reversible regulatory process catalyzed by the opposing reactions of protein kinases and phosphatases, which are central to the proper functioning of the cell. Dysfunction of members in either the protein kinase or phosphatase family can have wide-ranging deleterious effects in both metazoans and plants alike. Previously, three bacterial-like phosphoprotein phosphatase classes were uncovered in eukaryotes and named according to the bacterial sequences with which they have the greatest similarity: Shewanella-like (SLP), Rhizobiales-like (RLPH), and ApaH-like (ALPH) phosphatases. Utilizing the wealth of data resulting from recently sequenced complete eukaryotic genomes, we conducted database searching by hidden Markov models, multiple sequence alignment, and phylogenetic tree inference with Bayesian and maximum likelihood methods to elucidate the pattern of evolution of eukaryotic bacterial-like phosphoprotein phosphatase sequences, which are predominantly distributed in photosynthetic eukaryotes. We uncovered a pattern of ancestral mitochondrial (SLP and RLPH) or archaeal (ALPH) gene entry into eukaryotes, supplemented by possible instances of lateral gene transfer between bacteria and eukaryotes. In addition to the previously known green algal and plant SLP1 and SLP2 protein forms, a more ancestral third form (SLP3) was found in green algae. Data from in silico subcellular localization predictions revealed class-specific differences in plants likely to result in distinct functions, and for SLP sequences, distinctive and possibly functionally significant differences between plants and nonphotosynthetic eukaryotes. Conserved carboxyl-terminal sequence motifs with class-specific patterns of residue substitutions, most prominent in photosynthetic organisms, raise the possibility of complex interactions with regulatory proteins.

Green targeting predictor and ambiguous targeting predictor 2: the pitfalls of plant protein targeting prediction and of transient protein expression in heterologous systems

The challenges of plant protein targeting prediction are the existence of dual subcellular targets and the bias of experimentally confirmed data towards few and mostly nonplant model species. To assess whether training with proteins from evolutionarily distant species has a negative impact on prediction accuracy, we developed the Green Targeting Predictor tool, which was trained with a species-specific data set for Physcomitrella patens. Its performance was compared with that of the same tool trained with a mixed data set. In addition, we updated the Ambiguous Targeting Predictor. We found that predictions deviated from in vivo observations predominantly for proteins diverging within the green lineage, as well as for dual targeted proteins. To evaluate the usefulness of heterologous expression systems, selected proteins were subjected to localization studies in P. patens, Arabidopsis thaliana and Nicotiana tabacum. Four out of six proteins that show dual targeting in the original plant system were located only in a single compartment in one or both heterologous systems. We conclude that targeting signals of divergent plant species exhibit differences, calling for custom in silico and in vivo approaches when aiming to unravel the actual distribution patterns of proteins within a plant cell.

Is ABCC6 a genuine mitochondrial protein?

Background

A recent article in Circulation Research suggests that the protein ABCC6, which when defective is responsible for pseudoxanthoma elasticum, an inherited condition with skin, eye and cardiovascular manifestations, is associated with dysfunction in mitochondria – Martin et al.: ABCC6 Localizes to the Mitochondria-Associated Membrane.Circ Res 2012, 111:516–520. We present complementary information based on a bioinformatics analysis, which was not performed in the article cited, to examine the suggestion that ABCC6 is localized to mitochondria.

Results

All the computational strategies and integrative approaches that constitute references in the field indicate that ABCC6 is localized outside of mitochondria.

Conclusion

Our computational and integrative results, including both experimental and predictive data, show that there is no support in favor of the localization of ABCC6 in mitochondria.