Evolutionary dynamics of protein domain architecture in plants

Phylogenetic and Evolutionary Analysis of Plant Small RNA 2'-O-Methyltransferase (HEN1) Protein Family

HUA ENHANCER 1 (HEN1) is a pivotal mediator in protecting sRNAs from 3'-end uridylation and 3' to 5' exonuclease-mediated degradation in plants. Here, we investigated the pattern of the HEN1 protein family evolutionary history and possible relationships in the plant lineages using protein sequence analyses and conserved motifs composition, functional domain identification, architecture, and phylogenetic tree reconstruction and evolutionary history inference. According to our results, HEN1 protein sequences bear several highly conserved motifs in plant species retained during the evolution from their ancestor. However, several motifs are present only in Gymnosperms and Angiosperms. A similar trend showed for their domain architecture. At the same time, phylogenetic analysis revealed the grouping of the HEN1 proteins in the three main super clads. In addition, the Neighbor-net network analysis result provides some nodes have multiple parents indicating a few conflicting signals in the data, which is not the consequence of sampling error, the effect of the selected model, or the estimation method. By reconciling the protein and species tree, we considered the gene duplications in several given species and found 170 duplication events in the evolution of HEN1 in the plant lineages. According to our analysis, the main HEN1 superclass mostly showed orthologous sequences that illustrate the vertically transmitting of HEN1 to the main lines. However, in both orthologous and paralogs, we predicted insignificant structural deviations. Our analysis implies that small local structural changes that occur continuously during the folds can moderate the changes created in the sequence. According to our results, we proposed a hypothetical model and evolutionary trajectory for the HEN1 protein family in the plant kingdom.

Recent Insights into Plant miRNA Biogenesis: Multiple Layers of miRNA Level Regulation

MicroRNAs are small RNAs, 20-22 nt long, the main role of which is to downregulate gene expression at the level of mRNAs. MiRNAs are fundamental regulators of plant growth and development in response to internal signals as well as in response to abiotic and biotic factors. Therefore, the deficiency or excess of individual miRNAs is detrimental to particular aspects of a plant's life. In consequence, the miRNA levels must be appropriately adjusted. To obtain proper expression of each miRNA, their biogenesis is controlled at multiple regulatory layers. Here, we addressed processes discovered to influence miRNA steady-state levels, such as MIR transcription, co-transcriptional pri-miRNA processing (including splicing, polyadenylation, microprocessor assembly and activity) and miRNA-encoded peptides synthesis. MiRNA stability, RISC formation and miRNA export out of the nucleus and out of the plant cell also define the levels of miRNAs in various plant tissues. Moreover, we show the evolutionary conservation of miRNA biogenesis core proteins across the plant kingdom.

Comparative phylogeny and evolutionary analysis of Dicer-like protein family in two plant monophyletic lineages

BACKGROUND

Small RNAs (sRNAs) that do not get untranslated into proteins exhibit a pivotal role in the expression regulation of their cognate gene(s) in almost all eukaryotic lineages, including plants. Hitherto, numerous protein families such as Dicer, a unique class of Ribonuclease III, have been reported to be involved in sRNAs processing pathways and silencing. In this study, we aimed to investigate the phylogenetic relationship and evolutionary history of the DCL protein family.

RESULTS

Our results illustrated the DCL family of proteins grouped into four main subfamilies (DCLs 1-4) presented in either Eudicotyledons or Liliopsids. The accurate observation of the phylogenetic trees supports the independent expansion of DCL proteins among the Eudicotyledons and Liliopsids species. They share the common origin, and the main duplication events for the formation of the DCL subfamilies occurred before the Eudicotyledons/Liliopsids split from their ancestral DCL. In addition, shreds of evidence revealed that the divergence happened when multicellularization started and since the need for complex gene regulation considered being a necessity by organisms. At that time, they have evolved independently among the monophyletic lineages. The other finding was that the combination of DCL protein subfamilies bears several highly conserved functional domains in plant species that originated from their ancestor architecture. The conservation of these domains happens to be both lineage-specific and inter lineage-specific.

CONCLUSIONS

