Seasonal human coronavirus antibodies are boosted upon SARS-CoV-2 infection but not associated with protection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly spread within the human population. Although SARS-CoV-2 is a novel coronavirus, most humans had been previously exposed to other antigenically distinct common seasonal human coronaviruses (hCoVs) before the coronavirus disease 2019 (COVID-19) pandemic. Here, we quantified levels of SARS-CoV-2-reactive antibodies and hCoV-reactive antibodies in serum samples collected from 431 humans before the COVID-19 pandemic. We then quantified pre-pandemic antibody levels in serum from a separate cohort of 251 individuals who became PCR-confirmed infected with SARS-CoV-2. Finally, we longitudinally measured hCoV and SARS-CoV-2 antibodies in the serum of hospitalized COVID-19 patients. Our studies indicate that most individuals possessed hCoV-reactive antibodies before the COVID-19 pandemic. We determined that ∼20% of these individuals possessed non-neutralizing antibodies that cross-reacted with SARS-CoV-2 spike and nucleocapsid proteins. These antibodies were not associated with protection against SARS-CoV-2 infections or hospitalizations, but they were boosted upon SARS-CoV-2 infection.

Prior vaccination promotes early activation of memory T cells and enhances immune responses during SARS-CoV-2 breakthrough infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection of vaccinated individuals is increasingly common but rarely results in severe disease, likely due to the enhanced potency and accelerated kinetics of memory immune responses. However, there have been few opportunities to rigorously study early recall responses during human viral infection. To better understand human immune memory and identify potential mediators of lasting vaccine efficacy, we used high-dimensional flow cytometry and SARS-CoV-2 antigen probes to examine immune responses in longitudinal samples from vaccinated individuals infected during the Omicron wave. These studies revealed heightened spike-specific responses during infection of vaccinated compared to unvaccinated individuals. Spike-specific cluster of differentiation (CD)4 T cells and plasmablasts expanded and CD8 T cells were robustly activated during the first week. In contrast, memory B cell activation, neutralizing antibody production and primary responses to nonspike antigens occurred during the second week. Collectively, these data demonstrate the functionality of vaccine-primed immune memory and highlight memory T cells as rapid responders during SARS-CoV-2 infection.

Whole-genome sequencing reveals host factors underlying critical COVID-19.

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2-4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes-including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)-in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease.

Gene burden analysis identifies genes associated with increased risk and severity of adult-onset hearing loss in a diverse hospital-based cohort.

Loss or absence of hearing is common at both extremes of human lifespan, in the forms of congenital deafness and age-related hearing loss. While these are often studied separately, there is increasing evidence that their genetic basis is at least partially overlapping. In particular, both common and rare variants in genes associated with monogenic forms of hearing loss also contribute to the more polygenic basis of age-related hearing loss. Here, we directly test this model in the Penn Medicine BioBank-a healthcare system cohort of around 40,000 individuals with linked genetic and electronic health record data. We show that increased burden of predicted deleterious variants in Mendelian hearing loss genes is associated with increased risk and severity of adult-onset hearing loss. As a specific example, we identify one gene-TCOF1, responsible for a syndromic form of congenital hearing loss-in which deleterious variants are also associated with adult-onset hearing loss. We also identify four additional novel candidate genes (COL5A1, HMMR, RAPGEF3, and NNT) in which rare variant burden may be associated with hearing loss. Our results confirm that rare variants in Mendelian hearing loss genes contribute to polygenic risk of hearing loss, and emphasize the utility of healthcare system cohorts to study common complex traits and diseases.

Germline Mutations in CIDEB and Protection against Liver Disease.

