Candidate gene prioritization by network analysis of differential expression using machine learning approaches

Background  

Discovering novel disease genes is still challenging for diseases for which no prior knowledge - such as known disease genes or disease-related pathways - is available. Performing genetic studies frequently results in large lists of candidate genes of which only few can be followed up for further investigation. We have recently developed a computational method for constitutional genetic disorders that identifies the most promising candidate genes by replacing prior knowledge by experimental data of differential gene expression between affected and healthy individuals.     

CD4+ natural regulatory T cells prevent experimental cerebral malaria via CTLA-4 when expanded in vivo

Studies in malaria patients indicate that higher frequencies of peripheral blood CD4(+) Foxp3(+) CD25(+) regulatory T (Treg) cells correlate with increased blood parasitemia. This observation implies that Treg cells impair pathogen clearance and thus may be detrimental to the host during infection. In C57BL/6 mice infected with Plasmodium berghei ANKA, depletion of Foxp3(+) cells did not improve parasite control or disease outcome. In contrast, elevating frequencies of natural Treg cells in vivo using IL-2/anti-IL-2 complexes resulted in complete protection against severe disease. This protection was entirely dependent upon Foxp3(+) cells and resulted in lower parasite biomass, impaired antigen-specific CD4(+) T and CD8(+) T cell responses that would normally promote parasite tissue sequestration in this model, and reduced recruitment of conventional T cells to the brain. Furthermore, Foxp3(+) cell-mediated protection was dependent upon CTLA-4 but not IL-10. These data show that T cell-mediated parasite tissue sequestration can be reduced by regulatory T cells in a mouse model of malaria, thereby limiting malaria-induced immune pathology.     

Massively Parallel Signature Sequencing and Bioinformatics Analysis Identifies Up-Regulation of TGFBI and SOX4 in Human Glioblastoma

Background

A comprehensive network-based understanding of molecular pathways abnormally altered in glioblastoma multiforme (GBM) is essential for developing effective therapeutic approaches for this deadly disease.

Methodology/Principal Findings

Applying a next generation sequencing technology, massively parallel signature sequencing (MPSS), we identified a total of 4535 genes that are differentially expressed between normal brain and GBM tissue. The expression changes of three up-regulated genes, CHI3L1, CHI3L2, and FOXM1, and two down-regulated genes, neurogranin and L1CAM, were confirmed by quantitative PCR. Pathway analysis revealed that TGF- β pathway related genes were significantly up-regulated in GBM tumor samples. An integrative pathway analysis of the TGF β signaling network identified two alternative TGF−β signaling pathways mediated by SOX4 (sex determining region Y-box 4) and TGFBI (Transforming growth factor beta induced). Quantitative RT-PCR and immunohistochemistry staining demonstrated that SOX4 and TGFBI expression is elevated in GBM tissues compared with normal brain tissues at both the RNA and protein levels. In vitro functional studies confirmed that TGFBI and SOX4 expression is increased by TGF- β stimulation and decreased by a specific inhibitor of TGF- β receptor 1 kinase.

Conclusions/Significance

Our MPSS database for GBM and normal brain tissues provides a useful resource for the scientific community. The identification of non-SMAD mediated TGF−β signaling pathways acting through SOX4 and TGFBI (GENE ID:7045) in GBM indicates that these alternative pathways should be considered, in addition to the canonical SMAD mediated pathway, in the development of new therapeutic strategies targeting TGF−β signaling in GBM. Finally, the construction of an extended TGF- β signaling network with overlaid gene expression changes between GBM and normal brain extends our understanding of the biology of GBM.     

The structure of borders in a small world

Territorial subdivisions and geographic borders are essential for understanding phenomena in sociology, political science, history, and economics. They influence the interregional flow of information and cross-border trade and affect the diffusion of innovation and technology. However, it is unclear if existing administrative subdivisions that typically evolved decades ago still reflect the most plausible organizational structure of today. The complexity of modern human communication, the ease of long-distance movement, and increased interaction across political borders complicate the operational definition and assessment of geographic borders that optimally reflect the multi-scale nature of today's human connectivity patterns. What border structures emerge directly from the interplay of scales in human interactions is an open question. Based on a massive proxy dataset, we analyze a multi-scale human mobility network and compute effective geographic borders inherent to human mobility patterns in the United States. We propose two computational techniques for extracting these borders and for quantifying their strength. We find that effective borders only partially overlap with existing administrative borders, and show that some of the strongest mobility borders exist in unexpected regions. We show that the observed structures cannot be generated by gravity models for human traffic. Finally, we introduce the concept of link significance that clarifies the observed structure of effective borders. Our approach represents a novel type of quantitative, comparative analysis framework for spatially embedded multi-scale interaction networks in general and may yield important insight into a multitude of spatiotemporal phenomena generated by human activity.     

