Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells.

To better understand host-virus genetic dependencies and find potential therapeutic targets for COVID-19, we performed a genome-scale CRISPR loss-of-function screen to identify host factors required for SARS-CoV-2 viral infection of human alveolar epithelial cells. Top-ranked genes cluster into distinct pathways, including the vacuolar ATPase proton pump, Retromer, and Commander complexes. We validate these gene targets using several orthogonal methods such as CRISPR knockout, RNA interference knockdown, and small-molecule inhibitors. Using single-cell RNA-sequencing, we identify shared transcriptional changes in cholesterol biosynthesis upon loss of top-ranked genes. In addition, given the key role of the ACE2 receptor in the early stages of viral entry, we show that loss of RAB7A reduces viral entry by sequestering the ACE2 receptor inside cells. Overall, this work provides a genome-scale, quantitative resource of the impact of the loss of each host gene on fitness/response to viral infection.

Reconstructing Denisovan Anatomy Using DNA Methylation Maps.

Denisovans are an extinct group of humans whose morphology remains unknown. Here, we present a method for reconstructing skeletal morphology using DNA methylation patterns. Our method is based on linking unidirectional methylation changes to loss-of-function phenotypes. We tested performance by reconstructing Neanderthal and chimpanzee skeletal morphologies and obtained >85% precision in identifying divergent traits. We then applied this method to the Denisovan and offer a putative morphological profile. We suggest that Denisovans likely shared with Neanderthals traits such as an elongated face and a wide pelvis. We also identify Denisovan-derived changes, such as an increased dental arch and lateral cranial expansion. Our predictions match the only morphologically informative Denisovan bone to date, as well as the Xuchang skull, which was suggested by some to be a Denisovan. We conclude that DNA methylation can be used to reconstruct anatomical features, including some that do not survive in the fossil record.

Reshaping endoplasmic reticulum quality control through the unfolded protein response.

Endoplasmic reticulum quality control (ERQC) pathways comprising chaperones, folding enzymes, and degradation factors ensure the fidelity of ER protein folding and trafficking to downstream secretory environments. However, multiple factors, including tissue-specific secretory proteomes, environmental and genetic insults, and organismal aging, challenge ERQC. Thus, a key question is: how do cells adapt ERQC to match the diverse, ever-changing demands encountered during normal physiology and in disease? The answer lies in the unfolded protein response (UPR), a signaling mechanism activated by ER stress. In mammals, the UPR comprises three signaling pathways regulated downstream of the ER membrane proteins IRE1, ATF6, and PERK. Upon activation, these UPR pathways remodel ERQC to alleviate cellular stress and restore ER function. Here, we describe how UPR signaling pathways adapt ERQC, highlighting their importance for maintaining ER function across tissues and the potential for targeting the UPR to mitigate pathologies associated with protein misfolding diseases.

Immune surveillance of brain metastatic cancer cells is mediated by IFITM1.

Brain metastasis, most commonly originating from lung cancer, increases cancer morbidity and mortality. Although metastatic colonization is the rate-limiting and most complex step of the metastatic cascade, the underlying mechanisms are poorly understood. Here, in vivo genome-wide CRISPR-Cas9 screening revealed that loss of interferon-induced transmembrane protein 1 (IFITM1) promotes brain colonization of human lung cancer cells. Incipient brain metastatic cancer cells with high expression of IFITM1 secrete microglia-activating complement component 3 and enhance the cytolytic activity of CD8+ T cells by increasing the expression and membrane localization of major histocompatibility complex class I. After activation, microglia (of the innate immune system) and cytotoxic CD8+ T lymphocytes (of the adaptive immune system) were found to jointly eliminate cancer cells by releasing interferon-gamma and inducing phagocytosis and T-cell-mediated killing. In human cancer clinical trials, immune checkpoint blockade therapy response was significantly correlated with IFITM1 expression, and IFITM1 enhanced the brain metastasis suppression efficacy of PD-1 blockade in mice. Our results exemplify a novel mechanism through which metastatic cancer cells overcome the innate and adaptive immune responses to colonize the brain, and suggest that a combination therapy increasing IFITM1 expression in metastatic cells with PD-1 blockade may be a promising strategy to reduce metastasis.

