Combining elicited expert knowledge into Bayesian priors

Information that is elicited from experts can be treated as `data', so can be analysed using a Bayesian statistical model, to formulate a prior model. Typically methods for encoding a single expert's knowledge have been parametric, constrained by the extent of an expert's knowledge and energy regarding a target parameter. Interestingly these methods have often been deterministic, in that all elicited information is treated at `face value', without error. Here we sought a parametric and statistical approach for encoding assessments from multiple experts. Our recent work proposed and demonstrated the use of a flexible hierarchical model for this purpose. In contrast to previous mathematical approaches like linear or geometric pooling, our new approach accounts for several sources of variation: elicitation error, encoding error and expert diversity. Of interest are the practical, mathematical and philosophical interpretations of this form of hierarchical pooling (which is both statistical and parametric), and how it fits within the subjective Bayesian paradigm. Case studies from a bioassay and project management (on PhDs) are used to illustrate the approach.

In Plant Activation: An inducible, hyperexpression platform for recombinant protein production in plants

In this study, we describe a novel protein production platform that provides both activation and amplification of transgene expression in planta. The In Plant Activation (INPACT) system is based on the replication machinery of tobacco yellow dwarf mastrevirus (TYDV) and is essentially transient gene expression from a stably transformed plant, thus combining the advantages of both means of expression. The INPACT cassette is uniquely arranged such that the gene of interest is split and only reconstituted in the presence of the TYDV-encoded Rep/RepA proteins. Rep/RepA expression is placed under the control of the AlcA:AlcR gene switch, which is responsive to trace levels of ethanol. Transgenic tobacco (Nicotiana tabacum cv Samsun) plants containing an INPACT cassette encoding the b-glucuronidase (GUS) reporter had negligible background expression but accumulated very high GUS levels (up to 10% total soluble protein) throughout the plant, within 3 d of a 1% ethanol application. The GUS reporter was replaced with a gene encoding a lethal ribonuclease, barnase, demonstrating that the INPACT system provides exquisite control of transgene expression and can be adapted to potentially toxic or inhibitory compounds. The INPACT gene expression platform is scalable, not host-limited, and has been used to express both a therapeutic and an industrial protein.

High-resolution loss of heterozygosity screening implicates PTPRJ as a potential tumor suppressor gene that affects susceptibility to non-Hodgkin's Lymphoma

We employed a Hidden-Markov-Model (HMM) algorithm in loss of heterozygosity (LOH) analysis of high-density single nucleotide polymorphism (SNP) array data from Non-Hodgkin’s lymphoma (NHL) entities, follicular lymphoma (FL), and diffuse large B-cell lymphoma (DLBCL). This revealed a high frequency of LOH over the chromosomal region 11p11.2, containing the gene encoding the protein tyrosine phosphatase receptor type J (PTPRJ). Although PTPRJ regulates components of key survival pathways in B-cells (i.e., BCR, MAPK, and PI3K signaling), its role in B-cell development is poorly understood. LOH of PTPRJ has been described in several types of cancer but not in any hematological malignancy. Interestingly, FL cases with LOH exhibited down-regulation of PTPRJ, in contrast no significant variation of expression was shown in DLBCLs. In addition, sequence screening in Exons 5 and 13 of PTPRJ identified the G973A (rs2270993), T1054C (rs2270992), A1182C (rs1566734), and G2971C (rs4752904) coding SNPs (cSNPs). The A1182 allele was significantly more frequent in FLs and in NHLs with LOH. Significant over-representation of the C1054 (rs2270992) and the C2971 (rs4752904) alleles were also observed in LOH cases. A haplotype analysis also revealed a significant lower frequency of haplotype GTCG in NHL cases, but it was only detected in cases with retention. Conversely, haplotype GCAC was over-representated in cases with LOH. Altogether, these results indicate that the inactivation of PTPRJ may be a common lymphomagenic mechanism in these NHL subtypes and that haplotypes in PTPRJ gene may play a role in susceptibility to NHL, by affecting activation of PTPRJ in these B-cell lymphomas.

