SLC6A1 variant pathogenicity, molecular function and phenotype: a genetic and clinical analysis

Genetic variants in the SLC6A1 gene can cause a broad phenotypic disease spectrum by altering the protein function. Thus, systematically curated clinically relevant genotype-phenotype associations are needed to understand the disease mechanism and improve therapeutic decision-making. We aggregated genetic and clinical data from 172 individuals with likely pathogenic/pathogenic (lp/p) SLC6A1 variants and functional data for 184 variants (14.1% lp/p). Clinical and functional data were available for a subset of 126 individuals. We explored the potential associations of variant positions on the GAT1 3D structure with variant pathogenicity, altered molecular function and phenotype severity using bioinformatic approaches. The GAT1 transmembrane domains 1, 6 and extracellular loop 4 (EL4) were enriched for patient over population variants. Across functionally tested missense variants (n = 156), the spatial proximity from the ligand was associated with loss-of-function in the GAT1 transporter activity. For variants with complete loss of in vitro GABA uptake, we found a 4.6-fold enrichment in patients having severe disease versus non-severe disease (P = 2.9 × 10-3, 95% confidence interval: 1.5-15.3). In summary, we delineated associations between the 3D structure and variant pathogenicity, variant function and phenotype in SLC6A1-related disorders. This knowledge supports biology-informed variant interpretation and research on GAT1 function. All our data can be interactively explored in the SLC6A1 portal (https://slc6a1-portal.broadinstitute.org/).

Rare CIDEC coding variants enriched in age-related macular degeneration patients with small low-luminance deficit cause lipid droplet and fat storage defects

BACKGROUND

The basis of Age-related macular degeneration (AMD) genetic risk has been well documented; however, few studies have looked at genetic biomarkers of disease progression or treatment response within advanced AMD patients. Here we report the first genome-wide analysis of genetic determinants of low-luminance vision deficit (LLD), which is seen as predictive of visual acuity loss and anti-VEGF treatment response in neovascular AMD patients.

METHODS

AMD patients were separated into small- and large-LLD groups for comparison and whole genome sequencing was performed. Genetic determinants of LLD were assessed by common and rare variant genetic analysis. Follow-up functional analysis of rare coding variants identified by the burden test was then performed in vitro.

RESULTS

We identified four coding variants in the CIDEC gene. These rare variants were only present in patients with a small LLD, which has been previously shown to indicate better prognosis and better anti-VEGF treatment response. Our in vitro functional characterization of these CIDEC alleles revealed that all decrease the binding affinity between CIDEC and the lipid droplet fusion effectors PLIN1, RAB8A and AS160. The rare CIDEC alleles all cause a hypomorphic defect in lipid droplet fusion and enlargement, resulting in a decreased fat storage capability in adipocytes.

CONCLUSIONS

As we did not detect CIDEC expression in the ocular tissue affected by AMD, our results suggest that the CIDEC variants do not play a direct role in the eye and influence low-luminance vision deficit via an indirect and systemic effect related to fat storage capacity.

Identifying therapeutic drug targets using bidirectional effect genes

Prioritizing genes for translation to therapeutics for common diseases has been challenging. Here, we propose an approach to identify drug targets with high probability of success by focusing on genes with both gain of function (GoF) and loss of function (LoF) mutations associated with opposing effects on phenotype (Bidirectional Effect Selected Targets, BEST). We find 98 BEST genes for a variety of indications. Drugs targeting those genes are 3.8-fold more likely to be approved than non-BEST genes. We focus on five genes (IGF1R, NPPC, NPR2, FGFR3, and SHOX) with evidence for bidirectional effects on stature. Rare protein-altering variants in those genes result in significantly increased risk for idiopathic short stature (ISS) (OR = 2.75, p = 3.99 × 10-8). Finally, using functional experiments, we demonstrate that adding an exogenous CNP analog (encoded by NPPC) rescues the phenotype, thus validating its potential as a therapeutic treatment for ISS. Our results show the value of looking for bidirectional effects to identify and validate drug targets.

Influence of genetic copy number variants of the human GLUT3 glucose transporter gene SLC2A3 on protein expression, glycolysis and rheumatoid arthritis risk: A genetic replication study

Objectives

The gene encoding glucose transporter 3 (GLUT3, SLC2A3) is present in the human population at variable copy number. An overt disease phenotype of SLC2A3 copy number variants has not been reported; however, deletion of SLC2A3 has been previously reported to protect carriers from rheumatoid arthritis, implicating GLUT3 as a therapeutic target in rheumatoid arthritis. Here we aim to perform functional analysis of GLUT3 copy number variants in immune cells, and test the reported protective association of the GLUT3 copy number variants for rheumatoid arthritis in a genetic replication study.

Methods

Cells from genotyped healthy controls were analyzed for SLC2A3/GLUT3 expression and glycolysis capacity. We genotyped the SLC2A3 copy number variant in four independent cohorts of rheumatoid arthritis and controls and one cohort of multiple sclerosis and controls.

Results

Heterozygous deletion of SLC2A3 correlates directly with expression levels of GLUT3 and influences glycolysis rates in the human immune system. The frequency of the SLC2A3 copy number variant is not different between rheumatoid arthritis, multiple sclerosis and control groups.

Conclusions

Despite a robust SLC2A3 gene copy number dependent phenotype, our study of large groups of rheumatoid arthritis cases and controls provides no evidence for rheumatoid arthritis disease protection in deletion carriers. These data emphasize the importance of well powered replication studies to confirm or refute genetic associations, particularly for relatively rare variants.

•

T cell and macrophage expression of SLC2A3/GLUT3 correlates to SLC2A3 gene copy number in a dose dependent manner.

•

Glycolysis rates are reduced in individuals harboring a deletion of the GLUT3 gene SLC2A3

•

Deletion of SLC2A3 is not associated with protection from rheumatoid arthritis

•

Deletion of SLC2A3 is not associated with risk for multiple sclerosis

•

GLUT3 is not a viable therapeutic target for RA as previously proposed based on a protective association of SLC2A3 deletion.

LACC1 Regulates TNF and IL-17 in Mouse Models of Arthritis and Inflammation

Both common and rare genetic variants of laccase domain-containing 1 (LACC1, previously C13orf31) are associated with inflammatory bowel disease, leprosy, Behcet disease, and systemic juvenile idiopathic arthritis. However, the functional relevance of these variants is unclear. In this study, we use LACC1-deficient mice to gain insight into the role of LACC1 in regulating inflammation. Following oral administration of Citrobacter rodentium, LACC1 knockout (KO) mice had more severe colon lesions compared with wildtype (WT) controls. Immunization with collagen II, a collagen-induced arthritis (CIA) model, resulted in an accelerated onset of arthritis and significantly worse arthritis and inflammation in LACC1 KO mice. Similar results were obtained in a mannan-induced arthritis model. Serum and local TNF in CIA paws and C. rodentium colons were significantly increased in LACC1 KO mice compared with WT controls. The percentage of IL-17A-producing CD4⁺ T cells was elevated in LACC1 KO mice undergoing CIA as well as aged mice compared with WT controls. Neutralization of IL-17, but not TNF, prevented enhanced mannan-induced arthritis in LACC1 KO mice. These data provide new mechanistic insight into the function of LACC1 in regulating TNF and IL-17 during inflammatory responses. We hypothesize that these effects contribute to immune-driven pathologies observed in individuals carrying LACC1 variants.

Paired Immunoglobulin-like Type 2 Receptor Alpha G78R variant alters ligand binding and confers protection to Alzheimer's disease

Paired Immunoglobulin-like Type 2 Receptor Alpha (PILRA) is a cell surface inhibitory receptor that recognizes specific O-glycosylated proteins and is expressed on various innate immune cell types including microglia. We show here that a common missense variant (G78R, rs1859788) of PILRA is the likely causal allele for the confirmed Alzheimer’s disease risk locus at 7q21 (rs1476679). The G78R variant alters the interaction of residues essential for sialic acid engagement, resulting in >50% reduced binding for several PILRA ligands including a novel ligand, complement component 4A, and herpes simplex virus 1 (HSV-1) glycoprotein B. PILRA is an entry receptor for HSV-1 via glycoprotein B, and macrophages derived from R78 homozygous donors showed significantly decreased levels of HSV-1 infection at several multiplicities of infection compared to homozygous G78 macrophages. We propose that PILRA G78R protects individuals from Alzheimer’s disease risk via reduced inhibitory signaling in microglia and reduced microglial infection during HSV-1 recurrence.

Alzheimer’s disease (AD) is a devastating neurodegenerative disorder resulting from a complex interaction of environmental and genetic risk factors. Despite considerable progress in defining the genetic component of AD risk, understanding the biology of common variant associations is a challenge. We find that PILRA G78R (rs1859788) is the likely AD risk variant from the 7q21 locus (rs1476679) and PILRA G78R reduces PILRA endogenous and exogenous ligand binding. Our study highlights a new immune signaling axis in AD and suggests a role for exogenous ligands (HSV-1). Further, we have identified that reduced function of a negative regulator of microglia and neutrophils is protective from AD risk, providing a new candidate therapeutic target.

Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions

Decoding the information in mRNA during protein synthesis relies on tRNA adaptors, the abundance of which can affect the decoding rate and translation efficiency. To determine whether cells alter tRNA abundance to selectively regulate protein expression, we quantified changes in the abundance of individual tRNAs at different time points in response to diverse stress conditions in Saccharomyces cerevisiae. We found that the tRNA pool was dynamic and rearranged in a manner that facilitated selective translation of stress-related transcripts. Through genomic analysis of multiple data sets, stochastic simulations, and experiments with designed sequences of proteins with identical amino acids but altered codon usage, we showed that changes in tRNA abundance affected protein expression independently of factors such as mRNA abundance. We suggest that cells alter their tRNA abundance to selectively affect the translation rates of specific transcripts to increase the amounts of required proteins under diverse stress conditions.

Inflammatory Bowel Disease Susceptibility Gene C1ORF106 Regulates Intestinal Epithelial Permeability

Intestinal epithelial cells form a physical barrier that is tightly regulated to control intestinal permeability. Proinflammatory cytokines, such as TNF-α, increase epithelial permeability through disruption of epithelial junctions. The regulation of the epithelial barrier in inflammatory gastrointestinal disease remains to be fully characterized. In this article, we show that the human inflammatory bowel disease genetic susceptibility gene C1ORF106 plays a key role in regulating gut epithelial permeability. C1ORF106 directly interacts with cytohesins to maintain functional epithelial cell junctions. C1orf106-deficient mice are hypersensitive to TNF-α-induced increase in epithelial permeability, and this is associated with increased diarrhea. This study identifies C1ORF106 as an epithelial cell junction protein, and the loss of C1ORF106 augments TNF-α-induced intestinal epithelial leakage and diarrhea that may play a critical role in the development of inflammatory bowel disease.

Estimating the human mutation rate from autozygous segments reveals population differences in human mutational processes

Heterozygous mutations within homozygous sequences descended from a recent common ancestor offer a way to ascertain de novo mutations across multiple generations. Using exome sequences from 3222 British-Pakistani individuals with high parental relatedness, we estimate a mutation rate of 1.45 ± 0.05 × 10-8 per base pair per generation in autosomal coding sequence, with a corresponding non-crossover gene conversion rate of 8.75 ± 0.05 × 10-6 per base pair per generation. This is at the lower end of exome mutation rates previously estimated in parent-offspring trios, suggesting that post-zygotic mutations contribute little to the human germ-line mutation rate. We find frequent recurrence of mutations at polymorphic CpG sites, and an increase in C to T mutations in a 5' CCG 3' to 5' CTG 3' context in the Pakistani population compared to Europeans, suggesting that mutational processes have evolved rapidly between human populations.Estimates of human mutation rates differ substantially based on the approach. Here, the authors present a multi-generational estimate from the autozygous segment in a non-European population that gives insight into the contribution of post-zygotic mutations and population-specific mutational processes.

The Ox40/Ox40 Ligand Pathway Promotes Pathogenic Th Cell Responses, Plasmablast Accumulation, and Lupus Nephritis in NZB/W F1 Mice

Ox40 ligand (Ox40L) locus genetic variants are associated with the risk for systemic lupus erythematosus (SLE); however, it is unclear how Ox40L contributes to SLE pathogenesis. In this study, we evaluated the contribution of Ox40L and its cognate receptor, Ox40, using in vivo agonist and antagonist approaches in the NZB × NZW (NZB/W) F1 mouse model of SLE. Ox40 was highly expressed on several CD4 Th cell subsets in the spleen and kidney of diseased mice, and expression correlated with disease severity. Treatment of aged NZB/W F1 mice with agonist anti-Ox40 mAbs potently exacerbated renal disease, which was accompanied by activation of kidney-infiltrating T cells and cytokine production. The agonist mAbs also induced activation and inflammatory gene expression in splenic CD4 T cells, including IFN-regulated genes, increased the number of follicular helper T cells and plasmablasts in the spleen, and led to elevated levels of serum IgM and enhanced renal glomerular IgM deposition. In a type I IFN-accelerated lupus model, treatment with an antagonist Ox40:Fc fusion protein significantly delayed the onset of severe proteinuria and improved survival. These data support the hypothesis that the Ox40/Ox40L pathway drives cellular and humoral autoimmune responses during lupus nephritis in NZB/W F1 mice and emphasize the potential clinical value of targeting this pathway in human lupus.

An Organismal CNV Mutator Phenotype Restricted to Early Human Development

De novo copy number variants (dnCNVs) arising at multiple loci in a personal genome have usually been considered to reflect cancer somatic genomic instabilities. We describe a multiple dnCNV (MdnCNV) phenomenon in which individuals with genomic disorders carry five to ten constitutional dnCNVs. These CNVs originate from independent formation incidences, are predominantly tandem duplications or complex gains, exhibit breakpoint junction features reminiscent of replicative repair, and show increased de novo point mutations flanking the rearrangement junctions. The active CNV mutation shower appears to be restricted to a transient perizygotic period. We propose that a defect in the CNV formation process is responsible for the "CNV-mutator state," and this state is dampened after early embryogenesis. The constitutional MdnCNV phenomenon resembles chromosomal instability in various cancers. Investigations of this phenomenon may provide unique access to understanding genomic disorders, structural variant mutagenesis, human evolution, and cancer biology.

geneAttribution: trait agnostic identification of candidate genes associated with noncoding variation

Motivation

We have developed geneAttribution, an R package that assigns candidate causal gene(s) to a risk variant identified by a genetic association study such as a GWAS. The method combines user-supplied functional annotation such as expression quantitative trait loci (eQTL) or Hi-C genome conformation data and reports the most likely candidate genes. In the absence of annotation data, geneAttribution relies on the distances between the genes and the input variant.

Availability and Implementation

The package is freely available from http://www.bioconductor.org/. A quick-start vignette is included with the package.

Regulation of T-cell receptor downmodulation by DENND1B in TH2 cells and allergic disease

The DENN-domain is an evolutionary conserved protein module found in all eukaryotes and serves as an exchange factor for Rab-GTPases to regulate diverse cellular functions. Variants in DENND1B are associated with development of childhood asthma and other immune disorders. To understand how DENND1B may contribute to diseases, Dennd1b−/− mice were generated and exhibit hyper-allergic responses following antigen challenge. Dennd1b−/− TH2, but not other TH cells, exhibit delayed receptor-induced TCR downmodulation, enhanced TCR signaling, and increased production of effector cytokines. As DENND1B interacts with AP-2 and Rab35, TH2 cells deficient in AP-2 or Rab35 also exhibit enhanced TCR-mediated effector functions. Moreover, human TH2 cells carrying asthma-associated DENND1B variants express less DENND1B and phenocopy Dennd1b−/− TH2 cells. These results provide a molecular basis for how DENND1B, a previously unrecognized regulator of TCR downmodulation in TH2 cells, contributes to asthma pathogenesis and how DENN-domain containing proteins may contribute to other human disorders.

Thinking differently about lupus

A search for the genetic causes of an autoimmune disease called systemic lupus erythematosus reveals a new twist on an old story.

Regulation of T Cell Receptor Signaling by DENND1B in TH2 Cells and Allergic Disease

The DENN domain is an evolutionary conserved protein module found in all eukaryotes and serves as an exchange factor for Rab-GTPases to regulate diverse cellular functions. Variants in DENND1B are associated with development of childhood asthma and other immune disorders. To understand how DENND1B may contribute to human disease, Dennd1b(-/-) mice were generated and exhibit hyper-allergic responses following antigen challenge. Dennd1b(-/-) TH2, but not other TH cells, exhibit delayed receptor-induced T cell receptor (TCR) downmodulation, enhanced TCR signaling, and increased production of effector cytokines. As DENND1B interacts with AP-2 and Rab35, TH2 cells deficient in AP-2 or Rab35 also exhibit enhanced TCR-mediated effector functions. Moreover, human TH2 cells carrying asthma-associated DENND1B variants express less DENND1B and phenocopy Dennd1b(-/-) TH2 cells. These results provide a molecular basis for how DENND1B, a previously unrecognized regulator of TCR downmodulation in TH2 cells, contributes to asthma pathogenesis and how DENN-domain-containing proteins may contribute to other human disorders.

Timing, rates and spectra of human germline mutation

Germline mutations are a driving force behind genome evolution and genetic disease. We investigated genome-wide mutation rates and spectra in multi-sibling families. Mutation rate increased with paternal age in all families, but the number of additional mutations per year differed more than two-fold between families. Meta-analysis of 6,570 mutations showed that germline methylation influences mutation rates. In contrast to somatic mutations, we found remarkable consistency of germline mutation spectra between the sexes and at different paternal ages. 3.8% of mutations were mosaic in the parental germline, resulting in 1.3% of mutations being shared between siblings. The number of these shared mutations varied significantly between families. Our data suggest that the mutation rate per cell division is higher during both early embryogenesis and differentiation of primordial germ cells, but is reduced substantially during post-pubertal spermatogenesis. These findings have important consequences for the recurrence risks of disorders caused by de novo mutations.

Spial: analysis of subtype-specific features in multiple sequence alignments of proteins

Motivation: Spial (Specificity in alignments) is a tool for the comparative analysis of two alignments of evolutionarily related sequences that differ in their function, such as two receptor subtypes. It highlights functionally important residues that are either specific to one of the two alignments or conserved across both alignments. It permits visualization of this information in three complementary ways: by colour-coding alignment positions, by sequence logos and optionally by colour-coding the residues of a protein structure provided by the user. This can aid in the detection of residues that are involved in the subtype-specific interaction with a ligand, other proteins or nucleic acids. Spial may also be used to detect residues that may be post-translationally modified in one of the two sets of sequences.

Availability:http://www.mrc-lmb.cam.ac.uk/genomes/spial/; supplementary information is available at http://www.mrc-lmb.cam.ac.uk/genomes/spial/help.html

Contact:ajv@mrc-lmb.cam.ac.uk

Transcriptional control of the quorum sensing response in yeast

Quorum sensing is a process of intercellular communication. It allows individual cells to assess population density and to co-ordinate behaviour by secreting and sensing communication molecules. In the yeast Saccharomyces cerevisiae, the communication molecules are the aromatic alcohols tryptophol and phenylethanol, and quorum sensing regulates the transition between the solitary yeast form and the filamentous form. Though it is known that addition of these communication molecules to yeast cultures causes large changes in gene expression, how these changes are orchestrated and whether this system is conserved in related fungal species is still unknown. In this work, by employing an integrated computational approach that makes use of large-scale genomics datasets, such as ChIP-ChIP and expression analysis upon deletion and over-expression of transcriptional factors, we predict CAT8 and MIG1 as key transcriptional regulators that control the differential expression of the genes affected by aromatic alcohol communication. In addition, through a comparative genomic analysis involving 31 fungal species, we show that the S. cerevisiae quorum sensing system is a recent evolutionary innovation and that the genes which are differentially expressed upon treatment with these molecules are distributed across the genome in a highly non-random manner. The identified transcription factors will aid in further unravelling the molecular mechanisms of S. cerevisiae quorum sensing and may facilitate the engineering of regulatory circuits for applications such as the expression of heterologous proteins via aromatic alcohols.

Genomic analysis reveals a tight link between transcription factor dynamics and regulatory network architecture

Although several studies have provided important insights into the general principles of biological networks, the link between network organization and the genome-scale dynamics of the underlying entities (genes, mRNAs, and proteins) and its role in systems behavior remain unclear. Here we show that transcription factor (TF) dynamics and regulatory network organization are tightly linked. By classifying TFs in the yeast regulatory network into three hierarchical layers (top, core, and bottom) and integrating diverse genome-scale datasets, we find that the TFs have static and dynamic properties that are similar within a layer and different across layers. At the protein level, the top-layer TFs are relatively abundant, long-lived, and noisy compared with the core- and bottom-layer TFs. Although variability in expression of top-layer TFs might confer a selective advantage, as this permits at least some members in a clonal cell population to initiate a response to changing conditions, tight regulation of the core- and bottom-layer TFs may minimize noise propagation and ensure fidelity in regulation. We propose that the interplay between network organization and TF dynamics could permit differential utilization of the same underlying network by distinct members of a clonal cell population.

Chemogenomics and biotechnology

A robust knowledge of the interactions between small molecules and specific proteins aids the development of new biotechnological tools and the identification of new drug targets, and can lead to specific biological insights. Such knowledge can be obtained through chemogenomic screens. In these screens, each small molecule from a chemical library is applied to each cell type from a library of cells, and the resulting phenotypes are recorded. Chemogenomic screens have recently become very common and will continue to generate large amounts of data. The interpretation of this data will occupy biologists and chemists alike for some time to come. This review discusses methods for the acquisition and interpretation of chemogenomic data, in addition to possible applications of chemogenomics in biotechnology.