In silico identification of conserved intercoding sequences in Leishmania genomes: Unraveling putative cis-regulatory elements

In silico analyses of Leishmania spp. genome data are a powerful resource to improve the understanding of these pathogens' biology. Trypanosomatids such as Leishmania spp. have their protein-coding genes grouped in long polycistronic units of functionally unrelated genes. The control of gene expression happens by a variety of posttranscriptional mechanisms. The high degree of synteny among Leishmania species is accompanied by highly conserved coding sequences (CDS) and poorly conserved intercoding untranslated sequences. To identify the elements involved in the control of gene expression, we conducted an in silico investigation to find conserved intercoding sequences (CICS) in the genomes of L major, L infantum, and L braziliensis. We used a combination of computational tools, such as Linux-Shell, PERL and R languages, BLAST, MSPcrunch, SSAKE, and Pred-A-Term algorithms to construct a pipeline which was able to: (i) search for conservation in target-regions, (ii) eliminate CICS redundancy and mask repeat elements, (iii) predict the mRNA's extremities, (iv) analyze the distribution of orthologous genes within the generated LeishCICS-clusters, (v) assign GO terms to the LeishCICS-clusters. and (vi) provide statistical support for the gene-enrichment annotation. We associated the LeishCICS-cluster data, generated at the end of the pipeline, with the expression profile oft. donovani genes during promastigote-amastigote differentiation, as previously evaluated by others (GEO accession: GSE21936). A Pearson's correlation coefficient greater than 0.5 was observed for 730 LeishCICS-clusters containing from 2 to 17 genes. The designed computational pipeline is a useful tool and its application identified potential regulatory cis elements and putative regulons in Leishmania. (C) 2012 Elsevier B.V. All rights reserved.

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.

Computational identification and functional validation of regulatory motifs in cartilage-expressed genes.

Chondrocyte gene regulation is important for the generation and maintenance of cartilage tissues. Several regulatory factors have been identified that play a role in chondrogenesis, including the positive transacting factors of the SOX family such as SOX9, SOX5, and SOX6, as well as negative transacting factors such as C/EBP and delta EF1. However, a complete understanding of the intricate regulatory network that governs the tissue-specific expression of cartilage genes is not yet available. We have taken a computational approach to identify cis-regulatory, transcription factor (TF) binding motifs in a set of cartilage characteristic genes to better define the transcriptional regulatory networks that regulate chondrogenesis. Our computational methods have identified several TFs, whose binding profiles are available in the TRANSFAC database, as important to chondrogenesis. In addition, a cartilage-specific SOX-binding profile was constructed and used to identify both known, and novel, functional paired SOX-binding motifs in chondrocyte genes. Using DNA pattern-recognition algorithms, we have also identified cis-regulatory elements for unknown TFs. We have validated our computational predictions through mutational analyses in cell transfection experiments. One novel regulatory motif, N1, found at high frequency in the COL2A1 promoter, was found to bind to chondrocyte nuclear proteins. Mutational analyses suggest that this motif binds a repressive factor that regulates basal levels of the COL2A1 promoter.

Molecular cloning and characterization of a ripening-induced polygalacturonase related to papaya fruit softening

Pulp softening is one of the most remarkable changes during ripening of papaya (Carica papaya) fruit and it is a major cause for post-harvest losses. Although cell wall catabolism has a major influence on papaya fruit, quality information on the gene products involved in this process is limited. A full-length polygalacturonase cDNA (cpPG) was isolated from papaya pulp and used to study gene expression and enzyme activity during normal and ethylene-induced ripening and after exposure of the fruit to 1-MCP. Northern-blot analysis demonstrated that cpPG transcription was strongly induced during ripening and was highly ethylene-dependent. The accumulation of cpPG transcript was paralleled by enzyme activity, and inversely correlated to the pulp firmness. Preliminary in silica analysis of the cpPG genomic sequence revealed the occurrence of putative regulatory motifs in the promoter region that may help to explain the effects of plant hormones and non-abiotic stresses on papaya fruit firmness. This newly isolated cpPG is an important candidate for functional characterization and manipulation to control the process of pulp softening during papaya ripening. (C) 2009 Elsevier Masson SAS. All rights reserved.

