Bayesian latent variable models for mixed discrete outcomes

In studies of complex health conditions, mixtures of discrete outcomes (event time, count, binary, ordered categorical) are commonly collected. For example, studies of skin tumorigenesis record latency time prior to the first tumor, increases in the number of tumors at each week, and the occurrence of internal tumors at the time of death. Motivated by this application, we propose a general underlying Poisson variable framework for mixed discrete outcomes, accommodating dependency through an additive gamma frailty model for the Poisson means. The model has log-linear, complementary log-log, and proportional hazards forms for count, binary and discrete event time outcomes, respectively. Simple closed form expressions can be derived for the marginal expectations, variances, and correlations. Following a Bayesian approach to inference, conditionally-conjugate prior distributions are chosen that facilitate posterior computation via an MCMC algorithm. The methods are illustrated using data from a Tg.AC mouse bioassay study.     

Distinct mechanisms of integrin binding by Yersinia pseudotuberculosis adhesins determine the phagocytic response of host macrophages

The enteropathogenic yersiniae express two outer membrane adhesins, invasin and YadA, that contribute to pathogenesis. While invasin binds directly to beta1 integrin receptors with high affinity, YadA binds indirectly through extracellular matrix (ECM) components. In this study, Yersinia pseudotuberculosis inv and yadA mutants were used to investigate how these distinct binding mechanisms compare and potentially compete in activating signalling pathways and promoting bacterial uptake by host macrophages. The efficiency of adhesin-mediated phagocytic responses was found to be dependent on the relative expression of invasin and YadA on the bacterial surface as well as the expression of ECM proteins in the extracellular milieu. Under conditions of low concentrations of ECM, invasin was found to be the dominant adhesin, promoting high levels of phagocytosis coincident with robust and sustained activation of the protein tyrosine kinases Fak and Pyk2, phosphorylation of the adaptor molecule Cas and activation of the small GTPase Rac1. In the presence of higher concentrations of ECM, YadA became the dominant functional adhesin through its ability to engage integrin receptors via an ECM bridge. We propose a model whereby invasin promotes robust and prolonged activation of phagocytic signalling cascades by inducing a 'high-affinity' integrin conformation as well as integrin clustering. We postulate that YadA-ECM promotes phagocytosis through a more transient activation of signalling cascades that arises from integrin clustering in the context of a cross-linked fibrillar ECM network.     

An interactional network of genes involved in chitin synthesis in Saccharomyces cerevisiae

Background

The identification of disease-associated genes using single nucleotide polymorphisms (SNPs) has been increasingly reported. In particular, the Affymetrix Mapping 10 K SNP microarray platform uses one PCR primer to amplify the DNA samples and determine the genotype of more than 10,000 SNPs in the human genome. This provides the opportunity for large scale, rapid and cost-effective genotyping assays for linkage analysis. However, the analysis of such datasets is nontrivial because of the large number of markers, and visualizing the linkage scores in the context of genome maps remains less automated using the current linkage analysis software packages. For example, the haplotyping results are commonly represented in the text format.

Results

Here we report the development of a novel software tool called CompareLinkage for automated formatting of the Affymetrix Mapping 10 K genotype data into the "Linkage" format and the subsequent analysis with multi-point linkage software programs such as Merlin and Allegro. The new software has the ability to visualize the results for all these programs in dChip in the context of genome annotations and cytoband information. In addition we implemented a variant of the Lander-Green algorithm in the dChipLinkage module of dChip software (V1.3) to perform parametric linkage analysis and haplotyping of SNP array data. These functions are integrated with the existing modules of dChip to visualize SNP genotype data together with LOD score curves. We have analyzed three families with recessive and dominant diseases using the new software programs and the comparison results are presented and discussed.

Conclusions

