GeneTIER: prioritization of candidate disease genes using tissue-specific gene expression profiles.

Motivation In attempts to determine the genetic causes of human disease, researchers are often faced with a large number of candidate genes. Linkage studies can point to a genomic region containing hundreds of genes, while the high-throughput sequencing approach will often identify a great number of non-synonymous genetic variants. Since systematic experimental verification of each such candidate gene is not feasible, a method is needed to decide which genes are worth investigating further. Computational gene prioritization presents itself as a solution to this problem, systematically analyzing and sorting each gene from the most to least likely to be the disease-causing gene, in a fraction of the time it would take a researcher to perform such queries manually. Results Here we present GeneTIER (Gene TIssue Expression Ranker), a new web-based application for candidate gene prioritization. GeneTIER replaces knowledge-based inference traditionally used in candidate disease gene prioritization applications with experimental data from tissue-specific gene expression datasets and thus largely overcomes the bias towards the better characterized genes/diseases that commonly afflict other methods. We show that our approach is capable of accurate candidate gene prioritization and illustrate its strengths and weaknesses using case study examples. Availability and Implementation Freely available on the web at http://dna.leeds.ac.uk/GeneTIER/ Contact: [email protected]

A Chromosome 7 Pericentric Inversion Defined at Single-Nucleotide Resolution Using Diagnostic Whole Genome Sequencing in a Patient with Hand-Foot-Genital Syndrome.

Next generation sequencing methodologies are facilitating the rapid characterisation of novel structural variants at nucleotide resolution. These approaches are particularly applicable to variants initially identified using alternative molecular methods. We report a child born with bilateral postaxial syndactyly of the feet and bilateral fifth finger clinodactyly. This was presumed to be an autosomal recessive syndrome, due to the family history of consanguinity. Karyotype analysis revealed a homozygous pericentric inversion of chromosome 7 (46,XX,inv(7)(p15q21)x2) which was confirmed to be heterozygous in both unaffected parents. Since the resolution of the karyotype was insufficient to identify any putatively causative gene, we undertook medium-coverage whole genome sequencing using paired-end reads, in order to elucidate the molecular breakpoints. In a two-step analysis, we first narrowed down the region by identifying discordant read-pairs, and then determined the precise molecular breakpoint by analysing the mapping locations of "soft-clipped" breakpoint-spanning reads. PCR and Sanger sequencing confirmed the identified breakpoints, both of which were located in intergenic regions. Significantly, the 7p15 breakpoint was located 523 kb upstream of HOXA13, the locus for hand-foot-genital syndrome. By inference from studies of HOXA locus control in the mouse, we suggest that the inversion has delocalised a HOXA13 enhancer to produce the phenotype observed in our patient. This study demonstrates how modern genetic diagnostic approach can characterise structural variants at nucleotide resolution and provide potential insights into functional regulation

Robust diagnostic genetic testing using solution capture enrichment and a novel variant-filtering interface.

