Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures.

INTRODUCTION: Breast cancer subtyping and prognosis have been studied extensively by gene expression profiling, resulting in disparate signatures with little overlap in their constituent genes. Although a previous study demonstrated a prognostic concordance among gene expression signatures, it was limited to only one dataset and did not fully elucidate how the different genes were related to one another nor did it examine the contribution of well-known biological processes of breast cancer tumorigenesis to their prognostic performance. METHOD: To address the above issues and to further validate these initial findings, we performed the largest meta-analysis of publicly available breast cancer gene expression and clinical data, which are comprised of 2,833 breast tumors. Gene coexpression modules of three key biological processes in breast cancer (namely, proliferation, estrogen receptor [ER], and HER2 signaling) were used to dissect the role of constituent genes of nine prognostic signatures. RESULTS: Using a meta-analytical approach, we consolidated the signatures associated with ER signaling, ERBB2 amplification, and proliferation. Previously published expression-based nomenclature of breast cancer 'intrinsic' subtypes can be mapped to the three modules, namely, the ER-/HER2- (basal-like), the HER2+ (HER2-like), and the low- and high-proliferation ER+/HER2- subtypes (luminal A and B). We showed that all nine prognostic signatures exhibited a similar prognostic performance in the entire dataset. Their prognostic abilities are due mostly to the detection of proliferation activity. Although ER- status (basal-like) and ERBB2+ expression status correspond to bad outcome, they seem to act through elevated expression of proliferation genes and thus contain only indirect information about prognosis. Clinical variables measuring the extent of tumor progression, such as tumor size and nodal status, still add independent prognostic information to proliferation genes. CONCLUSION: This meta-analysis unifies various results of previous gene expression studies in breast cancer. It reveals connections between traditional prognostic factors, expression-based subtyping, and prognostic signatures, highlighting the important role of proliferation in breast cancer prognosis.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa.

PURPOSE: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. METHODS: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. RESULTS: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. CONCLUSIONS: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Alterations of the 5'untranslated region of SLC16A12 lead to age-related cataract.

PURPOSE. Knowledge of genetic factors predisposing to age-related cataract is very limited. The aim of this study was to identify DNA sequences that either lead to or predispose for this disease. METHODS. The candidate gene SLC16A12, which encodes a solute carrier of the monocarboxylate transporter family, was sequenced in 484 patients with cataract (134 with juvenile cataract, 350 with age-related cataract) and 190 control subjects. Expression studies included luciferase reporter assay and RT-PCR experiments. RESULTS. One patient with age-related cataract showed a novel heterozygous mutation (c.-17A>G) in the 5'untranslated region (5'UTR). This mutation is in cis with the minor G-allele of the single nucleotide polymorphism (SNP) rs3740030 (c.-42T/G), also within the 5'UTR. Using a luciferase reporter assay system, a construct with the patient's haplotype caused a significant upregulation of luciferase activity. In comparison, the SNP G-allele alone promoted less activity, but that amount was still significantly higher than the amount of the common T-allele. Analysis of SLC16A12 transcripts in surrogate tissue demonstrated striking allele-specific differences causing 5'UTR heterogeneity with respect to sequence and quantity. These differences in gene expression were mirrored in an allele-specific predisposition to age-related cataract, as determined in a Swiss population (odds ratio approximately 2.2; confidence intervals, 1.23-4.3). CONCLUSIONS. The monocarboxylate transporter SLC16A12 may contribute to age-related cataract. Sequences within the 5'UTR modulate translational efficiency with pathogenic consequences.

Probing the role of point mutations in the cyp51A gene from Aspergillus fumigatus in the model yeast Saccharomyces cerevisiae.

Azole-resistant strains of Aspergillus fumigatus have been detected and the underlying molecular mechanisms of resistance characterized. Point mutations in the cyp51A gene have been proved to be related to azole resistance in A. fumigatus clinical strains and with different resistance profiles depending on the amino acid change (G54E, G54V, G54R, G54W, M220V, M220K, M220T, M220I). The aim of this work was to express A. fumigatus cyp51A genes in the yeast Saccharomyces cerevisiae in order to better assess the contribution of each independent amino acid substitution to resistance. A tetracycline regulatable system allowing repression of the endogenous essential ERG11 gene was used. The expression of Aspergillus cyp51A alleles could efficiently restore the absence of ERG11 in S. cerevisiae. In general, S. cerevisiae clones expressing. A. fumigatus cyp51A alleles from azole-resistant isolates showed higher MICs to all azoles tested than those expressing alleles from susceptible isolates. The azole susceptibility profiles obtained in S. cerevisiae upon expression of specific cyp51A alleles recapitulated susceptibility profiles observed from their A. fumigatus origins. In conclusion this work supports the concept that characteristics of specific A. fumigatus cyp51A alleles could be investigated in the heterologous host S. cerevisiae.

Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