BACKGROUND: Exome sequencing in hundreds of thousands of persons may enable the identification of rare protein-coding genetic variants associated with protection from human diseases like liver cirrhosis, providing a strategy for the discovery of new therapeutic targets. METHODS: We performed a multistage exome sequencing and genetic association analysis to identify genes in which rare protein-coding variants were associated with liver phenotypes. We conducted in vitro experiments to further characterize associations. RESULTS: The multistage analysis involved 542,904 persons with available data on liver aminotransferase levels, 24,944 patients with various types of liver disease, and 490,636 controls without liver disease. We found that rare coding variants in APOB, ABCB4, SLC30A10, and TM6SF2 were associated with increased aminotransferase levels and an increased risk of liver disease. We also found that variants in CIDEB, which encodes a structural protein found in hepatic lipid droplets, had a protective effect. The burden of rare predicted loss-of-function variants plus missense variants in CIDEB (combined carrier frequency, 0.7%) was associated with decreased alanine aminotransferase levels (beta per allele, -1.24 U per liter; 95% confidence interval [CI], -1.66 to -0.83; P = 4.8×10-9) and with 33% lower odds of liver disease of any cause (odds ratio per allele, 0.67; 95% CI, 0.57 to 0.79; P = 9.9×10-7). Rare coding variants in CIDEB were associated with a decreased risk of liver disease across different underlying causes and different degrees of severity, including cirrhosis of any cause (odds ratio per allele, 0.50; 95% CI, 0.36 to 0.70). Among 3599 patients who had undergone bariatric surgery, rare coding variants in CIDEB were associated with a decreased nonalcoholic fatty liver disease activity score (beta per allele in score units, -0.98; 95% CI, -1.54 to -0.41 [scores range from 0 to 8, with higher scores indicating more severe disease]). In human hepatoma cell lines challenged with oleate, CIDEB small interfering RNA knockdown prevented the buildup of large lipid droplets. CONCLUSIONS: Rare germline mutations in CIDEB conferred substantial protection from liver disease. (Funded by Regeneron Pharmaceuticals.).

Impact of natural selection on global patterns of genetic variation and association with clinical phenotypes at genes involved in SARS-CoV-2 infection.

Human genomic diversity has been shaped by both ancient and ongoing challenges from viruses. The current coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a devastating impact on population health. However, genetic diversity and evolutionary forces impacting host genes related to SARS-CoV-2 infection are not well understood. We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (angiotensin converting enzyme 2 [ACE2], transmembrane protease serine 2 [TMPRSS2], dipeptidyl peptidase 4 [DPP4], and lymphocyte antigen 6 complex locus E [LY6E]). We analyzed data from 2,012 ethnically diverse Africans and 15,977 individuals of European and African ancestry with electronic health records and integrated with global data from the 1000 Genomes Project. At ACE2, we identified 41 nonsynonymous variants that were rare in most populations, several of which impact protein function. However, three nonsynonymous variants (rs138390800, rs147311723, and rs145437639) were common among central African hunter-gatherers from Cameroon (minor allele frequency 0.083 to 0.164) and are on haplotypes that exhibit signatures of positive selection. We identify signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage compared with the chimpanzee genome. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19. Our study provides insights into global variation at host genes related to SARS-CoV-2 infection, which have been shaped by natural selection in some populations, possibly due to prior viral infections.

Multiancestry Genome-Wide Association Study of Aortic Stenosis Identifies Multiple Novel Loci in the Million Veteran Program.

BACKGROUND: Calcific aortic stenosis (CAS) is the most common valvular heart disease in older adults and has no effective preventive therapies. Genome-wide association studies (GWAS) can identify genes influencing disease and may help prioritize therapeutic targets for CAS. METHODS: We performed a GWAS and gene association study of 14 451 patients with CAS and 398 544 controls in the Million Veteran Program. Replication was performed in the Million Veteran Program, Penn Medicine Biobank, Mass General Brigham Biobank, BioVU, and BioMe, totaling 12 889 cases and 348 094 controls. Causal genes were prioritized from genome-wide significant variants using polygenic priority score gene localization, expression quantitative trait locus colocalization, and nearest gene methods. CAS genetic architecture was compared with that of atherosclerotic cardiovascular disease. Causal inference for cardiometabolic biomarkers in CAS was performed using Mendelian randomization and genome-wide significant loci were characterized further through phenome-wide association study. RESULTS: We identified 23 genome-wide significant lead variants in our GWAS representing 17 unique genomic regions. Of the 23 lead variants, 14 were significant in replication, representing 11 unique genomic regions. Five replicated genomic regions were previously known risk loci for CAS (PALMD, TEX41, IL6, LPA, FADS) and 6 were novel (CEP85L, FTO, SLMAP, CELSR2, MECOM, CDAN1). Two novel lead variants were associated in non-White individuals (P<0.05): rs12740374 (CELSR2) in Black and Hispanic individuals and rs1522387 (SLMAP) in Black individuals. Of the 14 replicated lead variants, only 2 (rs10455872 [LPA], rs12740374 [CELSR2]) were also significant in atherosclerotic cardiovascular disease GWAS. In Mendelian randomization, lipoprotein(a) and low-density lipoprotein cholesterol were both associated with CAS, but the association between low-density lipoprotein cholesterol and CAS was attenuated when adjusting for lipoprotein(a). Phenome-wide association study highlighted varying degrees of pleiotropy, including between CAS and obesity at the FTO locus. However, the FTO locus remained associated with CAS after adjusting for body mass index and maintained a significant independent effect on CAS in mediation analysis. CONCLUSIONS: We performed a multiancestry GWAS in CAS and identified 6 novel genomic regions in the disease. Secondary analyses highlighted the roles of lipid metabolism, inflammation, cellular senescence, and adiposity in the pathobiology of CAS and clarified the shared and differential genetic architectures of CAS with atherosclerotic cardiovascular diseases.