Targeted hybridization enrichment prior to next-generation sequencing is a widespread method for characterizing sequence variation in a research setting, and is being adopted by diagnostic laboratories. However, the number of variants identified can overwhelm clinical laboratories with strict time constraints, the final interpretation of likely pathogenicity being a particular bottleneck. To address this, we have developed an approach in which, after automatic variant calling on a standard unix pipeline, subsequent variant filtering is performed interactively, using AgileExomeFilter and AgilePindelFilter (http://dna.leeds.ac.uk/agile), tools designed for clinical scientists with standard desktop computers. To demonstrate the method's diagnostic efficacy, we tested 128 patients using (1) a targeted capture of 36 cancer-predisposing genes or (2) whole-exome capture for diagnosis of the genetically heterogeneous disorder primary ciliary dyskinesia (PCD). In the cancer cohort, complete concordance with previous diagnostic data was achieved across 793 variant genotypes. A high yield (42%) was also achieved for exome-based PCD diagnosis, underscoring the scalability of our method. Simple adjustments to the variant filtering parameters further allowed the identification of a homozygous truncating mutation in a presumptive new PCD gene, DNAH8. These tools should allow diagnostic laboratories to expand their testing portfolios flexibly, using a standard set of reagents and techniques

Targeted nanopore sequencing enables complete characterisation of structural deletions initially identified using exon-based short-read sequencing strategies

Background The widespread adoption of exome sequencing has greatly increased the rate of genetic diagnosis for inherited conditions. However, the detection and validation of large deletions remains challenging. While numerous bioinformatics approaches have been developed to detect deletions from whole - exome sequencing and targeted panels, further work is typically required to define the physical breakpoints or integration sites. Accurate characterisation requires either expensive follow - up whole - genome sequencing or the time - consuming, laborious process of PCR walking, both of which are challenging when dealing with the repeat sequences which frequently intersect deletion breakpoints. The aim of this study was to develop a cost-effective, long-range sequencing method to characterise deletions. Methods Genomic DNA was amplified with primers spanning the deletion using long-range PCR and the products purified. Sequencing was performed on MinION flongle flowcells. The resulting fast5 files were basecalled using Guppy, trimmed using Porechop and aligned using Minimap2. Filtering was performed using NanoFilt. Nanopore sequencing results were verified by Sanger sequencing. Results Four cases with deletions detected following comparative read-depth analysis of targeted short-read sequencing were analysed. Nanopore sequencing defined breakpoints at the molecular level in all cases including homozygous breakpoints in EYS, CNGA1 and CNGB1 and a heterozygous deletion in PRPF31. All breakpoints were verified by Sanger sequencing. Conclusions In this study, a quick, accurate and cost - effective method is described to characterise deletions identified from exome, and similar data, using nanopore sequencing

Inherited CD19 Deficiency Does Not Impair Plasma Cell Formation or Response to CXCL12

Background The human CD19 antigen is expressed throughout B cell ontogeny with the exception of neoplastic plasma cells and a subset of normal plasma cells. CD19 plays a role in propagating signals from the B cell receptor and other receptors such as CXCR4 in mature B cells. Studies of CD19-deficient patients have confirmed its function during the initial stages of B cell activation and the production of memory B cells; however, its role in the later stages of B cell differentiation is unclear. Objective Using B cells from a newly identified CD19-deficient individual, we investigated the role of CD19 in the generation and function of plasma cells using an in vitro differentiation model. Methods Flow cytometry and long-read nanopore sequencing using locus-specific long-range amplification products were used to screen a patient with suspected primary immunodeficiency. Purified B cells from the patient and healthy controls were activated with CD40L, IL-21, IL-2, and anti-Ig, then transferred to different cytokine conditions to induce plasma cell differentiation. Subsequently, the cells were stimulated with CXCL12 to induce signalling through CXCR4. Phosphorylation of key downstream proteins including ERK and AKT was assessed by Western blotting. RNA-seq was also performed on in vitro differentiating cells. Results Long-read nanopore sequencing identified the homozygous pathogenic mutation c.622del (p.Ser208Profs*19) which was corroborated by the lack of CD19 cell surface staining. CD19-deficient B cells that are predominantly naïve generate phenotypically normal plasma cells with expected patterns of differentiation-associated genes and normal levels of CXCR4. Differentiated CD19-deficient cells were capable of responding to CXCL12; however, plasma cells derived from naïve B cells, both CD19-deficient and sufficient, had relatively diminished signaling compared to those generated from total B cells. Additionally, CD19 ligation on normal plasma cells results in AKT phosphorylation. Conclusion CD19 is not required for generation of antibody-secreting cells or the responses of these populations to CXCL12, but may alter the response other ligands that require CD19 potentially affecting localization, proliferation, or survival. The observed hypogammaglobulinemia in CD19-deficient individuals is therefore likely attributable to the lack of memory B cells

Second asymptomatic carotid surgery trial (ACST-2): a randomised comparison of carotid artery stenting versus carotid endarterectomy

Background: Among asymptomatic patients with severe carotid artery stenosis but no recent stroke or transient cerebral ischaemia, either carotid artery stenting (CAS) or carotid endarterectomy (CEA) can restore patency and reduce long-term stroke risks. However, from recent national registry data, each option causes about 1% procedural risk of disabling stroke or death. Comparison of their long-term protective effects requires large-scale randomised evidence. Methods: ACST-2 is an international multicentre randomised trial of CAS versus CEA among asymptomatic patients with severe stenosis thought to require intervention, interpreted with all other relevant trials. Patients were eligible if they had severe unilateral or bilateral carotid artery stenosis and both doctor and patient agreed that a carotid procedure should be undertaken, but they were substantially uncertain which one to choose. Patients were randomly allocated to CAS or CEA and followed up at 1 month and then annually, for a mean 5 years. Procedural events were those within 30 days of the intervention. Intention-to-treat analyses are provided. Analyses including procedural hazards use tabular methods. Analyses and meta-analyses of non-procedural strokes use Kaplan-Meier and log-rank methods. The trial is registered with the ISRCTN registry, ISRCTN21144362. Findings: Between Jan 15, 2008, and Dec 31, 2020, 3625 patients in 130 centres were randomly allocated, 1811 to CAS and 1814 to CEA, with good compliance, good medical therapy and a mean 5 years of follow-up. Overall, 1% had disabling stroke or death procedurally (15 allocated to CAS and 18 to CEA) and 2% had non-disabling procedural stroke (48 allocated to CAS and 29 to CEA). Kaplan-Meier estimates of 5-year non-procedural stroke were 2·5% in each group for fatal or disabling stroke, and 5·3% with CAS versus 4·5% with CEA for any stroke (rate ratio [RR] 1·16, 95% CI 0·86–1·57; p=0·33). Combining RRs for any non-procedural stroke in all CAS versus CEA trials, the RR was similar in symptomatic and asymptomatic patients (overall RR 1·11, 95% CI 0·91–1·32; p=0·21). Interpretation: Serious complications are similarly uncommon after competent CAS and CEA, and the long-term effects of these two carotid artery procedures on fatal or disabling stroke are comparable. Funding: UK Medical Research Council and Health Technology Assessment Programme

An alternative to array-based diagnostics: a prospectively recruited cohort, comparing arrayCGH to next-generation sequencing to evaluate foetal structural abnormalities

Molecular diagnostic investigations, following the identification of foetal abnormalities, are routinely performed using array comparative genomic hybridisation (aCGH). Despite the utility of this technique, contemporary approaches for the detection of copy number variation are typically based on next-generation sequencing (NGS). We sought to compare an in-house NGS-based workflow (CNVseq) with aCGH, for invasively obtained foetal samples from pregnancies complicated by foetal structural abnormality. DNA from 40 foetuses was screened using both 8 × 60 K aCGH oligoarrays and low-coverage whole genome sequencing. Sequencer-compatible libraries were combined in a ten–sample multiplex and sequenced using an Illumina HiSeq2500. The mean resolution of CNVseq was 29 kb, compared to 60 kb for aCGH analyses. Four clinically significant, concordant, copy number imbalances were detected using both techniques, however, genomic breakpoints were more precisely defined by CNVseq. This data indicates CNVseq is a robust and sensitive alternative to aCGH, for the prenatal investigation of foetuses with structural abnormalities. • Impact statement • What is already known about this subject? Copy number variant analysis using next-generation sequencing has been successfully applied to investigations of tumour specimens and patients with developmental delays. The application of our approach, to a prospective prenatal diagnosis cohort, has not hitherto been assessed. • What do the results of this study add? Next-generation sequencing has a comparable turnaround time and assay sensitivity to copy number variant analysis performed using array CGH. We demonstrate that having established a next-generation sequencing facility, high-throughput CNVseq sample processing and analysis can be undertaken within the framework of a regional diagnostic service. • What are the implications of these findings for clinical practice and/or further research? Array CGH is a legacy technology which is likely to be superseded by low-coverage whole genome sequencing, for the detection of copy number variants, in the prenatal diagnosis of structural abnormalities

Haplotyping Using Long-Range PCR and Nanopore Sequencing of Phase Variants; Lessons Learned From the ABCA4 Locus

Short-read next-generation sequencing has revolutionized our ability to identify variants underlying inherited diseases, but it does not allow the phasing of variants to clarify their diagnostic interpretation. The advent of widespread, increasingly accurate long-read sequencing has opened up new applications not currently available through short-read next-generation sequencing. One such use is the ability to phase variants to clarify their diagnostic interpretation and to investigate the increasingly prevalent role of cis-acting variants in the pathogenesis of the inherited disease, so-called complex alleles. Complex alleles are becoming an increasingly prevalent part of the study of genes associated with inherited disease, for example, in ABCA4-related disease. We sought to establish a cost-effective method to phase contiguous segments of the 130kb ABCA4 locus by long-read sequencing of overlapping amplification products. Using the comprehensively characterized CEPH sample, NA12878, we verified the accuracy and robustness of our assay. However, in-field assessment of its utility using clinical test cases was hampered by the paucity and distribution of identified variants and by PCR chimerism, particularly where the number of PCR cycles were high. Despite this, we were able to construct robust phase blocks of up to 94.9kb, representing 73% of the ABCA4 locus. We conclude that, while haplotype analysis of variants located within discrete amplification products was robust and informative, the stitching-together of larger phase blocks using overlapping single-molecule reads remained practically challenging