beta2-adrenergic receptor signaling and desensitization elucidated by quantitative modeling of real time cAMP dynamics

G protein-coupled receptor signaling is dynamically regulated by multiple feedback mechanisms, which rapidly attenuate signals elicited by ligand stimulation, causing desensitization. The individual contributions of these mechanisms, however, are poorly understood. Here, we use an improved fluorescent biosensor for cAMP to measure second messenger dynamics stimulated by endogenous beta(2)-adrenergic receptor (beta(2)AR) in living cells. beta(2)AR stimulation with isoproterenol results in a transient pulse of cAMP, reaching a maximal concentration of approximately 10 microm and persisting for less than 5 min. We investigated the contributions of cAMP-dependent kinase, G protein-coupled receptor kinases, and beta-arrestin to the regulation of beta(2)AR signal kinetics by using small molecule inhibitors, small interfering RNAs, and mouse embryonic fibroblasts. We found that the cAMP response is restricted in duration by two distinct mechanisms in HEK-293 cells: G protein-coupled receptor kinase (GRK6)-mediated receptor phosphorylation leading to beta-arrestin mediated receptor inactivation and cAMP-dependent kinase-mediated induction of cAMP metabolism by phosphodiesterases. A mathematical model of beta(2)AR signal kinetics, fit to these data, revealed that direct receptor inactivation by cAMP-dependent kinase is insignificant but that GRK6/beta-arrestin-mediated inactivation is rapid and profound, occurring with a half-time of 70 s. This quantitative system analysis represents an important advance toward quantifying mechanisms contributing to the physiological regulation of receptor signaling.     

Genetic determinants of acute hypoxic ventilation: patterns of inheritance in mice.

Acutely lowering ambient O(2) tension increases ventilation in many mammalian species, including humans and mice. Inheritance patterns among kinships and between mouse strains suggest that a robust genetic influence determines individual hypoxic ventilatory responses (HVR). Here, we tested specific genetic hypotheses to describe the inheritance patterns of HVR phenotypes among two inbred mouse strains and their segregant and nonsegregant progeny. Using whole body plethysmography, we assessed the magnitude and pattern of ventilation in C3H/HeJ (C3) and C57BL/6J (B6) progenitor strains at baseline and during acute (3-5 min) hypoxic [mild hypercapnic hypoxia, inspired O(2) fraction (FI(O(2))) = 0.10] and normoxic (mild hypercapnic normoxia, FI(O(2)) = 0.21) inspirate challenges in mild hypercapnia (inspired CO(2) fraction = 0.03). First- and second-filial generations and two backcross progeny were also studied to assess response distributions of HVR phenotypes relative to the parental strains. Although the minute ventilation (VE) during hypoxia was comparable between the parental strains, breathing frequency (f) and tidal volume were significantly different; C3 mice demonstrated a slow, deep HVR relative to a rapid, shallow phenotype of B6 mice. The HVR profile in B6C3F(1)/J mice suggested that this offspring class represented a third phenotype, distinguishable from the parental strains. The distribution of HVR among backcross and intercross offspring suggested that the inheritance patterns for f and VE during mild hypercapnic hypoxia are consistent with models that incorporate two genetic determinants. These results further suggest that the quantitative genetic expression of alleles derived from C3 and B6 parental strains interact to significantly attenuate individual HVR in the first- and second-filial generations. In conclusion, the genetic control of HVR in this model was shown to exhibit a relatively simple genetic basis in terms of respiratory timing characteristics.     

Regulation of cell signaling dynamics by the protein kinase-scaffold Ste5

Cell differentiation requires the ability to detect and respond appropriately to a variety of extracellular signals. Here we investigate a differentiation switch induced by changes in the concentration of a single stimulus. Yeast cells exposed to high doses of mating pheromone undergo cell division arrest. Cells at intermediate doses become elongated and divide in the direction of a pheromone gradient (chemotropic growth). Either of the pheromone-responsive MAP kinases, Fus3 and Kss1, promotes cell elongation, but only Fus3 promotes chemotropic growth. Whereas Kss1 is activated rapidly and with a graded dose-response profile, Fus3 is activated slowly and exhibits a steeper dose-response relationship (ultrasensitivity). Fus3 activity requires the scaffold protein Ste5; when binding to Ste5 is abrogated, Fus3 behaves like Kss1, and the cells no longer respond to a gradient or mate efficiently with distant partners. We propose that scaffold proteins serve to modulate the temporal and dose-response behavior of the MAP kinase.     

Using the optimal receiver operating characteristic curve to design a predictive genetic test, exemplified with type 2 diabetes

Current extensive genetic research into common complex diseases, especially with the completion of genome-wide association studies, is bringing to light many novel genetic risk loci. These new discoveries, along with previously known genetic risk variants, offer an important opportunity for researchers to improve health care. We describe a method of quick evaluation of these new findings for potential clinical practice by designing a new predictive genetic test, estimating its classification accuracy, and determining the sample size required for the verification of this accuracy. The proposed predictive test is asymptotically more powerful than tests built on any other existing method and can be extended to scenarios where loci are linked or interact. We illustrate the approach for the case of type 2 diabetes. We incorporate recently discovered risk factors into the proposed test and find a potentially better predictive genetic test. The area under the receiver operating characteristic (ROC) curve (AUC) of the proposed test is estimated to be higher (AUC = 0.671) than for the existing test (AUC = 0.580).     

Optimizing the evidence for linkage by permuting marker order