A genome-first approach to rare variants in hypertrophic cardiomyopathy genes MYBPC3 and MYH7 in a medical biobank.

'Genome-first' approaches to analyzing rare variants can reveal new insights into human biology and disease. Because pathogenic variants are often rare, new discovery requires aggregating rare coding variants into 'gene burdens' for sufficient power. However, a major challenge is deciding which variants to include in gene burden tests. Pathogenic variants in MYBPC3 and MYH7 are well-known causes of hypertrophic cardiomyopathy (HCM), and focusing on these 'positive control' genes in a genome-first approach could help inform variant selection methods and gene burdening strategies for other genes and diseases. Integrating exome sequences with electronic health records among 41 759 participants in the Penn Medicine BioBank, we evaluated the performance of aggregating predicted loss-of-function (pLOF) and/or predicted deleterious missense (pDM) variants in MYBPC3 and MYH7 for gene burden phenome-wide association studies (PheWAS). The approach to grouping rare variants for these two genes produced very different results: pLOFs but not pDM variants in MYBPC3 were strongly associated with HCM, whereas the opposite was true for MYH7. Detailed review of clinical charts revealed that only 38.5% of patients with HCM diagnoses carrying an HCM-associated variant in MYBPC3 or MYH7 had a clinical genetic test result. Additionally, 26.7% of MYBPC3 pLOF carriers without HCM diagnoses had clear evidence of left atrial enlargement and/or septal/LV hypertrophy on echocardiography. Our study shows the importance of evaluating both pLOF and pDM variants for gene burden testing in future studies to uncover novel gene-disease relationships and identify new pathogenic loss-of-function variants across the human genome through genome-first analyses of healthcare-based populations.

Bi-allelic variants in the ER quality-control mannosidase gene EDEM3 cause a congenital disorder of glycosylation.

EDEM3 encodes a protein that converts Man8GlcNAc2 isomer B to Man7-5GlcNAc2. It is involved in the endoplasmic reticulum-associated degradation pathway, responsible for the recognition of misfolded proteins that will be targeted and translocated to the cytosol and degraded by the proteasome. In this study, through a combination of exome sequencing and gene matching, we have identified seven independent families with 11 individuals with bi-allelic protein-truncating variants and one individual with a compound heterozygous missense variant in EDEM3. The affected individuals present with an inherited congenital disorder of glycosylation (CDG) consisting of neurodevelopmental delay and variable facial dysmorphisms. Experiments in human fibroblast cell lines, human plasma, and mouse plasma and brain tissue demonstrated decreased trimming of Man8GlcNAc2 isomer B to Man7GlcNAc2, consistent with loss of EDEM3 enzymatic activity. In human cells, Man5GlcNAc2 to Man4GlcNAc2 conversion is also diminished with an increase of Glc1Man5GlcNAc2. Furthermore, analysis of the unfolded protein response showed a reduced increase in EIF2AK3 (PERK) expression upon stimulation with tunicamycin as compared to controls, suggesting an impaired unfolded protein response. The aberrant plasma N-glycan profile provides a quick, clinically available test for validating variants of uncertain significance that may be identified by molecular genetic testing. We propose to call this deficiency EDEM3-CDG.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Evinacumab Reduces Triglyceride-Rich Lipoproteins in Patients with Hyperlipidemia: A Post-Hoc Analysis of Three Randomized Clinical Trials.

