One of the most important and challenging ``knowledge extraction' tasks in bioinformatics is the reverse engineering of gene regulatory networks (GRNs) from DNA microarray gene expression data. Indeed, as a result of the development of high-throughput data-collection techniques, biology is experiencing a data flood phenomenon that pushes biologists toward a new view of biology--systems biology--that aims at system-level understanding of biological systems.<p><p>Unfortunately, even for small model organisms such as the yeast Saccharomyces cerevisiae, the number p of genes is much larger than the number n of expression data samples. The dimensionality issue induced by this ``small n, large p' data setting renders standard statistical learning methods inadequate. Restricting the complexity of the models enables to deal with this serious impediment. Indeed, by introducing (a priori undesirable) bias in the model selection procedure, one reduces the variance of the selected model thereby increasing its accuracy.<p><p>Gaussian graphical models (GGMs) have proven to be a very powerful formalism to infer GRNs from expression data. Standard GGM selection techniques can unfortunately not be used in the ``small n, large p' data setting. One way to overcome this issue is to resort to regularization. In particular, shrinkage estimators of the covariance matrix--required to infer GGMs--have proven to be very effective. Our first contribution consists in a new shrinkage estimator that improves upon existing ones through the use of a Monte Carlo (parametric bootstrap) procedure.<p><p>Another approach to GGM selection in the ``small n, large p' data setting consists in reverse engineering limited-order partial correlation graphs (q-partial correlation graphs) to approximate GGMs. Our second contribution consists in an inference algorithm, the q-nested procedure, that builds a sequence of nested q-partial correlation graphs to take advantage of the smaller order graphs' topology to infer higher order graphs. This allows us to significantly speed up the inference of such graphs and to avoid problems related to multiple testing. Consequently, we are able to consider higher order graphs, thereby increasing the accuracy of the inferred graphs.<p><p>Another important challenge in bioinformatics is the prediction of gene function. An example of such a prediction task is the identification of genes that are targets of the nitrogen catabolite repression (NCR) selection mechanism in the yeast Saccharomyces cerevisiae. The study of model organisms such as Saccharomyces cerevisiae is indispensable for the understanding of more complex organisms. Our third contribution consists in extending the standard two-class classification approach by enriching the set of variables and comparing several feature selection techniques and classification algorithms.<p><p>Finally, our fourth contribution formulates the prediction of NCR target genes as a network inference task. We use GGM selection to infer multivariate dependencies between genes, and, starting from a set of genes known to be sensitive to NCR, we classify the remaining genes. We hence avoid problems related to the choice of a negative training set and take advantage of the robustness of GGM selection techniques in the ``small n, large p' data setting. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

La perméabilité membranaire joue un rôle important dans la résistance aux antibiotiques chez lesbactéries à Gram négatif.L’objectif de notre travail était de caractériser la fonction tenue par les deux régulateurs globaux dela perméabilité membranaire chez Enterobacer aerogenes: mar et ram. L’objectif initial futd’identifier le répresseur spécifique de RamA qui manquait en tant qu’élément de la cascade derégulation actuellement définie. La sélection de 60 souches nous a permis de confirmer le rôlecentral joué par RamA dans la régulation, ainsi qu’identifier des mutations pouvant être critiques, auniveau de RamR. Ainsi, l’absence variations observées dans le régulon marRAB et l’expressionmodérée des transcrits montrée par qRT-PCR laisse penser, que RamA a un rôle clef dans larégulation de l’expression des porines et des pompes d’efflux chez E. aerogenes.L’autre partie de notre travail reposait sur l’étude de la translocation des antibiotiques au traversdes porines. L’étude des interactions porine-carbapénèmes s’est faite sur la porine sauvage OmpFd'Escherichia coli et deux mutations. Les résultats indiquent également l'importance de l'aspartateen position 113 dans la sélectivité de translocation des carbapénèmes au sein de la porine OmpF.Ce travail montre ainsi que la translocation des pénicillines est aussi sous la dépendance desinteractions qui se créent entre le substrat et le résidu en position 113 de OmpF et limitent alorsleur passage au niveau du canal porine. Nous avons recherché la contribution attribuée à la porineOmp36 d'E. aerogenes dans la translocation de certaines béta-lactamines. Les mesures ont permisde conclure que les deux beta-lactamin / Genetic permeability plays an important role in antibiotic resistance of Gram-negative bacteria.Our work was to characterize and better understand of the genetic regulation of membranepermeability in E. aerogenes. We focused on two global regulators, mar and ram, in about 60clinical isolates. Alterations in the upstream region of ramA and in ramR but no mutations in marAnor marR were observed. Overexpression of ramA or ramR led to an altered antibiotic susceptibilityassociated to decrease of porins expression and over-expression of efflux-pumps. qRT-PCR pointedout the estimated importance of the ram-regulon in the regulation cascade.Another part of this work was to characterize the translocation of compounds through porins andthe role of porins in drug uptake in general. Measurement of the rate of antibiotic action of threecarbapenems in an E. coli strain solely expressing OmpF as porin clearly indicated the importanceof the aspartate at position 113 in antibiotic translocation. A multi-disciplinary three way approachof computer modeling, black-lipid-bilayer assays and measurement of antibiotic action, suggestedthat interactions with residue D113 of E. coli porin OmpF are rate-limiting for transport throughthe porin channel. Combination of biological and biophysical measurements with E. aerogenesporin Omp36 denoted that interactions between the porin channel and the antibiotic facilitate andaccelerate transport.

