Natural rubber content GWAS
GLM and MLM associations, model diagnostics and audited downloads.
GLM: PC1–PC3; MLM: PC1–PC3 and a standardized genomic relationship matrix. Interpretation is limited to this F1 family.
Genome-wide results for the selected model
Red dashed line: fixed exploratory screening at −log10(P)=4 (P≤1×10⁻⁴); gray dotted line: suggestive 1/M. Variants and candidate genes use the fixed cutoff; no empirical permutation calibration was performed.
View P≤1e-4 variants · View original Bonferroni reference plot
QQ plot and model statistics
λ = 0.978；Genome-wide threshold P ≤ 4.5e-08。
Filter association statistics
| Variant | Position (bp) | Effect allele | n | β | P | FDR q | Nearest gene |
|---|---|---|---|---|---|---|---|
| TK1151:8:70383246:A:C | Chr08:70,383,246 | C | 124 | 3.159 | 7.37e-07 | 0.0865 | TkA08G270160 |
Beta is the change in content percentage points per effect allele. Table filters do not alter the genome-wide plot above.
Methods, sources and scope
This is one F1 family. Phenotypes are published mean natural rubber contents; the VCF was generated by this project from CRA009020 reads. Sample correspondence is supported by unique exact read-count matches; consistency with the author’s original bin coding remains unresolved. Results support exploratory within-family association; candidates require independent validation.
Biallelic SNPs on8 chromosomes; nonvarying offspring GT dosage excluded; GT depth ≥5; variant missingness ≤10%; MAF ≥0.05; sample missingness ≤15%; no HWE filtering in the F1 family. GLM uses nonmissing genotypes per variant; MLM handles residual missing calls by its software defaults. Nearest genes indicate proximity; published marker intervals are not full QTL confidence intervals.
Yang et al., 2023Published QTLsMethods & versionsNumerical validationWorkflow & provenanceFurther analysis plan