A modular approach for integrative analysis of large-scale gene-expression and drug-response data

High-throughput technologies are now used to generate more than one type of data from the same biological samples. To properly integrate such data, we propose using co-modules, which describe coherent patterns across paired data sets, and conceive several modular methods for their identification. We first test these methods using in silico data, demonstrating that the integrative scheme of our Ping-Pong Algorithm uncovers drug-gene associations more accurately when considering noisy or complex data. Second, we provide an extensive comparative study using the gene-expression and drug-response data from the NCI-60 cell lines. Using information from the DrugBank and the Connectivity Map databases we show that the Ping-Pong Algorithm predicts drug-gene associations significantly better than other methods. Co-modules provide insights into possible mechanisms of action for a wide range of drugs and suggest new targets for therapy

Coevolutionary networks: a novel approach to understanding the relationships of humans with the infectious agents

Human organism is interpenetrated by the world of microorganisms, from the conception until the death. This interpenetration involves different levels of interactions between the partners including trophic exchanges, bi-directional cell signaling and gene activation, besides genetic and epigenetic phenomena, and tends towards mutual adaptation and coevolution. Since these processes are critical for the survival of individuals and species, they rely on the existence of a complex organization of adaptive systems aiming at two apparently conflicting purposes: the maintenance of the internal coherence of each partner, and a mutually advantageous coexistence and progressive adaptation between them. Humans possess three adaptive systems: the nervous, the endocrine and the immune system, each internally organized into subsystems functionally connected by intraconnections, to maintain the internal coherence of the system. The three adaptive systems aim at the maintenance of the internal coherence of the organism and are functionally linked by interconnections, in such way that what happens to one is immediately sensed by the others. The different communities of infectious agents that live within the organism are also organized into functional networks. The members of each community are linked by intraconnections, represented by the mutual trophic, metabolic and other influences, while the different infectious communities affect each other through interconnections. Furthermore, by means of its adaptive systems, the organism influences and is influenced by the microbial communities through the existence of transconnections. It is proposed that these highly complex and dynamic networks, involving gene exchange and epigenetic phenomena, represent major coevolutionary forces for humans and microorganisms.

PAP1 [poly(A) polymerase 1] homozygosity and hyperadenylation are major determinants of increased mRNA stability of CDR1 in azole-resistant clinical isolates of Candida albicans.

Using genetically matched azole-susceptible (AS) and azole-resistant (AR) clinical isolates of Candida albicans, we recently demonstrated that CDR1 overexpression in AR isolates is due to its enhanced transcriptional activation and mRNA stability. This study examines the molecular mechanisms underlying enhanced CDR1 mRNA stability in AR isolates. Mapping of the 3' untranslated region (3' UTR) of CDR1 revealed that it was rich in adenylate/uridylate (AU) elements, possessed heterogeneous polyadenylation sites, and had putative consensus sequences for RNA-binding proteins. Swapping of heterologous and chimeric lacZ-CDR1 3' UTR transcriptional reporter fusion constructs did not alter the reporter activity in AS and AR isolates, indicating that cis-acting sequences within the CDR1 3' UTR itself are not sufficient to confer the observed differential mRNA decay. Interestingly, the poly(A) tail of the CDR1 mRNA of AR isolates was approximately 35-50 % hyperadenylated as compared with AS isolates. C. albicans poly(A) polymerase (PAP1), responsible for mRNA adenylation, resides on chromosome 5 in close proximity to the mating type-like (MTL) locus. Two different PAP1 alleles, PAP1-a/PAP1-alpha, were recovered from AS (MTL-a/MTL-alpha), while a single type of PAP1 allele (PAP1-alpha) was recovered from AR isolates (MTL-alpha/MTL-alpha). Among the heterozygous deletions of PAP1-a (Deltapap1-a/PAP1-alpha) and PAP1-alpha (PAP1-a/Deltapap1-alpha), only the former led to relatively enhanced drug resistance, to polyadenylation and to transcript stability of CDR1 in the AS isolate. This suggests a dominant negative role of PAP1-a in CDR1 transcript polyadenylation and stability. Taken together, our study provides the first evidence, to our knowledge, that loss of heterozygosity at the PAP1 locus is linked to hyperadenylation and subsequent increased stability of CDR1 transcripts, thus contributing to enhanced drug resistance.