Skeletal muscle wasting is a consequence of numerous physiological conditions, including denervation, corticosteroid treatment, immobilization, and aging. The E3 ubiquitin ligases, MuRF1 and MAFbx, are induced under nearly all atrophy conditions and are believed to play a key role in protein degradation in atrophying muscle. However, the preliminary data described in this study provides new evidence that MuRF1 may also act as a transcriptional modulator of atrophy-induced gene activity, including the regulation of MAFbx and MuRF1 expression. To characterize the transcriptional regulation of MuRF1 and MAFbx, reporter gene constructs containing fragments of the proximal promoter regions of these genes were developed, transfected into C2C12 cells with or without a MuRF1 expression plasmid and monitored for differences in reporter gene activity. The MuRF1 and MAFbx reporters each showed repressed activity in cells ectopically expressing MuRF1 compared to cells that did not overexpress MuRF1. Furthermore, ectopic expression of the myogenic regulatory factors (MRFs), MyoD1 and myogenin, caused significant activation of the MuRF1 and MAFbx reporter constructs. However, co-overexpression of MuRF1 with MyoD1 or myogenin resulted in reversal of MRF induction of reporter gene activity, and synergistic repression of a constructed E-box reporter system. To further characterize the role of the MuRF1 gene product in repression of MuRF1 expression, a MuRF1 RING domain mutant and a MuRF1 c-terminal mutant were created. The mutant constructs were then co-transfected along with MRF expression plasmids and the MuRF1 reporter construct into C2C12 cells and reporter gene activity was assessed. The MuRF1 RING mutant failed to reverse MRF activation of the reporter gene, while the c-terminal mutant successfully reversed activation of the reporter gene. These findings suggest that ubiquitin ligase activity is required for MuRF1 transcriptional regulatory effects. These data offer exciting evidence of a potential new function for MuRF1 as a transcriptional modulator of atrophy-induced changes in gene expression.

Indiana University-Purdue University Indianapolis (IUPUI) / The LHX3 transcription factor is essential for pituitary gland and nervous system development in mammals. In humans, mutations in the LHX3 gene underlie combined pituitary hormone deficiency (CPHD) disease featuring deficits in anterior pituitary hormones and defects in the nervous system. The mechanisms that control temporal and spatial expression of the LHX3 gene are poorly understood. The proximal promoters of the human LHX3 gene are insufficient to guide expression in vivo and downstream elements including a conserved 7.9 kilobase (kb) enhancer region appear to play a role in tissue-specific expression in the pituitary and nervous system. In this study, I characterized the activity of this downstream enhancer region in regulating gene expression at the cellular level during development. Human LHX3 enhancer-driven Cre reporter transgenic mice were generated to facilitate studies of enhancer actions. The downstream LHX3 enhancer primarily guides gene transcription in αGSU-expressing cells secreting the TSHβ, LHβ or FSHβ hormones and expressing the GATA2 and SF1 transcription factors. In the developing nervous system, the enhancer serves as a targeting module for expression specifically in V2a interneurons. These results demonstrate that the downstream LHX3 enhancer is important in specific endocrine and neural cell types but also indicate that additional regulatory elements are likely involved in LHX3 gene expression in other cell types. Further, these studies demonstrate significant gonadotrope cell heterogeneity during pituitary development, providing insights into the cellular physiology of this key reproductive regulatory cell. The human LHX3 enhancer-driven Cre reporter transgenic mice provide a valuable tool for further developmental studies of cell determination and differentiation in the pituitary and nervous system. Furthermore understanding the regulation of human LHX3 gene will help develop tools to better diagnose and treat pituitary CPHD disease.

