Deep mutational scanning and massively parallel kinetics of plasminogen activator inhibitor-1 functional stability to probe its latency transition.

Plasminogen activator inhibitor-1 (PAI-1), a member of the serine protease inhibitor superfamily of proteins, is unique among serine protease inhibitors for exhibiting a spontaneous conformational change to a latent or inactive state. The functional half-life for this transition at physiologic temperature and pH is ∼1 to 2 h. To better understand the molecular mechanisms underlying this transition, we now report on the analysis of a comprehensive PAI-1 variant library expressed on filamentous phage and selected for functional stability after 48 h at 37 °C. Of the 7201 possible single amino acid substitutions in PAI-1, we identified 439 that increased the functional stability of PAI-1 beyond that of the WT protein. We also found 1549 single amino acid substitutions that retained inhibitory activity toward the canonical target protease of PAI-1 (urokinase-like plasminogen activator), whereas exhibiting functional stability less than or equal to that of WT PAI-1. Missense mutations that increase PAI-1 functional stability are concentrated in highly flexible regions within the PAI-1 structure. Finally, we developed a method for simultaneously measuring the functional half-lives of hundreds of PAI-1 variants in a multiplexed, massively parallel manner, quantifying the functional half-lives for 697 single missense variants of PAI-1 by this approach. Overall, these findings provide novel insight into the mechanisms underlying the latency transition of PAI-1 and provide a database for interpreting human PAI-1 genetic variants.

Hepatic inactivation of murine Surf4 results in marked reduction in plasma cholesterol.

PCSK9 negatively regulates low-density lipoprotein receptor (LDLR) abundance on the cell surface, leading to decreased hepatic clearance of LDL particles and increased levels of plasma cholesterol. We previously identified SURF4 as a cargo receptor that facilitates PCSK9 secretion in HEK293T cells (Emmer et al., 2018). Here, we generated hepatic SURF4-deficient mice (Surf4fl/fl Alb-Cre+) to investigate the physiologic role of SURF4 in vivo. Surf4fl/fl Alb-Cre+ mice exhibited normal viability, gross development, and fertility. Plasma PCSK9 levels were reduced by ~60% in Surf4fl/fl Alb-Cre+ mice, with a corresponding ~50% increase in steady state LDLR protein abundance in the liver, consistent with SURF4 functioning as a cargo receptor for PCSK9. Surprisingly, these mice exhibited a marked reduction in plasma cholesterol and triglyceride levels out of proportion to the partial increase in hepatic LDLR abundance. Detailed characterization of lipoprotein metabolism in these mice instead revealed a severe defect in hepatic lipoprotein secretion, consistent with prior reports of SURF4 also promoting the secretion of apolipoprotein B (APOB). Despite a small increase in liver mass and lipid content, histologic evaluation revealed no evidence of steatohepatitis or fibrosis in Surf4fl/fl Alb-Cre+ mice. Acute depletion of hepatic SURF4 by CRISPR/Cas9 or liver-targeted siRNA in adult mice confirms these findings. Together, these data support the physiologic significance of SURF4 in the hepatic secretion of PCSK9 and APOB-containing lipoproteins and its potential as a therapeutic target in atherosclerotic cardiovascular diseases.

Identification of novel γ-globin inducers among all potential erythroid druggable targets.

Human γ-globin is predominantly expressed in fetal liver erythroid cells during gestation from 2 nearly identical genes, HBG1 and HBG2, that are both perinatally silenced. Reactivation of these fetal genes in adult red blood cells can ameliorate many symptoms associated with the inherited ß-globinopathies, sickle cell disease, and Cooley anemia. Although promising genetic strategies to reactivate the γ-globin genes to treat these diseases have been explored, there are significant barriers to their effective implementation worldwide; alternatively, pharmacological induction of γ-globin synthesis could readily reach the majority of affected individuals. In this study, we generated a CRISPR knockout library that targeted all erythroid genes for which prospective or actual therapeutic compounds already exist. By probing this library for genes that repress fetal hemoglobin (HbF), we identified several novel, potentially druggable, γ-globin repressors, including VHL and PTEN. We demonstrate that deletion of VHL induces HbF through activation of the HIF1α pathway and that deletion of PTEN induces HbF through AKT pathway stimulation. Finally, we show that small-molecule inhibitors of PTEN and EZH induce HbF in both healthy and ß-thalassemic human primary erythroid cells.

