Two forms of Drosophila melanogaster Gs alpha are produced by alternate splicing involving an unusual splice site

SpliceFinder: ab initio prediction of splice sites using convolutional neural network

Background

Identifying splice sites is a necessary step to analyze the location and structure of genes. Two dinucleotides, GT and AG, are highly frequent on splice sites, and many other patterns are also on splice sites with important biological functions. Meanwhile, the dinucleotides occur frequently at the sequences without splice sites, which makes the prediction prone to generate false positives. Most existing tools select all the sequences with the two dimers and then focus on distinguishing the true splice sites from those pseudo ones. Such an approach will lead to a decrease in false positives; however, it will result in non-canonical splice sites missing.

Result

We have designed SpliceFinder based on convolutional neural network (CNN) to predict splice sites. To achieve the ab initio prediction, we used human genomic data to train our neural network. An iterative approach is adopted to reconstruct the dataset, which tackles the data unbalance problem and forces the model to learn more features of splice sites. The proposed CNN obtains the classification accuracy of 90.25%, which is 10% higher than the existing algorithms. The method outperforms other existing methods in terms of area under receiver operating characteristics (AUC), recall, precision, and F1 score. Furthermore, SpliceFinder can find the exact position of splice sites on long genomic sequences with a sliding window. Compared with other state-of-the-art splice site prediction tools, SpliceFinder generates results in about half lower false positive while keeping recall higher than 0.8. Also, SpliceFinder captures the non-canonical splice sites. In addition, SpliceFinder performs well on the genomic sequences of Drosophila melanogaster, Mus musculus, Rattus, and Danio rerio without retraining.

Conclusion

Based on CNN, we have proposed a new ab initio splice site prediction tool, SpliceFinder, which generates less false positives and can detect non-canonical splice sites. Additionally, SpliceFinder is transferable to other species without retraining. The source code and additional materials are available at https://gitlab.deepomics.org/wangruohan/SpliceFinder.

Molecular Cloning and Characterization of G Alpha Proteins from the Western Tarnished Plant Bug, Lygus hesperus

The Gα subunits of heterotrimeric G proteins play critical roles in the activation of diverse signal transduction cascades. However, the role of these genes in chemosensation remains to be fully elucidated. To initiate a comprehensive survey of signal transduction genes, we used homology-based cloning methods and transcriptome data mining to identity Gα subunits in the western tarnished plant bug (Lygus hesperus Knight). Among the nine sequences identified were single variants of the Gαi, Gαo, Gαs, and Gα12 subfamilies and five alternative splice variants of the Gαq subfamily. Sequence alignment and phylogenetic analyses of the putative L. hesperus Gα subunits support initial classifications and are consistent with established evolutionary relationships. End-point PCR-based profiling of the transcripts indicated head specific expression for LhGαq4, and largely ubiquitous expression, albeit at varying levels, for the other LhGα transcripts. All subfamilies were amplified from L. hesperus chemosensory tissues, suggesting potential roles in olfaction and/or gustation. Immunohistochemical staining of cultured insect cells transiently expressing recombinant His-tagged LhGαi, LhGαs, and LhGαq1 revealed plasma membrane targeting, suggesting the respective sequences encode functional G protein subunits.

A comprehensive survey of non-canonical splice sites in the human transcriptome

We uncovered the diversity of non-canonical splice sites at the human transcriptome using deep transcriptome profiling. We mapped a total of 3.7 billion human RNA-seq reads and developed a set of stringent filters to avoid false non-canonical splice site detections. We identified 184 splice sites with non-canonical dinucleotides and U2/U12-like consensus sequences. We selected 10 of the herein identified U2/U12-like non-canonical splice site events and successfully validated 9 of them via reverse transcriptase-polymerase chain reaction and Sanger sequencing. Analyses of the 184 U2/U12-like non-canonical splice sites indicate that 51% of them are not annotated in GENCODE. In addition, 28% of them are conserved in mouse and 76% are involved in alternative splicing events, some of them with tissue-specific alternative splicing patterns. Interestingly, our analysis identified some U2/U12-like non-canonical splice sites that are converted into canonical splice sites by RNA A-to-I editing. Moreover, the U2/U12-like non-canonical splice sites have a differential distribution of splicing regulatory sequences, which may contribute to their recognition and regulation. Our analysis provides a high-confidence group of U2/U12-like non-canonical splice sites, which exhibit distinctive features among the total human splice sites.

Cytoarchitecture and ultrastructure of neural stem cell niches and neurogenic complexes maintaining adult neurogenesis in the olfactory midbrain of spiny lobsters, Panulirus argus

