Genome writing and Targeted Delivery of the NKX6-3/ANK1 gene cluster and its Type 2 Diabetes GWAS Variants to Human iPSCs

Genome-wide association studies (GWAS) identified over 600 loci containing single-nucleotide polymorphisms (SNPs) associated with type 2 diabetes (T2D), most of which reside in non-coding regions. Among the set of T2D SNPs, linking causal genome variants to disease risk experimentally has remained a challenge; however, advances in synthetic mammalian genome writing techniques now enable the delivery of multiple haplotypes to human induced pluripotent stem cells (hiPSCs) to create a series of isogenic cell lines that can be differentiated and phenotyped in vitro. Here, to begin efforts in dissecting a T2D GWAS locus, we engineered an NKX6-3/ANK1 gene cluster knockout hiPSC line and introduced a landing pad facilitating the delivery of synthetic haplotype payloads. We built four haplotypes, including several that are not observed in nature, containing risk SNPs spanning the NKX6-3/ANK1 gene cluster using a method called "variant Switching Auxotrophic markers for Integration" (vSwAP-In), and integrated them precisely into hiPSCs. NKX6-3/ANK1 deletion blocked pancreatic progenitor and skeletal muscle differentiation, suggesting that NKX6-3 and ANK1 are required for early pancreatic and skeletal muscle development, and perhaps related to the existence of two nonoverlapping sets of SNPs in linkage disequilibrium that associate with the expression of the two adjacent genes. When NKX6-3/ANK1 T2D "Risk" haplotypes were reintroduced, skeletal muscle and pancreatic progenitor differentiation capabilities were restored. ANK1 expression was elevated in the ANK1 Risk and All-Risk haplotypes compared to the NKX6-3 Risk and Non-Risk haplotypes, establishing a functional experimental platform to examine risk SNP clusters in their native contexts. Overall, this work establishes a platform for the dissection of GWAS loci using synthetic haplotype genomics in hiPSCs. Significance StatementGenome-wide association studies have been used to identify disease-associated SNPs; however, most SNPs lie in non-coding regions, making functional experimentation difficult to perform. Using vSwAP-In, a yeast-based DNA variant-building method, and mSwAP-In, a mammalian genome engineering approach, we establish a platform for functional GWAS dissection in hiPSCs. This platform allows us to build DNA harboring virtually any combination of disease-risk SNPs, allowing for functional characterization of SNPs without the limitations of linkage disequilibrium. We demonstrate this approach using a Type 2 diabetes GWAS gene cluster, NKX6-3/ANK1.