Plant protein-coding gene families: emerging bioinformatics approaches

Transcriptome-Based Identification of the SaR2R3-MYB Gene Family in Sophora alopecuroides and Function Analysis of SaR2R3-MYB15 in Salt Stress Tolerance

As one of the most prominent gene families, R2R3-MYB transcription factors significantly regulate biochemical and physiological processes under salt stress. However, in Sophora alopecuroides, a perennial herb known for its exceptional saline alkali resistance, the comprehensive identification and characterization of SaR2R3-MYB genes and their potential functions in response to salt stress have yet to be determined. We investigated the expression profiles and biological functions of SaR2R3-MYB transcription factors in response to salt stress, utilizing a transcriptome-wide mining method. Our analysis identified 28 SaR2R3-MYB transcription factors, all sharing a highly conserved R2R3 domain, which were further divided into 28 subgroups through phylogenetic analysis. Some SaR2R3-MYB transcription factors showed induction under salt stress, with SaR2R3-MYB15 emerging as a potential regulator based on analysis of the protein-protein interaction network. Validation revealed the transcriptional activity and nuclear localization of SaR2R3-MYB15. Remarkably, overexpression of SaR2R3-MYB15 in transgenic plants could increase the activity of antioxidant enzymes and the accumulation of proline but decrease the content of malondialdehyde (MDA), compared with wild-type plants. Moreover, several salt stress-related genes showed higher expression levels in transgenic plants, implying their potential to enhance salt tolerance. Our findings shed light on the role of SaR2R3-MYB genes in salt tolerance in S. alopecuroides.

dEMBF v2.0: An Updated Database of Enzymes for Microalgal Biofuel Feedstock

In light of increasing algal genomics data and knowledge of biosynthetic pathways responsible for biofuel production, an integrated resource for easy access to all information is essential to improve our understanding of algal lipid metabolism. Against this backdrop, dEMBF v2.0, a significantly updated and improved version of our database of microalgae lipid biosynthetic enzymes for biofuel production, has been developed. dEMBF v2.0 now contains a comprehensive annotation of 2018 sequences encoding 35 enzymes, an increase of over 7-fold as compared with the first version. Other improved features include an increase in species coverage to 32 algal genomes, analysis of additional metabolic pathways, expanded annotation thoroughly detailing sequence and structural features, including enzyme-ligand interactions, and integration of supporting experimental evidence to demonstrate the role of enzymes in increasing lipid content. Along with a complete redesign of the interface, the updated database provides several inbuilt tools and user-friendly functionalities for more interactive and dynamic visualization of data.

Multispecies genome-wide analysis defines the MAP3K gene family in Gossypium hirsutum and reveals conserved family expansions

Background

Gene families are sets of structurally and evolutionarily related genes – in one or multiple species – that typically share a conserved biological function. As such, the identification and subsequent analyses of entire gene families are widely employed in the fields of evolutionary and functional genomics of both well established and newly sequenced plant genomes. Currently, plant gene families are typically identified using one of two major ways: 1) HMM-profile based searches using models built on Arabidopsis thaliana genes or 2) coding sequence homology searches using curated databases. Integrated databases containing functionally annotated genes and gene families have been developed for model organisms and several important crops; however, a comprehensive methodology for gene family annotation is currently lacking, preventing automated annotation of newly sequenced genomes.

Results

This paper proposes a combined measure of homology identification, motif conservation, phylogenomic and integrated gene expression analyses to define gene family structures in multiple plant species. The MAP3K gene families in seven plant species, including two currently unexamined species Gossypium hirsutum, and Zostera marina, were characterized to reveal new insights into their collective function and evolution and demonstrate the effectiveness of our novel methodology.

Conclusion

Compared with recent reports, this methodology performs significantly better for the identification and analysis of gene family members in several monocots/dicots, diploid as well as polyploid plant species.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2624-9) contains supplementary material, which is available to authorized users.

Insights into the molecular evolution of peptidase inhibitors in arthropods

