Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

The implementation of high-throughput sequencing (HTS) technologies in research and diagnostic laboratories has linked many new genes to rare bleeding, thrombotic, and platelet disorders (BTPD), and revealed multiple genetic variants linked to those disorders, many of them being of uncertain pathogenicity when considering the accepted evidence (variant consequence, frequency in control datasets, number of reported patients, prediction models, and functional assays). The sequencing effort has also resulted in resources for gathering disease-causing variants associated with specific genes, but for BTPD, such well-curated databases exist only for a few genes. On the other hand, submissions by individuals or diagnostic laboratories to the variant database ClinVar are hampered by the lack of a submission process tailored to capture the specific features of hemostatic diseases. As we move toward the implementation of HTS in the diagnosis of BTPD, the Scientific and Standardization Committee for Genetics in Thrombosis and Haemostasis has developed and tested a REDCap-based interface, aimed at the community, to submit curated genetic variants for diagnostic-grade BTPD genes. Here, we describe the use of the interface and the initial submission of 821 variants from 30 different centers covering 14 countries. This open-access variant resource will be shared with the community to improve variant classification and regular bulk data transfer to ClinVar.

Missense and truncating variants in the X-chromosome-linked CLCN4 gene, resulting in reduced or complete loss-of-function (LOF) of the encoded chloride/proton exchanger ClC-4, were recently demonstrated to cause a neurocognitive phenotype in both males and females. Through international clinical matchmaking and interrogation of public variant databases we assembled a database of 90 rare CLCN4 missense variants in 90 families: 41 unique and 18 recurrent variants in 49 families. For 43 families, including 22 males and 33 females, we collated detailed clinical and segregation data. To confirm causality of variants and to obtain insight into disease mechanisms, we investigated the effect on electrophysiological properties of 59 of the variants in Xenopus oocytes using extended voltage and pH ranges. Detailed analyses revealed new pathophysiological mechanisms: 25% (15/59) of variants demonstrated LOF, characterized by a "shift" of the voltage-dependent activation to more positive voltages, and nine variants resulted in a toxic gain-of-function, associated with a disrupted gate allowing inward transport at negative voltages. Functional results were not always in line with in silico pathogenicity scores, highlighting the complexity of pathogenicity assessment for accurate genetic counselling. The complex neurocognitive and psychiatric manifestations of this condition, and hitherto under-recognized impacts on growth, gastrointestinal function, and motor control are discussed. Including published cases, we summarize features in 122 individuals from 67 families with CLCN4-related neurodevelopmental condition and suggest future research directions with the aim of improving the integrated care for individuals with this diagnosis.

No abstract available

N-alpha-acetylation is a common co-translational protein modification that is essential for normal cell function in humans. We previously identified the genetic basis of an X-linked infantile lethal Mendelian disorder involving a c.109T>C (p.Ser37Pro) missense variant in NAA10, which encodes the catalytic subunit of the N-terminal acetyltransferase A (NatA) complex. The auxiliary subunit of the NatA complex, NAA15, is the dimeric binding partner for NAA10. Through a genotype-first approach with whole-exome or genome sequencing (WES/WGS) and targeted sequencing analysis, we identified and phenotypically characterized 38 individuals from 33 unrelated families with 25 different de novo or inherited, dominantly acting likely gene disrupting (LGD) variants in NAA15. Clinical features of affected individuals with LGD variants in NAA15 include variable levels of intellectual disability, delayed speech and motor milestones, and autism spectrum disorder. Additionally, mild craniofacial dysmorphology, congenital cardiac anomalies, and seizures are present in some subjects. RNA analysis in cell lines from two individuals showed degradation of the transcripts with LGD variants, probably as a result of nonsense-mediated decay. Functional assays in yeast confirmed a deleterious effect for two of the LGD variants in NAA15. Further supporting a mechanism of haploinsufficiency, individuals with copy-number variant (CNV) deletions involving NAA15 and surrounding genes can present with mild intellectual disability, mild dysmorphic features, motor delays, and decreased growth. We propose that defects in NatA-mediated N-terminal acetylation (NTA) lead to variable levels of neurodevelopmental disorders in humans, supporting the importance of the NatA complex in normal human development.

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.