Epistatic effects of potassium channel variation on cardiac repolarization and atrial fibrillation risk

Objectives The aim of this study was to evaluate the role of cardiac K+ channel gene variants in families with atrial fibrillation (AF). Background The K+ channels play a major role in atrial repolarization but single mutations in cardiac K+ channel genes are infrequently present in AF families. The collective effect of background K+ channel variants of varying prevalence and effect size on the atrial substrate for AF is largely unexplored. Methods Genes encoding the major cardiac K+ channels were resequenced in 80 AF probands. Nonsynonymous coding sequence variants identified in AF probands were evaluated in 240 control subjects. Novel variants were characterized using patch-clamp techniques and in silico modeling was performed using the Courtemanche atrial cell model. Results Nineteen nonsynonymous variants in 9 genes were found, including 11 rare variants. Rare variants were more frequent in AF probands (18.8% vs. 4.2%, p < 0.001), and the mean number of variants was greater (0.21 vs. 0.04, p < 0.001). The majority of K+ channel variants individually had modest functional effects. Modeling simulations to evaluate combinations of K+ channel variants of varying population frequency indicated that simultaneous small perturbations of multiple current densities had nonlinear interactions and could result in substantial (>30 ms) shortening or lengthening of action potential duration as well as increased dispersion of repolarization. Conclusions Families with AF show an excess of rare functional K+ channel gene variants of varying phenotypic effect size that may contribute to an atrial arrhythmogenic substrate. Atrial cell modeling is a useful tool to assess epistatic interactions between multiple variants.

Functional analysis of missense variants in the TRESK (KCNK18) K+ channel

A loss of function mutation in the TRESK K2P potassium channel (KCNK18), has recently been linked with typical familial migraine with aura. We now report the functional characterisation of additional TRESK channel missense variants identified in unrelated patients. Several variants either had no apparent functional effect, or they caused a reduction in channel activity. However, the C110R variant was found to cause a complete loss of TRESK function, yet is present in both sporadic migraine and control cohorts, and no variation in KCNK18 copy number was found. Thus despite the previously identified association between loss of TRESK channel activity and migraine in a large multigenerational pedigree, this finding indicates that a single non-functional TRESK variant is not alone sufficient to cause typical migraine and highlights the genetic complexity of this disorder. Migraine is a common, disabling neurological disorder with a genetic, environmental and in some cases hormonal component. It is characterized by attacks of severe, usually unilateral and throbbing headache, can be accompanied by nausea, vomiting and photophobia and is clinically divided into two main subtypes, migraine with aura (MA) when a migraine is accompanied by transient and reversible focal neurological symptoms and migraine without aura (MO)1. The multifactorial and clinical heterogeneity of the disorder have considerably hindered the identification of common migraine susceptibility genes and most of our current understanding comes from the studies of familial hemiplegic migraine (FHM), a rare monogenic autosomal dominant form of MA2. So far, the three susceptibility genes that have been convincingly identified in FHM families all encode ion channels or transporters: CACNA1A encoding the α1 subunit of the Cav2.1 calcium channel3, SCN1A encoding the Nav1.1 sodium channel4 and ATP1A2 encoding the α2 subunit of the Na+/K+ pump5. It is believed that mutations in these genes may lead to increased efflux of glutamate and potassium in the synapse and thereby cause migraine by rendering the brain more susceptible to cortical spreading depression (CSD)6 which is thought to play a role in initiating a migraine attack7,8. However, these genes have not to date been implicated in common forms of migraine9. Nevertheless, current opinion suggests that typical migraine, like FHM, is also disorder of neuronal excitability, ion homeostasis and neurotransmitter release10,11,12. Mutations in the SLC4A4 gene encoding the sodium-bicarbonate cotransporter NBCe1, have recently been implicated in several different forms of migraine13, and a variety of genes involved in glutamate homeostasis (PGCP, MTDH14 and LRP115) and a cation channel (TRPM8)15 have also recently been implicated in migraine via genome-wide association studies. Ion channels are therefore highly likely to play an important role in the pathogenesis of typical migraine. TRESK (KCNK18), is a member of the two-pore domain (K2P) family of potassium channels involved in the control of cellular electrical excitability16. Regulation of TRESK activity by the calcium-dependent phosphatase calcineurin17, as well as its expression in dorsal root ganglia (DRG)18 and trigeminal ganglia (TG)19,20 has led to a proposed role for this channel in a variety of pain pathways. In a recent study, a frameshift mutation (F139Wfsx24) in TRESK was identified in a large multigenerational pedigree where it co-segregated perfectly with typical MA and a significant genome-wide linkage LOD score of 3.0. Furthermore, functional analysis revealed that this mutation caused a complete loss of TRESK function and that the truncated subunit was also capable of down regulating wild-type channel function. This therefore highlighted KCNK18 as potentially important candidate gene and suggested that TRESK dysfunction might play a possible role in the pathogenesis of familial migraine with visual aura20. Additional screening for KCNK18 mutations in unrelated sporadic migraine and control cohorts also identified a number of other missense variants; R10G, A34V, C110R, S231P and A233V20. The A233V variant was found only in the control cohort, whilst A34V was identified in a single Australian migraine proband for which family samples were not available, but it was not detected in controls. By contrast, the R10G, C110R, and S231P variants were found in both migraineurs and controls in both cohorts. In this study, we have investigated the functional effect of these variants to further probe the potential association of TRESK dysfunction with typical migraine.

