Deleterious GRM1 mutations in schizophrenia

We analysed a phenotypically well-characterised sample of 450 schziophrenia patients and 605 controls for rare non-synonymous single nucleotide polymorphisms (nsSNPs) in the GRM1 gene, their functional effects and family segregation. GRM1 encodes the metabotropic glutamate receptor 1 (mGluR1), whose documented role as a modulator of neuronal signalling and synaptic plasticity makes it a plausible schizophrenia candidate. In a recent study, this gene was shown to harbour a cluster of deleterious nsSNPs within a functionally important domain of the receptor, in patients with schizophrenia and bipolar disorder. Our Sanger sequencing of the GRM1 coding regions detected equal numbers of nsSNPs in cases and controls, however the two groups differed in terms of the potential effects of the variants on receptor function: 6/6 case-specific and only 1/6 control-specific nsSNPs were predicted to be deleterious. Our in-vitro experimental follow-up of the case-specific mutants showed that 4/6 led to significantly reduced inositol phosphate production, indicating impaired function of the major mGluR1 signalling pathway; 1/6 had reduced cell membrane expression; inconclusive results were obtained in 1/6. Family segregation analysis indicated that these deleterious nsSNPs were inherited. Interestingly, four of the families were affected by multiple neuropsychiatric conditions, not limited to schizophrenia, and the mutations were detected in relatives with schizophrenia, depression and anxiety, drug and alcohol dependence, and epilepsy. Our findings suggest a possible mGluR1 contribution to diverse psychiatric conditions, supporting the modulatory role of the receptor in such conditions as proposed previously on the basis of in vitro experiments and animal studies.     

Two newly identified exons in human GRM1 express a novel splice variant of metabotropic glutamate 1 receptor

To date, five human metabotropic glutamate (mGlu) 1 receptor splice variants (1a, 1b, 1d, 1f, and 1g) have been described, all of which involve alternative C-terminal splicing. mGlu1a receptor contains a long C-terminal domain (341 amino acids), which has been shown to scaffold with several proteins and contribute to the structure of the post-synaptic density. However, several shorter mGlu1 receptor splice variants lack the sequence required for these interactions, and no major functional differences between these short splice variants have been described. By using RT-PCR we have shown that two human melanoma cell lines express both mGlu1a and mGlu1b receptors. In addition, using 3′RACE, we identified three previously unknown mGlu1 receptor mRNAs. Two differ in the length of their 3′ untranslated region (UTR), and encode the same predicted protein as mGlu1g receptor—the shortest of all mGlu1 receptor splice variants. The third mRNA, named mGlu1h, encodes a predicted C-terminal splice variant of 10 additional amino acids. mGlu1h mRNA was observed in two different melanoma cell lines and is overexpressed, compared with melanoma precursor cells, melanocytes. Most importantly, this new splice variant, mGlu1h receptor, is encoded by two previously unidentified exons located within the human GRM1 gene. Additionally, these new exons are found exclusively within the GRM1 genes of higher primates and are highly conserved. Therefore, we hypothesize that mGlu1h receptors play a distinct role in primate glutamatergic signaling.     

A Cytogenetic Abnormality and Rare Coding Variants Identify ABCA13 as a Candidate Gene in Schizophrenia, Bipolar Disorder, and Depression

Schizophrenia and bipolar disorder are leading causes of morbidity across all populations, with heritability estimates of ∼80% indicating a substantial genetic component. Population genetics and genome-wide association studies suggest an overlap of genetic risk factors between these illnesses but it is unclear how this genetic component is divided between common gene polymorphisms, rare genomic copy number variants, and rare gene sequence mutations. We report evidence that the lipid transporter gene ABCA13 is a susceptibility factor for both schizophrenia and bipolar disorder. After the initial discovery of its disruption by a chromosome abnormality in a person with schizophrenia, we resequenced ABCA13 exons in 100 cases with schizophrenia and 100 controls. Multiple rare coding variants were identified including one nonsense and nine missense mutations and compound heterozygosity/homozygosity in six cases. Variants were genotyped in additional schizophrenia, bipolar, depression (n > 1600), and control (n > 950) cohorts and the frequency of all rare variants combined was greater than controls in schizophrenia (OR = 1.93, p = 0.0057) and bipolar disorder (OR = 2.71, p = 0.00007). The population attributable risk of these mutations was 2.2% for schizophrenia and 4.0% for bipolar disorder. In a study of 21 families of mutation carriers, we genotyped affected and unaffected relatives and found significant linkage (LOD = 4.3) of rare variants with a phenotype including schizophrenia, bipolar disorder, and major depression. These data identify a candidate gene, highlight the genetic overlap between schizophrenia, bipolar disorder, and depression, and suggest that rare coding variants may contribute significantly to risk of these disorders.     