New interneurons are continuously generated in small proliferation zones within neuronal somata clusters in the olfactory deutocerebrum of adult decapod crustaceans. Each proliferation zone is connected to a clump of cells containing one neural stem cell (i.e., adult neuroblast), thus forming a "neurogenic complex." Here we provide a detailed analysis of the cytoarchitecture of neurogenic complexes in adult spiny lobsters, Panulirus argus, based on transmission electron microscopy and labeling with cell-type-selective markers. The clump of cells is composed of unique bipolar clump-forming cells that collectively completely envelop the adult neuroblast and are themselves ensheathed by a layer of processes of multipolar cell body glia. An arteriole is attached to the clump of cells, but dye perfusion experiments show that hemolymph has no access to the interior of the clump of cells. Thus, the clump of cells fulfills morphological criteria of a protective stem cell niche, with clump-forming cells constituting the adult neuroblast's microenvironment together with the cell body glia processes separating it from other tissue components. Bromodeoxyuridine pulse-chase experiments with short survival times suggest that adult neuroblasts are not quiescent but rather cycle actively during daytime. We propose a cell lineage model in which an asymmetrically dividing adult neuroblast repopulates the pool of neuronal progenitor cells in the associated proliferation zone. In conclusion, as in mammalian brains, adult neurogenesis in crustacean brains is fueled by neural stem cells that are maintained by stem cell niches that preserve elements of the embryonic microenvironment and contain glial and vascular elements.

Proteolytic activation and function of the cytokine Spätzle in the innate immune response of a lepidopteran insect, Manduca sexta

The innate immune response of insects includes induced expression of genes encoding a variety of antimicrobial peptides. The signaling pathways that stimulate this gene expression have been well characterized by genetic analysis in Drosophila melanogaster, but are not well understood in most other insect species. One such pathway involves proteolytic activation of a cytokine called Spätzle, which functions in dorsal-ventral patterning in early embryonic development and in the antimicrobial immune response in larvae and adults. We have investigated the function of Spätzle in a lepidopteran insect, Manduca sexta, in which hemolymph proteinases activated during immune responses have been characterized biochemically. Two cDNA isoforms for M. sexta Spätzle-1 differ because of alternative splicing, resulting in a 10 amino acid residue insertion in the pro-region of proSpätzle-1B that is not present in proSpätzle-1A. The proSpätzle-1A cDNA encodes a 32.7 kDa polypeptide that is 23% and 44% identical to D. melanogaster and Bombyx mori Spätzle-1, respectively. Recombinant proSpätzle-1A was a disulfide-linked homodimer. M. sexta hemolymph proteinase 8 cleaved proSpätzle-1A to release Spätzle-C108, a dimer of the C-terminal 108 residue cystine-knot domain. Injection of Spätzle-C108, but not proSpätzle-1A, into larvae stimulated expression of several antimicrobial peptides and proteins, including attacin-1, cecropin-6, moricin, lysozyme, and the immunoglobulin domain protein hemolin, but did not significantly affect the expression of two bacteria-inducible pattern recognition proteins, immulectin-2 and beta-1,3-glucan recognition protein-2. The results of this and other recent studies support a model for a pathway in which the clip-domain proteinase pro-hemolymph proteinase 6 becomes activated in plasma upon exposure to Gram-negative or Gram-positive bacteria or to beta-1,3-glucan. Hemolymph proteinase 6 then activates pro-hemolymph proteinase 8, which in turn activates Spätzle-1. The resulting Spätzle-C108 dimer is likely to function as a ligand to activate a Toll pathway in M. sexta as a response to a wide variety of microbial challenges, stimulating a broad response to infection. Structured digital abstract * MINT-7295125: Spätzle 1A (uniprotkb:C8BMD1) and Spätzle 1A (uniprotkb:C8BMD1) bind (MI:0407) by comigration in gel electrophoresis (MI:0807).

Violating the splicing rules: TG dinucleotides function as alternative 3' splice sites in U2-dependent introns

TG dinucleotides functioning as alternative 3' splice sites were identified and experimentally verified in 36 human genes.

Background

Despite some degeneracy of sequence signals that govern splicing of eukaryotic pre-mRNAs, it is an accepted rule that U2-dependent introns exhibit the 3' terminal dinucleotide AG. Intrigued by anecdotal evidence for functional non-AG 3' splice sites, we carried out a human genome-wide screen.

Results

