A widespread chromosomal inversion polymorphism contributes to a major life-history transition, local adaptation, and reproductive isolation.

The role of chromosomal inversions in adaptation and speciation is controversial. Historically, inversions were thought to contribute to these processes either by directly causing hybrid sterility or by facilitating the maintenance of co-adapted gene complexes. Because inversions suppress recombination when heterozygous, a recently proposed local adaptation mechanism predicts that they will spread if they capture alleles at multiple loci involved in divergent adaptation to contrasting environments. Many empirical studies have found inversion polymorphisms linked to putatively adaptive phenotypes or distributed along environmental clines. However, direct involvement of an inversion in local adaptation and consequent ecological reproductive isolation has not to our knowledge been demonstrated in nature. In this study, we discovered that a chromosomal inversion polymorphism is geographically widespread, and we test the extent to which it contributes to adaptation and reproductive isolation under natural field conditions. Replicated crosses between the prezygotically reproductively isolated annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus, revealed that alternative chromosomal inversion arrangements are associated with life-history divergence over thousands of kilometers across North America. The inversion polymorphism affected adaptive flowering time divergence and other morphological traits in all replicated crosses between four pairs of annual and perennial populations. To determine if the inversion contributes to adaptation and reproductive isolation in natural populations, we conducted a novel reciprocal transplant experiment involving outbred lines, where alternative arrangements of the inversion were reciprocally introgressed into the genetic backgrounds of each ecotype. Our results demonstrate for the first time in nature the contribution of an inversion to adaptation, an annual/perennial life-history shift, and multiple reproductive isolating barriers. These results are consistent with the local adaptation mechanism being responsible for the distribution of the two inversion arrangements across the geographic range of M. guttatus and that locally adaptive inversion effects contribute directly to reproductive isolation. Such a mechanism may be partially responsible for the observation that closely related species often differ by multiple chromosomal rearrangements.

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

Mitochondrial DNA polymorphism A4917G is independently associated with age-related macular degeneration.

The objective of this study was to determine if MTND2*LHON4917G (4917G), a specific non-synonymous polymorphism in the mitochondrial genome previously associated with neurodegenerative phenotypes, is associated with increased risk for age-related macular degeneration (AMD). A preliminary study of 393 individuals (293 cases and 100 controls) ascertained at Vanderbilt revealed an increased occurrence of 4917G in cases compared to controls (15.4% vs.9.0%, p = 0.11). Since there was a significant age difference between cases and controls in this initial analysis, we extended the study by selecting Caucasian pairs matched at the exact age at examination. From the 1547 individuals in the Vanderbilt/Duke AMD population association study (including 157 in the preliminary study), we were able to match 560 (280 cases and 280 unaffected) on exact age at examination. This study population was genotyped for 4917G plus specific AMD-associated nuclear genome polymorphisms in CFH, LOC387715 and ApoE. Following adjustment for the listed nuclear genome polymorphisms, 4917G independently predicts the presence of AMD (OR = 2.16, 95%CI 1.20-3.91, p = 0.01). In conclusion, a specific mitochondrial polymorphism previously implicated in other neurodegenerative phenotypes (4917G) appears to convey risk for AMD independent of recently discovered nuclear DNA polymorphisms.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Impact of economic, regulatory, and patent policies on innovation in cancer chemoprevention.

Chemoprevention agents are an emerging new scientific area that holds out the promise of delaying or avoiding a number of common cancers. These new agents face significant scientific, regulatory, and economic barriers, however, which have limited investment in their research and development (R&D). These barriers include above-average clinical trial scales, lengthy time frames between discovery and Food and Drug Administration approval, liability risks (because they are given to healthy individuals), and a growing funding gap for early-stage candidates. The longer time frames and risks associated with chemoprevention also cause exclusivity time on core patents to be limited or subject to significant uncertainties. We conclude that chemoprevention uniquely challenges the structure of incentives embodied in the economic, regulatory, and patent policies for the biopharmaceutical industry. Many of these policy issues are illustrated by the recently Food and Drug Administration-approved preventive agents Gardasil and raloxifene. Our recommendations to increase R&D investment in chemoprevention agents include (a) increased data exclusivity times on new biological and chemical drugs to compensate for longer gestation periods and increasing R&D costs; chemoprevention is at the far end of the distribution in this regard; (b) policies such as early-stage research grants and clinical development tax credits targeted specifically to chemoprevention agents (these are policies that have been very successful in increasing R&D investment for orphan drugs); and (c) a no-fault liability insurance program like that currently in place for children's vaccines.