Macrophage migration inhibitory factor is critically involved in basal and fluoxetine-stimulated adult hippocampal cell proliferation and in anxiety, depression, and memory-related behaviors.

Intensive research is devoted to unravel the neurobiological mechanisms mediating adult hippocampal neurogenesis, its regulation by antidepressants, and its behavioral consequences. Macrophage migration inhibitory factor (MIF) is a pro-inflammatory cytokine that is expressed in the CNS, where its function is unknown. Here, we show, for the first time, the relevance of MIF expression for adult hippocampal neurogenesis. We identify MIF expression in neurogenic cells (in stem cells, cells undergoing proliferation, and in newly proliferated cells undergoing maturation) in the subgranular zone of the rodent dentate gyrus. A causal function for MIF in cell proliferation was shown using genetic (MIF gene deletion) and pharmacological (treatment with the MIF antagonist Iso-1) approaches. Behaviorally, genetic deletion of MIF resulted in increased anxiety- and depression-like behaviors, as well as of impaired hippocampus-dependent memory. Together, our studies provide evidence supporting a pivotal function for MIF in both basal and antidepressant-stimulated adult hippocampal cell proliferation. Moreover, loss of MIF results in a behavioral phenotype that, to a large extent, corresponds with alterations predicted to arise from reduced hippocampal neurogenesis. These findings underscore MIF as a potentially relevant molecular target for the development of treatments linked to deficits in neurogenesis, as well as to problems related to anxiety, depression, and cognition.

Comparative genomics suggests that the human pathogenic fungus Pneumocystis jirovecii acquired obligate biotrophy through gene loss.

Pneumocystis jirovecii is a fungal parasite that colonizes specifically humans and turns into an opportunistic pathogen in immunodeficient individuals. The fungus is able to reproduce extracellularly in host lungs without eliciting massive cellular death. The molecular mechanisms that govern this process are poorly understood, in part because of the lack of an in vitro culture system for Pneumocystis spp. In this study, we explored the origin and evolution of the putative biotrophy of P. jirovecii through comparative genomics and reconstruction of ancestral gene repertoires. We used the maximum parsimony method and genomes of related fungi of the Taphrinomycotina subphylum. Our results suggest that the last common ancestor of Pneumocystis spp. lost 2,324 genes in relation to the acquisition of obligate biotrophy. These losses may result from neutral drift and affect the biosyntheses of amino acids and thiamine, the assimilation of inorganic nitrogen and sulfur, and the catabolism of purines. In addition, P. jirovecii shows a reduced panel of lytic proteases and has lost the RNA interference machinery, which might contribute to its genome plasticity. Together with other characteristics, that is, a sex life cycle within the host, the absence of massive destruction of host cells, difficult culturing, and the lack of virulence factors, these gene losses constitute a unique combination of characteristics which are hallmarks of both obligate biotrophs and animal parasites. These findings suggest that Pneumocystis spp. should be considered as the first described obligate biotrophs of animals, whose evolution has been marked by gene losses.

In vitro characterization of immune-related properties of human fetal bone cells for potential tissue engineering applications.

We describe herein some immunological properties of human fetal bone cells recently tested for bone tissue-engineering applications. Adult mesenchymal stem cells (MSCs) and osteoblasts were included in the study for comparison. Surface markers involved in bone metabolism and immune recognition were analyzed using flow cytometry before and after differentiation or treatment with cytokines. Immunomodulatory properties were studied on activated peripheral blood mononuclear cells (PBMCs). The immuno-profile of fetal bone cells was further investigated at the gene expression level. Fetal bone cells and adult MSCs were positive for Stro-1, alkaline phosphatase, CD10, CD44, CD54, and beta2-microglobulin, but human leukocyte antigen (HLA)-I and CD80 were less present than on adult osteoblasts. All cells were negative for HLA-II. Treatment with recombinant human interferon gamma increased the presence of HLA-I in adult cells much more than in fetal cells. In the presence of activated PBMCs, fetal cells had antiproliferative effects, although with patterns not always comparable with those of adult MSCs and osteoblasts. Because of the immunological profile, and with their more-differentiated phenotype than of stem cells, fetal bone cells present an interesting potential for allogeneic cell source in tissue-engineering applications.

