Study of cerebral gene expression densities using Voronoi analysis

As the available public cerebral gene expression image data increasingly grows, the demand for automated methods to analyze such large amount of data also increases. An important study that can be carried out on these data is related to the spatial relationship between gene expressions. Similar spatial density distribution of expression between genes may indicate they are functionally correlated, thus the identification of these similarities is useful in suggesting directions of investigation to discover gene interactions and their correlated functions. In this paper, we describe the use of a high-throughput methodology based on Voronoi diagrams to automatically analyze and search for possible local spatial density relationships between gene expression images. We tested this method using mouse brain section images from the Allen Mouse Brain Atlas public database. This methodology provided measurements able to characterize the similarity of the density distribution between gene expressions and allowed the visualization of the results through networks and Principal Component Analysis (PCA). These visualizations are useful to analyze the similarity level between gene expression patterns, as well as to compare connection patterns between region networks. Some genes were found to have the same type of function and to be near each other in the PCA visualizations. These results suggest cerebral density correlations between gene expressions that could be further explored. (C) 2011 Elsevier B.V. All rights reserved.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

The SUCEST-FUN regulatory network database: Designing an energy grass

Modern sugarcane cultivars are complex hybrids resulting from crosses among several Saccharum species. Traditional breeding methods have been employed extensively in different countries over the past decades to develop varieties with increased sucrose yield and resistance to pests and diseases. Conventional variety improvement, however, may be limited by the narrow pool of suitable genes. Thus, molecular genetics is seen as a promising tool to assist in the process of developing improved varieties. The SUCEST-FUN Project (http://sucest-fun.org) aims to associate function with sugarcane genes using a variety of tools, in particular those that enable the study of the sugarcane transcriptome. An extensive analysis has been conducted to characterise, phenotypically, sugarcane genotypes with regard to their sucrose content, biomass and drought responses. Through the analysis of different cultivars, genes associated with sucrose content, yield, lignin and drought have been identified. Currently, tools are being developed to determine signalling and regulatory networks in grasses, and to sequence the sugarcane genome, as well as to identify sugarcane promoters. This is being implemented through the SUCEST-FUN (http://sucest-fun.org) and GRASSIUS databases (http://grassius.org), the cloning of sugarcane promoters, the identification of cis-regulatory elements (CRE) using Chromatin Immunoprecipitation-sequencing (ChIP-Seq) and the generation of a comprehensive Signal Transduction and Transcription gene catalogue (SUCAST Catalogue).

Modeling the Emergence of Circadian Rhythms in a Clock Neuron Network

Circadian rhythms in pacemaker cells persist for weeks in constant darkness, while in other types of cells the molecular oscillations that underlie circadian rhythms damp rapidly under the same conditions. Although much progress has been made in understanding the biochemical and cellular basis of circadian rhythms, the mechanisms leading to damped or self-sustained oscillations remain largely unknown. There exist many mathematical models that reproduce the circadian rhythms in the case of a single cell of the Drosophila fly. However, not much is known about the mechanisms leading to coherent circadian oscillation in clock neuron networks. In this work we have implemented a model for a network of interacting clock neurons to describe the emergence (or damping) of circadian rhythms in Drosophila fly, in the absence of zeitgebers. Our model consists of an array of pacemakers that interact through the modulation of some parameters by a network feedback. The individual pacemakers are described by a well-known biochemical model for circadian oscillation, to which we have added degradation of PER protein by light and multiplicative noise. The network feedback is the PER protein level averaged over the whole network. In particular, we have investigated the effect of modulation of the parameters associated with (i) the control of net entrance of PER into the nucleus and (ii) the non-photic degradation of PER. Our results indicate that the modulation of PER entrance into the nucleus allows the synchronization of clock neurons, leading to coherent circadian oscillations under constant dark condition. On the other hand, the modulation of non-photic degradation cannot reset the phases of individual clocks subjected to intrinsic biochemical noise.

Sequence and Copy Number Analyses of HEXB Gene in Patients Affected by Sandhoff Disease: Functional Characterization of 9 Novel Sequence Variants

Sandhoff disease (SD) is a lysosomal disorder caused by mutations in the HEXB gene. To date, 43 mutations of HEXB have been described, including 3 large deletions. Here, we have characterized 14 unrelated SD patients and developed a Multiplex Ligation-dependent Probe Amplification (MLPA) assay to investigate the presence of large HEXB deletions. Overall, we identified 16 alleles, 9 of which were novel, including 4 sequence variation leading to aminoacid changes [c.626C>T (p.T209I), c.634C>A (p.H212N), c.926G>T (p.C309F), c.1451G>A (p.G484E)] 3 intronic mutations (c.1082+5G>A, c.1242+1G>A, c.1169+5G>A), 1 nonsense mutation c.146C>A (p.S49X) and 1 small in-frame deletion c.1260_1265delAGTTGA (p.V421_E422del). Using the new MLPA assay, 2 previously described deletions were identified. In vitro expression studies showed that proteins bearing aminoacid changes p.T209I and p.G484E presented a very low or absent activity, while proteins bearing the p.H212N and p.C309F changes retained a significant residual activity. The detrimental effect of the 3 novel intronic mutations on the HEXB mRNA processing was demonstrated using a minigene assay. Unprecedentedly, minigene studies revealed the presence of a novel alternative spliced HEXB mRNA variant also present in normal cells. In conclusion, we provided new insights into the molecular basis of SD and validated an MLPA assay for detecting large HEXB deletions.

Pore size regulates cell and tissue interactions with PLGA-CaP scaffolds used for bone engineering

A common subject in bone tissue engineering is the need for porous scaffolds to support cell and tissue interactions aiming at repairing bone tissue. As poly(lactide-co-glycolide)calcium phosphate (PLGACaP) scaffolds can be manufactured with different pore sizes, the aim of this study was to evaluate the effect of pore diameter on osteoblastic cell responses and bone tissue formation. Scaffolds were prepared with 85% porosity, with pore diameters in the ranges 470590, 590850 and 8501200 mu m. Rat bone marrow stem cells differentiated into osteoblasts were cultured on the scaffolds for up to 10 days to evaluate cell growth, alkaline phosphatase (ALP) activity and the gene expression of the osteoblast markers RUNX2, OSX, COL, MSX2, ALP, OC and BSP by real-time PCR. Scaffolds were implanted in critical size rat calvarial defects for 2, 4, and 8 weeks for histomorphometric analysis. Cell growth and ALP activity were not affected by the pore size; however, there was an increase in the gene expression of osteoblastic markers with the increase in the pore sizes. At 2 weeks all scaffolds displayed a similar amount of bone and blood vessels formation. At 4 and 8 weeks much more bone formation and an increased number of blood vessels were observed in scaffolds with pores of 470590 mu m. These results show that PLGACaP is a promising biomaterial for bone engineering. However, ideally, combinations of larger (similar to 1000 mu m) and smaller (similar to 500 mu m) pores in a single scaffold would optimize cellular and tissue responses during bone healing. Copyright (C) 2011 John Wiley & Sons, Ltd.

Interactions between the NO-Citrulline Cycle and Brain-derived Neurotrophic Factor in Differentiation of Neural Stem Cells

The diffusible messenger NO plays multiple roles in neuroprotection, neurodegeneration, and brain plasticity. Argininosuccinate synthase (AS) is a ubiquitous enzyme in mammals and the key enzyme of the NO-citrulline cycle, because it provides the substrate L-arginine for subsequent NO synthesis by inducible, endothelial, and neuronal NO synthase (NOS). Here, we provide evidence for the participation of AS and of the NO-citrulline cycle in the progress of differentiation of neural stem cells (NSC) into neurons, astrocytes, and oligodendrocytes. AS expression and activity and neuronal NOS expression, as well as L-arginine and NOx production, increased along neural differentiation, whereas endothelial NOS expression was augmented in conditions of chronic NOS inhibition during differentiation, indicating that this NOS isoform is amenable to modulation by extracellular cues. AS and NOS inhibition caused a delay in the progress of neural differentiation, as suggested by the decreased percentage of terminally differentiated cells. On the other hand, BDNF reversed the delay of neural differentiation of NSC caused by inhibition of NOx production. Alikely cause is the lack of NO, which up-regulated p75 neurotrophin receptor expression, a receptor required for BDNF-induced differentiation of NSC. We conclude that the NO-citrulline cycle acts together with BDNF for maintaining the progress of neural differentiation.

Sugarwin: A Sugarcane Insect-Induced Gene with Antipathogenic Activity

In sugarcane fields, colonization of the stalk by opportunistic fungi usually occurs after the caterpillar Diatraea saccharalis attacks the sugarcane plant. Plants respond to insect attack by inducing and accumulating a large set of defense proteins. Two homologues of a barley wound-inducible protein (BARWIN), sugarcane wound-inducible proteins SUGARWIN1 and SUGARWIN2, have been identified in sugarcane by an in silico analysis. Antifungal properties have been described for a number of BARWIN homologues. We report that a SUGARWIN:green fluorescent protein fusion protein is located in the endoplasmic reticulum and in the extracellular space of sugarcane plants. The induction of sugarwin transcripts occurs in response to mechanical wounding, D. saccharalis damage, and methyl jasmonate treatment. The accumulation of transcripts is late induced and is restricted to the site of the wound. Although the transcripts of sugarwin genes were strongly increased following insect attack, the protein itself did not show any effect on insect development; rather, it altered fungal morphology, leading to the apoptosis of the germlings. These results suggest that, in the course of evolution, sugarwin-encoding genes were recruited by sugarcane due to their antipathogenic activity. We rationalize that sugarcane is able to induce sugarwin gene expression in response to D. saccharalis feeding as a concerted plant response to the anticipated invasion by the fungi that typically penetrate the plant stalk after insect damage.

Extensive cross-talk and global regulators identified from an analysis of the integrated transcriptional and signaling network in Escherichia coli

To understand the regulatory dynamics of transcription factors (TFs) and their interplay with other cellular components we have integrated transcriptional, protein-protein and the allosteric or equivalent interactions which mediate the physiological activity of TFs in Escherichia coli. To study this integrated network we computed a set of network measurements followed by principal component analysis (PCA), investigated the correlations between network structure and dynamics, and carried out a procedure for motif detection. In particular, we show that outliers identified in the integrated network based on their network properties correspond to previously characterized global transcriptional regulators. Furthermore, outliers are highly and widely expressed across conditions, thus supporting their global nature in controlling many genes in the cell. Motifs revealed that TFs not only interact physically with each other but also obtain feedback from signals delivered by signaling proteins supporting the extensive cross-talk between different types of networks. Our analysis can lead to the development of a general framework for detecting and understanding global regulatory factors in regulatory networks and reinforces the importance of integrating multiple types of interactions in underpinning the interrelationships between them.

Structure and Ionic Conductivity in the Mixed-Network Former Chalcogenide Glass System [Na2S](2/3)[(B2S3)(x)(P2S5)(1-x)](1/3)

Glasses in the system [Na2S](2/3)[(B2S3)(x)(P2S5)(1-x)](1/3) (0.0 <= x <= 1.0) were prepared by the melt quenching technique, and their properties were characterized by thermal analysis and impedance spectroscopy. Their atomic-level structures were comprehensively characterized by Raman spectroscopy and B-11, P-31, and Na-23 high resolution solid state magic-angle spinning (MAS) NMR techniques. P-31 MAS NMR peak assignments were made by the presence or absence of homonuclear indirect P-31-P-31 spin-spin interactions as detected using homonuclear J-resolved and refocused INADEQUATE techniques. The extent of B-S-P connectivity in the glassy network was quantified by P-31{B-11} and B-11{P-31} rotational echo double resonance spectroscopy. The results clearly illustrate that the network modifier alkali sulfide, Na2S, is not proportionally shared between the two network former components, B and P. Rather, the thiophosphate (P) component tends to attract a larger concentration of network modifier species than predicted by the bulk composition, and this results in the conversion of P2S74-, pyrothiophosphate, Na/P = 2:1, units into PS43-, orthothiophosphate, Na/P = 3:1, groups. Charge balance is maintained by increasing the net degree of polymerization of the thioborate (B) units through the formation of covalent bridging sulfur (BS) units, B S B. Detailed inspection of the B-11 MAS NMR spectra reveals that multiple thioborate units are formed, ranging from neutral BS3/2 groups all the way to the fully depolymerized orthothioborate (BS33-) species. On the basis of these results, a comprehensive and quantitative structural model is developed for these glasses, on the basis of which the compositional trends in the glass transition temperatures (T-g) and ionic conductivities can be rationalized. Up to x = 0.4, the dominant process can be described in a simplified way by the net reaction equation P-1 + B-1 reversible arrow P-0 + B-4, where the superscripts denote the number of BS atoms for the respective network former species. Above x = 0.4, all of the thiophosphate units are of the P-0 type and both pyro-(B-1) and orthothioborate (B-0) species make increasing contributions to the network structure with increasing x. In sharp contrast to the situation in sodium borophosphate glasses, four-coordinated thioborate species are generally less abundant and heteroatomic B-S-P linkages appear to not exist. On the basis of this structural information, compositional trends in the ionic conductivities are discussed in relation to the nature of the charge-compensating anionic species and the spatial distribution of the charge carriers.

Gene Expression Profile in Response to Doxorubicin-Rapamycin Combined Treatment of HER-2-Overexpressing Human Mammary Epithelial Cell Lines

HER-2-positive breast cancers frequently sustain elevated AKT/mTOR signaling, which has been associated with resistance to doxorubicin treatment. Here, we investigated whether rapamycin, an mTOR inhibitor, increased the sensitivity to doxorubicin therapy in two HER-2-overexpressing cell lines: C5.2, which was derived from the parental HB4a by transfection with HER-2 and SKBR3, which exhibits HER-2 amplification. The epithelial mammary cell line HB4a was also analyzed. The combined treatment using 20 nmol/L of rapamycin and 30 nmol/L of doxorubicin arrested HB4a and C5.2 cells in S to G(2)-M, whereas SKBR3 cells showed an increase in the G(0)-G(1) phase. Rapamycin increased the sensitivity to doxorubicin in HER-2-overexpressing cells by approximately 2-fold, suggesting that the combination displayed a more effective antiproliferative action. Gene expression profiling showed that these results might reflect alterations in genes involved in canonical pathways related to purine metabolism, oxidative phosphorylation, protein ubiquitination, and mitochondrial dysfunction. A set of 122 genes modulated by the combined treatment and specifically related to HER-2 overexpression was determined by finding genes commonly regulated in both C5.2 and SKBR3 that were not affected in HB4a cells. Network analysis of this particular set showed a smaller subgroup of genes in which coexpression pattern in HB4a cells was disrupted in C5.2 and SKBR3. Altogether, our data showed a subset of genes that might be more robust than individual markers in predicting the response of HER-2-overexpressing breast cancers to doxorubicin and rapamycin combination. Mol Cancer Ther; 11(2); 464-74. (C) 2011 AACR.

Molecular determinants of caste differentiation in the highly eusocial honeybee Apis mellifera

Abstract Background In honeybees, differential feeding of female larvae promotes the occurrence of two different phenotypes, a queen and a worker, from identical genotypes, through incremental alterations, which affect general growth, and character state alterations that result in the presence or absence of specific structures. Although previous studies revealed a link between incremental alterations and differential expression of physiometabolic genes, the molecular changes accompanying character state alterations remain unknown. Results By using cDNA microarray analyses of >6,000 Apis mellifera ESTs, we found 240 differentially expressed genes (DEGs) between developing queens and workers. Many genes recorded as up-regulated in prospective workers appear to be unique to A. mellifera, suggesting that the workers' developmental pathway involves the participation of novel genes. Workers up-regulate more developmental genes than queens, whereas queens up-regulate a greater proportion of physiometabolic genes, including genes coding for metabolic enzymes and genes whose products are known to regulate the rate of mass-transforming processes and the general growth of the organism (e.g., tor). Many DEGs are likely to be involved in processes favoring the development of caste-biased structures, like brain, legs and ovaries, as well as genes that code for cytoskeleton constituents. Treatment of developing worker larvae with juvenile hormone (JH) revealed 52 JH responsive genes, specifically during the critical period of caste development. Using Gibbs sampling and Expectation Maximization algorithms, we discovered eight overrepresented cis-elements from four gene groups. Graph theory and complex networks concepts were adopted to attain powerful graphical representations of the interrelation between cis-elements and genes and objectively quantify the degree of relationship between these entities. Conclusion We suggest that clusters of functionally related DEGs are co-regulated during caste development in honeybees. This network of interactions is activated by nutrition-driven stimuli in early larval stages. Our data are consistent with the hypothesis that JH is a key component of the developmental determination of queen-like characters. Finally, we propose a conceptual model of caste differentiation in A. mellifera based on gene-regulatory networks.

Analysis of Schistosoma mansoni genes shared with Deuterostomia and with possible roles in host interactions

Abstract Background: Schistosoma mansoni is a blood helminth parasite that causes schistosomiasis, a disease that affects 200 million people in the world. Many orthologs of known mammalian genes have been discovered in this parasite and evidence is accumulating that some of these genes encode proteins linked to signaling pathways in the parasite that appear to be involved with growth or development, suggesting a complex co-evolutionary process. Results: In this work we found 427 genes conserved in the Deuterostomia group that have orthologs in S. mansoni and no members in any nematodes and insects so far sequenced. Among these genes we have identified Insulin Induced Gene (INSIG), Interferon Regulatory Factor (IRF) and vasohibin orthologs, known to be involved in mammals in mevalonate metabolism, immune response and angiogenesis control, respectively. We have chosen these three genes for a more detailed characterization, which included extension of their cloned messages to obtain full-length sequences. Interestingly, SmINSIG showed a 10-fold higher expression in adult females as opposed to males, in accordance with its possible role in regulating egg production. SmIRF has a DNA binding domain, a tryptophan-rich N-terminal region and several predicted phosphorylation sites, usually important for IRF activity. Fourteen different alternatively spliced forms of the S. mansoni vasohibin (SmVASL) gene were detected that encode seven different protein isoforms including one with a complete C-terminal end, and other isoforms with shorter C-terminal portions. Using S. mansoni homologs, we have employed a parsimonious rationale to compute the total gene losses/gains in nematodes, arthropods and deuterostomes under either the Coelomata or the Ecdysozoa evolutionary hypotheses; our results show a lower losses/gains number under the latter hypothesis. Conclusion: The genes discussed which are conserved between S. mansoni and deuterostomes, probably have an ancient origin and were lost in Ecdysozoa, being still present in Lophotrochozoa. Given their known functions in Deuterostomia, it is possible that some of them have been co-opted to perform functions related (directly or indirectly) to host adaptation or interaction with host signaling processes.

DNA vaccine containing the mycobacterial hsp65 gene prevented insulitis in MLD-STZ diabetes

Abstract Background Our group previously demonstrated that a DNA plasmid encoding the mycobacterial 65-kDa heat shock protein (DNA-HSP65) displayed prophylactic and therapeutic effect in a mice model for tuberculosis. This protection was attributed to induction of a strong cellular immunity against HSP65. As specific immunity to HSP60 family has been detected in arthritis, multiple sclerosis and diabetes, the vaccination procedure with DNA-HSP65 could induce a cross-reactive immune response that could trigger or worsen these autoimmune diseases. Methods In this investigation was evaluated the effect of a previous vaccination with DNA-HSP65 on diabetes development induced by Streptozotocin (STZ). C57BL/6 mice received three vaccine doses or the corresponding empty vector and were then injected with multiple low doses of STZ. Results DNA-HSP65 vaccination protected mice from STZ induced insulitis and this was associated with higher production of IL-10 in spleen and also in the islets. This protective effect was also concomitant with the appearance of a regulatory cell population in the spleen and a decreased infiltration of the islets by T CD8+ lymphocytes. The vector (DNAv) also determined immunomodulation but its protective effect against insulitis was very discrete. Conclusion The data presented in this study encourages a further investigation in the regulatory potential of the DNA-HSP65 construct. Our findings have important implications for the development of new immune therapy strategies to combat autoimmune diseases.

Multivariate gene expression analysis reveals functional connectivity changes between normal/tumoral prostates

Abstract Background Prostate cancer is a leading cause of death in the male population, therefore, a comprehensive study about the genes and the molecular networks involved in the tumoral prostate process becomes necessary. In order to understand the biological process behind potential biomarkers, we have analyzed a set of 57 cDNA microarrays containing ~25,000 genes. Results Principal Component Analysis (PCA) combined with the Maximum-entropy Linear Discriminant Analysis (MLDA) were applied in order to identify genes with the most discriminative information between normal and tumoral prostatic tissues. Data analysis was carried out using three different approaches, namely: (i) differences in gene expression levels between normal and tumoral conditions from an univariate point of view; (ii) in a multivariate fashion using MLDA; and (iii) with a dependence network approach. Our results show that malignant transformation in the prostatic tissue is more related to functional connectivity changes in their dependence networks than to differential gene expression. The MYLK, KLK2, KLK3, HAN11, LTF, CSRP1 and TGM4 genes presented significant changes in their functional connectivity between normal and tumoral conditions and were also classified as the top seven most informative genes for the prostate cancer genesis process by our discriminant analysis. Moreover, among the identified genes we found classically known biomarkers and genes which are closely related to tumoral prostate, such as KLK3 and KLK2 and several other potential ones. Conclusion We have demonstrated that changes in functional connectivity may be implicit in the biological process which renders some genes more informative to discriminate between normal and tumoral conditions. Using the proposed method, namely, MLDA, in order to analyze the multivariate characteristic of genes, it was possible to capture the changes in dependence networks which are related to cell transformation.