Quantifying single nucleotide variant detection sensitivity in exome sequencing

BACKGROUND: The targeted capture and sequencing of genomic regions has rapidly demonstrated its utility in genetic studies. Inherent in this technology is considerable heterogeneity of target coverage and this is expected to systematically impact our sensitivity to detect genuine polymorphisms. To fully interpret the polymorphisms identified in a genetic study it is often essential to both detect polymorphisms and to understand where and with what probability real polymorphisms may have been missed. RESULTS: Using down-sampling of 30 deeply sequenced exomes and a set of gold-standard single nucleotide variant (SNV) genotype calls for each sample, we developed an empirical model relating the read depth at a polymorphic site to the probability of calling the correct genotype at that site. We find that measured sensitivity in SNV detection is substantially worse than that predicted from the naive expectation of sampling from a binomial. This calibrated model allows us to produce single nucleotide resolution SNV sensitivity estimates which can be merged to give summary sensitivity measures for any arbitrary partition of the target sequences (nucleotide, exon, gene, pathway, exome). These metrics are directly comparable between platforms and can be combined between samples to give “power estimates” for an entire study. We estimate a local read depth of 13X is required to detect the alleles and genotype of a heterozygous SNV 95% of the time, but only 3X for a homozygous SNV. At a mean on-target read depth of 20X, commonly used for rare disease exome sequencing studies, we predict 5–15% of heterozygous and 1–4% of homozygous SNVs in the targeted regions will be missed. CONCLUSIONS: Non-reference alleles in the heterozygote state have a high chance of being missed when commonly applied read coverage thresholds are used despite the widely held assumption that there is good polymorphism detection at these coverage levels. Such alleles are likely to be of functional importance in population based studies of rare diseases, somatic mutations in cancer and explaining the “missing heritability” of quantitative traits

Using high angular resolution diffusion imaging data to discriminate cortical regions

Brodmann's 100-year-old summary map has been widely used for cortical localization in neuroscience. There is a pressing need to update this map using non-invasive, high-resolution and reproducible data, in a way that captures individual variability. We demonstrate here that standard HARDI data has sufficiently diverse directional variation among grey matter regions to inform parcellation into distinct functional regions, and that this variation is reproducible across scans. This characterization of the signal variation as non-random and reproducible is the critical condition for successful cortical parcellation using HARDI data. This paper is a first step towards an individual cortex-wide map of grey matter microstructure, The gray/white matter and pial boundaries were identified on the high-resolution structural MRI images. Two HARDI data sets were collected from each individual and aligned with the corresponding structural image. At each vertex point on the surface tessellation, the diffusion-weighted signal was extracted from each image in the HARDI data set at a point, half way between gray/white matter and pial boundaries. We then derived several features of the HARDI profile with respect to the local cortical normal direction, as well as several fully orientationally invariant features. These features were taken as a fingerprint of the underlying grey matter tissue, and used to distinguish separate cortical areas. A support-vector machine classifier, trained on three distinct areas in repeat 1 achieved 80-82% correct classification of the same three areas in the unseen data from repeat 2 in three volunteers. Though gray matter anisotropy has been mostly overlooked hitherto, this approach may eventually form the foundation of a new cortical parcellation method in living humans. Our approach allows for further studies on the consistency of HARDI based parcellation across subjects and comparison with independent microstructural measures such as ex-vivo histology

The inference of gray whale (Eschrichtius robustus) historical population attributes from whole-genome sequences

Commercial whaling caused extensive demographic declines in many great whale species, including gray whales that were extirpated from the Atlantic Ocean and dramatically reduced in the Pacific Ocean. The Eastern Pacific gray whale has recovered since the 1982 ban on commercial whaling, but the Western Pacific gray whale-once considered possibly extinct-consists of only about 200 individuals and is considered critically endangered by some international authorities. Herein, we use whole-genome sequencing to investigate the demographic history of gray whales from the Pacific and use environmental niche modelling to make predictions about future gene flow.Our sequencing efforts and habitat niche modelling indicate that: i) western gray whale effective population sizes have declined since the last glacial maximum; ii) contemporary gray whale genomes, both eastern and western, harbor less autosomal nucleotide diversity than most other marine mammals and megafauna; iii) the extent of inbreeding, as measured by autozygosity, is greater in the Western Pacific than in the Eastern Pacific populations; and iv) future climate change is expected to open new migratory routes for gray whales.Our results indicate that gray whale genomes contain low nucleotide diversity and have been subject to both historical and recent inbreeding. Population sizes over the last million years likely peaked about 25,000 years before present and have declined since then. Our niche modelling suggests that novel migratory routes may develop within the next century and if so this could help retain overall genetic diversity, which is essential for adaption and successful recovery in light of global environmental change and past exploitation