DCL subfamilies (i.e., DCL1-DCL4) distribute in their single clades after diverging from their common ancestor and before emerging into higher plants. Therefore, it seems that the main duplication events for the formation of the DCL subfamilies occurred before the Eudicotyledons/Liliopsida split and before the appearance of moss, and after the single-cell green algae. We also observed the same trends among the main DCL subfamilies from functional unit composition and architecture. Despite the long evolutionary course from the divergence of Liliopsida lineage from the Eudicotyledons, a significant diversifying force to domain composition and orientation was absent. The results of this study provide a deeper insight into DCL protein evolutionary history and possible sequence and structural relationships between DCL protein subfamilies in the main higher plant monophyletic lineages; i.e., Eudicotyledons and Liliopsida.

TIAMMAt: Leveraging biodiversity to revise protein domain models, evidence from innate immunity

Sequence annotation is fundamental for studying the evolution of protein families, particularly when working with nonmodel species. Given the rapid, ever-increasing number of species receiving high-quality genome sequencing, accurate domain modeling that is representative of species diversity is crucial for understanding protein family sequence evolution and their inferred function(s). Here, we describe a bioinformatic tool called Taxon-Informed Adjustment of Markov Model Attributes (TIAMMAt) which revises domain profile hidden Markov models (HMMs) by incorporating homologous domain sequences from underrepresented and nonmodel species. Using innate immunity pathways as a case study, we show that revising profile HMM parameters to directly account for variation in homologs among underrepresented species provides valuable insight into the evolution of protein families. Following adjustment by TIAMMAt, domain profile HMMs exhibit changes in their per-site amino acid state emission probabilities and insertion/deletion probabilities while maintaining the overall structure of the consensus sequence. Our results show that domain revision can heavily impact evolutionary interpretations for some families (i.e., NLR’s NACHT domain), whereas impact on other domains (e.g., rel homology domain and interferon regulatory factor domains) is minimal due to high levels of sequence conservation across the sampled phylogenetic depth (i.e., Metazoa). Importantly, TIAMMAt revises target domain models to reflect homologous sequence variation using the taxonomic distribution under consideration by the user. TIAMMAt’s flexibility to revise any subset of the Pfam database using a user-defined taxonomic pool will make it a valuable tool for future protein evolution studies, particularly when incorporating (or focusing) on nonmodel species.

CsubMADS1, a lag phase transcription factor, controls development of polar eukaryotic microalga Coccomyxa subellipsoidea C-169

MADS-box transcription factors (TFs) have not been functionally delineated in microalgae. In this study, the role of CsubMADS1 from microalga Coccomyxa subellipsoidea C-169 has been explored. Unlike Type II MADS-box proteins of seed plants with MADS, Intervening, K-box, and C domains, CsubMADS1 only has MADS and Intervening domains. It forms a group with MADS TFs from algae in the phylogenetic tree within the Type II MIKC^C clade. CsubMADS1 is expressed strongly in the lag phase of growth. The CsubMADS1 monomer does not have a specific localization in the nucleus, and it forms homodimers to localize exclusively in the nucleus. The monomer has two nuclear localization signals (NLSs): an N-terminal NLS and an internal NLS. The internal NLS is functional, and the homodimer requires two NLSs for specific nuclear localization. Overexpression (OX) of CsubMADS1 slows down the growth of the culture and leads to the creation of giant polyploid multinucleate cells, resembling autospore mother cells. This implies that the release of autospores from autospore mother cells may be delayed. Thus, in wild-type (WT) cells, CsubMADS1 may play a crucial role in slowing down growth during the lag phase. Due to starvation in 2-month-old colonies on solid media, the WT colonies produce mucilage, whereas OX colonies produce significantly less mucilage. Thus, CsubMADS1 also negatively regulates stress-induced mucilage production and probably plays a role in stress tolerance during the lag phase. Taken together, our results reveal that CsubMADS1 is a key TF involved in the development and stress tolerance of this polar microalga.

Mix and match: Patchwork domain evolution of the land plant-specific Ca2+-permeable mechanosensitive channel MCA