Peptidase inhibitors are key proteins involved in the control of peptidases. In arthropods, peptidase inhibitors modulate the activity of peptidases involved in endogenous physiological processes and peptidases of the organisms with which they interact. Exploring available arthropod genomic sequences is a powerful way to obtain the repertoire of peptidase inhibitors in every arthropod species and to understand the evolutionary mechanisms involved in the diversification of this kind of proteins. A genomic comparative analysis of peptidase inhibitors in species belonging to different arthropod taxonomic groups was performed. The results point out: i) species or clade-specific presence is shown for several families of peptidase inhibitors; ii) multidomain peptidase inhibitors are commonly found in many peptidase inhibitor families; iii) several families have a wide range of members in different arthropod species; iv) several peptidase inhibitor families show species-specific (or clade-specific) gene family expansions; v) functional divergence may be assumed for particular clades; vi) passive expansions may be used by natural selection to fix adaptations. In conclusion, conservation and divergence of duplicated genes and the potential recruitment as peptidase inhibitors of proteins from other families are the main mechanisms used by arthropods to fix diversity. This diversity would be associated to the control of target peptidases and, as consequence, to adapt to specific environments.

Expression analysis of the Arabidopsis thaliana AtSpen2 gene, and its relationship with other plant genes encoding Spen proteins

Proteins of the Split ends (Spen) family are characterized by an N-terminal domain, with one or more RNA recognition motifs and a SPOC domain. In Arabidopsis thaliana, the Spen protein FPA is involved in the control of flowering time as a component of an autonomous pathway independent of photoperiod. The A. thaliana genome encodes another gene for a putative Spen protein at the locus At4g12640, herein named AtSpen2. Bioinformatics analysis of the AtSPEN2 SPOC domain revealed low sequence similarity with the FPA SPOC domain, which was markedly lower than that found in other Spen proteins from unrelated plant species. To provide experimental information about the function of AtSpen2, A. thaliana plants were transformed with gene constructs of its promoter region with uidA::gfp reporter genes; the expression was observed in vascular tissues of leaves and roots, as well as in ovules and developing embryos. There was absence of a notable phenotype in knockout and overexpressing lines, suggesting that its function in plants might be specific to certain endogenous or environmental conditions. Our results suggest that the function of Atspen2 diverged from that of fpa due in part to their different transcription expression pattern and divergence of the regulatory SPOC domain.

Comparative Study of Lectin Domains in Model Species: New Insights into Evolutionary Dynamics

Lectins are present throughout the plant kingdom and are reported to be involved in diverse biological processes. In this study, we provide a comparative analysis of the lectin families from model species in a phylogenetic framework. The analysis focuses on the different plant lectin domains identified in five representative core angiosperm genomes (Arabidopsisthaliana, Glycine max, Cucumis sativus, Oryza sativa ssp. japonica and Oryza sativa ssp. indica). The genomes were screened for genes encoding lectin domains using a combination of Basic Local Alignment Search Tool (BLAST), hidden Markov models, and InterProScan analysis. Additionally, phylogenetic relationships were investigated by constructing maximum likelihood phylogenetic trees. The results demonstrate that the majority of the lectin families are present in each of the species under study. Domain organization analysis showed that most identified proteins are multi-domain proteins, owing to the modular rearrangement of protein domains during evolution. Most of these multi-domain proteins are widespread, while others display a lineage-specific distribution. Furthermore, the phylogenetic analyses reveal that some lectin families evolved to be similar to the phylogeny of the plant species, while others share a closer evolutionary history based on the corresponding protein domain architecture. Our results yield insights into the evolutionary relationships and functional divergence of plant lectins.

Misannotation Awareness: A Tale of Two Gene-Groups

Incorrectly or simply not annotated data is largely increasing in most public databases, undoubtedly caused by the rise in sequence data and the more recent boom of genomic projects. Molecular biologists and bioinformaticists should join efforts to tackle this issue. Practical challenges have been experienced when studying the alternative oxidase (AOX) gene family, and hence the motivation for the present work. Commonly used databases were screened for their capacity to distinguish AOX from the plastid terminal oxidase (also called plastoquinol terminal oxidase; PTOX) and we put forward a simple approach, based on amino acids signatures, that unequivocally distinguishes these gene families. Further, available sequence data on the AOX family in plants was carefully revised to: (1) confirm the classification as AOX and (2) identify to which AOX family member they belong to. We bring forward the urgent need of misannotation awareness and re-annotation of public AOX sequences by highlighting different types of misclassifications and the large under-estimation of data availability.

Bioinformatic landscapes for plant transcription factor system research

