Improved imputation of low-frequency and rare variants using the UK10K haplotype reference panel

Imputing genotypes from reference panels created by whole-genome sequencing (WGS) provides a cost-effective strategy for augmenting the single-nucleotide polymorphism (SNP) content of genome-wide arrays. The UK10K Cohorts project has generated a data set of 3,781 whole genomes sequenced at low depth (average 7x), aiming to exhaustively characterize genetic variation down to 0.1% minor allele frequency in the British population. Here we demonstrate the value of this resource for improving imputation accuracy at rare and low-frequency variants in both a UK and an Italian population. We show that large increases in imputation accuracy can be achieved by re-phasing WGS reference panels after initial genotype calling. We also present a method for combining WGS panels to improve variant coverage and downstream imputation accuracy, which we illustrate by integrating 7,562 WGS haplotypes from the UK10K project with 2,184 haplotypes from the 1000 Genomes Project. Finally, we introduce a novel approximation that maintains speed without sacrificing imputation accuracy for rare variants

Whole-genome sequence-based analysis of thyroid function

Tiina Paunio on työryhmän UK10K Consortium jäsen.Normal thyroid function is essential for health, but its genetic architecture remains poorly understood. Here, for the heritable thyroid traits thyrotropin (TSH) and free thyroxine (FT4), we analyse whole-genome sequence data from the UK10K project (N = 2,287). Using additional whole-genome sequence and deeply imputed data sets, we report meta-analysis results for common variants (MAF >= 1%) associated with TSH and FT4 (N = 16,335). For TSH, we identify a novel variant in SYN2 (MAF = 23.5%, P = 6.15 x 10(-9)) and a new independent variant in PDE8B (MAF = 10.4%, P = 5.94 x 10(-14)). For FT4, we report a low-frequency variant near B4GALT6/ SLC25A52 (MAF = 3.2%, P = 1.27 x 10(-9)) tagging a rare TTR variant (MAF = 0.4%, P = 2.14 x 10(-11)). All common variants explain >= 20% of the variance in TSH and FT4. Analysis of rare variants (MAFPeer reviewe

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia.

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Cav2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca2+ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Cav2.2 in normal human neurodevelopment.MAK is funded by an NIHR Research Professorship and receives funding from the Wellcome Trust, Great Ormond Street Children's Hospital Charity, and Rosetrees Trust. E.M. received funding from the Rosetrees Trust (CD-A53) and Great Ormond Street Hospital Children's Charity. K.G. received funding from Temple Street Foundation. A.M. is funded by Great Ormond Street Hospital, the National Institute for Health Research (NIHR), and Biomedical Research Centre. F.L.R. and D.G. are funded by Cambridge Biomedical Research Centre. K.C. and A.S.J. are funded by NIHR Bioresource for Rare Diseases. The DDD Study presents independent research commissioned by the Health Innovation Challenge Fund (grant number HICF-1009-003), a parallel funding partnership between the Wellcome Trust and the Department of Health, and the Wellcome Trust Sanger Institute (grant number WT098051). We acknowledge support from the UK Department of Health via the NIHR comprehensive Biomedical Research Centre award to Guy's and St. Thomas' National Health Service (NHS) Foundation Trust in partnership with King's College London. This research was also supported by the NIHR Great Ormond Street Hospital Biomedical Research Centre. J.H.C. is in receipt of an NIHR Senior Investigator Award. The research team acknowledges the support of the NIHR through the Comprehensive Clinical Research Network. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR, Department of Health, or Wellcome Trust. E.R.M. acknowledges support from NIHR Cambridge Biomedical Research Centre, an NIHR Senior Investigator Award, and the University of Cambridge has received salary support in respect of E.R.M. from the NHS in the East of England through the Clinical Academic Reserve. I.E.S. is supported by the National Health and Medical Research Council of Australia (Program Grant and Practitioner Fellowship)

Performance of PVP in variant prioritization in WGS data, separated by mode of inheritance of the disease.

<p>Performance of PVP in variant prioritization in WGS data, separated by mode of inheritance of the disease.</p

Overview of the performance of PVP, CADD, DANN, GWAVA and Exomiser in prioritizing causative variants in WGS data.

<p>We prioritize all variants in a VCF file resulting from WGS using the same models. Analysis is separated reflecting the performance of the various tools identifying exonic and non-exonic variants. For CADD, DANN, and GWAVA, we report only analysis results for which a prediction score is returned; consequently, total numbers are less than the total of 11,251 causative variants.</p

Performance of PVP in identifying causative variants in whole genome sequences using human phenotypes (PVP-Human), model organisms phenotypes (PVP-Model), and combined phenotypes (PVP), and comparison of PVP to CADD, DANN, GWAVA, and Genomiser.

<p>Performance of PVP in identifying causative variants in whole genome sequences using human phenotypes (PVP-Human), model organisms phenotypes (PVP-Model), and combined phenotypes (PVP), and comparison of PVP to CADD, DANN, GWAVA, and Genomiser.</p

Homozygous loss-of-function mutations in SLC26A7 cause goitrous congenital hypothyroidism

Abstract Defects in genes mediating thyroid hormone biosynthesis result in dyshormonogenic congenital hypothyroidism (CH). Here, we report homozygous truncating mutations in SLC26A7 in 6 unrelated families with goitrous CH and show that goitrous hypothyroidism also occurs in Slc26a7-null mice. In both species, the gene is expressed predominantly in the thyroid gland, and loss of function is associated with impaired availability of iodine for thyroid hormone synthesis, partially corrected in mice by iodine supplementation.SLC26A7 is a member of the same transporter family as SLC26A4 (pendrin), an anion exchanger with affinity for iodide and chloride (among others), whose gene mutations cause congenital deafness and dyshormonogenic goiter. However, in contrast to pendrin, SLC26A7 does not mediate cellular iodide efflux and hearing in affected individuals is normal. We delineate a hitherto unrecognized role for SLC26A7 in thyroid hormone biosynthesis, for which the mechanism remains unclear

The IGSF1 deficiency syndrome: Characteristics of male and female patients

Context: Ig superfamily member1 (IGSF1) deficiency was recently discovered as a novel X-linked cause of central hypothyroidism (CeH) and macro-orchidism. However, clinical and biochemical data regarding growth, puberty, and metabolic outcome, as well as features of female carriers, are scarce. Objective: Our objective was to investigate clinical and biochemical characteristics associated with IGSF1 deficiency in both sexes.Methods: All patients (n=42, 24 males) from 10 families examined in the university clinics of Leiden, Amsterdam, Cambridge, and Milan were included in this case series. Detailed clinical data were collected with an identical protocol, and biochemical measurements were performed in a central laboratory. Results: Male patients (age 0-87 years, 17 index cases and 7 from family studies) showed CeH (100%), hypoprolactinemia (n = 16, 67%), and transient partial GH deficiency (n = 3, 13%). Pubertal testosterone productionwasdelayed, aswerethe growth spurtandpubic hair development. However, testicular growth started at a normal age and attained macro-orchid size in all evaluable adults. Body mass index, percent fat, and waist circumference tended to be elevated. The metabolic syndrome was present in 4 of 5 patients over 55 years of age. Heterozygous female carriers (age 32-80 years) showed CeH in 6 of 18 cases (33%), hypoprolactinemia in 2 (11%), and GH deficiency in none. As in men, body mass index, percent fat, and waist circumference were relatively high, and the metabolic syndrome was present in 3 cases. Conclusion: In male patients, the X-linked IGSF1 deficiency syndrome is characterized by CeH, hypoprolactinemia, delayed puberty, macro-orchidism, and increased body weight. A subset of female carriers also exhibits CeH

Unpaved road verges as hotspots of fleshy-fruited shrub recruitment and establishment

Hypothetical low-quality habitats can hold an overlooked conservation value. Some frugivorous mam- mals such as the red fox (Vulpes vulpes) and the European rabbit (Oryctolagus cuniculus) disperse many viable seeds of fleshy-fruited shrubs along the verges of soft linear developments (SLD), such as trails and firebreaks. However, seed arrival does not guarantee plant recruitment, since several post-dispersal processes can alter seed rain. To examine whether SLD verges assist shrub recruitment and establish- ment, we compared the density and the structure of a community of Mediterranean shrubs between SLD verges and the adjacent scrubland. Both seedlings and adult fleshy-fruited shrubs dispersed by foxes and rabbits reached higher densities along SLD verges than in the scrubland, suggesting SLD verges can be suitable habitats for shrub recruit- ment and establishment. Bird-dispersed shrubs showed a similar pattern, whereas shrubs dispersed by ungulates and badgers (Meles meles) as well as rockroses (Cistaceae) showed similar densities in both hab- itats. Shrub species composition and diversity were similar between habitats. Due to a marked differential seed arrival, SLD verges housed higher densities of fleshy-fruited shrubs than the adjacent scrubland. Established shrubs may attract seed-dispersing wildlife, and create proper environments for plant recruitment, generating a reforestation feedback. Incipient shrub populations along roadsides may act as stepping stones with potential to connect isolated populations in fragmented landscapes, where SLD are pervasive. We recommend careful management of frugivore populations and SLD verges in order to favor the diversity and the structural complexity of native vegetation while preventing the spread of invasive species.Peer reviewe

Rare Variant Analysis of Human and Rodent Obesity Genes in Individuals with Severe Childhood Obesity

Obesity is a genetically heterogeneous disorder. Using targeted and whole-exome sequencing, we studied 32 human and 87 rodent obesity genes in 2,548 severely obese children and 1,117 controls. We identified 52 variants contributing to obesity in 2% of cases including multiple novel variants in GNAS, which were sometimes found with accelerated growth rather than short stature as describedw previously. Nominally significant associations were found for rare functional variants in BBS1, BBS9, GNAS, MKKS, CLOCK and ANGPTL6. The p.S284X variant in ANGPTL6 drives the association signal (rs201622589, MAF∼0.1%, odds ratio = 10.13, p-value = 0.042) and results in complete loss of secretion in cells. Further analysis including additional case-control studies and population controls (N = 260,642) did not support association of this variant with obesity (odds ratio = 2.34, p-value = 2.59 × 10-3), highlighting the challenges of testing rare variant associations and the need for very large sample sizes. Further validation in cohorts with severe obesity and engineering the variants in model organisms will be needed to explore whether human variants in ANGPTL6 and other genes that lead to obesity when deleted in mice, do contribute to obesity. Such studies may yield druggable targets for weight loss therapies