Protein Kinase C Quality Control by Phosphatase PHLPP1 Unveils Loss-of-Function Mechanism in Cancer.

Protein kinase C (PKC) isozymes function as tumor suppressors in increasing contexts. In contrast to oncogenic kinases, whose function is acutely regulated by transient phosphorylation, PKC is constitutively phosphorylated following biosynthesis to yield a stable, autoinhibited enzyme that is reversibly activated by second messengers. Here, we report that the phosphatase PHLPP1 opposes PKC phosphorylation during maturation, leading to the degradation of aberrantly active species that do not become autoinhibited. Cancer-associated hotspot mutations in the pseudosubstrate of PKCß that impair autoinhibition result in dephosphorylated and unstable enzymes. Protein-level analysis reveals that PKCα is fully phosphorylated at the PHLPP site in over 5,000 patient tumors, with higher PKC levels correlating (1) inversely with PHLPP1 levels and (2) positively with improved survival in pancreatic adenocarcinoma. Thus, PHLPP1 provides a proofreading step that maintains the fidelity of PKC autoinhibition and reveals a prominent loss-of-function mechanism in cancer by suppressing the steady-state levels of PKC.

Advanced variant classification framework reduces the false positive rate of predicted loss-of-function variants in population sequencing data.

Predicted loss of function (pLoF) variants are often highly deleterious and play an important role in disease biology, but many pLoF variants may not result in loss of function (LoF). Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines' PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 genes associated with autosomal-recessive disease from the Genome Aggregation Database (gnomAD v.2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in a low proportion expressed across transcripts (pext) scored region, or the presence of cryptic in-frame splice rescues. Variants predicted to evade LoF or to be potential artifacts were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of pLoF variants predicted as likely not LoF/not LoF, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

Natural selection has driven the recurrent loss of an immunity gene that protects Drosophila against a major natural parasite.

Polymorphisms in immunity genes can have large effects on susceptibility to infection. To understand the origins of this variation, we have investigated the genetic basis of resistance to the parasitoid wasp Leptopilina boulardi in Drosophila melanogaster. We found that increased expression of the gene lectin-24A after infection by parasitic wasps was associated with a faster cellular immune response and greatly increased rates of killing the parasite. lectin-24A encodes a protein that is strongly up-regulated in the fat body after infection and localizes to the surface of the parasite egg. In certain susceptible lines, a deletion upstream of the lectin-24A has largely abolished expression. Other mutations predicted to abolish the function of this gene have arisen recurrently in this gene, with multiple loss-of-expression alleles and premature stop codons segregating in natural populations. The frequency of these alleles varies greatly geographically, and in some southern African populations, natural selection has driven them near to fixation. We conclude that natural selection has favored the repeated loss of an important component of the immune system, suggesting that in some populations, a pleiotropic cost to lectin-24A expression outweighs the benefits of resistance.

Mutant analysis of Kcng4b reveals how the different functional states of the voltage-gated potassium channel regulate ear development.

The voltage gated (Kv) slow-inactivating delayed rectifier channel regulates the development of hollow organs of the zebrafish. The functional channel consists of the tetramer of electrically active Kcnb1 (Kv2.1) subunits and Kcng4b (Kv6.4) modulatory or electrically silent subunits. The two mutations in zebrafish kcng4b gene - kcng4b-C1 and kcng4b-C2 (Gasanov et al., 2021) - have been studied during ear development using electrophysiology, developmental biology and in silico structural modelling. kcng4b-C1 mutation causes a C-terminal truncation characterized by mild Kcng4b loss-of-function (LOF) manifested by failure of kinocilia to extend and formation of ectopic otoliths. In contrast, the kcng4b-C2-/- mutation causes the C-terminal domain to elongate and the ectopic seventh transmembrane (TM) domain to form, converting the intracellular C-terminus to an extracellular one. Kcng4b-C2 acts as a Kcng4b gain-of-function (GOF) allele. Otoliths fail to develop and kinocilia are reduced in kcng4b-C2-/-. These results show that different mutations of the silent subunit Kcng4 can affect the activity of the Kv channel and cause a wide range of developmental defects.

