Understanding transcriptional regulation quantitatively is a crucial step towards uncovering and ultimately utilizing the regulatory semantics encoded in the genome. Transcription of a gene is induced by the binding of site-specific transcription factors (TFs) to so-called cis-regulatory-modules (CRMs). The frequency and duration of the binding events are influenced by the concentrations of the TFs, the binding affinities and location of the transcription factor binding sites (TFBSs) in the CRM as well as the properties of the TFs themselves (e.g. effectiveness, competitive interaction with other TFs). Modeling these interactions using a mathematical approach, based on sound biochemical and thermodynamic foundations, enables the understanding and quantitative prediction of transcriptional output of a target gene. In the thesis I introduce the developed framework for modeling, visualizing, and predicting the regulation of the transcription rate of a target gene. Given the concentrations of a set of TFs, the TFBSs in a regulatory DNA region, and the transcription rate of the target gene, the method will optimize its parameters to generate a predictive model that incorporates the regulatory mechanism of the observed gene. I demonstrate the generalization ability of the model by subjecting it to standard machine learning and hypothesis testing procedures. Furthermore, I demonstrate the potential of the approach by training the method on a gene in D. melanogaster and predicting the output of the homologous genes in 12 other Drosophila species where the regulatory sequence has evolved substantially but the regulatory mechanism was conserved. Finally, I investigate the proposed role-switching behaviour of TFs regulating the development of D. melanogaster, which suggests that SUMOylation is the biological mechanism facilitating the switch.

Understanding transcriptional regulation quantitatively is a crucial step towards uncovering and ultimately utilizing the regulatory semantics encoded in the genome. Transcription of a gene is induced by the binding of site-specific transcription factors (TFs) to so-called cis-regulatory-modules (CRMs). The frequency and duration of the binding events are influenced by the concentrations of the TFs, the binding affinities and location of the transcription factor binding sites (TFBSs) in the CRM as well as the properties of the TFs themselves (e.g. effectiveness, competitive interaction with other TFs). Modeling these interactions using a mathematical approach, based on sound biochemical and thermodynamic foundations, enables the understanding and quantitative prediction of transcriptional output of a target gene. In the thesis I introduce the developed framework for modeling, visualizing, and predicting the regulation of the transcription rate of a target gene. Given the concentrations of a set of TFs, the TFBSs in a regulatory DNA region, and the transcription rate of the target gene, the method will optimize its parameters to generate a predictive model that incorporates the regulatory mechanism of the observed gene. I demonstrate the generalization ability of the model by subjecting it to standard machine learning and hypothesis testing procedures. Furthermore, I demonstrate the potential of the approach by training the method on a gene in D. melanogaster and predicting the output of the homologous genes in 12 other Drosophila species where the regulatory sequence has evolved substantially but the regulatory mechanism was conserved. Finally, I investigate the proposed role-switching behaviour of TFs regulating the development of D. melanogaster, which suggests that SUMOylation is the biological mechanism facilitating the switch.

Understanding transcriptional regulation quantitatively is a crucial step towards uncovering and ultimately utilizing the regulatory semantics encoded in the genome. Transcription of a gene is induced by the binding of site-specific transcription factors (TFs) to so-called cis-regulatory-modules (CRMs). The frequency and duration of the binding events are influenced by the concentrations of the TFs, the binding affinities and location of the transcription factor binding sites (TFBSs) in the CRM as well as the properties of the TFs themselves (e.g. effectiveness, competitive interaction with other TFs). Modeling these interactions using a mathematical approach, based on sound biochemical and thermodynamic foundations, enables the understanding and quantitative prediction of transcriptional output of a target gene. In the thesis I introduce the developed framework for modeling, visualizing, and predicting the regulation of the transcription rate of a target gene. Given the concentrations of a set of TFs, the TFBSs in a regulatory DNA region, and the transcription rate of the target gene, the method will optimize its parameters to generate a predictive model that incorporates the regulatory mechanism of the observed gene. I demonstrate the generalization ability of the model by subjecting it to standard machine learning and hypothesis testing procedures. Furthermore, I demonstrate the potential of the approach by training the method on a gene in D. melanogaster and predicting the output of the homologous genes in 12 other Drosophila species where the regulatory sequence has evolved substantially but the regulatory mechanism was conserved. Finally, I investigate the proposed role-switching behaviour of TFs regulating the development of D. melanogaster, which suggests that SUMOylation is the biological mechanism facilitating the switch.