Multidomain proteins can have a complex evolutionary history that may involve de novo domain evolution, recruitment and / or recombination of existing domains and domain losses. Here, the domain evolution of the plant-specific Ca²⁺-permeable mechanosensitive channel protein, MID1-COMPLEMENTING ACTIVITY (MCA), was investigated. MCA, a multidomain protein, possesses a Ca²⁺-influx-MCA^func domain and a PLAC8 domain. Profile Hidden Markov Models (HMMs) of domains were assessed in 25 viridiplantae proteomes. While PLAC8 was detected in plants, animals, and fungi, MCA^func was found in streptophytes but not in chlorophytes. Full MCA proteins were only found in embryophytes. We identified the MCA^func domain in all streptophytes including charophytes where it appeared in E3 ubiquitin ligase-like proteins. Our Maximum Likelihood (ML) analyses suggested that the MCA^func domain evolved early in the history of streptophytes. The PLAC8 domain showed similarity to Plant Cadmium Resistance (PCR) genes, and the coupling of MCA^func and PLAC8 seemed to represent a single evolutionary event. This combination is unique in MCA, and does not exist in other plant mechanosensitive channels. Within angiosperms, gene duplications increased the number of MCAs. Considering their role in mechanosensing in roots, MCA might be instrumental for the rise of land plants. This study provides a textbook example of de novo domain emergence, recombination, duplication, and losses, leading to the convergence of function of proteins in plants.

Using Cellular Automata to Simulate Domain Evolution in Proteins

Proteins play primary roles in important biological processes such as catalysis, physiological functions, and immune system functions. Thus, the research on how proteins evolved has been a nuclear question in the field of evolutionary biology. General models of protein evolution help to determine the baseline expectations for evolution of sequences, and these models have been extensively useful in sequence analysis as well as for the computer simulation of artificial sequence data sets. We have developed a new method of simulating multi-domain protein evolution, including fusions of domains, insertion, and deletion. It has been observed via the simulation test that the success rates achieved by the proposed predictor are remarkably high. For the convenience of the most experimental scientists, a user-friendly web server has been established at http://jci-bioinfo.cn/domainevo, by which users can easily get their desired results without having to go through the detailed mathematics. Through the simulation results of this website, users can predict the evolution trend of the protein domain architecture.

Evolutionary Analysis of JAZ Proteins in Plants: An Approach in Search of the Ancestral Sequence

Jasmonates are phytohormones that regulate development, metabolism and immunity. Signal transduction is critical to activate jasmonate responses, but the evolution of some key regulators such as jasmonate-ZIM domain (JAZ) repressors is not clear. Here, we identified 1065 JAZ sequence proteins in 66 lower and higher plants and analyzed their evolution by bioinformatics methods. We found that the TIFY and Jas domains are highly conserved along the evolutionary scale. Furthermore, the canonical degron sequence LPIAR(R/K) of the Jas domain is conserved in lower and higher plants. It is noteworthy that degron sequences showed a large number of alternatives from gymnosperms to dicots. In addition, ethylene-responsive element binding factor-associated amphiphilic repression (EAR) motifs are displayed in all plant lineages from liverworts to angiosperms. However, the cryptic MYC2-interacting domain (CMID) domain appeared in angiosperms for the first time. The phylogenetic analysis performed using the Maximum Likelihood method indicated that JAZ ortholog proteins are grouped according to their similarity and plant lineage. Moreover, ancestral JAZ sequences were constructed by PhyloBot software and showed specific changes in the TIFY and Jas domains during evolution from liverworts to dicots. Finally, we propose a model for the evolution of the ancestral sequences of the main eight JAZ protein subgroups. These findings contribute to the understanding of the JAZ family origin and expansion in land plants.

Domain-based Comparative Analysis of Bacterial Proteomes: Uniqueness, Interactions, and the Dark Matter

Background

Proteins may have none, single, double, or multiple domains, while a single domain may appear in multiple proteins. Their distribution patterns may have impacts on bacterial physi-ology and lifestyle.Objective: This study aims to understand how domains are distributed and duplicated in bacterial prote-omes, in order to better understand bacterial physiology and lifestyles.

Methods

In this study, we used 16712 Hidden Markov Models to screen 944 bacterial reference prote-omes versus a threshold E-value<0.001. The number of non-redundant domains and duplication rates of redundant domains for each species were calculated. The unique domains, if any, were also identified for each species. In addition, the properties of no-domain proteins were investigated in terms of physico-chemical properties.

Results

