Data set of Tanaka et al. "Hidden balance in sexual conflict over inbreeding: antagonistic coevolution between male control and female resistance"

1) Data description

Data description

This dataset contains experimental data collected to investigate sexual conflict over inbreeding in the West Indian sweet potato weevil, Euscepes postfasciatus. The study combined physiological and behavioral experiments to test how inbreeding depression, male accessory gland substances, and kinship influence female mating receptivity and reproductive output.

Specifically, the data include:
(1) measurements of female reproductive output from inbred (full-sib) and outbred (non-sib) matings to quantify inbreeding depression;
(2) results of microinjection experiments testing the effects of concentration and kinship origin of male accessory gland substances on female mating receptivity; and
(3) behavioral observations of female remating under natural copulation conditions following mating with kin or non-kin males.

All experiments were conducted under controlled laboratory conditions, and the data were used to test predictions derived from theoretical models of sexual conflict over inbreeding.

⸻

2) Files and variables

Files and variables

The Dryad repository includes the following data files and analysis scripts:

Data files
•TanakaData_Exp1.csv
Data for Experiment 1 (inbreeding depression). Each row represents a female–male pair.
•TanakaData_Exp2.csv
Data for Experiment 2 (microinjection assay of male accessory gland substances). Each row represents an individual female.
•TanakaData_Exp3.csv
Data for Experiment 3 (natural mating assay). Each row represents an individual female.

Common variables (may vary slightly among files)
•FemaleID: Unique identifier for each female
•Kin: Kinship between mating partners (Kin = full-sib, NonKin = non-sib)
•Fsize: Female body size (left elytron length, mm)
•Fno / OffspringNumber: Number of adult offspring produced
•Inj: Injection treatment (K = kin-derived suspension, cNK = non-kin suspension, NK2 = 2× concentration non-kin suspension)
•Period / Period2: Time (days) until mating or remating
•Fin / Fin3: Event indicator (1 = mating occurred, 0 = censored)

Missing values

Missing data are indicated as NA. Blank cells are not used.

All variable names are consistent with those used in the accompanying R scripts.

⸻

3) Code / software

Code and software

All analyses were conducted using R version 4.5.1 (R Core Team, 2025).
The following R packages were used:
•MASS (for negative binomial and Poisson GLMs)
•survival (for survival analysis and Cox proportional hazards models)
•survminer (for visualization and pairwise survival comparisons)
•ggplot2 (for plotting)

The file TanakaCode.R contains all scripts required to reproduce the statistical analyses and figures reported in the manuscript. The workflow is as follows:
1.Load the CSV data files for each experiment.
2.Perform GLM analyses to estimate inbreeding depression (Experiment 1).
3.Conduct survival analyses and pairwise Peto & Peto tests for mating receptivity (Experiments 2 and 3).
4.Generate figures using ggsurvplot() and related plotting functions.

All software used is free and open-source.

1) 数据集说明

数据集说明

本数据集收录的实验数据，旨在探究西印度群岛甘薯象甲（Euscepes postfasciatus）近交相关的性冲突。研究结合生理学与行为学实验，检验近交衰退、雄性附腺物质以及亲缘关系如何影响雌性交配接受度与繁殖产出。

具体而言，数据集包含以下内容：
（1）近交（全同胞交配）与远交（非同胞交配）组合下的雌性繁殖产出数据，用于量化近交衰退程度；
（2）显微注射实验结果，用于检验雄性附腺物质的浓度与亲缘来源对雌性交配接受度的影响；
（3）雌性与亲缘/非亲缘雄性交配后，在自然交配条件下的再次交配行为观测数据。

所有实验均在可控实验室条件下开展，所得数据用于验证基于近交性冲突理论模型提出的相关预测。

⸻

2) 文件与变量

文件与变量

本Dryad仓储库包含如下数据文件与分析脚本：

数据文件
• TanakaData_Exp1.csv：实验1（近交衰退检测）数据，每一行对应一组雌-雄交配对。
• TanakaData_Exp2.csv：实验2（雄性附腺物质显微注射检测）数据，每一行对应1只雌性个体。
• TanakaData_Exp3.csv：实验3（自然交配检测）数据，每一行对应1只雌性个体。

通用变量（各文件中变量可能略有差异）
• FemaleID：每只雌性个体的唯一标识符
• Kin：交配双方的亲缘关系（Kin=全同胞，NonKin=非同胞）
• Fsize：雌性个体体型（左侧鞘翅长度，单位：mm）
• Fno / OffspringNumber：产出的成虫子代数量
• Inj：注射处理组（K=亲缘来源悬浮液，cNK=非亲缘来源悬浮液，NK2=2倍浓度非亲缘悬浮液）
• Period / Period2：直至首次交配或再次交配的时长（单位：天）
• Fin / Fin3：事件指示变量（1=发生交配，0=删失数据）

缺失值处理
缺失数据以NA标记，空单元格未被使用。

所有变量命名均与配套R脚本中使用的命名保持一致。

⸻

3) 代码与软件

代码与软件

所有分析均基于R版本4.5.1（R核心团队，2025）开展。
所用R包包括：
• MASS：用于负二项分布与泊松广义线性模型（GLMs）分析
• survival：用于生存分析与Cox比例风险模型构建
• survminer：用于生存分析可视化与两两生存比较
• ggplot2：用于绘图可视化

TanakaCode.R文件包含复现论文中统计分析与图表所需的全部脚本，分析流程如下：
1. 加载各实验对应的CSV数据文件
2. 开展广义线性模型分析，估算实验1的近交衰退程度
3. 针对实验2与3的交配接受度数据，开展生存分析及两两Peto & Peto检验
4. 使用ggsurvplot()及相关绘图函数生成图表

所用软件均为免费开源软件。