Bipolar I disorder and schizophrenia: a 440-single-nucleotide polymorphism screen of 64 candidate genes among Ashkenazi Jewish case-parent trios

Bipolar, schizophrenia, and schizoaffective disorders are common, highly heritable psychiatric disorders, for which familial coaggregation, as well as epidemiological and genetic evidence, suggests overlapping etiologies. No definitive susceptibility genes have yet been identified for any of these disorders. Genetic heterogeneity, combined with phenotypic imprecision and poor marker coverage, has contributed to the difficulty in defining risk variants. We focused on families of Ashkenazi Jewish descent, to reduce genetic heterogeneity, and, as a precursor to genomewide association studies, we undertook a single-nucleotide polymorphism (SNP) genotyping screen of 64 candidate genes (440 SNPs) chosen on the basis of previous linkage or of association and/or biological relevance. We genotyped an average of 6.9 SNPs per gene, with an average density of 1 SNP per 11.9 kb in 323 bipolar I disorder and 274 schizophrenia or schizoaffective Ashkenazi case-parent trios. Using single-SNP and haplotype-based transmission/disequilibrium tests, we ranked genes on the basis of strength of association (P<.01). Six genes (DAO, GRM3, GRM4, GRIN2B, IL2RB, and TUBA8) met this criterion for bipolar I disorder; only DAO has been previously associated with bipolar disorder. Six genes (RGS4, SCA1, GRM4, DPYSL2, NOS1, and GRID1) met this criterion for schizophrenia or schizoaffective disorder; five replicate previous associations, and one, GRID1, shows a novel association with schizophrenia. In addition, six genes (DPYSL2, DTNBP1, G30/G72, GRID1, GRM4, and NOS1) showed overlapping suggestive evidence of association in both disorders. These results may help to prioritize candidate genes for future study from among the many suspected/proposed for schizophrenia and bipolar disorders. They provide further support for shared genetic susceptibility between these two disorders that involve glutamate-signaling pathways.     

Melanoma‐associated GRM3 variants dysregulate melanosome trafficking and cAMP signaling

Large‐scale sequencing studies have revealed several genes that are recurrently mutated in melanomas. To annotate the melanoma genome, we have expressed tumor‐associated variants of these genes in zebrafish and characterized their effects on melanocyte development and function. Here, we describe expression of tumor‐associated variants of the recurrently mutated metabotropic glutamate receptor 3 (GRM3) gene. Unlike wild‐type GRM3, tumor‐associated GRM3 variants disrupted trafficking of melanosomes, causing their aggregation in the cell body. Melanosomes are trafficked in a cAMP‐dependent manner, and drugs that directly or indirectly increased cAMP levels were able to suppress melanosome aggregation in mutant GRM3‐expressing melanocytes. Our data show that oncogenic GRM3 variants dysregulate cAMP signaling, a heretofore unknown role for these oncogenes. cAMP signaling has been implicated in melanoma progression and drug resistance, and our data show that oncogenic properties of GRM3 could be mediated, at least in part, by alterations in cAMP signaling.     

Cloning and characterization of a human orphan family C G-protein coupled receptor GPRC5D

Recently three orphan G-protein coupled receptors, RAIG1, GPRC5B and GPRC5C, with homology to members of family C (metabotropic glutamate receptor-like) have been identified. Using the protein sequences of these receptors as queries we identified overlapping expressed sequence tags which were predicted to encode an additional subtype. The full length coding regions of mouse mGprc5d and human GPRC5D were cloned and shown to contain predicted open reading frames of 300 and 345 amino acids, respectively. GPRC5D has seven putative transmembrane segments and is expressed in the cell membrane. The four human receptor subtypes, which we assign to group 5 of family C GPCRs, show 31-42% amino acid sequence identity to each other and 20-25% sequence identity to the transmembrane domains of metabotropic glutamate receptor subtypes 2 and 3 and other family C members. In contrast to the remaining family C members, the group 5 receptors have short amino terminal domains of some 30-50 amino acids. GPRC5D was shown to be clustered with RAIG1 on chromosome 12p13.3 and like RAIG1 and GPRC5B to consist of three exons, the first exon being the largest containing all seven transmembrane segments. GPRC5D mRNA is widely expressed in the peripheral system but all four receptors show distinct expression patterns. Interestingly, mRNA levels of all four group 5 receptors were found in medium to high levels in the kidney, pancreas and prostate and in low to medium levels in the colon and the small intestine, whereas other organs only express a subset of the genes. In an attempt to delineate the signal transduction pathway(s) of the orphan receptors, a series of chimeric receptors containing the amino terminal domain of the calcium sensing receptor or metabotropic glutamate receptor subtype 1, and the seven transmembrane domain of the orphan receptors were constructed and tested in binding and functional assays.     

