The putative Leishmania telomerase RNA (leishTER) undergoes trans-splicing and contains a conserved template sequence

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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.

Komparative Analyse der humanen Chromosomenregion 11p15.3 um die Gene LMO1,TUB und der orthologen Region in der Maus

Ziel der vorliegenden Arbeit war die vergleichende Sequenzierung und nachfolgende Analyse des syntänen chromosomalen Abschnitts auf dem kurzen Arm des humanen Chromosoms 11 in der Region 11p15.3 mit den Genen LMO1, TUB und dem orthologen Genomabschnitt der Maus auf Chromosom 7 F2. Die im Rahmen dieser Arbeit durchgeführte Kartierung dieser beiden chromosomalen Bereiche ermöglichte die Komplettierung einer genomischen Karte auf insgesamt über eine Megabase, die im Kooperationssequenzierprojekt der Universitäts-Kinderklinik und dem Institut für Molekulargenetik in Mainz erstellt wurde. Mit Hilfe von 28 PAC- und Cosmid-Klonen konnten in dieser Arbeit 383 kb an genomischer DNA des Menschen und mit sechs BAC- und PAC-Klonen 412 kb an genomischer DNA der Maus dargestellt werden. Dies ermöglichte erstmals die exakte Festlegung der Reihenfolge der in diesem chromosomalen Abschnitt enthaltenen Gene und die genaue Kartierung von acht STS-Markern des Menschen, bzw. vier STS-Sonden der Maus. Es zeigte sich dabei, dass die chromosomale Orientierung telomer-/centromerwärts des orthologen Bereichs in der Maus im Vergleich zum Menschen in invertierter Ausrichtung vorliegt. Die Sequenzierung von drei humanen Klonen ermöglichte die Bestimmung von 319.119 bp an zusammenhängender genomischer DNA. Dadurch konnte die genaue Lokalisation und Strukturaufklärung der Gene LMO1, ein putatives Tumorsuppressorgen, das mit der Entstehung von Leukämien assoziiert ist, und TUB, ein Transkriptionsmodulator, der in die Fettstoffwechselregulation involviert ist, vorgenommen werden. Für das murine Genom wurden 412.827 bp an neuer DNA-Sequenz durch Sequenzierung von ebenfalls drei Klonen generiert. Der im Vergleich zum Menschen ca. 100 kb größere Genombereich beinhaltete zudem die neuen Gene Stk33 und Eif3. Es handelte sich dabei um zwei Gene, die erst im Rahmen dieser Arbeit entdeckt und charakterisiert wurden. Die parallele Bearbeitung beider Genombereiche ermöglichte eine umfassende komparative Analyse nach kodierenden, funktionellen und strukturgebenden Sequenzabschnitten in beiden Spezies. Es konnten dabei für beide Organismen die Exon-Intron-Strukturen der Gene LMO1/Lmo1 und TUB/Tub geklärt. Zudem konnten vier neue Exons und zwei neue speziesspezifischer Spleißvarianten für TUB/Tub beschrieben werden. Die Identifizierung dieser neuen Spleißvarianten offenbart neue Möglichkeiten für alternative Regulation und Funktion, oder für eine veränderte Proteinstruktur, die weitere Erklärungsansätze für die Entstehung der mit diesen Genen assoziierten Erkrankungen zulässt. In der sequenzierten, größeren Genomsequenz der Maus konnte in den flankierenden, nicht mit der sequenzierten Humansequenz überlappenden Bereich das neue Gen Eif3 in seiner Exon-Intron-Struktur und die beiden letzten Exons 11 und 12 des Gens Stk33 kartiert und charakterisiert werden. Die umfangreiche Sequenzanalyse beider sequenzierter Genombereiche ergab für den Abschnitt des Menschen insgesamt 229 potentielle Exonsequenzen und für den Bereich der Maus 527 mögliche Exonbereiche. Davon konnten beim Menschen explizit 21 Exons und bei der Maus 31 Exons als exprimierte Bereiche identifiziert und experimentell mittels RT-PCR, bzw. durch cDNA-Sequenzierung verifiziert werden. Diese Abschnitte beschrieben nicht nur die Exonbereiche der oben genannten vier Gene, sondern konnten auch neuen nicht weiter definierten EST-Sequenzen zugeordnet werden. Mittels des Interspeziesvergleiches war darüber hinaus auch die Analyse der nichtkodierenden Intergen-Bereiche möglich. So konnten beispielsweise im ersten Intron des LMO1/Lmo1 sieben Sequenzbereiche mit Konservierungen von ca. 90% bestimmt werden. Auch die Charakterisierung von Promotor- und putativ regulatorischen Sequenzabschnitten konnte mit Hilfe unterschiedlicher bioinformatischer Analyse-Tools durchgeführt werden. Die konservierten Sequenzbereiche der DNA zeigen im Durchschnitt eine Homologie von mehr als 65% auf. Auch die Betrachtung der Genomorganisation zeigte Gemeinsamkeiten, die sich meist nur in ihrer graduellen Ausprägung unterschieden. So weist ein knapp 80 kb großer Bereich proximal zum humanen TUB-Gen einen deutlich erhöhten AT-Gehalt auf, der ebenso im murinen Genom nur in verkürzter Version und schwächer ausgeprägt in Erscheinung tritt. Die zusätzliche Vergleichsanalyse mit einer weiteren Spezies, den orthologen Genomabschnitten von Fugu, zeigte, dass es sich bei den untersuchten Genen LMO1 und TUB um sehr konservierte und evolutiv alte Gene handelt, deren genomisches Organisationsmuster sich auch bei den paralogen Genfamilienmitglieder innerhalb derselben Spezies wiederfindet. Insgesamt konnte durch die Kartierung, Sequenzierung und Analyse eine umfassende Datenbasis für die betrachtete Genomregion und die beschriebenen Gene generiert werden, die für zukünftige Untersuchungen und Fragestellungen wertvolle Informationen bereithält.