[Truncated abstract] Human complement receptor 2 (CR2/CD21) is a cell surface glycoprotein detected on specific cells involved in immunity, which binds complement C3 cleavage fragments, cellular ligands IFN-? and CD23 as well as the EBV coat protein, gp350/220. During the early stages of B-cell development CR2/CD21 is silenced. Expression is initiated on immature B-cells escaping negative selection. During peripheral maturation CR2/CD21 is up-regulated with B-cell sub-populations showing distinctive surface levels (comparatively low, intermediate or high). CR2/CD21 is silenced upon terminal plasmacytic differentiation. Appropriate timing and expression level of CR2/CD21 is important for the development of a healthy B-cell repertoire. Previous studies have identified sequences within the proximal promoter and first intron of CR2/CD21 that cooperate within native chromatin to control cell-specific silencing. Further, analysis of cultured human cells has revealed chromatin structural variation causing DNase I hypersensitivity at these regulatory sites in a CR2/CD21-expressing mature B-cell line (Raji) which are absent in a non-lymphoid cell type (K562). The primary focus of the present study involved characterising chromatin structural variation over previously recognized DNase I hypersensitive regions at the CR2/CD21 locus in human cells to understand how chromatin structure might regulate developmental expression of CR2/CD21. ... These studies provide evidence that notch signaling influences CR2/CD21 expression in human cell lines. First, in vivo binding of CBF1 to CR2/CD21 sequences in the proximal promoter and CRS implies that CR2/CD21 is a direct target of notch activation. Second, the effect of exogenous notch signalling molecules on CR2/CD21 proximal promoter activity was modulated by factors binding tandem E-boxes near the transcriptional start site suggesting that the notch pathway may also influence CR2/CD21 expression via control of HLH molecules. Third, initiation of CR2/CD21 expression was observed in a nonexpressing pre-B cell line (Reh) by co-culture with stromal cells expressing a notch ligand (OP9-DL) but not control stroma (OP9-GFP). Together, these findings support a role for notch regulation of B-cell maturation and invite speculation that initiation of CR2/CD21 expression following negative selection of immature B-cells involves crosstalk between HLH transcriptional regulators and the notch pathway. Furthermore, the Reh/OP9-DL co-culture system may provide a model to directly study the relationship between cell signalling molecules, transcription factor regulation, chromatin structural variation and differentiation of B-cells.

Historically, sperm has been considered merely as a carrier of genetic material at fertilisation. However, it is known that sperm supports embryonic development better than other cell types, suggesting that it might also have additional important, non-genetic contributions to embryonic development. The work described in this dissertation focuses on identifying the molecular determinants of developmental programming of sperm. First, the development of embryos derived from sperm and spermatids, immature precursors of sperm was compared. Sperm-derived embryos developed significantly better than spermatid-derived embryos. Further research aiming to identify the reasons for the developmental advantage of sperm led to the identification of proteins that are present specifically in sperm and not in spermatids. Moreover, egg factors which are preferentially incorporated into the sperm, but not into the spermatid chromatin were identified with the use of egg extracts, suggesting that the chromatin of sperm could be programmed to interact with the components of the egg. Subsequently, the reasons for developmental failure of spermatid-derived embryos were investigated. By comparing the sperm with spermatids it was shown that the programming of sperm to support efficient development is linked to its special ability to regulate expression of developmentally-important embryonic genes, and not to its ability to support DNA replication or rRNA production. Further characterisation of the sperm and spermatid chromatin with the use of genome-wide sequencing allowed me to link the correct regulation of gene expression in the embryo with a certain combination of epigenetic marks in the sperm, but not in the spermatid chromatin. Finally, it is shown that enzymatic removal of epigenetic modifications at fertilisation leads to misregulation of gene expression. This therefore suggests that epigenetic information contained in parental genomes at fertilisation is required for a proper regulation of embryonic transcription. My results support the hypothesis that the sperm is not only a carrier of genetic material, but also provides the embryo with epigenetic information for regulation of transcription after fertilisation. I believe that these findings advance our current understanding of the nature and mechanisms of sperm programming for embryonic development, and are important contributions to the emerging field of transgenerational inheritance of epigenetic traits in general.

