FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Feature architecture aware phylogenetic profiling indicates a functional diversification of type IVa pili in the nosocomial pathogen Acinetobacter baumannii.

The Gram-negative bacterial pathogen Acinetobacter baumannii is a major cause of hospital-acquired opportunistic infections. The increasing spread of pan-drug resistant strains makes A. baumannii top-ranking among the ESKAPE pathogens for which novel routes of treatment are urgently needed. Comparative genomics approaches have successfully identified genetic changes coinciding with the emergence of pathogenicity in Acinetobacter. Genes that are prevalent both in pathogenic and a-pathogenic Acinetobacter species were not considered ignoring that virulence factors may emerge by the modification of evolutionarily old and widespread proteins. Here, we increased the resolution of comparative genomics analyses to also include lineage-specific changes in protein feature architectures. Using type IVa pili (T4aP) as an example, we show that three pilus components, among them the pilus tip adhesin ComC, vary in their Pfam domain annotation within the genus Acinetobacter. In most pathogenic Acinetobacter isolates, ComC displays a von Willebrand Factor type A domain harboring a finger-like protrusion, and we provide experimental evidence that this finger conveys virulence-related functions in A. baumannii. All three genes are part of an evolutionary cassette, which has been replaced at least twice during A. baumannii diversification. The resulting strain-specific differences in T4aP layout suggests differences in the way how individual strains interact with their host. Our study underpins the hypothesis that A. baumannii uses T4aP for host infection as it was shown previously for other pathogens. It also indicates that many more functional complexes may exist whose precise functions have been adjusted by modifying individual components on the domain level.

ncOrtho: efficient and reliable identification of miRNA orthologs.

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

Evolutionarily stable gene clusters shed light on the common grounds of pathogenicity in the Acinetobacter calcoaceticus-baumannii complex.

Nosocomial pathogens of the Acinetobacter calcoaceticus-baumannii (ACB) complex are a cautionary example for the world-wide spread of multi- and pan-drug resistant bacteria. Aiding the urgent demand for novel therapeutic targets, comparative genomics studies between pathogens and their apathogenic relatives shed light on the genetic basis of human-pathogen interaction. Yet, existing studies are limited in taxonomic scope, sensing of the phylogenetic signal, and resolution by largely analyzing genes independent of their organization in functional gene clusters. Here, we explored more than 3,000 Acinetobacter genomes in a phylogenomic framework integrating orthology-based phylogenetic profiling and microsynteny conservation analyses. We delineate gene clusters in the type strain A. baumannii ATCC 19606 whose evolutionary conservation indicates a functional integration of the subsumed genes. These evolutionarily stable gene clusters (ESGCs) reveal metabolic pathways, transcriptional regulators residing next to their targets but also tie together sub-clusters with distinct functions to form higher-order functional modules. We shortlisted 150 ESGCs that either co-emerged with the pathogenic ACB clade or are preferentially found therein. They provide a high-resolution picture of genetic and functional changes that coincide with the manifestation of the pathogenic phenotype in the ACB clade. Key innovations are the remodeling of the regulatory-effector cascade connecting LuxR/LuxI quorum sensing via an intermediate messenger to biofilm formation, the extension of micronutrient scavenging systems, and the increase of metabolic flexibility by exploiting carbon sources that are provided by the human host. We could show experimentally that only members of the ACB clade use kynurenine as a sole carbon and energy source, a substance produced by humans to fine-tune the antimicrobial innate immune response. In summary, this study provides a rich and unbiased set of novel testable hypotheses on how pathogenic Acinetobacter interact with and ultimately infect their human host. It is a comprehensive resource for future research into novel therapeutic strategies.

FAS: assessing the similarity between proteins using multi-layered feature architectures.

MOTIVATION: Protein sequence comparison is a fundamental element in the bioinformatics toolkit. When sequences are annotated with features such as functional domains, transmembrane domains, low complexity regions or secondary structure elements, the resulting feature architectures allow better informed comparisons. However, many existing schemes for scoring architecture similarities cannot cope with features arising from multiple annotation sources. Those that do fall short in the resolution of overlapping and redundant feature annotations. RESULTS: Here, we introduce FAS, a scoring method that integrates features from multiple annotation sources in a directed acyclic architecture graph. Redundancies are resolved as part of the architecture comparison by finding the paths through the graphs that maximize the pair-wise architecture similarity. In a large-scale evaluation on more than 10 000 human-yeast ortholog pairs, architecture similarities assessed with FAS are consistently more plausible than those obtained using e-values to resolve overlaps or leaving overlaps unresolved. Three case studies demonstrate the utility of FAS on architecture comparison tasks: benchmarking of orthology assignment software, identification of functionally diverged orthologs, and diagnosing protein architecture changes stemming from faulty gene predictions. With the help of FAS, feature architecture comparisons can now be routinely integrated into these and many other applications. AVAILABILITY AND IMPLEMENTATION: FAS is available as python package: https://pypi.org/project/greedyFAS/.

