Iron dysregulation and inflammatory stress erythropoiesis associates with long-term outcome of COVID-19.

Persistent symptoms following SARS-CoV-2 infection are increasingly reported, although the drivers of post-acute sequelae (PASC) of COVID-19 are unclear. Here we assessed 214 individuals infected with SARS-CoV-2, with varying disease severity, for one year from COVID-19 symptom onset to determine the early correlates of PASC. A multivariate signature detected beyond two weeks of disease, encompassing unresolving inflammation, anemia, low serum iron, altered iron-homeostasis gene expression and emerging stress erythropoiesis; differentiated those who reported PASC months later, irrespective of COVID-19 severity. A whole-blood heme-metabolism signature, enriched in hospitalized patients at month 1-3 post onset, coincided with pronounced iron-deficient reticulocytosis. Lymphopenia and low numbers of dendritic cells persisted in those with PASC, and single-cell analysis reported iron maldistribution, suggesting monocyte iron loading and increased iron demand in proliferating lymphocytes. Thus, defects in iron homeostasis, dysregulated erythropoiesis and immune dysfunction due to COVID-19 possibly contribute to inefficient oxygen transport, inflammatory disequilibrium and persisting symptomatology, and may be therapeutically tractable.

Detection and annotation of transposable element insertions and deletions on the human genome using nanopore sequencing.

Repetitive sequences represent about 45% of the human genome. Some are transposable elements (TEs) with the ability to change their position in the genome, creating genetic variability both as insertions or deletions, with potential pathogenic consequences. We used long-read nanopore sequencing to identify TE variants in the genomes of 24 patients with antithrombin deficiency. We identified 7 344 TE insertions and 3 056 TE deletions, 2 926 were not previously described in publicly available databases. The insertions affected 3 955 genes, with 6 insertions located in exons, 3 929 in introns, and 147 in promoters. Potential functional impact was evaluated with gene annotation and enrichment analysis, which suggested a strong relationship with neuron-related functions and autism. We conclude that this study encourages the generation of a complete map of TEs in the human genome, which will be useful for identifying new TEs involved in genetic disorders.

Genome sequencing and comprehensive rare-variant analysis of 465 families with neurodevelopmental disorders.

Despite significant progress in unraveling the genetic causes of neurodevelopmental disorders (NDDs), a substantial proportion of individuals with NDDs remain without a genetic diagnosis after microarray and/or exome sequencing. Here, we aimed to assess the power of short-read genome sequencing (GS), complemented with long-read GS, to identify causal variants in participants with NDD from the National Institute for Health and Care Research (NIHR) BioResource project. Short-read GS was conducted on 692 individuals (489 affected and 203 unaffected relatives) from 465 families. Additionally, long-read GS was performed on five affected individuals who had structural variants (SVs) in technically challenging regions, had complex SVs, or required distal variant phasing. Causal variants were identified in 36% of affected individuals (177/489), and a further 23% (112/489) had a variant of uncertain significance after multiple rounds of re-analysis. Among all reported variants, 88% (333/380) were coding nuclear SNVs or insertions and deletions (indels), and the remainder were SVs, non-coding variants, and mitochondrial variants. Furthermore, long-read GS facilitated the resolution of challenging SVs and invalidated variants of difficult interpretation from short-read GS. This study demonstrates the value of short-read GS, complemented with long-read GS, in investigating the genetic causes of NDDs. GS provides a comprehensive and unbiased method of identifying all types of variants throughout the nuclear and mitochondrial genomes in individuals with NDD.

A signature of platelet reactivity in CBC scattergrams reveals genetic predictors of thrombotic disease risk.

Genetic studies of platelet reactivity (PR) phenotypes may identify novel antiplatelet drug targets. However, such studies have been limited by small sample sizes (n < 5000) because of the complexity of measuring PR. We trained a model to predict PR from complete blood count (CBC) scattergrams. A genome-wide association study of this phenotype in 29 806 blood donors identified 21 distinct associations implicating 20 genes, of which 6 have been identified previously. The effect size estimates were significantly correlated with estimates from a study of flow cytometry-measured PR and a study of a phenotype of in vitro thrombus formation. A genetic score of PR built from the 21 variants was associated with the incidence rates of myocardial infarction and pulmonary embolism. Mendelian randomization analyses showed that PR was causally associated with the risks of coronary artery disease, stroke, and venous thromboembolism. Our approach provides a blueprint for using phenotype imputation to study the determinants of hard-to-measure but biologically important hematological traits.