We identified TG dinucleotides functioning as alternative 3' splice sites in 36 human genes. The TG-derived splice variants were experimentally validated with a success rate of 92%. Interestingly, ratios of alternative splice variants are tissue-specific for several introns. TG splice sites and their flanking intron sequences are substantially conserved between orthologous vertebrate genes, even between human and frog, indicating functional relevance. Remarkably, TG splice sites are exclusively found as alternative 3' splice sites, never as the sole 3' splice site for an intron, and we observed a distance constraint for TG-AG splice site tandems.

Conclusion

Since TGs splice sites are exclusively found as alternative 3' splice sites, the U2 spliceosome apparently accomplishes perfect specificity for 3' AGs at an early splicing step, but may choose 3' TGs during later steps. Given the tiny fraction of TG 3' splice sites compared to the vast amount of non-viable TGs, cis-acting sequence signals must significantly contribute to splice site definition. Thus, we consider TG-AG 3' splice site tandems as promising subjects for studies on the mechanisms of 3' splice site selection.

Galpha encoding gene family of the malaria vector mosquito Anopheles gambiae: expression analysis and immunolocalization of AGalphaq and AGalphao in female antennae

To initiate a comprehensive investigation of chemosensory signal transduction downstream of odorant receptors, we identified and characterized the complete set of genes that encode G-protein alpha subunits in the genome of the malaria vector mosquito An. gambiae. Data are provided on the tissue-specific expression patterns of 10 corresponding aga-transcripts in adult mosquitoes and pre-imago developmental stages. Specific immunoreactivity in chemosensory hairs of female antennae provides evidence in support of the participation of a subset of AGalphaq isoforms in olfactory signal transduction in this mosquito. In contrast, AGalphao is localized along the flagellar axon bundle but is absent from chemosensory sensilla, which suggests that this G-protein alpha subunit does not participate in olfactory signal transduction.

Gsalpha is involved in sugar perception in Drosophila melanogaster

In Drosophila melanogaster, gustatory receptor genes (Grs) encode G-protein-coupled receptors (GPCRs) in gustatory receptor neurons (GRNs) and some olfactory receptor neurons. One of the Gr genes, Gr5a, encodes a sugar receptor that is expressed in a subset of GRNs and has been most extensively studied both molecularly and physiologically, but the G-protein alpha subunit (Galpha) that is coupled to this sugar receptor remains unknown. Here, we propose that Gs is the Galpha that is responsible for Gr5a-mediated sugar-taste transduction, based on the following findings: First, immunoreactivities against Gs were detected in a subset of GRNs including all Gr5a-expressing neurons. Second, trehalose-intake is reduced in flies heterozygous for null mutations in DGsalpha, a homolog of mammalian Gs, and trehalose-induced electrical activities in sugar-sensitive GRNs were depressed in those flies. Furthermore, expression of wild-type DGsalpha in sugar-sensitive GRNs in heterozygotic DGsalpha mutant flies rescued those impairments. Third, expression of double-stranded RNA for DGsalpha in sugar-sensitive GRNs depressed both behavioral and electrophysiological responses to trehalose. Together, these findings indicate that DGsalpha is involved in trehalose perception. We suggest that sugar-taste signals are processed through the Gsalpha-mediating signal transduction pathway in sugar-sensitive GRNs in Drosophila.

Presynaptic impairment of synaptic transmission in Drosophila embryos lacking Gs(alpha)

Gs(alpha) is a subunit of the heterotrimeric G-protein complex, expressed ubiquitously in all types of cells, including neurons. Drosophila larvae, which have mutations in the Gs(alpha) gene, are lethargic, suggesting an impairment of neuronal functions. In this study, we examined synaptic transmission at the neuromuscular synapse in Gs(alpha)-null (dgsR60) embryos shortly before they hatched. At low-frequency nerve stimulation, synaptic transmission in mutant embryos was not very different from that in controls. In contrast, facilitation during tetanic stimulation was minimal in dgsR60, and no post-tetanic potentiation was observed. Miniature synaptic currents (mSCs) were slightly smaller in amplitude and less frequent in dgsR60 embryos in normal-K+ saline. In high-K+ saline, mSCs with distinctly large amplitude occurred frequently in controls at late embryonic stages, whereas those mSCs were rarely observed in dgsR60 embryos, suggesting a developmental defect in the mutant. Using the Gal4-UAS expression system, we found that these phenotypes in dgsR60 were caused predominantly by lack of Gs(alpha) in presynaptic neurons and not in postsynaptic muscles. To test whether Gs(alpha) couples presynaptic modulator receptors to adenylyl cyclase (AC), we examined the responses of two known G-protein-coupled receptors in dgsR60 embryos. Both metabotropic glutamate and octopamine receptor responses were indistinguishable from those of controls, indicating that these receptors are not linked to AC by Gs(alpha). We therefore suggest that synaptic transmission is compromised in dgsR60 embryos because of presynaptic defects in two distinct processes; one is uncoupling between the yet-to-be-known modulator receptor and AC activation, and the other is a defect in synapse formation.