Genomic evidence for the widespread presence of GH45 cellulases among soil invertebrates.

Lignocellulose is a major component of vascular plant biomass. Its decomposition is crucial for the terrestrial carbon cycle. Microorganisms are considered primary decomposers, but evidence increases that some invertebrates may also decompose lignocellulose. We investigated the taxonomic distribution and evolutionary origins of GH45 hydrolases, important enzymes for the decomposition of cellulose and hemicellulose, in a collection of soil invertebrate genomes. We found that these genes are common in springtails and oribatid mites. Phylogenetic analysis revealed that cellulase genes were acquired early in the evolutionary history of these groups. Domain architectures and predicted 3D enzyme structures indicate that these cellulases are functional. Patterns of presence and absence of these genes across different lineages prompt further investigation into their evolutionary and ecological benefits. The ubiquity of cellulase genes suggests that soil invertebrates may play a role in lignocellulose decomposition, independently or in synergy with microorganisms. Understanding the ecological and evolutionary implications might be crucial for understanding soil food webs and the carbon cycle.

The structure of the Aquifex aeolicus MATE family multidrug resistance transporter and sequence comparisons suggest the existence of a new subfamily.

Multidrug and toxic compound extrusion (MATE) transporters are widespread in all domains of life. Bacterial MATE transporters confer multidrug resistance by utilizing an electrochemical gradient of H+ or Na+ to export xenobiotics across the membrane. Despite the availability of X-ray structures of several MATE transporters, a detailed understanding of the transport mechanism has remained elusive. Here we report the crystal structure of a MATE transporter from Aquifex aeolicus at 2.0-Å resolution. In light of its phylogenetic placement outside of the diversity of hitherto-described MATE transporters and the lack of conserved acidic residues, this protein may represent a subfamily of prokaryotic MATE transporters, which was proven by phylogenetic analysis. Furthermore, the crystal structure and substrate docking results indicate that the substrate binding site is located in the N bundle. The importance of residues surrounding this binding site was demonstrated by structure-based site-directed mutagenesis. We suggest that Aq_128 is functionally similar but structurally diverse from DinF subfamily transporters. Our results provide structural insights into the MATE transporter, which further advances our global understanding of this important transporter family.

Tek (Tie2) is not required for cardiovascular development in zebrafish.

Angiopoietin/TIE signalling plays a major role in blood and lymphatic vessel development. In mouse, Tek (previously known as Tie2) mutants die prenatally due to a severely underdeveloped cardiovascular system. In contrast, in zebrafish, previous studies have reported that although embryos injected with tek morpholinos (MOs) exhibit severe vascular defects, tek mutants display no obvious vascular malformations. To further investigate the function of zebrafish Tek, we generated a panel of loss-of-function tek mutants, including RNA-less alleles, an allele lacking the MO-binding site, an in-frame deletion allele and a premature termination codon-containing allele. Our data show that all these mutants survive to adulthood with no obvious cardiovascular defects. MO injections into tek mutants lacking the MO-binding site or the entire tek locus cause similar vascular defects to those observed in MO-injected +/+ siblings, indicating off-target effects of the MOs. Surprisingly, comprehensive phylogenetic profiling and synteny analyses reveal that Tek was lost in the largest teleost clade, suggesting a lineage-specific shift in the function of TEK during vertebrate evolution. Altogether, these data show that Tek is dispensable for zebrafish development, and probably dispensable in most teleost species.

Minimized combinatorial CRISPR screens identify genetic interactions in autophagy.