Linkage disequilibrium analysis in the genetically isolated Norfolk Island population

Norfolk Island is a human genetic isolate, possessing unique population characteristics that could be utilized for complex disease gene localization. Our intention was to evaluate the extent and strength of linkage disequilibrium (LD) in the Norfolk isolate by investigating markers within Xq13.3 and the NOS2A gene encoding the inducible nitric oxide synthase. A total of six microsatellite markers spanning approximately 11 Mb were assessed on chromosome Xq13.3 in a group of 56 men from Norfolk Island. Additionally, three single nucleotide polymorphisms (SNPs) localizing to the NOS2A gene were analyzed in a subset of the complex Norfolk pedigree. With the exception of two of the marker pairs, one of which is the most distantly spaced marker, all the Xq13.3 marker pairs were found to be in significant LD indicating that LD extends up to 9.5-11.5 Mb in the Norfolk Island population. Also, all SNPs studied showed significant LD in both Norfolk Islanders and Australian Caucasians, with two of the marker pairs in complete LD in the Norfolk population only. The Norfolk Island study population possesses a unique set of characteristics including founder effect, geographical isolation, exhaustive genealogical information and phenotypic data of use to cardiovascular disease risk traits. With LD extending up to 9.5-11 Mb, the Norfolk isolate should be a powerful resource for the localization of complex disease genes.

Association of the NuMA region on chromosome 11q13 with breast cancer susceptibility

The development of breast cancer is a complex process that involves multiple genes at many stages, from initial cell cycle dysregulation to disease progression. To identify genetic variations that influence this process, we conducted a large-scale association study using a collection of German cases and controls and >25,000 SNPs located within 16,000 genes. One of the loci identified was located on chromosome 11q13 [odds ratio (OR)=1.85, P=0.017]. The initial association was subsequently tested in two independent breast cancer collections. In both sample sets, the frequency of the susceptibility allele was increased in the cases (OR=1.6, P=0.01). The susceptibility allele was also associated with an increase in cancer family history (P=0.1). Fine mapping showed that the region of association extends approximately 300 kb and spans several genes, including the gene encoding the nuclear mitotic apparatus protein (NuMA). A nonsynonymous SNP (A794G) in NuMA was identified that showed a stronger association with breast cancer risk than the initial marker SNP (OR=2.8, P=0.005 initial sample; OR=2.1, P=0.002 combined). NuMA is a cell cycle-related protein essential for normal mitosis that is degraded in early apoptosis. NuMA-retinoic acid receptor alpha fusion proteins have been described in acute promyelocytic leukemia. Although the potential functional relevance of the A794G variation requires further biological validation, we conclude that variations in the NuMA gene are likely responsible for the observed increased breast cancer risk.