PURPOSE: Natural selection (Mendelian randomization) studies support a causal relationship between elevated triglyceride-rich lipoproteins (TRLs) and atherosclerotic cardiovascular disease (ASCVD). This post-hoc analysis assessed the efficacy of evinacumab in reducing TRLs in patient cohorts from three separate clinical trials with evinacumab. METHODS: Patients with homozygous familial hypercholesterolemia (HoFH) and low-density lipoprotein cholesterol (LDL-C) ≥ 70 mg/dL were enrolled in a phase III trial (R1500-CL-1629; NCT03399786). Patients diagnosed with refractory hypercholesterolemia, with LDL-C ≥ 70 mg/dL or ≥ 100 mg/dL for those with or without ASCVD, respectively, were enrolled in a phase II trial (R1500-CL-1643; NCT03175367). Patients with severe hypertriglyceridemia (fasting TGs ≥ 500 mg/dL) were enrolled in a phase II trial (R1500-HTG-1522; NCT03452228). Patients received evinacumab intravenously (5 or 15 mg/kg) every 4 weeks, or subcutaneously (300 or 450 mg) every week or every 2 weeks. Efficacy outcomes included change in TRLs (calculated as total cholesterol minus high-density lipoprotein cholesterol minus LDL-C) and other lipid parameters from baseline to 12, 16, or 24 weeks for trial 1522, 1643, and 1629, respectively. RESULTS: At baseline, TRL levels were higher for patients with severe hypertriglyceridemia entering the 1522 trial vs. other cohorts. Reductions in TRLs were observed across all studies with evinacumab, with > 50% reduction from baseline observed at the highest doses evaluated in patients with HoFH or refractory hypercholesterolemia. Within all three trials, evinacumab was generally well tolerated. CONCLUSIONS: Despite limitations in direct comparisons between study groups, these data indicate that TRL levels could be a future target for lipid-lowering therapies.

Endothelial lipase mediates efficient lipolysis of triglyceride-rich lipoproteins.

Triglyceride-rich lipoproteins (TRLs) are circulating reservoirs of fatty acids used as vital energy sources for peripheral tissues. Lipoprotein lipase (LPL) is a predominant enzyme mediating triglyceride (TG) lipolysis and TRL clearance to provide fatty acids to tissues in animals. Physiological and human genetic evidence support a primary role for LPL in hydrolyzing TRL TGs. We hypothesized that endothelial lipase (EL), another extracellular lipase that primarily hydrolyzes lipoprotein phospholipids may also contribute to TRL metabolism. To explore this, we studied the impact of genetic EL loss-of-function on TRL metabolism in humans and mice. Humans carrying a loss-of-function missense variant in LIPG, p.Asn396Ser (rs77960347), demonstrated elevated plasma TGs and elevated phospholipids in TRLs, among other lipoprotein classes. Mice with germline EL deficiency challenged with excess dietary TG through refeeding or a high-fat diet exhibited elevated TGs, delayed dietary TRL clearance, and impaired TRL TG lipolysis in vivo that was rescued by EL reconstitution in the liver. Lipidomic analyses of postprandial plasma from high-fat fed Lipg-/- mice demonstrated accumulation of phospholipids and TGs harboring long-chain polyunsaturated fatty acids (PUFAs), known substrates for EL lipolysis. In vitro and in vivo, EL and LPL together promoted greater TG lipolysis than either extracellular lipase alone. Our data positions EL as a key collaborator of LPL to mediate efficient lipolysis of TRLs in humans and mice.

Targeting the coronavirus nucleocapsid protein through GSK-3 inhibition.

