Search

Summary DeepPheWAS is an R package for phenome-wide association studies that creates clinically curated composite phenotypes and integrates quantitative phenotypes from primary care data, longitudinal trajectories of quantitative measures, disease progression and drug response phenotypes. Tools are provided for efficient analysis of association with any genetic input, under any genetic model, with optional sex-stratified analysis, and for developing novel phenotypes. Availability and implementation The DeepPheWAS R package is freely available under GNU general public licence v3.0 from at https://github.com/Richard-Packer/DeepPheWAS. Supplementary information Supplementary data are available at Bioinformatics online.

Lung function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry GWAS meta-analysis of lung function to date, comprising 580,869 participants, 1020 independent association signals identified 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score (GRS) showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies (PheWAS) for selected associated variants, and trait and pathway-specific GRS to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Background The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma. Methods We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping. Results Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11). Conclusions At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P-8) interactions with sex on lung function, and 21 showed suggestive interactions (P-6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV1) (P=3.15x10-15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV1 more in males (untransformed FEV1 β=0.028 [SE 0.0022] litres) than females (β=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein (HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10-6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

Objective\ud \ud Pathogenic variants in SCN3A , encoding the voltage‐gated sodium channel subunit Nav1.3, cause severe childhood‐onset epilepsy and malformation of cortical development. Here, we define the spectrum of clinical, genetic, and neuroimaging features of SCN3A ‐related neurodevelopmental disorder.\ud Methods\ud \ud Patients were ascertained via an international collaborative network. We compared sodium channels containing wild‐type vs. variant Nav1.3 subunits co‐expressed with β1 and β2 subunits using whole‐cell voltage clamp electrophysiological recordings in a heterologous mammalian system (HEK‐293 T cells).\ud Results\ud \ud Of 22 patients with pathogenic SCN3A variants, most had treatment‐resistant epilepsy beginning in the first year of life (16/21, 76%; median onset, 2 weeks), with severe or profound developmental delay (15/20; 75%). Many, but not all (15/19; 79%), exhibited malformations of cortical development. Pathogenic variants clustered in transmembrane segments 4–6 of domains II‐IV. Most pathogenic missense variants tested (10/11; 91%) displayed gain of channel function, with increased persistent current and/or a leftward shift in the voltage dependence of activation, and all variants associated with malformation of cortical development exhibited gain of channel function. One variant (p.Ile1468Arg) exhibited mixed effects, with gain and partial loss of function. Two variants demonstrated loss of channel function.\ud Interpretation\ud \ud Our study defines SCN3A‐ related neurodevelopmental disorder along a spectrum of severity, but typically including epilepsy and severe or profound developmental delay/intellectual disability. Malformations of cortical development are a characteristic feature of this unusual channelopathy syndrome, present in over 75% of affected individuals. Gain of function at the channel level in developing neurons is likely an important mechanism of disease pathogenesis.

Homozygosity for theSERPINA1Z allele causes alpha-1 antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on non-respiratory phenotypes, and on lung function in the general population, remain unclear. We conducted the largest population-based study to date to determine Z allele effects on >2,400 phenotypes using the UK Biobank study (N>303,353). We detected strong associations between heterozygosity and non-respiratory phenotypes including increased height, increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher lung function in non-smokers, but smoking appears to abolish this protective effect. Individuals heterozygous for the Z allele may therefore have altered risk of smoking-induced lung disease and other, non-respiratory conditions.

Reduced lung function predicts mortality and is key to the diagnosis of chronic obstructive pulmonary disease (COPD). In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, 139 of which are new. In combination, these variants strongly predict COPD in independent patient populations. Furthermore, the combined effect of these variants showed generalizability across smokers and never-smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD., Editorial summary: A genome-wide association study in over 400,000 individuals identifies 139 new signals for lung function. These variants can predict chronic obstructive pulmonary disease in independent, trans-ethnic cohorts.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P

PCGF2 encodes the polycomb group ring finger 2 protein, a transcriptional repressor involved in cell proliferation, differentiation, and embryogenesis. PCGF2 is a component of the polycomb repressive complex 1 (PRC1), a multiprotein complex which controls gene silencing through histone modification and chromatin remodelling. We report the phenotypic characterization of 13 patients (11 unrelated individuals and a pair of monozygotic twins) with missense mutations in PCGF2. All the mutations affected the same highly conserved proline in PCGF2 and were de novo, excepting maternal mosaicism in one. The patients demonstrated a recognizable facial gestalt, intellectual disability, feeding problems, impaired growth, and a range of brain, cardiovascular, and skeletal abnormalities. Computer structural modeling suggests the substitutions alter an N-terminal loop of PCGF2 critical for histone biding. Mutant PCGF2 may have dominant-negative effects, sequestering PRC1 components into complexes that lack the ability to interact efficiently with histones. These findings demonstrate the important role of PCGF2 in human development and confirm that heterozygous substitutions of the Pro65 residue of PCGF2 cause a recognizable syndrome characterized by distinctive craniofacial, neurological, cardiovascular, and skeletal features.

The redox state of the neural progenitors regulates physiological processes such as neuronal differentiation and dendritic and axonal growth. The relevance of endoplasmic reticulum (ER)-associated oxidoreductases in these processes is largely unexplored. We describe a severe neurological disorder caused by bi-allelic loss-of-function variants in thioredoxin (TRX)-related transmembrane-2 (TMX2); these variants were detected by exome sequencing in 14 affected individuals from ten unrelated families presenting with congenital microcephaly, cortical polymicrogyria, and other migration disorders. TMX2 encodes one of the five TMX proteins of the protein disulfide isomerase family, hitherto not linked to human developmental brain disease. Our mechanistic studies on protein function show that TMX2 localizes to the ER mitochondria-associated membranes (MAMs), is involved in posttranslational modification and protein folding, and undergoes physical interaction with the MAM-associated and ER folding chaperone calnexin and ER calcium pump SERCA2. These interactions are functionally relevant because TMX2-deficient fibroblasts show decreased mitochondrial respiratory reserve capacity and compensatory increased glycolytic activity. Intriguingly, under basal conditions TMX2 occurs in both reduced and oxidized monomeric form, while it forms a stable dimer under treatment with hydrogen peroxide, recently recognized as a signaling molecule in neural morphogenesis and axonal pathfinding. Exogenous expression of the pathogenic TMX2 variants or of variants with an in vitro mutagenized TRX domain induces a constitutive TMX2 polymerization, mimicking an increased oxidative state. Altogether these data uncover TMX2 as a sensor in the MAM-regulated redox signaling pathway and identify it as a key adaptive regulator of neuronal proliferation, migration, and organization in the developing brain.

We report on a white Afrikaans family from eastern South Africa with three members affected with North Sea progressive myoclonus epilepsy, resulting from a homozygous founder GOSR2 mutation (c.430G>T, p.Gly144Trp). The mutation was identified by exomic sequencing in a research study investigating childhood onset ataxias. All three subjects presented with ataxia, tremor, early gait difficulties, and myoclonic and generalized tonic clonic (GTC) epilepsy. Each patient underwent deep brain stimulation of the caudal Zona Incerta before coming to the attention of the authors. In each case there was a reduction in GTC seizures, and two patients exhibited a reduction in involuntary movements, as evaluated during long‐term follow‐up. In one case there was an improvement in gait and stance when assessed while the stimulation was on.

Adult lung function is highly heritable and 279 genetic loci were recently reported as associated with spirometry-based measures of lung function. Though lung development and function differ between males and females throughout life, there has been no genome-wide study to identify genetic variants with differential effects on lung function in males and females. Here, we present the first genome-wide genotype-by-sex interaction study on four lung function traits in 303,612 participants from the UK Biobank. We detected five SNPs showing genome-wide significant (P−8) interactions with sex on lung function, as well as 21 suggestively significant interactions (P−6). The strongest sex interaction signal came from rs7697189 at 4:145436894 on forced expiratory volume in 1 second (FEV1) (P = 3.15 × 10−15), and was replicated (P = 0.016) in 75,696 individuals in the SpiroMeta consortium. Sex-stratified analyses demonstrated that the minor (C) allele of rs7697189 increased lung function to a greater extent in males than females (untransformed FEV1β = 0.028 [SE 0.0022] litres in males vs β = 0.009 [SE 0.0014] litres in females), and this effect was not accounted for by differential effects on height, smoking or age at puberty. This SNP resides upstream of the gene encoding hedgehog-interacting protein (HHIP) and has previously been reported for association with lung function andHHIPexpression in lung tissue. In our analyses, whileHHIPexpression in lung tissue was significantly different between the sexes with females having higher expression (most significant probeset P=6.90 × 10−6) after adjusting for age and smoking, rs7697189 did not demonstrate sex differential effects on expression. Establishing the mechanism by whichHHIPSNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases, such as chronic obstructive pulmonary disease (COPD), in both sexes.

Reduced lung function predicts mortality and is key to the diagnosis of COPD. In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, one-half of which are new. In combination these variants strongly predict COPD in deeply-phenotyped patient populations. Furthermore, the combined effect of these variants showed generalisability across smokers and never-smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD.

NMDA receptor agonists have been used for many years to generate animal models of polymicrogyria, a malformation of cortical development. Fry et al. identify de novo GRIN1 mutations in eleven patients with severe bilateral polymicrogyria. Polymicrogyria-associated GRIN1 mutations cluster in specific protein domains and significantly alter NMDA receptor function., Polymicrogyria is a malformation of cortical development. The aetiology of polymicrogyria remains poorly understood. Using whole-exome sequencing we found de novo heterozygous missense GRIN1 mutations in 2 of 57 parent-offspring trios with polymicrogyria. We found nine further de novo missense GRIN1 mutations in additional cortical malformation patients. Shared features in the patients were extensive bilateral polymicrogyria associated with severe developmental delay, postnatal microcephaly, cortical visual impairment and intractable epilepsy. GRIN1 encodes GluN1, the essential subunit of the N-methyl-d-aspartate receptor. The polymicrogyria-associated GRIN1 mutations tended to cluster in the S2 region (part of the ligand-binding domain of GluN1) or the adjacent M3 helix. These regions are rarely mutated in the normal population or in GRIN1 patients without polymicrogyria. Using two-electrode and whole-cell voltage-clamp analysis, we showed that the polymicrogyria-associated GRIN1 mutations significantly alter the in vitro activity of the receptor. Three of the mutations increased agonist potency while one reduced proton inhibition of the receptor. These results are striking because previous GRIN1 mutations have generally caused loss of function, and because N-methyl-d-aspartate receptor agonists have been used for many years to generate animal models of polymicrogyria. Overall, our results expand the phenotypic spectrum associated with GRIN1 mutations and highlight the important role of N-methyl-d-aspartate receptor signalling in the pathogenesis of polymicrogyria.

We read with interest “Non-Cardiac Chest Pain: Management in the Ambulatory Clinic setting” (Acute Med 2019;18(3)165-70). It is useful to know about the various musculoskeletal causes of chest pain, especially where specific treatment may be available. We were alarmed, however, about the discussion of pleural pain. This is a common presenting complaint to acute medical teams both on the acute medical take and on ambulatory care units. The authors do not list a differential diagnosis for pain of pleural origin, apart from pulmonary embolism. A chest radiograph is recommended, as well as either a WELL’s (sic) score or D-dimer. Furthermore, bedside lung ultrasound is recommended if the chest radiograph is normal, despite limited availability and training in the UK. Ironically, the reference provided (Aydogdu M et al, Tuberk Toraks 2014;62:12-21) is much more reflective of the modern approach to suspected pulmonary embolism. It discusses the PERC (PE rule-out criteria) which perform very well in people who present with chest pain, and reduces reliance on D-dimer testing. D-dimer testing should never be performed without determining a pre-test probability using validated scores such as the Wells score. As the authors suggest, not all chest pain that worsens with inspiration is “pleuritic”. In addition, substernal pain with a respiratory component does not increase the likelihood of PE, whereas lateral and posterior pain does (Kline J, Thomb Res 2018;163:207-20). Kind regards, Jacob F. de Wolff FRCP Katherine M. Fawcett FRCP Consultant Acute Physicians Northwick Park Hospital

Spinal muscular atrophies (SMAs) are a heterogeneous group of inherited disorders characterized by degeneration of anterior horn cells and progressive muscle weakness. In two unrelated families affected by a distinct form of autosomal-dominant distal SMA initially manifesting with calf weakness, we identified by genetic linkage analysis and exome sequencing a heterozygous missense mutation, c.616T>C (p.Cys206Arg), in F-box protein 38 (FBXO38). FBXO38 is a known coactivator of the transcription factor Krüppel-like factor 7 (KLF7), which regulates genes required for neuronal axon outgrowth and repair. The p.Cys206Arg substitution did not alter the subcellular localization of FBXO38 but did impair KLF7-mediated transactivation of a KLF7-responsive promoter construct and endogenous KLF7 target genes in both heterologously expressing human embryonic kidney 293T cells and fibroblasts derived from individuals with the FBXO38 missense mutation. This transcriptional dysregulation was associated with an impairment of neurite outgrowth in primary motor neurons. Together, these results suggest that a transcriptional regulatory pathway that has a well-established role in axonal development could also be critical for neuronal maintenance and highlight the importance of FBXO38 and KLF7 activity in motor neurons.

Dystonia type 4 (DYT4) was first described in a large family from Heacham in Norfolk with an autosomal dominantly inherited whispering dysphonia, generalized dystonia, and a characteristic hobby horse ataxic gait. We carried out a genetic linkage analysis in the extended DYT4 family that spanned 7 generations from England and Australia, revealing a single LOD score peak of 6.33 on chromosome 19p13.12-13. Exome sequencing in 2 cousins identified a single cosegregating mutation (p.R2G) in the β-tubulin 4a (TUBB4a) gene that was absent in a large number of controls. The mutation is highly conserved in the β-tubulin autoregulatory MREI (methionine-arginine-glutamic acid-isoleucine) domain, highly expressed in the central nervous system, and extensive in vitro work has previously demonstrated that substitutions at residue 2, specifically R2G, disrupt the autoregulatory capability of the wild-type β-tubulin peptide, affirming the role of the cytoskeleton in dystonia pathogenesis. © 2013 American Neurological Association.

Kohlschutter-Tonz syndrome (KTS) is a rare autosomal recessive disorder characterized by amelogenesis imperfecta, psychomotor delay or regression and seizures starting early in childhood. KTS was established as a distinct clinical entity after the first report by Kohlschutter in 1974, and to date, only a total of 20 pedigrees have been reported. The genetic etiology of KTS remained elusive until recently when mutations in ROGDI were independently identified in three unrelated families and in five likely related Druze families. Herein, we report a clinical and genetic study of 10 KTS families. By using a combination of whole exome sequencing, linkage analysis, and Sanger sequencing, we identify novel homozygous or compound heterozygous ROGDI mutations in five families, all presenting with a typical KTS phenotype. The other families, mostly presenting with additional atypical features, were negative for ROGDI mutations, suggesting genetic heterogeneity of atypical forms of the disease. © 2012 Wiley Periodicals, Inc.

Spinocerebellar ataxias (SCA) are a genetically heterogeneous group of neurodegenerative diseases characterised by progressive cerebellar ataxia, dysarthria and oculomotor abnormalities. Recently the prodynorphin (PDYN) gene was identified as the cause of SCA23 in four Dutch families displaying progressive gait and limb ataxia. In this study we aimed to assess the frequency of PDYN gene defects and extend the phenotype of SCA23 patients in a UK ataxia series and also in patients from Greece, Egypt and India. We sequenced the coding and flanking intronic regions of the PDYN gene in a total of 852 ataxia patients, of which 356 were sporadic with no family history, 320 had a positive family history, and 176 probands had a positive family history and at least one family member had also been investigated. We also analysed 190 patients with multiple-system atrophy with cerebellar features (MSA-C), a phenocopy of SCA23. We identified a novel putative pathogenic heterozygous missense variant in the PDYN gene in an early onset SCA patient with an unknown family history. This variant was not present in 570 matched British controls. This is the first study to screen for SCA23 in UK patients and confirms that PDYN mutations are a very rare cause of spinocerebellar ataxia, accounting for ~ 0.1 % of ataxia cases but perhaps with a higher frequency in pure cerebellar ataxia. Given the rarity of PDYN mutations, front-line diagnostic evaluation of UK familial and early onset pure spinocerebellar ataxia patients should focus on other known ataxia genes.

Background Charcot–Marie–Tooth disease (CMT) is a clinically and genetically heterogeneous group of diseases with approximately 45 different causative genes described. The aims of this study were to determine the frequency of different genes in a large cohort of patients with CMT and devise guidelines for genetic testing in practice. Methods The genes known to cause CMT were sequenced in 1607 patients with CMT (425 patients attending an inherited neuropathy clinic and 1182 patients whose DNA was sent to the authors for genetic testing) to determine the proportion of different subtypes in a UK population. Results A molecular diagnosis was achieved in 62.6% of patients with CMT attending the inherited neuropathy clinic; in 80.4% of patients with CMT1 (demyelinating CMT) and in 25.2% of those with CMT2 (axonal CMT). Mutations or rearrangements in PMP22 , GJB1 , MPZ and MFN2 accounted for over 90% of the molecular diagnoses while mutations in all other genes tested were rare. Conclusion Four commonly available genes account for over 90% of all CMT molecular diagnoses; a diagnostic algorithm is proposed based on these results for use in clinical practice. Any patient with CMT without a mutation in these four genes or with an unusual phenotype should be considered for referral for an expert opinion to maximise the chance of reaching a molecular diagnosis.

In 2007, an association of single nucleotide polymorphisms (SNPs) in the fat mass and obesity-associated (FTO) gene region with body mass index (BMI) and risk of obesity was identified in multiple populations, making FTO the first locus unequivocally associated with adiposity. At the time, FTO was a gene of unknown function and it was not known whether these SNPs exerted their effect on adiposity by affecting FTO or neighboring genes. Therefore, this breakthrough association inspired a wealth of in silico, in vitro, and in vivo analyses in model organisms and humans to improve knowledge of FTO function. These studies suggested that FTO plays a role in controlling feeding behavior and energy expenditure. Here, we review the approaches taken that provide a blueprint for the study of other obesity-associated genes in the hope that this strategy will result in increased understanding of the biological mechanisms underlying body weight regulation.

OBJECTIVE— Loss of lipin 1 activity causes lipodystrophy and insulin resistance in the fld mouse, and LPIN1 expression and common genetic variation were recently suggested to influence adiposity and insulin sensitivity in humans. We aimed to conduct a comprehensive association study to clarify the influence of common LPIN1 variation on adiposity and insulin sensitivity in U.K. populations and to examine the role of LPIN1 mutations in insulin resistance syndromes. RESEARCH DESIGN AND METHOD— Twenty-two single nucleotide polymorphisms tagging common LPIN1 variation were genotyped in Medical Research Council (MRC) Ely (n = 1,709) and Hertfordshire (n = 2,901) population-based cohorts. LPIN1 exons, exon/intron boundaries, and 3′ untranslated region were sequenced in 158 patients with idiopathic severe insulin resistance (including 23 lipodystrophic patients) and 48 control subjects. RESULTS— We found no association between LPIN1 single nucleotide polymorphisms and fasting insulin but report a nominal association between rs13412852 and BMI (P = 0.042) in a meta-analysis of 8,504 samples from in-house and publicly available studies. Three rare nonsynonymous variants (A353T, R552K, and G582R) were detected in severely insulin-resistant patients. However, these did not cosegregate with disease in affected families, and Lipin1 protein expression and phosphorylation in patients with variants were indistinguishable from those in control subjects. CONCLUSIONS— Our data do not support a major effect of common LPIN1 variation on metabolic traits and suggest that mutations in LPIN1 are not a common cause of lipodystrophy in humans. The nominal associations with BMI and other metabolic traits in U.K. cohorts require replication in larger cohorts.

A 64 year-old male with metastatic prostate adenocarcinoma presented with bilateral hydronephrosis and renal impairment. Bilateral percutaneous nephrostomy drainage followed by ante-grade stenting was done. Shortly afterwards, the patient developed an extensive left-sided pleural effusion. His serum creatinine rose and he became anuric. Emergency pleural aspiration and later, pleural drainage were performed. Pleural aspirate was diagnostic of urinothorax and non contrast CT scan demonstrated a left reno-pleural fistula. The right stent was removed cystoscopically. The left stent could not be removed cystoscopically and was replaced in an ante grade manner through a fresh percutaneous renal approach. This led to cessation of pleural fluid accumulation. The patient was discharged with bilateral ureteric stents and normal renal function. A month later, he had normal renal function, no hydronephrosis and normal chest x-rays.

Summary: Computational genomics seeks to draw biological inferences from genomic datasets, often by integrating and contextualizing next-generation sequencing data. CGAT provides an extensive suite of tools designed to assist in the analysis of genome scale data from a range of standard file formats. The toolkit enables filtering, comparison, conversion, summarization and annotation of genomic intervals, gene sets and sequences. The tools can both be run from the Unix command line and installed into visual workflow builders, such as Galaxy. Availability: The toolkit is freely available from http://github.com/CGATOxford/cgat Contact: andreas.heger@dpag.ox.ac.uk

Background TRPV4 mutations have been identified in Charcot–Marie–Tooth type 2 (CMT2), scapuloperoneal spinal muscular atrophy and distal hereditary motor neuropathy (dHMN). Objective We aimed to screen the TRPV4 gene in 422 British patients with inherited neuropathy for potentially pathogenic mutations. Methods We sequenced TRPV4 coding regions and splice junctions in 271 patients with CMT2 and 151 patients with dHMN. Mutations were clinically and genetically characterised and screened in ≥345 matched controls. Results 13 missense and nonsense variants were identified, of which five were novel and absent from controls (G20R, E218K, N302Y, Y567X and T701I). N302Y and T701I mutations were present in typical CMT2 cases and are potentially pathogenic based on in silico analyses. G20R was detected in a patient with dHMN and her asymptomatic father and is possibly pathogenic with variable expressivity. The Y567X variant segregated with disease in a family with severe CMT2 but also with a MFN2 mutation reported to cause a mild CMT2 phenotype. Although Y567X caused nonsense mediated mRNA decay, the amount of TRPV4 protein on western blotting of patient lymphoblasts was no different to control. Y567X is therefore unlikely to be pathogenic. E218K is unlikely to be pathogenic based on segregation. Conclusions In this comprehensive analysis of the TRPV4 gene, we identified mutations in TRPV4 likely harbours many missense and nonsense non-pathogenic variants that should be analysed in detail to prove pathogenicity before results are given to patients.

Genome-wide association studies have revealed that common noncoding variants in MTNR1B (encoding melatonin receptor 1B, also known as MT 2) increase type 2 diabetes (T2D) risk. Although the strongest association signal was highly significant (P < 1 - 10 g 20), its contribution to T2D risk was modest (odds ratio (OR) of g1/41.10g1.15). We performed large-scale exon resequencing in 7,632 Europeans, including 2,186 individuals with T2D, and identified 40 nonsynonymous variants, including 36 very rare variants (minor allele frequency (MAF)

OBJECTIVE Wolfram syndrome 1 (WFS1) single nucleotide polymorphisms (SNPs) are associated with risk of type 2 diabetes. In this study we aimed to refine this association and investigate the role of low-frequency WFS1 variants in type 2 diabetes risk. RESEARCH DESIGN AND METHODS For fine-mapping, we sequenced WFS1 exons, splice junctions, and conserved noncoding sequences in samples from 24 type 2 diabetic case and 68 control subjects, selected tagging SNPs, and genotyped these in 959 U.K. type 2 diabetic case and 1,386 control subjects. The same genomic regions were sequenced in samples from 1,235 type 2 diabetic case and 1,668 control subjects to compare the frequency of rarer variants between case and control subjects. RESULTS Of 31 tagging SNPs, the strongest associated was the previously untested 3′ untranslated region rs1046320 (P = 0.008); odds ratio 0.84 and P = 6.59 × 10−7 on further replication in 3,753 case and 4,198 control subjects. High correlation between rs1046320 and the original strongest SNP (rs10010131) (r2 = 0.92) meant that we could not differentiate between their effects in our samples. There was no difference in the cumulative frequency of 82 rare (minor allele frequency [MAF] CONCLUSIONS We identified six highly correlated SNPs that show strong and comparable associations with risk of type 2 diabetes, but further refinement of these associations will require large sample sizes (>100,000) or studies in ethnically diverse populations. Low frequency variants in WFS1 are unlikely to have a large impact on type 2 diabetes risk in white U.K. populations, highlighting the complexities of undertaking association studies with low-frequency variants identified by resequencing.

To the Editor: Glucose uptake into muscle and fat is impaired in insulin-resistant states [1]. In response to insulin, vesicles containing GLUT4 are redistributed from the cell interior to the plasma membrane where they dock and fuse with the plasma membrane (GLUT4 translocation), enabling glucose uptake [2]. In vitro studies suggest that Tre-2 BUB2 CDC16, 1 domain family member 4 (TBC1D4), which is also known as Akt substrate of 160 kDa (AS160), plays an important role in GLUT4 translocation [3] making TBC1D4 an excellent candidate gene for insulin resistance.

Tre-2, BUB2, CDC16, 1 domain family member 4 (TBC1D4) (AS160) is a Rab-GTPase activating protein implicated in insulin-stimulated glucose transporter 4 (GLUT4) translocation in adipocytes and myotubes. To determine whether loss-of-function mutations in TBC1D4 might impair GLUT4 translocation and cause insulin resistance in humans, we screened the coding regions of this gene in 156 severely insulin-resistant patients. A female presenting at age 11 years with acanthosis nigricans and extreme postprandial hyperinsulinemia was heterozygous for a premature stop mutation (R363X) in TBC1D4 . After demonstrating reduced expression of wild-type TBC1D4 protein and expression of the truncated protein in lymphocytes from the proband, we further characterized the biological effects of the truncated protein in 3T3L1 adipocytes. Prematurely truncated TBC1D4 protein tended to increase basal cell membrane GLUT4 levels ( P = 0.053) and significantly reduced insulin-stimulated GLUT4 cell membrane translocation ( P < 0.05). When coexpressed with wild-type TBC1D4, the truncated protein dimerized with full-length TBC1D4, suggesting that the heterozygous truncated variant might interfere with its wild-type counterpart in a dominant negative fashion. Two overweight family members with the mutation also manifested normal fasting glucose and insulin levels but disproportionately elevated insulin levels following an oral glucose challenge. This family provides unique genetic evidence of TBC1D4 involvement in human insulin action.

Mutations at the gene encoding wolframin (WFS1) cause Wolfram syndrome, a rare neurological condition. Associations between single nucleotide polymorphisms (SNPs) at WFS1 and type 2 diabetes have recently been reported. Thus, our aim was to replicate those associations in a northern Swedish case-control study of type 2 diabetes. We also performed a meta-analysis of published and previously unpublished data from Sweden, Finland and France, to obtain updated summary effect estimates.Four WFS1 SNPs (rs10010131, rs6446482, rs752854 and rs734312 [H611R]) were genotyped in a type 2 diabetes case-control study (n = 1,296/1,412) of Swedish adults. Logistic regression was used to assess the association between each WFS1 SNP and type 2 diabetes, following adjustment for age, sex and BMI. We then performed a meta-analysis of 11 studies of type 2 diabetes, comprising up to 14,139 patients and 16,109 controls, to obtain a summary effect estimate for the WFS1 variants.In the northern Swedish study, the minor allele at rs752854 was associated with reduced type 2 diabetes risk [odds ratio (OR) 0.85, 95% CI 0.75-0.96, p=0.010]. Borderline statistical associations were observed for the remaining SNPs. The meta-analysis of the four independent replication studies for SNP rs10010131 and correlated variants showed evidence for statistical association (OR 0.87, 95% CI 0.82-0.93, p=4.5 x 10(-5)). In an updated meta-analysis of all 11 studies, strong evidence of statistical association was also observed (OR 0.89, 95% CI 0.86-0.92; p=4.9 x 10(-11)).In this study of WFS1 variants and type 2 diabetes risk, we have replicated the previously reported associations between SNPs at this locus and the risk of type 2 diabetes.

We studied genes involved in pancreatic β cell function and survival, identifying associations between SNPs in WFS1 and diabetes risk in UK populations that we replicated in an Ashkenazi population and in additional UK studies. In a pooled analysis comprising 9,533 cases and 11,389 controls, SNPs in WFS1 were strongly associated with diabetes risk. Rare mutations in WFS1 cause Wolfram syndrome; using a gene-centric approach, we show that variation in WFS1 also predisposes to common type 2 diabetes.

The inherited neuropathies are a genetically and clinically heterogeneous group of disorders encompassing Charcot-Marie-Tooth disease (CMT), hereditary motor neuropathy (HMN) and hereditary sensory and autonomic neuropathy (HSAN). Mutations in over 45 genes have been implicated in causing the inherited neuropathies, and many more remain unknown. In order to determine the proportion of patients with inherited neuropathy attending our clinic who obtained a genetic diagnosis, we screened all known genes that cause CMT, HMN and HSAN. Of 871 patients attending a specialist neuropathy clinic over a 5-year period, 613 had an inherited neuropathy of whom 409 had CMT. A genetic diagnosis was reached in 59% of patients with CMT; 54% accounted for by chromosome 17 rearrangements or mutations in the commonly available genes PMP22, GJB1, MPZ and MFN2. A genetic diagnosis was achieved in 79% of patients with CMT1 and 35% of CMT2. The chromosome 17 duplication only accounted for 38% of CMT, less than reported in other populations; this reflects a bias in our clinic towards patients with rarer forms of CMT. We give details of the individual contributions the rarer genes made to our cohort. These results suggest that the majority of patients with CMT1 should obtain a genetic diagnosis. Mutations in the rarer known genes account for a small proportion of patients with CMT, and many genes for CMT2 remain unknown. We suggest an algorithm for genetic testing in patients with CMT in the UK.