Generation of Human Isogenic Induced Pluripotent Stem Cell Lines with CRISPR Prime Editing.

We developed an efficient CRISPR prime editing protocol and generated isogenic-induced pluripotent stem cell (iPSC) lines carrying heterozygous or homozygous alleles for putatively causal single nucleotide variants at six type 2 diabetes loci (ABCC8, MTNR1B, TCF7L2, HNF4A, CAMK1D, and GCK). Our two-step sequence-based approach to first identify transfected cell pools with the highest fraction of edited cells significantly reduced the downstream efforts to isolate single clones of edited cells. We found that prime editing can make targeted genetic changes in iPSC and optimization of system components and guide RNA designs that were critical to achieve acceptable efficiency. Systems utilizing PEmax, epegRNA modifications, and MLH1dn provided significant benefit, producing editing efficiencies of 36-73%. Editing success and pegRNA design optimization required for each variant differed depending on the sequence at the target site. With attention to design, prime editing is a promising approach to generate isogenic iPSC lines, enabling the study of specific genetic changes in a common genetic background.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Subcutaneous adipose tissue splice quantitative trait loci reveal differences in isoform usage associated with cardiometabolic traits.

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r2 > 0.8) located within an intron spanned by the splice junction. We identified 203 sQTLs that colocalized with 141 genome-wide association study (GWAS) signals for cardiometabolic traits, including 25 signals for lipid traits, 24 signals for body mass index (BMI), and 12 signals for waist-hip ratio adjusted for BMI. Among all 141 GWAS signals colocalized with an sQTL, we detected 26 that also colocalized with an exon eQTL for an exon skipped by the sQTL splice junction. At a GWAS signal for high-density lipoprotein cholesterol colocalized with an NR1H3 sQTL splice junction, we show that the alternative splice product encodes an NR1H3 transcription factor that lacks a DNA binding domain and fails to activate transcription. Together, these results detect splicing events and candidate mechanisms that may contribute to gene function at GWAS loci.

Transcriptional and Cellular Diversity of the Human Heart.

BACKGROUND: The human heart requires a complex ensemble of specialized cell types to perform its essential function. A greater knowledge of the intricate cellular milieu of the heart is critical to increase our understanding of cardiac homeostasis and pathology. As recent advances in low-input RNA sequencing have allowed definitions of cellular transcriptomes at single-cell resolution at scale, we have applied these approaches to assess the cellular and transcriptional diversity of the nonfailing human heart. METHODS: Microfluidic encapsulation and barcoding was used to perform single nuclear RNA sequencing with samples from 7 human donors, selected for their absence of overt cardiac disease. Individual nuclear transcriptomes were then clustered based on transcriptional profiles of highly variable genes. These clusters were used as the basis for between-chamber and between-sex differential gene expression analyses and intersection with genetic and pharmacologic data. RESULTS: We sequenced the transcriptomes of 287 269 single cardiac nuclei, revealing 9 major cell types and 20 subclusters of cell types within the human heart. Cellular subclasses include 2 distinct groups of resident macrophages, 4 endothelial subtypes, and 2 fibroblast subsets. Comparisons of cellular transcriptomes by cardiac chamber or sex reveal diversity not only in cardiomyocyte transcriptional programs but also in subtypes involved in extracellular matrix remodeling and vascularization. Using genetic association data, we identified strong enrichment for the role of cell subtypes in cardiac traits and diseases. Intersection of our data set with genes on cardiac clinical testing panels and the druggable genome reveals striking patterns of cellular specificity. CONCLUSIONS: Using large-scale single nuclei RNA sequencing, we defined the transcriptional and cellular diversity in the normal human heart. Our identification of discrete cell subtypes and differentially expressed genes within the heart will ultimately facilitate the development of new therapeutics for cardiovascular diseases.

Genetic Reduction in Left Ventricular Protein Kinase C-α and Adverse Ventricular Remodeling in Human Subjects.

BACKGROUND: Inhibition of PKC-α (protein kinase C-α) enhances contractility and cardioprotection in animal models, but effects in humans are unknown. Genotypes at rs9912468 strongly associate with PRKCA expression in the left ventricle, enabling genetic approaches to measure effects of reduced PKC-α in human populations. METHODS AND RESULTS: We analyzed the cis expression quantitative trait locus for PRKCA marked by rs9912468 using 313 left ventricular specimens from European Ancestry patients. The forward strand minor allele (G) at rs9912468 is associated with reduced PKC-α transcript abundance (1.7-fold reduction in minor allele homozygotes, P=1×10-41). This association was cardiac specific in expression quantitative trait locus data sets that span 16 human tissues. Cardiac epigenomic data revealed a predicted enhancer in complete (R2=1.0) linkage disequilibrium with rs9912468 within intron 2 of PRKCA. We cloned this region and used reporter constructs to verify cardiac-specific enhancer activity in vitro in cardiac and noncardiac cells and in vivo in zebrafish. The PRKCA enhancer contains 2 common genetic variants and 4 haplotypes; the haplotype correlated with the rs9912468 PKC-α-lowering allele (G) showed lowest activity. In contrast to previous reports in animal models, the PKC-α-lowering allele is associated with adverse left ventricular remodeling (higher mass, larger diastolic dimension), reduced fractional shortening, and higher risk of dilated cardiomyopathy in human populations. CONCLUSIONS: These findings support PKC-α as a regulator of the human heart but suggest that PKC-α inhibition may adversely affect the left ventricle depending on timing and duration. Pharmacological studies in human subjects are required to discern potential benefits and harms of PKC-α inhibitors as an approach to treat heart disease.

Novel Mutation in FLNC (Filamin C) Causes Familial Restrictive Cardiomyopathy.

