APOE Molecular Spectrum in a French Cohort with Primary Dyslipidemia.

Primary hypercholesterolemia is characterized by elevated LDL-cholesterol (LDL-C) levels isolated in autosomal dominant hypercholesterolemia (ADH) or associated with elevated triglyceride levels in familial combined hyperlipidemia (FCHL). Rare APOE variants are known in ADH and FCHL. We explored the APOE molecular spectrum in a French ADH/FCHL cohort of 5743 unrelated probands. The sequencing of LDLR, PCSK9, APOB, and APOE revealed 76 carriers of a rare APOE variant, with no mutation in LDLR, PCSK9, or APOB. Among the 31 APOE variants identified here, 15 are described in ADH, 10 in FCHL, and 6 in both probands. Five were previously reported with dyslipidemia and 26 are novel, including 12 missense, 5 synonymous, 2 intronic, and 7 variants in regulatory regions. Sixteen variants were predicted as pathogenic or likely pathogenic, and their carriers had significantly lower polygenic risk scores (wPRS) than carriers of predicted benign variants. We observed no correlation between LDL-C levels and wPRS, suggesting a major effect of APOE variants. Carriers of p.Leu167del were associated with a severe phenotype. The analysis of 11 probands suggests that carriers of an APOE variant respond better to statins than carriers of a LDLR mutation. Altogether, we show that the APOE variants account for a significant contribution to ADH and FCHL.

Whole Exome/Genome Sequencing Joint Analysis of a Family with Oligogenic Familial Hypercholesterolemia.

Autosomal Dominant Hypercholesterolemia (ADH) is a genetic disorder caused by pathogenic variants in LDLR, APOB, PCSK9 and APOE genes. We sought to identify new candidate genes responsible for the ADH phenotype in patients without pathogenic variants in the known ADH-causing genes by focusing on a French family with affected and non-affected members who presented a high ADH polygenic risk score (wPRS). Linkage analysis, whole exome and whole genome sequencing resulted in the identification of variants p.(Pro398Ala) in CYP7A1, p.(Val1382Phe) in LRP6 and p.(Ser202His) in LDLRAP1. A total of 6 other variants were identified in 6 of 160 unrelated ADH probands: p.(Ala13Val) and p.(Aps347Asn) in CYP7A1; p.(Tyr972Cys), p.(Thr1479Ile) and p.(Ser1612Phe) in LRP6; and p.(Ser202LeufsTer19) in LDLRAP1. All six probands presented a moderate wPRS. Serum analyses of carriers of the p.(Pro398Ala) variant in CYP7A1 showed no differences in the synthesis of bile acids compared to the serums of non-carriers. Functional studies of the four LRP6 mutants in HEK293T cells resulted in contradictory results excluding a major effect of each variant alone. Within the family, none of the heterozygous for only the LDLRAP1 p.(Ser202His) variant presented ADH. Altogether, each variant individually does not result in elevated LDL-C; however, the oligogenic combination of two or three variants reveals the ADH phenotype.

Posttranscriptional Regulation of the Human LDL Receptor by the U2-Spliceosome.

BACKGROUND: The LDLR (low-density lipoprotein receptor) in the liver is the major determinant of LDL-cholesterol levels in human plasma. The discovery of genes that regulate the activity of LDLR helps to identify pathomechanisms of hypercholesterolemia and novel therapeutic targets against atherosclerotic cardiovascular disease. METHODS: We performed a genome-wide RNA interference screen for genes limiting the uptake of fluorescent LDL into Huh-7 hepatocarcinoma cells. Top hit genes were validated by in vitro experiments as well as analyses of data sets on gene expression and variants in human populations. RESULTS: The knockdown of 54 genes significantly inhibited LDL uptake. Fifteen of them encode for components or interactors of the U2-spliceosome. Knocking down any one of 11 out of 15 genes resulted in the selective retention of intron 3 of LDLR. The translated LDLR fragment lacks 88% of the full length LDLR and is detectable neither in nontransfected cells nor in human plasma. The hepatic expression of the intron 3 retention transcript is increased in nonalcoholic fatty liver disease as well as after bariatric surgery. Its expression in blood cells correlates with LDL-cholesterol and age. Single nucleotide polymorphisms and 3 rare variants of one spliceosome gene, RBM25, are associated with LDL-cholesterol in the population and familial hypercholesterolemia, respectively. Compared with overexpression of wild-type RBM25, overexpression of the 3 rare RBM25 mutants in Huh-7 cells led to lower LDL uptake. CONCLUSIONS: We identified a novel mechanism of posttranscriptional regulation of LDLR activity in humans and associations of genetic variants of RBM25 with LDL-cholesterol levels.