Diverse bioinformatic resources have been developed for plant transcription factor (TF) research. This review presents the bioinformatic resources and methodologies for the elucidation of plant TF-mediated biological events. Such information is helpful to dissect the transcriptional regulatory systems in the three reference plants Arabidopsis , rice, and maize and translation to other plants. Transcription factors (TFs) orchestrate diverse biological programs by the modulation of spatiotemporal patterns of gene expression via binding cis-regulatory elements. Advanced sequencing platforms accompanied by emerging bioinformatic tools revolutionize the scope and extent of TF research. The system-level integration of bioinformatic resources is beneficial to the decoding of TF-involved networks. Herein, we first briefly introduce general and specialized databases for TF research in three reference plants Arabidopsis, rice, and maize. Then, as proof of concept, we identified and characterized heat shock transcription factor (HSF) members through the TF databases. Finally, we present how the integration of bioinformatic resources at -omics layers can aid the dissection of TF-mediated pathways. We also suggest ways forward to improve the bioinformatic resources of plant TFs. Leveraging these bioinformatic resources and methodologies opens new avenues for the elucidation of transcriptional regulatory systems in the three model systems and translation to other plants.

Cell wall proteome of sugarcane stems: comparison of a destructive and a non-destructive extraction method showed differences in glycoside hydrolases and peroxidases

BACKGROUND

Sugarcane has been used as the main crop for ethanol production for more than 40 years in Brazil. Recently, the production of bioethanol from bagasse and straw, also called second generation (2G) ethanol, became a reality with the first commercial plants started in the USA and Brazil. However, the industrial processes still need to be improved to generate a low cost fuel. One possibility is the remodeling of cell walls, by means of genetic improvement or transgenesis, in order to make the bagasse more accessible to hydrolytic enzymes. We aimed at characterizing the cell wall proteome of young sugarcane culms, to identify proteins involved in cell wall biogenesis. Proteins were extracted from the cell walls of 2-month-old culms using two protocols, non-destructive by vacuum infiltration vs destructive. The proteins were identified by mass spectrometry and bioinformatics.

RESULTS

A predicted signal peptide was found in 84 different proteins, called cell wall proteins (CWPs). As expected, the non-destructive method showed a lower percentage of proteins predicted to be intracellular than the destructive one (33% vs 44%). About 19% of CWPs were identified with both methods, whilst the infiltration protocol could lead to the identification of 75% more CWPs. In both cases, the most populated protein functional classes were those of proteins related to lipid metabolism and oxido-reductases. Curiously, a single glycoside hydrolase (GH) was identified using the non-destructive method whereas 10 GHs were found with the destructive one. Quantitative data analysis allowed the identification of the most abundant proteins.

CONCLUSIONS

The results highlighted the importance of using different protocols to extract proteins from cell walls to expand the coverage of the cell wall proteome. Ten GHs were indicated as possible targets for further studies in order to obtain cell walls less recalcitrant to deconstruction. Therefore, this work contributed to two goals: enlarge the coverage of the sugarcane cell wall proteome, and provide target proteins that could be used in future research to facilitate 2G ethanol production.

Plant protein peptidase inhibitors: an evolutionary overview based on comparative genomics

Background

Peptidases are key proteins involved in essential plant physiological processes. Although protein peptidase inhibitors are essential molecules that modulate peptidase activity, their global presence in different plant species remains still unknown. Comparative genomic analyses are powerful tools to get advanced knowledge into the presence and evolution of both, peptidases and their inhibitors across the Viridiplantae kingdom.

Results

A genomic comparative analysis of peptidase inhibitors and several groups of peptidases in representative species of different plant taxonomic groups has been performed. The results point out: i) clade-specific presence is common to many families of peptidase inhibitors, being some families present in most land plants; ii) variability is a widespread feature for peptidase inhibitory families, with abundant species-specific (or clade-specific) gene family proliferations; iii) peptidases are more conserved in different plant clades, being C1A papain and S8 subtilisin families present in all species analyzed; and iv) a moderate correlation among peptidases and their inhibitors suggests that inhibitors proliferated to control both endogenous and exogenous peptidases.

Conclusions

Comparative genomics has provided valuable insights on plant peptidase inhibitor families and could explain the evolutionary reasons that lead to the current variable repertoire of peptidase inhibitors in specific plant clades.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-812) contains supplementary material, which is available to authorized users.