Relationship between estrogen receptor alpha location and gene induction reveals the importance of downstream sites and cofactors.

BACKGROUND: To understand cancer-related modifications to transcriptional programs requires detailed knowledge about the activation of signal-transduction pathways and gene expression programs. To investigate the mechanisms of target gene regulation by human estrogen receptor alpha (hERalpha), we combine extensive location and expression datasets with genomic sequence analysis. In particular, we study the influence of patterns of DNA occupancy by hERalpha on expression phenotypes. RESULTS: We find that strong ChIP-chip sites co-localize with strong hERalpha consensus sites and detect nucleotide bias near hERalpha sites. The localization of ChIP-chip sites relative to annotated genes shows that weak sites are enriched near transcription start sites, while stronger sites show no positional bias. Assessing the relationship between binding configurations and expression phenotypes, we find binding sites downstream of the transcription start site (TSS) to be equally good or better predictors of hERalpha-mediated expression as upstream sites. The study of FOX and SP1 cofactor sites near hERalpha ChIP sites shows that induced genes frequently have FOX or SP1 sites. Finally we integrate these multiple datasets to define a high confidence set of primary hERalpha target genes. CONCLUSION: Our results support the model of long-range interactions of hERalpha with the promoter-bound cofactor SP1 residing at the promoter of hERalpha target genes. FOX motifs co-occur with hERalpha motifs along responsive genes. Importantly we show that the spatial arrangement of sites near the start sites and within the full transcript is important in determining response to estrogen signaling.

Hypermethylation-mediated regulation of CD44 gene expression in human neuroblastoma

The CD44 adhesion receptor is silenced in highly malignant neuroblastomas (NBs) with MYCN amplification. Because its functional expression is associated with decreased tumorigenic properties, CD44 behaves as a tumor suppressor gene in NB and other cancers. Given that the precise mechanisms responsible for CD44 silencing are not elucidated, we investigated whether CD44 expression could be regulated by DNA hypermethylation. The methylation status of CD44 gene promoter and exon 1 regions was analyzed in 12 NB cell lines and 21 clinical samples after bisulfite genomic modification, followed by PCR and single-strand conformation polymorphism analysis and genomic sequencing. The results showed that almost all CD44-negative cell lines displayed hypermethylation in both regions, whereas all CD44-expressing cell lines were unmethylated. These observations correlated with the ability to restore CD44 mRNA and protein expression by treatment of CD44-negative cells with the 5-aza-2'-deoxycytidine demethylating agent. In contrast, no CD44 gene hypermethylation could be detected in 21 NB clinical samples of different stages, irrespective of CD44 expression. Although our results suggest that aberrant methylation of promoter and exon 1 regions is involved in CD44 silencing in NB cell lines, they also indicate that methylation of unidentified regulatory sequences or methylation-independent mechanisms also control the expression of CD44 in primary NB tumors and cell lines. We therefore conclude that CD44 silencing is controlled by complex and tumor cell-specific processes, including gene hypermethylation. Further investigation of other mechanisms and genes involved in CD44 regulation will be needed before demethylation-mediated reactivation of the CD44 gene can be considered as therapeutic strategy for neuroblastoma and perhaps other related cancers.

Perinuclear Mlp proteins down-regulate gene expression in response to a defect in mRNA export