Combinatorial CRISPR-Cas screens have advanced the mapping of genetic interactions, but their experimental scale limits the number of targetable gene combinations. Here, we describe 3Cs multiplexing, a rapid and scalable method to generate highly diverse and uniformly distributed combinatorial CRISPR libraries. We demonstrate that the library distribution skew is the critical determinant of its required screening coverage. By circumventing iterative cloning of PCR-amplified oligonucleotides, 3Cs multiplexing facilitates the generation of combinatorial CRISPR libraries with low distribution skews. We show that combinatorial 3Cs libraries can be screened with minimal coverages, reducing associated efforts and costs at least 10-fold. We apply a 3Cs multiplexing library targeting 12,736 autophagy gene combinations with 247,032 paired gRNAs in viability and reporter-based enrichment screens. In the viability screen, we identify, among others, the synthetic lethal WDR45B-PIK3R4 and the proliferation-enhancing ATG7-KEAP1 genetic interactions. In the reporter-based screen, we identify over 1,570 essential genetic interactions for autophagy flux, including interactions among paralogous genes, namely ATG2A-ATG2B, GABARAP-MAP1LC3B and GABARAP-GABARAPL2. However, we only observe few genetic interactions within paralogous gene families of more than two members, indicating functional compensation between them. This work establishes 3Cs multiplexing as a platform for genetic interaction screens at scale.

Ten Years of Collaborative Progress in the Quest for Orthologs.

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

Advances and Applications in the Quest for Orthologs.

Gene families evolve by the processes of speciation (creating orthologs), gene duplication (paralogs), and horizontal gene transfer (xenologs), in addition to sequence divergence and gene loss. Orthologs in particular play an essential role in comparative genomics and phylogenomic analyses. With the continued sequencing of organisms across the tree of life, the data are available to reconstruct the unique evolutionary histories of tens of thousands of gene families. Accurate reconstruction of these histories, however, is a challenging computational problem, and the focus of the Quest for Orthologs Consortium. We review the recent advances and outstanding challenges in this field, as revealed at a symposium and meeting held at the University of Southern California in 2017. Key advances have been made both at the level of orthology algorithm development and with respect to coordination across the community of algorithm developers and orthology end-users. Applications spanned a broad range, including gene function prediction, phylostratigraphy, genome evolution, and phylogenomics. The meetings highlighted the increasing use of meta-analyses integrating results from multiple different algorithms, and discussed ongoing challenges in orthology inference as well as the next steps toward improvement and integration of orthology resources.

Adjustment of photosynthetic activity to drought and fluctuating light in wheat.

Drought is a major cause of losses in crop yield. Under field conditions, plants exposed to drought are usually also experiencing rapid changes in light intensity. Accordingly, plants need to acclimate to both, drought and light stress. Two crucial mechanisms in plant acclimation to changes in light conditions comprise thylakoid protein phosphorylation and dissipation of light energy as heat by non-photochemical quenching (NPQ). Here, we analyzed the acclimation efficacy of two different wheat varieties, by applying fluctuating light for analysis of plants, which had been subjected to a slowly developing drought stress as it usually occurs in the field. This novel approach allowed us to distinguish four drought phases, which are critical for grain yield, and to discover acclimatory responses which are independent of photodamage. In short-term, under fluctuating light, the slowdown of NPQ relaxation adjusts the photosynthetic activity to the reduced metabolic capacity. In long-term, the photosynthetic machinery acquires a drought-specific configuration by changing the PSII-LHCII phosphorylation pattern together with protein stoichiometry. Therefore, the fine-tuning of NPQ relaxation and PSII-LHCII phosphorylation pattern represent promising traits for future crop breeding strategies.

Macrophage HIF-2α regulates tumor-suppressive Spint1 in the tumor microenvironment.

In solid tumors, tumor-associated macrophages (TAMs) commonly accumulate within hypoxic areas. Adaptations to such environments evoke transcriptional changes by the hypoxia-inducible factors (HIFs). While HIF-1α is ubiquitously expressed, HIF-2α appears tissue-specific with consequences of HIF-2α expression in TAMs only being poorly characterized. An E0771 allograft breast tumor model revealed faster tumor growth in myeloid HIF-2α knockout (HIF-2αLysM-/- ) compared with wildtype (wt) mice. In an RNA-sequencing approach of FACS sorted wt and HIF-2α LysM-/- TAMs, serine protease inhibitor, Kunitz type-1 ( Spint1) emerged as a promising candidate for HIF-2α-dependent regulation. We validated reduced Spint1 messenger RNA expression and concomitant Spint1 protein secretion under hypoxia in HIF-2α-deficient bone marrow-derived macrophages (BMDMs) compared with wt BMDMs. In line with the physiological function of Spint1 as an inhibitor of hepatocyte growth factor (HGF) activation, supernatants of hypoxic HIF-2α knockout BMDMs, not containing Spint1, were able to release proliferative properties of inactive pro-HGF on breast tumor cells. In contrast, hypoxic wt BMDM supernatants containing abundant Spint1 amounts failed to do so. We propose that Spint1 contributes to the tumor-suppressive function of HIF-2α in TAMs in breast tumor development.

