Researcher practices on returning genetic research results.

BACKGROUND/AIMS: as genetic and genomic research proliferates, debate has ensued about returning results to participants. In addition to consideration of the benefits and harms to participants, researchers must also consider the logistical and financial feasibility of returning research results. However, little data exist of actual researcher practices. METHODS: we conducted an online survey of 446 corresponding authors of genetic/genomic studies conducted in the United States and published in 2006-2007 to assess the frequency with which they considered, offered to, or actually returned research results, what factors influenced these decisions, and the method of communicating results. RESULTS: the response rate was 24% (105/446). Fifty-four percent of respondents considered the issue of returning research results to participants, 28% offered to return individual research results, and 24% actually returned individual research results. Of those who considered the issue of returning research results during the study planning phase, the most common factors considered were whether research results were deemed clinically useful (18%) and respect for participants (13%). Researchers who had a medical degree and conducted studies on children were significantly more likely to offer to return or actually return individual results compared to those with a Ph.D. only. CONCLUSIONS: we speculate that issues associated with clinical validity and respect for participants dominated concerns of time and expense given the prominent and continuing ethical debates surrounding genetics and genomics research. The substantial number of researchers who did not consider returning research results suggests that researchers and institutional review boards need to devote more attention to a topic about which research participants are interested.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

Diagnosis of partial body radiation exposure in mice using peripheral blood gene expression profiles.

In the event of a terrorist-mediated attack in the United States using radiological or improvised nuclear weapons, it is expected that hundreds of thousands of people could be exposed to life-threatening levels of ionizing radiation. We have recently shown that genome-wide expression analysis of the peripheral blood (PB) can generate gene expression profiles that can predict radiation exposure and distinguish the dose level of exposure following total body irradiation (TBI). However, in the event a radiation-mass casualty scenario, many victims will have heterogeneous exposure due to partial shielding and it is unknown whether PB gene expression profiles would be useful in predicting the status of partially irradiated individuals. Here, we identified gene expression profiles in the PB that were characteristic of anterior hemibody-, posterior hemibody- and single limb-irradiation at 0.5 Gy, 2 Gy and 10 Gy in C57Bl6 mice. These PB signatures predicted the radiation status of partially irradiated mice with a high level of accuracy (range 79-100%) compared to non-irradiated mice. Interestingly, PB signatures of partial body irradiation were poorly predictive of radiation status by site of injury (range 16-43%), suggesting that the PB molecular response to partial body irradiation was anatomic site specific. Importantly, PB gene signatures generated from TBI-treated mice failed completely to predict the radiation status of partially irradiated animals or non-irradiated controls. These data demonstrate that partial body irradiation, even to a single limb, generates a characteristic PB signature of radiation injury and thus may necessitate the use of multiple signatures, both partial body and total body, to accurately assess the status of an individual exposed to radiation.

Mechanical defenses in leaves eaten by Costa Rican howling monkeys (Alouatta palliata).

Primate species often eat foods of different physical properties. This may have implications for tooth structure and wear in those species. The purpose of this study was to examine the mechanical defenses of leaves eaten by Alouatta palliata from different social groups at Hacienda La Pacifica in Costa Rica. Leaves were sampled from the home-ranges of groups living in different microhabitats. Specimens were collected during the wet and dry seasons from the same tree, same plant part, and same degree of development as those eaten by the monkeys. The toughness of over 300 leaves was estimated using a scissors test on a Darvell mechanical tester. Toughness values were compared between social groups, seasons, and locations on the leaves using ANOVA. Representative samples of leaves were also sun-dried for subsequent scanning electron microscopy and energy dispersive x-ray (EDX) analyses in an attempt to locate silica on the leaves. Both forms of mechanical defense (toughness and silica) were found to be at work in the plants at La Pacifica. Fracture toughness varied significantly by location within single leaves, indicating that measures of fracture toughness must be standardized by location on food items. Monkeys made some food choices based on fracture toughness by avoiding the toughest parts of leaves and consuming the least tough portions. Intergroup and seasonal differences in the toughness of foods suggest that subtle differences in resource availability can have a significant impact on diet and feeding in Alouatta palliata. Intergroup differences in the incidence of silica on leaves raise the possibility of matching differences in the rates and patterns of tooth wear.

