A Relative Variation-Based Method to Unraveling Gene Regulatory Networks

Gene regulatory network (GRN) reconstruction is essential in understanding the functioning and pathology of a biological system. Extensive models and algorithms have been developed to unravel a GRN. The DREAM project aims to clarify both advantages and disadvantages of these methods from an application viewpoint. An interesting yet surprising observation is that compared with complicated methods like those based on nonlinear differential equations, etc., methods based on a simple statistics, such as the so-called -score, usually perform better. A fundamental problem with the -score, however, is that direct and indirect regulations can not be easily distinguished. To overcome this drawback, a relative expression level variation (RELV) based GRN inference algorithm is suggested in this paper, which consists of three major steps. Firstly, on the basis of wild type and single gene knockout/knockdown experimental data, the magnitude of RELV of a gene is estimated. Secondly, probability for the existence of a direct regulation from a perturbed gene to a measured gene is estimated, which is further utilized to estimate whether a gene can be regulated by other genes. Finally, the normalized RELVs are modified to make genes with an estimated zero in-degree have smaller RELVs in magnitude than the other genes, which is used afterwards in queuing possibilities of the existence of direct regulations among genes and therefore leads to an estimate on the GRN topology. This method can in principle avoid the so-called cascade errors under certain situations. Computational results with the Size 100 sub-challenges of DREAM3 and DREAM4 show that, compared with the -score based method, prediction performances can be substantially improved, especially the AUPR specification. Moreover, it can even outperform the best team of both DREAM3 and DREAM4. Furthermore, the high precision of the obtained most reliable predictions shows that the suggested algorithm may be very helpful in guiding biological experiment designs

From Africa to Europe and back: refugia and range shifts cause high genetic differentiation in the Marbled White butterfly Melanargia galathea

<p>Abstract</p> <p>Background</p> <p>The glacial-interglacial oscillations caused severe range modifications of biota. Thermophilic species became extinct in the North and survived in southern retreats, e.g. the Mediterranean Basin. These repeated extinction and (re)colonisation events led to long-term isolation and intermixing of populations and thus resulted in strong genetic imprints in many European species therefore being composed of several genetic lineages. To better understand these cycles of repeated expansion and retraction, we selected the Marbled White butterfly <it>Melanargia galathea</it>. Fourty-one populations scattered over Europe and the Maghreb and one population of the sibling taxon <it>M. lachesis </it>were analysed using allozyme electrophoresis.</p> <p>Results</p> <p>We obtained seven distinct lineages applying neighbour joining and STRUCTURE analyses: (i) Morocco, (ii) Tunisia, (iii) Sicily, (iv) Italy and southern France, (v) eastern Balkans extending to Central Europe, (vi) western Balkans with western Carpathian Basin as well as (vii) south-western Alps. The hierarchy of these splits is well matching the chronology of glacial and interglacial cycles since the Günz ice age starting with an initial split between the <it>galathea </it>group in North Africa and the <it>lachesis </it>group in Iberia. These genetic structures were compared with past distribution patterns during the last glacial stage calculated with distribution models.</p> <p>Conclusions</p> <p>Both methods suggest climatically suitable areas in the Maghreb and the southern European peninsulas with distinct refugia during the last glacial period and underpin strong range expansions to the North during the Postglacial. However, the allozyme patterns reveal biogeographical structures not detected by distribution modelling as two distinct refugia in the Maghreb, two or more distinct refugia at the Balkans and a close link between the eastern Maghreb and Sicily. Furthermore, the genetically highly diverse western Maghreb might have acted as source or speciation centre of this taxon, while the eastern, genetically impoverished Maghreb population might result from a relatively recent recolonisation from Europe via Sicily.</p

Coordinated Expression Domains in Mammalian Genomes

Gene order in eukaryotic genomes is not random. Genes showing similar expression (coexpression) patterns are often clustered along the genome. The goal of this study is to characterize coexpression clustering in mammalian genomes and to investigate the underlying mechanisms.We detect clustering of coexpressed genes across multiple scales, from neighboring genes to chromosomal domains that span tens of megabases and, in some cases, entire chromosomes. Coexpression domains may be positively or negatively correlated with other domains, within and between chromosomes. We find that long-range expression domains are associated with gene density, which in turn is related to physical organization of the chromosomes within the nucleus. We show that gene expression changes between healthy and diseased tissue samples occur in a gene density-dependent manner.We demonstrate that coexpression domains exist across multiple scales. We identify potential mechanisms for short-range as well as long-range coexpression domains. We provide evidence that the three-dimensional architecture of the chromosomes may underlie long-range coexpression domains. Chromosome territory reorganization may play a role in common human diseases such as Alzheimer's disease and psoriasis