The null allele of GSTM1 does not affect susceptibility to solar keratoses in the Australian white population

Solar keratoses (SKs) are induced by exposure to UV radiation and are capable of undergoing transformation to squamous cell carcinoma (SCC).1 The two main factors influencing the occurrence of SK are the sensitivity of the skin to sunlight and the total duration of solar exposure. These factors are responsible for the high incidence of SK in Australia. Although the influence of genetic factors is not defined, there is evidence that the gene encoding the enzyme, glutathione S-transferase, may be implicated in cancer predisposition and therefore SK. Glutathione S-transferase Mu-1 (GSTM1) is an isoenzyme involved in the detoxification of carcinogens. The GSTM1 protein is completely absent in approximately 50% of white persons. This absence is caused by a homozygous gene deletion on chromosome 1p resulting in a null genotype.2 Katoh3 showed that the frequency of the GSTM1 null genotype was significantly higher in 85 patients with urothelial cancer (61.2%; p < 0.05), suggesting that the null genotype may increase cancer susceptibility. This finding was supported by Lafuente et al.4 who found evidence that persons who lack the GSTM1 gene have approximately twice the chance of experiencing malignant melanoma. Further research in the United Kingdom found that patients with two or more skin tumors of different types, basal cell carcinoma (BCC) and SCC, had a significantly higher frequency of GSTM1 null genotypes than controls (71%; p = 0.033). However the GSTM1 genotype in patients with only SCC was not excessive in this population.5 Persons residing in northern Australia have the highest incidence of nonmelanoma skin cancer (SCC and BCC) in the world6 and receive far greater solar exposure than persons residing in the United Kingdom. It is possible that the GSTM1 null genotype may affect susceptibility to SK, which may act as SCC precursors, in Australians exposed to these high levels of solar radiation.

The importance of individual difference : examining anthropomorphic tendency and responses to spokes-characters

Anthropomorphism is a cognitive bias, which occurs when individuals see human characteristics in a non-human agent, object or animal. Anthropomorphism is especially interesting to marketers, because once anthropomorphic bias has been triggered, it can lead to a greater feeling of connectedness to a non-human agent (Tam, Lee and Chao, 2013), the emulation of behaviours (Aggarwal and McGill, 2012) or greater attribution of brand personality and brand liking (Delbaere, McQuarrie and Phillips, 2011). Importantly, research now shows that levels of this tendency vary between individuals (Waytz, Cacioppo and Epley, 2010), but research to date has failed to focus on how anthropomorphic tendency influences individual responses to marketing communications messages. Spokes-characters present an ideal context through which to examine this gap, given that they function as personified brands, designed to trigger consumer anthropomorphic tendency. Further, little is understood about how spokes-characters operate and which consumers will prefer them to their human counterparts. Like anthropomorphic research, much empirical work to date has focused on design and outcomes, examining the sender’s encoding process and the feedback generated, but ignoring the individual decoding process that is so important to understanding individual differences and message effectiveness. The current research employs three experiments using an online survey with stimulus exposure to show that anthropomorphic tendency, personality similarity and spokes-character type all have relevance to the understanding of this complex relationship. Study one and two indicate that while a human spokesperson is still preferred by many, higher levels of anthropomorphic tendency increase likeability of cartoon spokes characters. Study three highlights the importance of personality similarity, which further increases likability. Additional analyses provide key findings concerning the nature of anthropomorphic tendency as an individual difference and trait. This research contributes to a greater understanding of anthropomorphism theory and fills existing gaps in the consumer psychology and marketing communications literature.