Submitted by Karin Goebel (karing@fiocruz.br) on 2014-11-25T18:05:56Z No. of bitstreams: 1 Dissertação Bruna Hilzendeger Marcon.pdf: 6455497 bytes, checksum: 8ea632ce91cdf16edd8b86a624972dba (MD5) / Approved for entry into archive by Karin Goebel (karing@fiocruz.br) on 2014-11-25T18:06:29Z (GMT) No. of bitstreams: 1 Dissertação Bruna Hilzendeger Marcon.pdf: 6455497 bytes, checksum: 8ea632ce91cdf16edd8b86a624972dba (MD5) / Made available in DSpace on 2014-11-25T18:06:29Z (GMT). No. of bitstreams: 1 Dissertação Bruna Hilzendeger Marcon.pdf: 6455497 bytes, checksum: 8ea632ce91cdf16edd8b86a624972dba (MD5) Previous issue date: 2014 / Fundação Oswaldo Cruz. Instituto Carlos Chagas. Curitiba, PR, Brasil / As células-tronco (CTs) caracterizam-se por possuírem a capacidade de se autorrenovar e de dar origem a um ou mais tipos celulares diferenciados. Nos últimos anos, diversos trabalhos mostraram a existência de CTs em tecidos adultos, tornando-as uma alternativa interessante para uso em terapias celulares. Contudo, para melhor utilizar as CTs, é preciso primeiramente compreender como ocorre a diferenciação em um tipo celular específico e, principalmente, como é regulada a expressão gênica durante este processo. Em 2009, Ingolia e colaboradores apresentaram uma nova técnica conhecida como ribosome profiling, a qual consiste no isolamento e sequenciamento em larga escala dos fragmentos de RNA associados e protegidos pelos ribossomos, os quais têm um tamanho aproximado de 30 nucleotídeos (conhecido com footprint ribossomal). Ao mapear as sequências obtidas, é possível obter informações não apenas sobre quais sequências estão sendo traduzidas, mas também sobre a cinética da tradução e sua extensiva rede de regulação. Assim, o objetivo deste trabalho foi aplicar a técnica de ribosome profiling ao estudo do processo de diferenciação de CTs adultas. Como modelo de estudo, foram utilizadas CTs obtidas de tecido adiposo antes (t=0) e após a indução para diferenciação adipogênica por 3 dias (t=72h). O primeiro passo do trabalho foi a adaptação do protocolo de ribosome profiling para o estudo de CTs adultas, o qual consiste na lise celular, digestão do lisado com uma RNA nuclease (a qual irá degradar o RNA exposto, preservando os fragmentos protegidos pelo ribossomo), ultracentifrugação do homogenato sobre colchão de sacarose 1 M para sedimentação dos ribossomos, extração de RNA e isolamento dos fragmentos de 30 nucleotídeos. Também foi feita extração do RNA poliA. As amostras foram sequenciadas (SOLiD™) e os dados obtidos foram triados e mapeados contra um banco de dados de RNAm, utilizando-se a ferramenta CLC Genomics Workbench. Foram identificados mais de 8.000 transcritos para as amostras de ribosome profiling e mais de 17.000 para as de poliA. Ao calcular o fold change entre as condições t=0 e t=72h, foi possível verificar que mais de 50% dos genes foram detectados como diferencialmente expressos apenas por ribosome profiling. Observou-se que genes relacionados com vias de diferenciação adipogênica e de metabolismo de lipídeos encontravam-se regulados positivamente em ambas as amostras de RNA. Por outro lado, observou-se que vias de regulação do citoesqueleto de actina e de adesão focal estavam reguladas negativamente apenas nas amostras de ribosome profiling. Isso é interessante, uma vez que a inibição destas vias já foi descrita como importante para o processo de adipogênese. Além disso, foi observada uma forte redução na eficiência de tradução de genes relacionados com a tradução após 72 horas de indução para diferenciação. Os resultados obtidos no presente trabalho reforçam as evidências de que os mecanismos de regulação pós-transcricionais e traducionais têm um papel muito importante na regulação da diferenciação celular de CTs, sendo que a técnica de ribosome profiling permitiu obter informações mais detalhadas de como este processo pode estar acontecendo. / Stem cells (SC) are characterized by their capacity of both self-renewing and giving rise to new differentiated cells. SC are found in adult tissues, which are considered a putative source for cell therapy. However, little is known about the mechanisms involved in the trigger of SCs differentiation into a specific cell type. Understanding adult SCs differentiation process is a fundamental step to better use and to take advantage of their potential. In 2009, Ingolia and collaborators presented a new methodology of transcriptome analysis named ribosome profiling, which consists on the isolation and deep-sequencing of the mRNA fragments enclosed by ribosomes. When lysed cells are submitted to nuclease digestion, unprotected mRNA is degraded, while fragments within ribosomes are preserved and have a known footprint of 30 nucleotides. Sequencing these ribosome-protected fragments results in a high-precision measurement of in vivo translation, providing precise information about translation kinetics and its extensive regulation. The objective of this work was to apply the ribosome profiling methodology to the study of adipogenic differentiation in adult SCs. SCs were isolated from human adipose tissue from three donors and were cultured in a control medium (t=0) and induced to adipogenic differentiation for 72 hours (t=72h). The first step was to adapt and optimize the ribosome profiling protocol to the SC model, which consists in cell lysis, cell lysate digestion by nuclease (to degrade unprotected RNA, preserving ribosome-protected fragments), ultracentrifugation over a 1M sucrose cushion to pellet ribosomes, RNA extraction and 30 nucleotides fragments isolation. poliA RNA was also isolated. Samples were submitted to deep-sequencing (SOLiD™) and the reads obtained were trimmed and mapped onto the reference mRNA database using the CLC Genomics Workbench. Over 8000 transcripts were identified in ribosome profiling samples and over 17000 in poliA samples. Fold change analysis between t=0 and t=72h of both RNA samples showed that differential expression of more than 50% of the genes was identified only by ribosome profiling. Pathways related to adipogenesis and lipid metabolism were upregulated in both RNA samples. However, regulation of the actin cytoskeleton and focal adhesion proteins were downregulated only in ribosome profiling samples. Interestingly, downregulation of these pathways was already described as an important phenomenon to cell adipogenesis. Besides, we observed a strong reduction of translational efficiency of genes involved in translation at t=72h. Our results reinforce previous data, suggesting that posttranscriptional and translational regulation play a fundamental role in the regulation of SC differentiation process and that ribosome profiling is an important tool to better understand this process.

