A global analysis of low-complexity regions in the Trypanosoma brucei proteome reveals enrichment in the C-terminus of nucleic acid binding proteins providing potential targets of phosphorylation

Background: Low-complexity regions (LCRs) on proteins have attracted increasing attention recently due to their role in the assembly of membraneless organelles or granules by liquid-liquid phase separation. Several examples of such granules have been shown to sequester RNA and proteins in an inactive state, providing an important mechanism for dynamic post-transcriptional gene regulation. In trypanosome parasites, post-transcriptional control overwhelmingly dominates gene regulation due to the organisation of their genome into polycistronic transcription units. The purpose of the current study was to generate a substantially more comprehensive genome-wide survey of LCRs on trypanosome proteins than currently available . Methods: Using the Shannon's entropy method, provided in the R package 'entropy', we identified LCRs in the proteome of Trypanosoma brucei. Our analysis predicts LCRs and their positional enrichment in distinct protein cohorts and superimposes on this a range of post-translational modifications derived from available experimental datasets. Results: We have identified 8162 LCRs present on 4914 proteins, representing 42% of the proteome, placing Trypanosoma brucei among the eukaryotes with the highest percentage of LCRs . Our results highlight the enrichment of LCRs in the C-terminal region of predicted nucleic acid binding proteins, these acting as favoured sites for potential phosphorylation. Phosphorylation represents 51% of the post-translational modifications present on LCRs compared to 16% on the rest of the proteome. Conclusions: The post-translational modifications of LCRs, and in particular phosphorylation events, could contribute to post-transcriptional gene expression control and the dynamics of protein targeting to membraneless organelles in kinetoplastid parasites.

A global analysis of low-complexity regions in the Trypanosoma brucei proteome reveals enrichment in the C-terminus of nucleic acid binding proteins providing potential targets of phosphorylation

Background: Low-complexity regions (LCRs) on proteins have attracted increasing attention recently due to their role in the assembly of membraneless organelles or granules by liquid-liquid phase separation. Several examples of such granules have been shown to sequester RNA and proteins in an inactive state, providing an important mechanism for dynamic post-transcriptional gene regulation. In trypanosome parasites, post-transcriptional control overwhelmingly dominates gene regulation due to the organisation of their genome into polycistronic transcription units. The purpose of the current study was to generate a substantially more comprehensive genome-wide survey of LCRs on trypanosome proteins than currently available . Methods: Using the Shannon's entropy method, provided in the R package 'entropy', we identified LCRs in the proteome of Trypanosoma brucei. Our analysis predicts LCRs and their positional enrichment in distinct protein cohorts and superimposes on this a range of post-translational modifications derived from available experimental datasets. Results: Our results highlight the enrichment of LCRs in the C-terminal region of predicted nucleic acid binding proteins, these acting as favoured sites for potential phosphorylation. Conclusions: The post-translational modifications of LCRs, and in particular the phosphorylation events, could contribute to post-transcriptional gene expression control and the dynamics of protein targeting to membraneless organelles in kinetoplastid parasites.

The Methyl-CpG-Binding Protein Mbd2 Regulates Susceptibility to Experimental Colitis via Control of CD11c+ Cells and Colonic Epithelium

Methyl-CpG-binding domain-2 (Mbd2) acts as an epigenetic regulator of gene expression, by linking DNA methylation to repressive chromatin structure. Although Mbd2 is widely expressed in gastrointestinal immune cells and is implicated in regulating intestinal cancer, anti-helminth responses and colonic inflammation, the Mbd2-expressing cell types that control these responses are incompletely defined. Indeed, epigenetic control of gene expression in cells that regulate intestinal immunity is generally poorly understood, even though such mechanisms may explain the inability of standard genetic approaches to pinpoint the causes of conditions like inflammatory bowel disease. In this study we demonstrate a vital role for Mbd2 in regulating murine colonic inflammation. Mbd2^−/− mice displayed dramatically worse pathology than wild type controls during dextran sulfate sodium (DSS) induced colitis, with increased inflammatory (IL-1β⁺) monocytes. Profiling of mRNA from innate immune and epithelial cell (EC) populations suggested that Mbd2 suppresses inflammation and pathology via control of innate-epithelial cell crosstalk and T cell recruitment. Consequently, restriction of Mbd2 deficiency to CD11c⁺ dendritic cells and macrophages, or to ECs, resulted in increased DSS colitis severity. Our identification of this dual role for Mbd2 in regulating the inflammatory capacity of both CD11c⁺ cells and ECs highlights how epigenetic control mechanisms may limit intestinal inflammatory responses.

The lung environment controls alveolar macrophage metabolism and responsiveness in type 2 inflammation

Fine control of macrophage activation is needed to prevent inflammatory disease, particularly at barrier sites such as the lungs. However, the dominant mechanisms that regulate the activation of pulmonary macrophages during inflammation are poorly understood. We found that alveolar macrophages (AlvMs) were much less able to respond to the canonical type 2 cytokine IL-4, which underpins allergic disease and parasitic worm infections, than macrophages from lung tissue or the peritoneal cavity. We found that the hyporesponsiveness of AlvMs to IL-4 depended upon the lung environment but was independent of the host microbiota or the lung extracellular matrix components surfactant protein D (SP-D) and mucin 5b (Muc5b). AlvMs showed severely dysregulated metabolism relative to that of cavity macrophages. After removal from the lungs, AlvMs regained responsiveness to IL-4 in a glycolysis-dependent manner. Thus, impaired glycolysis in the pulmonary niche regulates AlvM responsiveness during type 2 inflammation.