IS-RT-PCR assay detection of MT-MMP in a human breast cancer cell line

The in situ-reverse transcription-polymerase chain reaction (IS-RT-PCR) is a method that allows the direct localisation of gene expression. The method utilises the dual buffer mediated activity of the enzyme rTth DNA polymerase enabling both reverse transcription and DNA amplification. Labelled nucleoside triphosphates allow the site of expression to be labelled, rather than the PCR primers themselves, giving a more accurate localisation of transcript expression and decreased background than standard in situ hybridisation (ISH) assays. The MDA-MB-231 human breast carcinoma (HBC) cell line was assayed via the IS-RT-PCR technique, using primers encoding MT-MMP (membrane-type matrix metalloproteinase) and human β-actin. Our results clearly indicate baseline expression of MT-MMP in the relatively invasive MDA-MB-231 cell line at a signal intensity similar to the housekeeping gene β-actin, and results following induction with Concanavalin A (Con A) are consistent with our previous results obtained via Northern blotting.

Epstein-Barr virus-positive diffuse large B-cell lymphoma of the elderly expresses EBNA3A with conserved CD8+ T-cell epitopes

Post-transplantation lymphoproliferative disorders (PTLD) arise in the immunosuppressed and are frequently Epstein-Barr virus (EBV) associated. The most common PTLD histological sub-type is diffuse large B-cell lymphoma (EBV+DLBCL-PTLD). Restoration of EBV-specific T-cell immunity can induce EBV+DLBCL-PTLD regression. The most frequent B-cell lymphoma in the immunocompetent is also DLBCL. ‘EBV-positive DLBCL of the elderly’ (EBV+DLBCL) is a rare but well-recognized DLBCL entity that occurs in the overtly immunocompetent, that has an adverse outcome relative to EBV-negative DLBCL. Unlike PTLD (which is classified as viral latency III), literature suggests EBV+DLBCL is typically latency II, i.e. expression is limited to the immuno-subdominant EBNA1, LMP1 and LMP2 EBV-proteins. If correct, this would be a major impediment for T-cell immunotherapeutic strategies. Unexpectedly we observed EBV+DLBCL-PTLD and EBV+DLBCL both shared features consistent with type III EBV-latency, including expression of the immuno-dominant EBNA3A protein. Extensive analysis showed frequent polymorphisms in EBNA1 and LMP1 functionally defined CD8+ T-cell epitope encoding regions, whereas EBNA3A polymorphisms were very rare making this an attractive immunotherapy target. As with EBV+DLBCL-PTLD, the antigen presenting machinery within lymphomatous nodes was intact. EBV+DLBCL express EBNA3A suggesting it is amenable to immunotherapeutic strategies.

Exome sequencing identifies recurrent somatic RAC1 mutations in melanoma

We characterized the mutational landscape of melanoma, the form of skin cancer with the highest mortality rate, by sequencing the exomes of 147 melanomas. Sun-exposed melanomas had markedly more ultraviolet (UV)-like C>T somatic mutations compared to sun-shielded acral, mucosal and uveal melanomas. Among the newly identified cancer genes was PPP6C, encoding a serine/threonine phosphatase, which harbored mutations that clustered in the active site in 12% of sun-exposed melanomas, exclusively in tumors with mutations in BRAF or NRAS. Notably, we identified a recurrent UV-signature, an activating mutation in RAC1 in 9.2% of sun-exposed melanomas. This activating mutation, the third most frequent in our cohort of sun-exposed melanoma after those of BRAF and NRAS, changes Pro29 to serine (RAC1P29S) in the highly conserved switch I domain. Crystal structures, and biochemical and functional studies of RAC1P29S showed that the alteration releases the conformational restraint conferred by the conserved proline, causes an increased binding of the protein to downstream effectors, and promotes melanocyte proliferation and migration. These findings raise the possibility that pharmacological inhibition of downstream effectors of RAC1 signaling could be of therapeutic benefit.