We developed a new marker-reordering algorithm to find the best order of fine-mapping markers for multipoint linkage analysis. The algorithm searches for the best order of fine-mapping markers such that the sum of the squared differences in identity-by-descent distribution between neighboring markers is minimized. To test this algorithm, we examined its effect on the evidence for linkage in the simulated and the Collaborative Studies on Genetics of Alcoholism (COGA) data. We found enhanced evidence for linkage with the reordered map at the true location in the simulated data (p-value decreased from 1.16 x 10(-9) to 9.70 x 10(-10)). Analysis of the White population from the COGA data with the reordered map for alcohol dependence led to a significant change of the linkage signal (p = 0.0365 decreased to p = 0.0039) on chromosome 1 between marker D1S1592 and D1S1598. Our results suggest that reordering fine-mapping markers in candidate regions when the genetic map is uncertain can be a critical step when considering a dense map.     

Herpes virus in juvenile Pacific oysters Crassostrea gigas from Tomales Bay, California, coincides with summer mortality episodes

Pacific Crassostrea gigas and eastern C. virginica oysters were examined between June 2002 and April 2003 from 8 locations along the east, west and south USA coasts for oyster herpes virus (OsHV) infections using the A primer set in a previously developed PCR test. Only surviving Pacific oysters from a mortality event in Tomales Bay, California, USA, where annual losses of oysters have occurred each summer since 1993, were infected with a herpes-like virus in 2002. PCR examination using template amounts of both 50 and 500 ng were essential for OsHV detection. Sequence analysis indicated that the Tomales Bay OsHV was similar to that identified in France with the exception of a single base pair substitution in a 917 bp fragment of the viral genome. However, unlike the French OsHV-1, the Tomales Bay OsHV did not amplify with the primer pair of a second OsHV-1 PCR assay, suggesting that further characterization of these viruses is warranted. No evidence of Cowdry type A viral infections characteristic of herpes virus infections or other pathogens were observed in OsHV-infected oysters. Hemocytosis, diapedesis and hemocyte degeneration characterized by nuclear pycnosis and fragmentation were observed in infected oysters, which is consistent with previous observations of OsHV infections in France. Together these data suggest that OsHV may be associated with the annual summer Pacific oyster seed mortality observed in Tomales Bay, but establishment of a causal relationship warrants further investigation.     

Fine mapping of autistic disorder to chromosome 15q11-q13 by use of phenotypic subtypes

Autistic disorder (AutD) is a complex genetic disease. Available evidence suggests that several genes contribute to the underlying genetic risk for the development of AutD. However, both etiologic heterogeneity and genetic heterogeneity confound the discovery of AutD-susceptibility genes. Chromosome 15q11-q13 has been identified as a strong candidate region on the basis of both the frequent occurrence of chromosomal abnormalities in that region and numerous suggestive linkage and association findings. Ordered-subset analysis (OSA) is a novel statistical method to identify a homogeneous subset of families that contribute to overall linkage at a given chromosomal location and thus to potentially help in the fine mapping and localization of the susceptibility gene within a chromosomal area. For the present analysis, a factor that represents insistence on sameness (IS)--derived from a principal-component factor analysis using data on 221 patients with AutD from the repetitive behaviors/stereotyped patterns domain in the Autism Diagnostic Interview-Revised--was used as a covariate in OSA. Analysis of families sharing high scores on the IS factor increased linkage evidence for the 15q11-q13 region, at the GABRB3 locus, from a LOD score of 1.45 to a LOD score of 4.71. These results narrow our region of interest on chromosome 15 to an area surrounding the gamma-aminobutyric acid-receptor subunit genes, in AutD, and support the hypothesis that the analysis of phenotypic homogeneous subtypes may be a powerful tool for the mapping of disease-susceptibility genes in complex traits.     

The structure of linkage disequilibrium at the DBH locus strongly influences the magnitude of association between diallelic markers and plasma dopamine beta-hydroxylase activity

There is currently a great deal of interest in using linkage disequilibrium (LD) mapping to locate both disease and quantitative-trait loci on a genomewide scale. Recent findings suggest that much of the human genome is organized in discrete "blocks" of low haplotype diversity, but the utility of such blocks in identifying genes influencing complex traits is not yet known and must ultimately be tested empirically through use of real data. We recently identified a putative functional polymorphism (-1021C-->T) in the 5' upstream region of the DBH gene that accounted for 35%-52% of the total phenotypic variance in plasma dopamine beta-hydroxylase (DBH) activity in samples from three distinct populations. In the present study, we genotyped 11 diallelic markers at the DBH locus surrounding -1021C-->T in 386 unrelated individuals of European origin. We identified a single 10-kb block containing -1021C-->T, in which four haplotypes comprised 93% of the observed chromosomes. Only markers within the block were highly associated with phenotype (P< or =2.2 x 10(-10)), with one exception. In general, association with phenotype was strongly correlated with the degree of LD between each marker and -1021C-->T. Of four LD measures assessed, d(2) was the best predictor of this relationship. Had one attempted to map quantitative-trait loci for plasma DBH activity on a genomewide basis without prior knowledge of candidate regions and not included (by chance) markers within this haplotype block, the DBH locus might have been missed entirely. These results provide a direct example of the potential value of constructing a haplotype map of the human genome prior to embarking on large-scale association studies.     

X-linked extension of the revised Haseman-Elston algorithm for linkage analysis in sib pairs

Haseman and Elston (H-E) proposed a regression-based robust test of linkage between a marker and an autosomal quantitative trait locus, using the squared sib pair trait difference as a dependent variable and the proportion of alleles shared identical by descent by the sib pair as an independent variable. Several authors have proposed improvement of the original H-E's seminal work by using an optimal linear combination of squared sum and squared difference as the dependent variable. In this paper, we extend Haseman and Elston's sib pair method to an X-linked locus. We give a general formulation of the complete regression model and details of the regression coefficients in terms of variance components. Simulation results are presented to describe the power of this technique for a theoretical best case scenario.