Ecological Diversification in an Adaptive Radiation of Plants: The Role of De Novo Mutation and Introgression.

Adaptive radiations are characterized by rapid ecological diversification and speciation events, leading to fuzzy species boundaries between ecologically differentiated species. Adaptive radiations are therefore key systems for understanding how species are formed and maintained, including the role of de novo mutations versus preexisting variation in ecological adaptation and the genome-wide consequences of hybridization events. For example, adaptive introgression, where beneficial alleles are transferred between lineages through hybridization, may fuel diversification in adaptive radiations and facilitate adaptation to new environments. In this study, we employed whole-genome resequencing data to investigate the evolutionary origin of hummingbird-pollinated flowers and to characterize genome-wide patterns of phylogenetic discordance and introgression in Penstemon subgenus Dasanthera, a small and diverse adaptive radiation of plants. We found that magenta hummingbird-adapted flowers have apparently evolved twice from ancestral blue-violet bee-pollinated flowers within this radiation. These shifts in flower color are accompanied by a variety of inactivating mutations to a key anthocyanin pathway enzyme, suggesting that independent de novo loss-of-function mutations underlie the parallel evolution of this trait. Although patterns of introgression and phylogenetic discordance were heterogenous across the genome, a strong effect of gene density suggests that, in general, natural selection opposes introgression and maintains genetic differentiation in gene-rich genomic regions. Our results highlight the importance of both de novo mutation and introgression as sources of evolutionary change and indicate a role for de novo mutation in driving parallel evolution in adaptive radiations.

Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites.

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

Analysis of missense variants in the human genome reveals widespread gene-specific clustering and improves prediction of pathogenicity.

We used a machine learning approach to analyze the within-gene distribution of missense variants observed in hereditary conditions and cancer. When applied to 840 genes from the ClinVar database, this approach detected a significant non-random distribution of pathogenic and benign variants in 387 (46%) and 172 (20%) genes, respectively, revealing that variant clustering is widespread across the human exome. This clustering likely occurs as a consequence of mechanisms shaping pathogenicity at the protein level, as illustrated by the overlap of some clusters with known functional domains. We then took advantage of these findings to develop a pathogenicity predictor, MutScore, that integrates qualitative features of DNA substitutions with the new additional information derived from this positional clustering. Using a random forest approach, MutScore was able to identify pathogenic missense mutations with very high accuracy, outperforming existing predictive tools, especially for variants associated with autosomal-dominant disease and cancer. Thus, the within-gene clustering of pathogenic and benign DNA changes is an important and previously underappreciated feature of the human exome, which can be harnessed to improve the prediction of pathogenicity and disambiguation of DNA variants of uncertain significance.

Response to Wyckelsma et al.: Loss of α-actinin-3 during human evolution provides superior cold resilience and muscle heat generation.

The common loss-of-function mutation R577X in the structural muscle protein ACTN3 emerged as a potential target of positive selection from early studies and has been the focus of insightful physiological work suggesting a significant impact on muscle metabolism. Adaptation to cold climates has been proposed as a key adaptive mechanism explaining its global allele frequency patterns. Here, we re-examine this hypothesis analyzing modern (n = 3,626) and ancient (n = 1,651) genomic data by using allele-frequency as well as haplotype homozygosity-based methods. The presented results are more consistent with genetic drift rather than selection in cold climates as the main driver of the ACTN3 R577X frequency distribution in human populations across the world. This Matters Arising paper is in response to Wyckelsma et al. (2021),1 published in The American Journal of Human Genetics. See also the response by Wyckelsma et al. (2022),2 published in this issue.

Biochemical characterization of the Drosophila insulin receptor kinase and longevity-associated mutants.