Leishmania cysteine proteinases: from gene to subunit vaccine.

Whole genome sequences of microbial pathogens present new opportunities for clinical application. Presently, genome sequencing of the human protozoan parasite Leishmania major is in progress. The driving forces behind the genome project are to identify genes with key cellular functions and new drug targets, to increase knowledge on mechanisms of drug resistance and to favor technology transfer to scientists from endemic countries. Sequencing of the genome is also aimed at the identification of genes that are expressed in the infectious stages of the parasite and in particular in the intracellular form of the parasite. Several protective antigens of Leishmania have been identified. In addition to these antigens, lysosomal cysteine proteinases (CPs) have been characterized in different strains of Leishmania and Trypanosoma, as new target molecules. Recently, we have isolated and characterized Type I (CPB) and Type II (CPA) cysteine proteinase encoding genes from L. major. The exact function of cysteine proteinases of Leishmania is not completely understood, although there are a few reports describing their role as virulence factors. One specific feature of CPB in Leishmania and other trypanosomatids, is the presence of a Cterminal extension (CTE) which is possibly indicative of conserved structure and function. Recently, we demonstrated that DNA immunization of genetically susceptible BALB / c mice, using a cocktail of CPB and CPA genes, induced long lasting protection against L. major infection. This review intends to give an overview of the current knowledge on genetic vaccination used against leishmaniasis and the importance of CP genes for such an approach.

Molecular Modeling Approaches for Determining Gene Function: application to a Putative Poly-A Binding Protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 aminoacids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.

Different epitopes of the LACK protein are recognized by V beta 4 V alpha 8 CD4+ T cells in H-2b and H-2d mice susceptible to Leishmania major.

After inoculation of Leishmania major, a rapid production of IL-4 by LACK-specific CD4+ T cells has been shown to drive Th2 cell development in susceptible mice i.e. BALB/c and C57BL/6 mice rendered susceptible by neutralization of IFN-gamma at the onset of infection. Here, we showed that peptide AA 156-173 induced an early IL-4 mRNA expression not only in BALB/c mice but also in resistant B10.D2 mice when IFN-gamma is neutralized. Epitope mapping of LACK protein demonstrated that peptide containing AA 293-305 induced early IL-4 mRNA transcripts in susceptible H-2b mice i.e. BALB/b and resistant C57BL/6 mice when IFN-gamma is neutralized. Stringently, the early IL-4 response to the H-2d (AA 156-173) or the H-2b (AA 293-305) epitopes occurred in V beta 4 V alpha 8 CD4+ T cells from either H-2d or H-2b susceptible mice, respectively.

Gene expression databases and data mining.

The DNA microarray technology has arguably caught the attention of the worldwide life science community and is now systematically supporting major discoveries in many fields of study. The majority of the initial technical challenges of conducting experiments are being resolved, only to be replaced with new informatics hurdles, including statistical analysis, data visualization, interpretation, and storage. Two systems of databases, one containing expression data and one containing annotation data are quickly becoming essential knowledge repositories of the research community. This present paper surveys several databases, which are considered "pillars" of research and important nodes in the network. This paper focuses on a generalized workflow scheme typical for microarray experiments using two examples related to cancer research. The workflow is used to reference appropriate databases and tools for each step in the process of array experimentation. Additionally, benefits and drawbacks of current array databases are addressed, and suggestions are made for their improvement.

Comparative study and identification of potent eukaryotic transcriptional repressors in gene switch systems.

In mammalian cells, proper gene regulation is achieved by the complex interplay of transcription factors that activate or repress gene expression by binding to the regulatory regions of target promoters. While transcriptional activators have been extensively characterised and classified into functional groups, relatively little is known about the comparative strength and cell type-specificity of transcriptional repressors. Here, we have compared the ability of a series of eukaryotic repression domains to silence basal and activated transcription. A series of the most potent repression domains was further tested in the context of a gene therapy gene-switch system in various cell types. The results indicate that the analysed repression domains exert varying silencing activities in different promoter contexts. Furthermore, their potential for gene silencing varies also depending on the cellular context. When multimerised within one chimeric repressor protein, particular combinations of repressor domains were found to display synergistic repressing effects and efficient repression in a panel of cell lines. This approach thus allowed the identification of transcriptional repressors that are both potent and versatile in terms of cellular specificity as a basis for gene switch systems.