Indiana University-Purdue University Indianapolis (IUPUI) / LIM homeodomain 3 (LHX3) is an important regulator of pituitary and nervous system development. To date, twelve LHX3 gene mutations have been identified in patients with combined pituitary hormone deficiency disease (CPHD). Understanding the molecular mechanisms governing LHX3/Lhx3 gene regulation will provide critical insights into organ development pathways and associated diseases. DNA methylation has been implicated in gene regulation in multiple physiological systems. This dissertation examines the role of DNA methylation in regulating the murine Lhx3 gene. To determine if demethylation of the Lhx3 gene promoter would induce its expression, murine pre-somatotrope pituitary cells that do not normally express Lhx3 (Pit-1/0 cells) were treated with the demethylating reagent, 5-Aza-2’-deoxycytidine. This treatment lead to activation of the Lhx3 gene and thus suggested that methylation contributes to Lhx3 gene regulation. Proteins that modify chromatin, such as histone deacetylases (HDACs) have also been shown to affect DNA methylation patterns and subsequent gene activation. Pit-1/0 pituitary cells treated with a combination of the demethylating reagent and the HDAC inhibitor, Trichostatin A led to activation of the Lhx3 gene, suggesting crosstalk between DNA methylation and histone modification processes. To assess DNA methylation levels, treated and untreated Pit-1/0 genomic DNA were subjected to bisulfite conversion and sequencing. Treated Pit-1/0 cells had decreased methylation compared to untreated cells. Chromatin immunoprecipitation assays demonstrated interactions between the methyl-binding protein, MeCP2 and the Lhx3 promoter regions in the Pit-1/0 cell line. Overall, the study demonstrates that DNA methylation patterns of the Lhx3 gene are associated with its expression status.

Indiana University-Purdue University Indianapolis (IUPUI) / Histone acetylation has been linked to developmental changes in gene expression and is a validated drug target of apicomplexan parasites, but little is known about the roles of individual histone modifying enzymes and how they are recruited to target genes. The protozoan parasite Toxoplasma gondii (phylum Apicomplexa) is unusual among invertebrates in possessing two GCN5-family lysine acetyltransferases (KATs). While GCN5a is required for gene expression in response to alkaline stress, this KAT is dispensable for parasite proliferation in normal culture conditions. In contrast, GCN5b cannot be disrupted, suggesting it is essential for Toxoplasma viability. To further explore the function of GCN5b, we generated clonal parasites expressing an inducible HA-tagged form of GCN5b containing a point mutation that ablates enzymatic activity (E703G). Stabilization of this dominant-negative form of GCN5b was mediated through ligand-binding to a destabilization domain (dd) fused to the protein. Induced accumulation of the ddHAGCN5b(E703G) protein led to a rapid arrest in parasite replication. Growth arrest was accompanied by a decrease in histone H3 acetylation at specific lysine residues as well as reduced expression of GCN5b target genes in GCN5b(E703G) parasites, which were identified using chromatin immunoprecipitation coupled with microarray hybridization (ChIP-chip). We also demonstrate that GCN5b interacts with AP2-domain proteins, which are plant-like transcription factors in Apicomplexa. The interactions between GCN5b, AP2IX-7, and AP2X-8 were confirmed by reciprocal co-immunoprecipitation and revealed a “core complex” that includes the co-activator ADA2-A, TFIID subunits, LEO1 polymerase-associated factor (Paf1) subunit, and RRM proteins. The dominant-negative phenotype of ddHAGCN5b(E703G) parasites, considered with the proteomics and ChIP-chip data, indicate that GCN5b plays a central role in transcriptional and chromatin remodeling complexes. We conclude that GCN5b has a non-redundant and indispensable role in regulating gene expression required during the Toxoplasma lytic cycle.