Résumé Une caractéristique des cellules eucaryotes est le confinement du matériel génétique (ADN/DNA) dans le noyau. Pour décoder cette information, un ARN messager (mRNA) est d'abord transcrit sous forme d'un ARN prémessager (pré-mRNA). Ce-dernier doit subir plusieurs étapes de maturation pour aboutir à une particule ribonucléoprotéique (mRNP) qui sera exportée vers le cytoplasme et traduite en protéine. La protéine de levure Mex67p et son homologue humain TAP sont des récepteurs d'export médiant la translocation du mRNP au travers des complexes du pore nucléaire (NPC). Mex67p/TAP ne se lient pas directement au mRNA, mais nécessitent la présence de protéines adaptatrices, telles que Yra1p et son homologue humain REF1. Afin d'identifier de nouveaux facteurs impliqués dans l'export des mRNPs ou de nouvelles fonctions pour Yra1p, nous avons effectué un crible génétique avec un mutant thermosensible de Yra1p, GFP-yra 1 -8. Ce mutant présente un défaut d'export des mRNAs et une diminution des niveaux de transcrits du gène rapporteur LacZ ainsi que de certains transcrits endogènes. Nous avons trouvé que la perte de Mlp2p, ou d'une protéine hautement similaire, Mlp1p, restaure la croissance du mutant GFP-yra1-8 à température restrictive. Mlp1p et Mlp2p sont des protéines nucléaires, dont l'homologue humain est TPR. Les Mlp (myosin¬like proteins) ainsi que TPR forment des structures filamenteuses ancrées aux NPC. Bien que la fonction des Mlp ne soit pas clairement définie, un rôle dans la biogenèse et la surveillance des mRNPs a été récemment proposé. Notre étude montre que la perte des Mlp, non seulement restaure la croissance de GFP-yra1-8, mais augmente aussi les niveaux des transcrits LacZ et facilite leur apparition dans le cytoplasme. Des expériences d'immunoprécipitations de la chromatine révèlent que Mlp2p diminue le taux de synthèse du transcrit LacZ dans GFP-yra1-8. Des analyses du transcriptome montrent que Mlp2p réduit aussi les niveaux d'une population de transcrits endogènes dans le mutant. Finalement, des localisations in situ suggèrent que la transcription du rapporteur LacZ a lieu à la périphérie du noyau, à proximité des Mlp. Ainsi, les protéines Mlp pourraient préférentiellement diminuer la transcription de gènes exprimés à la périphérie nucléaire. Nous montrons aussi que Yra1p interagit génétiquement avec Nab2p une protéine liée au mRNA et impliquée dans son export, mais non avec d'autres protéines également impliquées dans l'export des mRNAs. Les résultats obtenus soutiennent un modèle où les protéines Yra1p et Nab2p sont nécessaires à l'arrimage des mRNPs sur la plate-forme des Mlp. Si ces signaux manquent ou sont défectueux, les mRNPs ne peuvent pas poursuivre leur trajet vers le canal central du NPC. Ce bloc induirait par la suite une diminution de la transcription d'une population de gènes potentiellement localisée à la périphérie nucléaire. Dans son ensemble, cette étude suggère que les protéines Mlp établissent un lien entre la transcription de certains mRNAs et leur export au travers du pore nucléaire. Summary A hallmark of the eukaryotic cell is the packaging of DNA in the nucleus. To decode the genetic information, a messenger RNA (mRNA) is first synthesized as a pre-mRNA molecule, which undergoes different maturation steps resulting in an mRNP (messenger RNA ribonucleoprotein), which can be actively transported to the cytoplasm and translated into a protein. Yeast Mex67p and its human homologue TAP are export receptors mediating mRNP translocation through the nuclear pore complex (NPC). The recruitment of Mex67p/TAP to mRNA is mediated by mRNA export adaptors of the evolutionarily conserved REF (RNA and Export Factor binding) family: yeast Yra1p and human REF1. To uncover new functions of Yra1p or new factors implicated in mRNA export, we performed a genetic screen with a themiosensitive (ts) yra1 mutant, GFP-yra1-8. This mutant exhibits mRNA export defects and a decrease in the levels of LacZ reporter and certain endogenous transcripts. We found that the loss of Mlp2p, or the related Mlp1p protein, substantially rescues the growth defect of the GFP-yra1 -8 mutant. Mlp1p and M1p2p are large non-essential proteins, homologous to human TPR, proposed to form intra-nuclear filamentous structures anchored at the NPC. Their role is not clearly defined, but they have been implicated in mRNP biogenesis and surveillance. Our study shows that loss of Mlp proteins not only restores growth of GFP-yra1-8, but also rescues LacZ mRNA levels and increases their appearance in the cytoplasm. Chromatin immunoprecipitation and pulse chase experiments indicate that Mlp2p down-regulates LacZ mRNA synthesis in GFP-yra1-8. DNA micro- array analyses reveal that Mlp2p also reduces the levels of a subset of cellular transcripts in the yra1 mutant strain. In situ localizations suggest that LacZ transcription occurs at the nuclear periphery, in close proximity to Mlp proteins. Thus, Mlp proteins may preferentially down-regulate genes expressed at the nuclear periphery. Finally, we show that Yra1p genetically interacts with the shuttling mRNA-binding protein Nab2p and that loss of Mlp proteins rescues the growth defect of yra1 and nab2, but not other mRNA export mutants. The data support a model in which Nab2p and Yra1p are required for rnRNP docking to the Mlp platform. Lack of these signals prevents mRNPs from crossing the Mlp gate. This block may then negatively feed-back on the transcription of a subset of genes, potentially located at the nuclear envelope. Overall, this study suggests that perinuclear Mlp proteins establish a link between mRNA transcription and export.