Drosophila melanogaster (fruit fly) insulin receptor (D-IR) is highly homologous to the human counterpart. Like the human pathway, D-IR responds to numerous insulin-like peptides to activate cellular signals that regulate growth, development, and lipid metabolism in fruit flies. Allelic mutations in the D-IR kinase domain elevate life expectancy in fruit flies. We developed a robust heterologous expression system to express and purify wild-type and longevity-associated mutant D-IR kinase domains to investigate enzyme kinetics and substrate specificities. D-IR exhibits remarkable similarities to the human insulin receptor kinase domain but diverges in substrate preferences. We show that longevity-associated mutations reduce D-IR catalytic activity. Deletion of the unique kinase insert domain portion or mutations proximal to activating tyrosines do not influence kinase activity, suggesting their potential role in substrate recruitment and downstream signaling. Through biochemical investigations, this study enhances our comprehension of D-IR's role in Drosophila physiology, complementing genetic studies and expanding our knowledge on the catalytic functions of this conserved signaling pathway.

Heterozygous variants in USP25 cause genetic generalized epilepsy.

USP25 encodes ubiquitin-specific proteases 25, a key member of deubiquitinating enzyme family and is involved in neural fate determination. Although abnormal expression in Down's syndrome was reported previously, the specific role of USP25 in human diseases has not been defined. In this study, we performed trio-based whole exome sequencing in a cohort of 319 cases (families) with generalized epilepsy of unknown etiology. Five heterozygous USP25 variants including two de novo and three co-segregated variants were determined in eight individuals affected by generalized seizures and/or febrile seizures from five unrelated families. The frequency of USP25 variants showed a significantly high aggregation in this cohort compared to the East Asian population and all populations in the gnomAD database. The mean onset ages of febrile and afebrile seizures were 10 months (infancy) and 11.8 years (juvenile), respectively. The patients achieved seizure freedom except one had occasional nocturnal seizures at the last follow-up. Two patients exhibited intellectual disability. Usp25 was ubiquitously expressed in mouse brain with two peaks on embryonic days (E14‒E16) and postnatal day 21, respectively. Similarly, USP25 expressed in fetus/early childhood stage with a second peak at approximately 12‒20 years old in human brain, consistent with the seizure onset age at infancy and juvenile in the patients. To investigate the functional impact of USP25 deficiency in vivo, we established Usp25 knock-out mice, which showed increased seizure susceptibility compared to wild-type mice in pentylenetetrazol-induced seizure test. To explore the impact of USP25 variants, we employed multiple functional detections. In HEK293T cells, the severe phenotype associated variant (p.Gln889Ter) led to a significant reduction of mRNA and protein expressions but formed a stable truncated dimers with increment of deubiquitinating enzyme activities and abnormal cellular aggregations, indicating a gain-of-function effect. The p.Gln889Ter and p.Leu1045del increased neuronal excitability in mice brain, with a higher firing ability in p.Gln889Ter. These functional impairments align with the severity of the observed phenotypes, suggesting a genotype-phenotype correlation. Hence, a moderate association between USP25 and epilepsy was noted, indicating USP25 is potentially a predisposing gene for epilepsy. Our results from Usp25 null mice and the patient-derived variants indicated that USP25 would play epileptogenic role via loss-of-function or gain-of-function effects. The truncated variant p.Gln889Ter would have profoundly different effect on epilepsy. Together, our results underscore the significance of USP25 heterozygous variants in epilepsy, thereby highlighting the critical role of USP25 in the brain.

Rescuing tri-heteromeric NMDA receptor function: the potential of pregnenolone-sulfate in loss-of-function GRIN2B variants.