Indiana University-Purdue University Indianapolis (IUPUI) / In systems biology study, biological networks were used to gain insights into biological systems. While the traditional approach to studying biological networks is based on the identification of interactions among genes or the identification of a gene set ranking according to differentially expressed gene lists, little is known about interactions between higher order biological systems, a network of gene sets. Several types of gene set network have been proposed including co-membership, linkage, and co-enrichment human gene set networks. However, to our knowledge, none of them contains directionality information. Therefore, in this study we proposed a method to construct a regulatory gene set network, a directed network, which reveals novel relationships among gene sets. A regulatory gene set network was constructed by using publicly available gene regulation data. A directed edge in regulatory gene set networks represents a regulatory relationship from one gene set to the other gene set. A regulatory gene set network was compared with another type of gene set network to show that the regulatory network provides additional information. In order to show that a regulatory gene set network is useful for understand the underlying mechanism of a disease, an Alzheimer's disease (AD) regulatory gene set network was constructed. In addition, we developed Pathway and Annotated Gene-set Electronic Repository (PAGER), an online systems biology tool for constructing and visualizing gene and gene set networks from multiple gene set collections. PAGER is available at http://discern.uits.iu.edu:8340/PAGER/. Global regulatory and global co-membership gene set networks were pre-computed. PAGER contains 166,489 gene sets, 92,108,741 co-membership edges, 697,221,810 regulatory edges, 44,188 genes, 651,586 unique gene regulations, and 650,160 unique gene interactions. PAGER provided several unique features including constructing regulatory gene set networks, generating expanded gene set networks, and constructing gene networks within a gene set. However, tissue specific or disease specific information was not considered in the disease specific network constructing process, so it might not have high accuracy of presenting the high level relationship among gene sets in the disease context. Therefore, our framework can be improved by collecting higher resolution data, such as tissue specific and disease specific gene regulations and gene sets. In addition, experimental gene expression data can be applied to add more information to the gene set network. For the current version of PAGER, the size of gene and gene set networks are limited to 100 nodes due to browser memory constraint. Our future plans is integrating internal gene or proteins interactions inside pathways in order to support future systems biology study.

Indiana University-Purdue University Indianapolis (IUPUI) / RNA-binding proteins (RBPs) are vital post-transcriptional regulatory molecules in transcriptome of mammalian species. It necessitates studying their expression dynamics to extract how post-transcriptional networks work in various mammalian tissues. RNA binding proteins (RBPs) play important roles in controlling the post-transcriptional fate of RNA molecules, yet their evolutionary dynamics remains largely unknown. As expression profiles of genes encoding for RBPs can yield insights about their evolutionary trajectories on the post-transcriptional regulatory networks across species, we performed a comparative analyses of RBP expression profiles across 8 tissues (brain, cerebellum, heart, lung, liver, lung, skeletal muscle, testis) in 11 mammals (human, chimpanzee, gorilla, orangutan, macaque, rat, mouse, platypus, opossum, cow) and chicken & frog (evolutionary outgroups). Noticeably, orthologous gene expression profiles suggest a significantly higher expression level for RBPs than their non-RBP gene counterparts, which include other protein-coding and non-coding genes, across all the mammalian tissues studied here. This trend is significant irrespective of the tissue and species being compared, though RBP gene expression distribution patterns were found to be generally diverse in nature. Our analysis also shows that RBPs are expressed at a significantly lower level in human and mouse tissues compared to their expression levels in equivalent tissues in other mammals: chimpanzee, orangutan, rat, etc., which are all likely exposed to diverse natural habitats and ecological settings compared to more stable ecological environment humans and mice might have been exposed, thus reducing the need for complex and extensive post-transcriptional control. Further analysis of the similarity of orthologous RBP expression profiles between all pairs of tissue-mammal combinations clearly showed the grouping of RBP expression profiles across tissues in a given mammal, in contrast to the clustering of expression profiles for non-RBPs, which frequently grouped equivalent tissues across diverse mammalian species together, suggesting a significant evolution of RBPs expression after speciation events. Calculation of species specificity indices (SSIs) for RBPs across various tissues, to identify those that exhibited restricted expression to few mammals, revealed that about 30% of the RBPs are species-specific in at least one tissue studied here, with lung, liver, kidney & testis exhibiting a significantly higher proportion of species specifically expressed RBPs. We conducted a differential expression analysis of RBPs in human, mouse and chicken tissues to study the evolution of expression levels in recently evolved species (i.e., humans and mice) than evolutionarily-distant species (i.e., chickens). We identified more than 50% of the orthologous RBPs to be differentially expressed in at least one tissue, compared between human and mouse, but not so between human and an outgroup chicken, in which RBP expression levels are relatively conserved. Among the studied tissues (brain, liver and kidney) showed a higher fraction of differentially expressed RBPs, which may suggest hyper- regulatory activities by RBPs in these tissues with species evolution. Overall, this study forms a foundation for understanding the evolution of expression levels of RBPs in mammals, facilitating a snapshot of the wiring patterns of post-transcriptional regulatory networks in mammalian genomes. In our second study, we focused on elucidating novel features of post-transcriptional regulatory molecules called as circRNA from LongPolyA RNA-sequence data. The debate over presence of nonlinear exon splicing such as exon-shuffling or formation of circularized forms has finally come to an end as numerous repertoires have shown of their occurrence and presence through transcriptomic analyses. It is evident from previous studies that along with consensus-site splicing non-consensus site splicing is robustly occurring in the cell. Also, in spite of applying different high-throughput approaches (both computational and experimental) to determine their abundance, the signal is consistent and strongly conforming the plausible circularization mechanisms. Earlier studies hypothesized and hence focused on the ribo-minus non-polyA RNA-sequence data to identify circular RNA structures in cell and compared their abundance levels with their linear counterparts. Thus far, the studies show their conserved nature across tissues and species also that they are not translated and preferentially are without poly (A) tail, with one to five exons long. Much of this initial work has been performed using non-polyA sequencing thus probably underestimates the abundance of circular RNAs originating from long poly (A) RNA isoforms. Our hypothesis is if the circular RNA events are not the artifact of random events, but has a structured and defined mechanism for their formation, then there would not be biases on preferential selection / leaving of polyA tails, while forming the circularized isoforms. We have applied an existing computational pipeline from earlier studies by Memczack et. al., on ENCODE cell-lines long poly (A) RNA-sequence data. With the same pipeline, we achieve a significant number of circular RNA isoforms in the data, some of which are overlapping with known circular RNA isoforms from the literature. We identified an approach and worked upon to identify the precise structure of circular RNA, which is not plausible from the existing computational approaches. We aim to study their expression profiles in normal and cancer cell-lines, and see if there exists any pattern and functional significance based on their abundance levels in the cell.

