Protein subcellular relocalization of duplicated genes in Arabidopsis

Protein subcellular relocalization and function of duplicated flagellar calcium binding protein genes in honey bee trypanosomatid parasite

The honey bee trypanosomatid parasite, Lotmaria passim, contains two genes that encode the flagellar calcium binding protein (FCaBP) through tandem duplication in its genome. FCaBPs localize in the flagellum and entire body membrane of L. passim through specific N-terminal sorting sequences. This finding suggests that this is an example of protein subcellular relocalization resulting from gene duplication, altering the intracellular localization of FCaBP. However, this phenomenon may not have occurred in Leishmania, as one or both of the duplicated genes have become pseudogenes. Multiple copies of the FCaBP gene are present in several Trypanosoma species and Leptomonas pyrrhocoris, indicating rapid evolution of this gene in trypanosomatid parasites. The N-terminal flagellar sorting sequence of L. passim FCaBP1 is in close proximity to the BBSome complex, while that of Trypanosoma brucei FCaBP does not direct GFP to the flagellum in L. passim. Deletion of the two FCaBP genes in L. passim affected growth and impaired flagellar morphogenesis and motility, but it did not impact host infection. Therefore, FCaBP represents a duplicated gene with a rapid evolutionary history that is essential for flagellar structure and function in a trypanosomatid parasite.

A tyrosine phospho-switch within the Longin domain of VAMP721 modulates SNARE functionality

The final step in secretion is membrane fusion facilitated by SNARE proteins that reside in opposite membranes. The formation of a trans-SNARE complex between one R and three Q coiled-coiled SNARE domains drives the final approach of the membranes providing the mechanical energy for fusion. Biological control of this mechanism is exerted by additional domains within some SNAREs. For example, the N-terminal Longin domain (LD) of R-SNAREs (also called Vesicle-associated membrane proteins, VAMPs) can fold back onto the SNARE domain blocking interaction with other cognate SNAREs. The LD may also determine the subcellular localization via interaction with other trafficking-related proteins. Here, we provide cell-biological and genetic evidence that phosphorylation of the Tyrosine57 residue regulates the functionality of VAMP721. We found that an aspartate mutation mimics phosphorylation, leading to protein instability and subsequent degradation in lytic vacuoles. The mutant SNARE also fails to rescue the defects of vamp721vamp722 loss-of-function lines in spite of its wildtype-like localization within the secretory pathway and the ability to interact with cognate SNARE partners. Most importantly, it imposes a dominant negative phenotype interfering with root growth, normal secretion and cytokinesis in wildtype plants generating large aggregates that mainly contain secretory vesicles. Non-phosphorylatable VAMP721Y57F needs higher gene dosage to rescue double mutants in comparison to native VAMP721 underpinning that phosphorylation modulates SNARE function. We propose a model where short-lived phosphorylation of Y57 serves as a regulatory step to control VAMP721 activity, favoring its open state and interaction with cognate partners to ultimately drive membrane fusion.

Systematic investigation of the R2R3-MYB gene family in Aquilaria sinensis reveals a transcriptional repressor AsMYB054 involved in 2-(2-phenylethyl)chromone biosynthesis

Trees in the genus Aquilaria produce agarwood, a valuable resin used in medicine, perfumes, and incense. 2-(2-Phenethyl)chromones (PECs) are characteristic components of agarwood; however, molecular mechanisms underlying PEC biosynthesis and regulation remain largely unknown. The R2R3-MYB transcription factors play important regulatory roles in the biosynthesis of various secondary metabolites. In this study, 101 R2R3-MYB genes in Aquilaria sinensis were systematically identified and analyzed at the genome-wide level. Transcriptomic analysis revealed that 19 R2R3-MYB genes were significantly regulated by an agarwood inducer, and showed significant correlations with PEC accumulation. Expression and evolutionary analyses revealed that AsMYB054, a subgroup 4 R2R3-MYB, was negatively correlated with PEC accumulation. AsMYB054 was located in the nucleus and functioned as a transcriptional repressor. Moreover, AsMYB054 could bind to the promoters of the PEC biosynthesis related genes AsPKS02 and AsPKS09, and inhibit their transcriptional activity. These findings suggested that AsMYB054 functions as a negative regulator of PEC biosynthesis via the inhibition of AsPKS02 and AsPKS09 in A. sinensis. Our results provide a comprehensive understanding of the R2R3-MYB subfamily in A. sinensis and lay a foundation for further functional analyses of R2R3-MYB genes in PEC biosynthesis.