TRAPID: an efficient online tool for the functional and comparative analysis of de novo RNA-Seq transcriptomes

Transcriptome analysis through next-generation sequencing technologies allows the generation of detailed gene catalogs for non-model species, at the cost of new challenges with regards to computational requirements and bioinformatics expertise. Here, we present TRAPID, an online tool for the fast and efficient processing of assembled RNA-Seq transcriptome data, developed to mitigate these challenges. TRAPID offers high-throughput open reading frame detection, frameshift correction and includes a functional, comparative and phylogenetic toolbox, making use of 175 reference proteomes. Benchmarking and comparison against state-of-the-art transcript analysis tools reveals the efficiency and unique features of the TRAPID system. TRAPID is freely available at http://bioinformatics.psb.ugent.be/webtools/trapid/.

From plant genomes to protein families: computational tools

The development of new high-throughput sequencing technologies has increased dramatically the number of successful genomic projects. Thus, draft genomic sequences of more than 60 plant species are currently available. Suitable bioinformatics tools are being developed to assemble, annotate and analyze the enormous number of sequences produced. In this context, specific plant comparative genomic databases are become powerful tools for gene family annotation in plant clades. In this mini-review, the current state-of-art of genomic projects is glossed. Besides, the computational tools developed to compare genomic data are compiled.

Phylogenetically distant barley legumains have a role in both seed and vegetative tissues

Legumains or vacuolar processing enzymes are cysteine peptidases (C13 family, clan CD) with increasingly recognized physiological significance in plants. They have previously been classified as seed and vegetative legumains. In this work, the entire barley legumain family is described. The eight members of this family belong to the two phylogenetic clades in which the angiosperm legumains are distributed. An in-depth molecular and functional characterization of a barley legumain from each group, HvLeg-2 and HvLeg-4, was performed. Both legumains contained a signal peptide and were located in the endoplasmic reticulum, were expressed in seeds and vegetative tissues, and when expressed as recombinant proteins showed legumain and caspase proteolytic activities. However, the role of each protein seemed to be different in their target tissues. HvLeg-2 responded in leaves to biotic and abiotic stimuli, such as salicylic acid, jasmonic acid, nitric oxide, abscisic acid, and aphid infestation, and was induced by gibberellic acid in seeds, where the protein is able to degrade storage globulins. HvLeg-4 responded in leaves to wounding, nitric oxide, and abscisic acid treatments, and had an unknown role in the germinating seed. From these results, a multifunctional role was assumed for these two phylogenetically distant legumains, achieving different physiological functions in both seed and vegetative tissues.

MIPS PlantsDB: a database framework for comparative plant genome research