Une méta-analyse GWA de la glycémie après 2h d'HGPO révèle que GIPR est associé avec la sécrétion de l'insuline en réponse au glucose et que l'ADCY5 est un nouveau gène de suseptibilité au diabète de type 2. [A meta-analysis of GWA blood glucose after 2 h of OGTT revealed that GIPR is associated with the secretion of insulin based on glucose and ADCY5 is a new susceptibility gene for type 2 diabetes.]

Introduction: Les études GVvA (Genome-wide association ,-studies) ont identifié et confirmé plus de 20 gènes de susceptibilité au DT2 et ont contribué à mieux comprendre sa physiopathologie. L'hyperglycémie à jeun (GJ), et 2 heures après une HGPO (G2h) sont les deux mesures cliniques du diagnostic du DT2. Nous avons identifié récemment la G6P du pancréas (G6PC2) comme déterminant de la variabilité physiologique de la GJ puis Ie récepteur à la mélatonine (MTNRIB) qui de plus lie la régulation du rythme circadien au DT2. Dans ce travail nous avons étudié la génétique de la G2h à l'aide de l'approche GWA. Résultats: Nous avons réalisé une méta-analyse GWA dans le cadre de MAGIC (Meta-Analysis of Glucose and Insulin related traits Consortium) qui a inclus 9 études GWA (N=15'234). La réplication de 29 loci (N=6958-30 121, P < 10-5 ) a confirmé 5 nouveaux loci; 2 étant connus comme associés avec Ie DT2 (TCF7L2, P = 1,6 X 10-10 ) et la GJ (GCKR, p = 5,6 X 10-10 ); alors que GIPR (p= 5,2 X 10-12), VSP13C (p= 3,9 X 10-8) et ADCY5 (p = 1,11 X 10-15 ) sont inédits. GIPR code Ie récepteur au GIP (gastric inhibitory polypeptide) qui est sécrété par les ceIlules intestinales pour stimuler la sécrétion de l'insuline en réponse au glucose (l'effet incrétine). Les porteurs du variant GIPR qui augmente la G2h ont également un indice insulinogénique plus bas, (p= 1,0 X 10-17) mais ils ne présentent aucune modification de leur glycémie suite à une hyperglycémie provoquée par voie veineuse (p= 0,21). Ces résultats soutiennent un effet incrétine du locus GIPR qui expliquerait ~9,6 % de la variance total de ce trait. La biologie de ADCY5 et VPS13C et son lien avec l'homéostasie du glucose restent à élucider. GIPR n'est pas associé avec le risque de DT2 indiquant qu'il influence la variabilité physiologique de la G2h alors que le locus ADCY5 est associé avec le DT2 (OR = 1,11, P = 1,5 X 10-15). Conclusion: Notre étude démontre que l'étude de la G2h est une approche efficace d'une part pour la compréhension de la base génétique de la physiologie de ce trait clinique important et d'autre part pour identifier de nouveaux gènes de susceptibilité au DT2.

NLRC5 exclusively transactivates MHC class I and related genes through a distinctive SXY module.

