Imputation strategies for genomic prediction using nanopore sequencing.

BACKGROUND: Genomic prediction describes the use of SNP genotypes to predict complex traits and has been widely applied in humans and agricultural species. Genotyping-by-sequencing, a method which uses low-coverage sequence data paired with genotype imputation, is becoming an increasingly popular SNP genotyping method for genomic prediction. The development of Oxford Nanopore Technologies' (ONT) MinION sequencer has now made genotyping-by-sequencing portable and rapid. Here we evaluate the speed and accuracy of genomic predictions using low-coverage ONT sequence data in a population of cattle using four imputation approaches. We also investigate the effect of SNP reference panel size on imputation performance. RESULTS: SNP array genotypes and ONT sequence data for 62 beef heifers were used to calculate genomic estimated breeding values (GEBVs) from 641 k SNP for four traits. GEBV accuracy was much higher when genome-wide flanking SNP from sequence data were used to help impute the 641 k panel used for genomic predictions. Using the imputation package QUILT, correlations between ONT and low-density SNP array genomic breeding values were greater than 0.91 and up to 0.97 for sequencing coverages as low as 0.1 × using a reference panel of 48 million SNP. Imputation time was significantly reduced by decreasing the number of flanking sequence SNP used in imputation for all methods. When compared to high-density SNP arrays, genotyping accuracy and genomic breeding value correlations at 0.5 × coverage were also found to be higher than those imputed from low-density arrays. CONCLUSIONS: Here we demonstrated accurate genomic prediction is possible with ONT sequence data from sequencing coverages as low as 0.1 × , and imputation time can be as short as 10 min per sample. We also demonstrate that in this population, genotyping-by-sequencing at 0.1 × coverage can be more accurate than imputation from low-density SNP arrays.

Molecular hypotheses to explain the shared pathways and underlying pathobiological causes in catatonia and in catatonic presentations in neuropsychiatric disorders.

The pathobiological causes, the shared cellular and molecular pathways in catatonia and in catatonic presentation in neuropsychiatric disorders are yet to be determined. The hypotheses in this paper have been deduced from the latest scientific research findings and clinical observations of patients with genetic disorders, behavioral phenotypes and other family members suffering mental disorders. The first hypothesis postulates that catatonia and the heterogeneity of catatonic signs and symptoms involve nucleolar dysfunction arising from abnormalities of the brain-specific, non-coding micro-RNA, SNORD115 genes (either duplications or deletions) which result in pathobiological dysfunction of various combinations in the downstream pathways (possibly along with other genes in these shared pathways). SNORD115 controls five genes CRHR1, PBRM1, TAF1, DPM2, and RALGPS1 as well as the alternative splicing of serotonin 2C receptor. SNORD115 abnormalities with varying downstream multigene involvement would account for catatonia across the life span within some subtypes of autism spectrum disorders, schizophrenia, bipolar and major depressive disorder, psychosis, genetic disorders, and in immune disorders such as anti-N-methyl-d-aspartate receptor (NMDAR) antibody encephalitis as well as the susceptibility to the neuroleptic malignant syndrome (NMS) if environmentally triggered. Furthermore, SNORD115 genes may underlie a genetic vulnerability when environmental triggers result in excess serotonin producing the serotonin syndrome, a condition similar to NMS in which catatonia may occur. Dysfunction of SNORD115-PBRM1 connecting with SMARCA2 as well as other proven schizophrenia-associated genes might explain why traditionally catatonia has been classified with schizophrenia. SNORD115-TAF1 and SNORD-DPM2 dysfunction introduce possible clues to the parkinsonism and increased creatinine phosphokinase in NMS, while abnormalities of SNORD115-RALGPS1 suggest links to both anti-NMDAR encephalitis and the proven predisposing catatonic SHANK3 gene. The second hypothesis postulates that periodic catatonia (PC) on 15q15 involves abnormalities of vacuolar protein sorting 39 (VPS39), a proven de novo schizophrenic gene in this chromosomal locus and part of the HOPS complex. These will impact the autophagic and endocytic pathways, thereby lowering lysosomal degradation. VPS39 mutations may be considered also to disrupt lysosome-mitochondria tethering and transport of lipids and calcium through membrane contact sites (MCSs). To account for the periodicity in PC it is speculated that the mammalian equivalent of the vacuole and mitochondria patch (vCLAMP) would be altered by VPS39 mutations and subsequently followed by the mammalian equivalent of endoplasmic reticulum mitochondria encounter structure (ERMES) restoring mitochondrial homeostasis. Future precision psychiatry will require accurate pathophysiologically-defined psychiatric diagnoses to accelerate the discovery of specific molecular-targeted medications to improve therapeutic outcomes.