Indiana University-Purdue University Indianapolis (IUPUI) / In response to different environmental stresses, phosphorylation of eukaryotic initiation factor-2 (eIF2) rapidly reduces protein synthesis, which lowers energy expenditure and facilitates reprogramming of gene expression to remediate stress damage. Central to the changes in gene expression, eIF2 phosphorylation also enhances translation of ATF4, a transcriptional activator of genes subject to the Integrated Stress Response (ISR). The ISR increases the expression of genes important for alleviating stress, or alternatively triggering apoptosis. One ISR target gene encodes the transcriptional regulator CHOP whose accumulation is critical for stress-induced apoptosis. In this dissertation research, I show that eIF2 phosphorylation induces preferential translation of CHOP by a mechanism involving a single upstream ORF (uORF) located in the 5’-leader of the CHOP mRNA. In the absence of stress and low eIF2 phosphorylation, translation of the uORF serves as a barrier that prevents translation of the downstream CHOP coding region. Enhanced eIF2 phosphorylation during stress facilitates ribosome bypass of the uORF, and instead results in the translation of CHOP. Stable cell lines were also constructed that express CHOP transcript containing the wild type uORF or deleted for the uORF and each were analyzed for expression changes in response to the different stress conditions. Increased CHOP levels due to the absence of inhibitory uORF sensitized the cells to stress-induced apoptosis when compared to the cells that express CHOP mRNA containing the wild type uORF. This new mechanism of translational control explains how expression of CHOP and the fate of cells are tightly linked to the levels of phosphorylated eIF2 and stress damage.

Indiana University-Purdue University Indianapolis (IUPUI) / The epigenetic state of any cell is, in part, regulated by the interaction of DNA with nuclear histones. Histone tails can be modified in a number of ways that impact on the availability of DNA to interact with transcriptional complexes, including methylation, acetylation, phosphorylation, ubiquituination, and sumoylation. Histones are acetylated by a large family of enzymes, histone acetyl transferases (HATs), and deacetylated by the histone deacetylases (HDACs). Acetylated histones are generally considered markers of genomic regions that are actively being transcribed, whereas deacetylated and methylated histones are generally markers of regions that are inactive. The goal of the present study was to 1) study the epigenetic state with regard to the presence of euchromatin and heterochromatin in the developing chick and murine retina, 2) study and compare the localization patterns of the classical HDACs in the developing chick and murine retina with respect retinal progenitors and early differentiated cell types 3) to test the hypothesis that overall HDAC activity is required for dividing retinal progenitors to leave the cell cycle and differentiate. Our results showed that the classical HDACs were ubiquitously expressed in the developing chick and murine retinas. Species specific differences as well as stage dependent variations were observed in the localization of the HDACs in the cell types that were studied in the chick and murine retina. Our preliminary results also showed that HDAC inhibition may lead to the inability of the cell types to leave the cell cycle and a subsequent increase in the number of progenitor cells present in the developing chick retina.