The Medical Genome Reference Bank contains whole genome and phenotype data of 2570 healthy elderly

Population health research is increasingly focused on the genetic determinants of healthy ageing, but there is no public resource of whole genome sequences and phenotype data from healthy elderly individuals. Here we describe the first release of the Medical Genome Reference Bank (MGRB), comprising whole genome sequence and phenotype of 2570 elderly Australians depleted for cancer, cardiovascular disease, and dementia. We analyse the MGRB for single-nucleotide, indel and structural variation in the nuclear and mitochondrial genomes. MGRB individuals have fewer disease-associated common and rare germline variants, relative to both cancer cases and the gnomAD and UK Biobank cohorts, consistent with risk depletion. Age-related somatic changes are correlated with grip strength in men, suggesting blood-derived whole genomes may also provide a biologic measure of age-related functional deterioration. The MGRB provides a broadly applicable reference cohort for clinical genetics and genomic association studies, and for understanding the genetics of healthy ageing

Alternative pathway dysregulation in tissues drives sustained complement activation and predicts outcome across the disease course in COVID-19

Complement, a critical defence against pathogens, has been implicated as a driver of pathology in COVID-19. Complement activation products are detected in plasma and tissues and complement blockade considered for therapy. To delineate roles of complement in immunopathogenesis, we undertook the largest comprehensive study of complement in an COVID-19 to date, a comprehensive profiling of 16 complement biomarkers, including key components, regulators and activation products, in 966 plasma samples from 682 hospitalised COVID-19 patients collected across the hospitalisation period as part of the UK ISARIC4C study. Unsupervised clustering of complement biomarkers mapped to disease severity and supervised machine learning identified marker sets in early samples that predicted peak severity. Compared to heathy controls, complement proteins and activation products (Ba, iC3b, terminal complement complex) were significantly altered in COVID-19 admission samples in all severity groups. Elevated alternative pathway activation markers (Ba and iC3b) and decreased alternative pathway regulator (properdin) in admission samples associated with more severe disease and risk of death. Levels of most complement biomarkers were reduced in severe disease, consistent with consumption and tissue deposition. Latent class mixed modelling and cumulative incidence analysis identified the trajectory of increase of Ba to be a strong predictor of peak COVID-19 disease severity and death. The data demonstrate that early-onset, uncontrolled activation of complement, driven by sustained and progressive amplification through the alternative pathway amplification loop is a ubiquitous feature of COVID-19, further exacerbated in severe disease. These findings provide novel insights into COVID-19 immunopathogenesis and inform strategies for therapeutic intervention

Para-infectious brain injury in COVID-19 persists at follow-up despite attenuated cytokine and autoantibody responses

