Analysis of stop-gain and frameshift variants in human innate immunity genes.

Loss-of-function variants in innate immunity genes are associated with Mendelian disorders in the form of primary immunodeficiencies. Recent resequencing projects report that stop-gains and frameshifts are collectively prevalent in humans and could be responsible for some of the inter-individual variability in innate immune response. Current computational approaches evaluating loss-of-function in genes carrying these variants rely on gene-level characteristics such as evolutionary conservation and functional redundancy across the genome. However, innate immunity genes represent a particular case because they are more likely to be under positive selection and duplicated. To create a ranking of severity that would be applicable to innate immunity genes we evaluated 17,764 stop-gain and 13,915 frameshift variants from the NHLBI Exome Sequencing Project and 1,000 Genomes Project. Sequence-based features such as loss of functional domains, isoform-specific truncation and nonsense-mediated decay were found to correlate with variant allele frequency and validated with gene expression data. We integrated these features in a Bayesian classification scheme and benchmarked its use in predicting pathogenic variants against Online Mendelian Inheritance in Man (OMIM) disease stop-gains and frameshifts. The classification scheme was applied in the assessment of 335 stop-gains and 236 frameshifts affecting 227 interferon-stimulated genes. The sequence-based score ranks variants in innate immunity genes according to their potential to cause disease, and complements existing gene-based pathogenicity scores. Specifically, the sequence-based score improves measurement of functional gene impairment, discriminates across different variants in a given gene and appears particularly useful for analysis of less conserved genes.

Splice variant-specific stabilization of JNKs by IB1/JIP1.

Islet-Brain 1 (IB1) (also called JNK-interacting protein 1; JIP1) is a scaffold protein that tethers components of the JNK mitogen-activated protein kinase pathway inducing a modulation of the activity and the target specificity of the JNK kinases. Dysfunctions in IB1 have been associated with diseases such as early type II diabetes. To gain more insight in the functions of IB1, its ability to modulate the expression levels of the various JNK proteins was assessed. Each of the three JNK genes gives rise to several splice variants encoding short or long proteins. The expression levels of the short JNK proteins, but not of the long variants, were systematically higher in rat tissues and in transformed cell lines expressing high IB1 levels compared to tissues and cells with no or low IB1 expression. HEK293 cells bearing a tetracycline-inducible IB1 construct showed a specific increase of the short JNK endogenous splice variants in the presence of tetracycline. The augmented expression level of the short JNK splice variants induced by IB1 resulted from an increased stability towards degradation. Modulation of the stability of specific JNK splice variants represents therefore a newly identified mechanism used by IB1 to regulate the JNK MAPK pathway.

Association study of 182 candidate genes in anorexia nervosa.

We performed association studies with 5,151 SNPs that were judged as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN) based on location under reported linkage peaks, previous results in the literature (182 candidate genes), brain expression, biological plausibility, and estrogen responsivity. We employed a case-control design that tested each SNP individually as well as haplotypes derived from these SNPs in 1,085 case individuals with AN diagnoses and 677 control individuals. We also performed separate association analyses using three increasingly restrictive case definitions for AN: all individuals with any subtype of AN (All AN: n = 1,085); individuals with AN with no binge eating behavior (AN with No Binge Eating: n = 687); and individuals with the restricting subtype of AN (Restricting AN: n = 421). After accounting for multiple comparisons, there were no statistically significant associations for any individual SNP or haplotype block with any definition of illness. These results underscore the importance of large samples to yield appropriate power to detect genotypic differences in individuals with AN and also motivate complementary approaches involving Genome-Wide Association (GWA) studies, Copy Number Variation (CNV) analyses, sequencing-based rare variant discovery assays, and pathway-based analysis in order to make up for deficiencies in traditional candidate gene approaches to AN.

Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila.

