Network enrichment analysis: extension of gene-set enrichment analysis to gene networks

ABSTRACT: BACKGROUND: Gene-set enrichment analyses (GEA or GSEA) are commonly used for biological characterization of an experimental gene-set. This is done by finding known functional categories, such as pathways or Gene Ontology terms, that are over-represented in the experimental set; the assessment is based on an overlap statistic. Rich biological information in terms of gene interaction network is now widely available, but this topological information is not used by GEA, so there is a need for methods that exploit this type of information in high-throughput data analysis. RESULTS: We developed a method of network enrichment analysis (NEA) that extends the overlap statistic in GEA to network links between genes in the experimental set and those in the functional categories. For the crucial step in statistical inference, we developed a fast network randomization algorithm in order to obtain the distribution of any network statistic under the null hypothesis of no association between an experimental gene-set and a functional category. We illustrate the NEA method using gene and protein expression data from a lung cancer study. CONCLUSIONS: The results indicate that the NEA method is more powerful than the traditional GEA, primarily because the relationships between gene sets were more strongly captured by network connectivity rather than by simple overlaps.     

Equilibrium and structure of thallium(III)–ethylenediamine complexes in pyridine solution and in solid

The formation of three [Tl(en)_n]³⁺ complexes (n=1–3) in a pyridine solvent has been established by means of ²⁰⁵Tl and ¹H NMR. Their stepwise stability constants based on concentrations, K_n=[Tl(en)_n ³⁺]/{[Tl(en)_n−1 ³⁺]·[en]}, at 298 K in 0.5 M NaClO₄ ionic medium in pyridine, were calculated from ²⁰⁵Tl NMR integrals: log K₁=7.6±0.7; log K₂=5.2±0.5 and log K₃=2.64±0.05. Linear correlation between both the ²⁰⁵Tl NMR shifts and spin–spin coupling ²⁰⁵Tl–¹H versus the stability constants has been found and discussed. A single crystal with the composition [Tl(en)₃](ClO₄)₃ was synthesized and its structure determined by X-ray diffraction. The Tl³⁺ ion is coordinated by three ethylenediamine ligands via six N-donor atoms in a distorted octahedral fashion.     

Synthesis, cellular uptake and HIV-1 Tat-dependent trans-activation inhibition activity of oligonucleotide analogues disulphide-conjugated to cell-penetrating peptides

Oligonucleotides composed of 2′-O-methyl and locked nucleic acid residues complementary to HIV-1 trans-activation responsive element TAR block Tat-dependent trans-activation in a HeLa cell assay when delivered by cationic lipids. We describe an improved procedure for synthesis and purification under highly denaturing conditions of 5′-disulphide-linked conjugates of 3′-fluorescein labelled oligonucleotides with a range of cell-penetrating peptides and investigate their abilities to enter HeLa cells and block trans-activation. Free uptake of 12mer OMe/LNA oligonucleotide conjugates to Tat (48–58), Penetratin and R₉F₂ was observed in cytosolic compartments of HeLa cells. Uptake of the Tat conjugate was enhanced by N-terminal addition of four Lys or Arg residues or a second Tat peptide. None of the conjugates entered the nucleus or inhibited trans-activation when freely delivered, but inhibition was obtained in the presence of cationic lipids. Nuclear exclusion was seen for free delivery of Tat (48–58), Penetratin and R₉ conjugates of 16mer phosphorothioate OMe oligonucleotide. Uptake into human fibroblast cytosolic compartments was seen for Tat, Penetratin, R₉F₂ and Transportan conjugates. Large enhancements of HeLa cell uptake into cytosolic compartments were seen when free Tat peptide was added to Tat conjugate of 12mer OMe/LNA oligonucleotide or Penetratin peptide to Penetratin conjugate of the same oligonucleotide.     

Genome‐wide association study identifies loci associated with liability to alcohol and drug dependence that is associated with variability in reward‐related ventral striatum activity in African‐ and European‐Americans