Data Availability Statement: The individual-level data from these studies is not publicly available to main confidentiality. Data generated by the ISARIC4C consortium is available for collaborative analysis projects through an independent data and materials access committee at isaric4c.net/sample_access. Data and samples from the COVID-Clinical Neuroscience Study are available through collaborative research by application through the NIHR bioresource at https://bioresource.nihr.ac.uk/using-our-bioresource/apply-for-bioresource-data-access/. Brain injury marker and immune mediator data are present in the paper and in the source data file. Source data are provided with this paper.To understand neurological complications of COVID-19 better both acutely and for recovery, we measured markers of brain injury, inflammatory mediators, and autoantibodies in 203 hospitalised participants; 111 with acute sera (1–11 days post-admission) and 92 convalescent sera (56 with COVID-19-associated neurological diagnoses). Here we show that compared to 60 uninfected controls, tTau, GFAP, NfL, and UCH-L1 are increased with COVID-19 infection at acute timepoints and NfL and GFAP are significantly higher in participants with neurological complications. Inflammatory mediators (IL-6, IL-12p40, HGF, M-CSF, CCL2, and IL-1RA) are associated with both altered consciousness and markers of brain injury. Autoantibodies are more common in COVID-19 than controls and some (including against MYL7, UCH-L1, and GRIN3B) are more frequent with altered consciousness. Additionally, convalescent participants with neurological complications show elevated GFAP and NfL, unrelated to attenuated systemic inflammatory mediators and to autoantibody responses. Overall, neurological complications of COVID-19 are associated with evidence of neuroglial injury in both acute and late disease and these correlate with dysregulated innate and adaptive immune responses acutely.National Institute for Health and Care Research (NIHR) (CO-CIN-01) and jointly by NIHR and UK Research and Innovation (CV220-169, MC_PC_19059). B.D.M. is supported by the UKRI/MRC (MR/V03605X/1), the MRC/UKRI (MR/V007181/1), MRC (MR/T028750/1) and Wellcome (ISSF201902/3). C.D. is supported by MRC (MC_PC_19044). We would like to thank the University of Liverpool GCP laboratory facility team for Luminex assistance and the Liverpool University Biobank team for all their help, especially Dr. Victoria Shaw, Lara Lavelle-Langham, and Sue Holden. We would like to acknowledge the Liverpool Experimental Cancer Medicine Centre for providing infrastructure support for this research (Grant Reference: C18616/A25153). We acknowledge the Liverpool Centre for Cell Imaging (CCI) for provision of imaging equipment (Dragonfly confocal microscope) and excellent technical assistance (BBSRC grant number BB/R01390X/1). Tom Solomon is supported by The Pandemic Institute and the NIHR Health Protection Research Unit (HPRU) in Emerging and Zoonotic Infections at University of Liverpool. D.K.M. and E.N. are supported by the NIHR Cambridge Biomedical Centre and by NIHR funding to the NIHR BioResource (RG94028 and RG85445), and by funding from Brain Research UK 201819-20. We thank NIHR BioResource volunteers for their participation, and gratefully acknowledge NIHR BioResource centres, NHS Trusts and staff for their contribution. We thank the National Institute for Health and Care Research, NHS Blood and Transplant, and Health Data Research UK as part of the Digital Innovation Hub Programme. Support for title page creation and format was provided by AuthorArranger, a tool developed at the National Cancer Institute. The authors would like to acknowledge the eDRIS team (Public Health Scotland) for their support in obtaining approvals, the provisioning and linking of data and facilitating access to the National Safe Haven. The views expressed are those of the author(s) and not necessarily those of the UKRI, NHS, the NIHR or the Department of Health and Social Care

Implementation of corticosteroids in treatment of COVID-19 in the ISARIC WHO Clinical Characterisation Protocol UK: prospective, cohort study

Background: Dexamethasone was the first intervention proven to reduce mortality in patients with COVID-19 being treated in hospital. We aimed to evaluate the adoption of corticosteroids in the treatment of COVID-19 in the UK after the RECOVERY trial publication on June 16, 2020, and to identify discrepancies in care. Methods: We did an audit of clinical implementation of corticosteroids in a prospective, observational, cohort study in 237 UK acute care hospitals between March 16, 2020, and April 14, 2021, restricted to patients aged 18 years or older with proven or high likelihood of COVID-19, who received supplementary oxygen. The primary outcome was administration of dexamethasone, prednisolone, hydrocortisone, or methylprednisolone. This study is registered with ISRCTN, ISRCTN66726260. Findings: Between June 17, 2020, and April 14, 2021, 47 795 (75·2%) of 63 525 of patients on supplementary oxygen received corticosteroids, higher among patients requiring critical care than in those who received ward care (11 185 [86·6%] of 12 909 vs 36 415 [72·4%] of 50 278). Patients 50 years or older were significantly less likely to receive corticosteroids than those younger than 50 years (adjusted odds ratio 0·79 [95% CI 0·70–0·89], p=0·0001, for 70–79 years; 0·52 [0·46–0·58], p80 years), independent of patient demographics and illness severity. 84 (54·2%) of 155 pregnant women received corticosteroids. Rates of corticosteroid administration increased from 27·5% in the week before June 16, 2020, to 75–80% in January, 2021. Interpretation: Implementation of corticosteroids into clinical practice in the UK for patients with COVID-19 has been successful, but not universal. Patients older than 70 years, independent of illness severity, chronic neurological disease, and dementia, were less likely to receive corticosteroids than those who were younger, as were pregnant women. This could reflect appropriate clinical decision making, but the possibility of inequitable access to life-saving care should be considered. Funding: UK National Institute for Health Research and UK Medical Research Council