Ionotropic glutamate receptors (iGluRs) mediate neuronal communication at synapses throughout vertebrate and invertebrate nervous systems. We have characterized a family of iGluR-related genes in Drosophila, which we name ionotropic receptors (IRs). These receptors do not belong to the well-described kainate, AMPA, or NMDA classes of iGluRs, and they have divergent ligand-binding domains that lack their characteristic glutamate-interacting residues. IRs are expressed in a combinatorial fashion in sensory neurons that respond to many distinct odors but do not express either insect odorant receptors (ORs) or gustatory receptors (GRs). IR proteins accumulate in sensory dendrites and not at synapses. Misexpression of IRs in different olfactory neurons is sufficient to confer ectopic odor responsiveness. Together, these results lead us to propose that the IRs comprise a novel family of chemosensory receptors. Conservation of IR/iGluR-related proteins in bacteria, plants, and animals suggests that this receptor family represents an evolutionarily ancient mechanism for sensing both internal and external chemical cues.

Reduced IFNλ4 activity is associated with improved HCV clearance and reduced expression of interferon-stimulated genes

Hepatitis C virus (HCV) infections are the major cause of chronic liver disease, cirrhosis and hepatocellular carcinoma worldwide. Both spontaneous and treatment-induced clearance of HCV depend on genetic variation within the interferon-lambda locus, but until now no clear causal relationship has been established. Here we demonstrate that an amino-acid substitution in the IFNλ4 protein changing a proline at position 70 to a serine (P70S) substantially alters its antiviral activity. Patients harbouring the impaired IFNλ4-S70 variant display lower interferon-stimulated gene (ISG) expression levels, better treatment response rates and better spontaneous clearance rates, compared with patients coding for the fully active IFNλ4-P70 variant. Altogether, these data provide evidence supporting a role for the active IFNλ4 protein as the driver of high hepatic ISG expression as well as the cause of poor HCV clearance.

Transcription factor CTF1 acts as a chromatin domain boundary that shields human telomeric genes from silencing.

Telomeres are associated with chromatin-mediated silencing of genes in their vicinity. However, how epigenetic markers mediate mammalian telomeric silencing and whether specific proteins may counteract this effect are not known. We evaluated the ability of CTF1, a DNA- and histone-binding transcription factor, to prevent transgene silencing at human telomeres. CTF1 was found to protect a gene from silencing when its DNA-binding sites were interposed between the gene and the telomeric extremity, while it did not affect a gene adjacent to the telomere. Protein fusions containing the CTF1 histone-binding domain displayed similar activities, while mutants impaired in their ability to interact with the histone did not. Chromatin immunoprecipitation indicated the propagation of a hypoacetylated histone structure to various extents depending on the telomere. The CTF1 fusion protein was found to recruit the H2A.Z histone variant at the telomeric locus and to restore high histone acetylation levels to the insulated telomeric transgene. Histone lysine trimethylations were also increased on the insulated transgene, indicating that these modifications may mediate expression rather than silencing at human telomeres. Overall, these results indicate that transcription factors can act to delimit chromatin domain boundaries at mammalian telomeres, thereby blocking the propagation of a silent chromatin structure.

The novel ss469415590 variant predicts virological response to therapy in patients with chronic hepatitis C virus type 1 infection.

BACKGROUND: A novel dinucleotide variant TT/∆G (ss469415590) has been associated with hepatitis C virus clearance. AIM: To assess the role of the ss469415590 variant, compared with the known IL28B polymorphisms (rs8099917, rs12979860 and rs12980275) for predicting virological response to therapy in chronic hepatitis C, and its association with the CXCL10 chemokine serum levels - a surrogate marker of interferon-stimulated genes activation. METHODS: Multivariate analysis of factors predicting rapid and sustained virological response in 280 consecutive, treatment-naïve, nondiabetic, Caucasian patients with chronic hepatitis C treated with peginterferon alpha and ribavirin. RESULTS: In hepatitis C virus genotype 1, the OR (95% CI) for rapid and sustained virological response for the wild-type ss469415590 TT was 9.88 (1.99-48.99) and 7.25 (1.91-27.51), respectively, similar to those found for rs12979860 CC [9.55 (1.93-47.37) and 6.30 (1.71-23.13)] and for rs12980275 AA [9.62 (1.94-47.77] and 7.83 (2.02-30.34)], but higher than for rs8099917 TT [4.8 (1.73-13.33) and 4.75 (2.05-10.98)]. In hepatitis C virus genotype 1, mean (SD) CXCL10 levels in patients with the TT/TT, TT/∆G and ∆G/∆G variants were, respectively, 355.1 (240.6), 434.4 (247.4) and 569.9 (333.3) (P = 0.04). In patients with genotypes 2 and 3 no significant association was found for TT/∆G with viral response. The predictive value of ss469415590 was stronger in patients with advanced fibrosis. CONCLUSIONS: The novel IL28B variants at marker ss469415590 predict response to IFN alpha in chronic hepatitis C patients, especially in those with advanced fibrosis. Their determination may be superior to that of known IL28B variants for patient management using IFN-based regimens.