Deep mutational scanning of the plasminogen activator inhibitor-1 functional landscape.

The serine protease inhibitor (SERPIN) plasminogen activator inhibitor-1 (PAI-1) is a key regulator of the fibrinolytic system, inhibiting the serine proteases tissue- and urokinase-type plasminogen activator (tPA and uPA, respectively). Missense variants render PAI-1 non-functional through misfolding, leading to its turnover as a protease substrate, or to a more rapid transition to the latent/inactive state. Deep mutational scanning was performed to evaluate the impact of amino acid sequence variation on PAI-1 inhibition of uPA using an M13 filamentous phage display system. Error prone PCR was used to construct a mutagenized PAI-1 library encompassing ~ 70% of potential single amino acid substitutions. The relative effects of 27% of all possible missense variants on PAI-1 inhibition of uPA were determined using high-throughput DNA sequencing. 826 missense variants demonstrated conserved inhibitory activity while 1137 resulted in loss of PAI-1 inhibitory function. The least evolutionarily conserved regions of PAI-1 were also identified as being the most tolerant of missense mutations. The results of this screen confirm previous low-throughput mutational studies, including those of the reactive center loop. These data provide a powerful resource for explaining structure-function relationships for PAI-1 and for the interpretation of human genomic sequence variants.

Phage display broadly identifies inhibitor-reactive regions in von Willebrand factor.

BACKGROUND: Correction of von Willebrand factor (VWF) deficiency with replacement products containing VWF can lead to the development of anti-VWF alloantibodies (i.e., VWF inhibitors) in patients with severe von Willebrand disease (VWD). OBJECTIVE: Locate inhibitor-reactive regions within VWF using phage display. METHODS: We screened a phage library displaying random, overlapping fragments covering the full-length VWF protein sequence for binding to a commercial anti-VWF antibody or to immunoglobulins from three type 3 VWD patients who developed VWF inhibitors in response to treatment with plasma-derived VWF. Immunoreactive phage clones were identified and quantified by next-generation DNA sequencing (NGS). RESULTS: Next-generation DNA sequencing markedly increased the number of phages analyzed for locating immunoreactive regions within VWF following a single round of selection and identified regions not recognized in previous reports using standard phage display methods. Extending this approach to characterize VWF inhibitors from three type 3 VWD patients (including two siblings homozygous for the same VWF gene deletion) revealed patterns of immunoreactivity distinct from the commercial antibody and between unrelated patients, though with notable areas of overlap. Alloantibody reactivity against the VWF propeptide is consistent with incomplete removal of the propeptide from plasma-derived VWF replacement products. CONCLUSION: These results demonstrate the utility of phage display and NGS to characterize diverse anti-VWF antibody reactivities.

The Endoplasmic Reticulum Cargo Receptor SURF4 Facilitates Efficient Erythropoietin Secretion.

Erythropoietin (EPO) stimulates erythroid differentiation and maturation. Though the transcriptional regulation of EPO has been well studied, the molecular determinants of EPO secretion remain unknown. Here, we generated a HEK293T reporter cell line that provides a quantifiable and selectable readout of intracellular EPO levels and performed a genome-scale CRISPR screen that identified SURF4 as an important mediator of EPO secretion. Targeting SURF4 with multiple independent single guide RNAs (sgRNAs) resulted in intracellular accumulation and extracellular depletion of EPO. Both of these phenotypes were rescued by expression of SURF4 cDNA. Additionally, we found that disruption of SURF4 resulted in accumulation of EPO in the endoplasmic reticulum (ER) compartment and that SURF4 and EPO physically interact. Furthermore, SURF4 disruption in Hep3B cells also caused a defect in the secretion of endogenous EPO under conditions mimicking hypoxia, ruling out an artifact of heterologous overexpression. This work demonstrates that SURF4 functions as an ER cargo receptor that mediates the efficient secretion of EPO. Our findings also suggest that modulating SURF4 may be an effective treatment for disorders of erythropoiesis that are driven by aberrant EPO levels. Finally, we show that SURF4 overexpression results in increased secretion of EPO, suggesting a new strategy for more efficient production of recombinant EPO.