Genome-wide identification and comparative analysis of lipocalin families in Lepidoptera with an emphasis on Bombyx mori

Lipocalins exhibit functional diversity, including roles in retinol transport, invertebrate cryptic coloration, and stress response. However, genome-wide identification and characterization of lipocalin in the insect lineage have not been thoroughly explored. Here, we found that a lineage-specific expansion of the lipocalin genes in Lepidoptera occurred in large part due to tandem duplication events and several lipocalin genes involving insect coloration were expanded more via tandem duplication in butterflies. A comparative analysis of conserved motifs showed both conservation and divergence of lepidopteran lipocalin family protein structures during evolution. We observe dynamic changes in tissue expression preference of paralogs in Bombyx mori, suggesting differential contribution of paralogs to specific organ functions during evolution. Subcellular localization experiments revealed that lipocalins localize to the cytoplasm, nuclear membrane, or nucleus in BmN cells. Moreover, several lipocalin genes exhibited divergent responses to abiotic and biotic stresses, and 1 lipocalin gene was upregulated by 300 fold in B. mori. These results suggest that lipocalins act as signaling components in defense responses by mediating crosstalk between abiotic and biotic stress responses. This study deepens our understanding of the comprehensive characteristics of lipocalins in insects.

A metabolic, phylogenomic and environmental atlas of diatom plastid transporters from the model species Phaeodactylum

Diatoms are an important group of algae, contributing nearly 40% of total marine photosynthetic activity. However, the specific molecular agents and transporters underpinning the metabolic efficiency of the diatom plastid remain to be revealed. We performed in silico analyses of 70 predicted plastid transporters identified by genome-wide searches of Phaeodactylum tricornutum. We considered similarity with Arabidopsis thaliana plastid transporters, transcriptional co-regulation with genes encoding core plastid metabolic pathways and with genes encoded in the mitochondrial genomes, inferred evolutionary histories using single-gene phylogeny, and environmental expression trends using Tara Oceans meta-transcriptomics and meta-genomes data. Our data reveal diatoms conserve some of the ion, nucleotide and sugar plastid transporters associated with plants, such as non-specific triose phosphate transporters implicated in the transport of phosphorylated sugars, NTP/NDP and cation exchange transporters. However, our data also highlight the presence of diatom-specific transporter functions, such as carbon and amino acid transporters implicated in intricate plastid-mitochondria crosstalk events. These confirm previous observations that substrate non-specific triose phosphate transporters (TPT) may exist as principal transporters of phosphorylated sugars into and out of the diatom plastid, alongside suggesting probable agents of NTP exchange. Carbon and amino acid transport may be related to intricate metabolic plastid-mitochondria crosstalk. We additionally provide evidence from environmental meta-transcriptomic/meta- genomic data that plastid transporters may underpin diatom sensitivity to ocean warming, and identify a diatom plastid transporter (J43171) whose expression may be positively correlated with temperature.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Gene Duplication Accelerates the Pace of Protein Gain and Loss from Plant Organelles

Organelle biogenesis and function is dependent on the concerted action of both organellar-encoded (if present) and nuclear-encoded proteins. Differences between homologous organelles across the Plant Kingdom arise, in part, as a result of differences in the cohort of nuclear-encoded proteins that are targeted to them. However, neither the rate at which differences in protein targeting accumulate nor the evolutionary consequences of these changes are known. Using phylogenomic approaches coupled to ancestral state estimation, we show that the plant organellar proteome has diversified in proportion with molecular sequence evolution such that the proteomes of plant chloroplasts and mitochondria lose or gain on average 3.6 proteins per million years. We further demonstrate that changes in organellar protein targeting are associated with an increase in the rate of molecular sequence evolution and that such changes predominantly occur in genes with regulatory rather than metabolic functions. Finally, we show that gain and loss of protein target signals occurs at a higher rate following gene duplication, revealing that gene and genome duplication are a key facilitator of plant organelle evolution.