A natural structural variant of the mouse TCR beta-chain displays intrinsic receptor function and antigen specificity.

The Cbeta0 alternate cassette exon is located between the Jbeta1 and Cbeta1 genes in the mouse TCR beta-locus. In T cells with a VDJbeta1 rearrangement, the Cbeta0 exon may be included in TCRbeta transcripts (herein called TCRbeta-Cbeta0 transcripts), potentially inserting an additional 24 aa between the V and C domains of the TCR beta-chain. These TCRbeta splice isoforms may be differentially regulated after Ag activation, because we detected TCRbeta-Cbeta0 transcripts in a high proportion (>60%) of immature and mature T cells having VDJbeta1 rearrangements but found a substantially reduced frequency (<35%) of TCRbeta-Cbeta0 expression among CD8 T cells selected by Ag in vivo. To study the potential activity of the TCRbeta-Cbeta0 splice variant, we cloned full-length TCR cDNAs by single-cell RT-PCR into retroviral expression vectors. We found that the TCRbeta-Cbeta0 splice isoform can function during an early stage of T cell development normally dependent on TCR beta-chain expression. We also demonstrate that T hybridoma-derived cells expressing a TCRbeta-Cbeta0 isoform together with the clonally associated TCR alpha-chain recognize the same cognate peptide-MHC ligand as the corresponding normal alphabetaTCR. This maintenance of receptor function and specificity upon insertion of the Cbeta0 peptide cassette signifies a remarkable adaptability for the TCR beta-chain, and our findings open the possibility that this splice isoform may function in vivo.

Bioinformatics-driven identification and examination of candidate genes for non-alcoholic fatty liver disease.

ObjectiveCandidate genes for non-alcoholic fatty liver disease (NAFLD) identified by a bioinformatics approach were examined for variant associations to quantitative traits of NAFLD-related phenotypes.Research Design and MethodsBy integrating public database text mining, trans-organism protein-protein interaction transferal, and information on liver protein expression a protein-protein interaction network was constructed and from this a smaller isolated interactome was identified. Five genes from this interactome were selected for genetic analysis. Twenty-one tag single-nucleotide polymorphisms (SNPs) which captured all common variation in these genes were genotyped in 10,196 Danes, and analyzed for association with NAFLD-related quantitative traits, type 2 diabetes (T2D), central obesity, and WHO-defined metabolic syndrome (MetS).Results273 genes were included in the protein-protein interaction analysis and EHHADH, ECHS1, HADHA, HADHB, and ACADL were selected for further examination. A total of 10 nominal statistical significant associations (P<0.05) to quantitative metabolic traits were identified. Also, the case-control study showed associations between variation in the five genes and T2D, central obesity, and MetS, respectively. Bonferroni adjustments for multiple testing negated all associations.ConclusionsUsing a bioinformatics approach we identified five candidate genes for NAFLD. However, we failed to provide evidence of associations with major effects between SNPs in these five genes and NAFLD-related quantitative traits, T2D, central obesity, and MetS.

Copy number variation of KIR genes influences HIV-1 control.

A genome-wide screen for large structural variants showed that a copy number variant (CNV) in the region encoding killer cell immunoglobulin-like receptors (KIR) associates with HIV-1 control as measured by plasma viral load at set point in individuals of European ancestry. This CNV encompasses the KIR3DL1-KIR3DS1 locus, encoding receptors that interact with specific HLA-Bw4 molecules to regulate the activation of lymphocyte subsets including natural killer (NK) cells. We quantified the number of copies of KIR3DS1 and KIR3DL1 in a large HIV-1 positive cohort, and showed that an increase in KIR3DS1 count associates with a lower viral set point if its putative ligand is present (p = 0.00028), as does an increase in KIR3DL1 count in the presence of KIR3DS1 and appropriate ligands for both receptors (p = 0.0015). We further provide functional data that demonstrate that NK cells from individuals with multiple copies of KIR3DL1, in the presence of KIR3DS1 and the appropriate ligands, inhibit HIV-1 replication more robustly, and associated with a significant expansion in the frequency of KIR3DS1+, but not KIR3DL1+, NK cells in their peripheral blood. Our results suggest that the relative amounts of these activating and inhibitory KIR play a role in regulating the peripheral expansion of highly antiviral KIR3DS1+ NK cells, which may determine differences in HIV-1 control following infection.