Naturally occurring CCR5 extracellular and transmembrane domain variants affect HIV-1 Co-receptor and ligand binding function.

Analysis of CCR5 variants in human immunodeficiency virus, type 1 (HIV-1), high risk cohorts led to the identification of multiple single amino acid substitutions in the amino-terminal third of the HIV-1 co-receptor CCR5 suggesting the possibility of protective and permissive genotypes; unfortunately, the low frequency of these mutations did not led to correlation with function. Therefore, we used analytical methods to assess the functional and structural significance of six of these variant receptors in vitro. These studies showed three categories of effects on CCR5 function. 1) Mutations in the first extracellular domain of CCR5 severely reduce specific ligand binding and chemokine-induced chemotaxis. 2) An extracellular domain variant, A29S, when co-expressed with CD4, supported HIV-1 infection whereas the others do not. 3) The transmembrane region variants of CCR5 support monotropic HIV-1 infection that is blocked by addition of some receptor agonists. Mutations in the first and second transmembrane domains increase RANTES (regulated on activation normal T-cell expressed) binding affinity but did not affect MIP1beta binding affinity presumably based on differences in ligand-receptor interaction sites. Furthermore, the CCR5 transmembrane mutants do not respond to RANTES with the classical bell-shaped chemotactic response curve, suggesting that they are resistant to RANTES-induced desensitization. These data demonstrate that single amino acid changes in the extracellular domains of CCR5 can have profound effects on both HIV-1 co-receptor and specific ligand-induced functions, whereas mutations in the transmembrane domain only affect the response to chemokine ligands.     

Exome Sequencing of 75 Individuals from Multiply Affected Coeliac Families and Large Scale Resequencing Follow Up

Coeliac disease (CeD) is a highly heritable common autoimmune disease involving chronic small intestinal inflammation in response to dietary wheat. The human leukocyte antigen (HLA) region, and 40 newer regions identified by genome wide association studies (GWAS) and dense fine mapping, account for ∼40% of the disease heritability. We hypothesized that in pedigrees with multiple individuals with CeD rare [minor allele frequency (MAF) <0.5%] mutations of larger effect size (odds ratios of ∼ 2–5) might exist. We sequenced the exomes of 75 coeliac individuals of European ancestry from 55 multiply affected families. We selected interesting variants and genes for further follow up using a combination of: an assessment of shared variants between related subjects, a model-free linkage test, and gene burden tests for multiple, potentially causal, variants. We next performed highly multiplexed amplicon resequencing of all RefSeq exons from 24 candidate genes selected on the basis of the exome sequencing data in 2,248 unrelated coeliac cases and 2,230 controls. 1,335 variants with a 99.9% genotyping call rate were observed in 4,478 samples, of which 939 were present in coding regions of 24 genes (Ti/Tv 2.99). 91.7% of coding variants were rare (MAF <0.5%) and 60% were novel. Gene burden tests performed on rare functional variants identified no significant associations (p<1×10⁻³) in the resequenced candidate genes. Our strategy of sequencing multiply affected families with deep follow up of candidate genes has not identified any new CeD risk mutations.     

Update on molecular diagnosis of hereditary hemorrhagic telangiectasia

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant, vascular disease hallmarked by the development of arteriovenous malformations (AVMs). Germline mutations in two genes, endoglin (ENG) and activin receptor like kinase 1 (ACVRL1), have been implicated in this disease. This report describes molecular diagnosis in a consecutive series of 600 individuals with clinical features of HHT disease. Each coding exon and flanking intronic regions of ENG and ACVRL1 genes was sequenced. Exonic copy number was quantified in probands without a coding sequence mutation. Novel nonsynonymous variants were further analyzed to predict functional consequences. In addition, common single nucleotide polymorphisms genotypes and haplotypes for the two genes were compared between individuals with and without mutations. The highest mutation detection rate (87% [95% CI 80.2–91.5]) was observed in probands who met all four Curacao criteria (epistaxis, telangiectases, AVMs and family history). More than 30% of identified mutations were novel; however, only 6% were variants of unknown significance. Determining the significance of novel mutations as related to disease presents additional challenges. Detection of multiple alterations in the same proband also requires careful evaluation for disease association. In conclusion, the sensitivity of molecular diagnosis is highest in probands with a clinically confirmed diagnosis of HHT. However, a substantial fraction of probands in this group do not carry an identifiable mutation in the coding exons of these two genes. This suggests alternate mechanisms of gene inactivation or involvement of alternate loci, and it requires further investigation.     

