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:6:91279448:T:C | Chr06:91,279,448 | C | 116 | 2.952 | 4.49e-08 | 0.0498 | TkA06G368510 |
| TK1151:2:115973518:A:T | Chr02:115,973,518 | T | 126 | -4.33 | 3.68e-07 | 0.0865 | TkA02G455250 |
| TK1151:8:70383246:A:C | Chr08:70,383,246 | C | 124 | 3.159 | 7.37e-07 | 0.0865 | TkA08G270160 |
| TK1151:4:4042389:T:A | Chr04:4,042,389 | A | 115 | 1.753 | 8.07e-07 | 0.0865 | TkA04G016310 |
| TK1151:7:74767452:G:A | Chr07:74,767,452 | A | 125 | 7.586 | 1.19e-06 | 0.0865 | TkA07G297620 |
| TK1151:7:74767437:C:T | Chr07:74,767,437 | T | 126 | 7.563 | 1.36e-06 | 0.0865 | TkA07G297620 |
| TK1151:1:112651177:T:C | Chr01:112,651,177 | C | 122 | -10.53 | 1.51e-06 | 0.0865 | TkA01G449580 |
| TK1151:3:147542479:C:T | Chr03:147,542,479 | T | 116 | 2.201 | 1.65e-06 | 0.0865 | TkA03G581390 |
| TK1151:2:41556449:A:G | Chr02:41,556,449 | G | 118 | 2.043 | 1.83e-06 | 0.0865 | TkA02G162070 |
| TK1151:3:75945291:A:G | Chr03:75,945,291 | G | 122 | -10.51 | 1.9e-06 | 0.0865 | TkA03G308830 |
| TK1151:3:75945301:A:T | Chr03:75,945,301 | T | 122 | -10.51 | 1.9e-06 | 0.0865 | TkA03G308830 |
| TK1151:2:130028214:A:C | Chr02:130,028,214 | C | 121 | -10.53 | 1.92e-06 | 0.0865 | TkA02G506880 |
| TK1151:7:85280297:A:T | Chr07:85,280,297 | T | 124 | -10.52 | 2.02e-06 | 0.0865 | TkA07G341770 |
| TK1151:1:20765774:G:A | Chr01:20,765,774 | A | 122 | -3.867 | 2.78e-06 | 0.0865 | TkA01G081350 |
| TK1151:1:99932398:G:A | Chr01:99,932,398 | A | 124 | -10.43 | 2.87e-06 | 0.0865 | TkA01G400290 |
| TK1151:1:30661804:C:T | Chr01:30,661,804 | T | 126 | -10.4 | 2.92e-06 | 0.0865 | TkA01G121250 |
| TK1151:1:64093009:A:C | Chr01:64,093,009 | C | 126 | -10.4 | 2.92e-06 | 0.0865 | TkA01G254100 |
| TK1151:2:19562145:G:C | Chr02:19,562,145 | C | 126 | 10.4 | 2.92e-06 | 0.0865 | TkA02G075560 |
| TK1151:2:30564354:A:C | Chr02:30,564,354 | C | 126 | -10.4 | 2.92e-06 | 0.0865 | TkA02G117290 |
| TK1151:2:53377534:G:A | Chr02:53,377,534 | A | 126 | 10.4 | 2.92e-06 | 0.0865 | TkA02G209120 |
| TK1151:3:110564088:A:G | Chr03:110,564,088 | G | 126 | 10.4 | 2.92e-06 | 0.0865 | TkA03G443530 |
| TK1151:3:115558834:A:T | Chr03:115,558,834 | T | 126 | -10.4 | 2.92e-06 | 0.0865 | TkA03G462970 |
| TK1151:3:119778197:T:A | Chr03:119,778,197 | A | 126 | 10.4 | 2.92e-06 | 0.0865 | TkA03G479000 |
| TK1151:5:78561366:G:A | Chr05:78,561,366 | A | 126 | 10.4 | 2.92e-06 | 0.0865 | TkA05G318300 |
| TK1151:7:74707404:T:G | Chr07:74,707,404 | G | 126 | 10.4 | 2.92e-06 | 0.0865 | TkA07G297570 |
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