The in vivo endothelial cell translatome is highly heterogeneous across vascular beds.

Endothelial cells (ECs) are highly specialized across vascular beds. However, given their interspersed anatomic distribution, comprehensive characterization of the molecular basis for this heterogeneity in vivo has been limited. By applying endothelial-specific translating ribosome affinity purification (EC-TRAP) combined with high-throughput RNA sequencing analysis, we identified pan EC-enriched genes and tissue-specific EC transcripts, which include both established markers and genes previously unappreciated for their presence in ECs. In addition, EC-TRAP limits changes in gene expression after EC isolation and in vitro expansion, as well as rapid vascular bed-specific shifts in EC gene expression profiles as a result of the enzymatic tissue dissociation required to generate single-cell suspensions for fluorescence-activated cell sorting or single-cell RNA sequencing analysis. Comparison of our EC-TRAP with published single-cell RNA sequencing data further demonstrates considerably greater sensitivity of EC-TRAP for the detection of low abundant transcripts. Application of EC-TRAP to examine the in vivo host response to lipopolysaccharide (LPS) revealed the induction of gene expression programs associated with a native defense response, with marked differences across vascular beds. Furthermore, comparative analysis of whole-tissue and TRAP-selected mRNAs identified LPS-induced differences that would not have been detected by whole-tissue analysis alone. Together, these data provide a resource for the analysis of EC-specific gene expression programs across heterogeneous vascular beds under both physiologic and pathologic conditions.

Genome-wide linkage analysis and whole-exome sequencing identifies an ITGA2B mutation in a family with thrombocytopenia.

Hereditary thrombocytopenias can be subclassified based on mode of inheritance and platelet size. Here we report a family with autosomal dominant (AD) thrombocytopenia with normal platelet size. Linkage analysis and whole exome sequencing identified the R1026W substitution in ITGA2B as the causative defect. The same mutation has been previously reported in 7 Japanese families/patients with AD thrombocytopenia, but all of these patients had macrothrombocytopenia. This is the first report of a family with AD thrombocytopenia with normal platelet size resulting from mutation in ITGA2B. ITGA2B mutations should therefore be included in the differential diagnosis of this latter disorder.

Whole exome sequencing of ENU-induced thrombosis modifier mutations in the mouse.

Although the Factor V Leiden (FVL) gene variant is the most prevalent genetic risk factor for venous thrombosis, only 10% of FVL carriers will experience such an event in their lifetime. To identify potential FVL modifier genes contributing to this incomplete penetrance, we took advantage of a perinatal synthetic lethal thrombosis phenotype in mice homozygous for FVL (F5L/L) and haploinsufficient for tissue factor pathway inhibitor (Tfpi+/-) to perform a sensitized dominant ENU mutagenesis screen. Linkage analysis conducted in the 3 largest pedigrees generated from the surviving F5L/L Tfpi+/- mice ('rescues') using ENU-induced coding variants as genetic markers was unsuccessful in identifying major suppressor loci. Whole exome sequencing was applied to DNA from 107 rescue mice to identify candidate genes enriched for ENU mutations. A total of 3,481 potentially deleterious candidate ENU variants were identified in 2,984 genes. After correcting for gene size and multiple testing, Arl6ip5 was identified as the most enriched gene, though not reaching genome-wide significance. Evaluation of CRISPR/Cas9 induced loss of function in the top 6 genes failed to demonstrate a clear rescue phenotype. However, a maternally inherited (not ENU-induced) de novo mutation (Plcb4R335Q) exhibited significant co-segregation with the rescue phenotype (p = 0.003) in the corresponding pedigree. Thrombosis suppression by heterozygous Plcb4 loss of function was confirmed through analysis of an independent, CRISPR/Cas9-induced Plcb4 mutation (p = 0.01).

Sensitized mutagenesis screen in Factor V Leiden mice identifies thrombosis suppressor loci.