Genomic structure and ecdysone regulation of the prophenoloxidase 1 gene in the malaria vector Anopheles gambiae

Prophenoloxidase, a melanin-synthesizing enzyme, is considered to be an important arthropod immune protein. In mosquitoes, prophenoloxidase has been shown to be involved in refractory mechanisms against malaria parasites. In our study we used Anopheles gambiae, the most important human malaria vector, to characterize the first arthropod prophenoloxidase gene at the genomic level. The complete nucleotide sequence, including the immediate 5′ flanking sequence (−855 bp) of the prophenoloxidase 1 gene, was determined. The gene spans 10 kb and is composed of five exons and four introns coding for a 2.5-kb mRNA. In the 5′ flanking sequence, we found several putative regulatory motifs, two of which were identified as ecdysteroid regulatory elements. Electrophoretic mobility gel-shift assays and supershift assays demonstrated that the Aedes aegypti ecdysone receptor/Ultraspiracle nuclear receptor complex, and, seemingly, the endogenous Anopheles gambiae nuclear receptor complex, was able to bind one of the ecdysteroid response elements. Furthermore, 20-hydroxyecdysone stimulation was shown to up-regulate the transcription of the prophenoloxidase 1 gene in an A. gambiae cell line.

beta2 knockout mice develop parenchymal iron overload: A putative role for class I genes of the major histocompatibility complex in iron metabolism.

Hemochromatosis (HC) is an inherited disorder of iron absorption, mapping within the human major histocompatibility complex (MHC). We have identified a multigene system in the murine MHC that contains excellent candidates for the murine equivalent of the human HC locus and implicate nonclassical class I genes in the control of iron absorption. This gene system is characterized by multiple copies of two head-to-head genes encoded on opposite strands and driven by one common regulatory motif. This regulatory motif has a striking homology to the promoter region of the beta-globin gene, a gene obviously involved in iron metabolism and hence termed beta-globin analogous promoter (betaGAP). Upstream of the betaGAP sequence are nonclassical class I genes. At least one of these nonclassical class I genes, Q2, is expressed in the gastrointestinal tract, the primary site of iron absorption. Also expressed in the gastrointestinal tract and downstream of the betaGAP motif is a second set of putative genes, termed Hephaestus (HEPH). Based on these observations, we hypothesized that the genes that seem to be controlled by the betaGAP regulatory motifs would be responsible for the control of Fe absorption. As a test of this hypothesis, we predicted that mice which have altered expression of class I gene products, the beta2-microglobulin knockout mice, [beta2m(-/-)], would develop Fe overload. This prediction was confirmed, and these results indicate beta2m-associated proteins are involved in the control of intestinal Fe absorption.

Identification of the regulatory subunit of Arabidopsis thaliana acetohydroxyacid synthase and reconstitution with its catalytic subunit

Acetohydroxyacid synthase (EC 4.1.3.18; AHAS) catalyzes the initial step in the formation of the branched-chain amino acids. The enzyme from most bacteria is composed of a catalytic subunit, and a smaller regulatory subunit that is required for full activity and for sensitivity to feedback regulation by valine. A similar arrangement was demonstrated recently for yeast AHAS, and a putative regulatory subunit of tobacco AHAS has also been reported. In this latter case, the enzyme reconstituted from its purified subunits remained insensitive to feedback inhibition, unlike the enzyme extracted from native plant sources. Here we have cloned, expressed in Escherichia coil, and purified the AHAS regulatory subunit of Ambidopsis thaliana. Combining the protein with the purified A. thaliana catalytic subunit results in an activity stimulation that is sensitive to inhibition by valine, leucine, and isoleucine. Moreover, there is a strong synergy between the effects of leucine and valine, which closely mimics the properties of the native enzyme. The regulatory subunit contains a sequence repeat of approximately 180 residues, and we suggest that one repeat binds leucine while the second binds valine or isoleucine. This proposal is supported by reconstitution studies of the individual repeats, which were also cloned, expressed, and purified. The structure and properties of the regulatory subunit are reminiscent of the regulatory domain of threonine deaminase (EC 4.2.1.16), and it is suggested that the two proteins are evolutionarily related.

Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera

Background: The insect exoskeleton provides shape, waterproofing, and locomotion via attached somatic muscles. The exoskeleton is renewed during molting, a process regulated by ecdysteroid hormones. The holometabolous pupa transforms into an adult during the imaginal molt, when the epidermis synthe3sizes the definitive exoskeleton that then differentiates progressively. An important issue in insect development concerns how the exoskeletal regions are constructed to provide their morphological, physiological and mechanical functions. We used whole-genome oligonucleotide microarrays to screen for genes involved in exoskeletal formation in the honeybee thoracic dorsum. Our analysis included three sampling times during the pupal-to-adult molt, i.e., before, during and after the ecdysteroid-induced apolysis that triggers synthesis of the adult exoskeleton. Results: Gene ontology annotation based on orthologous relationships with Drosophila melanogaster genes placed the honeybee differentially expressed genes (DEGs) into distinct categories of Biological Process and Molecular Function, depending on developmental time, revealing the functional elements required for adult exoskeleton formation. Of the 1,253 unique DEGs, 547 were upregulated in the thoracic dorsum after apolysis, suggesting induction by the ecdysteroid pulse. The upregulated gene set included 20 of the 47 cuticular protein (CP) genes that were previously identified in the honeybee genome, and three novel putative CP genes that do not belong to a known CP family. In situ hybridization showed that two of the novel genes were abundantly expressed in the epidermis during adult exoskeleton formation, strongly implicating them as genuine CP genes. Conserved sequence motifs identified the CP genes as members of the CPR, Tweedle, Apidermin, CPF, CPLCP1 and Analogous-to-Peritrophins families. Furthermore, 28 of the 36 muscle-related DEGs were upregulated during the de novo formation of striated fibers attached to the exoskeleton. A search for cis-regulatory motifs in the 5′-untranslated region of the DEGs revealed potential binding sites for known transcription factors. Construction of a regulatory network showed that various upregulated CP- and muscle-related genes (15 and 21 genes, respectively) share common elements, suggesting co-regulation during thoracic exoskeleton formation. Conclusions: These findings help reveal molecular aspects of rigid thoracic exoskeleton formation during the ecdysteroid-coordinated pupal-to-adult molt in the honeybee.

Regulatory networks of Neisseria meningitidis and their implications for pathogenesis

Neisseria meningitidis, the leading cause of bacterial meningitis, can adapt to different host niches during human infection. Both transcriptional and post-transcriptional regulatory networks have been identified as playing a crucial role for bacterial stress responses and virulence. We investigated the N. meningitidis transcriptional landscape both by microarray and by RNA sequencing (RNAseq). Microarray analysis of N. meningitidis grown in the presence or absence of glucose allowed us to identify genes regulated by carbon source availability. In particular, we identified a glucose-responsive hexR-like transcriptional regulator in N. meningitidis. Deletion analysis showed that the hexR gene is accountable for a subset of the glucose-responsive regulation, and in vitro assays with the purified protein showed that HexR binds to the promoters of the central metabolic operons of meningococcus, by targeting a DNA region overlapping putative regulatory sequences. Our results indicate that HexR coordinates the central metabolism of meningococcus in response to the availability of glucose, and N. meningitidis strains lacking the hexR gene are also deficient in establishing successful bacteremia in a mouse model of infection. In parallel, RNAseq analysis of N. meningitidis cultured under standard or iron-limiting in vitro growth conditions allowed us to identify novel small non-coding RNAs (sRNAs) potentially involved in N. meningitidis regulatory networks. Manual curation of the RNAseq data generated a list of 51 sRNAs, 8 of which were validated by Northern blotting. Deletion of selected sRNAs caused attenuation of N. meningitidis infection in a murine model, leading to the identification of the first sRNAs influencing meningococcal bacteraemia. Furthermore, we describe the identification and initial characterization of a novel sRNA unique to meningococcus, closely associated to genes relevant for the intracellular survival of pathogenic Neisseriae. Taken together, our findings could help unravel the regulation of N. meningitidis adaptation to the host environment and its implications for pathogenesis.