Rationale: Despite the increased recognition of TBX4 (T-BOX transcription factor 4)-associated pulmonary arterial hypertension (PAH), genotype-phenotype associations are lacking and may provide important insights. Objectives: To compile and functionally characterize all TBX4 variants reported to date and undertake a comprehensive genotype-phenotype analysis. Methods: We assembled a multicenter cohort of 137 patients harboring monoallelic TBX4 variants and assessed the pathogenicity of missense variation (n = 42) using a novel luciferase reporter assay containing T-BOX binding motifs. We sought genotype-phenotype correlations and undertook a comparative analysis with patients with PAH with BMPR2 (Bone Morphogenetic Protein Receptor type 2) causal variants (n = 162) or no identified variants in PAH-associated genes (n = 741) genotyped via the National Institute for Health Research BioResource-Rare Diseases. Measurements and Main Results: Functional assessment of TBX4 missense variants led to the novel finding of gain-of-function effects associated with older age at diagnosis of lung disease compared with loss-of-function effects (P = 0.038). Variants located in the T-BOX and nuclear localization domains were associated with earlier presentation (P = 0.005) and increased incidence of interstitial lung disease (P = 0.003). Event-free survival (death or transplantation) was shorter in the T-BOX group (P = 0.022), although age had a significant effect in the hazard model (P = 0.0461). Carriers of TBX4 variants were diagnosed at a younger age (P < 0.001) and had worse baseline lung function (FEV1, FVC) (P = 0.009) than the BMPR2 and no identified causal variant groups. Conclusions: We demonstrated that TBX4 syndrome is not strictly the result of haploinsufficiency but can also be caused by gain of function. The pleiotropic effects of TBX4 in lung disease may be in part explained by the differential effect of pathogenic mutations located in critical protein domains.

Whole-genome sequencing (WGS) of cancers is becoming an accepted component of oncological care, and NHS England is currently rolling out WGS for all children with cancer. This approach was piloted during the 100,000 genomes (100 K) project. Here we share the experience of the East of England Genomic Medicine Centre (East-GMC), reporting the feasibility and clinical utility of centralised WGS for individual children locally.

Inherited retinal disease is a common cause of visual impairment and represents a highly heterogeneous group of conditions. Here, we present findings from a cohort of 722 individuals with inherited retinal disease, who have had whole-genome sequencing (n = 605), whole-exome sequencing (n = 72), or both (n = 45) performed, as part of the NIHR-BioResource Rare Diseases research study. We identified pathogenic variants (single-nucleotide variants, indels, or structural variants) for 404/722 (56%) individuals. Whole-genome sequencing gives unprecedented power to detect three categories of pathogenic variants in particular: structural variants, variants in GC-rich regions, which have significantly improved coverage compared to whole-exome sequencing, and variants in non-coding regulatory regions. In addition to previously reported pathogenic regulatory variants, we have identified a previously unreported pathogenic intronic variant in CHM in two males with choroideremia. We have also identified 19 genes not previously known to be associated with inherited retinal disease, which harbor biallelic predicted protein-truncating variants in unsolved cases. Whole-genome sequencing is an increasingly important comprehensive method with which to investigate the genetic causes of inherited retinal disease.

Although best practices have emerged on how to analyse and interpret personal genomes, the utility of whole genome screening remains underdeveloped. A large amount of information can be gathered from various types of analyses via whole genome sequencing including pathogenicity screening, genetic risk scoring, fitness, nutrition, and pharmacogenomic analysis. We recognize different levels of confidence when assessing the validity of genetic markers and apply rigorous standards for evaluation of phenotype associations. We illustrate the application of this approach on a family of five. By applying analyses of whole genomes from different methodological perspectives, we are able to build a more comprehensive picture to assist decision making in preventative healthcare and well-being management. Our interpretation and reporting outputs provide input for a clinician to develop a healthcare plan for the individual, based on genetic and other healthcare data.