Long-Read Sequencing Identifies the First Retrotransposon Insertion and Resolves Structural Variants Causing Antithrombin Deficiency.

The identification of inherited antithrombin deficiency (ATD) is critical to prevent potentially life-threatening thrombotic events. Causal variants in SERPINC1 are identified for up to 70% of cases, the majority being single-nucleotide variants and indels. The detection and characterization of structural variants (SVs) in ATD remain challenging due to the high number of repetitive elements in SERPINC1. Here, we performed long-read whole-genome sequencing on 10 familial and 9 singleton cases with type I ATD proven by functional and antigen assays, who were selected from a cohort of 340 patients with this rare disorder because genetic analyses were either negative, ambiguous, or not fully characterized. We developed an analysis workflow to identify disease-associated SVs. This approach resolved, independently of its size or type, all eight SVs detected by multiple ligation-dependent probe amplification, and identified for the first time a complex rearrangement previously misclassified as a deletion. Remarkably, we identified the mechanism explaining ATD in 2 out of 11 cases with previous unknown defect: the insertion of a novel 2.4 kb SINE-VNTR-Alu retroelement, which was characterized by de novo assembly and verified by specific polymerase chain reaction amplification and sequencing in the probands and affected relatives. The nucleotide-level resolution achieved for all SVs allowed breakpoint analysis, which revealed repetitive elements and microhomologies supporting a common replication-based mechanism for all the SVs. Our study underscores the utility of long-read sequencing technology as a complementary method to identify, characterize, and unveil the molecular mechanism of disease-causing SVs involved in ATD, and enlarges the catalogue of genetic disorders caused by retrotransposon insertions.

G protein-coupled receptor kinase 5 regulates thrombin signaling in platelets via PAR-1.

The interindividual variation in the functional response of platelets to activation by agonists is heritable. Genome-wide association studies (GWASs) of quantitative measures of platelet function have identified fewer than 20 distinctly associated variants, some with unknown mechanisms. Here, we report GWASs of pathway-specific functional responses to agonism by adenosine 5'-diphosphate, a glycoprotein VI-specific collagen mimetic, and thrombin receptor-agonist peptides, each specific to 1 of the G protein-coupled receptors PAR-1 and PAR-4, in subsets of 1562 individuals. We identified an association (P = 2.75 × 10-40) between a common intronic variant, rs10886430, in the G protein-coupled receptor kinase 5 gene (GRK5) and the sensitivity of platelets to activate through PAR-1. The variant resides in a megakaryocyte-specific enhancer that is bound by the transcription factors GATA1 and MEIS1. The minor allele (G) is associated with fewer GRK5 transcripts in platelets and the greater sensitivity of platelets to activate through PAR-1. We show that thrombin-mediated activation of human platelets causes binding of GRK5 to PAR-1 and that deletion of the mouse homolog Grk5 enhances thrombin-induced platelet activation sensitivity and increases platelet accumulation at the site of vascular injury. This corroborates evidence that the human G allele of rs10886430 is associated with a greater risk for cardiovascular disease. In summary, by combining the results of pathway-specific GWASs and expression quantitative trait locus studies in humans with the results from platelet function studies in Grk5-/- mice, we obtain evidence that GRK5 regulates the human platelet response to thrombin via the PAR-1 pathway.

Genetic determinants of ferritin, haemoglobin levels and haemoglobin trajectories: results from Donor InSight.