Mutations and variants of ONECUT1 in diabetes.

Genes involved in distinct diabetes types suggest shared disease mechanisms. Here we show that One Cut Homeobox 1 (ONECUT1) mutations cause monogenic recessive syndromic diabetes in two unrelated patients, characterized by intrauterine growth retardation, pancreas hypoplasia and gallbladder agenesis/hypoplasia, and early-onset diabetes in heterozygous relatives. Heterozygous carriers of rare coding variants of ONECUT1 define a distinctive subgroup of diabetic patients with early-onset, nonautoimmune diabetes, who respond well to diabetes treatment. In addition, common regulatory ONECUT1 variants are associated with multifactorial type 2 diabetes. Directed differentiation of human pluripotent stem cells revealed that loss of ONECUT1 impairs pancreatic progenitor formation and a subsequent endocrine program. Loss of ONECUT1 altered transcription factor binding and enhancer activity and NKX2.2/NKX6.1 expression in pancreatic progenitor cells. Collectively, we demonstrate that ONECUT1 controls a transcriptional and epigenetic machinery regulating endocrine development, involved in a spectrum of diabetes, encompassing monogenic (recessive and dominant) as well as multifactorial inheritance. Our findings highlight the broad contribution of ONECUT1 in diabetes pathogenesis, marking an important step toward precision diabetes medicine.

DNAJC3 deficiency induces ß-cell mitochondrial apoptosis and causes syndromic young-onset diabetes.

OBJECTIVE: DNAJC3, also known as P58IPK, is an Hsp40 family member that interacts with and inhibits PKR-like ER-localized eIF2α kinase (PERK). Dnajc3 deficiency in mice causes pancreatic ß-cell loss and diabetes. Loss-of-function mutations in DNAJC3 cause early-onset diabetes and multisystemic neurodegeneration. The aim of our study was to investigate the genetic cause of early-onset syndromic diabetes in two unrelated patients, and elucidate the mechanisms of ß-cell failure in this syndrome. METHODS: Whole exome sequencing was performed and identified variants were confirmed by Sanger sequencing. DNAJC3 was silenced by RNAi in INS-1E cells, primary rat ß-cells, human islets, and induced pluripotent stem cell-derived ß-cells. ß-cell function and apoptosis were assessed, and potential mediators of apoptosis examined. RESULTS: The two patients presented with juvenile-onset diabetes, short stature, hypothyroidism, neurodegeneration, facial dysmorphism, hypoacusis, microcephaly and skeletal bone deformities. They were heterozygous compound and homozygous for novel loss-of-function mutations in DNAJC3. DNAJC3 silencing did not impair insulin content or secretion. Instead, the knockdown induced rat and human ß-cell apoptosis and further sensitized cells to endoplasmic reticulum stress, triggering mitochondrial apoptosis via the pro-apoptototic Bcl-2 proteins BIM and PUMA. CONCLUSIONS: This report confirms previously described features and expands the clinical spectrum of syndromic DNAJC3 diabetes, one of the five monogenic forms of diabetes pertaining to the PERK pathway of the endoplasmic reticulum stress response. DNAJC3 deficiency may lead to ß-cell loss through BIM- and PUMA-dependent activation of the mitochondrial pathway of apoptosis.

Juvenile-Onset Diabetes and Congenital Cataract: "Double-Gene" Mutations Mimicking a Syndromic Diabetes Presentation.