The increasing number of non-redundant domains for a bacterial proteome follows the trend of an asymptotic function. The domain duplication rate is positively correlated with proteome size and in-creases more rapidly. The high percentage of single-domain proteins is more associated with small pro-teome size. For each proteome, unique domains were also obtained. Moreover, no-domain proteins show differences with the other three groups for several physicochemical properties analysed in this study.

Conclusion

The study confirmed that a low domain duplication rate and a high percentage of single-domain proteins are more likely to be associated with bacterial host-dependent or restricted niche-adapted lifestyle. In addition, the unique lifestyle and physiology were revealed based on the analysis of species-specific domains and core domain interactions or co-occurrences.

Evolution of Protein Domain Architectures

This chapter reviews current research on how protein domain architectures evolve. We begin by summarizing work on the phylogenetic distribution of proteins, as this will directly impact which domain architectures can be formed in different species. Studies relating domain family size to occurrence have shown that they generally follow power law distributions, both within genomes and larger evolutionary groups. These findings were subsequently extended to multi-domain architectures. Genome evolution models that have been suggested to explain the shape of these distributions are reviewed, as well as evidence for selective pressure to expand certain domain families more than others. Each domain has an intrinsic combinatorial propensity, and the effects of this have been studied using measures of domain versatility or promiscuity. Next, we study the principles of protein domain architecture evolution and how these have been inferred from distributions of extant domain arrangements. Following this, we review inferences of ancestral domain architecture and the conclusions concerning domain architecture evolution mechanisms that can be drawn from these. Finally, we examine whether all known cases of a given domain architecture can be assumed to have a single common origin (monophyly) or have evolved convergently (polyphyly). We end by a discussion of some available tools for computational analysis or exploitation of protein domain architectures and their evolution.

Order in Disorder as Observed by the "Hydrophobic Cluster Analysis" of Protein Sequences

Hydrophobic cluster analysis (HCA) is an original approach for protein sequence analysis, which provides access to the foldable repertoire of the protein universe, including yet unannotated protein segments ("dark proteome"). Foldable segments correspond to ordered regions, as well as to intrinsically disordered regions (IDRs) undergoing disorder to order transitions. In this review, how HCA can be used to give insight into this last category of foldable segments is illustrated, with examples matching known 3D structures. After reviewing the HCA principles, examples of short foldable segments are given, which often contain short linear motifs, typically matching hydrophobic clusters. These segments become ordered upon contact with partners, with secondary structure preferences generally corresponding to those observed in the 3D structures within the complexes. Such small foldable segments are sometimes larger than the segments of known 3D structures, including flanking hydrophobic clusters that may be critical for interaction specificity or regulation, as well as intervening sequences allowing fuzziness. Cases of larger conditionally disordered domains are also presented, with lower density in hydrophobic clusters than well-folded globular domains or with exposed hydrophobic patches, which are stabilized by interaction with partners.

Identification and characterization of candidates involved in production of OMEGAs in microalgae: a gene mining and phylogenomic approach

Optimizing the production of the high-value renewables such as OMEGAs through pathway engineering requires an in-depth understanding of the structure-function relationship of genes involved in the OMEGA biosynthetic pathways. In this preliminary study, our rationale is to identify and characterize the ∼221 putative genes involved in production of OMEGAs using bioinformatic analysis from the Streptophyte (plants), Chlorophyte (green algae), Rhodophyta (red algae), and Bacillariophyta (diatoms) lineages based on their phylogenomic profiling, conserved motif/domain organization and physico-chemical properties. The MEME suite predicted 12 distinct protein domains, which are conserved among these putative genes. The phylogenomic analysis of the putative candidate genes [such as FAD2 (delta-12 desaturase); ECR (enoyl-CoA reductase); FAD2 (delta-12 desaturase); ACOT (acyl CoA thioesterase); ECH (enoyl-CoA hydratase); and ACAT (acetyl-CoA acyltransferase)] with similar domains and motif patterns were remarkably well conserved. Furthermore, the subcellular network prediction of OMEGA biosynthetic pathway genes revealed a unique interaction between the light-dependent chlorophyll biosynthesis and glycerol-3-phosphate dehydrogenase, which predicts a major cross-talk between the key essential pathways. Such bioinformatic analysis will provide insights in finding the key regulatory genes to optimize the productivity of OMEGAs in microalgal cell factories.