Role of β-TrCP1, variant and homologue in HIV-1 life cycle

Summary The CD4 molecule plays a key role in AIDS pathogenesis, it is required for entry of the virus into permissive cells and its subsequent down-modulation of the cell surface is a hallmark of HN-1 infected cells. The virus encodes no less than three proteins that participate in this process: Nef, Vpu and Env. Vpu protein interacts with CD4 within the endoplasmic reticulum of infected cells, where it targets CD4 for degradation through the interaction with a cellular protein named ß-TrCP1. This F-box protein functions as the substrate recognition subunit of the SCF ß-Trcr E3 ubiquitin ligase, which normally induce the ubiquitination and subsequent degradation of various proteins such as ß-catenin and IxBa. Mammals possess a homologue of ß-TrCP1, HOS, also named ß-TrCP2 which has a cytoplasmic subcellular distribution. Structural analysis of the ligand-binding domain of both homologues shows striking surface similarities. Both F-box proteins have a redundant role in a number of cellular processes; however the potential role of ß-TrCP2 in HIV-1 infected cells has not been evaluated. In the present study, we assessed the existence of génetic variants of BRTC, encoding ß-TrCP1, and evaluated whether these variants would affect CD4 down-modulation. Additionally, we determined whether ß-TrCP2 shares with its homologue structural and functional properties that would allow it to bind Vpu, modulate CD4 expression, and thus participate in HN-1 pathogenesis. We identified a single nucleotide polymorphism present in the human population with an allelic frequency of 0.03 that leads to the substitution of alanine 507 by a serine. However, we showed by transient transfection in HeLa CD4+ cells that this variant behaves as ß-TrCP1 with respect to CD4 down-modulation. We established transient expression systems in HeLa CD4+ cells to test whether ß-TrCP2 is implicated in Vpu-mediated CD4 down-modulation. We show by coimmunoprecipitation experiments that ß-TrCP2 binds Vpu and is able to induce CD4 down-modulation as efficiently as ß-TrCP1. In two different cell lines, HeLa CD4+ and Jurkat, Vpu-mediated CD4 down-modulation could not be completely reversed through the silencing of endogenous ß-TrCP 1 or ß-TrCP2 individually, but required both genes to be silenced simultaneously. We evaluated the role of ß-TrCP1 and ß-TrCP2 in HIV-1 life cycle using silencing prior to actual viral infection. Both ß-TrCP1 and ß-TrCP2 contributed to CD4 down-modulation during aone-cycle viral infection iri Ghost cells. In addition, the combined silencing of both homologues in the absence of env and nef reversed CD4 down-modulation, showing that ß-TrCP 1 and ß-TrCP2 represent the main and additive effectors of HIV-1 encoded Vpu. In addition, we showed that silencing of ß-TrCPI but not ß-TrCP2 induced a decrease of HIV-1 LTR-driven expression. In a transient transfection system with Tat and a LTR luciferase reporter, both homologues modulated LTR-driven expression. The present study revealed that ß-TrCP2 represents a novel protein participating in HIV-1 cycle and complete comprehension of the complex interplay occurring between the two F-Box will improve our understanding of HIV-1 infection. Résumé La molécule CD4 joue un rôle clef dans la pathogenèse du SIDA ; elle est requise pour l'entrée du virus dans les cellules permissives et la diminution de sa concentration au niveau de la surface cellulaire est une importante caractéristique des cellules infectées par le VIH-1. Le virus encode pas moins de trois protéines qui participent à ce processus Nef, Vpu et Env. La protéine Vpu lie CD4 au niveau du réticulum endoplasmique et induit sa dégradation en interagissant avec une protéine cellulaire nommée ß-TrCP 1. Cette protéine de type F-Box est une sous unité du complexe ubiquitine-ligase E3 SCFß-TrCP. Elle permet la reconnaissance du substrat par le complexe qui induit l'ubiquitination et la subséquente dégradation de diverses protéines cellulaires comme la ß-catenin ou IκBα. Les mammifères possèdent un homologue à ß-TrCP1appelé ß-TrCP2 (ou HOS). L'analyse comparative du domaine permettant la reconnaissance des substrats des deux homologues montre de frappantes similarités. Le rôle de ß-TrCP2 dans le cycle viral du VIH-1 n'a pas encore été évalué. Lors de cette étude, nous avons recherché l'existence de variants génétique de BTRC (codant pour ß-TrCP1) et nous avons évalué si ces variants pourraient affecter la dégradation des molécules CD4 induite par le virus. Nous avons ainsi identifié un polymorphisme présent dans la population humaine avec une fréquence allélique de 0.03 qui consiste en une substitution de l'alanine 507 par une sérine. Nous avons cependant montré par transfection dans des cellules HeLa CD4+ que ce variant se comporte comme ß-TrCP 1 en ce qui concerne la modulation de CD4. De plus, nous avons déterminé si ß-TrCP2 partageait avec son homologue des propriétés structurelles et fonctionnelles qui lui permettraient de lier Vpu, moduler la concentration de CD4 et ainsi prendre part à la pathogenèse du SIDA. Pour ce faire, nous avons établi un système d'expression temporaire dans des cellules HeLa CD4+. Par co-immunoprécipitation, nous avons montré que ß-TrCP2 lie Vpu et est capable d'induire la dégradation de CD4 aussi efficacement que ß-TrCP1. Dans deux différentes lignées cellulaires, HeLa CD4+ et Jurkat, la dégradation de CD4 n'a pu être complètement inhibée par le silencing individuel de ß-TrCP 1 ou ß-TrCP2, mais nécessitait le silencing simultané des 2 gènes. Nous avons évalué le rôle des deux homologues dans le cycle viral du VIH-1 en infectant des cellules Ghost avec le virus après avoir effectué un silencing des deux protéines. Nous avons ainsi montré que ß-TrCP 1 et ß-TrCP2 contribuent de manière additive à la dégradation de CD4 induite par une infection du VIH-1. Le silencing combiné des deux homologues inhiba complètement cette dégradation en l'absence de env et nef, prouvant qu'aucune autre voie ne participe à ce processus: En outre, nous avons montré que le silencing de ß-TrCP 1 mais pas celui de ß-TrCP2 induisait une diminution de l'expression virale sous contrôle du LTR. Nous n'avons cependant pas été en mesure de reconstituer cet effet en exprimant Tat et un gène reporteur sous contrôle du LTR dans des cellules HeLa CD4+. Le présent travail révèle que ß-TrCP2 représente une nouvelle protéine participant dans le cycle viral du VIH-1. Une complète compréhension de l'effet de chacun des deux homologues sur le cycle viral permettra d'améliorer notre compréhension de l'infection par le VIH-1.

