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 (<= 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.

The genome of the blood fluke Schistosoma mansoni

Schistosoma mansoni is responsible for the neglected tropical disease schistosomiasis that affects 210 million people in 76 countries. Here we present analysis of the 363 megabase nuclear genome of the blood fluke. It encodes at least 11,809 genes, with an unusual intron size distribution, and new families of micro-exon genes that undergo frequent alternative splicing. As the first sequenced flatworm, and a representative of the Lophotrochozoa, it offers insights into early events in the evolution of the animals, including the development of a body pattern with bilateral symmetry, and the development of tissues into organs. Our analysis has been informed by the need to find new drug targets. The deficits in lipid metabolism that make schistosomes dependent on the host are revealed, and the identification of membrane receptors, ion channels and more than 300 proteases provide new insights into the biology of the life cycle and new targets. Bioinformatics approaches have identified metabolic chokepoints, and a chemogenomic screen has pinpointed schistosome proteins for which existing drugs may be active. The information generated provides an invaluable resource for the research community to develop much needed new control tools for the treatment and eradication of this important and neglected disease.

Exploring the Schistosoma mansoni adult male transcriptome using RNA-seq

Schistosoma mansoni is one of the agents of schistosomiasis, a chronic and debilitating disease. Here we, present a transcriptome-wide characterization of adult S. mansoni males by high-throughput RNA-sequencing. We obtained 1,620,432 high-quality ESTs from a directional strand-specific cDNA library, resulting in a 26% higher coverage of genome bases than that of the public ESTs available at NCBI. With a 15 x-deep coverage of transcribed genomic regions, our data were able to (i) confirm for the first time 990 predictions without previous evidence of transcription; (ii) correct gene predictions; (iii) discover 989 and 1196 RNA-seq contigs that map to intergenic and intronic genomic regions, respectively, where no gene had been predicted before. These contigs could represent new protein-coding genes or non-coding RNAs (ncRNAs). Interestingly, we identified 11 novel Micro-exon genes (MEGs). These data reveal new features of the S. mansoni transcriptional landscape and significantly advance our understanding of the parasite transcriptome. (c) 2011 Elsevier Inc. All rights reserved.

Transcriptional profiling of the oesophageal gland region of male worms of Schistosoma mansoni

The intestinal tract of schistosomes opens at the mouth and leads into the foregut or oesophageal region that is lined with syncytium continuous with the apical cytoplasm of the tegument. The oesophagus is surrounded by a specialised gland, the oesophageal gland. This gland releases materials into the lumen of the oesophagus and the region is thought to initiate the lysis of erythrocytes and neutralisation of immune effectors of the host. The oesophageal region is present in the early invasive schistosomulum, a stage potentially targetable by anti-schistosome vaccines. We used a 44k oligonucleotide microarray to identify highly up-regulated genes in microdissected frozen sections of the oesophageal gland of male worms of S. mansoni. We show that 122 genes were up-regulated 2-fold or higher in the oesophageal gland compared with a whole male worm tissue control. The enriched genes included several associated with lipid metabolism and transmembrane transport as well as some micro-exon genes. Since the oesophageal gland is important in the initiation of digestion and the fact that it develops early after invasion of the mammalian host, further study of selected highly up-regulated functionally important genes in this tissue may reveal new anti-schistosome intervention targets for schistosomiasis control.

Dopamine receptors in a songbird brain.

Dopamine is a key neuromodulatory transmitter in the brain. It acts through dopamine receptors to affect changes in neural activity, gene expression, and behavior. In songbirds, dopamine is released into the striatal song nucleus Area X, and the levels depend on social contexts of undirected and directed singing. This differential release is associated with differential expression of activity-dependent genes, such as egr1 (avian zenk), which in mammalian brain are modulated by dopamine receptors. Here we cloned from zebra finch brain cDNAs of all avian dopamine receptors: the D1 (D1A, D1B, D1D) and D2 (D2, D3, D4) families. Comparative sequence analyses of predicted proteins revealed expected phylogenetic relationships, in which the D1 family exists as single exon and the D2 family exists as spliced exon genes. In both zebra finch and chicken, the D1A, D1B, and D2 receptors were highly expressed in the striatum, the D1D and D3 throughout the pallium and within the mesopallium, respectively, and the D4 mainly in the cerebellum. Furthermore, within the zebra finch, all receptors, except for D4, showed differential expression in song nuclei relative to the surrounding regions and developmentally regulated expression that decreased for most receptors during the sensory acquisition and sensorimotor phases of song learning. Within Area X, half of the cells expressed both D1A and D2 receptors, and a higher proportion of the D1A-only-containing neurons expressed egr1 during undirected but not during directed singing. Our findings are consistent with hypotheses that dopamine receptors may be involved in song development and social context-dependent behaviors.

Alternative splicing of DNA polymerase beta in vertebrates

Alternative splicing (AS) is the predominant mechanism responsible for increasing eukaryotic transcriptome and proteome complexity. In this phenomenon, numerous mRNA transcripts are produced from a single pre-mRNA sequence. AS is reported to occur in 95% of human multi-exon genes; one specific gene that undergoes AS is DNA polymerase beta (POLB). POLB is the main DNA repair gene which performs short patch base excision repair (BER). In primate untransformed primary fibroblast cell lines, it was determined that the splice variant (SV) frequency of POLB correlates positively with species lifespan. To date, AS patterns of POLB have only been examined in mammals primarily through the use of cell lines. However, little attention has been devoted to investigating if such a relationship exists in non-mammals and whether cell lines reflect what is observed in vertebrate tissues. This idea was explored through cloning and characterization of 1,214 POLB transcripts from four non-mammalian species (Gallus gallus domesticus, Larus glaucescens, Xenopus laevis, and Pogona vitticeps) and two mammalian species (Sylvilagus floridanus and Homo sapiens) in two tissue types, liver and brain. POLB SV frequency occurred at low frequencies, < 3.2%, in non-mammalian tissues relative to mammalian (>20%). The highest POLB SV frequency was found in H. sapiens liver and brain tissues, occurring at 65.4% and 91.7%, respectively. Tissue specific AS of POLB was observed in L. glaucescens, P. vitticeps, and H. sapiens, but not G. gallus domesticus, X. laevis and S. floridanus.The AS patterns of a second gene, transient receptor potential cation channel subfamily V member 1 (TRPV1), were compared to those of POLB in liver and brain tissues of G. gallus domesticus, X. laevis and H. sapiens. This comparison was performed to investigate if any changes (either increase or decrease) observed in the AS of POLB were gene specific or if they were tissue specific, in which case similar changes in AS would be seen in POLB and TRPV1. Analysis did not reveal an increase or decrease in both the AS of POLB and TRPV1 in either the liver or brain tissues of G. gallus domesticus and H. sapiens. This result suggested that the AS patterns of POLB were not influenced by tissue specific rates of AS. Interestingly, an increase in the AS of both genes was only observed in X. laevis brain tissue. This result suggests that AS in general may be increased in the X. laevis brain as compared to liver tissue. No positive correlation between POLB SV frequency and species lifespan was found in non-mammalian tissues. The AS patterns of POLB in human primary untransformed fibroblast cell lines were representative of those seen in human liver tissue but not in brain tissue. Altogether, the AS patterns of POLB from vertebrate tissues and primate cell lines revealed a positive correlation between POLB SV frequency and lifespan in mammals, but not in non-mammals. It appears that this positive correlation does not exist in vertebrate species as a whole.

Functional consequences of new exon acquisition in mammalian chromodomain Y-like (CDYL) genes

The origin of new exons is an important mechanism for proteome diversity. Here, we report the recurrent origination of new exons in mammalian chromodomain Y-like (CDYL) genes and the functional consequences associated with the acquisition of the new exons

Transcriptional silencing of nonsense codon-containing immunoglobulin micro genes requires translation of its mRNA

Eukaryotes have evolved quality control mechanisms that prevent the expression of genes in which the protein coding potential is crippled by the presence of a premature translation-termination codon (PTC). In addition to nonsense-mediated mRNA decay (NMD), a well documented posttranscriptional consequence of the presence of a PTC in an mRNA, we recently reported the transcriptional silencing of PTC-containing immunoglobulin (Ig) mu and gamma minigenes when they are stably integrated into the genome of HeLa cells. Here we demonstrate that this transcriptional silencing of PTC-containing Ig-mu constructs requires active translation of the cognate mRNA, as it is not observed under conditions where translation of the PTC-containing mRNA is inhibited through an iron-responsive element in the 5'-untranslated region. Furthermore, RNA interference-mediated depletion of the essential NMD factor Upf1 not only abolishes NMD but also reduces the extent of nonsense-mediated transcriptional gene silencing (NMTGS). Collectively, our data indicate that NMTGS and NMD are linked, relying on the same mechanism for PTC recognition, and that the NMTGS pathway branches from the NMD pathway at a step after Upf1 function.

A RAB3GAP1 SINE Insertion in Alaskan Huskies with Polyneuropathy, Ocular Abnormalities and Neuronal Vacuolation (POANV) Resembling Human Warburg Micro Syndrome 1 (WARBM1)

We observed a hereditary phenotype in Alaskan Huskies, which was characterized by polyneuropathy with ocular abnormalities and neuronal vacuolation (POANV). The affected dogs developed a progressive severe ataxia, which led to euthanasia between 8 and 16 months of age. The pedigrees were consistent with a monogenic autosomal recessive inheritance. We localized the causative genetic defect to a 4 Mb interval on chromosome 19 by a combined linkage and homozygosity mapping approach. Whole genome sequencing of one affected dog, an obligate carrier and an unrelated control revealed a 218 bp SINE insertion into exon 7 of the RAB3GAP1 gene. The SINE insertion was perfectly associated with the disease phenotype in a cohort of 43 Alaskan Huskies and it was absent from 541 control dogs of diverse other breeds. The SINE insertion induced aberrant splicing and led to a transcript with a greatly altered exon 7. RAB3GAP1 loss-of-function variants in humans cause Warburg Micro Syndrome 1 (WARBM1), which is characterized by additional developmental defects compared to canine POANV, whereas Rab3gap1 deficient mice have a much milder phenotype than either humans or dogs. Thus the RAB3GAP1 mutant Alaskan Huskies provide an interesting intermediate phenotype that may help to better understand the function of RAB3GAP1 in development. Furthermore, the identification of the presumed causative genetic variant will enable genetic testing to avoid the non-intentional breeding of affected dogs.

Identification of micro-RNAs targeting genes of PI3K/AKT pathway in ovarian cancer

Ovarian cancer is the leading cause of cancer-related death for females due to lack of specific early detection method. It is of great interest to find molecular-based biomarkers which are sensitive and specific to ovarian cancer for early diagnosis, prognosis and therapeutics. miRNAs have been proposed to be potential biomarkers that could be used in cancer prevention and therapeutics. The current study analyzed the miRNA and mRNA expression data extracted from the Cancer Genome Atlas (TCGA) database. Using simple linear regression and multiple regression models, we found 71 miRNA-mRNA pairs which were negatively associated between 56 miRNAs and 24 genes of PI3K/AKT pathway. Among these miRNA and mRNA target pairs, 9 of them were in agreement with the predictions from the most commonly used target prediction programs including miRGen, miRDB, miRTarbase and miR2Disease. These shared miRNA-mRNA pairs were considered to be the most potential genes that were involved in ovarian cancer. Furthermore, 4 of the 9 target genes encode cell cycle or apoptosis related proteins including Cyclin D1, p21, FOXO1 and Bcl2, suggesting that their regulator miRNAs including miR-16, miR-96 and miR-21 most likely played important roles in promoting tumor growth through dysregulated cell cycle or apoptosis. miR-96 was also found to directly target IRS-1. In addition, the results showed that miR-17 and miR-9 may be involved in ovarian cancer through targeting JAK1. This study might provide evidence for using miRNA or miRNA profile as biomarker.^

Intron phase correlations and the evolution of the intron/exon structure of genes.

Two issues in the evolution of the intron/exon structure of genes are the role of exon shuffling and the origin of introns. Using a large data base of eukaryotic intron-containing genes, we have found that there are correlations between intron phases leading to an excess of symmetric exons and symmetric exon sets. We interpret these excesses as manifestations of exon shuffling and make a conservative estimate that at least 19% of the exons in the data base were involved in exon shuffling, suggesting an important role for exon shuffling in evolution. Furthermore, these excesses of symmetric exons appear also in those regions of eukaryotic genes that are homologous to prokaryotic genes: the ancient conserved regions. This last fact cannot be explained in terms of the insertional theory of introns but rather supports the concept that some of the introns were ancient, the exon theory of genes.

The 3'-terminal exon of the family of steroid and phenol sulfotransferase genes is spliced at the N-terminal glycine of the universally conserved GXXGXXK motif that forms the sulfonate donor binding site.

The guinea pig estrogen sulfotransferase gene has been cloned and compared to three other cloned steroid and phenol sulfotransferase genes (human estrogen sulfotransferase, human phenol sulfotransferase, and guinea pig 3 alpha-hydroxysteroid sulfotransferase). The four sulfotransferase genes demonstrate a common outstanding feature: the splice sites for their 3'-terminal exons are identically located. That is, the 3'-terminal exon splice sites involve a glycine that constitutes the N-terminal glycine of an invariably conserved GXXGXXK motif present in all steroid and phenol sulfotransferases for which primary structures are known. This consistency strongly suggests that all steroid and phenol sulfotransferase genes will be similarly spliced. The GXXGXXK motif forms the active binding site for the universal sulfonate donor 3'-phosphoadenosine 5'-phosphosulfate. Amino acid sequence alignment of 19 cloned steroid and phenol sulfotransferases starting with the GXXGXXK motif indicates that the 3'-terminal exon for each steroid and phenol sulfotransferase gene encodes a similarly sized C-terminal fragment of the protein. Interestingly, on further analysis of the alignment, three distinct amino acid sequence patterns emerge. The presence of the conserved functional GXXGXXK motif suggests that the protein domains encoded by steroid and phenol sulfotransferase 3'-terminal exons have evolved from a common ancestor. Furthermore, it is hypothesized that during the course of evolution, the 3'-terminal exon further diverged into at least three sulfotransferase subdivisions: a phenol or aryl group, an estrogen or phenolic steroid group, and a neutral steroid group.

Evidence against the exon theory of genes derived from the triose-phosphate isomerase gene.

The exon theory of genes proposes that the introns of protein-encoding nuclear genes are remnants of the DNA spacers between ancient minigenes. The discovery of an intron at a predicted position in the triose-phosphate isomerase (EC 5.3.1.1) gene of Culex mosquitoes has been hailed as an evidential pillar of the theory. We have found that that intron is also present in Aedes mosquitoes, which are closely related to Culex, but not in the phylogenetically more distant Anopheles, nor in the fly Calliphora vicina, nor in the moth Spodoptera littoralis. The presence of this intron in Culex and Aedes is parsimoniously explained as the result of an insertion in a recent common ancestor of these two species rather than as the remnant of an ancient intron. The absence of the intron in 19 species of very diverse organisms requires at least 10 independent evolutionary losses in order to be consistent with the exon theory.

A re-examination of the Ghrelin and Ghrelin receptor genes

The last few years have seen dramatic advances in genomics, including the discovery of a large number of non-coding and antisense transcripts. This has revolutionised our understanding of multifaceted transcript structures found within gene loci and their roles in the regulation of development, neurogenesis and other complex processes. The recent and continuing surge of knowledge has prompted researchers to reassess and further dissect gene loci. The ghrelin gene (GHRL) gives rise to preproghrelin, which in turn produces ghrelin, a 28 amino acid peptide hormone that acts via the ghrelin receptor (growth hormone secretagogue receptor/GHSR 1a). Ghrelin has many important physiological and pathophysiological roles, including the stimulation of growth hormone (GH) release, appetite regulation, and cancer development. A truncated receptor splice variant, GHSR 1b, does not bind ghrelin, but dimerises with GHSR 1a, and may act as a dominant negative receptor. The gene products of ghrelin and its receptor are frequently overexpressed in human cancer While it is well known that the ghrelin axis (ghrelin and its receptor) plays a range of important functional roles, little is known about the molecular structure and regulation of the ghrelin gene (GHRL) and ghrelin receptor gene (GHSR). This thesis reports the re-annotation of the ghrelin gene, discovery of alternative 5’ exons and transcription start sites, as well as the description of a number of novel splice variants, including isoforms with a putative signal peptide. We also describe the discovery and characterisation of a ghrelin antisense gene (GHRLOS), and the discovery and expression of a ghrelin receptor (growth hormone secretagogue receptor/GHSR) antisense gene (GHSR-OS). We have identified numerous ghrelin-derived transcripts, including variants with extended 5' untranslated regions and putative secreted obestatin and C-ghrelin transcripts. These transcripts initiate from novel first exons, exon -1, exon 0 and a 5' extended 1, with multiple transcription start sites. We used comparative genomics to identify, and RT-PCR to experimentally verify, that the proximal exon 0 and 5' extended exon 1 are transcribed in the mouse ghrelin gene, which suggests the mouse and human proximal first exon architecture is conserved. We have identified numerous novel antisense transcripts in the ghrelin locus. A candidate non-coding endogenous natural antisense gene (GHRLOS) was cloned and demonstrates very low expression levels in the stomach and high levels in the thymus, testis and brain - all major tissues of non-coding RNA expression. Next, we examined if transcription occurs in the antisense orientation to the ghrelin receptor gene, GHSR. A novel gene (GHSR-OS) on the opposite strand of intron 1 of the GHSR gene was identified and characterised using strand-specific RT-PCR and rapid amplification of cDNA ends (RACE). GHSR-OS is differentially expressed and a candidate non-coding RNA gene. In summary, this study has characterised the ghrelin and ghrelin receptor loci and demonstrated natural antisense transcripts to ghrelin and its receptor. Our preliminary work shows that the ghrelin axis generates a broad and complex transcriptional repertoire. This study provides the basis for detailed functional studies of the the ghrelin and GHSR loci and future studies will be needed to further unravel the function, diagnostic and therapeutic potential of the ghrelin axis.