Dynamics of Colon Monocyte and Macrophage Activation During Colitis

Background: Macrophages are pivotal in coordinating a range of important processes in the intestines, including controlling intracellular infections and limiting damaging inflammation against the microbiota. However, it is not clear how gut macrophages, relative to recruited blood monocytes and other myeloid cells, contribute to the intestinal inflammatory milieu, nor how macrophages and their monocyte precursors mediate recruitment of other immune cells to the inflamed intestine. Methods: Myeloid cell populations isolated from colonic inflammatory bowel disease (IBD) or murine dextran sulphate sodium (DSS) induced colitis were assessed using flow cytometry and compared to healthy controls. In addition, mRNA expression profiles in human and murine colon samples, and in macrophages and monocytes from healthy and inflamed murine colons, were analysed by quantitative PCR (qPCR) and mRNA microarray. Results: We show that the monocyte:macrophage balance is disrupted in colon inflammation to favour recruitment of CD14+HLA-DRInt cells in humans, and Ly6CHi monocytes in mice. In addition, we identify that murine blood monocytes receive systemic signals enabling increased release of IL-1β prior to egress from the blood into the colon. Further, once within the colon and relative to other myeloid cells, monocytes represent the dominant local source of both IL-1β and TNF. Finally, our data reveal that, independent of inflammation, murine colon macrophages act as a major source of Ccl7 and Ccl8 chemokines that trigger further recruitment of their pro-inflammatory monocyte precursors. Conclusions: Our work suggests that strategies targeting macrophage-mediated monocyte recruitment may represent a promising approach for limiting the chronic inflammation that characterises IBD.

Modulation of dendritic cell alternative activation and function by the vitamin A metabolite retinoic acid

Retinoic acid modulates the functions of IL-4 in alternatively activated DCs

The archetypal T_h2 cytokine IL-4 has previously been shown to alternatively activate murine macrophages and, more recently, dendritic cells (DCs) both in vitro and in vivo. IL-4 has also been shown to induce Aldh1a2 (aldehyde dehydrogenase 1a2) expression in murine macrophages recruited to the peritoneal cavity. However, the influence of IL-4 on DC Aldh1a2 induction in vivo has not yet been addressed. In this work, we found that DCs show enhanced aldehyde dehydrogenase enzyme activity in vivo, which led us to investigate the impact of the vitamin A metabolite all-trans retinoic acid (RA) on DC alternative activation and function. Antagonism of RA receptors reduced production of resistin-like molecule alpha by DCs responding to IL-4, while addition of exogenous RA enhanced production of this marker of alternative activation. Functionally, RA increased DC induction of CD4⁺ T-cell IL-10, while reducing CD4⁺ T-cell IL-4 and IL-13, revealing a previously unidentified role for RA in regulating the ability of alternatively activated DCs to influence T_h2 polarization.

Genomic analysis of the causative agents of coccidiosis in domestic chickens

Global production of chickens has trebled in the past two decades and they are now the most important source of dietary animal protein worldwide. Chickens are subject to many infectious diseases that reduce their performance and productivity. Coccidiosis, caused by apicomplexan protozoa of the genus Eimeria, is one of the most important poultry diseases. Understanding the biology of Eimeria parasites underpins development of new drugs and vaccines needed to improve global food security. We have produced annotated genome sequences of all seven species of Eimeria that infect domestic chickens, which reveal the full extent of previously described repeat-rich and repeat-poor regions and show that these parasites possess the most repeat-rich proteomes ever described. Furthermore, while no other apicomplexan has been found to possess retrotransposons, Eimeria is home to a family of chromoviruses. Analysis of Eimeria genes involved in basic biology and host-parasite interaction highlights adaptations to a relatively simple developmental life cycle and a complex array of co-expressed surface proteins involved in host cell binding.

Protein variation in blood-dwelling schistosome worms generated by differential splicing of micro-exon gene transcripts

Schistosoma mansoni is a well-adapted blood-dwelling parasitic helminth, persisting for decades in its human host despite being continually exposed to potential immune attack. Here, we describe in detail micro-exon genes (MEG) in S. mansoni, some present in multiple copies, which represent a novel molecular system for creating protein variation through the alternate splicing of short (< or =36 bp) symmetric exons organized in tandem. Analysis of three closely related copies of one MEG family allowed us to trace several evolutionary events and propose a mechanism for micro-exon generation and diversification. Microarray experiments show that the majority of MEGs are up-regulated in life cycle stages associated with establishment in the mammalian host after skin penetration. Sequencing of RT-PCR products allowed the description of several alternate splice forms of micro-exon genes, highlighting the potential use of these transcripts to generate a complex pool of protein variants. We obtained direct evidence for the existence of such pools by proteomic analysis of secretions from migrating schistosomula and mature eggs. Whole-mount in situ hybridization and immunolocalization showed that MEG transcripts and proteins were restricted to glands or epithelia exposed to the external environment. The ability of schistosomes to produce a complex pool of variant proteins aligns them with the other major groups of blood parasites, but using a completely different mechanism. We believe that our data open a new chapter in the study of immune evasion by schistosomes, and their ability to generate variant proteins could represent a significant obstacle to vaccine development.