The aims of our study were to use whole genome sequencing in a cross-sectional cohort of patients to identify new variants in genes implicated in neuropathic pain, to determine the prevalence of known pathogenic variants and to understand the relationship between pathogenic variants and clinical presentation. Patients with extreme neuropathic pain phenotypes (both sensory loss and gain) were recruited from secondary care clinics in the UK and underwent whole genome sequencing as part of the National Institute for Health and Care Research Bioresource Rare Diseases project. A multidisciplinary team assessed the pathogenicity of rare variants in genes previously known to cause neuropathic pain disorders and exploratory analysis of research candidate genes was completed. Association testing for genes carrying rare variants was completed using the gene-wise approach of the combined burden and variance-component test SKAT-O. Patch clamp analysis was performed on transfected HEK293T cells for research candidate variants of genes encoding ion channels. The results include the following: (i) Medically actionable variants were found in 12% of study participants (205 recruited), including known pathogenic variants: SCN9A(ENST00000409672.1): c.2544T>C, p.Ile848Thr that causes inherited erythromelalgia, and SPTLC1(ENST00000262554.2):c.340T>G, p.Cys133Tr variant that causes hereditary sensory neuropathy type-1. (ii) Clinically relevant variants were most common in voltage-gated sodium channels (Nav). (iii) SCN9A(ENST00000409672.1):c.554G>A, pArg185His variant was more common in non-freezing cold injury participants than controls and causes a gain of function of NaV1.7 after cooling (the environmental trigger for non-freezing cold injury). (iv) Rare variant association testing showed a significant difference in distribution for genes NGF, KIF1A, SCN8A, TRPM8, KIF1A, TRPA1 and the regulatory regions of genes SCN11A, FLVCR1, KIF1A and SCN9A between European participants with neuropathic pain and controls. (v) The TRPA1(ENST00000262209.4):c.515C>T, p.Ala172Val variant identified in participants with episodic somatic pain disorder demonstrated gain-of-channel function to agonist stimulation. Whole genome sequencing identified clinically relevant variants in over 10% of participants with extreme neuropathic pain phenotypes. The majority of these variants were found in ion channels. Combining genetic analysis with functional validation can lead to a better understanding as to how rare variants in ion channels lead to sensory neuron hyper-excitability, and how cold, as an environmental trigger, interacts with the gain-of-function NaV1.7 p.Arg185His variant. Our findings highlight the role of ion channel variants in the pathogenesis of extreme neuropathic pain disorders, likely mediated through changes in sensory neuron excitability and interaction with environmental triggers.

KAT6B sequence variants have been identified previously in both patients with the Say-Barber-Biesecker type of blepharophimosis mental retardation syndromes (SBBS) and in the more severe genitopatellar syndrome (GPS). We report on the findings in a previously unreported group of 57 individuals with suggestive features of SBBS or GPS. Likely causative variants have been identified in 34/57 patients and were commonly located in the terminal exons of KAT6B. Of those where parental samples could be tested, all occurred de novo. Thirty out of thirty-four had truncating variants, one had a missense variant and the remaining three had the same synonymous change predicted to affect splicing. Variants in GPS tended to occur more proximally to those in SBBS patients, and genotype/phenotype analysis demonstrated significant clinical overlap between SBBS and GPS. The de novo synonymous change seen in three patients with features of SBBS occurred more proximally in exon 16. Statistical analysis of clinical features demonstrated that KAT6B variant-positive patients were more likely to display hypotonia, feeding difficulties, long thumbs/great toes and dental, thyroid and patella abnormalities than KAT6B variant-negative patients. The few reported patients with KAT6B haploinsufficiency had a much milder phenotype, though with some features overlapping those of SBBS. We report the findings in a previously unreported patient with a deletion of the KAT6B gene to further delineate the haploinsufficiency phenotype. The molecular mechanisms giving rise to the SBBS and GPS phenotypes are discussed.

