Fig. 1
Genome-wide eQTL dynamic changes and characterization. a The proportion of static and dynamic distribution of genome-wide local and distant eQTLs at 20 and 40 DAF. b The number of eQTLs for each gene. The x-axis and y-axis represent the number of eQTLs mapped for each gene and the number of genes in each group, respectively. c Comparison of the explained variance (r2) of local and distant eQTLs for expression at 20 and 40 DAF. d Manhattan plot of BnaA07.WRI1 eGWAS at 20 DAF and 40 DAF. Gene names related to lipid metabolism are labeled on the plot. FUS3, FUSCA3; CD1, CUTIN DEFICIENT 1; SHN3, SHINE3; MYB30, MYB DOMAIN PROTEIN 30; LEC1, LEAFY COTYLEDON 1; PDAT, PHOSPHOLIPID: DIACYLGLYCEROL ACYLTRANSFERASE; ROD1, REDUCED OLEATE DESATURATION 1; TED4, TRACHEARY ELEMENT DIFFERENTIATION-RELATED 4; TGD4, TRIGALACTOSYLDIACYLGLYCEROL 4. e The distribution of overlapping ATAC peaks and eQTLs after randomly permuting the positions of ATAC peaks and eQTLs on the whole genome 1000 times at 20 DAF and 40 DAF. The inverted triangle represents the actual number of eQTLs in OCRs. The dashed lines indicate 95% confidence interval (CI) values