Unilateral focal polymicrogyria in a patient with classical Aarskog-Scott syndrome due to a novel missense mutation in an evolutionary conserved RhoGEF domain of the faciogenital dysplasia gene FGD1

Faciogenital dysplasia or Aarskog-Scott syndrome (AAS) is an X-linked disorder characterized by craniofacial, skeletal, and urogenital malformations and short stature. Mutations in the only known causative gene FGD1 are found in about one-fifth of the cases with the clinical diagnosis of AAS. FGD1 is a guanine nucleotide exchange factor (GEF) that specifically activates the Rho GTPase Cdc42 via its RhoGEF domain. The Cdc42 pathway is involved in skeletal formation and multiple aspects of neuronal development. We describe a boy with typical AAS and, in addition, unilateral focal polymicrogyria (PMG), a feature hitherto unreported in AAS. Sequencing of the FGD1 gene in the index case and his mother revealed the presence of a novel mutation (1396A>G; M466V), located in the evolutionary conserved alpha-helix 4 of the RhoGEF domain. M466V was not found in healthy family members, in >300 healthy controls and AAS patients, and has not been reported in the literature or mutation databases to date, indicating that this novel missense mutation causes AAS, and possibly PMG. Brain cortex malformations such as PMG could be initiated by mutations in the evolutionary conserved RhoGEF domain of FGD1, by perturbing the signaling via Rho GTPases such as Cdc42 known to cause brain malformation.

Characterization of the MLC box: A conserved {\it cis\/}-acting DNA element present in the chicken MLC1s/1c gene

The slow/cardiac alkali myosin light chain (MLC1s/1c) is a member of a multigene family whose protein products are essential for activation of the myosin ATPase. In the adult, the MLC1s/1c isoform is expressed in both cardiac and slow-twitch skeletal muscles, while it is expressed by all skeletal muscles during development.^ To elucidate the molecular mechanisms that underlie the transcriptional regulation of MLC1s/1c gene expression, the immediate 5$\sp\prime$ flanking region of the gene was isolated and shown to be capable of directing reporter gene expression. Analysis of this region revealed a 110 bp muscle-specific enhancer that includes a myocyte-specific enhancer-binding factor 2 (MEF-2) site, E-boxes, which are potential binding sites for the basic-helix-loop-helix proteins such as MyoD, and a MLC box. The focus of the thesis was to identify the role of the MLC box in expression of the MLC1s/1c gene.^ The MLC box is a member of the family of CArG box containing cis-acting DNA elements. Mutagenesis showed that the MLC box is necessary, but not sufficient, for the expression of a reporter gene linked to the 5$\sp\prime$ flanking region of the MLC1s/1c gene. Linker scanner and site-directed mutagenesis identified a number of potential sites within the 110 bp muscle-specific enhancer that may cooperate with the MLC box. These are the MEF-2 site, the E-box site, and a 10 bp element located upstream of the MEF-2 site that does not have sequence similarity with any known cis-acting element. The MLC box is capable of binding to factors present in muscle nuclear extracts, as well as to human recombinant serum response factor (SRF). Binding of SRF to the MLC box was correlated with the ability of the 5$\sp\prime$ flanking region of the MLC1s/1c gene to drive reporter gene expression. Results suggest a model in which binding of SRF to the MLC box activates expression of the MLC1s/1c gene while binding of the factors present in the nuclear extracts suppresses the expression of the gene. (Abstract shortened with permission of author.) ^

Distinct and conserved transcriptomic changes during nematode-induced giant cell development in tomato compared with Arabidopsis: a functional role for gene repression