Evolution of DDB1-binding WD40 (DWD) in the viridiplantae

Damaged DNA Binding 1 (DDB1)—binding WD40 (DWD) proteins are highly conserved and involved in a plethora of developmental and physiological processes such as flowering time control, photomorphogenesis, and abiotic stress responses. The phylogeny of this family of proteins in plants and algae of viridiplante is a critical area to understand the emergence of this family in such important and diverse functions. We aimed to investigate the putative homologs of DWD in the viridiplante and establish a deeper DWD evolutionary grasp. The advancement in publicly available genomic data allowed us to perform an extensive genome-wide DWD retrieval. Using annotated Arabidopsis thaliana DWDs as the reference, we generated and characterized a comprehensive DWD database for the studied photoautotrophs. Further, a generic DWD classification system (Type A to K), based on (i) position of DWD motifs, (ii) number of DWD motifs, and (iii) presence/absence of other domains, was adopted. About 72–80% DWDs have one DWD motif, whereas 17–24% DWDs have two and 0.5–4.7% DWDs have three DWD motifs. Neighbor-joining phylogenetic construction of A. thaliana DWDs facilitated us to tune these substrate receptors into 15 groups. Though the DWD count increases from microalgae to higher land plants, the ratio of DWD to WD40 remained constant throughout the viridiplante. The DWD expansion appeared to be the consequence of consistent DWD genetic flow accompanied by several gene duplication events. The network, phylogenetic, and statistical analysis delineated DWD evolutionary relevance in the viridiplante.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

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.

Analyses of Old "Prokaryotic" Proteins Indicate Functional Diversification in Arabidopsis and Oryza sativa

During evolution, various processes such as duplication, divergence, recombination, and many other events leads to the evolution of new genes with novel functions. These evolutionary events, thus significantly impact the evolution of cellular, physiological, morphological, and other phenotypic trait of organisms. While evolving, eukaryotes have acquired large number of genes from the earlier prokaryotes. This work is focused upon identification of old "prokaryotic" proteins in Arabidopsis and Oryza sativa genome, further highlighting their possible role(s) in the two genomes. Our results suggest that with respect to their genome size, the fraction of old "prokaryotic" proteins is higher in Arabidopsis than in Oryza sativa. The large fractions of such proteins encoding genes were found to be localized in various endo-symbiotic organelles. The domain architecture of the old "prokaryotic" proteins revealed similar distribution in both Arabidopsis and Oryza sativa genomes showing their conserved evolution. In Oryza sativa, the old "prokaryotic" proteins were more involved in developmental processes, might be due to constant man-made selection pressure for better agronomic traits/productivity. While in Arabidopsis, these proteins were involved in metabolic functions. Overall, the analysis indicates the distinct pattern of evolution of old "prokaryotic" proteins in Arabidopsis and Oryza sativa.

Presence of unique glyoxalase III proteins in plants indicates the existence of shorter route for methylglyoxal detoxification

Glyoxalase pathway, comprising glyoxalase I (GLY I) and glyoxalase II (GLY II) enzymes, is the major pathway for detoxification of methylglyoxal (MG) into D-lactate involving reduced glutathione (GSH). However, in bacteria, glyoxalase III (GLY III) with DJ-1/PfpI domain(s) can do the same conversion in a single step without GSH. Our investigations for the presence of DJ-1/PfpI domain containing proteins in plants have indicated the existence of GLY III-like proteins in monocots, dicots, lycopods, gymnosperm and bryophytes. A deeper in silico analysis of rice genome identified twelve DJ-1 proteins encoded by six genes. Detailed analysis has been carried out including their chromosomal distribution, genomic architecture and localization. Transcript profiling under multiple stress conditions indicated strong induction of OsDJ-1 in response to exogenous MG. A member of OsDJ-1 family, OsDJ-1C, showed high constitutive expression at all developmental stages and tissues of rice. MG depletion study complemented by simultaneous formation of D-lactate proved OsDJ-1C to be a GLY III enzyme that converts MG directly into D-lactate in a GSH-independent manner. Site directed mutagenesis of Cys-119 to Alanine significantly reduces its GLY III activity indicating towards the existence of functional GLY III enzyme in rice—a shorter route for MG detoxification.