An experimental study on the shallow-level migmatization of ferrogabbros from the Fuerteventura Basal Complex, Canary Islands

Spectacular shallow-level migmatization of ferrogabbroic rocks occurs in a metamorphic contact aureole of a gabbroic pluton of the Tierra Mala massif (TM) on Fuerteventura (Canary Islands). In order to improve our knowledge of the low pressure melting behavior of gabbroic rocks and to constrain the conditions of migmatization of the TM gabbros, we performed partial melting experiments on a natural ferrogabbro, which is assumed as protolith of the migmatites. The experiments were performed in an internally heated pressure vessel (IHPV) at 200 MPa, 930-1150 degreesC at relatively oxidizing conditions. Distinct amounts of water were added to the charge. From 930 to 1000 degreesC, the observed experimental phases are plagioclase (An(60-70)), clinopyroxene, amphibole (titanian magnesiohastingsites), two Fe-Ti oxides, and a basaltic, K-poor melt. Above 1000 degreesC, amphibole is no longer stable. The first melts are very rich in non-native plagioclase (>70 wt.%). This indicates that at the beginning of partial melting plagioclase is the major phase which is consumed to produce melt. In the experiments, plagioclase is stable up to high temperatures (1060 degreesC) showing increasing An content with temperature. This is not compatible with the natural migmatites, in which An-rich plagioclase is absent in the melanosomes, while amphibole is stable. Our results show that the partial melting of the natural rocks cannot be regarded as an ``in-situ'' process that occurred in a closed system. Considerable amounts of alkalis probably transported by water-rich fluids, derived from the mafic pluton underplating the TM gabbro, were necessary to drive the melting reaction out of the stability range of plagioclase. A partial melting experiment with a migmatite gabbro showing typical ``in-situ'' textures as starting material supports this assumption. Crystallization experiments performed at 1000 degreesC on a glass of the fitised ferrogabbro with different water contents added to the charge show that generally high water activities could be achieved (crystallization of amphibole), independently of the bulk water content, even in a system with very low initial bulk water content (0.3 wt.%). Increasing water contents produce plagioclase richer in An, reduces the modal proportion of plagioclase in the crystallizing assemblage and extends the melt fraction. High melt fractions of >30 wt.% could only be observed in systems with high bulk water contents (> - 2 wt.%). This indicates that the migmatites were generated under water-rich conditions (probably water-saturated), since those migmatites, which are characterized as ``in-situ'' formations, show generally high amounts of leucosomes (>30 wt.%). (C) 2003 Elsevier B.V. All rights reserved.

Late Jurassic oceanic crust and Upper Cretaceous Caribbean plateau picritic basalts exposed in the Duarte igneous complex, Hispaniola

Four distinct rock units have been recognized near El Aguacate, in the Janico-Juncalito-La Vega area of the Duarte complex (Dominican Republic): (1) serpentinites crosscut by numerous diabasic dikes, (2) basalts interbedded with Late Jurassic ribbon cherts, (3) picrites and ankaramites relatively enriched in incompatible trace elements, and (4) amphibolites and gneissic amphibolites chemically similar to Oceanic Plateau Basalts. Similar Ar-Ar ages of late magmatic amphibole from a picrite, and hornblende from an amphibolite (86.1 +/- 1.3 Ma and 86.7 +/- 1.6 Ma, respectively), suggest that the Duarte picrites are contemporaneous with the Deep Sea Drilling Program Leg 15 and Ocean Drilling Program Leg 126 basalts drilled from the Caribbean oceanic plateau. These basalts are associated with sediments containing Late Cretaceous faunas. Sr, Nd, and Pb data show that enriched picrites and amphibolites are isotopically similar to mafic lavas from previously described Caribbean plateau and Galapagos hotspot basalts. Major element, trace element, and lead isotopic features of Late Jurassic basalts and diabases are consistent with those of normal oceanic crust basalt. However, these basalts differ from typical N-MORB because they have lower epsilon Nd ratios that plot within the range of Ocean Island Basalts. These rocks appear to represent remnants of the Caribbean Jurassic oceanic crust formed from an oceanic ridge possibly close to a hotspot. Later, they were tectonically juxtaposed with Late Cretaceous slices of the Caribbean-Colombian plateau.

Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes

Mammalian sex chromosomes stem from ancestral autosomes and have substantially differentiated. It was shown that X-linked genes have generated duplicate intronless gene copies (retrogenes) on autosomes due to this differentiation. However, the precise driving forces for this out-of-X gene "movement" and its evolutionary onset are not known. Based on expression analyses of male germ-cell populations, we here substantiate and extend the hypothesis that autosomal retrogenes functionally compensate for the silencing of their X-linked housekeeping parental genes during, but also after, male meiotic sex chromosome inactivation (MSCI). Thus, sexually antagonistic forces have not played a major role for the selective fixation of X-derived gene copies in mammals. Our dating analyses reveal that although retrogenes were produced ever since the common mammalian ancestor, selectively driven retrogene export from the X only started later, on the placental mammal (eutherian) and marsupial (metatherian) lineages, respectively. Together, these observations suggest that chromosome-wide MSCI emerged close to the eutherian-marsupial split approximately 180 million years ago. Given that MSCI probably reflects the spread of the recombination barrier between the X and Y, crucial for their differentiation, our data imply that these chromosomes became more widely differentiated only late in the therian ancestor, well after the divergence of the monotreme lineage. Thus, our study also provides strong independent support for the recent notion that our sex chromosomes emerged, not in the common ancestor of all mammals, but rather in the therian ancestor, and therefore are much younger than previously thought

Arrested kinetic Li isotope fractionation at the margin of the llimaussaq complex, South Greenland: Evidence for open-system processes during final cooling of peralkaline igneous rocks

Li contents [Li] and isotopic composition (delta Li-7) of mafic minerals (mainly amphibole and clinopyroxene) from the alkaline to peralkaline Ilimaussaq plutonic complex, South Greenland, track the behavior of Li and its isotopes during magmatic differentiation and final cooling of an alkaline igneous system. [Li] in amphibole increase from < 10 ppm in Caamphiboles of the least differentiated unit to >3000 ppm in Na-amphiboles of the highly evolved units. In contrast, [Li] in clinopyroxene are comparatively low (<85 ppm) and do not vary systematically with differentiation. The distribution of Li between amphibole and pyroxene is controlled by the major element composition of the minerals (Ca-rich and Na-rich, respectively) and changes in oxygen fugacity (due to Li incorporation via coupled substitution with ferric iron) during magmatic differentiation. delta(7) Li values of all minerals span a wide range from + 17 to - 8 parts per thousand, with the different intrusive units of the complex having distinct Li isotopic systematics. Amphiboles, which dominate the Li budget of whole-rocks from the inner part of the complex, have constant delta Li-7 of + 1.8 +/- 2.2 parts per thousand (2 sigma, n = 15). This value reflects a homogeneous melt reservoir and is consistent with their mantle derivation, in agreement with published O and Nd isotopic data. Clinopyroxenes of these samples are consistently lighter, with Delta Li-7(amph-cpx). as large as 8 parts per thousand and are thus not in Li isotope equilibrium. These low values probably reflect late-stage diffusion of Li into clinopyroxene during final cooling of the rocks, thus enriching the clinopyroxene in 6 Li. At the margin of the complex delta(7) Li in the syenites increases systematically, from +2 to high values of + 14 parts per thousand. This, coupled with the observed Li isotope systematics of the granitic country rocks, reflects post-magmatic open-system processes occurring during final cooling of the intrusion. Although the shape and magnitude of the Li isotope and elemental profiles through syenite and country rock are suggestive of diffusion-driven isotope fractionation, they cannot be modeled by one-dimensional diffusive transport and point to circulation of a fluid having a high 67 Li value (possibly seawater) along the chilled contact. In all, this study demonstrates that Li isotopes can be used to identify complex fluid- and diffusion-governed processes taking place during the final cooling of such rocks. (c) 2007 Elsevier B.V All rights reserved.

