Bi-allelic genetic variants in the translational GTPases GTPBP1 and GTPBP2 cause a distinct identical neurodevelopmental syndrome.

The homologous genes GTPBP1 and GTPBP2 encode GTP-binding proteins 1 and 2, which are involved in ribosomal homeostasis. Pathogenic variants in GTPBP2 were recently shown to be an ultra-rare cause of neurodegenerative or neurodevelopmental disorders (NDDs). Until now, no human phenotype has been linked to GTPBP1. Here, we describe individuals carrying bi-allelic GTPBP1 variants that display an identical phenotype with GTPBP2 and characterize the overall spectrum of GTP-binding protein (1/2)-related disorders. In this study, 20 individuals from 16 families with distinct NDDs and syndromic facial features were investigated by whole-exome (WES) or whole-genome (WGS) sequencing. To assess the functional impact of the identified genetic variants, semi-quantitative PCR, western blot, and ribosome profiling assays were performed in fibroblasts from affected individuals. We also investigated the effect of reducing expression of CG2017, an ortholog of human GTPBP1/2, in the fruit fly Drosophila melanogaster. Individuals with bi-allelic GTPBP1 or GTPBP2 variants presented with microcephaly, profound neurodevelopmental impairment, pathognomonic craniofacial features, and ectodermal defects. Abnormal vision and/or hearing, progressive spasticity, choreoathetoid movements, refractory epilepsy, and brain atrophy were part of the core phenotype of this syndrome. Cell line studies identified a loss-of-function (LoF) impact of the disease-associated variants but no significant abnormalities on ribosome profiling. Reduced expression of CG2017 isoforms was associated with locomotor impairment in Drosophila. In conclusion, bi-allelic GTPBP1 and GTPBP2 LoF variants cause an identical, distinct neurodevelopmental syndrome. Mutant CG2017 knockout flies display motor impairment, highlighting the conserved role for GTP-binding proteins in CNS development across species.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 macrodomain of SARS-CoV-2.

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small-molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic, there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high-resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 153 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated conformational changes within the active site, and key inhibitor motifs that will template future drug development against Mac1.

Germline MBD4 deficiency causes a multi-tumor predisposition syndrome.

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

Non-coding region variants upstream of MEF2C cause severe developmental disorder through three distinct loss-of-function mechanisms.

Clinical genetic testing of protein-coding regions identifies a likely causative variant in only around half of developmental disorder (DD) cases. The contribution of regulatory variation in non-coding regions to rare disease, including DD, remains very poorly understood. We screened 9,858 probands from the Deciphering Developmental Disorders (DDD) study for de novo mutations in the 5' untranslated regions (5' UTRs) of genes within which variants have previously been shown to cause DD through a dominant haploinsufficient mechanism. We identified four single-nucleotide variants and two copy-number variants upstream of MEF2C in a total of ten individual probands. We developed multiple bespoke and orthogonal experimental approaches to demonstrate that these variants cause DD through three distinct loss-of-function mechanisms, disrupting transcription, translation, and/or protein function. These non-coding region variants represent 23% of likely diagnoses identified in MEF2C in the DDD cohort, but these would all be missed in standard clinical genetics approaches. Nonetheless, these variants are readily detectable in exome sequence data, with 30.7% of 5' UTR bases across all genes well covered in the DDD dataset. Our analyses show that non-coding variants upstream of genes within which coding variants are known to cause DD are an important cause of severe disease and demonstrate that analyzing 5' UTRs can increase diagnostic yield. We also show how non-coding variants can help inform both the disease-causing mechanism underlying protein-coding variants and dosage tolerance of the gene.

A human importin-ß-related disorder: Syndromic thoracic aortic aneurysm caused by bi-allelic loss-of-function variants in IPO8.