Genetic influences on alcohol and drug dependence partially overlap, however, specific loci underlying this overlap remain unclear. We conducted a genome‐wide association study (GWAS) of a phenotype representing alcohol or illicit drug dependence (ANYDEP) among 7291 European‐Americans (EA; 2927 cases) and 3132 African‐Americans (AA: 1315 cases) participating in the family‐based Collaborative Study on the Genetics of Alcoholism. ANYDEP was heritable (h ² in EA = 0.60, AA = 0.37). The AA GWAS identified three regions with genome‐wide significant (GWS; P < 5E‐08) single nucleotide polymorphisms (SNPs) on chromosomes 3 (rs34066662, rs58801820) and 13 (rs75168521, rs78886294), and an insertion‐deletion on chromosome 5 (chr5:141988181). No polymorphisms reached GWS in the EA. One GWS region (chromosome 1: rs1890881) emerged from a trans‐ancestral meta‐analysis (EA + AA) of ANYDEP, and was attributable to alcohol dependence in both samples. Four genes (AA: CRKL, DZIP3, SBK3; EA: P2RX6) and four sets of genes were significantly enriched within biological pathways for hemostasis and signal transduction. GWS signals did not replicate in two independent samples but there was weak evidence for association between rs1890881 and alcohol intake in the UK Biobank. Among 118 AA and 481 EA individuals from the Duke Neurogenetics Study, rs75168521 and rs1890881 genotypes were associated with variability in reward‐related ventral striatum activation. This study identified novel loci for substance dependence and provides preliminary evidence that these variants are also associated with individual differences in neural reward reactivity. Gene discovery efforts in non‐European samples with distinct patterns of substance use may lead to the identification of novel ancestry‐specific genetic markers of risk.     

Proteome analysis of rat liver mitochondria reveals a possible compensatory response to endotoxic shock

Organ failure induced by endotoxic shock has recently been associated with affected mitochondrial function. In this study, effects of in vivo lipopolysaccharide-challenge on protein patterns of rat liver mitochondria in treated animals versus controls were studied by two-dimensional electrophoresis (differential image gel electrophoresis). Significant upregulation was found for ATP-synthase alpha chain and superoxide dismutase [Mn]. Our data suggest that endotoxic shock mediated changes in the mitochondrial proteome contribute to a compensatory reaction (adaptation to endotoxic shock) rather than to a mechanism of cell damage.     

Protein-Protein Interactions in the β-Oxidation Part of the Phenylacetate Utilization Pathway: CRYSTAL STRUCTURE OF THE PaaF-PaaG HYDRATASE-ISOMERASE COMPLEX*

Microbial anaerobic and so-called hybrid pathways for degradation of aromatic compounds contain β-oxidation-like steps. These reactions convert the product of the opening of the aromatic ring to common metabolites. The hybrid phenylacetate degradation pathway is encoded in Escherichia coli by the paa operon containing genes for 10 enzymes. Previously, we have analyzed protein-protein interactions among the enzymes catalyzing the initial oxidation steps in the paa pathway (Grishin, A. M., Ajamian, E., Tao, L., Zhang, L., Menard, R., and Cygler, M. (2011) J. Biol. Chem. 286, 10735–10743). Here we report characterization of interactions between the remaining enzymes of this pathway and show another stable complex, PaaFG, an enoyl-CoA hydratase and enoyl-Coa isomerase, both belonging to the crotonase superfamily. These steps are biochemically similar to the well studied fatty acid β-oxidation, which can be catalyzed by individual monofunctional enzymes, multifunctional enzymes comprising several domains, or enzymatic complexes such as the bacterial fatty acid β-oxidation complex. We have determined the structure of the PaaFG complex and determined that although individually PaaF and PaaG are similar to enzymes from the fatty acid β-oxidation pathway, the structure of the complex is dissimilar from bacterial fatty acid β-oxidation complexes. The PaaFG complex has a four-layered structure composed of homotrimeric discs of PaaF and PaaG. The active sites of PaaF and PaaG are adapted to accept the intermediary components of the Paa pathway, different from those of the fatty acid β-oxidation. The association of PaaF and PaaG into a stable complex might serve to speed up the steps of the pathway following the conversion of phenylacetyl-CoA to a toxic and unstable epoxide-CoA by PaaABCE monooxygenase.     

The ecology of a late Visean forest at Catcraig (East Lothian,Scotland) based on multiproxy study of paleosol and root-casts

The Middle Longcraig Limestone (late Visean) at Catcraig, Scotland, is densely covered with large (about 1 m in diameter) hollows, which are surrounded by numerous fossilized roots. The latter represent the positions of tree-sized plants belonging to a fossil forest. This paper aims to reconstruct its ecology based on paleosol properties, size and spatial distribution of root-casts, rooting system morphology and characteristics of coalified plant remnants. The obtained data are consistent with the supposition that the Visean peat wetland forest of Catcraig represents a Cordaitalean dominated community composed of mono-sized, most probably, even-aged giant trees. The peaty paleosol (Histosol) contains pyrite and siderite, indicating reducing and acidic conditions, whilst the abundant presence of gypsum testifies periods of elevated salinity. The adaptation of trees to such conditions is supported by morphologically evident large root mass (typically more than 1000 individual roots per cast). Good preservation of plant tissues in peat, in addition to other paleosol and rooting system characteristics allow us to conclude that the trees occupied seacoast settings or lagoons, which were not permanently flooded.