Root-knot nematodes (RKNs) induce giant cells (GCs) from root vascular cells inside the galls. Accompanying molecular changes as a function of infection time and across different species, and their functional impact, are still poorly understood. Thus, the transcriptomes of tomato galls and laser capture microdissected (LCM) GCs over the course of parasitism were compared with those of Arabidopsis, and functional analysis of a repressed gene was performed. Microarray hybridization with RNA from galls and LCM GCs, infection-reproduction tests and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) transcriptional profiles in susceptible and resistant (Mi-1) lines were performed in tomato. Tomato GC-induced genes include some possibly contributing to the epigenetic control of GC identity. GC-repressed genes are conserved between tomato and Arabidopsis, notably those involved in lignin deposition. However, genes related to the regulation of gene expression diverge, suggesting that diverse transcriptional regulators mediate common responses leading to GC formation in different plant species. TPX1, a cell wall peroxidase specifically involved in lignification, was strongly repressed in GCs/galls, but induced in a nearly isogenic Mi-1 resistant line on nematode infection. TPX1 overexpression in susceptible plants hindered nematode reproduction and GC expansion. Time-course and cross-species comparisons of gall and GC transcriptomes provide novel insights pointing to the relevance of gene repression during RKN establishment.

Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

Autoregulation at the level of mRNA 3′ end formation of the suppressor of forked gene of Drosophila melanogaster is conserved in Drosophila virilis

The Drosophila melanogaster Suppressor of forked [Su(f)] protein shares homology with the yeast RNA14 protein and the 77-kDa subunit of human cleavage stimulation factor, which are proteins involved in mRNA 3′ end formation. This suggests a role for Su(f) in mRNA 3′ end formation in Drosophila. The su(f) gene produces three transcripts; two of them are polyadenylated at the end of the transcription unit, and one is a truncated transcript, polyadenylated in intron 4. Using temperature-sensitive su(f) mutants, we show that accumulation of the truncated transcript requires wild-type Su(f) protein. This suggests that the Su(f) protein autoregulates negatively its accumulation by stimulating 3′ end formation of the truncated su(f) RNA. Cloning of su(f) from Drosophila virilis and analysis of its RNA profile suggest that su(f) autoregulation is conserved in this species. Sequence comparison between su(f) from both species allows us to point out three conserved regions in intron 4 downstream of the truncated RNA poly(A) site. These conserved regions include the GU-rich downstream sequence involved in poly(A) site definition. Using transgenes truncated within intron 4, we show that sequence up to the conserved GU-rich domain is sufficient for production of the truncated RNA and for regulation of this production by su(f). Our results indicate a role of su(f) in the regulation of poly(A) site utilization and an important role of the GU-rich sequence for this regulation to occur.

The Drosophila LIM-only gene, dLMO, is mutated in Beadex alleles and might represent an evolutionarily conserved function in appendage development

The process of wing patterning involves precise molecular mechanisms to establish an organizing center at the dorsal–ventral boundary, which functions to direct the development of the Drosophila wing. We report that misexpression of dLMO, a Drosophila LIM-only protein, in specific patterns in the developing wing imaginal disc, disrupts the dorsal–ventral (D-V) boundary and causes errors in wing patterning. When dLMO is misexpressed along the anterior–posterior boundary, extra wing outgrowth occurs, similar to the phenotype seen when mutant clones lacking Apterous, a LIM homeodomain protein known to be essential for normal D-V patterning of the wing, are made in the wing disc. When dLMO is misexpressed along the D-V boundary in third instar larvae, loss of the wing margin is observed. This phenotype is very similar to the phenotype of Beadex, a long-studied dominant mutation that we show disrupts the dLMO transcript in the 3′ untranslated region. dLMO normally is expressed in the wing pouch of the third instar wing imaginal disc during patterning. A mammalian homolog of dLMO is expressed in the developing limb bud of the mouse. This indicates that LMO proteins might function in an evolutionarily conserved mechanism involved in patterning the appendages.

Independent deletions of a pathogen-resistance gene in Brassica and Arabidopsis