N-methyl-D-aspartate receptors (NMDARs emerging from GRIN genes) are tetrameric receptors that form diverse channel compositions in neurons, typically consisting of two GluN1 subunits combined with two GluN2(A-D) subunits. During prenatal stages, the predominant channels are di-heteromers with two GluN1 and two GluN2B subunits due to the high abundance of GluN2B subunits. Postnatally, the expression of GluN2A subunits increases, giving rise to additional subtypes, including GluN2A-containing di-heteromers and tri-heteromers with GluN1, GluN2A, and GluN2B subunits. The latter  emerge as the major receptor subtype at mature synapses in the hippocampus. Despite extensive research on purely di-heteromeric receptors containing two identical GRIN variants, the impact of a single variant on the function of other channel forms, notably tri-heteromers, is lagging. In this study, we systematically investigated the effects of two de novo GRIN2B variants (G689C and G689S) in pure, mixed di- and tri-heteromers. Our findings reveal that incorporating a single variant in mixed di-heteromers or tri-heteromers exerts a dominant negative effect on glutamate potency, although 'mixed' channels show improved potency compared to pure variant-containing di-heteromers. We show that a single variant within a receptor complex does not impair the response of all receptor subtypes to the positive allosteric modulator pregnenolone-sulfate (PS), whereas spermine completely fails to potentiate tri-heteromers containing GluN2A and -2B-subunits. We examined PS on primary cultured hippocampal neurons transfected with the variants, and observed a positive impact over current amplitudes and synaptic activity. Together, our study supports previous observations showing that mixed di-heteromers exhibit improved glutamate potency and extend these findings towards the exploration of the effect of Loss-of-Function variants over tri-heteromers. Notably, we provide an initial and crucial demonstration of the beneficial effects of GRIN2B-relevant potentiators on tri-heteromers. Our results underscore the significance of studying how different variants affect distinct receptor subtypes, as these effects cannot be inferred solely from observations made on pure di-heteromers. Overall, this study contributes to ongoing efforts to understand the pathophysiology of GRINopathies and provides insights into potential treatment strategies.

TBX20 loss-of-function variants in families with left ventricular non-compaction cardiomyopathy.

TBX20 encodes a cardiac transcription factor that is associated with atrial septal defects. Recent studies implicate loss-of-function TBX20 variants with left ventricular non-compaction cardiomyopathy (LVNC), although clinical and genetic data in families are limited. We report four families with TBX20 loss-of-function variants that segregate with LVNC. Genetic testing using genome or exome sequencing was performed in index cases, variants were validated with Sanger sequencing, and cascade genetic testing was performed in family members. A multi-exon deletion, small deletion, essential splice site variant and nonsense variant in TBX20 were found in four families. The index cases in two families were symptomatic children with identical congenital heart diseases and LVNC who developed different cardiomyopathy phenotypes with one developing heart failure requiring transplantation. In another family, the child index case had LVNC and congestive heart failure requiring heart transplantation. In the fourth family, the index case was a symptomatic adult with LVNC. In all families, the variants segregated in relatives with isolated LVNC, or with congenital heart disease or cardiomyopathy. Family members displayed a clinical spectrum from asymptomatic to severe presentations including heart failure. Our data strengthen TBX20 loss-of-function variants as a rare cause of LVNC and support TBX20 inclusion in genetic testing of LVNC.

Expanding the genetic and clinical spectrum of Tatton-Brown-Rahman syndrome in a series of 24 French patients.

BACKGROUND: Tatton-Brown-Rahman syndrome (TBRS; OMIM 615879), also known as DNA methyltransferase 3 alpha (DNMT3A)-overgrowth syndrome (DOS), was first described by Tatton-Brown in 2014. This syndrome is characterised by overgrowth, intellectual disability and distinctive facial features and is the consequence of germline loss-of-function variants in DNMT3A, which encodes a DNA methyltransferase involved in epigenetic regulation. Somatic variants of DNMT3A are frequently observed in haematological malignancies, including acute myeloid leukaemia (AML). To date, 100 individuals with TBRS with de novo germline variants have been described. We aimed to further characterise this disorder clinically and at the molecular level in a nationwide series of 24 French patients and to investigate the correlation between the severity of intellectual disability and the type of variant. METHODS: We collected genetic and medical information from 24 individuals with TBRS using a questionnaire released through the French National AnDDI-Rares Network. RESULTS: Here, we describe the first nationwide French cohort of 24 individuals with germline likely pathogenic/pathogenic variants in DNMT3A, including 17 novel variants. We confirmed that the main phenotypic features were intellectual disability (100% of individuals), distinctive facial features (96%) and overgrowth (87%). We highlighted novel clinical features, such as hypertrichosis, and further described the neurological features and EEG results. CONCLUSION: This study of a nationwide cohort of individuals with TBRS confirms previously published data and provides additional information and clarifies clinical features to facilitate diagnosis and improve care. This study adds value to the growing body of knowledge on TBRS and broadens its clinical and molecular spectrum.