Investigating the effect of two methane-mitigating diets on the rumen microbiome using massively parallel sequencing.

Variation in the composition of microorganisms in the rumen (the rumen microbiome) of dairy cattle (Bos taurus) is of great interest because of possible links to methane emission levels. Feed additives are one method being investigated to reduce enteric methane production by dairy cattle. Here we report the effect of 2 methane-mitigating feed additives (grapemarc and a combination of lipids and tannin) on rumen microbiome profiles of Holstein dairy cattle. We used untargeted (shotgun) massively parallel sequencing of microbes present in rumen fluid to generate quantitative rumen microbiome profiles. We observed large effects of the feed additives on the rumen microbiome profiles using multiple approaches, including linear mixed modeling, hierarchical clustering, and metagenomic predictions. The effect on the fecal microbiome profiles was not detectable using hierarchical clustering, but was significant in the linear mixed model and when metagenomic predictions were used, suggesting a more subtle effect of the diets on the lower gastrointestinal microbiome. A differential representation analysis (analogous to differential expression in RNA sequencing) showed significant overlap in the contigs (which are genome fragments representing different microorganism species) that were differentially represented between experiments. These similarities suggest that, despite the different additives used, the 2 diets assessed in this investigation altered the microbiomes of the samples in similar ways. Contigs that were differentially represented in both experiments were tested for associations with methane production in an independent set of animals. These animals were not treated with a methane-mitigating diet, but did show substantial natural variation in methane emission levels. The contigs that were significantly differentially represented in response to both dietary additives showed a significant enrichment for associations with methane production. This suggests that these methane-mitigating diets have altered the rumen microbiome toward naturally low methane-emitting microbial profiles. The contig sequences are predominantly new and include Faecalibacterium spp. The contigs we have identified here are potential biomarkers for low-methane-emitting cattle.

Human herpesviruses-6 and -7 each cause significant neurological morbidity in Britain and Ireland.

BACKGROUND: Primary human herpesvirus-6 and -7 (HHV-6/-7) infections cause febrile illness sometimes complicated by convulsions and rarely encephalopathy. AIMS: To explore the extent of such HHV-6 and -7 induced disease in young children. METHODS: In a three year prospective study in Britain and Ireland, 205 children (2-35 months old) hospitalised with suspected encephalitis and/or severe illness with fever and convulsions were reported via the British Paediatric Surveillance Unit network. Blood samples were tested for primary HHV-6 and -7 infections. RESULTS: 26/156 (17%) of children aged 2-23 months had primary infection (11 HHV-6; 13 HHV-7; two with both viruses) coinciding with the acute illness; this was much higher than the about three cases expected by chance. All 26 were pyrexial; 25 had convulsions (18 status epilepticus), 11 requiring ventilation. Median hospital stay was 7.5 days. For HHV-6 primary infection the median age was 53 weeks (range 42-94) and the distribution differed from that of uninfected children; for HHV-7, the median was 60 weeks (range 17-102) and the distribution did not differ for the uninfected. Fewer (5/15) children with primary HHV-7 infection had previously been infected with HHV-6 than expected. CONCLUSIONS: Primary HHV-6 and HHV-7 infections accounted for a significant proportion of cases in those <2 years old of severe illness with fever and convulsions requiring hospital admission; each virus contributed equally. Predisposing factors are age for HHV-6 and no previous infection with HHV-6 for HHV-7. Children with such neurological disease should be investigated for primary HHV-6/-7 infections, especially in rare cases coinciding by chance with immunisation to exclude misdiagnosis as vaccine reactions.