Identification of the novel class D ß-lactamase OXA-679 involved in carbapenem resistance in Acinetobacter calcoaceticus.

OBJECTIVES: The aim of this study was to characterize the Acinetobacter calcoaceticus clinical isolate AC_2117 with the novel carbapenem-hydrolysing class D ß-lactamase (CHDL) OXA-679. METHODS: Identification of the species and ß-lactamases was verified by genome sequencing (PacBio) and phylogenetic analyses. Antibiotic susceptibility of AC_2117 and transformants harbouring cloned blaOXA-679 was evaluated using antibiotic gradient strips and microbroth dilution. OXA-679 was purified heterologously and kinetic parameters were determined using spectrometry or isothermal titration calorimetry. The impact of OXA-679 production during imipenem therapy was evaluated in the Galleria mellonella infection model. RESULTS: Sequencing of the complete genome of the clinical A. calcoaceticus isolate AC_2117 identified a novel CHDL, termed OXA-679. This enzyme shared sequence similarity of 71% to each of the families OXA-143 and OXA-24/40. Phylogenetic analyses revealed that OXA-679 represents a member of a new OXA family. Cloning and expression of blaOXA-679 as well as measurement of kinetic parameters revealed the effective hydrolysis of carbapenems which resulted in reduced susceptibility to carbapenems in Escherichia coli and A. calcoaceticus, and high-level carbapenem resistance in Acinetobacter baumannii. Infection of larvae of G. mellonella with a sublethal dose of blaOXA-679-expressing A. baumannii could not be cured by high-dose imipenem therapy, indicating carbapenem resistance in vivo. CONCLUSIONS: We identified blaOXA-679 in a clinical A. calcoaceticus isolate that represents a member of the new OXA-679 family and that conferred high-level carbapenem resistance in vitro and in vivo.

PhyloProfile: dynamic visualization and exploration of multi-layered phylogenetic profiles.

Summary: Phylogenetic profiles form the basis for tracing proteins and their functions across species and through time. Novel genome sequences nowadays often represent species from the remotest corner of the tree of life. Thus, phylogenetic profiling becomes increasingly important for functionally annotating this data and to integrate it into a comprehensive view on organismal evolution. To strengthen the link between the sharing of a gene across species and of the corresponding function, it is meanwhile common to complement phylogenetic profiles with additional information, such as domain architecture similarities between orthologs, or pairwise similarities of other protein features. However, there are few visualization tools that facilitate an intuitive integration of these various information layers. Here, we present PhyloProfile, an R-based tool to visualize, explore and analyze multi-layered phylogenetic profiles. Availability and implementation: PhyloProfile is available as open source code under the MIT license at https://github.com/BIONF/phyloprofile. An online version for testing PhyloProfile and for small to medium-scale analyses is available at http://applbio.biologie.uni-frankfurt.de/phyloprofile.

Gearing up to handle the mosaic nature of life in the quest for orthologs.

The Quest for Orthologs (QfO) is an open collaboration framework for experts in comparative phylogenomics and related research areas who have an interest in highly accurate orthology predictions and their applications. We here report highlights and discussion points from the QfO meeting 2015 held in Barcelona. Achievements in recent years have established a basis to support developments for improved orthology prediction and to explore new approaches. Central to the QfO effort is proper benchmarking of methods and services, as well as design of standardized datasets and standardized formats to allow sharing and comparison of results. Simultaneously, analysis pipelines have been improved, evaluated and adapted to handle large datasets. All this would not have occurred without the long-term collaboration of Consortium members. Meeting regularly to review and coordinate complementary activities from a broad spectrum of innovative researchers clearly benefits the community. Highlights of the meeting include addressing sources of and legitimacy of disagreements between orthology calls, the context dependency of orthology definitions, special challenges encountered when analyzing very anciently rooted orthologies, orthology in the light of whole-genome duplications, and the concept of orthologous versus paralogous relationships at different levels, including domain-level orthology. Furthermore, particular needs for different applications (e.g. plant genomics, ancient gene families and others) and the infrastructure for making orthology inferences available (e.g. interfaces with model organism databases) were discussed, with several ongoing efforts that are expected to be reported on during the upcoming 2017 QfO meeting.

Recombinant production of A1S_0222 from Acinetobacter baumannii ATCC 17978 and confirmation of its DNA-(adenine N6)-methyltransferase activity.