Genome-Wide Comparative Analysis of R2R3 MYB Gene Family in Populus and Salix and Identification of Male Flower Bud Development-Related Genes

The MYB transcription factor (TF) family is one of the largest plant transcription factor gene family playing vital roles in plant growth and development, including defense, cell differentiation, secondary metabolism, and responses to biotic and abiotic stresses. As a model tree species of woody plants, in recent years, the identification and functional prediction of certain MYB family members in the poplar genome have been reported. However, to date, the characterization of the gene family in the genome of the poplar's sister species willow has not been done, nor are the differences and similarities between the poplar and willow genomes understood. In this study, we conducted the first genome-wide investigation of the R2R3 MYB subfamily in the willow, identifying 216 R2R3 MYB gene members, and combined with the poplar R2R3 MYB genes, performed the first comparative analysis of R2R3 MYB genes between the poplar and willow. We identified 81 and 86 pairs of R2R3 MYB paralogs in the poplar and willow, respectively. There were 17 pairs of tandem repeat genes in the willow, indicating active duplication of willow R2R3 MYB genes. A further 166 pairs of poplar and willow orthologs were identified by collinear and synonymous analysis. The findings support the duplication of R2R3 MYB genes in the ancestral species, with most of the R2R3 MYB genes being retained during the evolutionary process. The phylogenetic trees of the R2R3 MYB genes of 10 different species were drawn. The functions of the poplar and willow R2R3 MYB genes were predicted using reported functional groupings and clustering by OrthoFinder. Identified 5 subgroups in general expanded in woody species, three subgroups were predicted to be related to lignin synthesis, and we further speculate that the other two subgroups also play a role in wood formation. We analyzed the expression patterns of the GAMYB gene of subgroup 18 (S18) related to pollen development in the male flower buds of poplar and willow at different developmental stages by qRT-PCR. The results showed that the GAMYB gene was specifically expressed in the male flower bud from pollen formation to maturity, and that the expression first increased and then decreased. Both the specificity of tissue expression specificity and conservation indicated that GAMYB played an important role in pollen development in both poplar and willow and was an ideal candidate gene for the analysis of male flower development-related functions of the two species.

Divergence of duplicated genes by repeated partitioning of splice forms and subcellular localization

Gene duplication is a prominent and recurrent process in plant genomes. Among the possible fates of duplicated genes, subfunctionalization refers to duplicates taking on different parts of the function or expression pattern of the ancestral gene. This partitioning could be accompanied by changes in subcellular localization of the protein products. When alternative splicing of gene products leads to protein products with different subcellular localizations, we propose that after gene duplication there will be partitioning of the alternatively spliced forms such that the products of each duplicate are localized to only one of the original locations, which we refer to as sublocalization. We identified the plastid ascorbate peroxidase (cpAPX) genes across angiosperms and analyzed their duplication history, alternative splicing, and subcellular targeting patterns to identify cases of sublocalization. We found angiosperms typically have one cpAPX gene that generates both thylakoidal APX (tAPX) and stromal APX (sAPX) through alternative splicing. We identified several independent lineage-specific sublocalization cases with specialized paralogues of tAPX and sAPX. We determined that the sublocalization happened through two types of sequence evolution patterns. Our findings suggest that the divergence through sublocalization is key to the retention of paralogous cpAPX genes in angiosperms.

Location, location! cellular relocalization primes specialized metabolic diversification

Specialized metabolites are structurally diverse and cell- or tissue-specific molecules produced in restricted plant lineages. In contrast, primary metabolic pathways are highly conserved in plants and produce metabolites essential for all of life, such as amino acids and nucleotides. Substrate promiscuity - the capacity to accept non-native substrates - is a common characteristic of enzymes, and its impact is especially apparent in generating specialized metabolite variation. However, promiscuity only leads to metabolic diversity when alternative substrates are available; thus, enzyme cellular and subcellular localization directly influence chemical phenotypes. We review a variety of mechanisms that modulate substrate availability for promiscuous plant enzymes. We focus on examples where evolution led to modification of the 'cellular context' through changes in cell-type expression, subcellular relocalization, pathway sequestration, and cellular mixing via tissue damage. These varied mechanisms contributed to the emergence of structurally diverse plant specialized metabolites and inform future metabolic engineering approaches.