Biochemical characterization of individual components of the Allochromatium vinosum DsrMKJOP transmembrane complex aids understanding of complex function in vivo

The DsrMKJOP transmembrane complex has a most important function in dissimilatory sulfur metabolism and consists of cytoplasmic, periplasmic, and membrane integral proteins carrying FeS centers and b- and c-type cytochromes as cofactors. In this study, the complex was isolated from the purple sulfur bacterium Allochromatium vinosum and individual components were characterized as recombinant proteins. The two integral membrane proteins DsrM and DsrP were successfully produced in Escherichia coli C43(DE3) and C41(DE3), respectively. DsrM was identified as a diheme cytochrome b, and the two hemes were found to be in low-spin state. Their midpoint redox potentials were determined to be +60 and +110 mV. Although no hemes were predicted for DsrP, it was also clearly identified as a b-type cytochrome. To the best of our knowledge, this is the first time that heme binding has been experimentally proven for a member of the NrfD protein family. Both cytochromes were partly reduced after addition of a menaquinol analogue, suggesting interaction with quinones in vivo. DsrO and DsrK were both experimentally proven to be FeS-containing proteins. In addition, DsrK was shown to be membrane associated, and we propose a monotopic membrane anchoring for this protein. Coelution assays provide support for the proposed interaction of DsrK with the soluble cytoplasmic protein DsrC, which might be its substrate. A model for the function of DsrMKJOP in the purple sulfur bacterium A. vinosum is presented.     

Aminoguanidine administered during the induction of oral tolerance alters the systemic response of the tolerised rats

Herewith we investigated the role of nitric oxide synthase (NOS)-II in the establishment of oral tolerance induced by low antigen dose. To accomplish this, we used a rat model of oral tolerance induced by intragastric administration of low doses of ovalbumin (OVA). NOS-II was inhibited in vivo during the onset of tolerance by intraperitoneal (i.p.) treatment with aminoguanidine (AMG), a selective NOS-II inhibitor. Four experimental groups were generated: (TOL), tolerised rats, receiving OVA but no AMG; (TAG), rats tolerised with OVA and simultaneously receiving AMG i.p.; (CAG), controls treated with AMG but no oral antigen; and (CONT), controls receiving neither OVA nor AMG treatment. The state of oral tolerance was evaluated in all groups by analysing several immune parameters upon subcutaneous administration of OVA in Freund’s complete adjuvant. First, we were able to determine that NOS-II inhibition altered the TH1/TH2 balance in tolerised rats, driving the TH2 anti-OVA response in TOL rats towards TH1 in TAG animals, which showed enhanced delayed hypersensitivity responses. Second, splenocyte cultures from TAG rats showed lower levels of IL-10 production compared to TOL samples as determined by ELISA analysis. Last, we detected the presence of a functional distinct Tr1 regulatory T cell population in spleen samples recovered from TAG animals. Contrary to what happened with TOL Tr1 cells, the levels of Tr1 cells in TAG samples were modified by in vitro stimulation with OVA.All together, these data indicate a preponderant role for NOS-II in the process of oral tolerance induced by low antigen dose.     

Inferring meaningful pathways in weighted metabolic networks

An approach is presented for computing meaningful pathways in the network of small molecule metabolism comprising the chemical reactions characterized in all organisms. The metabolic network is described as a weighted graph in which all the compounds are included, but each compound is assigned a weight equal to the number of reactions in which it participates. Path finding is performed in this graph by searching for one or more paths with lowest weight. Performance is evaluated systematically by computing paths between the first and last reactions in annotated metabolic pathways, and comparing the intermediate reactions in the computed pathways to those in the annotated ones. For the sake of comparison, paths are computed also in the un-weighted raw (all compounds and reactions) and filtered (highly connected pool metabolites removed) metabolic graphs, respectively. The correspondence between the computed and annotated pathways is very poor (<30%) in the raw graph; increasing to approximately 65% in the filtered graph; reaching approximately 85% in the weighted graph. Considering the best-matching path among the five lightest paths increases the correspondence to 92%, on average. We then show that the average distance between pairs of metabolites is significantly larger in the weighted graph than in the raw unfiltered graph, suggesting that the small-world properties previously reported for metabolic networks probably result from irrelevant shortcuts through pool metabolites. In addition, we provide evidence that the length of the shortest path in the weighted graph represents a valid measure of the "metabolic distance" between enzymes. We suggest that the success of our simplistic approach is rooted in the high degree of specificity of the reactions in metabolic pathways, presumably reflecting thermodynamic constraints operating in these pathways. We expect our approach to find useful applications in inferring metabolic pathways in newly sequenced genomes.     