Importin 8, encoded by IPO8, is a ubiquitously expressed member of the importin-ß protein family that translocates cargo molecules such as proteins, RNAs, and ribonucleoprotein complexes into the nucleus in a RanGTP-dependent manner. Current knowledge of the cargoes of importin 8 is limited, but TGF-ß signaling components such as SMAD1-4 have been suggested to be among them. Here, we report that bi-allelic loss-of-function variants in IPO8 cause a syndromic form of thoracic aortic aneurysm (TAA) with clinical overlap with Loeys-Dietz and Shprintzen-Goldberg syndromes. Seven individuals from six unrelated families showed a consistent phenotype with early-onset TAA, motor developmental delay, connective tissue findings, and craniofacial dysmorphic features. A C57BL/6N Ipo8 knockout mouse model recapitulates TAA development from 8-12 weeks onward in both sexes but most prominently shows ascending aorta dilatation with a propensity for dissection in males. Compliance assays suggest augmented passive stiffness of the ascending aorta in male Ipo8-/- mice throughout life. Immunohistological investigation of mutant aortic walls reveals elastic fiber disorganization and fragmentation along with a signature of increased TGF-ß signaling, as evidenced by nuclear pSmad2 accumulation. RT-qPCR assays of the aortic wall in male Ipo8-/- mice demonstrate decreased Smad6/7 and increased Mmp2 and Ccn2 (Ctgf) expression, reinforcing a role for dysregulation of the TGF-ß signaling pathway in TAA development. Because importin 8 is the most downstream TGF-ß-related effector implicated in TAA pathogenesis so far, it offers opportunities for future mechanistic studies and represents a candidate drug target for TAA.

SCUBE3 loss-of-function causes a recognizable recessive developmental disorder due to defective bone morphogenetic protein signaling.

Signal peptide-CUB-EGF domain-containing protein 3 (SCUBE3) is a member of a small family of multifunctional cell surface-anchored glycoproteins functioning as co-receptors for a variety of growth factors. Here we report that bi-allelic inactivating variants in SCUBE3 have pleiotropic consequences on development and cause a previously unrecognized syndromic disorder. Eighteen affected individuals from nine unrelated families showed a consistent phenotype characterized by reduced growth, skeletal features, distinctive craniofacial appearance, and dental anomalies. In vitro functional validation studies demonstrated a variable impact of disease-causing variants on transcript processing, protein secretion and function, and their dysregulating effect on bone morphogenetic protein (BMP) signaling. We show that SCUBE3 acts as a BMP2/BMP4 co-receptor, recruits the BMP receptor complexes into raft microdomains, and positively modulates signaling possibly by augmenting the specific interactions between BMPs and BMP type I receptors. Scube3-/- mice showed craniofacial and dental defects, reduced body size, and defective endochondral bone growth due to impaired BMP-mediated chondrogenesis and osteogenesis, recapitulating the human disorder. Our findings identify a human disease caused by defective function of a member of the SCUBE family, and link SCUBE3 to processes controlling growth, morphogenesis, and bone and teeth development through modulation of BMP signaling.

Conclusion of diagnostic odysseys due to inversions disrupting GLI3 and FBN1.

Many genetic testing methodologies are biased towards picking up structural variants (SVs) that alter copy number. Copy-neutral rearrangements such as inversions are therefore likely to suffer from underascertainment. In this study, manual review prompted by a virtual multidisciplinary team meeting and subsequent bioinformatic prioritisation of data from the 100K Genomes Project was performed across 43 genes linked to well-characterised skeletal disorders. Ten individuals from three independent families were found to harbour diagnostic inversions. In two families, inverted segments of 1.2/14.8 Mb unequivocally disrupted GLI3 and segregated with skeletal features consistent with Greig cephalopolysyndactyly syndrome. For one family, phenotypic blending was due to the opposing breakpoint lying ~45 kb from HOXA13 In the third family, long suspected to have Marfan syndrome, a 2.0 Mb inversion disrupting FBN1 was identified. These findings resolved lengthy diagnostic odysseys of 9-20 years and highlight the importance of direct interaction between clinicians and data-analysts. These exemplars of a rare mutational class inform future SV prioritisation strategies within the NHS Genomic Medicine Service and similar genome sequencing initiatives. In over 30 years since these two disease-gene associations were identified, large inversions have yet to be described and so our results extend the mutational spectra linked to these conditions.