Dominant negative variants and cotranslational assembly of macromolecular complexes.

Pathogenic variants occurring in protein-coding regions underlie human genetic disease through various mechanisms. They can lead to a loss of function (LOF) such as in recessive conditions or in dominant conditions due to haploinsufficiency. Dominant-negative (DN) effects, counteracting the activity of the normal gene-product, and gain of function (GOF) are also mechanisms driving dominance. Here, I discuss a few papers on these specific mechanisms. In short, there is accumulating evidence pointing to differences between LOF versus non-LOF variants (DN and GOF). The latter are thought to have milder effects on protein structure and, as expected, DN variants are enriched at protein interfaces. This tendency to cluster in 3D space can help improve the ability of computational tools to predict the pathogenicity of DN variants, which is currently a challenging issue. More recent results support the hypothesis whereby cotranslational assembly of macromolecular complexes can buffer deleterious consequences of variants that would otherwise lead to DN effects (DNEs). Indeed, subunits the variants of which are responsible for DNEs tend to elude cotranslational assembly, thus poisoning complexes involving wild-type subunits. The constraints explaining why the buffering of DNEs is not universal require further investigation.

Loss of the TRPM4 channel in humans causes immune dysregulation with defective monocyte migration.

BACKGROUND: TRPM4 is a broadly expressed, calcium-activated, monovalent cation channel that regulates immune cell function in mice and cell lines. Clinically, however, partial loss- or gain-of-function mutations in TRPM4 lead to arrhythmia and heart disease, with no documentation of immunologic disorders. OBJECTIVE: To characterize functional cellular mechanisms underlying the immune dysregulation phenotype in a proband with a mutated TRPM4 gene. METHODS: We employed a combination of biochemical, cell biological, imaging, omics analyses, flow cytometry, and gene editing approaches. RESULTS: We report the first human cases to our knowledge with complete loss of the TRPM4 channel, leading to immune dysregulation with frequent bacterial and fungal infections. Single-cell and bulk RNA sequencing point to altered expression of genes affecting cell migration, specifically in monocytes. Inhibition of TRPM4 in T cells and the THP-1 monocyte cell line reduces migration. More importantly, primary T cells and monocytes from TRPM4 patients migrate poorly. Finally, CRISPR knockout of TRPM4 in THP-1 cells greatly reduces their migration potential. CONCLUSION: Our results demonstrate that TRPM4 plays a critical role in regulating immune cell migration, leading to increased susceptibility to infections.

Efficient identification of trait-associated loss-of-function variants in the UK Biobank cohort by exome-sequencing based genotype imputation.

The large-scale open access whole-exome sequencing (WES) data of the UK Biobank ~200,000 participants is accelerating a new wave of genetic association studies aiming to identify rare and functional loss-of-function (LoF) variants associated with complex traits and diseases. We proposed to merge the WES genotypes and the genome-wide genotyping (GWAS) genotypes of 167,000 UKB homogeneous European participants into a combined reference panel, and then to impute 241,911 UKB homogeneous European participants who had the GWAS genotypes only. We then used the imputed data to replicate association identified in the discovery WES sample. The average imputation accuracy measure r2 is modest to high for LoF variants at all minor allele frequency intervals: 0.942 at MAF interval (0.01, 0.5), 0.807 at (1.0 × 10-3 , 0.01), 0.805 at (1.0 × 10-4 , 1.0 × 10-3 ), 0.664 at (1.0 × 10-5 , 1.0 × 10-4 ) and 0.410 at (0, 1.0 × 10-5 ). As applications, we studied associations of LoF variants with estimated heel BMD and four lipid traits. In addition to replicating dozens of previously reported genes, we also identified three novel associations, two genes PLIN1 and ANGPTL3 for high-density-lipoprotein cholesterol and one gene PDE3B for triglycerides. Our results highlighted the strength of WES based genotype imputation as well as provided useful imputed data within the UKB cohort.