Functional analysis of rice bidirectional promoters

Bidirectional promoters regulate adjacent genes organized in a divergent fashion (head to head orientation). Several Reports pertaining to bidirectional promoters on a genomic scale exists in mammals. This work provides the essential background on theoretical and experimental work to carry out a genomic scale analysis of bidirectional promoters in plants. A computational study was performed to identify putative bidirectional promoters and the over-represented cis-regulatory motifs from three sequenced plant genomes: rice (Oryza sativa), Arabidopsis thaliana, and Populus trichocarpa using the Plant Cis-acting Regulatory DNA Elements (PLACE) and PLANT CARE databases. Over-represented motifs along with their possible function were described with the help of a few conserved representative putative bidirectional promoters from the three model plants. By doing so a foundation was laid for the experimental evaluation of bidirectional promoters in plants. A novel Agrobacterium tumefaciens mediated transient expression assay (AmTEA) was developed for young plants of different cereal species and the model dicot Arabidopsis thaliana. AmTEA was evaluated using five promoters (six constructs) and two reporter genes, gus and egfp. Efficacy and stability of AmTEA was compared with stable transgenics using the Arabidopsis DEAD-box RNA helicase family gene promoter. AmTEA was primarily developed to overcome the many problems associated with the development of transgenics and expression studies in plants. Finally a possible mechanism for the bidirectional activity of bidirectional promoters was highlighted. Deletion analysis using promoter-reporter gene constructs identified three rice promoters to be bidirectional. Regulatory elements located in the 5’- untranslated regions (UTR) of one of the genes of the divergent gene pair were found to be responsible for their bidirectional ctivity

Naturally occurring variation in copia expression is due to both element (cis) and host (trans) regulatory variation.

Significant differences in levels of copia [Drosophila long terminal repeat (LTR) retrotransposon] expression exist among six species representing the Drosophila melanogaster species complex (D. melanogaster, Drosophila mauritiana, Drosophila simulans, Drosophila sechellia, Drosophila yakuba, and Drosophila erecta) and a more distantly related species (Drosophila willistoni). These differences in expression are correlated with major size variation mapping to putative regulatory regions of the copia 5' LTR and adjacent untranslated leader region (ULR). Sequence analysis indicates that these size variants were derived from a series of regional duplication events. The ability of the copia LTR-ULR size variants to drive expression of a bacterial chloramphenicol acetyltransferase reporter gene was tested in each of the seven species. The results indicate that both element-encoded (cis) and host-genome-encoded (trans) genetic differences are responsible for the variability in copia expression within and between Drosophila species. This finding indicates that models purporting to explain the dynamics and distribution of retrotransposons in natural populations must consider the potential impact of both element-encoded and host-genome-encoded regulatory variation to be valid. We propose that interelement selection among retrotransposons may provide a molecular drive mechanism for the evolution of eukaryotic enhancers which can be subsequently distributed throughout the genome by retrotransposition.

The effect of male -male competition and its underlying regulatory mechanisms on the electric signal of the gymnotiform fish Brachyhypopomus gauderio

Sexually-selected communication signals can be used by competing males to settle contests without incurring the costs of fighting. The ability to dynamically regulate the signal in a context-dependent manner can further minimize the costs of male aggressive interactions. Such is the case in the gymnotiform fish Brachyhypopomus gauderio, which, by coupling its electric organ discharge (EOD) waveform to endocrine systems with circadian, seasonal, and behavioral drivers, can regulate its signal to derive the greatest reproductive benefit. My dissertation research examined the functional role of the EOD plasticity observed in male B. gauderio and the physiological mechanisms that regulate the enhanced male EOD. To evaluate whether social competition drives the EOD changes observed during male-male interactions, I manipulated the number of males in breeding groups to create conditions that exemplified low and high competition and measured their EOD and steroid hormone levels. My results showed that social competition drives the enhancement of the EOD amplitude of male B. gauderio. In addition, changes in the EOD of males due to changes in their social environment were paralleled by changes in the levels of androgens and cortisol. I also examined the relationship between body size asymmetry, EOD waveform parameters, and aggressive physical behaviors during male-male interactions in B. gauderio, in order to understand more fully the role of EOD waveforms as reliable signals. While body size was the best determinant of dominance in male B. gauderio, EOD amplitude reliably predicted body condition, a composite of length and weight, for fish in good body condition. To further characterize the mechanisms underlying the relationship between male-male interactions and EOD plasticity, I identified the expression of the serotonin receptor 1A, a key player in the regulation of aggressive behavior, in the brains of B. gauderio. I also identified putative regulatory regions in this receptor in B. gauderio and other teleost fish, highlighting the presence of additional plasticity. In conclusion, male-male competition seems to be a strong selective driver in the evolution of the male EOD plasticity in B. gauderio via the regulatory control of steroid hormones and the serotonergic system.