RGSZ1 and Ret RGS: two of several splice variants from the gene RGS20.

RGSZ1 and Ret RGS, members of the regulator of G-protein signaling (RGS) family, are GTPase-activating proteins (GAPs) with high selectivity for G alpha(z). We show here that RGSZ1 and Ret RGSZ1 are products of two of several splice variants of one gene, RGS20. RGS20 spans approximately 107 kb and contains at least seven exons. Five exons account for RGSZ1, including a single exon distinct to RGSZ1 that encodes a newly identified amino-terminal region. The previously described open reading frame (ORF) and 3' untranslated region are encoded by four downstream exons that also encode about half of Ret RGS. The 5' end of the RGSZ1 ORF contains several in-frame ATG codons (3-5 depending on the species), and multiple translational start sites may help explain the molecular weight heterogeneity of purified bovine brain RGSZ. Ret RGS replaces the 24 N-terminal amino acid residues of RGSZ1 with a large, N-terminal region that initially distinguished the bovine Ret RGS from human and mouse RGSZ1. This N-terminal domain is encoded by two distinct 5' exons that are variably combined with the four downstream exons shared with RGSZ1 to produce at least six mRNAs. They encode proteins with N termini that vary in size, hydrophobicity, and the presence of a cysteine string. At least two mRNAs that include the exon that encodes the N-terminal region unique to RGSZ1 were found in brain and a few other tissues, but not retina. RGS20 thus can account for multiple G(z)-selective GAPs in different tissues.

Binding of regulator of G protein signaling (RGS) proteins to phospholipid bilayers. Contribution of location and/or orientation to Gtpase-activating protein activity.

Regulator of G protein signaling (RGS) proteins must bind membranes in an orientation that permits the protein-protein interactions necessary for regulatory activity. RGS4 binds to phospholipid surfaces in a slow, multistep process that leads to maximal GTPase-activating protein (GAP) activity. When RGS4 is added to phospholipid vesicles that contain m2 or m1 muscarinic receptor and G(i), G(z), or G(q), GAP activity increases approximately 3-fold over 4 h at 30 degrees C and more slowly at 20 degrees C. This increase in GAP activity is preceded by several other events that suggest that, after binding, optimal interaction with G protein and receptor requires reorientation of RGS4 on the membrane surface, a conformational change, or both. Binding of RGS4 is initially reversible but becomes irreversible within 5 min. Onset of irreversibility parallels initial quenching of tryptophan fluorescence (t(12) approximately 30 s). Further quenching occurs after binding has become irreversible (t(12) approximately 6 min) but is complete well before maximal GAP activity is attained. These processes all appear to be energetically driven by the amphipathic N-terminal domain of RGS4 and are accelerated by palmitoylation of cysteine residues in this region. The RGS4 N-terminal domain confers similar membrane binding behavior on the RGS domains of either RGS10 or RGSZ1.

GTPase-activating proteins for heterotrimeric G proteins: regulators of G protein signaling (RGS) and RGS-like proteins.

GTPase-activating proteins (GAPs) regulate heterotrimeric G proteins by increasing the rates at which their subunits hydrolyze bound GTP and thus return to the inactive state. G protein GAPs act allosterically on G subunits, in contrast to GAPs for the Ras-like monomeric GTP-binding proteins. Although they do not contribute directly to the chemistry of GTP hydrolysis, G protein GAPs can accelerate hydrolysis >2000-fold. G protein GAPs include both effector proteins (phospholipase C-¿, p115RhoGEF) and a growing family of regulators of G protein signaling (RGS proteins) that are found throughout the animal and fungal kingdoms. GAP activity can sharpen the termination of a signal upon removal of stimulus, attenuate a signal either as a feedback inhibitor or in response to a second input, promote regulatory association of other proteins, or redirect signaling within a G protein signaling network. GAPs are regulated by various controls of their cellular concentrations, by complex interactions with G¿ or with G¿5 through an endogenous G-like domain, and by interaction with multiple other proteins.

Palmitoylation of a conserved cysteine in the regulator of G protein signaling (RGS) domain modulates the GTPase-activating activity of RGS4 and RGS10.

