There is a high unmet need for developing treatments for nonalcoholic fatty liver disease (NAFLD), for which there are no approved drugs today. Here, we used a human in vitro disease model to understand mechanisms linked to genetic risk variants associated with NAFLD. The model is based on 3D spheroids from primary human hepatocytes from five different donors. Across these donors, we observed highly reproducible differences in the extent of steatosis induction, demonstrating that inter-donor variability is reflected in the in vitro model. Importantly, our data indicates that the genetic variant TM6SF2 E167K, previously associated with increased risk for NAFLD, induces increased hepatocyte fat content by reducing APOB particle secretion. Finally, differences in gene expression pathways involved in cholesterol, fatty acid and glucose metabolism between wild type and TM6SF2 E167K mutation carriers (N = 125) were confirmed in the in vitro model. Our data suggest that the 3D in vitro spheroids can be used to investigate the mechanisms underlying the association of human genetic variants associated with NAFLD. This model may also be suitable to discover new treatments against NAFLD.