The CompareLinkage and dChipLinkage software packages are freely available. They provide the visualization tools for high-density oligonucleotide SNP array data, as well as the automated functions for formatting SNP array data for the linkage analysis programs Merlin and Allegro and calling these programs for linkage analysis. The results can be visualized in dChip in the context of genes and cytobands. In addition, a variant of the Lander-Green algorithm is provided that allows parametric linkage analysis and haplotyping.     

Protein kinase A as a therapeutic target for memory disorders: rationale and challenges

cAMP-dependent protein kinase A (PKA) signaling has a key role in memory processes and has been identified as a potential therapeutic target for memory disorders. The activation of PKA signaling is crucial for the consolidation of long-term memories dependent on the hippocampus and/or the amygdala, By contrast, initial studies indicate that cAMP-PKA activation might impair the working memory and executive functions of the prefrontal cortex. Furthermore, PKA activation in the nucleus accumbens might increase sensitivity to addiction. These complexities must be heeded when designing medications aimed at altering PKA activity. PKA might be most practical as a therapeutic target in disorders with global alterations in cAMP-PKA activity due to genetic or environmental factors.     

A simplified high-throughput method for pyrethroid knock-down resistance (kdr) detection in Anopheles gambiae

Background

A single base pair mutation in the sodium channel confers knock-down resistance to pyrethroids in many insect species. Its occurrence in Anopheles mosquitoes may have important implications for malaria vector control especially considering the current trend for large scale pyrethroid-treated bednet programmes. Screening Anopheles gambiae populations for the kdr mutation has become one of the mainstays of programmes that monitor the development of insecticide resistance. The screening is commonly performed using a multiplex Polymerase Chain Reaction (PCR) which, since it is reliant on a single nucleotide polymorphism, can be unreliable. Here we present a reliable and potentially high throughput method for screening An. gambiae for the kdr mutation.

Methods

A Hot Ligation Oligonucleotide Assay (HOLA) was developed to detect both the East and West African kdr alleles in the homozygous and heterozygous states, and was optimized for use in low-tech developing world laboratories. Results from the HOLA were compared to results from the multiplex PCR for field and laboratory mosquito specimens to provide verification of the robustness and sensitivity of the technique.

Results and Discussion

The HOLA assay, developed for detection of the kdr mutation, gives a bright blue colouration for a positive result whilst negative reactions remain colourless. The results are apparent within a few minutes of adding the final substrate and can be scored by eye. Heterozygotes are scored when a sample gives a positive reaction to the susceptible probe and the kdr probe. The technique uses only basic laboratory equipment and skills and can be carried out by anyone familiar with the Enzyme-linked immunosorbent assay (ELISA) technique. A comparison to the multiplex PCR method showed that the HOLA assay was more reliable, and scoring of the plates was less ambiguous.

Conclusion

The method is capable of detecting both the East and West African kdr alleles in the homozygous and heterozygous states from fresh or dried material using several DNA extraction methods. It is more reliable than the traditional PCR method and may be more sensitive for the detection of heterozygotes. It is inexpensive, simple and relatively safe making it suitable for use in resource-poor countries.     

Developing liquid chromatography ion mobility mass spectometry techniques

When a packet of ions in a buffer gas is exposed to a weak electric field, the ions will separate according to differences in their mobilities through the gas. This separation forms the basis of the analytical method known as ion mobility spectroscopy and is highly efficient, in that it can be carried out in a very short time frame (micro- to milliseconds). Recently, efforts have been made to couple the approach with liquid-phase separations and mass spectrometry in order to create a high-throughput and high-coverage approach for analyzing complex mixtures. This article reviews recent work to develop this approach for proteomics analyses. The instrumentation is described briefly. Several multidimensional data sets obtained upon analyzing complex mixtures are shown in order to illustrate the approach as well as provide a view of the limitations and required future work.     

Global whole-cell FTICR mass spectrometric proteomics analysis of the heat shock response in the radioresistant bacterium Deinococcus radiodurans