Anti-schistosomal intervention targets identified by lifecycle transcriptomic analyses

Background

Novel methods to identify anthelmintic drug and vaccine targets are urgently needed, especially for those parasite species currently being controlled by singular, often limited strategies. A clearer understanding of the transcriptional components underpinning helminth development will enable identification of exploitable molecules essential for successful parasite/host interactions. Towards this end, we present a combinatorial, bioinformatics-led approach, employing both statistical and network analyses of transcriptomic data, for identifying new immunoprophylactic and therapeutic lead targets to combat schistosomiasis.

Methodology/Principal Findings

Utilisation of a Schistosoma mansoni oligonucleotide DNA microarray consisting of 37,632 elements enabled gene expression profiling from 15 distinct parasite lifecycle stages, spanning three unique ecological niches. Statistical approaches of data analysis revealed differential expression of 973 gene products that minimally describe the three major characteristics of schistosome development: asexual processes within intermediate snail hosts, sexual maturation within definitive vertebrate hosts and sexual dimorphism amongst adult male and female worms. Furthermore, we identified a group of 338 constitutively expressed schistosome gene products (including 41 transcripts sharing no sequence similarity outside the Platyhelminthes), which are likely to be essential for schistosome lifecycle progression. While highly informative, statistics-led bioinformatics mining of the transcriptional dataset has limitations, including the inability to identify higher order relationships between differentially expressed transcripts and lifecycle stages. Network analysis, coupled to Gene Ontology enrichment investigations, facilitated a re-examination of the dataset and identified 387 clusters (containing 12,132 gene products) displaying novel examples of developmentally regulated classes (including 294 schistosomula and/or adult transcripts with no known sequence similarity outside the Platyhelminthes), which were undetectable by the statistical comparisons.

Conclusions/Significance

Collectively, statistical and network-based exploratory analyses of transcriptomic datasets have led to a thorough characterisation of schistosome development. Information obtained from these experiments highlighted key transcriptional programs associated with lifecycle progression and identified numerous anti-schistosomal candidate molecules including G-protein coupled receptors, tetraspanins, Dyp-type peroxidases, fucosyltransferases, leishmanolysins and the netrin/netrin receptor complex.

Despite the implementation of focused and well-funded chemotherapeutic control initiatives over the last decade, schistosomiasis remains a significant cause of morbidity and mortality within countries of the developing world. There is, therefore, an urgent need for the rapid translation of genomic information into viable vaccines or new drug classes capable of eradicating the parasitic schistosome worms responsible for this neglected tropical disease. In our effort to identify potential targets for novel chemotherapeutic and immunoprophylactic interventions, we detail a combined bioinformatics approach, comprising exploratory statistical and network analyses, to thoroughly describe the transcriptional progression of Schistosoma mansoni across three environmental niches. Our results indicate that, although schistosomes are masters at host deception and survival, there are numerous exploitable candidate molecules displaying either differential or constitutive expression throughout the parasite's lifecycle. Importantly, some of these transcripts represent gene families not commonly found outside—or expanded within—the phylum Platyhelminthes, and thus represent priority targets. Many of the candidates identified herein will be subjected to ongoing and future hypothesis-led functional investigations. The completion of such specific examinations ultimately will contribute to the successful development of novel control strategies useful in the alleviation of schistosome-induced immunopathologies, morbidities and mortalities.

Comparative expression profiling of Leishmania: modulation in gene expression between species and in different host genetic backgrounds

Background

Genome sequencing of Leishmania species that give rise to a range of disease phenotypes in the host has revealed highly conserved gene content and synteny across the genus. Only a small number of genes are differentially distributed between the three species sequenced to date, L. major, L. infantum and L. braziliensis. It is not yet known how many of these genes are expressed in the disease-promoting intracellular amastigotes of these species or whether genes conserved between the species are differentially expressed in the host.

Methods/Principal Findings

We have used customised oligonucleotide microarrays to confirm that all of the differentially distributed genes identified by genome comparisons are expressed in intracellular amastigotes, with only a few of these subject to regulation at the RNA level. In the first large-scale study of gene expression in L. braziliensis, we show that only ∼9% of the genes analysed are regulated in their RNA expression during the L. braziliensis life cycle, a figure consistent with that observed in other Leishmania species. Comparing amastigote gene expression profiles between species confirms the proposal that Leishmania transcriptomes undergo little regulation but also identifies conserved genes that are regulated differently between species in the host. We have also investigated whether host immune competence influences parasite gene expression, by comparing RNA expression profiles in L. major amastigotes derived from either wild-type (BALB/c) or immunologically compromised (Rag2^−/− γ_c^−/−) mice. While parasite dissemination from the site of infection is enhanced in the Rag2^−/− γ_c^−/− genetic background, parasite RNA expression profiles are unperturbed.