SCRIB and PUF60 Are Primary Drivers of the Multisystemic Phenotypes of the 8q24.3 Copy-Number Variant.

Copy-number variants (CNVs) represent a significant interpretative challenge, given that each CNV typically affects the dosage of multiple genes. Here we report on five individuals with coloboma, microcephaly, developmental delay, short stature, and craniofacial, cardiac, and renal defects who harbor overlapping microdeletions on 8q24.3. Fine mapping localized a commonly deleted 78 kb region that contains three genes: SCRIB, NRBP2, and PUF60. In vivo dissection of the CNV showed discrete contributions of the planar cell polarity effector SCRIB and the splicing factor PUF60 to the syndromic phenotype, and the combinatorial suppression of both genes exacerbated some, but not all, phenotypic components. Consistent with these findings, we identified an individual with microcephaly, short stature, intellectual disability, and heart defects with a de novo c.505C>T variant leading to a p.His169Tyr change in PUF60. Functional testing of this allele in vivo and in vitro showed that the mutation perturbs the relative dosage of two PUF60 isoforms and, subsequently, the splicing efficiency of downstream PUF60 targets. These data inform the functions of two genes not associated previously with human genetic disease and demonstrate how CNVs can exhibit complex genetic architecture, with the phenotype being the amalgam of both discrete dosage dysfunction of single transcripts and also of binary genetic interactions.