L’inhibition, par p16INK4a, de la progression du cycle cellulaire est considérée comme liée à un arrêt de la progression en phase G1 du à l’inhibition de l’activité de CDK4/6. Nous montrons que l’expression ectopique de p16INK4a dans trois lignées cellulaires malignes, p16-/- et pRb+/+, issues de tissus différents, provoque un allongement de la durée de la phase S et du cycle cellulaire total. L’ensemble de nos travaux sur p16INK4a sauvage et son mutant p16G101W indique que p16INK4a induit un allongement de la phase S i) indépendamment de l’origine tissulaire des cellules analysées et ii) en partie lié aux conséquences de l’inhibition de l’activité de CDK4/6 et peut-être des MAP-kinases. Dans sa localisation nucléaire, p16INK4a interviendrait dans la régulation du cycle cellulaire indépendamment de sa liaison à CDK4. L’expression de CDK1 est inhibée par p16INK4a dans les trois lignées analysées. Dans les cellules MCF7 et U87, cette inhibition est post-transcriptionnelle, médiée par la région 3’non traduite de l’ARNm de CDK1, et est associée à une modification de l’équilibre d’expression, tissu-spécifique, des microARN régulant potentiellement CDK1. Nous démontrons que CDK1 est une cible de miR- 410 et miR-650 induits par p16INK4a et le rôle de l’inhibition de la voie pRb/E2F par p16INK4a dans l’induction de miR-410. Ainsi, p16INK4a régule l’expression des gènes à différents niveaux en modifiant l’équilibre fonctionnel des facteurs de transcription et, en conséquence, des miARN. / The inhibition of cell cycle progression by p16INK4a, have been considered to result from arrest in G1 phase due to inhibition of CDK4/6 activity. We show that ectopic expression of p16INK4a in three human malignant cell lines, p16-/- and pRb +/ +, derived from different tissues, led to an increase in the length of S phase and of the entire cell cycle. Our studies using wild-type p16INK4a and p16G101W mutant indicated that p16INK4a induces a lengthening of S phase i) independently of tissue origins and ii) partly linked to the inhibition of CDK4/6 activity and possibly MAP-kinases. In the nucleus, p16INK4a may intervene in regulating the cell cycle independently of binding to CDK4. The expression of CDK1 is inhibited by p16INK4a in three cell lines analyzed. In MCF7 and U87cells, this inhibition is post-transcriptional, mediated by the 3' non translated region of CDK1 mRNA, and is associated with changes in the balance of the expression of microRNAs, which regulate potentially CDK1. We demonstrate that CDK1 is a target of miR-410 and miR-650, both induced by p16INK4a and the role of the inhibition of pRb/E2F pathway by p16INK4a in the induction of miR-410. Thus, p16INK4a regulates gene expression at different levels by modifying the functional balance of transcription factors and, consequently, the microRNAs