The rapidly increasing amount of plant genome (sequence) data enables powerful comparative analyses and integrative approaches and also requires structured and comprehensive information resources. Databases are needed for both model and crop plant organisms and both intuitive search/browse views and comparative genomics tools should communicate the data to researchers and help them interpret it. MIPS PlantsDB (http://mips.helmholtz-muenchen.de/plant/genomes.jsp) was initially described in NAR in 2007 [Spannagl,M., Noubibou,O., Haase,D., Yang,L., Gundlach,H., Hindemitt, T., Klee,K., Haberer,G., Schoof,H. and Mayer,K.F. (2007) MIPSPlantsDB–plant database resource for integrative and comparative plant genome research. Nucleic Acids Res., 35, D834–D840] and was set up from the start to provide data and information resources for individual plant species as well as a framework for integrative and comparative plant genome research. PlantsDB comprises database instances for tomato, Medicago, Arabidopsis, Brachypodium, Sorghum, maize, rice, barley and wheat. Building up on that, state-of-the-art comparative genomics tools such as CrowsNest are integrated to visualize and investigate syntenic relationships between monocot genomes. Results from novel genome analysis strategies targeting the complex and repetitive genomes of triticeae species (wheat and barley) are provided and cross-linked with model species. The MIPS Repeat Element Database (mips-REdat) and Catalog (mips-REcat) as well as tight connections to other databases, e.g. via web services, are further important components of PlantsDB.

Comparative co-expression analysis in plant biology

The analysis of gene expression data generated by high-throughput microarray transcript profiling experiments has shown that transcriptionally coordinated genes are often functionally related. Based on large-scale expression compendia grouping multiple experiments, this guilt-by-association principle has been applied to study modular gene programmes, identify cis-regulatory elements or predict functions for unknown genes in different model plants. Recently, several studies have demonstrated how, through the integration of gene homology and expression information, correlated gene expression patterns can be compared between species. The incorporation of detailed functional annotations as well as experimental data describing protein-protein interactions, phenotypes or tissue specific expression, provides an invaluable source of information to identify conserved gene modules and translate biological knowledge from model organisms to crops. In this review, we describe the different steps required to systematically compare expression data across species. Apart from the technical challenges to compute and display expression networks from multiple species, some future applications of plant comparative transcriptomics are highlighted.

Cysteine peptidases and their inhibitors in Tetranychus urticae: a comparative genomic approach

Background

Cysteine peptidases in the two-spotted spider mite Tetranychus urticae are involved in essential physiological processes, including proteolytic digestion. Cystatins and thyropins are inhibitors of cysteine peptidases that modulate their activity, although their function in this species has yet to be investigated. Comparative genomic analyses are powerful tools to obtain advanced knowledge into the presence and evolution of both, peptidases and their inhibitors, and could aid to elucidate issues concerning the function of these proteins.

Results

We have performed a genomic comparative analysis of cysteine peptidases and their inhibitors in T. urticae and representative species of different arthropod taxonomic groups. The results indicate: i) clade-specific proliferations are common to C1A papain-like peptidases and for the I25B cystatin family of inhibitors, whereas the C1A inhibitors thyropins are evolutionarily more conserved among arthropod clades; ii) an unprecedented extensive expansion for C13 legumain-like peptidases is found in T. urticae; iii) a sequence-structure analysis of the spider mite cystatins suggests that diversification may be related to an expansion of their inhibitory range; and iv) an in silico transcriptomic analysis shows that most cathepsin B and L cysteine peptidases, legumains and several members of the cystatin family are expressed at a higher rate in T. urticae feeding stages than in embryos.

Conclusion

Comparative genomics has provided valuable insights on the spider mite cysteine peptidases and their inhibitors. Mite-specific proliferations of C1A and C13 peptidase and I25 cystatin families and their over-expression in feeding stages of mites fit with a putative role in mite’s feeding and could have a key role in its broad host feeding range.

Update on sumoylation: defining core components of the plant SUMO conjugation system by phylogenetic comparison

The conjugation of the small ubiquitin-related modifier, SUMO, to substrate proteins is a reversible and dynamic process, and an important response of plants to environmental challenges. Nevertheless, reliable data have so far been restricted largely to the model plant Arabidopsis thaliana. The increasing availability of genome information for other plant species offers the possibility to identify a core set of indispensable components, and to discover species-specific features of the sumoylation pathway. We analyzed the enzymes responsible for the conjugation of SUMO to substrates for their conservation between dicots and monocots. We thus assembled gene sets that relate the Arabidopsis SUMO conjugation system to that of the dicot species tomato, grapevine and poplar, and to four plant species from the monocot class: rice, Brachypodium distachyon, Sorghum bicolor and maize. We found that a core set of genes with clear assignment in Arabidopsis had highly conserved homologs in all tested plants. However, we also observed a variation in the copy number of homologous genes, and sequence variations that suggested monocot-specific variants. Generally, SUMO ligases and proteases showed the most pronounced differences. Finally, we identified potential SUMO chain-binding ubiquitin ligases, pointing to an in vivo function of SUMO chains as degradation signals in plants.

Birth, death and subfunctionalization in the Arabidopsis genome

Arabidopsis thaliana is now a model system, not just for plant biology but also for comparative genomics. The completion of the sequences of two closely related species, Arabidopsis lyrata and Brassica rapa, is complemented by genomic comparisons among A. thaliana accessions and mutation accumulation lines. Together these genomic data document the birth of new genes via gene duplication, transposon exaptation and de novo formation of new genes from noncoding sequence. Most novel loci exhibit low expression, and are undergoing pseudogenization or subfunctionalization. Comparatively, A. thaliana has lost large amounts of sequence through deletion, particularly of transposable elements. Intraspecific genomic variation indicates high rates of deletion mutations and deletion polymorphisms across accessions, shedding light on the history of Arabidopsis genome architecture.