Witteveen-Kolk syndrome (OMIM 613406) is a recently defined neurodevelopmental syndrome caused by heterozygous loss-of-function variants in SIN3A. We define the clinical and neurodevelopmental phenotypes related to SIN3A-haploinsufficiency in 28 unreported patients. Patients with SIN3A variants adversely affecting protein function have mild intellectual disability, growth and feeding difficulties. Involvement of a multidisciplinary team including a geneticist, paediatrician and neurologist should be considered in managing these patients. Patients described here were identified through a combination of clinical evaluation and gene matching strategies (GeneMatcher and Decipher). All patients consented to participate in this study. Mean age of this cohort was 8.2 years (17 males, 11 females). Out of 16 patients ≥ 8 years old assessed, eight (50%) had mild intellectual disability (ID), four had moderate ID (22%), and one had severe ID (6%). Four (25%) did not have any cognitive impairment. Other neurological symptoms such as seizures (4/28) and hypotonia (12/28) were common. Behaviour problems were reported in a minority. In patients ≥2 years, three were diagnosed with Autism Spectrum Disorder (ASD) and four with Attention Deficit Hyperactivity Disorder (ADHD). We report 27 novel variants and one previously reported variant. 24 were truncating variants; three were missense variants and one large in-frame gain including exons 10-12.

ACTB encodes β-actin, an abundant cytoskeletal housekeeping protein. In humans, postulated gain-of-function missense mutations cause Baraitser-Winter syndrome (BRWS), characterized by intellectual disability, cortical malformations, coloboma, sensorineural deafness, and typical facial features. To date, the consequences of loss-of-function ACTB mutations have not been proven conclusively. We describe heterozygous ACTB deletions and nonsense and frameshift mutations in 33 individuals with developmental delay, apparent intellectual disability, increased frequency of internal organ malformations (including those of the heart and the renal tract), growth retardation, and a recognizable facial gestalt (interrupted wavy eyebrows, dense eyelashes, wide nose, wide mouth, and a prominent chin) that is distinct from characteristics of individuals with BRWS. Strikingly, this spectrum overlaps with that of several chromatin-remodeling developmental disorders. In wild-type mouse embryos, β-actin expression was prominent in the kidney, heart, and brain. ACTB mRNA expression levels in lymphoblastic lines and fibroblasts derived from affected individuals were decreased in comparison to those in control cells. Fibroblasts derived from an affected individual and ACTB siRNA knockdown in wild-type fibroblasts showed altered cell shape and migration, consistent with known roles of cytoplasmic β-actin. We also demonstrate that ACTB haploinsufficiency leads to reduced cell proliferation, altered expression of cell-cycle genes, and decreased amounts of nuclear, but not cytoplasmic, β-actin. In conclusion, we show that heterozygous loss-of-function ACTB mutations cause a distinct pleiotropic malformation syndrome with intellectual disability. Our biological studies suggest that a critically reduced amount of this protein alters cell shape, migration, proliferation, and gene expression to the detriment of brain, heart, and kidney development.

To determine the frequency of rare and pertinent disease-causing variants in small vessel disease (SVD)-associated genes (such as NOTCH3, HTRA1, COL4A1, COL4A2, FOXC1, TREX1, and GLA) in cerebral SVD, we performed targeted gene sequencing in 950 patients with younger-onset apparently sporadic SVD stroke using a targeted sequencing panel.

Bromodomain PHD finger transcription factor (BPTF) is the largest subunit of nucleosome remodeling factor (NURF), a member of the ISWI chromatin-remodeling complex. However, the clinical consequences of disruption of this complex remain largely uncharacterized. BPTF is required for anterior-posterior axis formation of the mouse embryo and was shown to promote posterior neuroectodermal fate by enhancing Smad2-activated wnt8 expression in zebrafish. Here, we report eight loss-of-function and two missense variants (eight de novo and two of unknown origin) in BPTF on 17q24.2. The BPTF variants were found in unrelated individuals aged between 2.1 and 13 years, who manifest variable degrees of developmental delay/intellectual disability (10/10), speech delay (10/10), postnatal microcephaly (7/9), and dysmorphic features (9/10). Using CRISPR-Cas9 genome editing of bptf in zebrafish to induce a loss of gene function, we observed a significant reduction in head size of F0 mutants compared to control larvae. Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) and phospho-histone H3 (PH3) staining to assess apoptosis and cell proliferation, respectively, showed a significant increase in cell death in F0 mutants compared to controls. Additionally, we observed a substantial increase of the ceratohyal angle of the craniofacial skeleton in bptf F0 mutants, indicating abnormal craniofacial patterning. Taken together, our data demonstrate the pathogenic role of BPTF haploinsufficiency in syndromic neurodevelopmental anomalies and extend the clinical spectrum of human disorders caused by ablation of chromatin remodeling complexes.