RGS4 and RGS10 expressed in Sf9 cells are palmitoylated at a conserved Cys residue (Cys(95) in RGS4, Cys(66) in RGS10) in the regulator of G protein signaling (RGS) domain that is also autopalmitoylated when the purified proteins are incubated with palmitoyl-CoA. RGS4 also autopalmitoylates at a previously identified cellular palmitoylation site, either Cys(2) or Cys(12). The C2A/C12A mutation essentially eliminates both autopalmitoylation and cellular [(3)H]palmitate labeling of Cys(95). Membrane-bound RGS4 is palmitoylated both at Cys(95) and Cys(2/12), but cytosolic RGS4 is not palmitoylated. RGS4 and RGS10 are GTPase-activating proteins (GAPs) for the G(i) and G(q) families of G proteins. Palmitoylation of Cys(95) on RGS4 or Cys(66) on RGS10 inhibits GAP activity 80-100% toward either Galpha(i) or Galpha(z) in a single-turnover, solution-based assay. In contrast, when GAP activity was assayed as acceleration of steady-state GTPase in receptor-G protein proteoliposomes, palmitoylation of RGS10 potentiated GAP activity >/=20-fold. Palmitoylation near the N terminus of C95V RGS4 did not alter GAP activity toward soluble Galpha(z) and increased G(z) GAP activity about 2-fold in the vesicle-based assay. Dual palmitoylation of wild-type RGS4 remained inhibitory. RGS protein palmitoylation is thus multi-site, complex in its control, and either inhibitory or stimulatory depending on the RGS protein and its sites of palmitoylation.

Reciprocal signaling between heterotrimeric G proteins and the p21-stimulated protein kinase.

p21-activated protein kinase (PAK)-1 phosphorylated Galpha(z), a member of the Galpha(i) family that is found in the brain, platelets, and adrenal medulla. Phosphorylation approached 1 mol of phosphate/mol of Galpha(z) in vitro. In transfected cells, Galpha(z) was phosphorylated both by wild-type PAK1 when stimulated by the GTP-binding protein Rac1 and by constitutively active PAK1 mutants. In vitro, phosphorylation occurred only at Ser(16), one of two Ser residues that are the major substrate sites for protein kinase C (PKC). PAK1 did not phosphorylate other Galpha subunits (i1, i2, i3, o, s, or q). PAK1-phosphorylated Galpha(z) was resistant both to RGSZ1, a G(z)-selective GTPase-activating protein (GAP), and to RGS4, a relatively nonselective GAP for the G(i) and G(q) families of G proteins. Phosphorylation of Ser(27) by PKC did not alter sensitivity to either GAP. The previously described inhibition of G(z) GAPs by PKC is therefore mediated by phosphorylation of Ser(16). Phosphorylation of either Ser(16) by PAK1 or Ser(27) by PKC decreased the affinity of Galpha(z) for Gbetagamma; phosphorylation of both residues by PKC caused no further effect. PAK1 thus regulates Galpha(z) function by attenuating the inhibitory effects of both GAPs and Gbetagamma. In this context, the kinase activity of PAK1 toward several protein substrates was directly inhibited by Gbetagamma, suggesting that PAK1 acts as a Gbetagamma-regulated effector protein. This inhibition of mammalian PAK1 by Gbetagamma contrasts with the stimulation of the PAK homolog Ste20p in Saccharomyces cerevisiae by the Gbetagamma homolog Ste4p/Ste18p.

Rapid GTP binding and hydrolysis by G(q) promoted by receptor and GTPase-activating proteins.