Distribution and evolution of the lectin family in soybean (Glycine max)

Lectins are a diverse group of proteins that bind specific carbohydrates and are found throughout all kingdoms. In plants, lectins are involved in a range of important processes such as plant defense and stress signaling. Although the genome sequence of Glycine max (soybean) has been published, little is known about the abundance and expansion patterns of lectin genes in soybean. Using BLAST and hidden Markov models, a total of 359 putative lectin genes have been identified. Furthermore, these sequences could be classified in nine of the twelve plant lectin families identified today. Analysis of the domain organization demonstrated that most of the identified lectin genes encode chimerolectins, consisting of one or multiple lectin domains combined with other known protein domains. Both tandem and segmental duplication events have contributed to the expansion of the lectin gene family. These data provide a detailed understanding of the domain architecture and molecular evolution of the lectin gene family in soybean.

P³DB 3.0: From plant phosphorylation sites to protein networks

In the past few years, the Plant Protein Phosphorylation Database (P³DB, http://p3db.org) has become one of the most significant in vivo data resources for studying plant phosphoproteomics. We have substantially updated P³DB with respect to format, new datasets and analytic tools. In the P³DB 3.0, there are altogether 47 923 phosphosites in 16 477 phosphoproteins curated across nine plant organisms from 32 studies, which have met our multiple quality standards for acquisition of in vivo phosphorylation site data. Centralized by these phosphorylation data, multiple related data and annotations are provided, including protein–protein interaction (PPI), gene ontology, protein tertiary structures, orthologous sequences, kinase/phosphatase classification and Kinase Client Assay (KiC Assay) data—all of which provides context for the phosphorylation event. In addition, P³DB 3.0 incorporates multiple network viewers for the above features, such as PPI network, kinase-substrate network, phosphatase-substrate network, and domain co-occurrence network to help study phosphorylation from a systems point of view. Furthermore, the new P³DB reflects a community-based design through which users can share datasets and automate data depository processes for publication purposes. Each of these new features supports the goal of making P³DB a comprehensive, systematic and interactive platform for phosphoproteomics research.

SAPs as novel regulators of abiotic stress response in plants

Stress associated proteins (SAPs), novel A20/AN1 zinc-finger domain-containing proteins, are fast emerging as potential candidates for biotechnological approaches in order to improve abiotic stress tolerance in plants - the ultimate aim of which is crop-yield protection. Until relatively recently, such proteins had only been identified in humans, where they had been shown to be key regulators of innate immunity. Their phylogenetic relationship and recruitment of diverse protein domains reflect an architectural and mechanistic diversity. Emerging evidence suggests that SAPs may act as ubiquitin ligase, redox sensor, and regulator of gene expression during stress. Here, we evaluate the new knowledge on SAPs with a view to understand their mechanism of action. Furthermore, we set an agenda for investigating hitherto unexplored roles of these proteins.

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the 'core' proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome-wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.

Variation of presence/absence genes among Arabidopsis populations

BACKGROUND

Gene presence/absence (P/A) polymorphisms are commonly observed in plants and are important in individual adaptation and species differentiation. Detecting their abundance, distribution and variation among individuals would help to understand the role played by these polymorphisms in a given species. The recently sequenced 80 Arabidopsis genomes provide an opportunity to address these questions.

RESULTS

By systematically investigating these accessions, we identified 2,407 P/A genes (or 8.9%) absent in one or more genomes, averaging 444 absent genes per accession. 50.6% of P/A genes belonged to multi-copy gene families, or 31.0% to clustered genes. However, the highest proportion of P/A genes, outnumbered in singleton genes, was observed in the regions near centromeres. In addition, a significant correlation was observed between the P/A gene frequency among the 80 accessions and the diversity level at P/A loci. Furthermore, the proportion of P/A genes was different among functional gene categories. Finally, a P/A gene tree showed a diversified population structure in the worldwide Arabidopsis accessions.

CONCLUSIONS

An estimate of P/A genes and their frequency distribution in the worldwide Arabidopsis accessions was obtained. Our results suggest that there are diverse mechanisms to generate or maintain P/A genes, by which individuals and functionally different genes can selectively maintain P/A polymorphisms for a specific adaptation.