Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Multiple imputations applied to the DREAM3 phosphoproteomics challenge: a winning strategy.

DREAM is an initiative that allows researchers to assess how well their methods or approaches can describe and predict networks of interacting molecules [1]. Each year, recently acquired datasets are released to predictors ahead of publication. Researchers typically have about three months to predict the masked data or network of interactions, using any predictive method. Predictions are assessed prior to an annual conference where the best predictions are unveiled and discussed. Here we present the strategy we used to make a winning prediction for the DREAM3 phosphoproteomics challenge. We used Amelia II, a multiple imputation software method developed by Gary King, James Honaker and Matthew Blackwell[2] in the context of social sciences to predict the 476 out of 4624 measurements that had been masked for the challenge. To chose the best possible multiple imputation parameters to apply for the challenge, we evaluated how transforming the data and varying the imputation parameters affected the ability to predict additionally masked data. We discuss the accuracy of our findings and show that multiple imputations applied to this dataset is a powerful method to accurately estimate the missing data. We postulate that multiple imputations methods might become an integral part of experimental design as a mean to achieve cost savings in experimental design or to increase the quantity of samples that could be handled for a given cost.

Systematic review and meta-analysis of risk factors for survival after lung metastasectomy in colorectal cancer patients

Objective: Resection of lung metastases (LM) from colorectal cancer (CRC)¦is increasingly performed with a curative intent.Most series report small groups¦of patients, and it is currently not possible to identify those CRC patients who¦may benefit the most of surgical management. It is clinically relevant to assess¦risk factors for prolonged survival after this type of procedures.¦Methods: A meta analysis of 24 series published between 2000 and 2011¦which focused on surgical management of LM from CRC and included more¦than 40 patients each, with or without prior resection of in transit liver¦metastases. Random effects were calculated for five variables considered as¦potential prognostic factors.¦Results: A total of 2815 patients who underwent surgery with a curative¦intent were considered in this analysis. Four parameters were associated with¦a decreased survival: 1) a short disease-free interval between primary tumor¦resection and development of LM (HR = 1·59, 95% CI 1·27-1·98); 2) multiple¦LM (HR = 2·04, 95%CI 1·72-2·41); 3) positive hilar/mediastinal lymph nodes¦(HR = 1·65, 95% CI 1·35-2·02); and 4) a high prethoracotomy CEA value (HR¦=1·91, 95% CI 1·57-2·32). By comparison, a history of resected liver metastases¦(HR = 1·36, 95% CI 0·92-2·03) did not achieve statistical significance.¦Conclusion: Risk factors for poor clinical outcome after surgery for lung¦metastases in CRC patients include: 1) synchronous lung metastases; 2) high¦pre-thoracotomy CEA; 3) hilar nodes involvement; and 4) multiple pulmonary¦lesions.

Insulinoma associated with a case of multiple endocrine neoplasia type I: Functional somatostatin receptors and abnormal glucose-induced insulin secretion.

A sporadic case of multiple endocrine neoplasia type I with coexisting insulinoma and hyperparathyroidism was investigated in vivo and in vitro. The insulinoma was localized by somatostatin receptor scintigraphy and these receptors were functionally active. Octreotide administration decreased the basal insulin and glucagon secretion by 90 and 46%, respectively. Immunocytochemistry of the insulinoma tissue was positive for insulin, chromogranin A and neuropeptide Y. The insulinoma cells were also isolated and cultured in vitro. Incubation experiments revealed that a low glucose concentration (1 mmol/l) was sufficient to increase cytosolic free calcium and to produce a maximal glucose-induced insulin release. Northern blot analysis of RNA obtained from the tumor showed a high abundance of the low Km glucose transporter GLUT1 but no transcript for the high Km glucose transporter GLUT2. The abnormal distribution of glucose transporters probably relates to the abnormal glucose sensing of insulinoma cells, and explains their sustained insulin secretion at low glucose concentrations. Whether these abnormalities share a pathogenetic link with the presence of functionally active somatostatin receptors remains to be elucidated.

Multiple refugia and barriers explain the phylogeography of the Valais shrew, Sorex antinorii (Mammalia: Soricomorpha), Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)-based hindcasting approach. We also used landscape genetics (including isolation-by-resistance) to infer the determinants of current intra-specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median-joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north-eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi.

Multiple adenosine 3',5'-cyclic [corrected] monophosphate response element DNA-binding proteins generated by gene diversification and alternative exon splicing.

