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.

Regulatory variation in a TBX5 enhancer leads to isolated congenital heart disease

Recent studies have identified the genetic underpinnings of a growing number of diseases through targeted exome sequencing. However, this strategy ignores the large component of the genome that does not code for proteins, but is nonetheless biologically functional. To address the possible involvement of regulatory variation in congenital heart diseases (CHDs), we searched for regulatory mutations impacting the activity of TBX5, a dosage-dependent transcription factor with well-defined roles in the heart and limb development that has been associated with the HoltOram syndrome (hearthand syndrome), a condition that affects 1/100 000 newborns. Using a combination of genomics, bioinformatics and mouse genetic engineering, we scanned approximate to 700 kb of the TBX5 locus in search of cis-regulatory elements. We uncovered three enhancers that collectively recapitulate the endogenous expression pattern of TBX5 in the developing heart. We re-sequenced these enhancer elements in a cohort of non-syndromic patients with isolated atrial and/or ventricular septal defects, the predominant cardiac defects of the HoltOram syndrome, and identified a patient with a homozygous mutation in an enhancer approximate to 90 kb downstream of TBX5. Notably, we demonstrate that this single-base-pair mutation abrogates the ability of the enhancer to drive expression within the heart in vivo using both mouse and zebrafish transgenic models. Given the population-wide frequency of this variant, we estimate that 1/100 000 individuals would be homozygous for this variant, highlighting that a significant number of CHD associated with TBX5 dysfunction might arise from non-coding mutations in TBX5 heart enhancers, effectively decoupling the heart and hand phenotypes of the HoltOram syndrome.

The SUCEST-FUN regulatory network database: Designing an energy grass