Mice selectively bred for High and Low fear behavior show differences in the number of pMAPK (p44/42 ERK) expressing neurons in lateral amygdala following Pavlovian fear conditioning

Individual variability in the acquisition, consolidation and extinction of conditioned fear potentially contributes to the development of fear pathology including posttraumatic stress disorder (PTSD). Pavlovian fear conditioning is a key tool for the study of fundamental aspects of fear learning. Here, we used a selected mouse line of High and Low Pavlovian conditioned fear created from an advanced intercrossed line (AIL) in order to begin to identify the cellular basis of phenotypic divergence in Pavlovian fear conditioning. We investigated whether phosphorylated MAPK (p44/42 ERK/MAPK), a protein kinase required in the amygdala for the acquisition and consolidation of Pavlovian fear memory, is differentially expressed following Pavlovian fear learning in the High and Low fear lines. We found that following Pavlovian auditory fear conditioning, High and Low line mice differ in the number of pMAPK-expressing neurons in the dorsal sub nucleus of the lateral amygdala (LAd). In contrast, this difference was not detected in the ventral medial (LAvm) or ventral lateral (LAvl) amygdala sub nuclei or in control animals. We propose that this apparent increase in plasticity at a known locus of fear memory acquisition and consolidation relates to intrinsic differences between the two fear phenotypes. These data provide important insights into the micronetwork mechanisms encoding phenotypic differences in fear. Understanding the circuit level cellular and molecular mechanisms that underlie individual variability in fear learning is critical for the development of effective treatment of fear-related illnesses such as PTSD.

The structure of Pavlovian fear conditioning in the amygdala

Do different brains forming a specific memory allocate the same groups of neurons to encode it? One way to test this question is to map neurons encoding the same memory and quantitatively compare their locations across individual brains. In a previous study, we used this strategy to uncover a common topography of neurons in the dorsolateral amygdala (LAd) that expressed a learning-induced and plasticity-related kinase (p42/44 mitogen-activated protein kinase; pMAPK), following auditory Pavlovian fear conditioning. In this series of experiments, we extend our initial findings to ask to what extent this functional topography depends upon intrinsic neuronal structure. We first showed that the majority (87 %) of pMAPK expression in the lateral amygdala was restricted to principal-type neurons. Next, we verified a neuroanatomical reference point for amygdala alignment using in vivo magnetic resonance imaging and in vitro morphometrics. We then determined that the topography of neurons encoding auditory fear conditioning was not exclusively governed by principal neuron cytoarchitecture. These data suggest that functional patterning of neurons undergoing plasticity in the amygdala following Pavlovian fear conditioning is specific to memory formation itself. Further, the spatial allocation of activated neurons in the LAd was specific to cued (auditory), but not contextual, fear conditioning. Spatial analyses conducted at another coronal plane revealed another spatial map unique to fear conditioning, providing additional evidence that the functional topography of fear memory storing cells in the LAd is non-random and stable. Overall, these data provide evidence for a spatial organizing principle governing the functional allocation of fear memory in the amygdala.

Pavlovian Fear Conditioning Activates a Common Pattern of Neurons in the Lateral Amygdala of Individual Brains

Understanding the physical encoding of a memory (the engram) is a fundamental question in neuroscience. Although it has been established that the lateral amygdala is a key site for encoding associative fear memory, it is currently unclear whether the spatial distribution of neurons encoding a given memory is random or stable. Here we used spatial principal components analysis to quantify the topography of activated neurons, in a select region of the lateral amygdala, from rat brains encoding a Pavlovian conditioned fear memory. Our results demonstrate a stable, spatially patterned organization of amygdala neurons are activated during the formation of a Pavlovian conditioned fear memory. We suggest that this stable neuronal assembly constitutes a spatial dimension of the engram. © 2011 This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.