The coronaviruses responsible for severe acute respiratory syndrome (SARS-CoV), COVID-19 (SARS-CoV-2), Middle East respiratory syndrome-CoV, and other coronavirus infections express a nucleocapsid protein (N) that is essential for viral replication, transcription, and virion assembly. Phosphorylation of N from SARS-CoV by glycogen synthase kinase 3 (GSK-3) is required for its function and inhibition of GSK-3 with lithium impairs N phosphorylation, viral transcription, and replication. Here we report that the SARS-CoV-2 N protein contains GSK-3 consensus sequences and that this motif is conserved in diverse coronaviruses, raising the possibility that SARS-CoV-2 may be sensitive to GSK-3 inhibitors, including lithium. We conducted a retrospective analysis of lithium use in patients from three major health systems who were PCR-tested for SARS-CoV-2. We found that patients taking lithium have a significantly reduced risk of COVID-19 (odds ratio = 0.51 [0.35-0.74], P = 0.005). We also show that the SARS-CoV-2 N protein is phosphorylated by GSK-3. Knockout of GSK3A and GSK3B demonstrates that GSK-3 is essential for N phosphorylation. Alternative GSK-3 inhibitors block N phosphorylation and impair replication in SARS-CoV-2 infected lung epithelial cells in a cell-type-dependent manner. Targeting GSK-3 may therefore provide an approach to treat COVID-19 and future coronavirus outbreaks.

Novel therapeutic opportunities for familial lecithin:cholesterol acyltransferase deficiency: promises and challenges.

PURPOSE OF REVIEW: Genetic lecithin:cholesterol acyltransferase (LCAT) deficiency is a rare, inherited, recessive disease, which manifests as two different syndromes: Familial LCAT deficiency (FLD) and Fish-eye disease (FED), characterized by low HDL-C and corneal opacity. FLD patients also develop anaemia and renal disease. There is currently no therapy for FLD, but novel therapeutics are at different stages of development. Here, we summarize the most recent advances and the opportunities for and barriers to the further development of such therapies. RECENT FINDINGS: Recent publications highlight the heterogeneous phenotype of FLD and the uncertainty over the natural history of disease and the factors contributing to disease progression. Therapies that restore LCAT function (protein and gene replacement therapies and LCAT activators) showed promising effects on markers of LCAT activity. Although they do not restore LCAT function, HDL mimetics may slow renal disease progression. SUMMARY: The further development of novel therapeutics requires the identification of efficacy endpoints, which include quantitative biomarkers of disease progression. Because of the heterogeneity of renal disease progression among FLD individuals, future treatments for FLD will have to be tailored based on the specific clinical characteristics of the patient. Extensive studies of the natural history and biomarkers of the disease will be required to achieve this goal.

Web of Science's Citation Median Metrics Overcome the Major Constraints of the Journal Impact Factor.

There are many metrics to evaluate the performance and status of journals. Among these, the journal impact factor (JIF) has become the dominant metric. The influence of JIF is illustrated by its widespread use to evaluate academic status, compensation, and funding decisions. However, as noted by Clarivate Analytics, the parent company of the Web of Science (WoS), the JIF should not be used without careful attention to the many phenomena that influence citation rates. To facilitate transparency, Clarivate Analytics provides all data used to determine the JIF. In addition, WoS provides other metrics for journal evaluation, including the article citation median and the review citation median. These metrics are represented as medians to minimize the confounding influence of a small number of highly cited articles that may occur when data are represented as means. Another feature of these WoS metrics is that data are separated according to different publication types of article (original research and review). To systematically compare these selected metrics, we used the data provided on the WoS web site to analyze 25 top ranked cardiovascular journals in the same mode as represented in the WoS citation distribution window. The results indicate that the article citation median and review citation median overcome several concerns that have been raised about the JIF and seem to provide enhanced objectivity as an indicator of journal impact in publishing original research and reviews. Therefore, we advocate that these additional WoS metrics might be preferentially considered as indicators of journal performance.

ILRUN Promotes Atherosclerosis Through Lipid-Dependent and Lipid-Independent Factors.