Monogenic forms of diabetes may account for 1-5% of all cases of diabetes, and may occur in the context of syndromic presentations. We investigated the case of a girl affected by insulin-dependent diabetes, diagnosed at 6 years old, associated with congenital cataract. Her consanguineous parents and her four other siblings did not have diabetes or cataract, suggesting a recessive syndrome. Using whole exome sequencing of the affected proband, we identified a heterozygous p.R825Q ABCC8 mutation, located at the exact same amino-acid position as the p.R825W recurring diabetes mutation, hence likely responsible for the diabetes condition, and a homozygous p.G71S mutation in CRYBB1, a gene known to be responsible for congenital cataract. Both mutations were predicted to be damaging and were absent or extremely rare in public databases. Unexpectedly, we found that the mother was also homozygous for the CRYBB1 mutation, and both the mother and one unaffected sibling were heterozygous for the ABCC8 mutation, suggesting incomplete penetrance of both mutations. Incomplete penetrance of ABCC8 mutations is well documented, but this is the first report of an incomplete penetrance of a CRYBB1 mutation, manifesting between susceptible subjects (unaffected mother vs. affected child) and to some extent within the patient herself, who had distinct cataract severities in both eyes. Our finding illustrates the importance of family studies to unmask the role of confounding factors such as double-gene mutations and incomplete penetrance that may mimic monogenic syndromes including in the case of strongly evocative family structure with consanguinity.

dUTPase (DUT) Is Mutated in a Novel Monogenic Syndrome With Diabetes and Bone Marrow Failure.

We describe a new syndrome characterized by early-onset diabetes associated with bone marrow failure, affecting mostly the erythrocytic lineage. Using whole-exome sequencing in a remotely consanguineous patient from a family with two affected siblings, we identified a single homozygous missense mutation (chr15.hg19:g.48,626,619A>G) located in the dUTPase (DUT) gene (National Center for Biotechnology Information Gene ID 1854), affecting both the mitochondrial (DUT-M p.Y142C) and the nuclear (DUT-N p.Y54C) isoforms. We found the same homozygous mutation in an unrelated consanguineous patient with diabetes and bone marrow aplasia from a family with two affected siblings, whereas none of the >60,000 subjects from the Exome Aggregation Consortium (ExAC) was homozygous for this mutation. This replicated observation probability was highly significant, thus confirming the role of this DUT mutation in this syndrome. DUT is a key enzyme for maintaining DNA integrity by preventing misincorporation of uracil into DNA, which results in DNA toxicity and cell death. We showed that DUT silencing in human and rat pancreatic ß-cells results in apoptosis via the intrinsic cell death pathway. Our findings support the importance of tight control of DNA metabolism for ß-cell integrity and warrant close metabolic monitoring of patients treated by drugs affecting dUTP balance.

A Missense Mutation in PPP1R15B Causes a Syndrome Including Diabetes, Short Stature, and Microcephaly.

Dysregulated endoplasmic reticulum stress and phosphorylation of eukaryotic translation initiation factor 2α (eIF2α) are associated with pancreatic ß-cell failure and diabetes. Here, we report the first homozygous mutation in the PPP1R15B gene (also known as constitutive repressor of eIF2α phosphorylation [CReP]) encoding the regulatory subunit of an eIF2α-specific phosphatase in two siblings affected by a novel syndrome of diabetes of youth with short stature, intellectual disability, and microcephaly. The R658C mutation in PPP1R15B affects a conserved amino acid within the domain important for protein phosphatase 1 (PP1) binding. The R658C mutation decreases PP1 binding and eIF2α dephosphorylation and results in ß-cell apoptosis. Our findings support the concept that dysregulated eIF2α phosphorylation, whether decreased by mutation of the kinase (EIF2AK3) in Wolcott-Rallison syndrome or increased by mutation of the phosphatase (PPP1R15B), is deleterious to ß-cells and other secretory tissues, resulting in diabetes associated with multisystem abnormalities.

Heterozygous HTRA1 mutations are associated with autosomal dominant cerebral small vessel disease.