BACKGROUND: Restrictive cardiomyopathy (RCM) is a rare cardiomyopathy characterized by impaired diastolic ventricular function resulting in a poor clinical prognosis. Rarely, heritable forms of RCM have been reported, and mutations underlying RCM have been identified in genes that govern the contractile function of the cardiomyocytes. METHODS AND RESULTS: We evaluated 8 family members across 4 generations by history, physical examination, electrocardiography, and echocardiography. Affected individuals presented with a pleitropic syndrome of progressive RCM, atrioventricular septal defects, and a high prevalence of atrial fibrillation. Exome sequencing of 5 affected members identified a single novel missense variant in a highly conserved residue of FLNC (filamin C; p.V2297M). FLNC encodes filamin C-a protein that acts as both a scaffold for the assembly and organization of the central contractile unit of striated muscle and also as a mechanosensitive signaling molecule during cell migration and shear stress. Immunohistochemical analysis of FLNC localization in cardiac tissue from an affected family member revealed a diminished localization at the z disk, whereas traditional localization at the intercalated disk was preserved. Stem cell-derived cardiomyocytes mutated to carry the effect allele had diminished contractile activity when compared with controls. CONCLUSION: We have identified a novel variant in FLNC as pathogenic variant for familial RCM-a finding that further expands on the genetic basis of this rare and morbid cardiomyopathy.

Diminished PRRX1 Expression Is Associated With Increased Risk of Atrial Fibrillation and Shortening of the Cardiac Action Potential.

BACKGROUND: Atrial fibrillation (AF) affects over 33 million individuals worldwide. Genome-wide association studies have identified at least 30 AF loci, but the mechanisms through which individual variants lead to altered disease risk have remained unclear for the majority of these loci. At the 1q24 locus, we hypothesized that the transcription factor PRRX1 could be a strong candidate gene as it is expressed in the pulmonary veins, a source of AF in many individuals. We sought to identify the molecular mechanism, whereby variation at 1q24 may lead to AF susceptibility. METHODS AND RESULTS: We sequenced a ≈158 kb region encompassing PRRX1 in 962 individuals with and without AF. We identified a broad region of association with AF at the 1q24 locus. Using in silico prediction and functional validation, we identified an enhancer that interacts with the promoter of PRRX1 in cells of cardiac lineage. Within this enhancer, we identified a single-nucleotide polymorphism, rs577676, which alters enhancer activity in a mouse atrial cell line and in embryonic zebrafish and differentially regulates PRRX1 expression in human left atria. We found that suppression of PRRX1 in human embryonic stem cell-derived cardiomyocytes and embryonic zebrafish resulted in shortening of the atrial action potential duration, a hallmark of AF. CONCLUSIONS: We have identified a functional genetic variant that alters PRRX1 expression, ultimately resulting in electrophysiological alterations in atrial myocytes that may promote AF.

Gain-of-function mutations in GATA6 lead to atrial fibrillation.

BACKGROUND: The genetic basis of atrial fibrillation (AF) and congenital heart disease remains incompletely understood. OBJECTIVE: We sought to determine the causative mutation in a family with AF, atrial septal defects, and ventricular septal defects. METHODS: We evaluated a pedigree with 16 family members, 1 with an atrial septal defect, 1 with a ventricular septal defect, and 3 with AF; we performed whole exome sequencing in 3 affected family members. Given that early-onset AF was prominent in the family, we then screened individuals with early-onset AF, defined as an age of onset <66 years, for mutations in GATA6. Variants were functionally characterized using reporter assays in a mammalian cell line. RESULTS: Exome sequencing in 3 affected individuals identified a conserved mutation, R585L, in the transcription factor gene GATA6. In the Massachusetts General Hospital Atrial Fibrillation (MGH AF) Study, the mean age of AF onset was 47.1 ± 10.9 years; 79% of the participants were men; and there was no evidence of structural heart disease. We identified 3 GATA6 variants (P91S, A177T, and A543G). Using wild-type and mutant GATA6 constructs driving atrial natriuretic peptide promoter reporter, we found that 3 of the 4 variants had a marked upregulation of luciferase activity (R585L: 4.1-fold, P < .0001; P91S: 2.5-fold, P = .0002; A177T; 1.7-fold, P = .03). In addition, when co-overexpressed with GATA4 and MEF2C, GATA6 variants exhibited upregulation of the alpha myosin heavy chain and atrial natriuretic peptide reporter activity. CONCLUSION: Overall, we found gain-of-function mutations in GATA6 in both a family with early-onset AF and atrioventricular septal defects as well as in a family with sporadic, early-onset AF.

A Multifaceted Role for Myd88-Dependent Signaling in Progression of Murine Mammary Carcinoma.

Previous data obtained in our laboratory suggested that there may be constitutive signaling through the myeloid differentiation primary response gene 88 (Myd88)-dependent signaling cascade in murine mammary carcinoma. Here, we extended these findings by showing that, in the absence of an added Toll-like receptor (TLR) agonist, the myddosome complex was preformed in 4T1 tumor cells, and that Myd88 influenced cytoplasmic extracellular signal-regulated kinase (Erk)1/Erk2 levels, nuclear levels of nuclear factor-kappaB (NFκB) and signal transducer and activator of transcription 5 (STAT5), tumor-derived chemokine (C-C motif) ligand 2 (CCL2) expression, and in vitro and in vivo tumor growth. In addition, RNA-sequencing revealed that Myd88-dependent signaling enhanced the expression of genes that could contribute to breast cancer progression and genes previously associated with poor outcome for patients with breast cancer, in addition to suppressing the expression of genes capable of inhibiting breast cancer progression. Yet, Myd88-dependent signaling in tumor cells also suppressed expression of genes that could contribute to tumor progression. Collectively, these data revealed a multifaceted role for Myd88-dependent signaling in murine mammary carcinoma.