DeepPheWAS: an R package for phenotype generation and association analysis for phenome-wide association studies.

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.

TMX2 Is a Crucial Regulator of Cellular Redox State, and Its Dysfunction Causes Severe Brain Developmental Abnormalities.

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.

Missense Mutations of the Pro65 Residue of PCGF2 Cause a Recognizable Syndrome Associated with Craniofacial, Neurological, Cardiovascular, and Skeletal Features.

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.

SCN3A-Related Neurodevelopmental Disorder: A Spectrum of Epilepsy and Brain Malformation.

OBJECTIVE: 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. METHODS: Patients were ascertained via an international collaborative network. We compared sodium channels containing wild-type versus variant Nav1.3 subunits coexpressed with ß1 and ß2 subunits using whole-cell voltage clamp electrophysiological recordings in a heterologous mammalian system (HEK-293T cells). RESULTS: 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 to 6 of domains II to 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. INTERPRETATION: 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 >75% of affected individuals. Gain of function at the channel level in developing neurons is likely an important mechanism of disease pathogenesis. ANN NEUROL 2020;88:348-362.

De novo mutations in GRIN1 cause extensive bilateral polymicrogyria.

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.

A dominant mutation in FBXO38 causes distal spinal muscular atrophy with calf predominance.

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.

De novo point mutations in patients diagnosed with ataxic cerebral palsy.

Cerebral palsy is a sporadic disorder with multiple likely aetiologies, but frequently considered to be caused by birth asphyxia. Genetic investigations are rarely performed in patients with cerebral palsy and there is little proven evidence of genetic causes. As part of a large project investigating children with ataxia, we identified four patients in our cohort with a diagnosis of ataxic cerebral palsy. They were investigated using either targeted next generation sequencing or trio-based exome sequencing and were found to have mutations in three different genes, KCNC3, ITPR1 and SPTBN2. All the mutations were de novo and associated with increased paternal age. The mutations were shown to be pathogenic using a combination of bioinformatics analysis and in vitro model systems. This work is the first to report that the ataxic subtype of cerebral palsy can be caused by de novo dominant point mutations, which explains the sporadic nature of these cases. We conclude that at least some subtypes of cerebral palsy may be caused by de novo genetic mutations and patients with a clinical diagnosis of cerebral palsy should be genetically investigated before causation is ascribed to perinatal asphyxia or other aetiologies.

Common variants in WFS1 confer risk of type 2 diabetes.

We studied genes involved in pancreatic beta 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.

CGAT: computational genomics analysis toolkit.

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.

The genetics of obesity: FTO leads the way.

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.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk.

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 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 showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores 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.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation.

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Comprehensive analysis of the TRPV4 gene in a large series of inherited neuropathies and controls.

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 <1% of patients with CMT2/dHMN. We found that 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.

A truncation mutation in TBC1D4 in a family with acanthosis nigricans and postprandial hyperinsulinemia.

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.

Exome-wide analysis of copy number variation shows association of the human leukocyte antigen region with asthma in UK Biobank.

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.

Pleiotropic associations of heterozygosity for the SERPINA1 Z allele in the UK Biobank.

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (ß=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including 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 height-adjusted forced expiratory volume in 1 s (FEV1) (ß=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (ß=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

Variants associated with HHIP expression have sex-differential effects on lung function.

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<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 ß=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.

New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries.

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 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.