relationships that were reproducible > 80% across runs. Similarly, to aggregate the TF-gene relationships across individuals, we aggregated TF-gene relationships that were detected in 80% of 12 individuals. Across these aggregation steps, we averaged the importance scores for selected TF-gene relationships. We selected likely direct TF-gene relationships by filtering relationships where the gene’s promoter contain the TF’s binding motifs using the cisTarget method implement in pySCENIC package using the pre-built TF-motif promoter binding databases for human accessed from (https://resources.aertslab.org/cistarget/databases/ in January 2023)47. We report the list of discovered TF-gene relationships in the Supplementary Data 1–S17.