The Multiscale Systems Immunology project: software for cell-based immunological simulation.

BACKGROUND: Computer simulations are of increasing importance in modeling biological phenomena. Their purpose is to predict behavior and guide future experiments. The aim of this project is to model the early immune response to vaccination by an agent based immune response simulation that incorporates realistic biophysics and intracellular dynamics, and which is sufficiently flexible to accurately model the multi-scale nature and complexity of the immune system, while maintaining the high performance critical to scientific computing. RESULTS: The Multiscale Systems Immunology (MSI) simulation framework is an object-oriented, modular simulation framework written in C++ and Python. The software implements a modular design that allows for flexible configuration of components and initialization of parameters, thus allowing simulations to be run that model processes occurring over different temporal and spatial scales. CONCLUSION: MSI addresses the need for a flexible and high-performing agent based model of the immune system.

A network of substrates of the E3 ubiquitin ligases MDM2 and HUWE1 control apoptosis independently of p53.

In the intrinsic pathway of apoptosis, cell-damaging signals promote the release of cytochrome c from mitochondria, triggering activation of the Apaf-1 and caspase-9 apoptosome. The ubiquitin E3 ligase MDM2 decreases the stability of the proapoptotic factor p53. We show that it also coordinated apoptotic events in a p53-independent manner by ubiquitylating the apoptosome activator CAS and the ubiquitin E3 ligase HUWE1. HUWE1 ubiquitylates the antiapoptotic factor Mcl-1, and we found that HUWE1 also ubiquitylated PP5 (protein phosphatase 5), which indirectly inhibited apoptosome activation. Breast cancers that are positive for the tyrosine receptor kinase HER2 (human epidermal growth factor receptor 2) tend to be highly aggressive. In HER2-positive breast cancer cells treated with the HER2 tyrosine kinase inhibitor lapatinib, MDM2 was degraded and HUWE1 was stabilized. In contrast, in breast cancer cells that acquired resistance to lapatinib, the abundance of MDM2 was not decreased and HUWE1 was degraded, which inhibited apoptosis, regardless of p53 status. MDM2 inhibition overcame lapatinib resistance in cells with either wild-type or mutant p53 and in xenograft models. These findings demonstrate broader, p53-independent roles for MDM2 and HUWE1 in apoptosis and specifically suggest the potential for therapy directed against MDM2 to overcome lapatinib resistance.

Bayesian Gaussian Copula Factor Models for Mixed Data.

Gaussian factor models have proven widely useful for parsimoniously characterizing dependence in multivariate data. There is a rich literature on their extension to mixed categorical and continuous variables, using latent Gaussian variables or through generalized latent trait models acommodating measurements in the exponential family. However, when generalizing to non-Gaussian measured variables the latent variables typically influence both the dependence structure and the form of the marginal distributions, complicating interpretation and introducing artifacts. To address this problem we propose a novel class of Bayesian Gaussian copula factor models which decouple the latent factors from the marginal distributions. A semiparametric specification for the marginals based on the extended rank likelihood yields straightforward implementation and substantial computational gains. We provide new theoretical and empirical justifications for using this likelihood in Bayesian inference. We propose new default priors for the factor loadings and develop efficient parameter-expanded Gibbs sampling for posterior computation. The methods are evaluated through simulations and applied to a dataset in political science. The models in this paper are implemented in the R package bfa.

A host transcriptional signature for presymptomatic detection of infection in humans exposed to influenza H1N1 or H3N2.