Conclusion/Significance

These findings support the hypothesis that Leishmania amastigotes are pre-adapted for intracellular survival and undergo little dynamic modulation of gene expression at the RNA level. Species-specific parasite factors contributing to virulence and pathogenicity in the host may be limited to the products of a small number of differentially distributed genes or the differential regulation of conserved genes, either of which are subject to translational and/or post-translational controls.

The single-celled parasite Leishmania, transmitted by sand flies in more than 88 tropical and sub-tropical countries globally, infects man and other mammals, causing a spectrum of diseases called the leishmaniases. Over 12 million people are currently infected worldwide with 2 million new cases reported each year. The type of leishmaniasis that develops in the mammalian host is dependent on the species of infecting parasite and the immune response to infection (that can be influenced by host genetic variation). Our research is focused on identifying parasite factors that contribute to pathogenicity in the host and understanding how these might differ between parasite species that give rise to the different clinical forms of leishmaniasis. Molecules of this type might lead to new therapeutic tools in the longer term. In this paper, we report a comparative analysis of gene expression profiles in three Leishmania species that give rise to different types of disease, focusing on the intracellular stages that reside in mammalian macrophages. Our results show that there are only a small number of differences between these parasite species, with host genetics playing only a minor role in influencing the parasites' response to their intracellular habitat. These small changes may be significant, however, in determining the clinical outcome of infection.

Altered retinal microRNA expression profile in a mouse model of retinitis pigmentosa

MicroRNA expression profiling showed that the retina of mice carrying a rhodopsin mutation that leads to retinitis pigmentosa have notably different microRNA profiles from wildtype mice; further in silico analyses identified potential retinal targets for differentially regulated microRNAs.

Background

The role played by microRNAs (miRs) as common regulators in physiologic processes such as development and various disease states was recently highlighted. Retinitis pigmentosa (RP) linked to RHO (which encodes rhodopsin) is the most frequent form of inherited retinal degeneration that leads to blindness, for which there are no current therapies. Little is known about the cellular mechanisms that connect mutations within RHO to eventual photoreceptor cell death by apoptosis.

Results

Global miR expression profiling using miR microarray technology and quantitative real-time RT-PCR (qPCR) was performed in mouse retinas. RNA samples from retina of a mouse model of RP carrying a mutant Pro347Ser RHO transgene and from wild-type retina, brain and a whole-body representation (prepared by pooling total RNA from eight different mouse organs) exhibited notably different miR profiles. Expression of retina-specific and recently described retinal miRs was semi-quantitatively demonstrated in wild-type mouse retina. Alterations greater than twofold were found in the expression of nine miRs in Pro347Ser as compared with wild-type retina (P < 0.05). Expression of miR-1 and miR-133 decreased by more than 2.5-fold (P < 0.001), whereas expression of miR-96 and miR-183 increased by more than 3-fold (P < 0.001) in Pro347Ser retinas, as validated by qPCR. Potential retinal targets for these miRs were predicted in silico.

Conclusion

This is the first miR microarray study to focus on evaluating altered miR expression in retinal disease. Additionally, novel retinal preference for miR-376a and miR-691 was identified. The results obtained contribute toward elucidating the function of miRs in normal and diseased retina. Modulation of expression of retinal miRs may represent a future therapeutic strategy for retinopathies such as RP.

Altered patterns of gene expression underlying the enhanced immunogenicity of radiation-attenuated schistosomes

BACKGROUND

Schistosome cercariae only elicit high levels of protective immunity against a challenge infection if they are optimally attenuated by exposure to ionising radiation that truncates their migration in the lungs. However, the underlying molecular mechanisms responsible for the altered phenotype of the irradiated parasite that primes for protection have yet to be identified.

METHODOLOGY/PRINCIPAL FINDINGS

We have used a custom microarray comprising probes derived from lung-stage parasites to compare patterns of gene expression in schistosomula derived from normal and irradiated cercariae. These were transformed in vitro and cultured for four, seven, and ten days to correspond in development to the priming parasites, before RNA extraction. At these late times after the radiation insult, transcript suppression was the principal feature of the irradiated larvae. Individual gene analysis revealed that only seven were significantly down-regulated in the irradiated versus normal larvae at the three time-points; notably, four of the protein products are present in the tegument or associated with its membranes, perhaps indicating a perturbed function. Grouping of transcripts using Gene Ontology (GO) and subsequent Gene Set Enrichment Analysis (GSEA) proved more informative in teasing out subtle differences. Deficiencies in signalling pathways involving G-protein-coupled receptors suggest the parasite is less able to sense its environment. Reduction of cytoskeleton transcripts could indicate compromised structure which, coupled with a paucity of neuroreceptor transcripts, may mean the parasite is also unable to respond correctly to external stimuli.

CONCLUSIONS/SIGNIFICANCE

The transcriptional differences observed are concordant with the known extended transit of attenuated parasites through skin-draining lymph nodes and the lungs: prolonged priming of the immune system by the parasite, rather than over-expression of novel antigens, could thus explain the efficacy of the irradiated vaccine.