Cerebral small vessel disease represents a heterogeneous group of disorders leading to stroke and cognitive impairment. While most small vessel diseases appear sporadic and related to age and hypertension, several early-onset monogenic forms have also been reported. However, only a minority of patients with familial small vessel disease carry mutations in one of known small vessel disease genes. We used whole exome sequencing to identify candidate genes in an autosomal dominant small vessel disease family in which known small vessel disease genes had been excluded, and subsequently screened all candidate genes in 201 unrelated probands with a familial small vessel disease of unknown aetiology, using high throughput multiplex polymerase chain reaction and next generation sequencing. A heterozygous HTRA1 variant (R166L), absent from 1000 Genomes and Exome Variant Server databases and predicted to be deleterious by in silico tools, was identified in all affected members of the index family. Ten probands of 201 additional unrelated and affected probands (4.97%) harboured a heterozygous HTRA1 mutation predicted to be damaging. There was a highly significant difference in the number of likely deleterious variants in cases compared to controls (P = 4.2 × 10(-6); odds ratio = 15.4; 95% confidence interval = 4.9-45.5), strongly suggesting causality. Seven of these variants were located within or close to the HTRA1 protease domain, three were in the N-terminal domain of unknown function and one in the C-terminal PDZ domain. In vitro activity analysis of HTRA1 mutants demonstrated a loss of function effect. Clinical features of this autosomal dominant small vessel disease differ from those of CARASIL and CADASIL by a later age of onset and the absence of the typical extraneurological features of CARASIL. They are similar to those of sporadic small vessel disease, except for their familial nature. Our data demonstrate that heterozygous HTRA1 mutations are an important cause of familial small vessel disease, and that screening of HTRA1 should be considered in all patients with a hereditary small vessel disease of unknown aetiology.

Loss of α1ß1 soluble guanylate cyclase, the major nitric oxide receptor, leads to moyamoya and achalasia.

Moyamoya is a cerebrovascular condition characterized by a progressive stenosis of the terminal part of the internal carotid arteries (ICAs) and the compensatory development of abnormal "moyamoya" vessels. The pathophysiological mechanisms of this condition, which leads to ischemic and hemorrhagic stroke, remain unknown. It can occur as an isolated cerebral angiopathy (so-called moyamoya disease) or in association with various conditions (moyamoya syndromes). Here, we describe an autosomal-recessive disease leading to severe moyamoya and early-onset achalasia in three unrelated families. This syndrome is associated in all three families with homozygous mutations in GUCY1A3, which encodes the α1 subunit of soluble guanylate cyclase (sGC), the major receptor for nitric oxide (NO). Platelet analysis showed a complete loss of the soluble α1ß1 guanylate cyclase and showed an unexpected stimulatory role of sGC within platelets. The NO-sGC-cGMP pathway is a major pathway controlling vascular smooth-muscle relaxation, vascular tone, and vascular remodeling. Our data suggest that alterations of this pathway might lead to an abnormal vascular-remodeling process in sensitive vascular areas such as ICA bifurcations. These data provide treatment options for affected individuals and strongly suggest that investigation of GUCY1A3 and other members of the NO-sGC-cGMP pathway is warranted in both isolated early-onset achalasia and nonsyndromic moyamoya.

Robust physical methods that enrich genomic regions identical by descent for linkage studies: confirmation of a locus for osteogenesis imperfecta.