There is great potential for host-based gene expression analysis to impact the early diagnosis of infectious diseases. In particular, the influenza pandemic of 2009 highlighted the challenges and limitations of traditional pathogen-based testing for suspected upper respiratory viral infection. We inoculated human volunteers with either influenza A (A/Brisbane/59/2007 (H1N1) or A/Wisconsin/67/2005 (H3N2)), and assayed the peripheral blood transcriptome every 8 hours for 7 days. Of 41 inoculated volunteers, 18 (44%) developed symptomatic infection. Using unbiased sparse latent factor regression analysis, we generated a gene signature (or factor) for symptomatic influenza capable of detecting 94% of infected cases. This gene signature is detectable as early as 29 hours post-exposure and achieves maximal accuracy on average 43 hours (p = 0.003, H1N1) and 38 hours (p-value = 0.005, H3N2) before peak clinical symptoms. In order to test the relevance of these findings in naturally acquired disease, a composite influenza A signature built from these challenge studies was applied to Emergency Department patients where it discriminates between swine-origin influenza A/H1N1 (2009) infected and non-infected individuals with 92% accuracy. The host genomic response to Influenza infection is robust and may provide the means for detection before typical clinical symptoms are apparent.

Bayesian inference of the number of factors in gene-expression analysis: application to human virus challenge studies.

BACKGROUND: Nonparametric Bayesian techniques have been developed recently to extend the sophistication of factor models, allowing one to infer the number of appropriate factors from the observed data. We consider such techniques for sparse factor analysis, with application to gene-expression data from three virus challenge studies. Particular attention is placed on employing the Beta Process (BP), the Indian Buffet Process (IBP), and related sparseness-promoting techniques to infer a proper number of factors. The posterior density function on the model parameters is computed using Gibbs sampling and variational Bayesian (VB) analysis. RESULTS: Time-evolving gene-expression data are considered for respiratory syncytial virus (RSV), Rhino virus, and influenza, using blood samples from healthy human subjects. These data were acquired in three challenge studies, each executed after receiving institutional review board (IRB) approval from Duke University. Comparisons are made between several alternative means of per-forming nonparametric factor analysis on these data, with comparisons as well to sparse-PCA and Penalized Matrix Decomposition (PMD), closely related non-Bayesian approaches. CONCLUSIONS: Applying the Beta Process to the factor scores, or to the singular values of a pseudo-SVD construction, the proposed algorithms infer the number of factors in gene-expression data. For real data the "true" number of factors is unknown; in our simulations we consider a range of noise variances, and the proposed Bayesian models inferred the number of factors accurately relative to other methods in the literature, such as sparse-PCA and PMD. We have also identified a "pan-viral" factor of importance for each of the three viruses considered in this study. We have identified a set of genes associated with this pan-viral factor, of interest for early detection of such viruses based upon the host response, as quantified via gene-expression data.

Dusp3 and Psme3 are associated with murine susceptibility to Staphylococcus aureus infection and human sepsis.

Using A/J mice, which are susceptible to Staphylococcus aureus, we sought to identify genetic determinants of susceptibility to S. aureus, and evaluate their function with regard to S. aureus infection. One QTL region on chromosome 11 containing 422 genes was found to be significantly associated with susceptibility to S. aureus infection. Of these 422 genes, whole genome transcription profiling identified five genes (Dcaf7, Dusp3, Fam134c, Psme3, and Slc4a1) that were significantly differentially expressed in a) S. aureus -infected susceptible (A/J) vs. resistant (C57BL/6J) mice and b) humans with S. aureus blood stream infection vs. healthy subjects. Three of these genes (Dcaf7, Dusp3, and Psme3) were down-regulated in susceptible vs. resistant mice at both pre- and post-infection time points by qPCR. siRNA-mediated knockdown of Dusp3 and Psme3 induced significant increases of cytokine production in S. aureus-challenged RAW264.7 macrophages and bone marrow derived macrophages (BMDMs) through enhancing NF-κB signaling activity. Similar increases in cytokine production and NF-κB activity were also seen in BMDMs from CSS11 (C57BL/6J background with chromosome 11 from A/J), but not C57BL/6J. These findings suggest that Dusp3 and Psme3 contribute to S. aureus infection susceptibility in A/J mice and play a role in human S. aureus infection.