Plant disease resistance (R) genes confer race-specific resistance to pathogens and are genetically defined on the basis of intra-specific functional polymorphism. Little is known about the evolutionary mechanisms that generate this polymorphism. Most R loci examined to date contain alternate alleles and/or linked homologs even in disease-susceptible plant genotypes. In contrast, the resistance to Pseudomonas syringae pathovar maculicola (RPM1) bacterial resistance gene is completely absent (rpm1-null) in 5/5 Arabidopsis thaliana accessions that lack RPM1 function. The rpm1-null locus contains a 98-bp segment of unknown origin in place of the RPM1 gene. We undertook comparative mapping of RPM1 and flanking genes in Brassica napus to determine the ancestral state of the RPM1 locus. We cloned two B. napus RPM1 homologs encoding hypothetical proteins with ≈81% amino acid identity to Arabidopsis RPM1. Collinearity of genes flanking RPM1 is conserved between B. napus and Arabidopsis. Surprisingly, we found four additional B. napus loci in which the flanking marker synteny is maintained but RPM1 is absent. These B. napus rpm1-null loci have no detectable nucleotide similarity to the Arabidopsis rpm1-null allele. We conclude that RPM1 evolved before the divergence of the Brassicaceae and has been deleted independently in the Brassica and Arabidopsis lineages. These results suggest that functional polymorphism at R gene loci can arise from gene deletions.

Conserved regulatory elements of the promoter sequence of the gene rpoH of enteric bacteria

The rpoH regulatory region of different members of the enteric bacteria family was sequenced or downloaded from GenBank and compared. In addition, the transcriptional start sites of rpoH of Yersinia frederiksenii and Proteus mirabilis, two distant members of this family, were determined. Sequences similar to the σ70 promoters P1, P4 and P5, to the σE promoter P3 and to boxes DnaA1, DnaA2, cAMP receptor protein (CRP) boxes CRP1, CRP2 and box CytR present in Escherichia coli K12, were identified in sequences of closely related bacteria such as: E.coli, Shigella flexneri, Salmonella enterica serovar Typhimurium, Citrobacter freundii, Enterobacter cloacae and Klebsiella pneumoniae. In more distant bacteria, Y.frederiksenii and P.mirabilis, the rpoH regulatory region has a distal P1-like σ70 promoter and two proximal promoters: a heat-induced σE-like promoter and a σ70 promoter. Sequences similar to the regulatory boxes were not identified in these bacteria. This study suggests that the general pattern of transcription of the rpoH gene in enteric bacteria includes a distal σ70 promoter, >200 nt upstream of the initiation codon, and two proximal promoters: a heat-induced σE-like promoter and a σ70 promoter. A second proximal σ70 promoter under catabolite-regulation is probably present only in bacteria closely related to E.coli.

Divergent and conserved features in the spatial expression of the Drosophila pseudoobscura esterase-5B gene and the esterase-6 gene of Drosophila melanogaster

The regulatory regions of homologous genes encoding esterase 6 (Est-6) of Drosophila melanogaster and esterase 5B (Est-5B) of Drosophila pseudoobscura show very little similarity. We have undertaken a comparative study of the pattern of expression directed by the Est-5B and Est-6 5′-flanking DNA to attempt to reveal conserved elements regulating tissue-specific expression in adults. Esterase regulatory sequences were linked to a lacZ reporter gene and transformed into D. melanogaster embryos. Est-5B, 5′ upstream elements, give rise to a β-galactosidase expression pattern that coincides with the wild-type expression of Est-5B in D. pseudoobscura. The expression patterns of the Est-5B/lacZ construct are different from those of a fusion gene containing the upstream region of Est-6. Common sites of expression for both kinds of constructs are the third segment of antenna, the maxillary palps, and salivary glands. In vitro deletion mutagenesis has shown that the two genes have a different organization of regulatory elements controlling expression in both the third segment of antenna and maxillary palps. The results suggest that the conservation of the expression pattern in genes that evolved from a common ancestor may not be accompanied by preservation of the corresponding cis-regulatory elements.

A gene family required for human germ cell development evolved from an ancient meiotic gene conserved in metazoans

The Deleted in AZoospermia (DAZ) genes encode potential RNA-binding proteins that are expressed exclusively in prenatal and postnatal germ cells and are strong candidates for human fertility factors. Here we report the identification of an additional member of the DAZ gene family, which we have called BOULE. With the identification of this gene, it is clear that the human DAZ gene family contains at least three members: DAZ, a Y-chromosome gene cluster that arose 30–40 million years ago and whose deletion is linked to infertility in men; DAZL, the “father” of DAZ, a gene that maps to human chromosome 3 and has homologs required for both female and male germ cell development in other organisms; and BOULE, a gene that we propose is the “grandfather” of DAZ and maps to human chromosome 2. Human and mouse BOULE resemble the invertebrate meiotic regulator Boule, the proposed ortholog of DAZ, in sequence and expression pattern and hence likely perform a similar meiotic function. In contrast, the previously identified human DAZ and DAZL are expressed much earlier than BOULE in prenatal germ stem cells and spermatogonia; DAZL also is expressed in female germ cells. These data suggest that homologs of the DAZ gene family can be grouped into two subfamilies (BOULE and DAZL) and that members of the DAZ family evolved from an ancestral meiotic regulator, Boule, to assume distinct, yet overlapping, functions in germ cell development.