BACKGROUND: The monogenic disease osteogenesis imperfecta (OI) is due to single mutations in either of the collagen genes ColA1 or ColA2, but within the same family a given mutation is accompanied by a wide range of disease severity. Although this phenotypic variability implies the existence of modifier gene variants, genome wide scanning of DNA from OI patients has not been reported. Promising genome wide marker-independent physical methods for identifying disease-related loci have lacked robustness for widespread applicability. Therefore we sought to improve these methods and demonstrate their performance to identify known and novel loci relevant to OI. RESULTS: We have improved methods for enriching regions of identity-by-descent (IBD) shared between related, afflicted individuals. The extent of enrichment exceeds 10- to 50-fold for some loci. The efficiency of the new process is shown by confirmation of the identification of the Col1A2 locus in osteogenesis imperfecta patients from Amish families. Moreover the analysis revealed additional candidate linkage loci that may harbour modifier genes for OI; a locus on chromosome 1q includes COX-2, a gene implicated in osteogenesis. CONCLUSION: Technology for physical enrichment of IBD loci is now robust and applicable for finding genes for monogenic diseases and genes for complex diseases. The data support the further investigation of genetic loci other than collagen gene loci to identify genes affecting the clinical expression of osteogenesis imperfecta. The discrimination of IBD mapping will be enhanced when the IBD enrichment procedure is coupled with deep resequencing.

Convergent evidence identifying MAP/microtubule affinity-regulating kinase 1 (MARK1) as a susceptibility gene for autism.

Autism spectrum disorders (ASDs) are common, heritable, but genetically heterogeneous neurodevelopmental conditions. We recently defined a susceptibility locus for ASDs on chromosome 1q41-q42. High-resolution single-nucleotide polymorphisms (126 SNPs) genotyping across the chromosome 1q41-q42 region, followed by a MARK1 (microtubule affinity-regulating kinase 1)-tagged-SNP association study in 276 families with autism from the Autism Genetic Research Exchange, showed that several SNPs within the MARK1 gene were significantly associated with ASDs by transmission disequilibrium tests. Haplotype rs12740310*C-rs3737296*G-rs12410279*A was overtransmitted (P(corrected)= 0.0016), with a relative risk for autism of 1.8 in homozygous carriers. Furthermore, ASD-associated SNP rs12410279 modulates the level of transcription of MARK1. We found that MARK1 was overexpressed in the prefrontal cortex (BA46) but not in cerebellar granule cells, on postmortem brain tissues from patients. MARK1 displayed an accelerated evolution along the lineage leading to humans, suggesting possible involvement of this gene in cognition. MARK1 encodes a kinase-regulating microtubule-dependent transport in axons and dendrites. Both overexpression and silencing of MARK1 resulted in significantly shorter dendrite length in mouse neocortical neurons and modified dendritic transport speed. As expected for a gene encoding a key polarity determinant Par-1 protein kinase, MARK1 is involved in axon-dendrite specification. Thus, MARK1 overexpression in humans may be responsible for subtle changes in dendritic functioning.

Association of autism with polymorphisms in the paired-like homeodomain transcription factor 1 (PITX1) on chromosome 5q31: a candidate gene analysis.

BACKGROUND: Autism is a complex, heterogeneous, behaviorally-defined disorder characterized by disruptions of the nervous system and of other systems such as the pituitary-hypothalamic axis. In a previous genome wide screen, we reported linkage of autism with a 1.2 Megabase interval on chromosome 5q31. For the current study, we hypothesized that 3 of the genes in this region could be involved in the development of autism: 1) paired-like homeodomain transcription factor 1 (PITX1), which is a key regulator of hormones within the pituitary-hypothalamic axis, 2) neurogenin 1, a transcription factor involved in neurogenesis, and 3) histone family member Y (H2AFY), which is involved in X-chromosome inactivation in females and could explain the 4:1 male:female gender distortion present in autism. METHODS: A total of 276 families from the Autism Genetic Resource Exchange (AGRE) repository composed of 1086 individuals including 530 affected children were included in the study. Single nucleotide polymorphisms tagging the three candidate genes were genotyped on the initial linkage sample of 116 families. A second step of analysis was performed using tightly linked SNPs covering the PITX1 gene. Association was evaluated using the FBAT software version 1.7.3 for single SNP analysis and the HBAT command from the same package for haplotype analysis respectively. RESULTS: Association between SNPs and autism was only detected for PITX1. Haplotype analysis within PITX1 showed evidence for overtransmission of the A-C haplotype of markers rs11959298 - rs6596189 (p = 0.0004). Individuals homozygous or heterozygous for the A-C haplotype risk allele were 2.54 and 1.59 fold more likely to be autistic than individuals who were not carrying the allele, respectively. CONCLUSION: Strong and consistent association was observed between a 2 SNPs within PITX1 and autism. Our data suggest that PITX1, a key regulator of hormones within the pituitary-hypothalamic axis, may be implicated in the etiology of autism.

