Plasmodium evolutionary rate covariation matrix
Data files
Jul 07, 2026 version files 114.35 MB
-
Plasmodium_ftERC.txt
98.30 MB
-
Plasmodium_OGs_singlecopy.xlsx
858.28 KB
-
Plasmodium_OGs.xlsx
15.18 MB
-
README.md
4.14 KB
-
Species_name_key.xlsx
9.48 KB
Jul 22, 2026 version files 114.40 MB
-
Plasmodium_ftERC.txt
98.36 MB
-
Plasmodium_OGs_singlecopy.xlsx
858.28 KB
-
Plasmodium_OGs.xlsx
15.18 MB
-
README.md
4.14 KB
-
Species_name_key.xlsx
9.48 KB
Abstract
Despite publication of the Plasmodium falciparum reference genome over 20 years ago, one-third of its genes remain functionally unannotated and information is limited for many others. Proteins that act in the same pathway or complex tend to experience similar shifts in evolutionary pressure, such that correlated constraints across species can indicate co-functional proteins. To investigate this connection, we calculated relative evolutionary rates across the genome on a phylogeny of 22 Plasmodium species and assigned each protein-protein pair a score representing the strength of evolutionary rate covariation (ERC). We show that known pathways and interacting proteins in Plasmodium have strong ERC, including proteins involved in red blood cell invasion and mosquito-specific pathways. By scanning genome-wide for additional proteins showing high ERC with established interacting proteins, we find enrichment of stage expression and physically interacting protein pairs supporting new candidate functions for proteins. More generally, we demonstrate the utility of ERC to prioritize proteins for hypothesis-driven functional follow up by showing that a protein with little functional characterization (PF3D7_0811600) shows high ERC and co-localizes with high molecular weight rhoptry proteins 2 and 3 (RhopH2 and RhopH3), which form an ion channel that enhances parasite permeability of infected red blood cells. The ERC matrix can be queried to extract Plasmodium proteins showing high ERC with any protein of interest to prioritize candidate genes and accelerate discovery of novel functional connections.
Dataset DOI: 10.5061/dryad.cfxpnvxmw
Description of the data and file structure
This dataset contains files for the analyses in the manuscript titled "Evolutionary rate covariation across malaria parasite species enables inference of protein interactions" These files include the estimated ortholog groups used to calculate the evolutionary rate covariation (ERC) matrix.
We assigned groups of paralogous and orthologous proteins across protein sequences from 22 Plasmodium species. Each group is named with the prefix "clusters" followed by a number (e.g., "clusters1000"). Missing values are denoted by "N/A".
Plasmodium_OGs.xlsx (single sheet): Each of the 8,085 rows represents one group of paralogous and orthologous proteins (column "V1"). Columns V2 through V714 represent individual protein members of each group, with cell values formatted as species name abbreviation | gene ID (for example, mal|PmUG01_00010100). This range corresponds to the largest group size in the dataset (713 proteins). Groups with fewer than 713 member proteins have "N/A" in the remaining columns. Group sizes vary because not all species are represented in every group, and some groups contain multiple proteins from a single species.
Plasmodium_OGs_singlecopy.xlsx (single sheet): Ortholog groups excluding paralogous genes and groups with fewer than 10 single copy orthologs, resulting in 4,410 groups (rows). The first column ("clusterID") indicates the group, and the 22 remaining columns are labeled by species (see Species_name_key.xlsx for abbreviations), each containing the single copy ortholog gene ID for that species, where one exists. "N/A" indicates species that do not have a single copy ortholog for a particular group.
Plasmodium_ftERC.txt (tab-delimited text file): Full genome evolutionary rate covariation (ERC) matrix containing ERC values across 4,360 genes, resulting in a 4,360 x 4,360 matrix. ERC values represent the correlation between two genes' relative evolutionary rates, which reflect how much faster or slower a gene is evolving in one species relative to all species across the tree. ERC values are Fisher transformed correlation coefficients (which corrects for number of datapoints and makes values more comparable across genes) [1] and thus are unitless and range from -7.2 to 8.9. ERC was calculated from multiple sequence alignments of proteins in the Plasmodium_OGs_singlecopy.xlsx file. Each row and column is named by cluster ID. The upper triangle and diagonal are set to "N/A". Groups in the lower triangle for which an ERC value could not be calculated are also set to "N/A".
Species_name_key.xlsx (single sheet): Key for species name abbreviations used in files in this dataset and reference genomes used to create the dataset. Column "Species abbreviation in Plasmodium_OGs": species abbreviation as used in Plasmodium_OGs.xlsx (e.g., berg) Column "Species abbreviation in Plasmodium_OGs_singlecopy": corresponding abbreviation used as column headers in Plasmodium_OGs_singlecopy.xlsx (e.g., berghei).
Species_name_key.xlsx (single sheet): Key linking species abbreviations used across files in this dataset to full species and strain names.
- Column "Species abbreviation in Plasmodium_OGs": species abbreviation as used in Plasmodium_OGs.xlsx (e.g., "berg").
- Column "Species abbreviation in Plasmodium_OGs_singlecopy": corresponding abbreviation used as column headers in Plasmodium_OGs_singlecopy.xlsx (e.g., "berghei").
- Column "Species reference genome": full species and strain name corresponding to each abbreviation (e.g., "Plasmodium berghei ANKA").
Sharing/Access information
We used the amino acid sequences from 22 Plasmodium species reference genomes from PlasmoDB release 54/September 2021 (https://plasmodb.org/plasmo/app/downloads).
For more information on ERC:
Little et al. 2025 ERC2.0 evolutionary rate covariation update improves inference of functional interactions across large phylogenies. Cold Spring Harbor Laboratory Press. DOI: 10.1101/gr.280586.125
Changes after Jul 7, 2026: Added explicit row-name labels to Plasmodium_ftERC.txt to match the existing column names. Data values are unchanged.