Microarray analysis identifies genes preferentially expressed in the lung schistosomulum of Schistosoma mansoni

The lung schistosomulum of Schistosoma mansoni is a validated target of protective immunity elicited in vaccinated mice. To identify genes expressed at this stage we constructed a microarray, representing 3088 contigs and singlets, with cDNA derived from in vitro cultured larvae and used it to screen RNA from seven life-cycle stages. Clustering of genes by expression profile across the life cycle revealed a number of membrane, membrane-associated and secreted proteins up-regulated at the lung stage, that may represent potential immune targets. Two promising secreted molecules have homology to antigens with vaccine and/or immunomodulatory potential in other helminths.

Integration of tools and resources for display and analysis of genomic data for protozoan parasites

Centralisation of tools for analysis of genomic data is paramount in ensuring that research is always carried out on the latest currently available data. As such, World Wide Web sites providing a range of online analyses and displays of data can play a crucial role in guaranteeing consistency of in silico work. In this respect, the protozoan parasite research community is served by several resources, either focussing on data and tools for one species or taking a broader view and providing tools for analysis of data from many species, thereby facilitating comparative studies. In this paper, we give a broad overview of the online resources available. We then focus on the GeneDB project, detailing the features and tools currently available through it. Finally, we discuss data curation and its importance in keeping genomic data 'relevant' to the research community.

Comparative genomics of trypanosomatid parasitic protozoa

A comparison of gene content and genome architecture of Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major, three related pathogens with different life cycles and disease pathology, revealed a conserved core proteome of about 6200 genes in large syntenic polycistronic gene clusters. Many species-specific genes, especially large surface antigen families, occur at nonsyntenic chromosome-internal and subtelomeric regions. Retroelements, structural RNAs, and gene family expansion are often associated with syntenic discontinuities that-along with gene divergence, acquisition and loss, and rearrangement within the syntenic regions-have shaped the genomes of each parasite. Contrary to recent reports, our analyses reveal no evidence that these species are descended from an ancestor that contained a photosynthetic endosymbiont.

The genome of the kinetoplastid parasite, Leishmania major

Leishmania species cause a spectrum of human diseases in tropical and subtropical regions of the world. We have sequenced the 36 chromosomes of the 32.8-megabase haploid genome of Leishmania major (Friedlin strain) and predict 911 RNA genes, 39 pseudogenes, and 8272 protein-coding genes, of which 36% can be ascribed a putative function. These include genes involved in host-pathogen interactions, such as proteolytic enzymes, and extensive machinery for synthesis of complex surface glycoconjugates. The organization of protein-coding genes into long, strand-specific, polycistronic clusters and lack of general transcription factors in the L. major, Trypanosoma brucei, and Trypanosoma cruzi (Tritryp) genomes suggest that the mechanisms regulating RNA polymerase II-directed transcription are distinct from those operating in other eukaryotes, although the trypanosomatids appear capable of chromatin remodeling. Abundant RNA-binding proteins are encoded in the Tritryp genomes, consistent with active posttranscriptional regulation of gene expression.

Expression profiling of the Leishmania life cycle: cDNA arrays identify developmentally regulated genes present but not annotated in the genome

As genomic sequencing of Leishmania nears completion, functional analyses that provide a global genetic perspective on biological processes are important. Despite polycistronic transcription, RNA transcript abundance can be measured using microarrays. To provide a resource to evaluate cDNA arrays, we undertook 5' expressed sequence tag analysis of 2183 full-length randomly selected cDNAs from Leishmania major promastigote (days 3, 7, 10 of culture in vitro), and lesion-derived amastigote libraries. PCR-amplified inserts from 1830 of these cDNA representing 1001 unique genes were spotted onto microarrays, and compared internally with PCR-amplified open reading frames (ORFs) from 904 genes representing 842 unique genes annotated in the L. major genome. Microarrays were screened with RNA from procyclic, metacyclic and amastigote populations of L. major. Redundant clones on the array gave highly reproducible results, providing confidence in identification of stage-specific gene expression. Four hundred and thirty unique (i.e. non-redundant) stage-specific genes were identified. A higher percentage of stage-specific gene expression was observed in amastigotes ( approximately 35%) compared to metacyclics ( approximately 12%) for both cDNAs and ORFs, but cDNAs provided a richer source of regulated genes than currently annotated ORFs from the Leishmania genome. In mapping cDNAs onto the Leishmania genome, we noted that approximately 42% aligned to regions not recognised as genes using current predictive annotation tools. These genes are highly represented in our stage-specific genes, and therefore represent important drug targets and vaccine candidates. Careful annotation of cDNAs onto the Leishmania genome will be important before producing the next generation of oligonucleotide arrays based on annotated genes of the genomic sequencing project.

GeneDB: a resource for prokaryotic and eukaryotic organisms