Receptor-promoted GTP binding and GTPase-activating protein (GAP)-promoted GTP hydrolysis determine the onset and termination of G protein signaling; they coordinately control signal amplitude. The mechanisms whereby cells independently regulate signal kinetics and signal amplitude are therefore central to understanding G protein function. We have used quench-flow kinetic methods to measure the rates of the individual reactions of the agonist-stimulated GTPase cycle for G(q) during steady-state signaling. G(q) and m1 muscarinic cholinergic receptor were co-reconstituted into proteoliposomes with one of two GAPs: phospholipase C (PLC)-beta1, the major G(q)-regulated effector protein, and RGS4, a GAP commonly thought to be an inhibitor of G(q) signaling. In this system, the rate constant for GAP-stimulated hydrolysis of Galpha(q)-bound GTP at 30 degrees C was 9-12 s(-1) for PLC-beta1 and 22-27 s(-1) for RGS4. These rates are 1,000- to 2,000-fold faster than in the absence of a GAP and far faster than measured previously. G(q) can thus hydrolyze bound GTP with deactivation half-times of 25-75 ms at 30 degrees C, commensurate with physiological rates of signal termination. GDP/GTP exchange, which reactivates G(q), was the principal rate-limiting step for the GTPase cycle and was also faster than previously thought. At physiological concentrations of GTP, exchange was limited by the rate of dissociation of GDP from the receptor-G(q) complex, with a maximal rate of 1.8 s(-1) at 30 degrees C. Comparison of activation and deactivation rates help explain how GDP/GTP exchange balance rapid GTP hydrolysis to maintain steady-state signal amplitude.

Modulation of the affinity and selectivity of RGS protein interaction with G alpha subunits by a conserved asparagine/serine residue.

The crystal structure of the complex between a G protein alpha subunit (Gi alpha 1) and its GTPase-activating protein (RGS4) demonstrated that RGS4 acts predominantly by stabilization of the transition state for GTP hydrolysis [Tesmer, J. J., et al. (1997) Cell 89, 251]. However, attention was called to a conserved Asn residue (Asn128) that could play a catalytic role by interacting, directly or indirectly, with the hydrolytic water molecule. We have analyzed the effects of several disparate substitutions for Asn128 on the GAP activity of RGS4 toward four G alpha substrates (Go, Gi, Gq, and Gz) using two assay formats. The results substantiate the importance of this residue but indicate that it is largely involved in substrate binding and that its function may vary with different G alpha targets. Various mutations decreased the apparent affinity of RGS4 for substrate G alpha proteins by several orders of magnitude, but had variable and modest effects on maximal rates of GTP hydrolysis when tested with different G alpha subunits. One mutation, N128F, that differentially decreased the GAP activity toward G alpha i compared with that toward G alpha q could be partially suppressed by mutation of the nearby residue in G alpha i to that found in G alpha q (K180P). Detection of GAP activities of the mutants was enhanced in sensitivity up to 100-fold by assay at steady state in proteoliposomes that contain heterotrimeric G protein and receptor.

Individualizing therapy to prevent long-term consequences of estrogen deficiency in postmenopausal women.

BACKGROUND: Alendronate sodium and raloxifene hydrochloride were recently approved for the prevention of postmenopausal osteoporosis, but data on their clinical efficacy are limited. We compared these drugs with hormone replacement therapy (HRT) to help women and physicians guide postmenopausal treatment decisions. OBJECTIVE: To help physicians understand how they can best help women choose the most beneficial therapy after menopause based on their individual risk profile. METHODS: We developed a decision analytic Markov model to compare the effects of alendronate therapy, raloxifene therapy, and HRT on risks of hip fracture, coronary heart disease (CHD), breast cancer, and life expectancy. Regression models linked individual risk factors to future disease risks and were modified by drug effects on bone density, lipid levels, and associated breast cancer effects. RESULTS: Hormone replacement therapy, alendronate therapy, and raloxifene therapy have similar predicted efficacies in preventing hip fractures (estimated relative risk, 0.57, 0.54, and 0.58, respectively). Hormone replacement therapy should be more than 10 times more effective than raloxifene therapy in preventing CHD, but raloxifene therapy may not induce breast cancer. Women at low risk for hip fracture, CHD, and breast cancer do not benefit significantly from any treatment. Among women at average risk, HRT was preferred unless raloxifene therapy could reduce the risk of breast cancer by at least 66%, compared with a 47% increase for HRT. Women at high risk for CHD benefit most from HRT; women at high risk for breast cancer but low risk for CHD benefit most from raloxifene therapy, but only if it lowers the risk of breast cancer. CONCLUSION: Because of significant differences in the impact of these drugs, treatment choice depends on an individual woman's risk for hip fracture, CHD, and breast cancer.