Analysis of canonical and non-canonical splice sites in mammalian genomes

A set of 43 337 splice junction pairs was extracted from mammalian GenBank annotated genes. Expressed sequence tag (EST) sequences support 22 489 of them. Of these, 98.71% contain canonical dinucleotides GT and AG for donor and acceptor sites, respectively; 0.56% hold non-canonical GC-AG splice site pairs; and the remaining 0.73% occurs in a lot of small groups (with a maximum size of 0.05%). Studying these groups we observe that many of them contain splicing dinucleotides shifted from the annotated splice junction by one position. After close examination of such cases we present a new classification consisting of only eight observed types of splice site pairs (out of 256 a priori possible combinations). EST alignments allow us to verify the exonic part of the splice sites, but many non-canonical cases may be due to intron sequencing errors. This idea is given substantial support when we compare the sequences of human genes having non-canonical splice sites deposited in GenBank by high throughput genome sequencing projects (HTG). A high proportion (156 out of 171) of the human non-canonical and EST-supported splice site sequences had a clear match in the human HTG. They can be classified after corrections as: 79 GC-AG pairs (of which one was an error that corrected to GC-AG), 61 errors that were corrected to GT-AG canonical pairs, six AT-AC pairs (of which two were errors that corrected to AT-AC), one case was produced from non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two cases left of supported non-canonical splice sites. If we assume that approximately the same situation is true for the whole set of annotated mammalian non-canonical splice sites, then the 99.24% of splice site pairs should be GT-AG, 0.69% GC-AG, 0.05% AT-AC and finally only 0.02% could consist of other types of non-canonical splice sites. We analyze several characteristics of EST-verified splice sites and build weight matrices for the major groups, which can be incorporated into gene prediction programs. We also present a set of EST-verified canonical splice sites larger by two orders of magnitude than the current one (22 199 entries versus approximately 600) and finally, a set of 290 EST-supported non-canonical splice sites. Both sets should be significant for future investigations of the splicing mechanism.

The heterotrimeric protein Go is required for the formation of heart epithelium in Drosophila

The gene encoding the alpha subunit of the Drosophila Go protein is expressed early in embryogenesis in the precursor cells of the heart tube, of the visceral muscles, and of the nervous system. This early expression coincides with the onset of the mesenchymal-epithelial transition to which are subjected the cardial cells and the precursor cells of the visceral musculature. This gene constitutes an appropriate marker to follow this transition. In addition, a detailed analysis of its expression suggests that the cardioblasts originate from two subpopulations of cells in each parasegment of the dorsal mesoderm that might depend on the wingless and hedgehog signaling pathways for both their determination and specification. In the nervous system, the expression of Goalpha shortly precedes the beginning of axonogenesis. Mutants produced in the Goalpha gene harbor abnormalities in the three tissues in which the gene is expressed. In particular, the heart does not form properly and interruptions in the heart epithelium are repeatedly observed, henceforth the brokenheart (bkh) name. Furthermore, in the bkh mutant embryos, the epithelial polarity of cardial cells was not acquired (or maintained) in various places of the cardiac tube. We predict that bkh might be involved in vesicular traffic of membrane proteins that is responsible for the acquisition of polarity.

Alternative splicing of the latency-related transcript of bovine herpesvirus 1 yields RNAs containing unique open reading frames

The latency-related transcript (LRT) of bovine herpesvirus 1 (BHV-1) is the only abundant viral RNA detected during latency. A previous study (A. Hossain, L. M. Schang, and C. Jones, J. Virol. 69:5345-5352, 1995) concluded that splicing of polyadenylated [poly(A)+] and splicing of nonpolyadenylated [poly(A)-] LRT are different. In this study, splice junction sites of LRT were identified. In trigeminal ganglia of acutely infected calves (1, 7, or 15 days postinfection [p.i.]) or in latently infected calves (60 days p.i.), alternative splicing of poly(A)+ LRT occurred. Productive viral gene expression in trigeminal ganglia is readily detected from 2 to 7 days p.i. but not at 15 days p.i. (L. M. Schang and C. Jones, J. Virol. 71:6786-6795, 1997), suggesting that certain aspects of a lytic infection occur in neurons and that these factors influence LRT splicing. Splicing of poly(A)- LRT was also detected in transfected COS-7 cells or infected MDBK cells. DNA sequence analysis of spliced LRT cDNAs, poly(A)+ or poly(A)-, revealed nonconsensus splice signals at exon/intron and intron/exon boundaries. The GC-AG splicing signal utilized by the herpes simplex virus type 1 latency-associated transcript in latently infected mice is also used by LRT in latently infected calves. Taken together, these results led us to hypothesize that (i) poly(A)+ LRT is spliced in trigeminal ganglia by neuron-specific factors, (ii) viral or virus-induced factors participate in splicing, and (iii) alternative splicing of LRT may result in protein isoforms which have novel biological properties.