Nesta tese propomos um modelo estocástico multimodal para regulação da expressão gênica em nível de transcrição. A definição de um espaço de parâmetros que contém o conteúdo biológico do sistema aliada à escolha apropriada de uma base para construir a matriz de acoplamento entre os estados do sistema levaram à obtenção de soluções exatas do modelo. Tais soluções são obtidas transformando as equações mestras em equações diferenciais parciais usando a técnica das funções geradoras e escrevendo os coeficientes das equações parciais em termos dos parâmetros biológicos do modelo. No regime estacionário obtivemos uma relação de recorrência para os coeficientes das séries de potências que definem as funções geradoras e a especificação das configurações de equilíbrio do sistema permite que estas séries sejam calculadas exatamente. Com as soluções exatas calculadas não só as distribuições de probabilidade foram obtidas como os momentos das distribuições. As distribuições de probabilidade de equilíbrio apresentam estruturas multimodais com vários picos e a análise do ruído (flutuação) mostra que a existência de um estado intermediário de eficiência transcricional leva a redução do ruído global do sistema. A inspeção dos autovalores da matriz de acoplamento mostrou que existem regiões onde a dinâmica dos momentos é de caráter oscilatória com amortecimento. Diferentes esquemas de acoplamento levam à diferentes regimes transientes, tal característica revela que o sistema multimodal apresentam maior flexibilidade adaptativa quando comparado com sistemas de um ou dois estados. / In this thesis we propose a stochastic model for multimodal regulation of gene expression at the transcriptional level. The definition of a parameter space that contains the contents of the biological system coupled with the appropriate choice of a base to build the coupling matrix between the states of the system led to the exact solutions of the model. Such solutions are obtained by transforming master equations in partial differential equations using the technique of generating functions and writing the coefficients of partial equations in terms of biological parameters of the model. In the steady a recurrence relation for the coefficients of power series defining the generating functions was obtained and specification of the equilibrium configurations of the system allows the exact calculation of these series. With the exact solutions calculated not only the probability distributions were obtained but also the moments of the distributions. The equilibrium distributions probability is multimodal and presents several peaks. Analysis of the noise (fluctuation) shows that the existence of an state with intermediate transcriptional efficiency leads to a reduction of the overall system noise. Inspection of the eigenvalues of the coupling matrix showed that there are regions where the dynamics of the moments is damped oscillating. Different coupling schemes lead to different transient regimes, this feature reveals that the multimodal system have greater adaptive flexibility when compared to systems of one or two states.