The results of previous studies indicated that D. radiodurans mounts a regulated protective response to heat shock, and that expression of more than 130 genes, including classical chaperones such as the groESL and dnaKJ operons and proteases such as clpB are induced in response to elevated temperature. In addition, previous qualitative whole-cell mass spectrometric studies conducted under heat shock conditions indicated global changes in the D. radiodurans proteome. To enable the discovery of novel heat shock inducible proteins as well as gain greater biological insight into the classical heat shock response at the protein level, we undertook the global whole-cell FTICR mass spectrometric proteomics study reported here. We have greatly increased the power of this approach by conducting a large number of replicate experiments in addition to taking a semiquantitative approach to data analysis, finding good reproducibility between replicates. Through this analysis, we have identified with high confidence a core set of classical heat shock proteins whose expression increases dramatically and reproducibly in response to elevated temperature. In addition, we have found that the heat shock proteome includes a large number of induced proteins that have not been identified previously as heat responsive, and have therefore been designated as candidate responders. Finally, our results are consistent with the hypothesis that elevated temperature stress could lead to cross-protection against other related stresses.     

Rat epidermal keratinocytes as an organotypic model for examining the role of Cx43 and Cx26 in skin differentiation

In order to characterize connexin expression and regulation in the epidermis, we have characterized a rat epidermal keratinocyte (REK) cell line that is phenotypically similar to basal keratinocytes in that they have the ability to differentiate into organotypic epidermis consisting of a basal cell layer, 2-3 suprabasal cell layers, and a cornified layer. RT-PCR revealed that REK cells express mRNA for Cx26, Cx31, Cx31.1, Cx37, and Cx43, which mimics the reported connexin profile for rat tissue. In addition, we report the expression of Cx30, Cx30.3, Cx40, and Cx45 in rat keratinocytes, highlighting the complexity of the connexin complement in rat epidermis. Furthermore, 3-dimensional analysis of organotypic skin revealed that Cx26 and Cx43 are exquisitely regulated during the differentiation process. The life-cycle of these connexins including their expression, transport, assembly into gap junctions, internalization, and degradation are elegantly depicted in organotypic epidermis as keratinocytes proceed from differentiation to programmed cell death.     

Changes in the metabolome of Saccharomyces cerevisiae associated with evolution in aerobic glucose-limited chemostats

The effect of culture age on intra- and extracellular metabolite levels as well as on in vitro determined specific activities of enzymes of central carbon metabolism was investigated during evolution for over 90 generations of Saccharomyces cerevisiae CEN.PK 113-7D in an aerobic glucose/ethanol-limited chemostat at a specific dilution rate of 0.052 h(-1). It was found that the fluxes of consumed (O2, glucose/ethanol) and secreted compounds (CO2) did not change significantly during the entire cultivation period. However, morphological changes were observed, leading to an increased cellular surface area. During 90 generations of chemostat growth not only the residual glucose concentration decreased, also the intracellular concentrations of trehalose, glycolytic intermediates, TCA cycle intermediates and amino acids were found to have decreased with a factor 5-10. The only exception was glyoxylate which showed a fivefold increase in concentration. In addition to this the specific activities of most glycolytic enzymes also decreased by a factor 5-10 during long-term cultivation. Exceptions to this were hexokinase, phosphofructokinase, pyruvate kinase and 6-phosphogluconate dehydrogenase of which the activities remained unchanged. Furthermore, the concentrations of the adenylate nucleotides as well as the energy charge of the cells did not change in a significant manner. Surprisingly, the specific activities of glucose-6-phosphate dehydrogenase (G6PDH), malate synthase (MS) and isocitrate lyase (ICL) increased significantly during 90 generations of chemostat cultivation. These changes seem to indicate a pattern where metabolic overcapacities (for reversible reactions) and storage pools (trehalose, high levels of amino acids and excess protein in enzymes) are lost during the evolution period. The driving force is proposed to be a growth advantage in the absence of these metabolic overcapacities.