BACKGROUND: Common genetic variation in close proximity to the ILRUN gene are significantly associated with coronary artery disease as well as with plasma lipid traits. We recently demonstrated that hepatic inflammation and lipid regulator with ubiquitin-associated domain-like and NBR1-like domains (ILRUN) regulates lipoprotein metabolism in vivo in mice. However, whether ILRUN, which is expressed in vascular cells, directly impacts atherogenesis remains unclear. We sought to determine the role of ILRUN in atherosclerosis development in mice. METHODS: For our study, we generated global Ilrun-deficient (IlrunKO) male and female mice on 2 hyperlipidemic backgrounds: low density lipoprotein receptor knockout (LdlrKO) and apolipoprotein E knockout (ApoeKO; double knockout [DKO]). RESULTS: Compared with littermate control mice (single LdlrKO or ApoeKO), deletion of Ilrun in DKO mice resulted in significantly attenuated both early and advanced atherosclerotic lesion development, as well as reduced necrotic area. DKO mice also had significantly decreased plasma cholesterol levels, primarily attributable to non-HDL (high-density lipoprotein) cholesterol. Hepatic-specific reconstitution of ILRUN in DKO mice on the ApoeKO background normalized plasma lipids, but atherosclerotic lesion area and necrotic area remained reduced in DKO mice. Further analysis showed that loss of Ilrun increased efferocytosis receptor MerTK expression in macrophages, enhanced in vitro efferocytosis, and significantly improved in situ efferocytosis in advanced lesions. CONCLUSIONS: Our results support ILRUN as an important novel regulator of atherogenesis that promotes lesion progression and necrosis. It influences atherosclerosis through both plasma lipid-dependent and lipid-independent mechanisms. These findings support ILRUN as the likely causal gene responsible for genetic association of variants with coronary artery disease at this locus and suggest that suppression of ILRUN activity might be expected to reduce atherosclerosis.

Human knockouts and phenotypic analysis in a cohort with a high rate of consanguinity.

A major goal of biomedicine is to understand the function of every gene in the human genome. Loss-of-function mutations can disrupt both copies of a given gene in humans and phenotypic analysis of such 'human knockouts' can provide insight into gene function. Consanguineous unions are more likely to result in offspring carrying homozygous loss-of-function mutations. In Pakistan, consanguinity rates are notably high. Here we sequence the protein-coding regions of 10,503 adult participants in the Pakistan Risk of Myocardial Infarction Study (PROMIS), designed to understand the determinants of cardiometabolic diseases in individuals from South Asia. We identified individuals carrying homozygous predicted loss-of-function (pLoF) mutations, and performed phenotypic analysis involving more than 200 biochemical and disease traits. We enumerated 49,138 rare (<1% minor allele frequency) pLoF mutations. These pLoF mutations are estimated to knock out 1,317 genes, each in at least one participant. Homozygosity for pLoF mutations at PLA2G7 was associated with absent enzymatic activity of soluble lipoprotein-associated phospholipase A2; at CYP2F1, with higher plasma interleukin-8 concentrations; at TREH, with lower concentrations of apoB-containing lipoprotein subfractions; at either A3GALT2 or NRG4, with markedly reduced plasma insulin C-peptide concentrations; and at SLC9A3R1, with mediators of calcium and phosphate signalling. Heterozygous deficiency of APOC3 has been shown to protect against coronary heart disease; we identified APOC3 homozygous pLoF carriers in our cohort. We recruited these human knockouts and challenged them with an oral fat load. Compared with family members lacking the mutation, individuals with APOC3 knocked out displayed marked blunting of the usual post-prandial rise in plasma triglycerides. Overall, these observations provide a roadmap for a 'human knockout project', a systematic effort to understand the phenotypic consequences of complete disruption of genes in humans.

Genomic profiling of human vascular cells identifies TWIST1 as a causal gene for common vascular diseases.

Genome-wide association studies have identified multiple novel genomic loci associated with vascular diseases. Many of these loci are common non-coding variants that affect the expression of disease-relevant genes within coronary vascular cells. To identify such genes on a genome-wide level, we performed deep transcriptomic analysis of genotyped primary human coronary artery smooth muscle cells (HCASMCs) and coronary endothelial cells (HCAECs) from the same subjects, including splicing Quantitative Trait Loci (sQTL), allele-specific expression (ASE), and colocalization analyses. We identified sQTLs for TARS2, YAP1, CFDP1, and STAT6 in HCASMCs and HCAECs, and 233 ASE genes, a subset of which are also GTEx eGenes in arterial tissues. Colocalization of GWAS association signals for coronary artery disease (CAD), migraine, stroke and abdominal aortic aneurysm with GTEx eGenes in aorta, coronary artery and tibial artery discovered novel candidate risk genes for these diseases. At the CAD and stroke locus tagged by rs2107595 we demonstrate colocalization with expression of the proximal gene TWIST1. We show that disrupting the rs2107595 locus alters TWIST1 expression and that the risk allele has increased binding of the NOTCH signaling protein RBPJ. Finally, we provide data that TWIST1 expression influences vascular SMC phenotypes, including proliferation and calcification, as a potential mechanism supporting a role for TWIST1 in CAD.