Piscine glycoprotein hormone (gonadotropin and thyrotropin) receptors: a review of recent developments

Similar to the higher vertebrates, the pituitary in bony fishes express three glycoprotein hormones: thyroid-stimulating hormone (TSH), follicle-stimulating hormone (FSH) and luteinizing hormone (LH). In addition to the appropriate secretion of these hormones, the timely and quantitative expression of their specific receptors (TSHR, FSHR and LHR) in the target tissues is an essential requirement for their physiological action. In fishes that constitute a very diverse group of vertebrates, there are only a few published reports of primary structure of these receptors although other examples have been communicated briefly. This review will summarize these reports as well as to describe the insights gained from what is known about the mammalian receptors. The structural organization of the fish receptors (as deduced from the encoding cDNAs) is highly homologous to the higher vertebrate receptors in that there is a 7-pass transmembrane region and an N-terminal extracellular domain, which contributes to ligand specificity. In mammals, the FSHR and the TSHR genes are composed of 10 exons whereas the LHR gene is composed of 11 exons. The position of the 'extra intron' is conserved in the catfish LHR gene. In the mammals, the transmembrane domain of each of the three glycoprotein hormone receptors is encoded by a single exon, however, in the salmon genes and homologous invertebrate genes, this portion of the receptor is encoded by multiple exons. In general, the tissue-specific expression of these receptors is similar to that seen in mammals, however, the gonadal expression of TSHR in the striped bass and sunrise sculpin and the renal expression of LHR in the channel catfish are unique.     

Mitochondrial Mutations in Subjects with Psychiatric Disorders

A considerable body of evidence supports the role of mitochondrial dysfunction in psychiatric disorders and mitochondrial DNA (mtDNA) mutations are known to alter brain energy metabolism, neurotransmission, and cause neurodegenerative disorders. Genetic studies focusing on common nuclear genome variants associated with these disorders have produced genome wide significant results but those studies have not directly studied mtDNA variants. The purpose of this study is to investigate, using next generation sequencing, the involvement of mtDNA variation in bipolar disorder, schizophrenia, major depressive disorder, and methamphetamine use. MtDNA extracted from multiple brain regions and blood were sequenced (121 mtDNA samples with an average of 8,800x coverage) and compared to an electronic database containing 26,850 mtDNA genomes. We confirmed novel and rare variants, and confirmed next generation sequencing error hotspots by traditional sequencing and genotyping methods. We observed a significant increase of non-synonymous mutations found in individuals with schizophrenia. Novel and rare non-synonymous mutations were found in psychiatric cases in mtDNA genes: ND6, ATP6, CYTB, and ND2. We also observed mtDNA heteroplasmy in brain at a locus previously associated with schizophrenia (T16519C). Large differences in heteroplasmy levels across brain regions within subjects suggest that somatic mutations accumulate differentially in brain regions. Finally, multiplasmy, a heteroplasmic measure of repeat length, was observed in brain from selective cases at a higher frequency than controls. These results offer support for increased rates of mtDNA substitutions in schizophrenia shown in our prior results. The variable levels of heteroplasmic/multiplasmic somatic mutations that occur in brain may be indicators of genetic instability in mtDNA.     

Genome Distribution of Replication-independent Histone H1 Variants Shows H1.0 Associated with Nucleolar Domains and H1X Associated with RNA Polymerase II-enriched Regions