On confidence intervals for genotype relative risks and attributable risks from case parent trio designs for candidate-gene studies.

Scherag et al. [Hum Hered 2002;54:210-217] recently proposed point estimates and asymptotic as well as exact confidence intervals for genotype relative risks (GRRs) and the attributable risk (AR) in case parent trio designs using single nucleotide polymorphism (SNP) data. The aim of this study was the investigation of coverage probabilities and bias in estimates if the marker locus is not identical to the disease locus. Using a variety of parameter constellations, including marker allele frequencies identical to and different from the SNP at the disease locus, we performed an analytical study to quantify the bias and a Monte-Carlo simulation study for quantifying both bias and coverage probabilities. No bias was observed if marker and trait locus coincided. Two parameters had a strong impact on coverage probabilities of confidence intervals and bias in point estimates if they did not coincide: the linkage disequilibrium (LD) parameter delta and the allele frequency at the marker SNP. If marker allele frequencies were different from the allele frequencies at the functional SNP, substantial biases occurred. Further, if delta between the marker and the disease locus was lower than the maximum possible delta, estimates were also biased. In general, biases were towards the null hypothesis for both GRRs and AR. If one GRR was not increased, as e.g. in a recessive genetic model, biases away from the null could be observed. If both GRRs were in identical directions and if both were substantially larger than 1, the bias always was towards the null. When applying point estimates and confidence intervals for GRRs and AR in candidate gene studies, great care is needed. Effect estimates are substantially biased towards the null if either the allele frequencies at the marker SNP and the true disease locus are different or if the LD between the marker SNP and the disease locus is not at its maximum. A bias away from the null occurs only in uncommon study situations; it is small and can therefore be ignored for applications.

Evolution of neural precursor selection: functional divergence of proneural proteins.

How conserved pathways are differentially regulated to produce diverse outcomes is a fundamental question of developmental and evolutionary biology. The conserved process of neural precursor cell (NPC) selection by basic helix-loop-helix (bHLH) proneural transcription factors in the peripheral nervous system (PNS) by atonal related proteins (ARPs) presents an excellent model in which to address this issue. Proneural ARPs belong to two highly related groups: the ATONAL (ATO) group and the NEUROGENIN (NGN) group. We used a cross-species approach to demonstrate that the genetic and molecular mechanisms by which ATO proteins and NGN proteins select NPCs are different. Specifically, ATO group genes efficiently induce neurogenesis in Drosophila but very weakly in Xenopus, while the reverse is true for NGN group proteins. This divergence in proneural activity is encoded by three residues in the basic domain of ATO proteins. In NGN proteins, proneural capacity is encoded by the equivalent three residues in the basic domain and a novel motif in the second Helix (H2) domain. Differential interactions with different types of zinc (Zn)-finger proteins mediate the divergence of ATO and NGN activities: Senseless is required for ATO group activity, whereas MyT1 is required for NGN group function. These data suggest an evolutionary divergence in the mechanisms of NPC selection between protostomes and deuterostomes.

Evolutionary trace of G protein-coupled receptors reveals clusters of residues that determine global and class-specific functions.

G protein-coupled receptor (GPCR) activation mediated by ligand-induced structural reorganization of its helices is poorly understood. To determine the universal elements of this conformational switch, we used evolutionary tracing (ET) to identify residue positions commonly important in diverse GPCRs. When mapped onto the rhodopsin structure, these trace residues cluster into a network of contacts from the retinal binding site to the G protein-coupling loops. Their roles in a generic transduction mechanism were verified by 211 of 239 published mutations that caused functional defects. When grouped according to the nature of the defects, these residues sub-divided into three striking sub-clusters: a trigger region, where mutations mostly affect ligand binding, a coupling region near the cytoplasmic interface to the G protein, where mutations affect G protein activation, and a linking core in between where mutations cause constitutive activity and other defects. Differential ET analysis of the opsin family revealed an additional set of opsin-specific residues, several of which form part of the retinal binding pocket, and are known to cause functional defects upon mutation. To test the predictive power of ET, we introduced novel mutations in bovine rhodopsin at a globally important position, Leu-79, and at an opsin-specific position, Trp-175. Both were functionally critical, causing constitutive G protein activation of the mutants and rapid loss of regeneration after photobleaching. These results define in GPCRs a canonical signal transduction mechanism where ligand binding induces conformational changes propagated through adjacent trigger, linking core, and coupling regions.