BACKGROUND AND OBJECTIVES: Blood donors might develop iron deficiency as approximately 250 mg of iron is lost with every donation. Susceptibility to iron deficiency and low haemoglobin levels differs between individuals, which might be due to genetic variation. Therefore, the aim of this study was to investigate associations between single nucleotide polymorphisms (SNPs) and haemoglobin trajectories, haemoglobin levels and ferritin levels in blood donors. MATERIALS AND METHODS: In 2655 donors participating in the observational cohort study Donor InSight-III (2015-2017), haemoglobin and ferritin levels were measured in venous EDTA whole blood and plasma samples, respectively. Haemoglobin trajectories (stable/declining) were determined by fitting growth-mixture models on repeated pre-donation capillary haemoglobin measurements. Genotyping was done using the UK Biobank - version 2 Axiom Array. Single SNP analyses adopting an additive genetic model on imputed genetic variants were performed for haemoglobin trajectories, haemoglobin levels and ferritin levels. Conditional analyses identified independent SNPs. RESULTS: Twelve, twenty and twenty-four independent SNPs were associated with haemoglobin trajectories, haemoglobin levels and ferritin levels respectively (P < 1 x 10-5 ). Rs112016443 reached genome-wide significance for ferritin levels, which influences WDSUB1 expression. CONCLUSION: Rs112016443 was genome-wide significantly associated with ferritin levels in Dutch donors. Further validation studies are needed, as well as studies towards underlying mechanisms and predicting iron deficiency using SNPs.

Development and validation of a universal blood donor genotyping platform: a multinational prospective study.

Each year, blood transfusions save millions of lives. However, under current blood-matching practices, sensitization to non-self-antigens is an unavoidable adverse side effect of transfusion. We describe a universal donor typing platform that could be adopted by blood services worldwide to facilitate a universal extended blood-matching policy and reduce sensitization rates. This DNA-based test is capable of simultaneously typing most clinically relevant red blood cell (RBC), human platelet (HPA), and human leukocyte (HLA) antigens. Validation was performed, using samples from 7927 European, 27 South Asian, 21 East Asian, and 9 African blood donors enrolled in 2 national biobanks. We illustrated the usefulness of the platform by analyzing antibody data from patients sensitized with multiple RBC alloantibodies. Genotyping results demonstrated concordance of 99.91%, 99.97%, and 99.03% with RBC, HPA, and HLA clinically validated typing results in 89 371, 3016, and 9289 comparisons, respectively. Genotyping increased the total number of antigen typing results available from 110 980 to >1 200 000. Dense donor typing allowed identification of 2 to 6 times more compatible donors to serve 3146 patients with multiple RBC alloantibodies, providing at least 1 match for 176 individuals for whom previously no blood could be found among the same donors. This genotyping technology is already being used to type thousands of donors taking part in national genotyping studies. Extraction of dense antigen-typing data from these cohorts provides blood supply organizations with the opportunity to implement a policy of genomics-based precision matching of blood.

Author Correction: Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Publisher Correction: Whole-genome sequencing of a sporadic primary immunodeficiency cohort.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Whole-genome sequencing of a sporadic primary immunodeficiency cohort.

Primary immunodeficiency (PID) is characterized by recurrent and often life-threatening infections, autoimmunity and cancer, and it poses major diagnostic and therapeutic challenges. Although the most severe forms of PID are identified in early childhood, most patients present in adulthood, typically with no apparent family history and a variable clinical phenotype of widespread immune dysregulation: about 25% of patients have autoimmune disease, allergy is prevalent and up to 10% develop lymphoid malignancies1-3. Consequently, in sporadic (or non-familial) PID genetic diagnosis is difficult and the role of genetics is not well defined. Here we address these challenges by performing whole-genome sequencing in a large PID cohort of 1,318 participants. An analysis of the coding regions of the genome in 886 index cases of PID found that disease-causing mutations in known genes that are implicated in monogenic PID occurred in 10.3% of these patients, and a Bayesian approach (BeviMed4) identified multiple new candidate PID-associated genes, including IVNS1ABP. We also examined the noncoding genome, and found deletions in regulatory regions that contribute to disease causation. In addition, we used a genome-wide association study to identify loci that are associated with PID, and found evidence for the colocalization of-and interplay between-novel high-penetrance monogenic variants and common variants (at the PTPN2 and SOCS1 loci). This begins to explain the contribution of common variants to the variable penetrance and phenotypic complexity that are observed in PID. Thus, using a cohort-based whole-genome-sequencing approach in the diagnosis of PID can increase diagnostic yield and further our understanding of the key pathways that influence immune responsiveness in humans.