Submitted by Cláudia Bueno (claudiamoura18@gmail.com) on 2015-11-12T18:06:58Z No. of bitstreams: 2 Dissertação - Iane de Oliveira Pires Porto - 2014.pdf: 1352493 bytes, checksum: 3e4c1213c96035cbae32d6eeccdd14e5 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2015-11-13T10:33:50Z (GMT) No. of bitstreams: 2 Dissertação - Iane de Oliveira Pires Porto - 2014.pdf: 1352493 bytes, checksum: 3e4c1213c96035cbae32d6eeccdd14e5 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2015-11-13T10:33:50Z (GMT). No. of bitstreams: 2 Dissertação - Iane de Oliveira Pires Porto - 2014.pdf: 1352493 bytes, checksum: 3e4c1213c96035cbae32d6eeccdd14e5 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2014-02-19 / Conselho Nacional de Pesquisa e Desenvolvimento Científico e Tecnológico - CNPq / MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional expression regulation by inducing mRNA degradation or translation inhibition. Some miRNAs are known to regulate HLA-G expression, an important immunemodulatory molecule that inhibits both Natural Killer and cytotoxic T cells through interaction with inhibitory receptors. The HLA-G is associated with maternal-fetal tolerance, tissue acceptance in transplants and the progression of tumors. The mechanisms underlying HLA-G expression control are not completely understood, however, its 3’untranslated region (3’UTR) is reported to play an important role on gene regulation influencing mRNA stability and interacting with miRNAs such as miR-148a-3p. In this study, we performed a systematic analysis of all miRNAs that are good candidates to act as HLA-G regulators. In order to determine the miRNAs with the highest potential to influence HLA-G expression, we compared the outputs of three distinct algorithms - miRanda, RNAhybrid and Pita. For this purpose, a method of miRNA inference was developed using Perl scripts to compare and filter results and a scoring system was created in order to evaluate both the binding stability of the miRNA/mRNA interaction and the miRNA specificity to its target sequence. Then, a panel of miRNAs with great potential of controlling HLA-G expression was generated. / MicroRNAs (miRNAs) são pequenos RNAs não codificantes envolvidos na regulação gênica pós-transcricional por meio da degradação da molécula de RNA mensageiro ou da inibição da tradução. Alguns miRNAs foram relatados como sendo responsáveis pela regulação da expressão do gene HLA-G, um importante imunomodulador que inibe a ação de células Natural Killer e células T citotóxicas ao interagir com receptores inibitórios. Este gene está associado à tolerância maternofetal, aceitação de tecidos após transplantes e progressão de tumores. Os mecanismos subjacentes à regulação da expressão de HLA-G não foram completamente elucidados, mas sabe-se que sua região 3’ não traduzida (3’NT) possui um papel importante na regulação gênica tanto por manter a estabilidade da molécula de mRNA quanto por interagir com miRNAs como miR-148a-3p. Neste estudo, foram inferidos miRNAs que são bons candidatos para atuarem como reguladores do gene HLA-G. Para determinar os miRNAs com o maior potencial de operarem no controle pós-transcricional dos níveis de HLA-G, comparamos os resultados de três algoritmos distintos – miRanda, RNAhybrid e Pita. Para tanto, foi desenvolvida uma estratégia de inferência de miRNAs que utiliza scripts em Perl para comparação e filtragem dos dados e um sistema de pontuação que permite avaliar tanto a estabilidade da interação miRNA/mRNA quanto a especificidade do miRNA à sua sequência alvo. Assim, um painel confiável de miRNAs com grande possibilidade de influenciar a expressão de HLA-G foi gerado considerando as regiões polimórficas e não polimórficas da região 3’NT do gene HLA-G individualmente.

