Data for: Natural variation reveals functional and genetic integration of a polyphenism
Data files
Mar 09, 2026 version files 44.60 MB
-
final.all.control.8.19.25.yaml
286 B
-
final.all.geno.8.19.25.csv
16.95 MB
-
final.all.gmap.csv
6.70 MB
-
final.all.ptypes.8.26.25.csv
3.40 KB
-
GM_RILs_all_RENAMED.txt
204.27 KB
-
GM_RILs.8.28.25.rmd
5.10 KB
-
onemap.rils.REAL.raw
20.70 MB
-
original_screening_phenotyping.txt
2.92 KB
-
ppa.cM.map.8.19.25.R
3.88 KB
-
PS312.RS5200B.map.30.samples.500k_8.19.25.txt
10.39 KB
-
qtl.with.cM.8.19.25.R
8.96 KB
-
README.md
5.11 KB
-
RIL.raw.corpses.txt
2.17 KB
-
RIL.raw.morph.ratio.txt
773 B
-
screening.ptypes.8.28.25.R
5.35 KB
Abstract
Integration and modularity can have a profound impact on the function and evolution of environmentally responsive traits, especially when they result in discrete, alternative forms—that is, developmental polyphenism. An unresolved issue for understanding this impact is the degree to which the genetic architectures of the individual components of a plastic trait permit independent versus coordinated evolution. The association of trait variation with genomic variation can provide a test of whether the same loci influence different components of the same integrated phenotype. An example of a coordinated, plastic trait is in the shark-tooth nematode Pristionchus pacificus, which develops into either a bacterial-feeding or a predatory adult morph, depending on its perception of local food availability. Moreover, this polyphenism, when measured as morph induction in response to a common set of cues, differs across natural isolates of the species. By creating recombinant inbred lines (RILs) from natural isolates that have diverged in their morph-induction bias, followed by quantitative trait locus analysis, we tested whether and to what extent component traits of this resource polyphenism are linked. We found that RILs with more frequent induction of the predatory morph also produced Eu individuals that were more effective predators. We also found that these two traits are associated with the same major-effect locus, suggesting that their causal genes are physically linked, if not the same, and are therefore likely to experience coordinated selection. In contrast, we found that morphological variation was not linked to these two traits and that such variation within each morph was even independent of variation in the other. Our findings show that the same coordinated plastic trait exhibits a blend of genetic correlation and independence, whose balance shapes the trait’s evolutionary potential
Dataset DOI: 10.5061/dryad.5qfttdzkm
Description of the data and file structure
These data primarily include phenotypic data--mouth form ratio, predatory ability, and morphometrics--of recombinant inbred lines (RILs) generated by crossing two natural isolates of the shark-tooth nematode Pristionchus pacificus that exhibit divergent ratios of mouth form production. In addition, code and genome-wide SNP calls for each RIL are included as input for trait associations and quantitative trait locus analysis. In most cases, data are organized into tab-delimited data frames easily read into R, but a few files have a specialized structure to accommodate particular packages used in the analyses.
Recombinant inbred lines were generated by crossing two strains of Pristionchus pacificus and then allowing the F1 progeny to self-fertilize until the F12 generation. RILs were then genotyped and phenotyped for quantitative trait locus mapping.
Files and variables
File: GM_RILs.8.28.25.rmd
Description: R markdown file for analyzing geometric morphometric data and generating values for each RIL that are used as input for correlation and QTL analyses
File: PS312.RS5200B.map.30.samples.500k_8.19.25.txt
Description: Map generated by the script ppa.cM.map.8.19.25.R and used as input for the script qtl.with.cM.8.19.25.R
Variables
- Chromosome: The chromosome of the marker
- Marker: The marker ID in the form of Chromosome.base position
- cM.Position: The marker position in cM
File: GM_RILs_all_RENAMED.txt
Description: Landmark data for 10 Eu and St individuals per RIL, with names that can be parsed in the associated R markdown file
Variables
- LM=20 denotes the number of landmarks for that individual
- Next are the x-y coordinates for the landmarks
- ID= denotes the sample ID
File: original_screening_phenotyping.txt
Description:
Variables
- RIL: The recombinant inbred line
- Eu: The number of Eu individuals
- St: The number of St individuals
- %Eu: The number of Eu individuals out of the total number of individuals evaluated for that RIL
File: ppa.cM.map.8.19.25.R
Description: R script for converting genome-wide SNP data to a genetic map with cM values for downstream QTL analysis.
File: RIL.raw.corpses.txt
Description: The number of corpses produced by each of the sequenced RILs.
Variables
- RIL: RIL identity
- Corpse: The number of corpses counted for that replicate
File: qtl.with.cM.8.19.25.R
Description: R script for performing QTL analysis of 10 phenotypic traits across 30 RILs
File: RIL.raw.morph.ratio.txt
Description: The Eu and St values for each sequenced RIL
Variables
- RIL: RIL identity
- Percent Eu: The percent of individuals that were Eu out of 50 individuals evaluated for that replicate
File: screening.ptypes.8.28.25.R
Description: R script for generating histograms and scatterplots of phenotypes for screened and sequenced RILs
File: final.all.control.8.19.25.yaml
Description: YAML format control file needed for QTL analysis
File: final.all.ptypes.8.26.25.csv
Description: Phenotype data for each RIL for QTL analyses
Variables
- id: RIL identity
- Mean. percent.Eu: Mean percent Eu for that RIL
- Median corpses: Median count of corpses for that RIL
- Distance.bn.morphs: Distance between Eu and St morphs according to geomorph trajectory analysis for that RIL
- All. disparity: Total disparity for that RIL based on geomorph analysis
- Eu.sub.disparity: Total disparity of Eu morph for that RIL
- St.sub.disparity: Total disparity of St morph for that RIL
- Eu.sub.PC1: Mean PC1 position of Eu morph for that RIL
- Eu.sub.PC2: Mean PC2 position of Eu morph for that RIL
- St.sub.PC1: Mean PC1 position of St morph for that RIL
- St.sub.PC2: Mean PC2 position of St morph for that RIL
File: onemap.rils.REAL.raw
Description: Genotype calls in 'raw' format used by the onemap package in R for generating a genetic map
File: final.all.gmap.csv
Description: The full set of markers that were genotyped. Used as input for QTL analysis.
Variables
- marker: The name of the marker
- chr: The chromosome the marker is located on
- pos: The physical position of the marker on the chromosome
File: final.all.geno.8.19.25.csv
Description: The full set of SNP calls for all sequenced RILs at all the makers in the final. all.gmap.csv. Is an input file for QTL analysis.
Variables
- Columns correspond to the RIL ID, then every marker, with values of E or S depending on whether the marker corresponds to the reference strain PS312 (E) or not.
Code/software
Code was generated using R v 4.4.1
Packages loaded in scripts include:
geomorph v 4.0.10
Morpho v2 .12
mclust v 6.1.1
factoextra v 1.0.7
ggplot2 v 3.5.1
dplyr v 1.1.4
corrplot v 0.95
qtl2 v 0.36
tidyr v 1.3.1
onemap v 3.0.0