Multiple GYPB gene deletions associated with the U- phenotype in those of African ancestry.

BACKGROUND: The MNS blood group system is defined by three homologous genes: GYPA, GYPB, and GYPE. GYPB encodes for glycophorin B (GPB) carrying S/s and the "universal" antigen U. RBCs of approximately 1% of individuals of African ancestry are U- due to absence of GPB. The U- phenotype has long been attributed to a deletion encompassing GYPB exons 2 to 5 and GYPE exon 1 (GYPB*01N). STUDY DESIGN AND METHODS: Samples from two U-individuals underwent Illumina short read whole genome sequencing (WGS) and Nanopore long read WGS. In addition, two existing WGS datasets, MedSeq (n = 110) and 1000 Genomes (1000G, n = 2535), were analyzed for GYPB deletions. Deletions were confirmed by Sanger sequencing. Twenty known U- donor samples were tested by a PCR assay to determine the specific deletion alleles present in African Americans. RESULTS: Two large GYPB deletions in U- samples of African ancestry were identified: a 110 kb deletion extending left of GYPB (DEL_B_LEFT) and a 103 kb deletion extending right (DEL_B_RIGHT). DEL_B_LEFT and DEL_B_RIGHT were the most common GYPB deletions in the 1000 Genomes Project 669 African genomes (allele frequencies 0.04 and 0.02). Seven additional deletions involving GYPB were seen in African, Admixed American, and South Asian samples. No samples analyzed had GYPB*01N. CONCLUSIONS: The U- phenotype in those of African ancestry is primarily associated with two different complete deletions of GYPB (with intact GYPE). Seven additional less common GYPB deletion backgrounds were found. GYPB*01N, long assumed to be the allele commonly encoding U- phenotypes, appears to be rare.

Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

Several strands of evidence question the dogma that human mitochondrial DNA (mtDNA) is inherited exclusively down the maternal line, most recently in three families where several individuals harbored a 'heteroplasmic haplotype' consistent with biparental transmission. Here we report a similar genetic signature in 7 of 11,035 trios, with allelic fractions of 5-25%, implying biparental inheritance of mtDNA in 0.06% of offspring. However, analysing the nuclear whole genome sequence, we observe likely large rare or unique nuclear-mitochondrial DNA segments (mega-NUMTs) transmitted from the father in all 7 families. Independently detecting mega-NUMTs in 0.13% of fathers, we see autosomal transmission of the haplotype. Finally, we show the haplotype allele fraction can be explained by complex concatenated mtDNA-derived sequences rearranged within the nuclear genome. We conclude that rare cryptic mega-NUMTs can resemble paternally mtDNA heteroplasmy, but find no evidence of paternal transmission of mtDNA in humans.

Contrasting patterns of longitudinal population dynamics and antimicrobial resistance mechanisms in two priority bacterial pathogens over 7 years in a single center.

BACKGROUND: Two of the most important pathogens contributing to the global rise in antimicrobial resistance (AMR) are Klebsiella pneumoniae and Enterobacter cloacae. Despite this, most of our knowledge about the changing patterns of disease caused by these two pathogens is based on studies with limited timeframes that provide few insights into their population dynamics or the dynamics in AMR elements that they can carry. RESULTS: We investigate the population dynamics of two priority AMR pathogens over 7 years between 2007 and 2012 in a major UK hospital, spanning changes made to UK national antimicrobial prescribing policy in 2007. Between 2006 and 2012, K. pneumoniae showed epidemiological cycles of multi-drug-resistant (MDR) lineages being replaced approximately every 2 years. This contrasted E. cloacae where there was no temporally changing pattern, but a continuous presence of the mixed population. CONCLUSIONS: The differing patterns of clonal replacement and acquisition of mobile elements shows that the flux in the K. pneumoniae population was linked to the introduction of globally recognized MDR clones carrying drug resistance markers on mobile elements. However, E. cloacae carries a chromosomally encoded ampC conferring resistance to front-line treatments and shows that MDR plasmid acquisition in E. cloacae was not indicative of success in the hospital. This led to markedly different dynamics in the AMR populations of these two pathogens and shows that the mechanism of the resistance and its location in the genome or mobile elements is crucial to predict population dynamics of opportunistic pathogens in clinical settings.