GeneDB (http://www.genedb.org/) is a genome database for prokaryotic and eukaryotic organisms. The resource provides a portal through which data generated by the Pathogen Sequencing Unit at the Wellcome Trust Sanger Institute and other collaborating sequencing centres can be made publicly available. It combines data from finished and ongoing genome and expressed sequence tag (EST) projects with curated annotation, that can be searched, sorted and downloaded, using a single web based resource. The current release stores 11 datasets of which six are curated and maintained by biologists, who review and incorporate information from the scientific literature, public databases and the respective research communities.

Leishmania major chromosome 3 contains two long convergent polycistronic gene clusters separated by a tRNA gene

Leishmania parasites (order Kinetoplastida, family Trypanosomatidae) cause a spectrum of human diseases ranging from asymptomatic to lethal. The approximately 33.6 Mb genome is distributed among 36 chromosome pairs that range in size from approximately 0.3 to 2.8 Mb. The complete nucleotide sequence of Leishmania major Friedlin chromosome 1 revealed 79 protein-coding genes organized into two divergent polycistronic gene clusters with the mRNAs transcribed towards the telomeres. We report here the complete nucleotide sequence of chromosome 3 (384 518 bp) and an analysis revealing 95 putative protein-coding ORFs. The ORFs are primarily organized into two large convergent polycistronic gene clusters (i.e. transcribed from the telomeres). In addition, a single gene at the left end is transcribed divergently towards the telomere, and a tRNA gene separates the two convergent gene clusters. Numerous genes have been identified, including those for metabolic enzymes, kinases, transporters, ribosomal proteins, spliceosome components, helicases, an RNA-binding protein and a DNA primase subunit.

Boudicca, a retrovirus-like long terminal repeat retrotransposon from the genome of the human blood fluke Schistosoma mansoni

The genome of Schistosoma mansoni contains a proviral form of a retrovirus-like long terminal repeat (LTR) retrotransposon, designated BOUDICCA: Sequence and structural characterization of the new mobile genetic element, which was found in bacterial artificial chromosomes prepared from S. mansoni genomic DNA, revealed the presence of three putative open reading frames (ORFs) bounded by direct LTRs of 328 bp in length. ORF1 encoded a retrovirus-like major homology region and a Cys/His box motif, also present in Gag polyproteins of related retrotransposons and retroviruses. ORF2 encoded enzymatic domains and motifs characteristic of a retrovirus-like polyprotein, including aspartic protease, reverse transcriptase, RNase H, and integrase, in that order, a domain order similar to that of the gypsy/Ty3 retrotransposons. An additional ORF at the 3' end of the retrotransposon may encode an envelope protein. Phylogenetic comparison based on the reverse transcriptase domain of ORF2 confirmed that Boudicca was a gypsy-like retrotransposon and showed that it was most closely related to CsRn1 from the Oriental liver fluke Clonorchis sinensis and to kabuki from Bombyx mori. Bioinformatics approaches together with Southern hybridization analysis of genomic DNA of S. mansoni and the screening of a bacterial artificial chromosome library representing approximately 8-fold coverage of the S. mansoni genome revealed that numerous copies of Boudicca were interspersed throughout the schistosome genome. By reverse transcription-PCR, mRNA transcripts were detected in the sporocyst, cercaria, and adult developmental stages of S. mansoni, indicating that Boudicca is actively transcribed in this trematode.

Leishmania RAB7: characterisation of terminal endocytic stages in an intracellular parasite

Leishmania species are intracellular parasites that inhabit a parasitophorous vacuole (PV) within host macrophages and engage with the host endo-membrane network to avoid clearance from the cell. Intracellular Leishmania amastigotes exhibit a high degree of proteolytic/lysosomal activity that may assist degradation of MHC class II molecules and subsequent interruption of antigen presentation. As an aid to further analysis of the endosomal/lysosomal events that could facilitate this process, we have characterised a Leishmania homologue of the late endosomal marker, Rab7, thought to be involved in the terminal steps of endocytosis and lysosomal delivery. The Leishmania major Rab7 (LmRAB7) protein is expressed throughout the life-cycle, shows 73 and 64% identity to Trypanosoma cruzi and Trypanosoma brucei Rab7s (TcRAB7 and TbRAB7), respectively, and includes a kinetoplastid-specific insertion. The recombinant protein binds GTP and polyclonal antibodies raised against this antigen recognise structures in the region of the cell between the nucleus and kinetoplast. By immunoelectron microscopy of axenic amastigotes, Leishmania mexicana Rab7 (LmexRAB7) is found juxtaposed to and overlapping membrane structures labelled for the megasomal marker, cysteine proteinase B, confirming a late-endosomal/lysosomal localisation.

Molecular cloning and characterization of two new isoforms of the protein kinase A catalytic subunit from the human parasite Leishmania

Leishmania are protozoan parasites that cause extensive morbidity and mortality in humans. Genes for two new isoforms of the protein kinase A catalytic subunit (PKAC) in Leishmania, Lmpkac2a and Lmpkac2b, were cloned and characterized. The predicted open reading frames for these isoforms are 93.4% identical over 338 amino acids (aa). The conserved PK catalytic cores (subdomains I-XI) are identical, while the carboxy-terminal extensions differ by only two aa. However, LmPKAC2 shares only 62% identity over the 255 aa catalytic core region with the previously described LmPKAC1 (c-lpk2). Unlike LmPKAC1, the location of the FXXF motif at the carboxy-terminus is conserved in both LmPKAC2 isoforms; however, the aa sequence, LXXF, in isoform-2a is unusual. The leishmanial isoforms can be distinguished by their NH(2)-terminal extensions, which show minimal similarity at the primary sequence level. Structural analysis of the three enzymes based on the crystal structure of mammalian PKAs predicts that both LmPKAC2 isoforms, unlike LmPKAC1, have identical alpha-helix structures in the NH(2)-terminal extension. Lmpkac2 genes are located on chromosome 35 just downstream from the leishmanial prp8 gene. This genomic organization is conserved in two species of Leishmania and Crithidia fasciculata and allowed for the partial analysis of Cfpkac2a. Phylogenetic analysis groups the two LmPKAC2 isoforms together and separately from LmPKAC1, which is more similar to the Euglena gracilis PKAC, EPK2. These findings provide the basis for additional studies on the role of the PKA family in parasite differentiation and virulence.

The Leishmania genome project: new insights into gene organization and function

The sequencing of Leishmania major Friedlin chromosome 1 (Chr1), Chr3, and Chr4 has been completed. and several other chromosomes are well underway. The complete genome sequence should be available by 2003. Over 1,000 full-length new genes have been identified, with the majority (approximately 75%) having unknown function. Many of these may be Leishmania (or kinetoplastid) specific. Most interestingly, the genes are organized into large (> 100-500 kb) polycistronic clusters of adjacent genes on the same DNA strand. Chr1 contains two such clusters organized in a "divergent" manner, i.e., the mRNAs for the two sets of genes are both transcribed towards the telomeres. Nuclear run-on analysis suggests that transcription is initiated in both directions within the "divergent" region. Chr3 and Chr4 contain two "convergent" clusters, with a single "divergent" gene at one telomere of Chr3. Sequence analysis of several genes from the LD1 region of Chr35 indicates a high degree of sequence conservation between L. major and L. donovani/L. infantum within protein-coding open reading frames (ORFs), with a lower degree of conservation within the non-coding regions. Immunization of mice with recombinant antigen from two of these genes, BTI (formerly ORFG) and ORFF, results in significant reduction in parasite burden following Leishmania challenge. Recombinant ORFF antigen shows promise as a serodiagnostic. We have also developed a tetracycline-regulated promoter system, which allows us to modulate gene expression in Leishmania.

Secondary DNA structure analysis of the coding strand switch regions of five Leishmania major Friedlin chromosomes

As part of the EULEISH international genome project, a region of 74,674 nucleotides from chromosome 21 of Leishmania major Friedlin was subcloned and sequenced; and 31 new coding sequences were predicted. Of particular interest was a unique coding strand switching region covering 1.6 kb of DNA; and this was subjected to further investigation. Bioinformatic analysis of this region revealed an unusually high AT composition, a lack of putative hairpins and a strong curvature of the DNA in agreement with the structural characteristics of similar regions of other Leishmania chromosomes. These observations and a comparison with the secondary DNA structure of four other Leishmania chromosomes and chromosomes of different organisms could suggest a functional role of this region in transcription and mitotic division.

Identification and cloning of Lmairk, a member of the Aurora/Ipl1p protein kinase family, from the human protozoan parasite Leishmania

Lmairk, a gene encoding a member of the Aurora/Ipl1p family of protein kinases (AIRK), was cloned from the protozoan parasite Leishmania major. Aurora kinases are key enzymes involved in the regulation of normal chromosome segregation during mitosis and cytokenesis of eukaryotic cells. This single-copy gene located on L. major chromosome 28 encodes a 301 amino acid polypeptide. All 11 conserved eukaryotic protein kinase catalytic subdomains are present and the proposed AIRK signature sequence was identified in the activation loop between subdomains VII and VIII. Lmairk is expressed, as an approximately 2.4 kb message, in at least three different species of Leishmania. This report represents the first identification of an AIRK from the trypanosomatid family of early divergent eukaryotes.

The Leishmania genome comes of Age

The Leishmania Genome Network (LGN) was born in Rio de Janeiro, Brazil in 1994. In the short period that has elapsed since then, the LGN has focused solely on the acquisition of the resources, and hence data, that have enabled a rational approach to genomic sequencing of the reference strain, Leishmania major Friedlin. This has now been achieved. In this review, Alasdair Ivens and Jennie Blackwell, secretary and chairman of the LGN, respectively, re-examine the approaches that were adopted, comment on some of the interesting data that have been obtained and introduce some genome-wide approaches that will facilitate functional studies of the parasite.

Genomics and the biology of parasites

Despite the advances of modern medicine, the threat of chronic illness, disfigurement, or death that can result from parasitic infection still affects the majority of the world population, retarding economic development. For most parasitic diseases, current therapeutics often leave much to be desired in terms of administration regime, toxicity, or effectiveness and potential vaccines are a long way from market. Our best prospects for identifying new targets for drug, vaccine, and diagnostics development and for dissecting the biological basis of drug resistance, antigenic diversity, infectivity and pathology lie in parasite genome analysis, and international mapping and gene discovery initiatives are under way for a variety of protozoan and helminth parasites. These are far from ideal experimental organisms, and the influence of biological and genomic characteristics on experimental approaches is discussed, progress is reviewed and future prospects are examined.

Differentially expressed Leishmania major gp63 genes encode cell surface leishmanolysin with distinct signals for glycosylphosphatidylinositol attachment

The Leishmania cell surface metalloproteinase, leishmanolysin or GP63, is expressed in all stages of Leishmania major. Initial studies reported that in L. major the gp63 genes were arranged as five homologous, tandemly repeated genes (gp63 genes 1-5) and a sixth, less conserved gp63 gene located 8 kb downstream of gp63 gene 5. This study compared the sequences of L. major gp63 gene 1 and gp63 gene 6 and identified a seventh L. major gp63 gene located downstream from gp63 gene 6. The L. major gp63 genes exhibited stage-specific differences in their expression: gp63 genes 1-5 were expressed in promastigotes only, gp63 gene 6 was expressed in promastigotes and amastigotes, while gp63 gene 7 was expressed predominantly in stationary phase promastigotes and in amastigotes. Analysis of the predicted protein sequence of gp63 gene 6 (GP63-6) and gp63 gene 1 (GP63-1) showed that these two proteins were homologous in terms of overall predicted domain structure. L. major GP63-1 has been reported to contain a glycosylphosphatidylinositol (GPI) membrane anchor while sequence analysis predicted that GP63-6 contained a different hydrophobic C-terminus that may act as a transmembrane region. Transfection studies using L. major gp63 gene 1 and gp63 gene 6 expressed in L. donovani promastigotes showed that GP63-6 was expressed at the cell surface and that the distinct GP63-6 C-terminus was capable of mediating GPI anchor attachment.

A physical map of the Leishmania major Friedlin genome

An extensive physical map of the Leishmania major Friedlin genome has been assembled by the combination of fingerprint analysis of a shuttle vector cosmid library and probe hybridization. The integrated data obtained for 9004 fingerprinted clones and 974 probes have placed 91.2% of the 33.58-Mb genome into contigs representing each of the 36 chromosomes. This first-generation map has already provided a suitable framework for both high-throughput DNA sequencing and functional studies of the L. major parasite.

Parasite genome analysis. A global map of the Leishmania major genome: prelude to genomic sequencing

In 1994, the World Health Organization (TDR) launched a new strategic initiative in parasite genome analysis, establishing international genome networks for filariae, Schistosoma, Leishmania, Trypanosoma brucei and T. cruzi. For Leishmania, a number of different but complementary approaches have been adopted by members of the Leishmania Genome Network. Our laboratory has been using cosmid clone fingerprinting to produce a physical map of the genome. Progress towards the completion of an integrated physical and biological map of L. major, and the preparations for genomic sequencing, are described.

Unravelling the Leishmania genome

The past few years have been significant advances in our understanding of eukaryotic genomes. In the field of parasitology, this is best exemplified by the application of genome mapping techniques to the study of genome structure and function in the protozoan parasite, Leishmania. Although much is known about the organism and the diseases it causes, molecular genetics has only recently begun to play a major part in elucidating some of the unusual characteristics of this interesting parasite. Mapping of the small (35 Mb) genome and determination of the functional role of genes by the application of in vitro homologous gene targeting techniques are revealing novel avenues for the development of prophylactic measures.

An integrated physical map of 210 markers assigned to the short arm of human chromosome 11

Using a panel of patient cell lines with chromosomal breakpoints, we constructed a physical map for the short arm of human chromosome 11. We focused on 11p15, a chromosome band harboring at least 25 known genes and associated with the Beckwith-Wiedemann syndrome, several childhood tumors, and genomic imprinting. This underlines the need for a physical map for this region. We divided the short arm of chromosome 11 into 18 breakpoint regions, and a large series of new and previously described genes and markers was mapped within these intervals using fluorescence in situ hybridization. Cosmid fingerprint analysis showed that 19 of these markers were included in cosmid contigs. A detailed 10-Mb pulsed-field physical map of the region 11p15.3-pter was constructed. These three different approaches enabled the high-resolution mapping of 210 markers, including 22 known genes.

The generation of ordered sets of cosmid DNA clones from human chromosome region 11p

We describe progress in a continuing project aimed at the generation of an overlapping cosmid DNA clone map of the short arm of human chromosome 11. The automated procedures used to prepare DNA samples and the computerized data collection and recording systems are described. We also demonstrate the use of the clones as reagents for the rapid isolation of genomic DNAs containing smaller probed regions. We have isolated approximately 4700 human cosmid DNA clones from mouse/human hybrid cell lines that contain predominantly human chromosomal region 11p. Of the DNA in the cell lines, 60% is derived from this chromosomal region, and the remaining 40% is derived from regions of chromosomes 3, 19, and 20. A total of 4159 clones have been fingerprinted to identify potential overlaps, and we have developed 535 sets ("contigs"). Using random modeling, it is estimated that 65% of 11p must be contained in the analyzed cosmids. The database of clones has been used to identify single or overlapping clones from noncosmid DNA probes. Examples are presented. It is proposed that cosmid reference filters be distributed to requesting laboratories.