Lysine 270 in the third intracellular domain of the oxytocin receptor is an important determinant for G alpha(q) coupling specificity.

To identify structural elements important to specific G alpha(q) coupling in the oxytocin receptor (OTR), intracellular domains were exchanged between OTR and G alpha(s)-coupled vasopressin V(2) receptors (V(2)Rs). Substitution of sequence from the second (2i) and third (3i) intracellular domains of V(2)R into comparable positions in OTR markedly reduced ligand affinity and resulted in a loss of G alpha(q) coupling. Substitution of the 2i domain of OTR into V(2)R decreased ligand affinity and vasopressin-stimulated adenylyl cyclase activity and only slightly increased phosphatidylinositide turnover. In contrast, substitution of the OTR3i domain into V(2)R produced a receptor chimera with high ligand affinity, decreased vasopressin-stimulated adenylyl cyclase activity, and markedly enhanced ligand-stimulated phosphatidylinositide turnover. The C-terminal 36 amino acids, but not the N-terminal 13 amino acids, of the OTR3i domain contained the determinants critical for enhanced activation of PLC. Mutation of a single lysine in the C-terminal OTR3i sequence to the corresponding V(2)R residue (valine) eliminated the enhanced ability of the V(2)R chimera to stimulate PLC but did not affect maximal adenylyl cyclase stimulation. Furthermore, mutation of this residue (K270) in wild-type OTR completely abolished the ability of the receptor to stimulate phosphatidylinositide turnover, with only a small reduction in ligand affinity. These data demonstrate that OTR K270 is critically important in the stimulation by OTR of phosphatidylinositide turnover and that this determinant can also increase this activity in the V(2)R chimera. Mutation of K270 also adversely affects the ability of OTR to stimulate ERK1/2 phosphorylation. Therefore, this residue plays an important role in the specificity of OTR/G alpha(q)/PLC coupling.

Structural clusters of evolutionary trace residues are statistically significant and common in proteins.

Given the massive increase in the number of new sequences and structures, a critical problem is how to integrate these raw data into meaningful biological information. One approach, the Evolutionary Trace, or ET, uses phylogenetic information to rank the residues in a protein sequence by evolutionary importance and then maps those ranked at the top onto a representative structure. If these residues form structural clusters, they can identify functional surfaces such as those involved in molecular recognition. Now that a number of examples have shown that ET can identify binding sites and focus mutational studies on their relevant functional determinants, we ask whether the method can be improved so as to be applicable on a large scale. To address this question, we introduce a new treatment of gaps resulting from insertions and deletions, which streamlines the selection of sequences used as input. We also introduce objective statistics to assess the significance of the total number of clusters and of the size of the largest one. As a result of the novel treatment of gaps, ET performance improves measurably. We find evolutionarily privileged clusters that are significant at the 5% level in 45 out of 46 (98%) proteins drawn from a variety of structural classes and biological functions. In 37 of the 38 proteins for which a protein-ligand complex is available, the dominant cluster contacts the ligand. We conclude that spatial clustering of evolutionarily important residues is a general phenomenon, consistent with the cooperative nature of residues that determine structure and function. In practice, these results suggest that ET can be applied on a large scale to identify functional sites in a significant fraction of the structures in the protein databank (PDB). This approach to combining raw sequences and structure to obtain detailed insights into the molecular basis of function should prove valuable in the context of the Structural Genomics Initiative.