Post-genomic update on a classical candidate gene for coronary artery disease: ESR1.

BACKGROUND: After age, sex is the most important risk factor for coronary artery disease (CAD). The mechanism through which women are protected from CAD is still largely unknown, but the observed sex difference suggests the involvement of the reproductive steroid hormone signaling system. Genetic association studies of the gene-encoding Estrogen Receptor α (ESR1) have shown conflicting results, although only a limited range of variation in the gene has been investigated. METHODS AND RESULTS: We exploited information made available by advanced new methods and resources in complex disease genetics to revisit the question of ESR1's role in risk of CAD. We performed a meta-analysis of 14 genome-wide association studies (CARDIoGRAM discovery analysis, N=≈87,000) to search for population-wide and sex-specific associations between CAD risk and common genetic variants throughout the coding, noncoding, and flanking regions of ESR1. In addition to samples from the MIGen (N=≈6000), WTCCC (N=≈7400), and Framingham (N=≈3700) studies, we extended this search to a larger number of common and uncommon variants by imputation into a panel of haplotypes constructed using data from the 1000 Genomes Project. Despite the widespread expression of ERα in vascular tissues, we found no evidence for involvement of common or low-frequency genetic variation throughout the ESR1 gene in modifying risk of CAD, either in the general population or as a function of sex. CONCLUSIONS: We suggest that future research on the genetic basis of sex-related differences in CAD risk should initially prioritize other genes in the reproductive steroid hormone biosynthesis system.

The fou2 gain-of-function allele and the wild-type allele of Two Pore Channel 1 contribute to different extents or by different mechanisms to defense gene expression in Arabidopsis.

The fatty acid oxygenation up-regulated 2 (fou2) mutant in Arabidopsis thaliana creates a gain-of-function allele in a non-selective cation channel encoded by the Two Pore Channel 1 (TPC1) gene. This mutant genetically implicates cation fluxes in the control of the positive feedback loop whereby jasmonic acid (JA) stimulates its own synthesis. In this study we observed extensive transcriptome reprogramming in healthy fou2 leaves closely resembling that induced by treatment with methyl jasmonate, biotic stresses and the potassium starvation response. Proteomic analysis of fou2 leaves identified increased levels of seven biotic stress- and JA-inducible proteins. In agreement with these analyses, epistasis studies performed by crossing fou2 with aos indicated that elevated levels of JA in fou2 are the major determinant of the mutant phenotype. In addition, generation of fou2 aba1-5, fou2 etr1-1 and fou2 npr1-1 double mutants showed that the fou2 phenotype was only weakly affected by ABA levels and unaffected by mutations in NPR1 and ETR1. The results now suggest possible mechanisms whereby fou2 could induce JA synthesis/signaling early in the wound response. In contrast to fou2, transcriptome analysis of a loss-of-function allele of TPC1, tpc1-2, revealed no differential expression of JA biosynthesis genes in resting leaves. However, the analysis disclosed reduced mRNA levels of the pathogenesis-related genes PDF1.2a and THI2.1 in healthy and diseased tpc1-2 leaves. The results suggest that wild-type TPC1 contributes to their expression by mechanisms somewhat different from those affecting their expression in fou2.

A mucin like gene different from the previously reported members of the mucin like gene families is transcribed in Trypanosoma cruzi but not in Trypanosoma rangeli

Trypanosoma cruzi expresses mucin like glycoproteins encoded by a complex multigene family. In this work, we report the transcription in T. cruzi but not in T. rangeli of a mucin type gene automatically annotated by the T. cruzi genome project. The gene showed no nucleotide similarities with the previously reported T. cruzi mucin like genes, although the computational analysis of the deduced protein showed that it has the characteristic features of mucins: a signal peptide sequence, O-glycosylation sites, and glycosylphosphatidylinositol (GPI) anchor sequence. The presence in this gene of N- terminal and C- terminal coding sequences common to other annotated mucin like genes suggests the existence of a new mucin like gene family.