Use of genome sequencing to hunt for cryptic second-hit variants: analysis of 31 cases recruited to the 100 000 Genomes Project.

BACKGROUND: Current clinical testing methods used to uncover the genetic basis of rare disease have inherent limitations, which can lead to causative pathogenic variants being missed. Within the rare disease arm of the 100 000 Genomes Project (100kGP), families were recruited under the clinical indication 'single autosomal recessive mutation in rare disease'. These participants presented with strong clinical suspicion for a specific autosomal recessive disorder, but only one suspected pathogenic variant had been identified through standard-of-care testing. Whole genome sequencing (WGS) aimed to identify cryptic 'second-hit' variants. METHODS: To investigate the 31 families with available data that remained unsolved following formal review within the 100kGP, SVRare was used to aggregate structural variants present in <1% of 100kGP participants. Small variants were assessed using population allele frequency data and SpliceAI. Literature searches and publicly available online tools were used for further annotation of pathogenicity. RESULTS: Using these strategies, 8/31 cases were solved, increasing the overall diagnostic yield of this cohort from 10/41 (24.4%) to 18/41 (43.9%). Exemplar cases include a patient with cystic fibrosis harbouring a novel exonic LINE1 insertion in CFTR and a patient with generalised arterial calcification of infancy with complex interlinked duplications involving exons 2-6 of ENPP1. Although ambiguous by short-read WGS, the ENPP1 variant structure was resolved using optical genome mapping and RNA analysis. CONCLUSION: Systematic examination of cryptic variants across a multi-disease cohort successfully identifies additional pathogenic variants. WGS data analysis in autosomal recessive rare disease should consider complex structural and small intronic variants as potentially pathogenic second hits.

GREEN-DB: a framework for the annotation and prioritization of non-coding regulatory variants from whole-genome sequencing data.

Non-coding variants have long been recognized as important contributors to common disease risks, but with the expansion of clinical whole genome sequencing, examples of rare, high-impact non-coding variants are also accumulating. Despite recent advances in the study of regulatory elements and the availability of specialized data collections, the systematic annotation of non-coding variants from genome sequencing remains challenging. Here, we propose a new framework for the prioritization of non-coding regulatory variants that integrates information about regulatory regions with prediction scores and HPO-based prioritization. Firstly, we created a comprehensive collection of annotations for regulatory regions including a database of 2.4 million regulatory elements (GREEN-DB) annotated with controlled gene(s), tissue(s) and associated phenotype(s) where available. Secondly, we calculated a variation constraint metric and showed that constrained regulatory regions associate with disease-associated genes and essential genes from mouse knock-outs. Thirdly, we compared 19 non-coding impact prediction scores providing suggestions for variant prioritization. Finally, we developed a VCF annotation tool (GREEN-VARAN) that can integrate all these elements to annotate variants for their potential regulatory impact. In our evaluation, we show that GREEN-DB can capture previously published disease-associated non-coding variants as well as identify additional candidate disease genes in trio analyses.

Biallelic PRMT7 pathogenic variants are associated with a recognizable syndromic neurodevelopmental disorder with short stature, obesity, and craniofacial and digital abnormalities.

