Two-marker association tests yield new disease associations for coronary artery disease and hypertension

It has been postulated that multiple-marker methods may have added ability, over single-marker methods, to detect genetic variants associated with disease. The Wellcome Trust Case Control Consortium (WTCCC) provided the first successful large genome-wide association studies (GWAS) which included single-marker association analyses for seven common complex diseases. Of those signals detected, only one was associated with coronary artery disease (CAD), and none were identified for hypertension (HTN). Our objective was to find additional genetic associations and pathways for cardiovascular disease by examining the WTCCC data for variants associated with CAD and HTN using two-marker testing methods. We applied two-marker association testing to the WTCCC dataset, which includes ~2,000 affected individuals with each disorder, and a shared pool of ~3,000 controls, all genotyped using Affymetrix GeneChip 500 K arrays. For CAD, we detected single nucleotide polymorphisms (SNP) pairs in three genes showing genome-wide significance: HFE2, STK32B, and DIPC2. The most notable SNP pairs in a non-protein-coding region were at 9p21, a known major CAD-associated region. For HTN, we detected SNP pairs in five genes: GPR39, XRCC4, MYO6, ZFAT, and MACROD2. Four further associated SNP pair regions were at least 70 kb from any known gene. We have shown that novel, multiple-marker, statistical methods can be of use in finding variants in GWAS. We describe many new, associated variants for both CAD and HTN and describe their known genetic mechanisms.     

Identifying when weather influences life-history traits of grazing herbivores

1. There is increasing evidence that density-independent weather effects influence life-history traits and hence the dynamics of populations of animals. Here, we present a novel statistical approach to estimate when such influences are strongest. The method is demonstrated by analyses investigating the timing of the influence of weather on the birth weight of sheep and deer. 2. The statistical technique allowed for the pattern of temporal correlation in the weather data enabling the effects of weather in many fine-scale time intervals to be investigated simultaneously. Thus, while previous studies have typically considered weather averaged across a single broad time interval during pregnancy, our approach enabled examination simultaneously of the relationships with weekly and fortnightly averages throughout the whole of pregnancy. 3. We detected a positive effect of temperature on the birth weight of deer, which is strongest in late pregnancy (mid-March to mid-April), and a negative effect of rainfall on the birthweight of sheep, which is strongest during mid-pregnancy (late January to early February). The possible mechanisms underlying these weather-birth weight relationships are discussed. 4. This study enhances our insight into the pattern of the timing of influence of weather on early development. The method is of much more general application and could provide valuable insights in other areas of ecology in which sequences of intercorrelated explanatory variables have been collected in space or in time.     

Improving power in contrasting linkage-disequilibrium patterns between cases and controls

Genetic association studies offer an opportunity to find genetic variants underlying complex human diseases. The success of this approach depends on the linkage disequilibrium (LD) between markers and the disease variant(s) in a local region of the genome. Because, in the region with a disease mutation, the LD pattern among markers may differ between cases and controls, in some scenarios, it is useful to compare a measure of this LD, to map disease mutations. For example, using the composite correlation to characterize the LD among markers, Zaykin et al. recently suggested an "LD contrast" test and showed that it has high power under certain haplotype-driven disease models. Furthermore, it is likely that individual variants observed at different positions in a gene act jointly with each other to influence the phenotype, and the LD contrast test is also a useful method to detect such joint action. However, the LD among markers introduced by mutations and their joint action is usually confounded by background LD, which is measured at the population level, especially in a local region with disease mutations. Because the measures of LD that are usually used, such as the composite correlation, represent both effects, they may not be optimal for the purpose of detecting association when high background LD exists. Here, we describe a test that improves the LD contrast test by taking into account the background LD. Because the proposed test is developed in a regression framework, it is very flexible and can be extended to continuous traits and to incorporate covariates. Our simulation results demonstrate the validity and substantially higher power of the proposed method over current methods. Finally, we illustrate our new method by applying it to real data from the International Collaborative Study on Hypertension in Blacks.     

CellGeo: A computational platform for the analysis of shape changes in cells with complex geometries