Factor V Leiden (F5L ) is a common genetic risk factor for venous thromboembolism in humans. We conducted a sensitized N-ethyl-N-nitrosourea (ENU) mutagenesis screen for dominant thrombosuppressor genes based on perinatal lethal thrombosis in mice homozygous for F5L (F5L/L ) and haploinsufficient for tissue factor pathway inhibitor (Tfpi+/- ). F8 deficiency enhanced the survival of F5L/LTfpi+/- mice, demonstrating that F5L/LTfpi+/- lethality is genetically suppressible. ENU-mutagenized F5L/L males and F5L/+Tfpi+/- females were crossed to generate 6,729 progeny, with 98 F5L/LTfpi+/- offspring surviving until weaning. Sixteen lines, referred to as "modifier of Factor 5 Leiden (MF5L1-16)," exhibited transmission of a putative thrombosuppressor to subsequent generations. Linkage analysis in MF5L6 identified a chromosome 3 locus containing the tissue factor gene (F3). Although no ENU-induced F3 mutation was identified, haploinsufficiency for F3 (F3+/- ) suppressed F5L/LTfpi+/- lethality. Whole-exome sequencing in MF5L12 identified an Actr2 gene point mutation (p.R258G) as the sole candidate. Inheritance of this variant is associated with suppression of F5L/LTfpi+/- lethality (P = 1.7 × 10-6), suggesting that Actr2p.R258G is thrombosuppressive. CRISPR/Cas9 experiments to generate an independent Actr2 knockin/knockout demonstrated that Actr2 haploinsufficiency is lethal, supporting a hypomorphic or gain-of-function mechanism of action for Actr2p.R258G Our findings identify F8 and the Tfpi/F3 axis as key regulators in determining thrombosis balance in the setting of F5L and also suggest a role for Actr2 in this process.

Massively parallel enzyme kinetics reveals the substrate recognition landscape of the metalloprotease ADAMTS13.

Proteases play important roles in many biologic processes and are key mediators of cancer, inflammation, and thrombosis. However, comprehensive and quantitative techniques to define the substrate specificity profile of proteases are lacking. The metalloprotease ADAMTS13 regulates blood coagulation by cleaving von Willebrand factor (VWF), reducing its procoagulant activity. A mutagenized substrate phage display library based on a 73-amino acid fragment of VWF was constructed, and the ADAMTS13-dependent change in library complexity was evaluated over reaction time points, using high-throughput sequencing. Reaction rate constants (kcat/KM) were calculated for nearly every possible single amino acid substitution within this fragment. This massively parallel enzyme kinetics analysis detailed the specificity of ADAMTS13 and demonstrated the critical importance of the P1-P1' substrate residues while defining exosite binding domains. These data provided empirical evidence for the propensity for epistasis within VWF and showed strong correlation to conservation across orthologs, highlighting evolutionary selective pressures for VWF.

Genetic variants in PLG, LPA, and SIGLEC 14 as well as smoking contribute to plasma plasminogen levels.

Plasminogen is the precursor of the serine protease plasmin, a central enzyme of the fibrinolytic system. Plasma levels of plasminogen vary by almost 2-fold among healthy individuals, yet little is known about its heritability or genetic determinants in the general population. In order to identify genetic factors affecting the natural variation of plasminogen levels, we performed a genome-wide association study and linkage analysis in a sample of 3456 young healthy individuals who participated in the Genes and Blood Clotting Study (GABC) or the Trinity Student Study (TSS). Heritability of plasminogen levels was 48.1% to 60.0%. Tobacco smoking and female sex were associated with higher levels of plasminogen. In the meta-analysis, 11 single-nucleotide polymorphisms (SNPs) in 2 regions reached genome-wide significance (P < 5.0E-8). Of these, 9 SNPs were near the PLG or LPA genes on Chr6q26, whereas 2 were on Chr19q13 and 5' upstream of SIGLEC14. These 11 SNPs represented 4 independent signals and collectively explained 6.8% of plasminogen level variation in the study populations. The strongest association was observed for a nonsynonymous SNP in the PLG gene (R523W). Individuals bearing an additional copy of this allele had an average decrease of 13.4% in plasma plasminogen level.

Absence of a red blood cell phenotype in mice with hematopoietic deficiency of SEC23B.