Genomic identification of cryptic susceptibility to penicillins and ß-lactamase inhibitors in methicillin-resistant Staphylococcus aureus.

Antibiotic resistance in bacterial pathogens threatens the future of modern medicine. One such resistant pathogen is methicillin-resistant Staphylococcus aureus (MRSA), which is resistant to nearly all ß-lactam antibiotics, limiting treatment options. Here, we show that a significant proportion of MRSA isolates from different lineages, including the epidemic USA300 lineage, are susceptible to penicillins when used in combination with ß-lactamase inhibitors such as clavulanic acid. Susceptibility is mediated by a combination of two different mutations in the mecA promoter region that lowers mecA-encoded penicillin-binding protein 2a (PBP2a) expression, and in the majority of isolates by either one of two substitutions in PBP2a (E246G or M122I) that increase the affinity of PBP2a for penicillin in the presence of clavulanic acid. Treatment of S. aureus infections in wax moth and mouse models shows that penicillin/ß-lactamase inhibitor susceptibility can be exploited as an effective therapeutic choice for 'susceptible' MRSA infection. Finally, we show that isolates with the PBP2a E246G substitution have a growth advantage in the presence of penicillin but the absence of clavulanic acid, which suggests that penicillin/ß-lactamase susceptibility is an example of collateral sensitivity (resistance to one antibiotic increases sensitivity to another). Our findings suggest that widely available and currently disregarded antibiotics could be effective in a significant proportion of MRSA infections.

Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

To identify novel causes of hereditary thrombocytopenia, we performed a genetic association analysis of whole-genome sequencing data from 13 037 individuals enrolled in the National Institute for Health Research (NIHR) BioResource, including 233 cases with isolated thrombocytopenia. We found an association between rare variants in the transcription factor-encoding gene IKZF5 and thrombocytopenia. We report 5 causal missense variants in or near IKZF5 zinc fingers, of which 2 occurred de novo and 3 co-segregated in 3 pedigrees. A canonical DNA-zinc finger binding model predicts that 3 of the variants alter DNA recognition. Expression studies showed that chromatin binding was disrupted in mutant compared with wild-type IKZF5, and electron microscopy revealed a reduced quantity of α granules in normally sized platelets. Proplatelet formation was reduced in megakaryocytes from 7 cases relative to 6 controls. Comparison of RNA-sequencing data from platelets, monocytes, neutrophils, and CD4+ T cells from 3 cases and 14 healthy controls showed 1194 differentially expressed genes in platelets but only 4 differentially expressed genes in each of the other blood cell types. In conclusion, IKZF5 is a novel transcriptional regulator of megakaryopoiesis and the eighth transcription factor associated with dominant thrombocytopenia in humans.

Complex structural variants in Mendelian disorders: identification and breakpoint resolution using short- and long-read genome sequencing.

BACKGROUND: Studies have shown that complex structural variants (cxSVs) contribute to human genomic variation and can cause Mendelian disease. We aimed to identify cxSVs relevant to Mendelian disease using short-read whole-genome sequencing (WGS), resolve the precise variant configuration and investigate possible mechanisms of cxSV formation. METHODS: We performed short-read WGS and analysis of breakpoint junctions to identify cxSVs in a cohort of 1324 undiagnosed rare disease patients. Long-read WGS and gene expression analysis were used to resolve one case. RESULTS: We identified three pathogenic cxSVs: a de novo duplication-inversion-inversion-deletion affecting ARID1B, a de novo deletion-inversion-duplication affecting HNRNPU and a homozygous deletion-inversion-deletion affecting CEP78. Additionally, a de novo duplication-inversion-duplication overlapping CDKL5 was resolved by long-read WGS demonstrating the presence of both a disrupted and an intact copy of CDKL5 on the same allele, and gene expression analysis showed both parental alleles of CDKL5 were expressed. Breakpoint analysis in all the cxSVs revealed both microhomology and longer repetitive elements. CONCLUSIONS: Our results corroborate that cxSVs cause Mendelian disease, and we recommend their consideration during clinical investigations. We show that resolution of breakpoints can be critical to interpret pathogenicity and present evidence of replication-based mechanisms in cxSV formation.