MICU2, a Paralog of MICU1, Resides within the Mitochondrial Uniporter Complex to Regulate Calcium Handling

Mitochondrial calcium uptake is present in nearly all vertebrate tissues and is believed to be critical in shaping calcium signaling, regulating ATP synthesis and controlling cell death. Calcium uptake occurs through a channel called the uniporter that resides in the inner mitochondrial membrane. Recently, we used comparative genomics to identify MICU1 and MCU as the key regulatory and putative pore-forming subunits of this channel, respectively. Using bioinformatics, we now report that the human genome encodes two additional paralogs of MICU1, which we call MICU2 and MICU3, each of which likely arose by gene duplication and exhibits distinct patterns of organ expression. We demonstrate that MICU1 and MICU2 are expressed in HeLa and HEK293T cells, and provide multiple lines of biochemical evidence that MCU, MICU1 and MICU2 reside within a complex and cross-stabilize each other's protein expression in a cell-type dependent manner. Using in vivo RNAi technology to silence MICU1, MICU2 or both proteins in mouse liver, we observe an additive impairment in calcium handling without adversely impacting mitochondrial respiration or membrane potential. The results identify MICU2 as a new component of the uniporter complex that may contribute to the tissue-specific regulation of this channel.National Institutes of Health (U.S.) (GM0077465)National Institutes of Health (U.S.) (DK080261

Distribution of misfolded prion protein seeding activity alone does not predict regions of neurodegeneration

Protein misfolding is common across many neurodegenerative diseases, with misfolded proteins acting as seeds for "prion-like" conversion of normally folded protein to abnormal conformations. A central hypothesis is that misfolded protein accumulation, spread and distribution is restricted to specific neuronal populations of the central nervous system and thus predict regions of neurodegeneration. We examined this hypothesis using a highly sensitive assay system for detection of misfolded protein seeds in a murine model of prion disease. Misfolded prion protein seeds were observed widespread throughout the brain accumulating in all brain regions examined irrespective of neurodegeneration. Importantly neither time of exposure nor amount of misfolded protein seeds present determined regions of neurodegeneration. We further demonstrate two distinct microglia responses in prion infected brains, a 11 novel homeostatic response in all regions and an innate immune response restricted to sites of 12 neurodegeneration. Therefore accumulation of misfolded prion protein alone does not define targeting 13 of neurodegeneration which instead results only when misfolded prion protein accompanies a specific 14 innate immune response

Models Analyses for Allelopathic Effects of Chicory at Equivalent Coupling of Nitrogen Supply and pH Level on F. arundinacea, T. repens and M. sativa

Alllelopathic potential of chicory was investigated by evaluating its effect on seed germination, soluble sugar, malondialdehyde (MDA) and the chlorophyll content of three target plants species (Festuca arundinacea, Trifolium repens and Medicago sativa). The secretion of allelochemicals was regulated by keeping the donor plant (chicory) separate from the three target plant species and using different pH and nitrogen levels. Leachates from donor pots with different pH levels and nitrogen concentrations continuously irrigated the target pots containing the seedlings. The allelopathic effects of the chicory at equivalent coupling of nitrogen supply and pH level on the three target plants species were explored via models analyses. The results suggested a positive effect of nitrogen supply and pH level on allelochemical secretion from chicory plants. The nitrogen supply and pH level were located at a rectangular area defined by 149 to 168 mg/l nitrogen supply combining 4.95 to 7.0 pH value and point located at nitrogen supply 177 mg/l, pH 6.33 when they were in equivalent coupling effects; whereas the inhibitory effects of equivalent coupling nitrogen supply and pH level were located at rectangular area defined by 125 to 131 mg/l nitrogen supply combining 6.71 to 6.88 pH value and two points respectively located at nitrogen supply 180 mg/l with pH 6.38 and nitrogen supply 166 mg/l with pH 7.59. Aqueous extracts of chicory fleshy roots and leaves accompanied by treatment at different sand pH values and nitrogen concentrations influenced germination, seedling growth, soluble sugar, MDA and chlorophyll of F. arundinacea, T. repens and M. sativa. Additionally, we determined the phenolics contents of root and leaf aqueous extracts, which were 0.104% and 0.044% on average, respectively