Modern sugarcane cultivars are complex hybrids resulting from crosses among several Saccharum species. Traditional breeding methods have been employed extensively in different countries over the past decades to develop varieties with increased sucrose yield and resistance to pests and diseases. Conventional variety improvement, however, may be limited by the narrow pool of suitable genes. Thus, molecular genetics is seen as a promising tool to assist in the process of developing improved varieties. The SUCEST-FUN Project (http://sucest-fun.org) aims to associate function with sugarcane genes using a variety of tools, in particular those that enable the study of the sugarcane transcriptome. An extensive analysis has been conducted to characterise, phenotypically, sugarcane genotypes with regard to their sucrose content, biomass and drought responses. Through the analysis of different cultivars, genes associated with sucrose content, yield, lignin and drought have been identified. Currently, tools are being developed to determine signalling and regulatory networks in grasses, and to sequence the sugarcane genome, as well as to identify sugarcane promoters. This is being implemented through the SUCEST-FUN (http://sucest-fun.org) and GRASSIUS databases (http://grassius.org), the cloning of sugarcane promoters, the identification of cis-regulatory elements (CRE) using Chromatin Immunoprecipitation-sequencing (ChIP-Seq) and the generation of a comprehensive Signal Transduction and Transcription gene catalogue (SUCAST Catalogue).

Whole-genome sequencing of the efficient industrial fuel-ethanol fermentative Saccharomyces cerevisiae strain CAT-1

The Saccharomyces cerevisiae strains widely used for industrial fuel-ethanol production have been developed by selection, but their underlying beneficial genetic polymorphisms remain unknown. Here, we report the draft whole-genome sequence of the S. cerevisiae strain CAT-1, which is a dominant fuel-ethanol fermentative strain from the sugarcane industry in Brazil. Our results indicate that strain CAT-1 is a highly heterozygous diploid yeast strain, and the similar to 12-Mb genome of CAT-1, when compared with the reference S228c genome, contains similar to 36,000 homozygous and similar to 30,000 heterozygous single nucleotide polymorphisms, exhibiting an uneven distribution among chromosomes due to large genomic regions of loss of heterozygosity (LOH). In total, 58 % of the 6,652 predicted protein-coding genes of the CAT-1 genome constitute different alleles when compared with the genes present in the reference S288c genome. The CAT-1 genome contains a reduced number of transposable elements, as well as several gene deletions and duplications, especially at telomeric regions, some correlated with several of the physiological characteristics of this industrial fuel-ethanol strain. Phylogenetic analyses revealed that some genes were likely associated with traits important for bioethanol production. Identifying and characterizing the allelic variations controlling traits relevant to industrial fermentation should provide the basis for a forward genetics approach for developing better fermenting yeast strains.

Distinct mechanisms underlie adaptation of proximal tubule Na+/H+ exchanger isoform 3 in response to chronic metabolic and respiratory acidosis

The Na+/H+ exchanger isoform 3 (NHE3) is essential for HCO3- reabsorption in renal proximal tubules. The expression and function of NHE3 must adapt to acid-base conditions. The goal of this study was to elucidate the mechanisms responsible for higher proton secretion in proximal tubules during acidosis and to evaluate whether there are differences between metabolic and respiratory acidosis with regard to NHE3 modulation and, if so, to identify the relevant parameters that may trigger these distinct adaptive responses. We achieved metabolic acidosis by lowering HCO3- concentration in the cell culture medium and respiratory acidosis by increasing CO2 tension in the incubator chamber. We found that cell-surface NHE3 expression was increased in response to both forms of acidosis. Mild (pH 7.21 +/- 0.02) and severe (6.95 +/- 0.07) metabolic acidosis increased mRNA levels, at least in part due to up-regulation of transcription, whilst mild (7.11 +/- 0.03) and severe (6.86 +/- 0.01) respiratory acidosis did not up-regulate NHE3 expression. Analyses of the Nhe3 promoter region suggested that the regulatory elements sensitive to metabolic acidosis are located between -466 and -153 bp, where two consensus binding sites for SP1, a transcription factor up-regulated in metabolic acidosis, were localised. We conclude that metabolic acidosis induces Nhe3 promoter activation, which results in higher mRNA and total protein level. At the plasma membrane surface, NHE3 expression was increased in metabolic and respiratory acidosis alike, suggesting that low pH is responsible for NHE3 displacement to the cell surface.

Sequence features responsible for intron retention in human

Abstract Background One of the least common types of alternative splicing is the complete retention of an intron in a mature transcript. Intron retention (IR) is believed to be the result of intron, rather than exon, definition associated with failure of the recognition of weak splice sites flanking short introns. Although studies on individual retained introns have been published, few systematic surveys of large amounts of data have been conducted on the mechanisms that lead to IR. Results TTo understand how sequence features are associated with or control IR, and to produce a generalized model that could reveal previously unknown signals that regulate this type of alternative splicing, we partitioned intron retention events observed in human cDNAs into two groups based on the relative abundance of both isoforms and compared relevant features. We found that a higher frequency of IR in human is associated with individual introns that have weaker splice sites, genes with shorter intron lengths, higher expression levels and lower density of both a set of exon splicing silencers (ESSs) and the intronic splicing enhancer GGG. Both groups of retained introns presented events conserved in mouse, in which the retained introns were also short and presented weaker splice sites. Conclusion Although our results confirmed that weaker splice sites are associated with IR, they showed that this feature alone cannot explain a non-negligible fraction of events. Our analysis suggests that cis-regulatory elements are likely to play a crucial role in regulating IR and also reveals previously unknown features that seem to influence its occurrence. These results highlight the importance of considering the interplay among these features in the regulation of the relative frequency of IR.

Promoter regions of Plasmodium vivax are poorly or not recognized by Plasmodium falciparum

Abstract Background Heterologous promoter analysis in Plasmodium has revealed the existence of conserved cis regulatory elements as promoters from different species can drive expression of reporter genes in heterologous transfection assays. Here, the functional characterization of different Plasmodium vivax promoters in Plasmodium falciparum using luciferase as the reporter gene is presented. Methods Luciferase reporter plasmids harboring the upstream regions of the msp1, dhfr, and vir3 genes as well as the full-length intergenic regions of the vir23/24 and ef-1�� genes of P. vivax were constructed and transiently transfected in P. falciparum. Results Only the constructs with the full-length intergenic regions of the vir23/24 and ef-1�� genes were recognized by the P. falciparum transcription machinery albeit to values approximately two orders of magnitude lower than those reported by luc plasmids harbouring promoter regions from P. falciparum and Plasmodium berghei. A bioinformatics approach allowed the identification of a motif (GCATAT) in the ef-1�� intergenic region that is conserved in five Plasmodium species but is degenerate (GCANAN) in P. vivax. Mutations of this motif in the P. berghei ef-1�� promoter region decreased reporter expression indicating it is active in gene expression in Plasmodium. Conclusion Together, this data indicates that promoter regions of P. vivax are poorly or not recognized by the P. falciparum transcription machinery suggesting the existence of P. vivax-specific transcription regulatory elements.

The clinical impact of chromosomal rearrangements with breakpoints upstream of the SOX9 gene: two novel de novo balanced translocations associated with acampomelic campomelic dysplasia

Abstract Background The association of balanced rearrangements with breakpoints near SOX9 [SRY (sex determining region Y)-box 9] with skeletal abnormalities has been ascribed to the presumptive altering of SOX9 expression by the direct disruption of regulatory elements, their separation from SOX9 or the effect of juxtaposed sequences. Case presentation We report on two sporadic apparently balanced translocations, t(7;17)(p13;q24) and t(17;20)(q24.3;q11.2), whose carriers have skeletal abnormalities that led to the diagnosis of acampomelic campomelic dysplasia (ACD; MIM 114290). No pathogenic chromosomal imbalances were detected by a-CGH. The chromosome 17 breakpoints were mapped, respectively, 917���855 kb and 601���585 kb upstream of the SOX9 gene. A distal cluster of balanced rearrangements breakpoints on chromosome 17 associated with SOX9-related skeletal disorders has been mapped to a segment 932���789 kb upstream of SOX9. In this cluster, the breakpoint of the herein described t(17;20) is the most telomeric to SOX9, thus allowing the redefining of the telomeric boundary of the distal breakpoint cluster region related to skeletal disorders to 601���585 kb upstream of SOX9. Although both patients have skeletal abnormalities, the t(7;17) carrier presents with relatively mild clinical features, whereas the t(17;20) was detected in a boy with severe broncheomalacia, depending on mechanical ventilation. Balanced and unbalanced rearrangements associated with disorders of sex determination led to the mapping of a regulatory region of SOX9 function on testicular differentiation to a 517���595 kb interval upstream of SOX9, in addition to TESCO (Testis-specific enhancer of SOX9 core). As the carrier of t(17;20) has an XY sex-chromosome constitution and normal male development for his age, the segment of chromosome 17 distal to the translocation breakpoint should contain the regulatory elements for normal testis development. Conclusions These two novel translocations illustrate the clinical variability in carriers of balanced translocations with breakpoints near SOX9. The translocation t(17;20) breakpoint provides further evidence for an additional testis-specific SOX9 enhancer 517 to 595 kb upstream of the SOX9 gene.

Leptin signaling in Kiss1 neurons arises after pubertal development

The adipocyte-derived hormone leptin is required for normal pubertal maturation in mice and humans and, therefore, leptin has been recognized as a crucial metabolic cue linking energy stores and the onset of puberty. Several lines of evidence have suggested that leptin acts via kisspeptin expressing neurons of the arcuate nucleus to exert its effects. Using conditional knockout mice, we have previously demonstrated that deletion of leptin receptors (LepR) from kisspeptin cells cause no puberty or fertility deficits. However, developmental adaptations and system redundancies may have obscured the physiologic relevance of direct leptin signaling in kisspeptin neurons. To overcome these putative effects, we re-expressed endogenous LepR selectively in kisspeptin cells of mice otherwise null for LepR, using the Cre-loxP system. Kiss1-Cre LepR null mice showed no pubertal development and no improvement of the metabolic phenotype, remaining obese, diabetic and infertile. These mice displayed decreased numbers of neurons expressing Kiss1 gene, similar to prepubertal control mice, and an unexpected lack of re-expression of functional LepR. To further assess the temporal coexpression of Kiss1 and Lepr genes, we generated mice with the human renilla green fluorescent protein (hrGFP) driven by Kiss1 regulatory elements and crossed them with mice that express Cre recombinase from the Lepr locus and the R26-tdTomato reporter gene. No coexpression of Kiss1 and LepR was observed in prepubertal mice. Our findings unequivocally demonstrate that kisspeptin neurons are not the direct target of leptin in the onset of puberty. Leptin signaling in kisspeptin neurons arises only after completion of sexual maturation.

Functional characterization of sugarcane mustang domesticated transposases and comparative diversity in sugarcane, rice, maize and sorghum

Transposable elements (TEs) account for a large portion of plant genomes, particularly in grasses, in which they correspond to 50%-80% of the genomic content. TEs have recently been shown to be a source of new genes and new regulatory networks. The most striking contribution of TEs is referred as "molecular domestication", by which the element coding sequence loses its movement capacity and acquires cellular function. Recently, domesticated transposases known as mustang and derived from the Mutator element have been described in sugarcane. In order to improve our understanding of the function of these proteins, we identified mustang genes from Sorghum bicolor and Zea mays and performed a phenetic analysis to assess the diversity and evolutionary history of this gene family. This analysis identified orthologous groups and showed that mustang genes are highly conserved in grass genomes. We also explored the transcriptional activity of sugarcane mustang genes in heterologous and homologous systems. These genes were found to be ubiquitously transcribed, with shoot apical meristem having the highest expression levels, and were downregulated by phytohormones. Together, these findings suggest the possible involvement of mustang proteins in the maintenance of hormonal homeostasis.

Analysis of plant LTR-retrotransposons at the fine-scale family level reveals individual molecular patterns

Background: Sugarcane is an important crop worldwide for sugar production and increasingly, as a renewable energy source. Modern cultivars have polyploid, large complex genomes, with highly unequal contributions from ancestral genomes. Long Terminal Repeat retrotransposons (LTR-RTs) are the single largest components of most plant genomes and can substantially impact the genome in many ways. It is therefore crucial to understand their contribution to the genome and transcriptome, however a detailed study of LTR-RTs in sugarcane has not been previously carried out. Results: Sixty complete LTR-RT elements were classified into 35 families within four Copia and three Gypsy lineages. Structurally, within lineages elements were similar, between lineages there were large size differences. FISH analysis resulted in the expected pattern of Gypsy/heterochromatin, Copia/euchromatin, but in two lineages there was localized clustering on some chromosomes. Analysis of related ESTs and RT-PCR showed transcriptional variation between tissues and families. Four distinct patterns were observed in sRNA mapping, the most unusual of which was that of Ale1, with very large numbers of 24nt sRNAs in the coding region. The results presented support the conclusion that distinct small RNA-regulated pathways in sugarcane target the lineages of LTR-RT elements. Conclusions: Individual LTR-RT sugarcane families have distinct structures, and transcriptional and regulatory signatures. Our results indicate that in sugarcane individual LTR-RT families have distinct behaviors and can potentially impact the genome in diverse ways. For instance, these transposable elements may affect nearby genes by generating a diverse set of small RNA's that trigger gene silencing mechanisms. There is also some evidence that ancestral genomes contribute significantly different element numbers from particular LTR-RT lineages to the modern sugarcane cultivar genome.

Preimplantation development regulatory pathway construction through a text-mining approach

Background: The integration of sequencing and gene interaction data and subsequent generation of pathways and networks contained in databases such as KEGG Pathway is essential for the comprehension of complex biological processes. We noticed the absence of a chart or pathway describing the well-studied preimplantation development stages; furthermore, not all genes involved in the process have entries in KEGG Orthology, important information for knowledge application with relation to other organisms. Results: In this work we sought to develop the regulatory pathway for the preimplantation development stage using text-mining tools such as Medline Ranker and PESCADOR to reveal biointeractions among the genes involved in this process. The genes present in the resulting pathway were also used as seeds for software developed by our group called SeedServer to create clusters of homologous genes. These homologues allowed the determination of the last common ancestor for each gene and revealed that the preimplantation development pathway consists of a conserved ancient core of genes with the addition of modern elements. Conclusions: The generation of regulatory pathways through text-mining tools allows the integration of data generated by several studies for a more complete visualization of complex biological processes. Using the genes in this pathway as ���seeds��� for the generation of clusters of homologues, the pathway can be visualized for other organisms. The clustering of homologous genes together with determination of the ancestry leads to a better understanding of the evolution of such process.

Modularity and evolutionary constraints in a baculovirus gene regulatory network

Abstract Background The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression. Results We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization. Conclusions Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.