Catharanthus roseus (Madagascar periwinkle), is a well-known medicinal plant that produces a vast array of terpenoid indole alkaloids (TIAs), including two anticancer compounds vinblastine and vincristine. Industrial scale production of TIAs is hampered by the difficulties of total chemical synthesis of these compounds and the fragmented knowledge on TIA pathway. Transcriptional regulation of the TIA biosynthetic pathway has not been thoroughly investigated in Catharanthus and only a few structural genes have been identified as the targets of two master regulators: the basic helix-loop-helix (bHLH) transcription factor (TF) CrMYC2 and APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF), ORCA3. Next generation sequencing (NGS) has been used as a tool to isolate novel genes encoding enzymes and regulators of TIA pathway in Catharanthus. In this dissertation, I have performed the transcriptomic analysis of transgenic Catharanthus hairy roots ectopically expressing a dominant repressive form of CrMYC2 or ORCA3 in order to understand their potential impact on the TIA transcriptional regulatory network and to identify and characterize novel target(s) of these two key TFs. MYC2 acts as regulatory hub involved in diverse aspects of plant growth, development, and specialized metabolite biosynthesis by coordinating the crosstalk among different phytohormone signals. CrMYC2 was initially identified in Catharanthus as a regulator of ORCA3. CrMYC2 transactivates ORCA3 by binding to the T/G-box in jasmonate-responsive element (JRE) of ORCA3 promoter. RNA interference (RNAi) mediated knockdown of CrMYC2 strongly reduced TIA accumulation in Catharanthus cell suspension culture. However, the potential influence of CrMYC2 on the expression of other regulatory and structural genes in the TIA pathway remains poorly understood. Transcriptomic analyses revealed that CrMYC2 plays an essential role in JA-induced gene expression and the differentially expressed genes are involved in diverse aspects of growth and development as well as abiotic and biotic stress responses in Catharanthus. Additionally, the expression of genes related to auxin, ethylene, and abscisic acid signaling cascades were affected in hairy roots with modified CrMYC2 expression, suggesting this TF mediates cross-talk between JA and other phytohormones. Surprisingly, overexpression of CrMYC2 resulted in repressed expression of TIA pathway genes in transgenic hairy roots. Expressions of key activators of indole and iridoid pathway were downregulated whereas expression of repressors were upregulated in CrMYC2 hairy roots. Activators (i.e. CrMYC2 and ORCA3) and repressors (i.e. G-box binding factors; GBFs) have been isolated and characterized for their role in regulation of TIA pathway. However, the interconnection between those regulators and the underlying molecular mechanism has not been throughly studied. I identified (i) the interaction of CrMYC2 with CrGBFs and (ii) how this cross-family transcription factor interactions fine-tunes TIA biosynthesis in Catharanthus. The expression profiles of CrMYC2 and CrGBFs were highly correlated in different tissues and in response JA. Moreover, CrMYC2 interacted with CrGBF1 and CrGBF2 in both yeast and plant cells. CrGBF1 and CrGBF2 could form homo- and hetero-dimer which bound T/G-box elements of TIA pathway gene promoters. In plant cells, CrGBF1 antagonizes the activity of CrMYC2 on target promoters in a dosage dependent-manner. Similarly, CrMYC2 can overcome CrGBF1-mediated repression of target promoters in a dosage dependent manner. ORCA3 is another major regulator of TIA biosynthesis in Catharanthus. The transcriptomic analysis of ORCA3 transgenic hairy roots revealed (i) the effect of ORCA3 on newly identified TIA pathway biosynthetic enzymes; (ii) identify the potential effect of ORCA3 on three biological processes: abiotic stress response, plant secondary metabolic process, and response to hormonal stimulus; and (iii) the identification of potential regulator(s) of TIA biosynthesis using ORCA3 based co-expression analysis.