A systematic review of the natural history and biomarkers of primary lecithin:cholesterol acyltransferase deficiency.

Syndromes associated with LCAT deficiency, a rare autosomal recessive condition, include fish-eye disease (FED) and familial LCAT deficiency (FLD). FLD is more severe and characterized by early and progressive chronic kidney disease (CKD). No treatment is currently available for FLD, but novel therapeutics are under development. Furthermore, although biomarkers of LCAT deficiency have been identified, their suitability to monitor disease progression and therapeutic efficacy is unclear, as little data exist on the rate of progression of renal disease. Here, we systematically review observational studies of FLD, FED, and heterozygous subjects, which summarize available evidence on the natural history and biomarkers of LCAT deficiency, in order to guide the development of novel therapeutics. We identified 146 FLD and 53 FED patients from 219 publications, showing that both syndromes are characterized by early corneal opacity and markedly reduced HDL-C levels. Proteinuria/hematuria were the first signs of renal impairment in FLD, followed by rapid decline of renal function. Furthermore, LCAT activity toward endogenous substrates and the percentage of circulating esterified cholesterol (EC%) were the best discriminators between these two syndromes. In FLD, higher levels of total, non-HDL, and unesterified cholesterol were associated with severe CKD. We reveal a nonlinear association between LCAT activity and EC% levels, in which subnormal levels of LCAT activity were associated with normal EC%. This review provides the first step toward the identification of disease biomarkers to be used in clinical trials and suggests that restoring LCAT activity to subnormal levels may be sufficient to prevent renal disease progression.

Prioritizing the Role of Major Lipoproteins and Subfractions as Risk Factors for Peripheral Artery Disease.

BACKGROUND: Lipoprotein-related traits have been consistently identified as risk factors for atherosclerotic cardiovascular disease, largely on the basis of studies of coronary artery disease (CAD). The relative contributions of specific lipoproteins to the risk of peripheral artery disease (PAD) have not been well defined. We leveraged large-scale genetic association data to investigate the effects of circulating lipoprotein-related traits on PAD risk. METHODS: Genome-wide association study summary statistics for circulating lipoprotein-related traits were used in the mendelian randomization bayesian model averaging framework to prioritize the most likely causal major lipoprotein and subfraction risk factors for PAD and CAD. Mendelian randomization was used to estimate the effect of apolipoprotein B (ApoB) lowering on PAD risk using gene regions proxying lipid-lowering drug targets. Genes relevant to prioritized lipoprotein subfractions were identified with transcriptome-wide association studies. RESULTS: ApoB was identified as the most likely causal lipoprotein-related risk factor for both PAD (marginal inclusion probability, 0.86; P=0.003) and CAD (marginal inclusion probability, 0.92; P=0.005). Genetic proxies for ApoB-lowering medications were associated with reduced risk of both PAD (odds ratio,0.87 per 1-SD decrease in ApoB [95% CI, 0.84-0.91]; P=9×10-10) and CAD (odds ratio,0.66 [95% CI, 0.63-0.69]; P=4×10-73), with a stronger predicted effect of ApoB lowering on CAD (ratio of effects, 3.09 [95% CI, 2.29-4.60]; P<1×10-6). Extra-small very-low-density lipoprotein particle concentration was identified as the most likely subfraction associated with PAD risk (marginal inclusion probability, 0.91; P=2.3×10-4), whereas large low-density lipoprotein particle concentration was the most likely subfraction associated with CAD risk (marginal inclusion probability, 0.95; P=0.011). Genes associated with extra-small very-low-density lipoprotein particle and large low-density lipoprotein particle concentration included canonical ApoB pathway components, although gene-specific effects were variable. Lipoprotein(a) was associated with increased risk of PAD independently of ApoB (odds ratio, 1.04 [95% CI, 1.03-1.04]; P=1.0×10-33). CONCLUSIONS: ApoB was prioritized as the major lipoprotein fraction causally responsible for both PAD and CAD risk. However, ApoB-lowering drug targets and ApoB-containing lipoprotein subfractions had diverse associations with atherosclerotic cardiovascular disease, and distinct subfraction-associated genes suggest possible differences in the role of lipoproteins in the pathogenesis of PAD and CAD.