Using the institutional grammar tool to understand regulatory compliance: The case of Colorado aquaculture

What is the relationship between the design of regulations and levels of individual compliance? To answer this question, Crawford and Ostrom's institutional grammar tool is used to deconstruct regulations governing the aquaculture industry in Colorado, USA. Compliance with the deconstructed regulatory components is then assessed based on the perceptions of the appropriateness of the regulations, involvement in designing the regulations, and intrinsic and extrinsic motivations. The findings suggest that levels of compliance with regulations vary across and within individuals regarding various aspects of the regulatory components. As expected, the level of compliance is affected by the perceived appropriateness of regulations, participation in designing the regulations, and feelings of guilt and fear of social disapproval. Furthermore, there is a strong degree of interdependence among the written components, as identified by the institutional grammar tool, in affecting compliance levels. The paper contributes to the regulation and compliance literature by illustrating the utility of the institutional grammar tool in understanding regulatory content, applying a new Q-Sort technique for measuring individual levels of compliance, and providing a rare exploration into feelings of guilt and fear outside of the laboratory setting. © 2012 Blackwell Publishing Asia Pty Ltd.

A C-terminal motif found in the beta2-adrenergic receptor, P2Y1 receptor and cystic fibrosis transmembrane conductance regulator determines binding to the Na+/H+ exchanger regulatory factor family of PDZ proteins.

The Na+/H+ exchanger regulatory factor (NHERF) binds to the tail of the beta2-adrenergic receptor and plays a role in adrenergic regulation of Na+/H+ exchange. NHERF contains two PDZ domains, the first of which is required for its interaction with the beta2 receptor. Mutagenesis studies of the beta2 receptor tail revealed that the optimal C-terminal motif for binding to the first PDZ domain of NHERF is D-S/T-x-L, a motif distinct from those recognized by other PDZ domains. The first PDZ domain of NHERF-2, a protein that is 52% identical to NHERF and also known as E3KARP, SIP-1, and TKA-1, exhibits binding preferences very similar to those of the first PDZ domain of NHERF. The delineation of the preferred binding motif for the first PDZ domain of the NHERF family of proteins allows for predictions for other proteins that may interact with NHERF or NHERF-2. For example, as would be predicted from the beta2 receptor tail mutagenesis studies, NHERF binds to the tail of the purinergic P2Y1 receptor, a seven-transmembrane receptor with an intracellular C-terminal tail ending in D-T-S-L. NHERF also binds to the tail of the cystic fibrosis transmembrane conductance regulator, which ends in D-T-R-L. Because the preferred binding motif of the first PDZ domain of the NHERF family of proteins is found at the C termini of a variety of intracellular proteins, NHERF and NHERF-2 may be multifunctional adaptor proteins involved in many previously unsuspected aspects of intracellular signaling.

Impaired formation of beta-adrenergic receptor-nucleotide regulatory protein complexes in pseudohypoparathyroidism.

Decreased activity of the guanine nucleotide regulatory protein (N) of the adenylate cyclase system is present in cell membranes of some patients with pseudohypoparathyrodism (PHP-Ia) whereas others have normal activity of N (PHP-Ib). Low N activity in PHP-Ia results in a decrease in hormone (H)-stimulatable adenylate cyclase in various tissues, which might be due to decreased ability to form an agonist-specific high affinity complex composed of H, receptor (R), and N. To test this hypothesis, we compared beta-adrenergic agonist-specific binding properties in erythrocyte membranes from five patients with PHP-Ia (N = 45% of control), five patients with PHP-Ib (N = 97%), and five control subjects. Competition curves that were generated by increasing concentrations of the beta-agonist isoproterenol competing with [125I]pindolol were shallow (slope factors less than 1) and were computer fit to a two-state model with corresponding high and low affinity for the agonist. The agonist competition curves from the PHP-Ia patients were shifted significantly (P less than 0.02) to the right as a result of a significant (P less than 0.01) decrease in the percent of beta-adrenergic receptors in the high affinity state from 64 +/- 22% in PHP-Ib and 56 +/- 5% in controls to 10 +/- 8% in PHP-Ia. The agonist competition curves were computer fit to a "ternary complex" model for the two-step reaction: H + R + N in equilibrium HR + N in equilibrium HRN. The modeling was consistent with a 60% decrease in the functional concentration of N, and was in good agreement with the biochemically determined decrease in erythrocyte N protein activity. These in vitro findings in erythrocytes taken together with the recent observations that in vivo isoproterenol-stimulated adenylate cyclase activity is decreased in patients with PHP (Carlson, H. E., and A. S. Brickman, 1983, J. Clin. Endocrinol. Metab. 56:1323-1326) are consistent with the notion that N is a bifunctional protein interacting with both R and the adenylate cyclase. It may be that in patients with PHP-Ia a single molecular and genetic defect accounts for both decreased HRN formation and decreased adenylate cyclase activity, whereas in PHP-Ib the biochemical lesion(s) appear not to affect HRN complex formation.