PURPOSE: Protein arginine methyltransferase 7 (PRMT7) is a member of a family of enzymes that catalyzes the methylation of arginine residues on several protein substrates. Biallelic pathogenic PRMT7 variants have previously been associated with a syndromic neurodevelopmental disorder characterized by short stature, brachydactyly, intellectual developmental disability, and seizures. To our knowledge, no comprehensive study describes the detailed clinical characteristics of this syndrome. Thus, we aim to delineate the phenotypic spectrum of PRMT7-related disorder. METHODS: We assembled a cohort of 51 affected individuals from 39 different families, gathering clinical information from 36 newly described affected individuals and reviewing data of 15 individuals from the literature. RESULTS: The main clinical characteristics of the PRMT7-related syndrome are short stature, mild to severe developmental delay/intellectual disability, hypotonia, brachydactyly, and distinct facial morphology, including bifrontal narrowing, prominent supraorbital ridges, sparse eyebrows, short nose with full/broad nasal tip, thin upper lip, full and everted lower lip, and a prominent or squared-off jaw. Additional variable findings include seizures, obesity, nonspecific magnetic resonance imaging abnormalities, eye abnormalities (i.e., strabismus or nystagmus), and hearing loss. CONCLUSION: This study further delineates and expands the molecular, phenotypic spectrum and natural history of PRMT7-related syndrome characterized by a neurodevelopmental disorder with skeletal, growth, and endocrine abnormalities.

Biallelic variants in PIGN cause Fryns syndrome, multiple congenital anomalies-hypotonia-seizures syndrome, and neurologic phenotypes: A genotype-phenotype correlation study.

PURPOSE: Biallelic PIGN variants have been described in Fryns syndrome, multiple congenital anomalies-hypotonia-seizure syndrome (MCAHS), and neurologic phenotypes. The full spectrum of clinical manifestations in relation to the genotypes is yet to be reported. METHODS: Genotype and phenotype data were collated and analyzed for 61 biallelic PIGN cases: 21 new and 40 previously published cases. Functional analysis was performed for 2 recurrent variants (c.2679C>G p.Ser893Arg and c.932T>G p.Leu311Trp). RESULTS: Biallelic-truncating variants were detected in 16 patients-10 with Fryns syndrome, 1 with MCAHS1, 2 with Fryns syndrome/MCAHS1, and 3 with neurologic phenotype. There was an increased risk of prenatal or neonatal death within this group (6 deaths were in utero or within 2 months of life; 6 pregnancies were terminated). Incidence of polyhydramnios, congenital anomalies (eg, diaphragmatic hernia), and dysmorphism was significantly increased. Biallelic missense or mixed genotype were reported in the remaining 45 cases-32 showed a neurologic phenotype and 12 had MCAHS1. No cases of diaphragmatic hernia or abdominal wall defects were seen in this group except patient 1 in which we found the missense variant p.Ser893Arg to result in functionally null alleles, suggesting the possibility of an undescribed functionally important region in the final exon. For all genotypes, there was complete penetrance for developmental delay and near-complete penetrance for seizures and hypotonia in patients surviving the neonatal period. CONCLUSION: We have expanded the described spectrum of phenotypes and natural history associated with biallelic PIGN variants. Our study shows that biallelic-truncating variants usually result in the more severe Fryns syndrome phenotype, but neurologic problems, such as developmental delay, seizures, and hypotonia, present across all genotypes. Functional analysis should be considered when the genotypes do not correlate with the predicted phenotype because there may be other functionally important regions in PIGN that are yet to be discovered.

Genome sequencing identifies KMT2E-disrupting cryptic structural variant in a female with O'Donnell-Luria-Rodan syndrome.

We describe a patient from the 100,000 Genomes Project with a complex de novo structural variant within KMT2E leading to O'Donnell-Luria-Rodan syndrome. This case expands the mutational spectrum for this syndrome and highlights the importance of revisiting unsolved cases using better SV prioritisation tools and updated gene panels.

The prevalence and phenotypic range associated with biallelic PKDCC variants.