The protective effect of the obesity-associated rs9939609 A variant in fat mass- and obesity-associated gene on depression.

Candidate gene and genome-wide association studies have not identified common variants, which are reliably associated with depression. The recent identification of obesity predisposing genes that are highly expressed in the brain raises the possibility of their genetic contribution to depression. As variation in the intron 1 of the fat mass- and obesity-associated (FTO) gene contributes to polygenic obesity, we assessed the possibility that FTO gene may contribute to depression in a cross-sectional multi-ethnic sample of 6561 depression cases and 21 932 controls selected from the EpiDREAM, INTERHEART, DeCC (depression case-control study) and Cohorte Lausannoise (CoLaus) studies. Major depression was defined according to DSM IV diagnostic criteria. Association analyses were performed under the additive genetic model. A meta-analysis of the four studies showed a significant inverse association between the obesity risk FTO rs9939609 A variant and depression (odds ratio=0.92 (0.89, 0.97), P=3 × 10(-4)) adjusted for age, sex, ethnicity/population structure and body-mass index (BMI) with no significant between-study heterogeneity (I(2)=0%, P=0.63). The FTO rs9939609 A variant was also associated with increased BMI in the four studies (β 0.30 (0.08, 0.51), P=0.0064) adjusted for age, sex and ethnicity/population structure. In conclusion, we provide the first evidence that the FTO rs9939609 A variant may be associated with a lower risk of depression independently of its effect on BMI. This study highlights the potential importance of obesity predisposing genes on depression.

KCTD13 is a major driver of mirrored neuroanatomical phenotypes of the 16p11.2 copy number variant.

Copy number variants (CNVs) are major contributors to genetic disorders. We have dissected a region of the 16p11.2 chromosome--which encompasses 29 genes--that confers susceptibility to neurocognitive defects when deleted or duplicated. Overexpression of each human transcript in zebrafish embryos identified KCTD13 as the sole message capable of inducing the microcephaly phenotype associated with the 16p11.2 duplication, whereas suppression of the same locus yielded the macrocephalic phenotype associated with the 16p11.2 deletion, capturing the mirror phenotypes of humans. Analyses of zebrafish and mouse embryos suggest that microcephaly is caused by decreased proliferation of neuronal progenitors with concomitant increase in apoptosis in the developing brain, whereas macrocephaly arises by increased proliferation and no changes in apoptosis. A role for KCTD13 dosage changes is consistent with autism in both a recently reported family with a reduced 16p11.2 deletion and a subject reported here with a complex 16p11.2 rearrangement involving de novo structural alteration of KCTD13. Our data suggest that KCTD13 is a major driver for the neurodevelopmental phenotypes associated with the 16p11.2 CNV, reinforce the idea that one or a small number of transcripts within a CNV can underpin clinical phenotypes, and offer an efficient route to identifying dosage-sensitive loci.

Orosomucoid (alpha-1 acid glycoprotein) phenotyping by use of immobilized pH gradients with 8 M urea and immunoblotting. A new variant encountered in a population study.

Orosomucoid (ORM) phenotyping has been performed on 329 unrelated Swiss subjects, using immobilized pH gradients with 8 M urea and 2% v/v 2-mercaptoethanol followed by immunoblotting. After desialylation the band patterns of ORM confirmed that the polymorphism of the structural locus ORM1 is controlled by three codominant autosomal alleles (ORM1*F1, ORM1*S and ORM1*F2). One rare and one new allele were detected. The rare variant, tentatively assigned to the second structural locus ORM2, is observed in a cathodal position and named ORM2 B1. The new variant, tentatively assigned to the first structural locus ORM1, is observed in a region located between ORM1 S and ORM1 F2, and named ORM1 F3. Moreover, the pI values of the ORM variants have been measured accurately with Immobiline Dry Plates (LKB): they were found to be within the pH range 4.93-5.14.