A functional polymorphism in the 5HTR2C gene associated with stress responses also predicts incident cardiovascular events.

Previously we have shown that a functional nonsynonymous single nucleotide polymorphism (rs6318) of the 5HTR2C gene located on the X-chromosome is associated with hypothalamic-pituitary-adrenal axis response to a stress recall task, and with endophenotypes associated with cardiovascular disease (CVD). These findings suggest that individuals carrying the rs6318 Ser23 C allele will be at higher risk for CVD compared to Cys23 G allele carriers. The present study examined allelic variation in rs6318 as a predictor of coronary artery disease (CAD) severity and a composite endpoint of all-cause mortality or myocardial infarction (MI) among Caucasian participants consecutively recruited through the cardiac catheterization laboratory at Duke University Hospital (Durham, NC) as part of the CATHGEN biorepository. Study population consisted of 6,126 Caucasian participants (4,036 [65.9%] males and 2,090 [34.1%] females). A total of 1,769 events occurred (1,544 deaths and 225 MIs; median follow-up time = 5.3 years, interquartile range = 3.3-8.2). Unadjusted Cox time-to-event regression models showed, compared to Cys23 G carriers, males hemizygous for Ser23 C and females homozygous for Ser23C were at increased risk for the composite endpoint of all-cause death or MI: Hazard Ratio (HR) = 1.47, 95% confidence interval (CI) = 1.17, 1.84, p = .0008. Adjusting for age, rs6318 genotype was not related to body mass index, diabetes, hypertension, dyslipidemia, smoking history, number of diseased coronary arteries, or left ventricular ejection fraction in either males or females. After adjustment for these covariates the estimate for the two Ser23 C groups was modestly attenuated, but remained statistically significant: HR = 1.38, 95% CI = 1.10, 1.73, p = .005. These findings suggest that this functional polymorphism of the 5HTR2C gene is associated with increased risk for CVD mortality and morbidity, but this association is apparently not explained by the association of rs6318 with traditional risk factors or conventional markers of atherosclerotic disease.

Functional annotation signatures of disease susceptibility loci improve SNP association analysis.

BACKGROUND: Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study's analysis plan. RESULTS: We developed a Bayesian statistical model for the prior probability of phenotype-genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super-track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen-2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non-informative predictors and evaluated the model's ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP's presence in the GC. Further, using data from a genome-wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome-wide scale and improves power to detect associations. CONCLUSIONS: We show how diverse functional annotations can be efficiently combined to create 'functional signatures' that predict the a priori odds of a variant's association to a trait and how these signatures can be integrated into a standard genome-wide-scale association analysis, resulting in improved power to detect truly associated variants.

Pharmacogenomics in early-phase clinical development.

Pharmacogenomics (PGx) offers the promise of utilizing genetic fingerprints to predict individual responses to drugs in terms of safety, efficacy and pharmacokinetics. Early-phase clinical trial PGx applications can identify human genome variations that are meaningful to study design, selection of participants, allocation of resources and clinical research ethics. Results can inform later-phase study design and pipeline developmental decisions. Nevertheless, our review of the clinicaltrials.gov database demonstrates that PGx is rarely used by drug developers. Of the total 323 trials that included PGx as an outcome, 80% have been conducted by academic institutions after initial regulatory approval. Barriers for the application of PGx are discussed. We propose a framework for the role of PGx in early-phase drug development and recommend PGx be universally considered in study design, result interpretation and hypothesis generation for later-phase studies, but PGx results from underpowered studies should not be used by themselves to terminate drug-development programs.