PKDCC encodes a component of Hedgehog signalling required for normal chondrogenesis and skeletal development. Although biallelic PKDCC variants have been implicated in rhizomelic shortening of limbs with variable dysmorphic features, this association was based on just two patients. In this study, data from the 100 000 Genomes Project was used in conjunction with exome sequencing and panel-testing results accessed via international collaboration to assemble a cohort of eight individuals from seven independent families with biallelic PKDCC variants. The allelic series included six frameshifts, a previously described splice-donor site variant and a likely pathogenic missense variant observed in two families that was supported by in silico structural modelling. Database queries suggested that the prevalence of this condition is between 1 of 127 and 1 of 721 in clinical cohorts with skeletal dysplasia of unknown aetiology. Clinical assessments, combined with data from previously published cases, indicate a predominantly upper limb involvement. Micrognathia, hypertelorism and hearing loss appear to be commonly co-occurring features. In conclusion, this study strengthens the link between biallelic inactivation of PKDCC and rhizomelic limb-shortening and will enable clinical testing laboratories to better interpret variants in this gene.

Genomic and transcriptomic correlates of Richter transformation in chronic lymphocytic leukemia.

The transformation of chronic lymphocytic leukemia (CLL) to high-grade B-cell lymphoma is known as Richter syndrome (RS), a rare event with dismal prognosis. In this study, we conducted whole-genome sequencing (WGS) of paired circulating CLL (PB-CLL) and RS biopsies (tissue-RS) from 17 patients recruited into a clinical trial (CHOP-O). We found that tissue-RS was enriched for mutations in poor-risk CLL drivers and genes in the DNA damage response (DDR) pathway. In addition, we identified genomic aberrations not previously implicated in RS, including the protein tyrosine phosphatase receptor (PTPRD) and tumor necrosis factor receptor-associated factor 3 (TRAF3). In the noncoding genome, we discovered activation-induced cytidine deaminase-related and unrelated kataegis in tissue-RS affecting regulatory regions of key immune-regulatory genes. These include BTG2, CXCR4, NFATC1, PAX5, NOTCH-1, SLC44A5, FCRL3, SELL, TNIP2, and TRIM13. Furthermore, differences between the global mutation signatures of pairs of PB-CLL and tissue-RS samples implicate DDR as the dominant mechanism driving transformation. Pathway-based clonal deconvolution analysis showed that genes in the MAPK and DDR pathways demonstrate high clonal-expansion probability. Direct comparison of nodal-CLL and tissue-RS pairs from an independent cohort confirmed differential expression of the same pathways by RNA expression profiling. Our integrated analysis of WGS and RNA expression data significantly extends previous targeted approaches, which were limited by the lack of germline samples, and it facilitates the identification of novel genomic correlates implicated in RS transformation, which could be targeted therapeutically. Our results inform the future selection of investigative agents for a UK clinical platform study. This trial was registered at www.clinicaltrials.gov as #NCT03899337.

Biallelic variants in SLC38A3 encoding a glutamine transporter cause epileptic encephalopathy.