A gene expression atlas of the domestic pig

<p>Abstract</p> <p>Background</p> <p>This work describes the first genome-wide analysis of the transcriptional landscape of the pig. A new porcine Affymetrix expression array was designed in order to provide comprehensive coverage of the known pig transcriptome. The new array was used to generate a genome-wide expression atlas of pig tissues derived from 62 tissue/cell types. These data were subjected to network correlation analysis and clustering.</p> <p>Results</p> <p>The analysis presented here provides a detailed functional clustering of the pig transcriptome where transcripts are grouped according to their expression pattern, so one can infer the function of an uncharacterized gene from the company it keeps and the locations in which it is expressed. We describe the overall transcriptional signatures present in the tissue atlas, where possible assigning those signatures to specific cell populations or pathways. In particular, we discuss the expression signatures associated with the gastrointestinal tract, an organ that was sampled at 15 sites along its length and whose biology in the pig is similar to human. We identify sets of genes that define specialized cellular compartments and region-specific digestive functions. Finally, we performed a network analysis of the transcription factors expressed in the gastrointestinal tract and demonstrate how they sub-divide into functional groups that may control cellular gastrointestinal development.</p> <p>Conclusions</p> <p>As an important livestock animal with a physiology that is more similar than mouse to man, we provide a major new resource for understanding gene expression with respect to the known physiology of mammalian tissues and cells. The data and analyses are available on the websites <url>http://biogps.org and http://www.macrophages.com/pig-atlas</url>.</p

Bioinformatic analyses identifies novel protein-coding pharmacogenomic markers associated with paclitaxel sensitivity in NCI60 cancer cell lines

<p>Abstract</p> <p>Background</p> <p>Paclitaxel is a microtubule-stabilizing drug that has been commonly used in treating cancer. Due to genetic heterogeneity within patient populations, therapeutic response rates often vary. Here we used the NCI60 panel to identify SNPs associated with paclitaxel sensitivity. Using the panel's GI50 response data available from Developmental Therapeutics Program, cell lines were categorized as either sensitive or resistant. PLINK software was used to perform a genome-wide association analysis of the cellular response to paclitaxel with the panel's SNP-genotype data on the Affymetrix 125 k SNP array. FastSNP software helped predict each SNP's potential impact on their gene product. mRNA expression differences between sensitive and resistant cell lines was examined using data from BioGPS. Using Haploview software, we investigated for haplotypes that were more strongly associated with the cellular response to paclitaxel. Ingenuity Pathway Analysis software helped us understand how our identified genes may alter the cellular response to paclitaxel.</p> <p>Results</p> <p>43 SNPs were found significantly associated (FDR < 0.005) with paclitaxel response, with 10 belonging to protein-coding genes (<it>CFTR</it>, <it>ROBO1</it>, <it>PTPRD</it>, <it>BTBD12</it>, <it>DCT</it>, <it>SNTG1</it>, <it>SGCD</it>, <it>LPHN2</it>, <it>GRIK1</it>, <it>ZNF607</it>). SNPs in <it>GRIK1</it>, <it>DCT</it>, <it>SGCD </it>and <it>CFTR </it>were predicted to be intronic enhancers, altering gene expression, while SNPs in <it>ZNF607 </it>and <it>BTBD12 </it>cause conservative missense mutations. mRNA expression analysis supported these findings as <it>GRIK1</it>, <it>DCT</it>, <it>SNTG1</it>, <it>SGCD </it>and <it>CFTR </it>showed significantly (p < 0.05) increased expression among sensitive cell lines. Haplotypes found in <it>GRIK1, SGCD, ROBO1, LPHN2</it>, and <it>PTPRD </it>were more strongly associated with response than their individual SNPs.</p> <p>Conclusions</p> <p>Our study has taken advantage of available genotypic data and its integration with drug response data obtained from the NCI60 panel. We identified 10 SNPs located within protein-coding genes that were not previously shown to be associated with paclitaxel response. As only five genes showed differential mRNA expression, the remainder would not have been detected solely based on expression data. The identified haplotypes highlight the role of utilizing SNP combinations within genomic loci of interest to improve the risk determination associated with drug response. These genetic variants represent promising biomarkers for predicting paclitaxel response and may play a significant role in the cellular response to paclitaxel.</p