Activation of protein kinase A-independent pathways by Gs alpha in Drosophila

One of the best-described transmembrane signal transduction mechanisms is based on receptor activation of the alpha subunit of the heterotrimeric G protein Gs, leading to stimulation of adenylyl cyclase and the production of cAMP. Intracellular cAMP is then thought to mediate its effects largely, if not entirely, by activation of protein kinase A and the subsequent phosphorylation of substrates which in turn control diverse cellular phenomena. In this report we demonstrate, by two different methods, that reduction or elimination of protein kinase A activity had no effect on phenotypes generated by activation of Gs alpha pathways in Drosophila wing epithelial cells. These genetic studies show that the Gs alpha pathway mediates its primary effects by a novel pathway in differentiating, wing epithelial cells. This novel pathway may in part be responsible for some of the complex, cell-specific responses observed following activation of this pathway in different cell types.

Drosophila stimulatory G protein alpha subunit activates mammalian adenylyl cyclase but interacts poorly with mammalian receptors: implications for receptor-G protein interaction

Heterotrimeric guanine nucleotide binding proteins (G proteins) transduce signals from cell-surface receptors to intracellular effector proteins. Two forms of stimulatory G protein (Gs) alpha-like subunit have been described in Drosophila melanogaster. To examine the function of these subunits we have used vaccinia virus vectors to express both proteins in cyc- cells, a murine S49 cell line deficient for Gs alpha activity. Receptor-independent activation of each Drosophila Gs alpha has demonstrated that both forms are capable of activating mammalian adenylyl cyclase and thus have the activity expected of stimulatory G proteins. However, the Drosophila Gs alpha subunits interact poorly with mammalian Gs-coupled receptors. These observations have helped to identify a region of high variability in Gs alpha proteins that may be important for receptor interactions.

Molecular characterization of the Drosophila melanogaster urate oxidase gene, an ecdysone-repressible gene expressed only in the malpighian tubules

The urate oxidase (UO) gene of Drosophila melanogaster is expressed during the third-instar larval and adult stages, exclusively within a subset of cells of the Malpighian tubules. The UO gene contains a 69-base-pair intron and encodes mature mRNAs of 1,224, 1,227, and 1,244 nucleotides, depending on the site of 3' endonucleolytic cleavage prior to polyadenylation. A direct repeat, 5'-AAGTGAGAGTGAT-3', is the proposed cis-regulatory element involved in 20-hydroxyecdysone repression of the UO gene. The deduced amino acid sequences of UO of D. melanogaster, rat, mouse, and pig and uricase II of soybean show 32 to 38% identity, with 22% of amino acid residues identical in all species. With use of P-element-mediated germ line transformation, 826 base pairs 5' and approximately 1,200 base pairs 3' of the D. melanogaster UO transcribed region contain all of the cis elements allowing for appropriate temporal regulation and Malpighian tubule-specific expression of the UO gene.

Molecular characterization of the Drosophila melanogaster urate oxidase gene, an ecdysone-repressible gene expressed only in the malpighian tubules

The urate oxidase (UO) gene of Drosophila melanogaster is expressed during the third-instar larval and adult stages, exclusively within a subset of cells of the Malpighian tubules. The UO gene contains a 69-base-pair intron and encodes mature mRNAs of 1,224, 1,227, and 1,244 nucleotides, depending on the site of 3' endonucleolytic cleavage prior to polyadenylation. A direct repeat, 5'-AAGTGAGAGTGAT-3', is the proposed cis-regulatory element involved in 20-hydroxyecdysone repression of the UO gene. The deduced amino acid sequences of UO of D. melanogaster, rat, mouse, and pig and uricase II of soybean show 32 to 38% identity, with 22% of amino acid residues identical in all species. With use of P-element-mediated germ line transformation, 826 base pairs 5' and approximately 1,200 base pairs 3' of the D. melanogaster UO transcribed region contain all of the cis elements allowing for appropriate temporal regulation and Malpighian tubule-specific expression of the UO gene.