Cell biologists increasingly rely on computer-aided image analysis, allowing them to collect precise, unbiased quantitative results. However, despite great progress in image processing and computer vision, current computational approaches fail to address many key aspects of cell behavior, including the cell protrusions that guide cell migration and drive morphogenesis. We developed the open source MATLAB application CellGeo, a user-friendly computational platform to allow simultaneous, automated tracking and analysis of dynamic changes in cell shape, including protrusions ranging from filopodia to lamellipodia. Our method maps an arbitrary cell shape onto a tree graph that, unlike traditional skeletonization algorithms, preserves complex boundary features. CellGeo allows rigorous but flexible definition and accurate automated detection and tracking of geometric features of interest. We demonstrate CellGeo’s utility by deriving new insights into (a) the roles of Diaphanous, Enabled, and Capping protein in regulating filopodia and lamellipodia dynamics in Drosophila melanogaster cells and (b) the dynamic properties of growth cones in catecholaminergic a–differentiated neuroblastoma cells.     

Patients with RET D631Y mutations most commonly present with pheochromocytoma and not medullary thyroid carcinoma

Multiple endocrine neoplasia type 2a results from an activating germline mutation in the RET proto-oncogene. Carriers of a RET mutation are at risk of medullary thyroid carcinoma, pheochromocytoma, and primary hyperparathyroidism. Most individuals with multiple endocrine neoplasia type 2a eventually develop medullary thyroid carcinoma and as there is a strong genotype-phenotype correlation, guidelines have been established as to the age recommended for prophylactic thyroidectomy. However for rare mutations in the RET proto-oncogene there is insufficient evidence to provide guidance as to the risk of medullary thyroid carcinoma. We present a family with the rare RET mutation, D631Y in which the proband initially presented with a pheochromocytoma, and review the available literature pertaining to this mutation. In 83% of index cases, pheochromocytoma was the presenting feature and only 37% of adult germline mutation carriers have developed medullary thyroid carcinoma, none of whom have been reported to have nodal or metastatic disease. Patients with a D631Y RET mutation typically present with pheochromocytoma and medullary thyroid carcinoma appears to occur with a later onset than reported with other RET mutations. Based on the current literature we recommend performing prophylactic total thyroidectomy by age 12 years for D631Y carriers although this recommendation may need to be reviewed as additional data becomes available.     

A novel stochastic simulation approach enables exploration of mechanisms for regulating polarity site movement

Cells polarize their movement or growth toward external directional cues in many different contexts. For example, budding yeast cells grow toward potential mating partners in response to pheromone gradients. Directed growth is controlled by polarity factors that assemble into clusters at the cell membrane. The clusters assemble, disassemble, and move between different regions of the membrane before eventually forming a stable polarity site directed toward the pheromone source. Pathways that regulate clustering have been identified but the molecular mechanisms that regulate cluster mobility are not well understood. To gain insight into the contribution of chemical noise to cluster behavior we simulated clustering using the reaction-diffusion master equation (RDME) framework to account for molecular-level fluctuations. RDME simulations are a computationally efficient approximation, but their results can diverge from the underlying microscopic dynamics. We implemented novel concentration-dependent rate constants that improved the accuracy of RDME-based simulations, allowing us to efficiently investigate how cluster dynamics might be regulated. Molecular noise was effective in relocating clusters when the clusters contained low numbers of limiting polarity factors, and when Cdc42, the central polarity regulator, exhibited short dwell times at the polarity site. Cluster stabilization occurred when abundances or binding rates were altered to either lengthen dwell times or increase the number of polarity molecules in the cluster. We validated key results using full 3D particle-based simulations. Understanding the mechanisms cells use to regulate the dynamics of polarity clusters should provide insights into how cells dynamically track external directional cues.     

Mathematical model of nucleotide regulation on airway epithelia. Implications for airway homeostasis