Fitting nutrition into the medical model: the role of decision analytic cost-effectiveness techniques.

Physicians are accustomed to making decisions based on information regarding the prevalence of disease, symptoms, physical signs, laboratory test results, and the risks and benefits of alternative treatments. If nutritional assessment and therapeutics are to become more common components of medical practice, significant barriers in each of these areas must be overcome. Even rudimentary dietary assessment is often missing from physician education. Dietary assessment tools that are readily available and that have demonstrated usefulness are largely unknown. In addition, many nutritional interventions have not been formally investigated in randomized, controlled trials, and thus their cost-effectiveness remains unknown. We present one approach to these issues by discussing the construction of a decision model examining strategies for vitamin D and calcium screening. The application of medical decision making techniques to problems in clinical nutrition illustrates how findings from research studies may be used to determine the risks, benefits and costs of alternative population based health related nutrition policies which can then be applied by physicians in their daily interactions with patients.

Kinetic control of guanine nucleotide binding to soluble Galpha(q).

Binding of guanine nucleotides to heterotrimeric G proteins is controlled primarily by kinetic factors, such as the release of bound GDP, rather than by affinity alone. Detergent-solubilized Galpha(q) displays unusual guanine nucleotide binding properties in comparison with other G protein alpha subunits. Under conditions where most G proteins bind nearly stoichiometric GTPgammaS in 5-30 min at micromolar nucleotide concentrations, GTPgammaS binding to Galpha(q) is slow (>1 hr to completion), markedly substoichiometric, and dependent upon high concentrations of nucleotide (0.1 to 0.2 mM). Although the latter two properties suggest low affinity, GTPgammaS dissociation is immeasurably slow under commonly used conditions. We found that purified Galpha(q) can bind stoichiometric GTPgammaS, but that binding is controlled kinetically by a combination of factors. GDP (or IDP) dissociated slowly from Galpha(q), but the dissociation rate increased linearly with the concentration of (NH4)2SO4 up to 0.75 M (approximately 20-fold acceleration). The resulting GDP-free Galpha(q) was labile to rapid and irreversible denaturation, however (rate constant > or = 1 min(-1) at 20 degrees). Denaturation competed kinetically with relatively slow GTPgammaS association, such that stoichiometric binding was only attained at 100 microM GTPgammaS. These findings reconcile the slowly reversible binding of GTPgammaS to Galpha(q) with the other behaviors that suggested lower affinity, and point out that events subsequent to GDP dissociation can markedly influence the rates and extents of guanine nucleotide binding to G protein alpha subunits. Understanding these interactions allowed the direct, accurate quantitation of active Galpha(q) by a simple GTPgammaS binding assay in the presence of (NH4)2SO4, and similarly can prevent underestimation of the concentrations of other G proteins.

Phospholipase C-beta1 directly accelerates GTP hydrolysis by Galphaq and acceleration is inhibited by Gbeta gamma subunits.

Phospholipase C-beta, the principal effector protein regulated by Galphaq, has been shown to increase the agonist-stimulated, steady-state GTPase activity of Gq in proteoliposomes that contain both heterotrimeric Gq and m1 muscarinic receptor. We now use a moderately stable complex of R183C Galphaq bound to GTP to show that PLC-beta1 acts directly as a GTPase-activating protein (GAP) for isolated Galphaq in a membrane-free system. PLC-beta1 accelerated the hydrolysis of GalphaqR183C.GTP up to 20-fold. The Km was 1.5 nM, which is similar both to the EC50 with which R183C and wild type Galphaq activate PLC-beta1 and to the EC50 with which PLC-beta1 acts as a Gq GAP in the vesicle-based assay. The Galphaq GAP activity of RGS4 can also be quantitated by this assay; it accelerated hydrolysis of bound GTP about 100-fold. The Gq GAP activities of both PLC-beta1 and RGS4 are blocked by Gbeta gamma subunits, probably by a competitive mechanism. These data suggest either that the Gbeta gamma subunits are not continuously required for receptor-catalyzed GDP/GTP exchange during steady-state GTP hydrolysis or that GAPs, either PLC-beta or RGS proteins, can substitute for Gbeta gamma in this set of reactions.