Unlike core histones, the linker histone H1 family is more evolutionarily diverse, and many organisms have multiple H1 variants or subtypes. In mammals, the H1 family includes seven somatic H1 variants; H1.1 to H1.5 are expressed in a replication-dependent manner, whereas H1.0 and H1X are replication-independent. Using ChIP-sequencing data and cell fractionation, we have compared the genomic distribution of H1.0 and H1X in human breast cancer cells, in which we previously observed differential distribution of H1.2 compared with the other subtypes. We have found H1.0 to be enriched at nucleolus-associated DNA repeats and chromatin domains, whereas H1X is associated with coding regions, RNA polymerase II-enriched regions, and hypomethylated CpG islands. Further, H1X accumulates within constitutive or included exons and retained introns and toward the 3′ end of expressed genes. Inducible H1X knockdown does not affect cell proliferation but dysregulates a subset of genes related to cell movement and transport. In H1X-depleted cells, the promoters of up-regulated genes are not occupied specifically by this variant, have a lower than average H1 content, and, unexpectedly, do not form an H1 valley upon induction. We conclude that H1 variants are not distributed evenly across the genome and may participate with some specificity in chromatin domain organization or gene regulation.     

A Hybrid Likelihood Model for Sequence-Based Disease Association Studies

In the past few years, case-control studies of common diseases have shifted their focus from single genes to whole exomes. New sequencing technologies now routinely detect hundreds of thousands of sequence variants in a single study, many of which are rare or even novel. The limitation of classical single-marker association analysis for rare variants has been a challenge in such studies. A new generation of statistical methods for case-control association studies has been developed to meet this challenge. A common approach to association analysis of rare variants is the burden-style collapsing methods to combine rare variant data within individuals across or within genes. Here, we propose a new hybrid likelihood model that combines a burden test with a test of the position distribution of variants. In extensive simulations and on empirical data from the Dallas Heart Study, the new model demonstrates consistently good power, in particular when applied to a gene set (e.g., multiple candidate genes with shared biological function or pathway), when rare variants cluster in key functional regions of a gene, and when protective variants are present. When applied to data from an ongoing sequencing study of bipolar disorder (191 cases, 107 controls), the model identifies seven gene sets with nominal p-values0.05, of which one MAPK signaling pathway (KEGG) reaches trend-level significance after correcting for multiple testing.     

TRPM1 Is Mutated in Patients with Autosomal-Recessive Complete Congenital Stationary Night Blindness

Night vision requires signaling from rod photoreceptors to adjacent bipolar cells in the retina. Mutations in the genes NYX and GRM6, expressed in ON bipolar cells, lead to a disruption of the ON bipolar cell response. This dysfunction is present in patients with complete X-linked and autosomal-recessive congenital stationary night blindness (CSNB) and can be assessed by standard full-field electroretinography (ERG), showing severely reduced rod b-wave amplitude and slightly altered cone responses. Although many cases of complete CSNB (cCSNB) are caused by mutations in NYX and GRM6, in ∼60% of the patients the gene defect remains unknown. Animal models of human diseases are a good source for candidate genes, and we noted that a cCSNB phenotype present in homozygous Appaloosa horses is associated with downregulation of TRPM1. TRPM1, belonging to the family of transient receptor potential channels, is expressed in ON bipolar cells and therefore qualifies as an excellent candidate. Indeed, mutation analysis of 38 patients with CSNB identified ten unrelated cCSNB patients with 14 different mutations in this gene. The mutation spectrum comprises missense, splice-site, deletion, and nonsense mutations. We propose that the cCSNB phenotype in these patients is due to the absence of functional TRPM1 in retinal ON bipolar cells.     

Identification of Novel GRM1 Mutations and Single Nucleotide Polymorphisms in Prostate Cancer Cell Lines and Tissues

Metabotropic glutamate receptor 1 (GRM1) signaling has been implicated in benign and malignant disorders including prostate cancer (PCa). To further explore the role of genetic alterations of GRM1 in PCa, we screened the entire human GRM1 gene including coding sequence, exon-intron junctions, and flanking untranslated regions (UTRs) for the presence of mutations and single nucleotide polymorphisms (SNPs) in several PCa cell lines and matched tumor-normal tissues from Caucasian Americans (CAs) and African Americans (AAs). We used bidirectional sequencing, allele-specific PCR, and bioinformatics to identify the genetic changes in GRM1 and to predict their functional role. A novel missense mutation identified at C1744T (582 Pro > Ser) position of GRM1 gene in a primary AA-PCa cell line (E006AA) was predicted to affect the protein stability and functions. Another novel mutation identified at exon-intron junction of exon-8 in C4-2B cell line resulted in alteration of the GRM1 splicing donor site. In addition, we found missense SNP at T2977C (993 Ser > Pro) position and multiple non-coding mutations and SNPs in 3′-UTR of GRM1 gene in PCa cell lines and tissues. These novel mutations may contribute to the disease by alterations in GRM1 gene splicing, receptor activation, and post-receptor downstream signaling.