The Effect of male-male competition and its Underlying Regulatory Mechanisms on the Electric Signal of the Gymnotiform fish Brachyhypopomus gauderio

Sexually-selected communication signals can be used by competing males to settle contests without incurring the costs of fighting. The ability to dynamically regulate the signal in a context-dependent manner can further minimize the costs of male aggressive interactions. Such is the case in the gymnotiform fish Brachyhypopomus gauderio, which, by coupling its electric organ discharge (EOD) waveform to endocrine systems with circadian, seasonal, and behavioral drivers, can regulate its signal to derive the greatest reproductive benefit. My dissertation research examined the functional role of the EOD plasticity observed in male B. gauderio and the physiological mechanisms that regulate the enhanced male EOD. To evaluate whether social competition drives the EOD changes observed during male-male interactions, I manipulated the number of males in breeding groups to create conditions that exemplified low and high competition and measured their EOD and steroid hormone levels. My results showed that social competition drives the enhancement of the EOD amplitude of male B. gauderio. In addition, changes in the EOD of males due to changes in their social environment were paralleled by changes in the levels of androgens and cortisol. I also examined the relationship between body size asymmetry, EOD waveform parameters, and aggressive physical behaviors during male-male interactions in B. gauderio, in order to understand more fully the role of EOD waveforms as reliable signals. While body size was the best determinant of dominance in male B. gauderio, EOD amplitude reliably predicted body condition, a composite of length and weight, for fish in good body condition. To further characterize the mechanisms underlying the relationship between male-male interactions and EOD plasticity, I identified the expression of the serotonin receptor 1A, a key player in the regulation of aggressive behavior, in the brains of B. gauderio. I also identified putative regulatory regions in this receptor in B. gauderio and other teleost fish, highlighting the presence of additional plasticity. In conclusion, male-male competition seems to be a strong selective driver in the evolution of the male EOD plasticity in B. gauderio via the regulatory control of steroid hormones and the serotonergic system.

TFEC is a macrophage-restricted member of the microphthalmia-TFE subfamily of basic helix-loop-helix leucine zipper transcription factors

The murine homologue of the TFEC was cloned as part of an analysis of the expression of the microphthalmia-TFE (MiT) subfamily of transcription factors in macrophages. TFEC, which most likely acts as a transcriptional repressor in heterodimers with other MiT family members, was identified in cells of the mononuclear phagocyte lineage, coexpressed,vith all other known MiT subfamily members (Mitf, TFE3, TFEB), Northern blot analysis of several different cell lineages indicated that the expression of murine TFEC (mTFEC) was restricted to macrophages. A 600-bp fragment of the TATA-less putative proximal promoter of TFEC shares features with many known macrophage-specific promoters and preferentially directs luciferase expression in the RAW264.7 macrophage cell line in transient transfection assays. Five of six putative Ets motifs identified in the TFEC promoter bind the macrophage-restricted transcription factor PU,I under in vitro conditions and in transfected 3T3 fibroblasts; the minimal luciferase activity of the TFEC promoter could be induced by coexpression of PU.1 or the related transcription factor Ets-2. The functional importance of the tissue-restricted expression of TFEC and a possible role in macrophage-specific gene regulation require further investigation, but are likely to be linked to the role of the other MiT family members in this lineage.

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.