The N-terminal domain of RGS4 confers receptor-selective inhibition of G protein signaling.

Regulators of heterotrimeric G protein signaling (RGS) proteins are GTPase-activating proteins (GAPs) that accelerate GTP hydrolysis by Gq and Gi alpha subunits, thus attenuating signaling. Mechanisms that provide more precise regulatory specificity have been elusive. We report here that an N-terminal domain of RGS4 discriminated among receptor signaling complexes coupled via Gq. Accordingly, deletion of the N-terminal domain of RGS4 eliminated receptor selectivity and reduced potency by 10(4)-fold. Receptor selectivity and potency of inhibition were partially restored when the RGS4 box was added together with an N-terminal peptide. In vitro reconstitution experiments also indicated that sequences flanking the RGS4 box were essential for high potency GAP activity. Thus, RGS4 regulates Gq class signaling by the combined action of two domains: 1) the RGS box accelerates GTP hydrolysis by Galphaq and 2) the N terminus conveys high affinity and receptor-selective inhibition. These activities are each required for receptor selectivity and high potency inhibition of receptor-coupled Gq signaling.

RGSZ1, a Gz-selective RGS protein in brain. Structure, membrane association, regulation by Galphaz phosphorylation, and relationship to a Gz gtpase-activating protein subfamily.

We cloned the cDNA for human RGSZ1, the major Gz-selective GTPase-activating protein (GAP) in brain (Wang, J., Tu, Y., Woodson, J., Song, X., and Ross, E. M. (1997) J. Biol. Chem. 272, 5732-5740) and a member of the RGS family of G protein GAPs. Its sequence is 83% identical to RET-RGS1 (except its N-terminal extension) and 56% identical to GAIP. Purified, recombinant RGSZ1, RET-RGS1, and GAIP each accelerated the hydrolysis of Galphaz-GTP over 400-fold with Km values of approximately 2 nM. RGSZ1 was 100-fold selective for Galphaz over Galphai, unusually specific among RGS proteins. Other enzymological properties of RGSZ1, brain Gz GAP, and RET-RGS1 were identical; GAIP differed only in Mg2+ dependence and in its slightly lower selectivity for Galphaz. RGSZ1, RET-RGS1, and GAIP thus define a subfamily of Gz GAPs within the RGS proteins. RGSZ1 has no obvious membrane-spanning region but is tightly membrane-bound in brain. Its regulatory activity in membranes depends on stable bilayer association. When co-reconstituted into phospholipid vesicles with Gz and m2 muscarinic receptors, RGSZ1 increased agonist-stimulated GTPase >15-fold with EC50 <12 nM, but RGSZ1 added to the vesicle suspension was <0.1% as active. RGSZ1, RET-RGS1, and GAIP share a cysteine string sequence, perhaps targeting them to secretory vesicles and allowing them to participate in the proposed control of secretion by Gz. Phosphorylation of Galphaz by protein kinase C inhibited the GAP activity of RGSZ1 and other RGS proteins, providing a mechanism for potentiation of Gz signaling by protein kinase C.

Dynamic regulation of RGS2 suggests a novel mechanism in G-protein signaling and neuronal plasticity.

Long-term neuronal plasticity is known to be dependent on rapid de novo synthesis of mRNA and protein, and recent studies provide insight into the molecules involved in this response. Here, we demonstrate that mRNA encoding a member of the regulator of G-protein signaling (RGS) family, RGS2, is rapidly induced in neurons of the hippocampus, cortex, and striatum in response to stimuli that evoke plasticity. Although several members of the RGS family are expressed in brain with discrete neuronal localizations, RGS2 appears unique in that its expression is dynamically responsive to neuronal activity. In biochemical assays, RGS2 stimulates the GTPase activity of the alpha subunit of Gq and Gi1. The effect on Gi1 was observed only after reconstitution of the protein in phospholipid vesicles containing M2 muscarinic acetylcholine receptors. RGS2 also inhibits both Gq- and Gi-dependent responses in transfected cells. These studies suggest a novel mechanism linking neuronal activity and signal transduction.