Congenital dyserythropoietic anemia type II (CDAII) is an autosomal recessive disease of ineffective erythropoiesis characterized by increased bi/multinucleated erythroid precursors in the bone marrow. CDAII results from mutations in SEC23B. The SEC23 protein is a core component of coat protein complex II-coated vesicles, which transport secretory proteins from the endoplasmic reticulum to the Golgi apparatus. Though the genetic defect underlying CDAII has been identified, the pathophysiology of this disease remains unknown. We previously reported that SEC23B-deficient mice die perinatally, exhibiting massive pancreatic degeneration, with this early mortality limiting evaluation of the adult hematopoietic compartment. We now report that mice with SEC23B deficiency restricted to the hematopoietic compartment survive normally and do not exhibit anemia or other CDAII characteristics. We also demonstrate that SEC23B-deficient hematopoietic stem cells (HSC) do not exhibit a disadvantage at reconstituting hematopoiesis when compared directly to wild-type HSC in a competitive repopulation assay. Secondary bone marrow transplants demonstrated continued equivalence of SEC23B-deficient and WT HSC in their hematopoietic reconstitution potential. The surprising discordance in phenotypes between SEC23B-deficient mice and humans may reflect an evolutionary shift in SEC23 paralog function and/or expression, or a change in a specific COPII cargo critical for erythropoiesis.

Linkage analysis identifies a locus for plasma von Willebrand factor undetected by genome-wide association.

The plasma glycoprotein von Willebrand factor (VWF) exhibits fivefold antigen level variation across the normal human population determined by both genetic and environmental factors. Low levels of VWF are associated with bleeding and elevated levels with increased risk for thrombosis, myocardial infarction, and stroke. To identify additional genetic determinants of VWF antigen levels and to minimize the impact of age and illness-related environmental factors, we performed genome-wide association analysis in two young and healthy cohorts (n = 1,152 and n = 2,310) and identified signals at ABO (P < 7.9E-139) and VWF (P < 5.5E-16), consistent with previous reports. Additionally, linkage analysis based on sibling structure within the cohorts, identified significant signals at chromosome 2q12-2p13 (LOD score 5.3) and at the ABO locus on chromosome 9q34 (LOD score 2.9) that explained 19.2% and 24.5% of the variance in VWF levels, respectively. Given its strong effect, the linkage region on chromosome 2 could harbor a potentially important determinant of bleeding and thrombosis risk. The absence of a chromosome 2 association signal in this or previous association studies suggests a causative gene harboring many genetic variants that are individually rare, but in aggregate common. These results raise the possibility that similar loci could explain a significant portion of the "missing heritability" for other complex genetic traits.

Analysis of informed consent document utilization in a minimal-risk genetic study.

BACKGROUND: The signed informed consent document certifies that the process of informed consent has taken place and provides research participants with comprehensive information about their role in the study. Despite efforts to optimize the informed consent document, only limited data are available about the actual use of consent documents by participants in biomedical research. OBJECTIVE: To examine the use of online consent documents in a minimal-risk genetic study. DESIGN: Prospective sibling cohort enrolled as part of a genetic study of hematologic and common human traits. SETTING: University of Michigan Campus, Ann Arbor, Michigan. PARTICIPANTS: Volunteer sample of healthy persons with 1 or more eligible siblings aged 14 to 35 years. Enrollment was through targeted e-mail to student lists. A total of 1209 persons completed the study. MEASUREMENTS: Time taken by participants to review a 2833-word online consent document before indicating consent and identification of a masked hyperlink near the end of the document. RESULTS: The minimum predicted reading time was 566 seconds. The median time to consent was 53 seconds. A total of 23% of participants consented within 10 seconds, and 93% of participants consented in less than the minimum predicted reading time. A total of 2.5% of participants identified the masked hyperlink. LIMITATION: The online consent process was not observed directly by study investigators, and some participants may have viewed the consent document more than once. CONCLUSION: Few research participants thoroughly read the consent document before agreeing to participate in this genetic study. These data suggest that current informed consent documents, particularly for low-risk studies, may no longer serve the intended purpose of protecting human participants, and the role of these documents should be reassessed. PRIMARY FUNDING SOURCE: National Institutes of Health.