The solute carrier (SLC) superfamily encompasses >400 transmembrane transporters involved in the exchange of amino acids, nutrients, ions, metals, neurotransmitters and metabolites across biological membranes. SLCs are highly expressed in the mammalian brain; defects in nearly 100 unique SLC-encoding genes (OMIM: https://www.omim.org) are associated with rare Mendelian disorders including developmental and epileptic encephalopathy and severe neurodevelopmental disorders. Exome sequencing and family-based rare variant analyses on a cohort with neurodevelopmental disorders identified two siblings with developmental and epileptic encephalopathy and a shared deleterious homozygous splicing variant in SLC38A3. The gene encodes SNAT3, a sodium-coupled neutral amino acid transporter and a principal transporter of the amino acids asparagine, histidine, and glutamine, the latter being the precursor for the neurotransmitters GABA and glutamate. Additional subjects with a similar developmental and epileptic encephalopathy phenotype and biallelic predicted-damaging SLC38A3 variants were ascertained through GeneMatcher and collaborations with research and clinical molecular diagnostic laboratories. Untargeted metabolomic analysis was performed to identify novel metabolic biomarkers. Ten individuals from seven unrelated families from six different countries with deleterious biallelic variants in SLC38A3 were identified. Global developmental delay, intellectual disability, hypotonia, and absent speech were common features while microcephaly, epilepsy, and visual impairment were present in the majority. Epilepsy was drug-resistant in half. Metabolomic analysis revealed perturbations of glutamate, histidine, and nitrogen metabolism in plasma, urine, and CSF of selected subjects, potentially representing biomarkers of disease. Our data support the contention that SLC38A3 is a novel disease gene for developmental and epileptic encephalopathy and illuminate the likely pathophysiology of the disease as perturbations in glutamine homeostasis.

Recurrent de novo missense variants in GNB2 can cause syndromic intellectual disability.

PURPOSE: Binding proteins (G-proteins) mediate signalling pathways involved in diverse cellular functions and comprise Gα and Gßγ units. Human diseases have been reported for all five Gß proteins. A de novo missense variant in GNB2 was recently reported in one individual with developmental delay/intellectual disability (DD/ID) and dysmorphism. We aim to confirm GNB2 as a neurodevelopmental disease gene, and elucidate the GNB2-associated neurodevelopmental phenotype in a patient cohort. METHODS: We discovered a GNB2 variant in the index case via exome sequencing and sought individuals with GNB2 variants via international data-sharing initiatives. In silico modelling of the variants was assessed, along with multiple lines of evidence in keeping with American College of Medical Genetics and Genomics guidelines for interpretation of sequence variants. RESULTS: We identified 12 unrelated individuals with five de novo missense variants in GNB2, four of which are recurrent: p.(Ala73Thr), p.(Gly77Arg), p.(Lys89Glu) and p.(Lys89Thr). All individuals have DD/ID with variable dysmorphism and extraneurologic features. The variants are located at the universally conserved shared interface with the Gα subunit, which modelling suggests weaken this interaction. CONCLUSION: Missense variants in GNB2 cause a congenital neurodevelopmental disorder with variable syndromic features, broadening the spectrum of multisystem phenotypes associated with variants in genes encoding G-proteins.

De Novo Heterozygous POLR2A Variants Cause a Neurodevelopmental Syndrome with Profound Infantile-Onset Hypotonia.

The RNA polymerase II complex (pol II) is responsible for transcription of all ∼21,000 human protein-encoding genes. Here, we describe sixteen individuals harboring de novo heterozygous variants in POLR2A, encoding RPB1, the largest subunit of pol II. An iterative approach combining structural evaluation and mass spectrometry analyses, the use of S. cerevisiae as a model system, and the assessment of cell viability in HeLa cells allowed us to classify eleven variants as probably disease-causing and four variants as possibly disease-causing. The significance of one variant remains unresolved. By quantification of phenotypic severity, we could distinguish mild and severe phenotypic consequences of the disease-causing variants. Missense variants expected to exert only mild structural effects led to a malfunctioning pol II enzyme, thereby inducing a dominant-negative effect on gene transcription. Intriguingly, individuals carrying these variants presented with a severe phenotype dominated by profound infantile-onset hypotonia and developmental delay. Conversely, individuals carrying variants expected to result in complete loss of function, thus reduced levels of functional pol II from the normal allele, exhibited the mildest phenotypes. We conclude that subtle variants that are central in functionally important domains of POLR2A cause a neurodevelopmental syndrome characterized by profound infantile-onset hypotonia and developmental delay through a dominant-negative effect on pol-II-mediated transcription of DNA.

Genome sequencing reveals underdiagnosis of primary ciliary dyskinesia in bronchiectasis.

BACKGROUND: Bronchiectasis can result from infectious, genetic, immunological and allergic causes. 60-80% of cases are idiopathic, but a well-recognised genetic cause is the motile ciliopathy, primary ciliary dyskinesia (PCD). Diagnosis of PCD has management implications including addressing comorbidities, implementing genetic and fertility counselling and future access to PCD-specific treatments. Diagnostic testing can be complex; however, PCD genetic testing is moving rapidly from research into clinical diagnostics and would confirm the cause of bronchiectasis. METHODS: This observational study used genetic data from severe bronchiectasis patients recruited to the UK 100,000 Genomes Project and patients referred for gene panel testing within a tertiary respiratory hospital. Patients referred for genetic testing due to clinical suspicion of PCD were excluded from both analyses. Data were accessed from the British Thoracic Society audit, to investigate whether motile ciliopathies are underdiagnosed in people with bronchiectasis in the UK. RESULTS: Pathogenic or likely pathogenic variants were identified in motile ciliopathy genes in 17 (12%) out of 142 individuals by whole-genome sequencing. Similarly, in a single centre with access to pathological diagnostic facilities, 5-10% of patients received a PCD diagnosis by gene panel, often linked to normal/inconclusive nasal nitric oxide and cilia functional test results. In 4898 audited patients with bronchiectasis, <2% were tested for PCD and <1% received genetic testing. CONCLUSIONS: PCD is underdiagnosed as a cause of bronchiectasis. Increased uptake of genetic testing may help to identify bronchiectasis due to motile ciliopathies and ensure appropriate management.

A statistical approach for tracking clonal dynamics in cancer using longitudinal next-generation sequencing data.

MOTIVATION: Tumours are composed of distinct cancer cell populations (clones), which continuously adapt to their local micro-environment. Standard methods for clonal deconvolution seek to identify groups of mutations and estimate the prevalence of each group in the tumour, while considering its purity and copy number profile. These methods have been applied on cross-sectional data and on longitudinal data after discarding information on the timing of sample collection. Two key questions are how can we incorporate such information in our analyses and is there any benefit in doing so? RESULTS: We developed a clonal deconvolution method, which incorporates explicitly the temporal spacing of longitudinally sampled tumours. By merging a Dirichlet Process Mixture Model with Gaussian Process priors and using as input a sequence of several sparsely collected samples, our method can reconstruct the temporal profile of the abundance of any mutation cluster supported by the data as a continuous function of time. We benchmarked our method on whole genome, whole exome and targeted sequencing data from patients with chronic lymphocytic leukaemia, on liquid biopsy data from a patient with melanoma and on synthetic data and we found that incorporating information on the timing of tissue collection improves model performance, as long as data of sufficient volume and complexity are available for estimating free model parameters. Thus, our approach is particularly useful when collecting a relatively long sequence of tumour samples is feasible, as in liquid cancers (e.g. leukaemia) and liquid biopsies. AVAILABILITY AND IMPLEMENTATION: The statistical methodology presented in this paper is freely available at github.com/dvav/clonosGP. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Can polygenic risk scores contribute to cost-effective cancer screening? A systematic review.

PURPOSE: Polygenic risk influences susceptibility to cancer. We assessed whether polygenic risk scores could be used in conjunction with other predictors of future disease status in cost-effective risk-stratified screening for cancer. METHODS: We undertook a systematic review of papers that evaluated the cost-effectiveness of screening interventions informed by polygenic risk scores compared with more conventional screening modalities. We included papers reporting cost-effectiveness outcomes with no restriction on type of cancer or form of polygenic risk modeled. We evaluated studies using the Quality of Health Economic Studies checklist. RESULTS: A total of 10 studies were included in the review, which investigated 3 cancers: prostate (n = 5), colorectal (n = 3), and breast (n = 2). Of the 10 papers, 9 scored highly (score >75 on a 0-100 scale) when assessed using the quality checklist. Of the 10 studies, 8 concluded that polygenic risk-informed cancer screening was likely to be more cost-effective than alternatives. CONCLUSION: Despite the positive conclusions of the included studies, it is unclear if polygenic risk stratification will contribute to cost-effective cancer screening given the absence of robust evidence on the costs of polygenic risk stratification, the effects of differential ancestry, potential downstream economic sequalae, and how large volumes of polygenic risk data would be collected and used.