Cutting edge: distinct glycolytic and lipid oxidative metabolic programs are essential for effector and regulatory CD4+ T cell subsets.

Stimulated CD4(+) T lymphocytes can differentiate into effector T cell (Teff) or inducible regulatory T cell (Treg) subsets with specific immunological roles. We show that Teff and Treg require distinct metabolic programs to support these functions. Th1, Th2, and Th17 cells expressed high surface levels of the glucose transporter Glut1 and were highly glycolytic. Treg, in contrast, expressed low levels of Glut1 and had high lipid oxidation rates. Consistent with glycolysis and lipid oxidation promoting Teff and Treg, respectively, Teff were selectively increased in Glut1 transgenic mice and reliant on glucose metabolism, whereas Treg had activated AMP-activated protein kinase and were dependent on lipid oxidation. Importantly, AMP-activated protein kinase stimulation was sufficient to decrease Glut1 and increase Treg generation in an asthma model. These data demonstrate that CD4(+) T cell subsets require distinct metabolic programs that can be manipulated in vivo to control Treg and Teff development in inflammatory diseases.

Comparative analyses of seven algorithms for copy number variant identification from single nucleotide polymorphism arrays.

Determination of copy number variants (CNVs) inferred in genome wide single nucleotide polymorphism arrays has shown increasing utility in genetic variant disease associations. Several CNV detection methods are available, but differences in CNV call thresholds and characteristics exist. We evaluated the relative performance of seven methods: circular binary segmentation, CNVFinder, cnvPartition, gain and loss of DNA, Nexus algorithms, PennCNV and QuantiSNP. Tested data included real and simulated Illumina HumHap 550 data from the Singapore cohort study of the risk factors for Myopia (SCORM) and simulated data from Affymetrix 6.0 and platform-independent distributions. The normalized singleton ratio (NSR) is proposed as a metric for parameter optimization before enacting full analysis. We used 10 SCORM samples for optimizing parameter settings for each method and then evaluated method performance at optimal parameters using 100 SCORM samples. The statistical power, false positive rates, and receiver operating characteristic (ROC) curve residuals were evaluated by simulation studies. Optimal parameters, as determined by NSR and ROC curve residuals, were consistent across datasets. QuantiSNP outperformed other methods based on ROC curve residuals over most datasets. Nexus Rank and SNPRank have low specificity and high power. Nexus Rank calls oversized CNVs. PennCNV detects one of the fewest numbers of CNVs.

Rapamycin Interferes With Postdepletion Regulatory T Cell Homeostasis and Enhances DSA Formation Corrected by CTLA4-Ig.

Previously, we demonstrated that alemtuzumab induction with rapamycin as sole maintenance therapy is associated with an increased incidence of humoral rejection in human kidney transplant patients. To investigate the role of rapamycin in posttransplant humoral responses after T cell depletion, fully MHC mismatched hearts were transplanted into hCD52Tg mice, followed by alemtuzumab treatment with or without a short course of rapamycin. While untreated hCD52Tg recipients acutely rejected B6 hearts (n = 12), hCD52Tg recipients treated with alemtuzumab alone or in conjunction with rapamycin showed a lack of acute rejection (MST > 100). However, additional rapamycin showed a reduced beating quality over time and increased incidence of vasculopathy. Furthermore, rapamycin supplementation showed an increased serum donor-specific antibodies (DSA) level compared to alemtuzumab alone at postoperation days 50 and 100. Surprisingly, additional rapamycin treatment significantly reduced CD4(+) CD25(+) FoxP3(+) T reg cell numbers during treatment. On the contrary, ICOS(+) PD-1(+) CD4 follicular helper T cells in the lymph nodes were significantly increased. Interestingly, CTLA4-Ig supplementation in conjunction with rapamycin corrected rapamycin-induced accelerated posttransplant humoral response by directly modulating Tfh cells but not Treg cells. This suggests that rapamycin after T cell depletion could affect Treg cells leading to an increase of Tfh cells and DSA production that can be reversed by CTLA4-Ig.