Background & AimAn unbiased genome-first approach can expand the molecular understanding of specific genes in disease-agnostic biobanks for deeper phenotyping. TM6SF2 represents a good candidate for this approach due to its known association with steatotic liver disease (SLD). MethodsWe screened participants with whole-exome sequences in the Penn Medicine Biobank (N>40,000) and the UK Biobank (UKB, N>200,000) for protein-altering variants in TM6SF2 and evaluated their association with liver phenotypes and clinical outcomes. ResultsThree missense variants in TM6SF2 (E167K, L156P, P216L) were associated with an increased risk of clinically-diagnosed and imaging-proven steatosis, independent of PNPLA3 I48M risk allele and hepatitis B/C (p<0.001). E167K homozygotes had significantly increased risk of SLD (OR:5.38, p<0.001), steatohepatitis (OR:5.76, p<0.05) and hepatocellular carcinoma (OR:11.22, p<0.0001) whilst nominal risk of steatohepatitis were present in L156P (OR:2.3, p<0.001) and P216L (OR:4.5, p<0.05) heterozygotes. In addition, carriers of E167K are at a 3-folded increase of at-risk MASH (OR:2.75, p<0.001). CT-derived liver fat scores were higher in E167K and L156P in an allele-dose manner (p<0.05). This corresponded with the UKB nuclear magnetic resonance-derived lipidomic analyses (N=105,348), revealing all carriers to exhibit lower total cholesterol, triglycerides and total choline. In silico predictions suggested that these missense variants cause structural disruptions in the EXPERA domain, leading to reduced protein function. This hypothesis was strengthened with a gene-burden of rare loss-of-function variants in TM6SF2 demonstrating an increased risk of SLD (OR:4.9, p<0.05), primarily driven by a novel rare stop-gain variant (W35X) in the same directionality. ConclusionThe functional genetic study of protein-altering variants provides insights on the association between loss of TM6SF2 function and steatotic liver disease and poses the bases for future mechanistic studies. Impact and implicationsThe genome-first approach expands insights into genetic risk factors for steatotic liver disease with TM6SF2 being a focal point due to its known association with plasma lipid traits. Our findings validated two missense variants (E167K and L156P) to be associated with increased risk of hepatic steatosis on CT and MRI scans, as well as the risk of clinically diagnosed hepatocellular carcinoma independent of the common PNPLA3 I48M risk variant. Notably, we also identified a predicted deleterious missense variant P216L to drive steatotic risk and demonstrated that an aggregated gene burden of rare variants demonstrating complete loss-of-function to be linked to the risk of hepatic steatosis. Combined, this study sets the stage for future mechanistic investigations into the functional consequences of TM6SF2 variants in metabolic dysfunction-associated steatotic liver disease to determine future therapeutic targets.