Automated typing of red blood cell and platelet antigens: a whole-genome sequencing study.

BACKGROUND: There are more than 300 known red blood cell (RBC) antigens and 33 platelet antigens that differ between individuals. Sensitisation to antigens is a serious complication that can occur in prenatal medicine and after blood transfusion, particularly for patients who require multiple transfusions. Although pre-transfusion compatibility testing largely relies on serological methods, reagents are not available for many antigens. Methods based on single-nucleotide polymorphism (SNP) arrays have been used, but typing for ABO and Rh-the most important blood groups-cannot be done with SNP typing alone. We aimed to develop a novel method based on whole-genome sequencing to identify RBC and platelet antigens. METHODS: This whole-genome sequencing study is a subanalysis of data from patients in the whole-genome sequencing arm of the MedSeq Project randomised controlled trial (NCT01736566) with no measured patient outcomes. We created a database of molecular changes in RBC and platelet antigens and developed an automated antigen-typing algorithm based on whole-genome sequencing (bloodTyper). This algorithm was iteratively improved to address cis-trans haplotype ambiguities and homologous gene alignments. Whole-genome sequencing data from 110 MedSeq participants (30 × depth) were used to initially validate bloodTyper through comparison with conventional serology and SNP methods for typing of 38 RBC antigens in 12 blood-group systems and 22 human platelet antigens. bloodTyper was further validated with whole-genome sequencing data from 200 INTERVAL trial participants (15 × depth) with serological comparisons. FINDINGS: We iteratively improved bloodTyper by comparing its typing results with conventional serological and SNP typing in three rounds of testing. The initial whole-genome sequencing typing algorithm was 99·5% concordant across the first 20 MedSeq genomes. Addressing discordances led to development of an improved algorithm that was 99·8% concordant for the remaining 90 MedSeq genomes. Additional modifications led to the final algorithm, which was 99·2% concordant across 200 INTERVAL genomes (or 99·9% after adjustment for the lower depth of coverage). INTERPRETATION: By enabling more precise antigen-matching of patients with blood donors, antigen typing based on whole-genome sequencing provides a novel approach to improve transfusion outcomes with the potential to transform the practice of transfusion medicine. FUNDING: National Human Genome Research Institute, Doris Duke Charitable Foundation, National Health Service Blood and Transplant, National Institute for Health Research, and Wellcome Trust.

Identification of LukPQ, a novel, equid-adapted leukocidin of Staphylococcus aureus.

Bicomponent pore-forming leukocidins are a family of potent toxins secreted by Staphylococcus aureus, which target white blood cells preferentially and consist of an S- and an F-component. The S-component recognizes a receptor on the host cell, enabling high-affinity binding to the cell surface, after which the toxins form a pore that penetrates the cell lipid bilayer. Until now, six different leukocidins have been described, some of which are host and cell specific. Here, we identify and characterise a novel S. aureus leukocidin; LukPQ. LukPQ is encoded on a 45 kb prophage (ΦSaeq1) found in six different clonal lineages, almost exclusively in strains cultured from equids. We show that LukPQ is a potent and specific killer of equine neutrophils and identify equine-CXCRA and CXCR2 as its target receptors. Although the S-component (LukP) is highly similar to the S-component of LukED, the species specificity of LukPQ and LukED differs. By forming non-canonical toxin pairs, we identify that the F-component contributes to the observed host tropism of LukPQ, thereby challenging the current paradigm that leukocidin specificity is driven solely by the S-component.