Phylogenomic analyses reveal introgression and cryptic speciation in the globally distributed, vector-transmitted pathogen Plasmodium relictum
Data files
May 04, 2026 version files 244.10 KB
-
a-read_quality.slurm
612 B
-
A01.ctrl
2.07 KB
-
A10_grw4_s_3.ctrl
2.32 KB
-
A10.ctrl
2.27 KB
-
b-remove_adaptors_quality_trim.sh
1.01 KB
-
c-map_to_reference.sh
813 B
-
codeml_model_selection.R
3.51 KB
-
d-make_consensus.sh
1.03 KB
-
e-prepare_probes.slurm
1.54 KB
-
extra_intro.csv
5.26 KB
-
extra_mig.csv
5.24 KB
-
f-get_gene.sh
3.14 KB
-
g-get_orthologues.sh
835 B
-
h_samples.txt
70 B
-
h-align_orthologues.sh
385 B
-
helpers.sh
2.06 KB
-
i-gene_tree_inference.sh
1.31 KB
-
j-generate_species_tree.sh
2.32 KB
-
k-filter_orthologue_output.sh
1.55 KB
-
l-paml_prep.sh
3.06 KB
-
m-bpp_prep.sh
935 B
-
model17.ctl
2.04 KB
-
model18.ctl
2.21 KB
-
model19.ctl
2.04 KB
-
model20.ctl
2.21 KB
-
model21.ctl
2.04 KB
-
model22.ctl
2.21 KB
-
model23.ctl
2.18 KB
-
model24.ctl
2.18 KB
-
model25.ctl
2.18 KB
-
n-bpp.sh
333 B
-
nh_samples.txt
403 B
-
o-bpp_ht.sh
1.48 KB
-
orthogroup_subset.txt
15.41 KB
-
plasphylo_figures.R
16.29 KB
-
q-branch_prep.sh
1.69 KB
-
r-codeml_loop.sh
3.56 KB
-
README.md
15.18 KB
-
s-stat_loop.sh
2.32 KB
-
sample_lineage_location.csv
700 B
-
SpeciesTreeParameterExplorationUpdated.txt
92.46 KB
-
strains_3.txt
821 B
-
strains.txt
812 B
-
sub_D00_loop.sh
20.71 KB
-
sub_D00_template_introgression.ctl
1.90 KB
-
sub_D00_template_migration.ctl
1.90 KB
-
t-param_loop.sh
846 B
-
template_M0.ctl
1.01 KB
-
template_M1.ctl
1.02 KB
-
template_M2.ctl
1.01 KB
-
template_M3.ctl
1.01 KB
-
template_M4.ctl
1.01 KB
-
test_pipeline.slurm
1.55 KB
Abstract
Establishing species limits is challenging, particularly for pathogens of wildlife. These pathogens can be difficult to sample and culture, and their genome sequencing must often be conducted in the presence of high levels of host DNA. Plasmodium relictum is a mosquito-vectored avian malaria pathogen that is a globally distributed host generalist, comprised of several genetic lineages. We used sequence capture data from 52 P. relictum infections originating from multiple continents to generate a genomic dataset of the pathogen. With this data, we established a robust phylogeny and determined species limits among P. relictum lineages. We generated phylogenomic trees by maximum likelihood and Bayesian methods with multi-species coalescent models and confirmed robustness of the topology by varying the amount of missing data in the analyses. Our results suggest the existence of two cryptic species among the infections we analyzed and provide evidence of genetic introgression between these species. One of the cryptic species, GRW4, devastated the endemic and immunologically naïve avifauna of Hawaii after its introduction to the islands ca. 100 years ago, and so was tested for positive selection in the GRW4 Hawaiian clade. Although we hypothesized it would be released from host selective pressures, we did not find evidence of positive selection in the Hawaiian GRW4 clade, and we discuss possible explanations. Overall, our results underscore the importance of genomic analyses for resolving pathogen species limits and understanding pathogen evolution.
Dataset DOI: 10.5061/dryad.jq2bvq8r8
Description of the data and file structure
We used sequence capture data from 52 P. relictum infections originating from multiple continents to generate a genomic dataset of the pathogen.
Files and variables
File: c-map_to_reference.sh
Description: Maps trimmed reads from the output of b-remove_adaptors_quality_trim.sh to the P. relictum reference genome GCA_900005765.1 and then records mapping statistics.
File: d-make_consensus.sh
Description: Generates a consensus genome for each sample from the c-map_to_reference.sh output.
File: extra_intro.csv
Description: Summary of output from additional introgression tests performed on groups from different chromosomal regions (orthogroup_subsets.txt) in sub_D00_loop.sh script.
Variables
- sample: Chromosomal region label used to manage samples in analysis. Corresponds to orthogroups in orthogroup_subsets.txt.
- label: Chromosome of each sample for color grouping in plot.
- location: Position for plotting result along genome; integer value corresponds to the chromosome and then the results from the first, second, and and last third of the chromosome assigned to +0.25, +0.5, and +0.75.
- phy: Probability of gene flow.
- lower: Lower end of 95 % HPD (highest posterior density) confidence intervals from MSC-I model.
- upper: Upper end of 95 % HPD (highest posterior density) confidence intervals from MSC-I model.
- count: Number of aligned sequences in chromosomal region used for analysis
- origin: Direction of gene flow, entry = x indicates flow in the direction from point X (GRW4 branch) to point Y (DONANA05/SGS1/GRW11 branch).
- avl: Average length of sequences in samples for analysis.
File: i-gene_tree_inference.sh
Description: Takes alignments from h-align_orthologues.sh that contain at least 4 sequences and checks with RAxML to remove any duplicate sequences or other errors. Then calculates the maximum likelihood best tree from 200 bootstrap replicates.
File: a-read_quality.slurm
Description: Runs FastQC and then MultiQC for quality assessment of the original sequence files, NCBI BioProject (Accession: PRJNA1330847).
File: e-prepare_probes.slurm
Description: Maps the original sequencing probes, (Ellis et al., 2022) Additional file 3: File S1, to both reference genomes (GCA_900005765.1, GCA_900005855.1) and generates two bed files of probe coordinates, one for the P.relictum reference and samples, and one for the P.gallinaceum reference.
File: f-get_gene.sh
Description: Applies the coordinates in the bed file from e-prepare_probes.slurm to generate a FASTA file of sequences within those coordinates for each sample from the consensus genome output from d-make_consensus.sh as well as for both reference genomes (GCA_900005765.1, GCA_900005855.1). Sequences are labeled with sample name and coordinates, then filtered based on the allowed percentage ambiguity in helpers.sh. Samples with no remaining sequences are removed.
File: g-get_orthologues.sh
Description: Applies Orthofinder to generate groups of orthologues. The groups are then filtered so that they contain only a single representative sequence from each sample. Groups that contain only a single sequence are removed.
File: o-bpp_ht.sh
Description: Generates the nested directories for 8-point calculation of each of the 9 models: 3 without introgression (model17.ctl, model19.ctl, model21.ctl), 3 with introgression (model18.ctl, model20.ctl, model22.ctl), and 3 with migration (model23.ctl, model24.ctl, model25.ctl). Copies the main template file for each model into the head folder, then uses BPP to generate control files for each of the 8 beta-weights for each model. Then runs BPP over all of the directories and control files (total 48). It refers to strains.txt for assigning the lineage names (except GALLUS01, marked as “Outgroup”) to the samples.
File: r-codeml_loop.sh
Description: Contains all of the loops for generating and then running CODEML control files for analysis, including checking for the best reading frame for the non-branch models. First two loops for base (template_M0.ctl) and site model (template_M1.ctl) and then for branch (template_M3.ctl) and branch-site estimate omega (template_M3.ctl) and branch-site with fixed omega for comparison (template_M4.ctl).
File: q-branch_prep.sh
Description: Prepares files for CODEML branch-site analysis by first calculating all of the stop codons (recorded as ???) in the three reading frame sequences alignment files, and the proceeding with the alignment with the fewest stop codons. Trees and alignments are filtered to contain at least one Hawaiian sample (h_samples.txt) and one non-Hawaiian sample (nh_samples.txt) and the occurrences of Hawaiian samples in the gene trees are labeled as the forward branch.
File: template_M1.ctl
Description: CODEML file. A control file for calculating several site models.
File: plasphylo_figures.R
Description: R-script that makes the figures, combined into a single script for simplicity.
File: sub_D00_loop.sh
Description: Bash script for preparing the control files and running the additional introgression and migration models based on orthogroups organized by coordinates described in orthogroup_subset.txt.
File: sample_lineage_location.csv
Description: A key for colors and shapes of the samples to the lineages with the location names for plot in Figure 4 and Supplementary Figure 1 in the associated manuscript.
Variables
- Sample_ID: Unique sample id corresponding to original sequence files NCBI BioProject (Accession: PRJNA1330847).
- Lineage: Sample lineage identified by mitochondrial gene cytochrome b.
File: A01.ctrl
Description: BPP script. A control file for species tree estimation, associated with the result presented in Figure 1, best trees in the MCMC sampling.
File: extra_mig.csv
Description: Summary of output from additional migration tests performed on groups from different chromosomal regions (orthogroup_subsets.txt) in sub_D00_loop.sh script.
Variables
- sample: Chromosomal region label used to manage samples in analysis. Corresponds to orthogroups in orthogroup_subsets.txt.
- label: Chromosome of each sample for color grouping in plot.
- location: Position for plotting result along genome; integer value corresponds to the chromosome and then the results from the first, second, and and last third of the chromosome assigned to +0.25, +0.5, and +0.75.
- M: mean migration rate.
- lower: Lower end of error bar for plotting. Calculated by subtracting the standard deviation from the mean.
- upper: Upper end of error bar for plotting. Calculated by adding the standard deviation to the mean.
- count: Number of aligned sequences in chromosomal region used for analysis
- origin: Direction of gene flow, entry = x indicates flow in the direction from point X (GRW4 branch) to point Y (DONANA05/SGS1/GRW11 branch).
- SD: Standard deviation of the migration rate.
File: codeml_model_selection.R
Description: Tests between likelihoods of models containing both Hawaiian and non-Hawaiian GRW4.
File: nh_samples.txt
Description: List of non-Hawaiian samples, including both reference genomes.
Variables
- List of samples corresponding to original sequence files NCBI BioProject (Accession: PRJNA1330847).
File: h_samples.txt
Description: List of Hawaiian samples
Variables
- List of sample IDs corresponding to original sequence files NCBI BioProject (Accession: PRJNA1330847)
File: h-align_orthologues.sh
Description: Applies mafft, set to global alignment, to the groups of orthologues output from g-get_orthologues.sh.
File: A10.ctrl
Description: BPP script. A control file for species delimitation analysis, associated with the result presented in Figure 2.
File: strains_3.txt
Description: Tab-separated list of sample IDs and corresponding lineage for BPP analysis, with GRW4 split geographically: GRW4_H for Hawaiian samples, GRW4_US for continental US samples, GRW4_B for Bulgaria sample.
Variables
- First column (no header): sample ID corresponding to original sequence files NCBI BioProject (Accession: PRJNA1330847).
- Second column (no header): sample lineage identified by mitochondrial gene cytochrome b, further separated by geography for GRW4.
File: template_M4.ctl
Description: CODEML file. A control file for branch-site model with fixed omega.
File: SpeciesTreeParameterExplorationUpdated.txt
Description: Collection of ASTRAL sample tree outputs from different parameter selections.
Variables
- Allowed_N: Maximum percentage of allowed N bases in sequence.
- Gene_Tree_Representation: Minimum required percent representation of sample in gene trees.
- Species_Tree: ASTRAL output tree with lengths for non-terminal branches.
- Species_Tree_LPP: ASTRAL output tree with local posterior probability for nodes.
- Species_Tree_QS: ASTRAL output tree with quartet support for nodes.
File: model19.ctl
Description: BPP file. A control file for MCS-I model, random seed = 40, associated with the results presented in Table 2 and Table 3 in the associated manuscript.
File: model20.ctl
Description: BPP file. A control file for MCS-I model, random seed = 40, associated with the results presented in Table 2 and Table 3.
File: m-bpp_prep.sh
Description: Takes aligned sequences from i-gene_tree_inference.sh and converts to PHYLIP format for BPP input. Then all of the aligned sequences are concatenated into a single file.
File: model21.ctl
Description: BPP file. A control file for non-introgression model, random seed = 50, associated with the results presented in Table 2 and Table 3.
File: model23.ctl
Description: BPP file. A control file for MCS-M model, random seed = 35, associated with the results presented in Table 2 and Table 3.
File: orthogroup_subset.txt
Description: List of Orthogroups sorted by chromosomal positions, corresponding to extra_intro.csv file.
Variables
- List of samples following chromosome labels: ch00 (99/00/MIT), chr01_1 (up t0 240,639), chr01_2 (between 240,639 and 481,278), etc.
File: helpers.sh
Description: Contains directory information for intermediate pipeline files and installed software. Includes variable information such as number of threads for some software, percentage allowed N for sequences, and minimum gene tree representation for final sample tree calculation.
File: b-remove_adaptors_quality_trim.sh
Description: Trims adaptor and poor-quality sequences from original sequence files, NCBI BioProject (Accession: PRJNA1330847), then calculates remaining quality statistics.
File: model22.ctl
Description: BPP file. A control file for MCS-I model, random seed = 50, associated with the results presented in Table 2 and Table 3.
File: model18.ctl
Description: BPP file. A control file for MCS-I model, random seed = 35, associated with the results presented in Table 2 and Table 3.
File: strains.txt
Description: Tab-separated list of sample IDs and corresponding lineage for BPP analysis
Variables
- First column (no header): sample ID corresponding to original sequence files NCBI BioProject (Accession: PRJNA1330847).
- Second column (no header): sample lineage identified by mitochondrial gene cytochrome b
File: model25.ctl
Description: BPP file. A control file for MCS-M model, random seed = 50, associated with the results presented in Table 2 and Table 3.
File: s-stat_loop.sh
Description: Collects model likelihood statistics from CODEML output files to be tested for significant difference in R.
File: j-generate_species_tree.sh
Description: Combines the sequence trees from the i-gene_tree_inference.sh output and then filters out trees that contain poorly represented samples, based on the helpers.sh input. Then applies ASTRAL to calculate a sample tree from the remaining gene trees.
File: model24.ctl
Description: BPP file. A control file for MCS-M model, random seed = 40, associated with the results presented in Table 2 and Table 3.
File: model17.ctl
Description: BPP file. A control file for non-introgression model, random seed = 35, associated with the results presented in Table 2 and Table 3.
File: sub_D00_template_introgression.ctl
Description: Control and bash file for running the additional introgression models based on orthogroups organized by coordinates described in orthogroup_subset.txt.
File: sub_D00_template_migration.ctl
Description: Control and bash file for running the additional migration models based on orthogroups organized by coordinates described in orthogroup_subset.txt.
File: template_M0.ctl
Description: CODEML file. A control file for homogenous model.
File: template_M3.ctl
Description: CODEML file. A control file for branch-site model.
File: test_pipeline.slurm
Description: Overall controlling file for SLURM batch system. Requests total time and resources for the job, calls each step while tracking the start and finish time for all of them, and emails the user when finished (failed, completed, or timed out).
File: A10_grw4_s_3.ctrl
Description: BPP script. An additional control file for the species delimitation analysis separating the GRW4 samples between Bulgaria, the continental US, and Hawaii.
File: n-bpp.sh
Description: Runs BPP analysis from A01.ctrl control file, associated with Figure 1 result (output A: best trees in the sample) and analysis A10.ctrl control file, associated with Figure 2 result (guide tree with posterior probability for presence of nodes). Also the additional species delimitation analysis separating the GRW4 samples geographically.
File: t-param_loop.sh
Description: Collects model parameters for Supplementary Table 6.
File: template_M2.ctl
Description: CODEML file. A control file for branch model.
File: k-filter_orthologue_output.sh
Description: Filters the Orthofinder output table of sequences per group based on the remaining groups used in the j-generate_species_tree.sh to generate content in Supplementary Table 2.
File: l-paml_prep.sh
Description: Prepares aligned sequence files for CODEML analysis by converting them to phylip format and generates the two other reading frame aligned sequence files. Also takes the gene tree outputs from i-gene_tree_interence.sh, adding the sequence length and number of samples to the top, and removing the branch length information.
Code/software
All code and scripts are included in the submitted files.