Acinetobacter baumannii appears as an often multidrug-resistant nosocomial pathogen in hospitals worldwide. Its remarkable persistence in the hospital environment is probably due to intrinsic and acquired resistance to disinfectants and antibiotics, tolerance to desiccation stress, capability to form biofilms, and is possibly facilitated by surface-associated motility. Our attempts to elucidate surface-associated motility in A. baumannii revealed a mutant inactivated in a putative DNA-(adenine N6)-methyltransferase, designated A1S_0222 in strain ATCC 17978. We recombinantly produced A1S_0222 as a glutathione S-transferase (GST) fusion protein and purified it to near homogeneity through a combination of GST affinity chromatography, cation exchange chromatography and PD-10 desalting column. Furthermore we demonstrate A1S_0222-dependent adenine methylation at a GAATTC site. We propose the name AamA (Acinetobacteradenine methyltransferase A) in addition to the formal names M.AbaBGORF222P/M.Aba17978ORF8565P. Small angle X-ray scattering (SAXS) revealed that the protein is monomeric and has an extended and likely two-domain shape in solution.

The association of late-acting snoRNPs with human pre-ribosomal complexes requires the RNA helicase DDX21.

Translation fidelity and efficiency require multiple ribosomal (r)RNA modifications that are mostly mediated by small nucleolar (sno)RNPs during ribosome production. Overlapping basepairing of snoRNAs with pre-rRNAs often necessitates sequential and efficient association and dissociation of the snoRNPs, however, how such hierarchy is established has remained unknown so far. Here, we identify several late-acting snoRNAs that bind pre-40S particles in human cells and show that their association and function in pre-40S complexes is regulated by the RNA helicase DDX21. We map DDX21 crosslinking sites on pre-rRNAs and show their overlap with the basepairing sites of the affected snoRNAs. While DDX21 activity is required for recruitment of the late-acting snoRNAs SNORD56 and SNORD68, earlier snoRNAs are not affected by DDX21 depletion. Together, these observations provide an understanding of the timing and ordered hierarchy of snoRNP action in pre-40S maturation and reveal a novel mode of regulation of snoRNP function by an RNA helicase in human cells.

An evolutionary perspective of AMPK-TOR signaling in the three domains of life.

AMPK and TOR protein kinases are the major control points of energy signaling in eukaryotic cells and organisms. They form the core of a complex regulatory network to co-ordinate metabolic activities in the cytosol with those in the mitochondria and plastids. Despite its relevance, it is still unclear when and how this regulatory pathway was formed during evolution, and to what extent its representations in the major eukaryotic lineages resemble each other. Here we have traced 153 essential proteins forming the human AMPK-TOR pathways across 412 species representing all three domains of life-prokaryotes (bacteria, archaea) and eukaryotes-and reconstructed their evolutionary history. The resulting phylogenetic profiles indicate the presence of primordial core pathways including seven proto-kinases in the last eukaryotic common ancestor. The evolutionary origins of the oldest components of the AMPK pathway, however, extend into the pre-eukaryotic era, and descendants of these ancient proteins can still be found in contemporary prokaryotes. The TOR complex in turn appears as a eukaryotic invention, possibly to aid in retrograde signaling between the mitochondria and the remainder of the cell. Within the eukaryotes, AMPK/TOR showed both a highly conserved core structure and a considerable plasticity. Most notably, KING1, a protein originally assigned as the γ subunit of AMPK in plants, is more closely related to the yeast SDS23 gene family than to the γ subunits in animals or fungi. This suggests its functional difference from a canonical AMPK γ subunit.

An evolutionary view on thylakoid protein phosphorylation uncovers novel phosphorylation hotspots with potential functional implications.

The regulation of photosynthetic light reactions by reversible protein phosphorylation is well established today, but functional studies have so far mostly been restricted to processes affecting light-harvesting complex II and the core proteins of photosystem II. Virtually no functional data are available on regulatory effects at the other photosynthetic complexes despite the identification of multiple phosphorylation sites. Therefore we summarize the available data from 50 published phospho-proteomics studies covering the main complexes involved in photosynthetic light reactions in the 'green lineage' (i.e. green algae and land plants) as well as its cyanobacterial counterparts. In addition, we performed an extensive orthologue search for the major photosynthetic thylakoid proteins in 41 sequenced genomes and generated sequence alignments to survey the phylogenetic distribution of phosphorylation sites and their evolutionary conservation from green algae to higher plants. We observed a number of uncharacterized phosphorylation hotspots at photosystem I and the ATP synthase with potential functional relevance as well as an unexpected divergence of phosphosites. Although technical limitations might account for a number of those differences, we think that many of these phosphosites have important functions. This is particularly important for mono- and dicot plants, where these sites might be involved in regulatory processes such as stress acclimation.