In the airways, adenine nucleotides support a complex signaling network mediating host defenses. Released by the epithelium into the airway surface liquid (ASL) layer, they regulate mucus clearance through P2 (ATP) receptors, and following surface metabolism through P1 (adenosine; Ado) receptors. The complexity of ASL nucleotide regulation provides an ideal subject for biochemical network modeling. A mathematical model was developed to integrate nucleotide release, the ectoenzymes supporting the dephosphorylation of ATP into Ado, Ado deamination into inosine (Ino), and nucleoside uptake. The model also includes ecto-adenylate kinase activity and feed-forward inhibition of Ado production by ATP and ADP. The parameters were optimized by fitting the model to experimental data for the steady-state and transient concentration profiles generated by adding ATP to polarized primary cultures of human bronchial epithelial (HBE) cells. The model captures major aspects of ATP and Ado regulation, including their >4-fold increase in concentration induced by mechanical stress mimicking normal breathing. The model also confirmed the independence of steady-state nucleotide concentrations on the ASL volume, an important regulator of airway clearance. An interactive approach between simulations and assays revealed that feed-forward inhibition is mediated by selective inhibition of ecto-5'-nucleotidase. Importantly, the model identifies ecto-adenylate kinase as a key regulator of ASL ATP and proposes novel strategies for the treatment of airway diseases characterized by impaired nucleotide-mediated clearance. These new insights into the biochemical processes supporting ASL nucleotide regulation illustrate the potential of this mathematical model for fundamental and clinical research.     

Inhibition of papillomavirus protein function in cervical cancer cells by intrabody targeting

Papillomaviruses (HPVs) are a major cause of human disease, and are responsible for approximately half a million cases of cervical cancer each year. HPVs also cause genital warts, and are the most common sexually transmitted disease in many countries. Despite their importance, there are currently no specific antivirals that are active against HPVs. Papillomavirus protein function is mediated largely by protein-protein interactions, which are difficult to inhibit using conventional approaches. To circumvent these problems, we have prepared an scFv library, and have used this to isolate high-affinity binding molecules that may stearically hinder the association of E6 with p53 and prevent E6-mediated p53 degradation in cervical cancer cells. One of the molecules isolated from the library (GTE6-1), had an affinity for 16E6 of 60nM, and bound within the first zinc finger of the protein. GTE6-1 was able to associate with non-denatured E6 following expression in mammalian cells and could inhibit E6-mediated p53 degradation in in vitro assays. E6-mediated p53 degradation is essential for the continuous growth of cervical cancer cells caused by HPV16. To examine the potential of GTE6-1 as an inhibitor of E6 function in such cells, the molecule was expressed in scFv, diabody and triabody formats in a number of cell lines that are driven to proliferate by the HPV16 oncogenes E6 and E7, including the cervical cancer cell line SiHa. In contrast to small E6-binding peptides containing the ELLG E6-binding motif, GTE6-1 expression lead to changes in nuclear structure, the appearance of apoptosis markers, and an elevation in the levels of p53. No effects were seen with a control scFv molecule, or when GTE6-1 was expressed in cells that are driven to proliferate by simian virus 40 (SV40) T-antigen. Given the accessibility of HPV-associated lesions to topical therapy, our results suggest that large interfering molecules such as intrabodies may be useful inhibitors of viral protein-protein interactions and be particularly appropriate for the treatment of HPV-associated disease.     

Expression of AMPA receptor subunits at synapses in laminae I–III of the rodent spinal dorsal horn

Neuropathic pain may arise following peripheral nerve injury though the molecular mechanisms associated with this are unclear. We used proteomic profiling to examine changes in protein expression associated with the formation of hyper-excitable neuromas derived from rodent saphenous nerves. A two-dimensional difference gel electrophoresis (2D-DIGE) profiling strategy was employed to examine protein expression changes between developing neuromas and normal nerves in whole tissue lysates. We found around 200 proteins which displayed a >1.75-fold change in expression between neuroma and normal nerve and identified 55 of these proteins using mass spectrometry. We also used immunoblotting to examine the expression of low-abundance ion channels Nav1.3, Nav1.8 and calcium channel α2δ-1 subunit in this model, since they have previously been implicated in neuronal hyperexcitability associated with neuropathic pain. Finally, S³⁵methionine in vitro labelling of neuroma and control samples was used to demonstrate local protein synthesis of neuron-specific genes. A number of cytoskeletal proteins, enzymes and proteins associated with oxidative stress were up-regulated in neuromas, whilst overall levels of voltage-gated ion channel proteins were unaffected. We conclude that altered mRNA levels reported in the somata of damaged DRG neurons do not necessarily reflect levels of altered proteins in hyper-excitable damaged nerve endings. An altered repertoire of protein expression, local protein synthesis and topological re-arrangements of ion channels may all play important roles in neuroma hyper-excitability.