Comprehensive analysis of Lon proteases in plants highlights independent gene duplication events

The degradation of damaged proteins is essential for cell viability. Lon is a highly conserved ATP-dependent serine-lysine protease that maintains proteostasis. We performed a comparative genome-wide analysis to determine the evolutionary history of Lon proteases. Prokaryotes and unicellular eukaryotes retained a single Lon copy, whereas multicellular eukaryotes acquired a peroxisomal copy, in addition to the mitochondrial gene, to sustain the evolution of higher order organ structures. Land plants developed small Lon gene families. Despite the Lon2 peroxisomal paralog, Lon genes triplicated in the Arabidopsis lineage through sequential evolutionary events including whole-genome and tandem duplications. The retention of Lon1, Lon4, and Lon3 triplicates relied on their differential and even contrasting expression patterns, distinct subcellular targeting mechanisms, and functional divergence. Lon1 seems similar to the pre-duplication ancestral gene unit, whereas the duplication of Lon3 and Lon4 is evolutionarily recent. In the wider context of plant evolution, papaya is the only genome with a single ancestral Lon1-type gene. The evolutionary trend among plants is to acquire Lon copies with ambiguous pre-sequences for dual-targeting to mitochondria and chloroplasts, and a substrate recognition domain that deviates from the ancestral Lon1 type. Lon genes constitute a paradigm of dynamic evolution contributing to understanding the functional fate of gene duplicates.

Landscape of gene transposition-duplication within the Brassicaceae family

We developed the CLfinder-OrthNet pipeline that detects co-linearity among multiple closely related genomes, finds orthologous gene groups, and encodes the evolutionary history of each orthologue group into a representative network (OrthNet). Using a search based on network topology, we identified 1,394 OrthNets that included gene transposition-duplication (tr-d) events, out of 17,432 identified in six Brassicaceae genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. These were derived from older tr-d events and had a higher chance of being expressed. We also found an enrichment of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3 that relate to their adaptive evolution. CLfinder-OrthNet provides a flexible toolkit to compare gene order, visualize evolutionary paths among orthologues as networks, and identify gene loci that share an evolutionary history.

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.

Atypical Protein Phosphatase 2A Gene Families Do Not Expand via Paleopolyploidization

Protein phosphatase 2A (PP2A) presents unique opportunities for analyzing molecular mechanisms of functional divergence between gene family members. The canonical PP2A holoenzyme regulates multiple eukaryotic signaling pathways by dephosphorylating target proteins and contains a catalytic (C) subunit, a structural/scaffolding (A) subunit, and a regulatory (B) subunit. Genes encoding PP2A subunits have expanded into multigene families in both flowering plants and mammals, and the extent to which different isoform functions may overlap is not clearly understood. To gain insight into the diversification of PP2A subunits, we used phylogenetic analyses to reconstruct the evolutionary histories of PP2A gene families in Arabidopsis (Arabidopsis thaliana). Genes encoding PP2A subunits in mammals represent ancient lineages that expanded early in vertebrate evolution, while flowering plant PP2A subunit lineages evolved much more recently. Despite this temporal difference, our data indicate that the expansion of PP2A subunit gene families in both flowering plants and animals was driven by whole-genome duplications followed by nonrandom gene loss. Selection analysis suggests that the expansion of one B subunit gene family (B56/PPP2R5) was driven by functional diversification rather than by the maintenance of gene dosage. We also observed reduced expansion rates in three distinct B subunit subclades. One of these subclades plays a highly conserved role in cell division, while the distribution of a second subclade suggests a specialized function in supporting beneficial microbial associations. Thus, while whole-genome duplications have driven the expansion and diversification of most PP2A gene families, members of functionally specialized subclades quickly revert to singleton status after duplication events.