Variants in genes encoding nuclear pore complex (NPC) proteins are a newly identified cause of paediatric steroid-resistant nephrotic syndrome (SRNS). Recent reports describing NUP93 variants suggest these could be a significant cause of paediatric onset SRNS. We report NUP93 cases in the UK and demonstrate in vivo functional effects of Nup93 depletion in a fly (Drosophila melanogaster) nephrocyte model.

Genomic rearrangements involving AUTS2 (7q11.22) are associated with autism and intellectual disability (ID), although evidence for causality is limited. By combining the results of diagnostic testing of 49,684 individuals, we identified 24 microdeletions that affect at least one exon of AUTS2, as well as one translocation and one inversion each with a breakpoint within the AUTS2 locus. Comparison of 17 well-characterized individuals enabled identification of a variable syndromic phenotype including ID, autism, short stature, microcephaly, cerebral palsy, and facial dysmorphisms. The dysmorphic features were more pronounced in persons with 3'AUTS2 deletions. This part of the gene is shown to encode a C-terminal isoform (with an alternative transcription start site) expressed in the human brain. Consistent with our genetic data, suppression of auts2 in zebrafish embryos caused microcephaly that could be rescued by either the full-length or the C-terminal isoform of AUTS2. Our observations demonstrate a causal role of AUTS2 in neurocognitive disorders, establish a hitherto unappreciated syndromic phenotype at this locus, and show how transcriptional complexity can underpin human pathology. The zebrafish model provides a valuable tool for investigating the etiology of AUTS2 syndrome and facilitating gene-function analysis in the future.

The ArrayExpress Archive of Functional Genomics Data (http://www.ebi.ac.uk/arrayexpress) is an international functional genomics database at the European Bioinformatics Institute (EMBL-EBI) recommended by most journals as a repository for data supporting peer-reviewed publications. It contains data from over 7000 public sequencing and 42,000 array-based studies comprising over 1.5 million assays in total. The proportion of sequencing-based submissions has grown significantly over the last few years and has doubled in the last 18 months, whilst the rate of microarray submissions is growing slightly. All data in ArrayExpress are available in the MAGE-TAB format, which allows robust linking to data analysis and visualization tools and standardized analysis. The main development over the last two years has been the release of a new data submission tool Annotare, which has reduced the average submission time almost 3-fold. In the near future, Annotare will become the only submission route into ArrayExpress, alongside MAGE-TAB format-based pipelines. ArrayExpress is a stable and highly accessed resource. Our future tasks include automation of data flows and further integration with other EMBL-EBI resources for the representation of multi-omics data.

Multiple primary tumors (MPTs) affect a substantial proportion of cancer survivors and can result from various causes, including inherited predisposition. Currently, germline genetic testing of MPT-affected individuals for variants in cancer-predisposition genes (CPGs) is mostly targeted by tumor type. We ascertained pre-assessed MPT individuals (with at least two primary tumors by age 60 years or at least three by 70 years) from genetics centers and performed whole-genome sequencing (WGS) on 460 individuals from 440 families. Despite previous negative genetic assessment and molecular investigations, pathogenic variants in moderate- and high-risk CPGs were detected in 67/440 (15.2%) probands. WGS detected variants that would not be (or were not) detected by targeted resequencing strategies, including low-frequency structural variants (6/440 [1.4%] probands). In most individuals with a germline variant assessed as pathogenic or likely pathogenic (P/LP), at least one of their tumor types was characteristic of variants in the relevant CPG. However, in 29 probands (42.2% of those with a P/LP variant), the tumor phenotype appeared discordant. The frequency of individuals with truncating or splice-site CPG variants and at least one discordant tumor type was significantly higher than in a control population (χ2 = 43.642; p ≤ 0.0001). 2/67 (3%) probands with P/LP variants had evidence of multiple inherited neoplasia allele syndrome (MINAS) with deleterious variants in two CPGs. Together with variant detection rates from a previous series of similarly ascertained MPT-affected individuals, the present results suggest that first-line comprehensive CPG analysis in an MPT cohort referred to clinical genetics services would detect a deleterious variant in about a third of individuals.