Asymmetric expression of the BMP antagonists chordin and gremlin in the sea anemone Nematostella vectensis: implications for the evolution of axial patterning

The evolutionary origin of the anterior-posterior and the dorsoventral body axes of Bilateria is a long-standing question. It is unclear how the main body axis of Cnidaria, the sister group to the Bilateria, is related to the two body axes of Bilateria. The conserved antagonism between two secreted factors, BMP2/4 (Dpp in Drosophila) and its antagonist Chordin (Short gastrulation in Drosophila) is a crucial component in the establishment of the dorsoventral body axis of Bilateria and could therefore provide important insight into the evolutionary origin of bilaterian axes. Here, we cloned and characterized two BMP ligands, dpp and GDF5-like as well as two secreted antagonists, chordin and gremlin, from the basal cnidarian Nematostella vectensis. Injection experiments in zebrafish show that the ventralizing activity of NvDpp mRNA is counteracted by NvGremlin and NvChordin, suggesting that Gremlin and Chordin proteins can function as endogenous antagonists of NvDpp. Expression analysis during embryonic and larval development of Nematostella reveals asymmetric expression of all four genes along both the oral-aboral body axis and along an axis perpendicular to this one, the directive axis. Unexpectedly, NvDpp and NvChordin show complex and overlapping expression on the same side of the embryo, whereas NvGDF5-like and NvGremlin are both expressed on the opposite side. Yet, the two pairs of ligands and antagonists only partially overlap, suggesting complex gradients of BMP activity along the directive axis but also along the oral-aboral axis. We conclude that a molecular interaction between BMP-like molecules and their secreted antagonists was already employed in the common ancestor of Cnidaria and Bilateria to create axial asymmetries, but that there is no simple relationship between the oral-aboral body axis of Nematostella and one particular body axis of Bilateria.     

Biochemical fates of alpha hemoglobin bound to alpha hemoglobin-stabilizing protein AHSP

Alpha hemoglobin-stabilizing protein (AHSP) is an erythroid protein that binds free alpha hemoglobin (alphaHb) to maintain its structure and limit its pro-oxidant activity. Prior studies have defined two different alphaHb.AHSP complexes. Binding of AHSP to Fe(II) alphaHb induces an unusual configuration in which the F helix of the globin becomes disordered and the heme ring becomes solvent-exposed. Over time, this intermediate oxidizes to form a stable hemichrome in which the proximal (F8) and distal (E7) histidines are coordinated to the heme iron atom. The addition of betaHb to either Fe(II) or Fe(III) alphaHb.AHSP displaces AHSP to generate tetrameric (alpha(2)beta(2)) HbA species. The biochemical properties and in vivo significance of the two alphaHb.AHSP complexes are poorly understood. Here we show that Fe(III) alphaHb.AHSP forms from auto-oxidation of oxygenated alphaHb bound to AHSP and that this process is greatly accelerated at physiologic temperature and oxygen pressures. In contrast to free Fe(III) alphaHb hemichromes, AHSP-bound Fe(III) alphaHb does not precipitate and can be recycled into functional HbA. This requires enzymatic reduction of AHSP-bound alphaHb, either prior to or after extraction by beta subunits. In contrast, reaction of Fe(II) alphaHb-AHSP with betaHb generates functional HbA directly. Our findings support a model in which AHSP can either stabilize alphaHb transiently en route to HbA formation during normal erythropoiesis or convert excessive free alphaHb into a more chemically inert state from which recovery of alphaHb is possible by redox cycling.     

Is the ApoE polymorphism associated with dilated cardiomyopathy?

Apolipoprotein E (ApoE) is 34 kDa protein involved in the modulation of cholesterol transport and homeostasis. Polymorphism of the ApoE gene has been implicated in many chronic cardiovascular and neuronal diseases. ApoE epsilon4 allele has been reported to be associated with increased risk of cardiovascular diseases such as myocardial infarction, hypertension, coronary heart disease, etc. Fifty patients with the end-stage dilated cardiomyopathy (DCM) and advanced congestive heart failure were examined in our study. For evaluation of ApoE polymorphism, novel approach of fast screening of ApoE gene polymorphism by combination of PCR and blotting (CVD StripAssay) was used. Individual genotypes were correlated with basic cardiologic clinical parameters. The reported frequency of this allele in Caucasian population is 14.7 %. Our results showed that in patients with DCM frequency of the ApoE epsilon4 allele is 40 %. Frequency of the genotype epsilon2/4 was 58 % and epsilon3/4 was 22 %. Comparison with control Caucasian groups monitored by others clearly revealed that frequency of epsilon4 alelle is increased in patients with advanced stages of DCM. This observation suggests association of ApoE polymorphism with severe form of DCM. Physiological consequences of this observation remain to be clarified.