The protein sequence deduced from the open reading frame of a human placental cDNA encoding a cAMP-responsive enhancer (CRE)-binding protein (CREB-327) has structural features characteristic of several other transcriptional transactivator proteins including jun, fos, C/EBP, myc, and CRE-BP1. Results of Southwestern analysis of nuclear extracts from several different cell lines show that there are multiple CRE-binding proteins, which vary in size in cell lines derived from different tissues and animal species. To examine the molecular diversity of CREB-327 and related proteins at the nucleic acid level, we used labeled cDNAs from human placenta that encode two different CRE-binding proteins (CREB-327 and CRE-BP1) to probe Northern and Southern blots. Both probes hybridized to multiple fragments on Southern blots of genomic DNA from various species. Alternatively, when a human placental c-jun probe was hybridized to the same blot, a single fragment was detected in most cases, consistent with the intronless nature of the human c-jun gene. The CREB-327 probe hybridized to multiple mRNAs, derived from human placenta, ranging in size from 2-9 kilobases. In contrast, the CRE-BP1 probe identified a single 4-kilobase mRNA. Sequence analyses of several overlapping human genomic cosmid clones containing CREB-327 sequences in conjunction with polymerase chain reaction indicates that the CREB-327/341 cDNAs are composed of at least eight or nine exons, and analyses of human placental cDNAs provide direct evidence for at least one alternatively spliced exon. Analyses of mouse/hamster-human hybridoma DNAs by Southern blotting and polymerase chain reaction localizes the CREB-327/341 gene to human chromosome 2. The results indicate that there is a dichotomy of CREB-like proteins, those that are related by overall structure and DNA-binding specificity as well as those that are related by close similarities of primary sequences.

Genetic variability of Aedes aegypti in the Americas using a mitochondrial gene: evidence of multiple introductions

To analyze the genetic relatedness and phylogeographic structure of Aedes aegypti, we collected samples from 36 localities throughout the Americas (Brazil, Peru, Venezuela, Guatemala, US), three from Africa (Guinea, Senegal, Uganda), and three from Asia (Singapore, Cambodia, Tahiti). Amplification and sequencing of a fragment of the mitochondrial NADH dehydrogenase subunit 4 gene identified 20 distinct haplotypes, of which 14 are exclusive to the Americas, four to African/Asian countries, one is common to the Americas and Africa, and one to the Americas and Asia. Nested clade analysis (NCA), pairwise distribution, statistical parsimony, and maximum parsimony analyses were used to infer evolutionary and historic processes, and to estimate phylogenetic relationships among haplotypes. Two clusters were found in all the analyses. Haplotypes clustered in the two clades were separated by eight mutational steps. Phylogeographic structure detected by the NCA was consistent with distant colonization within one clade and fragmentation followed by range expansion via long distance dispersal in the other. Three percent of nucleotide divergence between these two clades is suggestive of a gene pool division that may support the hypothesis of occurrence of two subspecies of Ae. aegypti in the Americas.

An in vitro model for dengue virus infection that exhibits human monocyte infection, multiple cytokine production and dexamethasone immunomodulation

An important cytokine role in dengue fever pathogenesis has been described. These molecules can be associated with haemorrhagic manifestations, coagulation disorders, hypotension and shock, all symptoms implicated in vascular permeability and disease worsening conditions. Several immunological diseases have been treated by cytokine modulation and dexamethasone is utilized clinically to treat pathologies with inflammatory and autoimmune ethiologies. We established an in vitro model with human monocytes infected by dengue virus-2 for evaluating immunomodulatory and antiviral activities of potential pharmaceutical products. Flow cytometry analysis demonstrated significant dengue antigen detection in target cells two days after infection. TNF-alpha, IFN-alpha, IL-6 and IL-10 are produced by in vitro infected monocytes and are significantly detected in cell culture supernatants by multiplex microbead immunoassay. Dexamethasone action was tested for the first time for its modulation in dengue infection, presenting optimistic results in both decreasing cell infection rates and inhibiting TNF-alpha, IFN-alpha and IL-10 production. This model is proposed for novel drug trials yet to be applyed for dengue fever.

Analysis of hepatitis C virus resistance to silibinin in vitro and in vivo points to a novel mechanism involving nonstructural protein 4B.

Intravenous silibinin (SIL) is an approved therapeutic that has recently been applied to patients with chronic hepatitis C, successfully clearing hepatitis C virus (HCV) infection in some patients even in monotherapy. Previous studies suggested multiple antiviral mechanisms of SIL; however, the dominant mode of action has not been determined. We first analyzed the impact of SIL on replication of subgenomic replicons from different HCV genotypes in vitro and found a strong inhibition of RNA replication for genotype 1a and genotype 1b. In contrast, RNA replication and infection of genotype 2a were minimally affected by SIL. To identify the viral target of SIL we analyzed resistance to SIL in vitro and in vivo. Selection for drug resistance in cell culture identified a mutation in HCV nonstructural protein (NS) 4B conferring partial resistance to SIL. This was corroborated by sequence analyses of HCV from a liver transplant recipient experiencing viral breakthrough under SIL monotherapy. Again, we identified distinct mutations affecting highly conserved amino acid residues within NS4B, which mediated phenotypic SIL resistance also in vitro. Analyses of chimeric viral genomes suggest that SIL might target an interaction between NS4B and NS3/4A. Ultrastructural studies revealed changes in the morphology of viral membrane alterations upon SIL treatment of a susceptible genotype 1b isolate, but not of a resistant NS4B mutant or genotype 2a, indicating that SIL might interfere with the formation of HCV replication sites. CONCLUSION: Mutations conferring partial resistance to SIL treatment in vivo and in cell culture argue for a mechanism involving NS4B. This novel mode of action renders SIL an attractive candidate for combination therapies with other directly acting antiviral drugs, particularly in difficult-to-treat patient cohorts.