MHC class II (MHCII) genes are transactivated by the NOD-like receptor (NLR) family member CIITA, which is recruited to SXY enhancers of MHCII promoters via a DNA-binding "enhanceosome" complex. NLRC5, another NLR protein, was recently found to control transcription of MHC class I (MHCI) genes. However, detailed understanding of NLRC5's target gene specificity and mechanism of action remained lacking. We performed ChIP-sequencing experiments to gain comprehensive information on NLRC5-regulated genes. In addition to classical MHCI genes, we exclusively identified novel targets encoding non-classical MHCI molecules having important functions in immunity and tolerance. ChIP-sequencing performed with Rfx5(-/-) cells, which lack the pivotal enhanceosome factor RFX5, demonstrated its strict requirement for NLRC5 recruitment. Accordingly, Rfx5-knockout mice phenocopy Nlrc5 deficiency with respect to defective MHCI expression. Analysis of B cell lines lacking RFX5, RFXAP, or RFXANK further corroborated the importance of the enhanceosome for MHCI expression. Although recruited by common DNA-binding factors, CIITA and NLRC5 exhibit non-redundant functions, shown here using double-deficient Nlrc5(-/-)CIIta(-/-) mice. These paradoxical findings were resolved by using a "de novo" motif-discovery approach showing that the SXY consensus sequence occupied by NLRC5 in vivo diverges significantly from that occupied by CIITA. These sequence differences were sufficient to determine preferential occupation and transactivation by NLRC5 or CIITA, respectively, and the S box was found to be the essential feature conferring NLRC5 specificity. These results broaden our knowledge on the transcriptional activities of NLRC5 and CIITA, revealing their dependence on shared enhanceosome factors but their recruitment to distinct enhancer motifs in vivo. Furthermore, we demonstrated selectivity of NLRC5 for genes encoding MHCI or related proteins, rendering it an attractive target for therapeutic intervention. NLRC5 and CIITA thus emerge as paradigms for a novel class of transcriptional regulators dedicated for transactivating extremely few, phylogenetically related genes.

Positional bias of general and tissue-specific regulatory motifs in mouse gene promoters

Background: The arrangement of regulatory motifs in gene promoters, or promoterarchitecture, is the result of mutation and selection processes that have operated over manymillions of years. In mammals, tissue-specific transcriptional regulation is related to the presence ofspecific protein-interacting DNA motifs in gene promoters. However, little is known about therelative location and spacing of these motifs. To fill this gap, we have performed a systematic searchfor motifs that show significant bias at specific promoter locations in a large collection ofhousekeeping and tissue-specific genes.Results: We observe that promoters driving housekeeping gene expression are enriched inparticular motifs with strong positional bias, such as YY1, which are of little relevance in promotersdriving tissue-specific expression. We also identify a large number of motifs that show positionalbias in genes expressed in a highly tissue-specific manner. They include well-known tissue-specificmotifs, such as HNF1 and HNF4 motifs in liver, kidney and small intestine, or RFX motifs in testis,as well as many potentially novel regulatory motifs. Based on this analysis, we provide predictionsfor 559 tissue-specific motifs in mouse gene promoters.Conclusion: The study shows that motif positional bias is an important feature of mammalianproximal promoters and that it affects both general and tissue-specific motifs. Motif positionalconstraints define very distinct promoter architectures depending on breadth of expression andtype of tissue.

Gene-disease network analysis reveals functional modules in mendelian, complex and environmental diseases

AbstractBACKGROUND: Scientists have been trying to understand the molecular mechanisms of diseases to design preventive and therapeutic strategies for a long time. For some diseases, it has become evident that it is not enough to obtain a catalogue of the disease-related genes but to uncover how disruptions of molecular networks in the cell give rise to disease phenotypes. Moreover, with the unprecedented wealth of information available, even obtaining such catalogue is extremely difficult.PRINCIPAL FINDINGS: We developed a comprehensive gene-disease association database by integrating associations from several sources that cover different biomedical aspects of diseases. In particular, we focus on the current knowledge of human genetic diseases including mendelian, complex and environmental diseases. To assess the concept of modularity of human diseases, we performed a systematic study of the emergent properties of human gene-disease networks by means of network topology and functional annotation analysis. The results indicate a highly shared genetic origin of human diseases and show that for most diseases, including mendelian, complex and environmental diseases, functional modules exist. Moreover, a core set of biological pathways is found to be associated with most human diseases. We obtained similar results when studying clusters of diseases, suggesting that related diseases might arise due to dysfunction of common biological processes in the cell.CONCLUSIONS: For the first time, we include mendelian, complex and environmental diseases in an integrated gene-disease association database and show that the concept of modularity applies for all of them. We furthermore provide a functional analysis of disease-related modules providing important new biological insights, which might not be discovered when considering each of the gene-disease association repositories independently. Hence, we present a suitable framework for the study of how genetic and environmental factors, such as drugs, contribute to diseases.AVAILABILITY: The gene-disease networks used in this study and part of the analysis are available at http://ibi.imim.es/DisGeNET/DisGeNETweb.html#Download