Massive expansions of Dscam splicing diversity via staggered homologous recombination during arthropod evolution

Following the Evolutionary Paths of Dscam1 Proteins toward Highly Specific Homophilic Interactions

Many adhesion proteins, evolutionarily related through gene duplication, exhibit distinct and precise interaction preferences and affinities crucial for cell patterning. Yet, the evolutionary paths by which these proteins acquire new specificities and prevent cross-interactions within their family members remain unknown. To bridge this gap, this study focuses on Drosophila Down syndrome cell adhesion molecule-1 (Dscam1) proteins, which are cell adhesion proteins that have undergone extensive gene duplication. Dscam1 evolved under strong selective pressure to achieve strict homophilic recognition, essential for neuronal self-avoidance and patterning. Through a combination of phylogenetic analyses, ancestral sequence reconstruction, and cell aggregation assays, we studied the evolutionary trajectory of Dscam1 exon 4 across various insect lineages. We demonstrated that recent Dscam1 duplications in the mosquito lineage bind with strict homophilic specificities without any cross-interactions. We found that ancestral and intermediate Dscam1 isoforms maintained their homophilic binding capabilities, with some intermediate isoforms also engaging in promiscuous interactions with other paralogs. Our results highlight the robust selective pressure for homophilic specificity integral to the Dscam1 function within the process of neuronal self-avoidance. Importantly, our study suggests that the path to achieving such selective specificity does not introduce disruptive mutations that prevent self-binding but includes evolutionary intermediates that demonstrate promiscuous heterophilic interactions. Overall, these results offer insights into evolutionary strategies that underlie adhesion protein interaction specificities.

Structure and evolution of neuronal wiring receptors and ligands

One of the fundamental properties of a neuronal circuit is the map of its connections. The cellular and developmental processes that allow for the growth of axons and dendrites, selection of synaptic targets, and formation of functional synapses use neuronal surface receptors and their interactions with other surface receptors, secreted ligands, and matrix molecules. Spatiotemporal regulation of the expression of these receptors and cues allows for specificity in the developmental pathways that wire stereotyped circuits. The families of molecules controlling axon guidance and synapse formation are generally conserved across animals, with some important exceptions, which have consequences for neuronal connectivity. Here, we summarize the distribution of such molecules across multiple taxa, with a focus on model organisms, evolutionary processes that led to the multitude of such molecules, and functional consequences for the diversification or loss of these receptors.

Hidden RNA pairings counteract the "first-come, first-served" splicing principle to drive stochastic choice in Dscam1 splice variants

Drosophila melanogaster Dscam1 encodes 38,016 isoforms via mutually exclusive splicing; however, the regulatory mechanism behind this is not fully understood. Here, we found a set of hidden RNA secondary structures that balance the stochastic choice of Dscam1 splice variants (designated balancer RNA secondary structures). In vivo mutational analyses revealed the dual function of these balancer interactions in driving the stochastic choice of splice variants, through enhancement of the inclusion of distal exon 6s by cooperating with docking site–selector pairing to form a stronger multidomain pre-mRNA structure and through simultaneous repression of the inclusion of proximal exon 6s by antagonizing their docking site–selector pairings. Thus, we provide an elegant molecular model based on competition and cooperation between two sets of docking site–selector and balancer pairings, which counteracts the “first-come, first-served” principle. Our findings provide conceptual and mechanistic insight into the dynamics and functions of long-range RNA secondary structures.

Complex RNA Secondary Structures Mediate Mutually Exclusive Splicing of Coleoptera Dscam1

Mutually exclusive splicing is an important mechanism for expanding protein diversity. An extreme example is the Down syndrome cell adhesion molecular (Dscam1) gene of insects, containing four clusters of variable exons (exons 4, 6, 9, and 17), which potentially generates tens of thousands of protein isoforms through mutually exclusive splicing, of which regulatory mechanisms are still elusive. Here, we systematically analyzed the variable exon 4, 6, and 9 clusters of Dscam1 in Coleoptera species. Through comparative genomics and RNA secondary structure prediction, we found apparent evidence that the evolutionarily conserved RNA base pairing mediates mutually exclusive splicing in the Dscam1 exon 4 cluster. In contrast to the fly exon 6, most exon 6 selector sequences in Coleoptera species are partially located in the variable exon region. Besides, bidirectional RNA–RNA interactions are predicted to regulate the mutually exclusive splicing of variable exon 9 of Dscam1. Although the docking sites in exon 4 and 9 clusters are clade specific, the docking sites-selector base pairing is conserved in secondary structure level. In short, our result provided a mechanistic framework for the application of long-range RNA base pairings in regulating the mutually exclusive splicing of Coleoptera Dscam1.

RNA secondary structures in Dscam1 mutually exclusive splicing: unique evolutionary signature from the midge

The Drosophila melanogaster gene Dscam1 potentially generates 38,016 distinct isoforms via mutually exclusive splicing, which are required for both nervous and immune functions. However, the mechanism underlying splicing regulation remains obscure. Here we show apparent evolutionary signatures characteristic of competing RNA secondary structures in exon clusters 6 and 9 of Dscam1 in the two midge species (Belgica antarctica and Clunio marinus). Surprisingly, midge Dscam1 encodes only ∼6000 different isoforms through mutually exclusive splicing. Strikingly, the docking site of the exon 6 cluster is conserved in almost all insects and crustaceans but is specific in the midge; however, the docking site-selector base-pairings are conserved. Moreover, the docking site is complementary to all predicted selector sequences downstream from every variable exon 9 of the midge Dscam1, which is in accordance with the broad spectrum of their isoform expression. This suggests that these cis-elements mainly function through the formation of long-range base-pairings. This study provides a vital insight into the evolution and mechanism of Dscam1 alternative splicing.

Dscam in immunity: A question of diversity in insects and crustaceans

In insects and crustaceans, thousands of Down syndrome cell adhesion molecules (Dscam) can be generated by alternative splicing of variable exons from a single-locus gene, Dscam-hv. This extraordinarily versatile gene (38,016 protein isoforms produced in Drosophila) was first proposed to be involved in exon guidance and subsequently implicated in immunity as a hypervariable immune molecule. Almost 20 y after discovery of Dscam-hv, there have been many studies in insects and crustaceans regarding roles of Dscam in immunity, with many similarities and concurrently, many differences. Here, we review the current status of Dscam-hv, presented as a comparison of similarities and differences in insects and crustaceans and discuss hypotheses of Dscam functions in immunity.

Acute thiamethoxam toxicity in honeybees is not enhanced by common fungicide and herbicide and lacks stress-induced changes in mRNA splicing

Securing food supply for a growing population is a major challenge and heavily relies on the use of agrochemicals to maximize crop yield. It is increasingly recognized, that some neonicotinoid insecticides have a negative impact on non-target organisms, including important pollinators such as the European honeybee Apis mellifera. Toxicity of neonicotinoids may be enhanced through simultaneous exposure with additional pesticides, which could help explain, in part, the global decline of honeybee colonies. Here we examined whether exposure effects of the neonicotinoid thiamethoxam on bee viability are enhanced by the commonly used fungicide carbendazim and the herbicide glyphosate. We also analysed alternative splicing changes upon pesticide exposure in the honeybee. In particular, we examined transcripts of three genes: (i) the stress sensor gene X box binding protein-1 (Xbp1), (ii) the Down Syndrome Cell Adhesion Molecule (Dscam) gene and iii) the embryonic lethal/abnormal visual system (elav) gene, which are important for neuronal function. Our results showed that acute thiamethoxam exposure is not enhanced by carbendazim, nor glyphosate. Toxicity of the compounds did not trigger stress-induced, alternative splicing in the analysed mRNAs, thereby leaving dormant a cellular response pathway to these man-made environmental perturbations.

Plasmid-based gap-repair recombineered transgenes reveal a central role for introns in mutually exclusive alternative splicing in Down Syndrome Cell Adhesion Molecule exon 4

Alternative splicing is a key feature of human genes, yet studying its regulation is often complicated by large introns. The Down Syndrome Cell Adhesion Molecule (Dscam) gene from Drosophila is one of the most complex genes generating vast molecular diversity by mutually exclusive alternative splicing. To resolve how alternative splicing in Dscam is regulated, we first developed plasmid-based UAS reporter genes for the Dscam variable exon 4 cluster and show that its alternative splicing is recapitulated by GAL4-mediated expression in neurons. We then developed gap-repair recombineering to very efficiently manipulate these large reporter plasmids in Escherichia coli using restriction enzymes or sgRNA/Cas9 DNA scission to capitalize on the many benefits of plasmids in phiC31 integrase-mediated transgenesis. Using these novel tools, we show that inclusion of Dscam exon 4 variables differs little in development and individual flies, and is robustly determined by sequences harbored in variable exons. We further show that introns drive selection of both proximal and distal variable exons. Since exon 4 cluster introns lack conserved sequences that could mediate robust long-range base-pairing to bring exons into proximity for splicing, our data argue for a central role of introns in mutually exclusive alternative splicing of Dscam exon 4 cluster.

B52 promotes alternative splicing of Dscam in Chinese mitten crab, Eriocheir sinensis

B52 is a member of the classical serine/arginine (SR)-rich proteins, which are phylogenetically conserved and play significant roles in mRNA maturation, including alternative splicing. In the present study, the docking site, selector sequences and locus control region of the Chinese mitten crab (Eriocheir sinensis) Down syndrome cell adhesion molecule (EsDscam) were identified. Alternative splicing of Dscam is essential to generate different isoforms. We also isolated and characterised the B52 gene from E. sinensis (EsB52). The 876 bp open reading frame of EsB52 encodes a 291 amino acid residue polypeptide, and EsB52 has two RNA recognition motifs (RRMs) at the N-terminus and an arginine/serine-rich domain at the C-terminus. Each RRM contains two degenerate short submotifs, RNP-1 and RNP2. Analysis of tissue distribution revealed that EsB52 mRNA expression was widespread in all tested tissues, and especially high in brain and hemocytes. In hemocytes, EsB52 was upregulated significantly after stimulation with pathogen-associated molecular patterns and bacteria. Furthermore, EsB52 RNAi decreased the number of Ig7 inclusion in mRNA rather than Ig2 or Ig3. Taken together, these findings suggest that EsB52 acts as an alternative splicing activator of EsDscam.

Revisiting Dscam diversity: lessons from clustered protocadherins

The complexity of neuronal wiring relies on the extraordinary recognition diversity of cell surface molecules. Drosophila Dscam1 and vertebrate clustered protocadherins (Pcdhs) are two classic examples of the striking diversity from a complex genomic locus, wherein the former encodes more than 10,000 distinct isoforms via alternative splicing, while the latter employs alternative promoters to attain isoform diversity. These structurally unrelated families show remarkably striking molecular parallels and even similar functions. Recent studies revealed a novel Dscam gene family with tandemly arrayed 5' cassettes in Chelicerata (e.g., the scorpion Mesobuthus martensii and the tick Ixodes scapularis), similar to vertebrate clustered Pcdhs. Likewise, octopus shows a more remarkable expansion of the Pcdh isoform repertoire than human. These discoveries of Dscam and Pcdh diversification reshape the evolutionary landscape of recognition molecule diversity and provide a greater understanding of convergent molecular strategies for isoform diversity. This article reviews new insights into the evolution, regulatory mechanisms, and functions of Dscam and Pcdh isoform diversity. In particular, the convergence of clustered Dscams and Pcdhs is highlighted.

Mutually exclusive alternative splicing of pre-mRNAs

Pre-mRNA alternative splicing is an important mechanism used to expand protein diversity in higher eukaryotes, and mutually exclusive splicing is a specific type of alternative splicing in which only one of the exons in a cluster is included in functional transcripts. The most extraordinary example of this is the Drosophila melanogaster Down's syndrome cell adhesion molecule gene (Dscam), which potentially encodes 38,016 different isoforms through mutually exclusive splicing. Mutually exclusive splicing is a unique and challenging model that can be used to elucidate the evolution, regulatory mechanism, and function of alternative splicing. The use of new approaches has not only greatly expanded the mutually exclusive exome, but has also enabled the systematic analyses of single-cell alternative splicing during development. Furthermore, the identification of long-range RNA secondary structures provides a mechanistic framework for the regulation of mutually exclusive splicing (i.e., Dscam splicing). This article reviews recent insights into the identification, underlying mechanism, and roles of mutually exclusive splicing. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

The landscape of human mutually exclusive splicing

Mutually exclusive splicing of exons is a mechanism of functional gene and protein diversification with pivotal roles in organismal development and diseases such as Timothy syndrome, cardiomyopathy and cancer in humans. In order to obtain a first genomewide estimate of the extent and biological role of mutually exclusive splicing in humans, we predicted and subsequently validated mutually exclusive exons (MXEs) using 515 publically available RNA‐Seq datasets. Here, we provide evidence for the expression of over 855 MXEs, 42% of which represent novel exons, increasing the annotated human mutually exclusive exome more than fivefold. The data provide strong evidence for the existence of large and multi‐cluster MXEs in higher vertebrates and offer new insights into MXE evolution. More than 82% of the MXE clusters are conserved in mammals, and five clusters have homologous clusters in Drosophila. Finally, MXEs are significantly enriched in pathogenic mutations and their spatio‐temporal expression might predict human disease pathology.

Role and convergent evolution of competing RNA secondary structures in mutually exclusive splicing

Exon or cassette duplication is an important means of expanding protein and functional diversity through mutually exclusive splicing. However, the mechanistic basis of this process in non-arthropod species remains poorly understood. Here, we demonstrate that MRP1 genes underwent tandem exon duplication in Nematoda, Platyhelminthes, Annelida, Mollusca, Arthropoda, Echinodermata, and early-diverging Chordata but not in late-diverging vertebrates. Interestingly, these events were of independent origin in different phyla, suggesting convergent evolution of alternative splicing. Furthermore, we showed that multiple sets of clade-conserved RNA pairings evolved to guide species-specific mutually exclusive splicing in Arthropoda. Importantly, we also identified a similar structural code in MRP exon clusters of the annelid, Capitella teleta, and chordate, Branchiostoma belcheri, suggesting an evolutionarily conserved competing pairing-guided mechanism in bilaterians. Taken together, these data reveal the molecular determinants and RNA pairing-guided evolution of species-specific mutually exclusive splicing spanning more than 600 million years of bilaterian evolution. These findings have a significant impact on our understanding of the evolution of and mechanism underpinning isoform diversity and complex gene structure.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

Assessing long-distance RNA sequence connectivity via RNA-templated DNA-DNA ligation

Many RNAs, including pre-mRNAs and long non-coding RNAs, can be thousands of nucleotides long and undergo complex post-transcriptional processing. Multiple sites of alternative splicing within a single gene exponentially increase the number of possible spliced isoforms, with most human genes currently estimated to express at least ten. To understand the mechanisms underlying these complex isoform expression patterns, methods are needed that faithfully maintain long-range exon connectivity information in individual RNA molecules. In this study, we describe SeqZip, a methodology that uses RNA-templated DNA–DNA ligation to retain and compress connectivity between distant sequences within single RNA molecules. Using this assay, we test proposed coordination between distant sites of alternative exon utilization in mouse Fn1, and we characterize the extraordinary exon diversity of Drosophila melanogaster Dscam1.

DOI:http://dx.doi.org/10.7554/eLife.03700.001

A flow chart can show how an outcome can be achieved from a particular start point by breaking down an activity into a list of possible steps. Often, a flow chart contains several alternative steps, not all of which are taken every time the flow chart is used. The same can be said of genes, which are biological instructions that often contain many options within their DNA sequences.

Proteins—which perform many roles in cells—are built following the instructions contained in genes. First, the DNA sequence of the gene is copied. This produces a molecule of ribonucleic acid (RNA), which is able to move around the cell to find the machinery that can use the genetic information to make a protein. Genes and their RNA copies contain instructions with more steps—called exons—than are necessary to make a working protein, so extra exons are removed (‘spliced’) from the RNA copies. Different combinations of exons can be removed, so splicing can make different versions of the RNA called isoforms. These allow a single gene to build many different proteins. In fruit flies, for example, the different exons of the gene Dscam1 can be spliced into one of 38,016 unique RNA isoforms.

Current technology only allows researchers to deduce the sequence of RNA molecules by combining sequences recorded from short fragments of the molecule. However, before splicing, RNA molecules tend to be much longer than this, so this restricts our understanding of the RNA isoforms found in cells. Here, Roy et al. devised and tested a new method called SeqZip to solve this problem.

SeqZip uses short fragments of DNA called ligamers that can only stick to the sections of RNA that will remain after the molecule has been spliced. After splicing, the ligamers can be stuck together to make a DNA replica of the spliced RNA. The end product is at least 49 times shorter than the original RNA, so it is easier to sequence. In addition, the combinations of the ligamers in the DNA replica show which exons of a specific gene are kept and which ones are spliced out.

To test the method, Roy et al. studied a mouse gene that has six RNA isoforms. SeqZip reduced the length of the RNA by five times and made it possible to measure how frequently the different isoforms naturally arise. Roy et al. also used SeqZip to work out which isoforms of the Dscam1 gene are used at different stages in the life of fruit fly larvae. SeqZip can provide insights into how complex organisms like flies, mice, and humans have evolved with relatively few—a little over 20,000—genes in their genomes.

DOI:http://dx.doi.org/10.7554/eLife.03700.002

Improvement of Dscam homophilic binding affinity throughout Drosophila evolution

Background

Drosophila Dscam1 is a cell-surface protein that plays important roles in neural development and axon tiling of neurons. It is known that thousands of isoforms bind themselves through specific homophilic interactions, a process which provides the basis for cellular self-recognition. Detailed biochemical studies of specific isoforms strongly suggest that homophilic binding, i.e. the formation of homodimers by identical Dscam1 isomers, is of great importance for the self-avoidance of neurons. Due to experimental limitations, it is currently impossible to measure the homophilic binding affinities for all 19,000 potential isoforms.

Results

Here we reconstructed the DNA sequences of an ancestral Dscam form (which likely existed approximately 40 ~ 50 million years ago) using a comparative genomic approach. On the basis of this sequence, we established a working model to predict the self-binding affinities of all isoforms in both the current and the ancestral genome, using machine-learning methods. Detailed computational analysis was performed to compare the self-binding affinities of all isoforms present in these two genomes. Our results revealed that 1) isoforms containing newly derived variable domains exhibit higher self-binding affinities than those with conserved domains, and 2) current isoforms display higher self-binding affinities than their counterparts in the ancient genome. As thousands of Dscam isoforms are needed for the self-avoidance of the neuron, we propose that an increase in self-binding affinity provides the basis for the successful evolution of the arthropod brain.

Conclusions

Our data presented here provide an excellent model for future experimental studies of the binding behavior of Dscam isoforms. The results of our analysis indicate that evolution favored the rise of novel variable domains thanks to their higher self-binding affinities, rather than selection merely on the basis of simple expansion of isoform diversity, as that this particular selection process would have established the powerful mechanisms required for neuronal self-avoidance. Thus, we reveal here a new molecular mechanism for the successful evolution of arthropod brains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0186-z) contains supplementary material, which is available to authorized users.

More than one way to produce protein diversity: duplication and limited alternative splicing of an adhesion molecule gene in basal arthropods

Exon duplication and alternative splicing evolved multiple times in metazoa and are of overall importance in shaping genomes and allowing organisms to produce many fold more proteins than there are genes in the genome. No other example is as striking as the one of the Down syndrome cell adhesion molecule (Dscam) of insects and crustaceans (pancrustaceans) involved in the nervous system differentiation and in the immune system. To elucidate the evolutionary history of this extraordinary gene, we investigated Dscam homologs in two basal arthropods, the myriapod Strigamia maritima and the chelicerate Ixodes scapularis. In both, Dscam diversified extensively by whole gene duplications resulting in multigene expansions. Within some of the S. maritima genes, exons coding for one of the immunoglobulin domains (Ig7) duplicated and are mutually exclusively alternatively spliced. Our results suggest that Dscam diversification was selected independently in chelicerates, myriapods, and pancrustaceans and that the usage of Dscam diversity by immune cells evolved for the first time in basal arthropods. We propose an evolutionary scenario for the appearance of the highly variable Dscam gene of pancrustaceans, adding to the understanding of how alternative splicing, exon, and gene duplication contribute to create molecular diversity associated with potentially new cellular functions.

Mechanisms of Drosophila Dscam mutually exclusive splicing regulation

Alternative splicing of pre-mRNA is a major mechanism to increase protein diversity in higher eukaryotes. Dscam, the Drosophila homologue of human DSCAM (Down's syndrome cell adhesion molecule), generates up to 38016 isoforms through mutually exclusive splicing in four variable exon clusters. This enormous molecular diversity is functionally important for wiring of the nervous system and phagocytosis of invading pathogens. Current models explaining this complex splicing regulation include a default repressed state of the variable exon clusters to prevent the splicing together of adjacent exons, the presence of RNA secondary structures important for the release of one specific variable exon from the repressed state and combinatorial interaction of RNA-binding proteins for choosing a specific exon.

The evolution of Dscam genes across the arthropods

BACKGROUND

One way of creating phenotypic diversity is through alternative splicing of precursor mRNAs. A gene that has evolved a hypervariable form is Down syndrome cell adhesion molecule (Dscam-hv), which in Drosophila melanogaster can produce thousands of isoforms via mutually exclusive alternative splicing. The extracellular region of this protein is encoded by three variable exon clusters, each containing multiple exon variants. The protein is vital for neuronal wiring where the extreme variability at the somatic level is required for axonal guidance, and it plays a role in immunity where the variability has been hypothesised to relate to recognition of different antigens. Dscam-hv has been found across the Pancrustacea. Additionally, three paralogous non-hypervariable Dscam-like genes have also been described for D. melanogaster. Here we took a bioinformatics approach, building profile Hidden Markov Models to search across species for putative orthologs to the Dscam genes and for hypervariable alternatively spliced exons, and inferring the phylogenetic relationships among them. Our aims were to examine whether Dscam orthologs exist outside the Bilateria, whether the origin of Dscam-hv could lie outside the Pancrustacea, when the Dscam-like orthologs arose, how many alternatively spliced exons of each exon cluster were present in the most common recent ancestor, and how these clusters evolved.

RESULTS

Our results suggest that the origin of Dscam genes may lie after the split between the Cnidaria and the Bilateria and supports the hypothesis that Dscam-hv originated in the common ancestor of the Pancrustacea. Our phylogeny of Dscam gene family members shows six well-supported clades: five containing Dscam-like genes and one containing all the Dscam-hv genes, a seventh clade contains arachnid putative Dscam genes. Furthermore, the exon clusters appear to have experienced different evolutionary histories.

CONCLUSIONS

Dscam genes have undergone independent duplication events in the insects and in an arachnid genome, which adds to the more well-known tandem duplications that have taken place within Dscam-hv genes. Therefore, two forms of gene expansion seem to be active within this gene family. The evolutionary history of this dynamic gene family will be further unfolded as genomes of species from more disparate groups become available.

Innate immune complexity in the purple sea urchin: diversity of the sp185/333 system

The California purple sea urchin, Strongylocentrotus purpuratus, is a long-lived echinoderm with a complex and sophisticated innate immune system. There are several large gene families that function in immunity in this species including the Sp185/333 gene family that has ∼50 (±10) members. The family shows intriguing sequence diversity and encodes a broad array of diverse yet similar proteins. The genes have two exons of which the second encodes the mature protein and has repeats and blocks of sequence called elements. Mosaics of element patterns plus single nucleotide polymorphisms-based variants of the elements result in significant sequence diversity among the genes yet maintains similar structure among the members of the family. Sequence of a bacterial artificial chromosome insert shows a cluster of six, tightly linked Sp185/333 genes that are flanked by GA microsatellites. The sequences between the GA microsatellites in which the Sp185/333 genes and flanking regions are located, are much more similar to each other than are the sequences outside the microsatellites suggesting processes such as gene conversion, recombination, or duplication. However, close linkage does not correspond with greater sequence similarity compared to randomly cloned and sequenced genes that are unlikely to be linked. There are three segmental duplications that are bounded by GAT microsatellites and include three almost identical genes plus flanking regions. RNA editing is detectible throughout the mRNAs based on comparisons to the genes, which, in combination with putative post-translational modifications to the proteins, results in broad arrays of Sp185/333 proteins that differ among individuals. The mature proteins have an N-terminal glycine-rich region, a central RGD motif, and a C-terminal histidine-rich region. The Sp185/333 proteins are localized to the cell surface and are found within vesicles in subsets of polygonal and small phagocytes. The coelomocyte proteome shows full-length and truncated proteins, including some with missense sequence. Current results suggest that both native Sp185/333 proteins and a recombinant protein bind bacteria and are likely important in sea urchin innate immunity.

Widespread recurrent evolution of genomic features

The recent explosion of genome sequences from all major phylogenetic groups has unveiled an unexpected wealth of cases of recurrent evolution of strikingly similar genomic features in different lineages. Here, we review the diverse known types of recurrent evolution in eukaryotic genomes, with a special focus on metazoans, ranging from reductive genome evolution to origins of splice-leader trans-splicing, from tandem exon duplications to gene family expansions. We first propose a general classification scheme for evolutionary recurrence at the genomic level, based on the type of driving force-mutation or selection-and the environmental and genomic circumstances underlying these forces. We then discuss various cases of recurrent genomic evolution under this scheme. Finally, we provide a broader context for repeated genomic evolution, including the unique relationship of genomic recurrence with the genotype-phenotype map, and the ways in which the study of recurrent genomic evolution can be used to understand fundamental evolutionary processes.

Population genetics of duplicated alternatively spliced exons of the Dscam gene in Daphnia and Drosophila

In insects and crustaceans, the Down syndrome cell adhesion molecule (Dscam) occurs in many different isoforms. These are produced by mutually exclusive alternative splicing of dozens of tandem duplicated exons coding for parts or whole immunoglobulin (Ig) domains of the Dscam protein. This diversity plays a role in the development of the nervous system and also in the immune system. Structural analysis of the protein suggested candidate epitopes where binding to pathogens could occur. These epitopes are coded by regions of the duplicated exons and are therefore diverse within individuals. Here we apply molecular population genetics and molecular evolution analyses using Daphnia magna and several Drosophila species to investigate the potential role of natural selection in the divergence between orthologs of these duplicated exons among species, as well as between paralogous exons within species. We found no evidence for a role of positive selection in the divergence of these paralogous exons. However, the power of this test was low, and the fact that no signs of gene conversion between paralogous exons were found suggests that paralog diversity may nonetheless be maintained by selection. The analysis of orthologous exons in Drosophila and in Daphnia revealed an excess of non-synonymous polymorphisms in the epitopes putatively involved in pathogen binding. This may be a sign of balancing selection. Indeed, in Dr. melanogaster the same derived non-synonymous alleles segregate in several populations around the world. Yet other hallmarks of balancing selection were not found. Hence, we cannot rule out that the excess of non-synonymous polymorphisms is caused by segregating slightly deleterious alleles, thus potentially indicating reduced selective constraints in the putative pathogen binding epitopes of Dscam.

Introns within ribosomal protein genes regulate the production and function of yeast ribosomes

In budding yeast, the most abundantly spliced pre-mRNAs encode ribosomal proteins (RPs). To investigate the contribution of splicing to ribosome production and function, we systematically eliminated introns from all RP genes to evaluate their impact on RNA expression, pre-rRNA processing, cell growth, and response to stress. The majority of introns were required for optimal cell fitness or growth under stress. Most introns are found in duplicated RP genes, and surprisingly, in the majority of cases, deleting the intron from one gene copy affected the expression of the other in a nonreciprocal manner. Consistently, 70% of all duplicated genes were asymmetrically expressed, and both introns and gene deletions displayed copy-specific phenotypic effects. Together, our results indicate that splicing in yeast RP genes mediates intergene regulation and implicate the expression ratio of duplicated RP genes in modulating ribosome function.

Invertebrate immune diversity

The arms race between hosts and pathogens (and other non-self) drives the molecular diversification of immune response genes in the host. Over long periods of evolutionary time, many different defense strategies have been employed by a wide variety of invertebrates. We review here penaeidins and crustins in crustaceans, the allorecognition system encoded by fuhc, fester and Uncle fester in a colonial tunicate, Dscam and PGRPs in arthropods, FREPs in snails, VCBPs in protochordates, and the Sp185/333 system in the purple sea urchin. Comparisons among immune systems, including those reviewed here have not identified an immune specific regulatory "genetic toolkit", however, repeatedly identified sequences (or "building materials" on which the tools act) are present in a broad range of immune systems. These include a Toll/TLR system, a primitive complement system, an LPS binding protein, and a RAG core/Transib element. Repeatedly identified domains and motifs that function in immune proteins include NACHT, LRR, Ig, death, TIR, lectin domains, and a thioester motif. In addition, there are repeatedly identified mechanisms (or "construction methods") that generate sequence diversity in genes with immune function. These include genomic instability, duplications and/or deletions of sequences and the generation of clusters of similar genes or exons that appear as families, gene recombination, gene conversion, retrotransposition, alternative splicing, multiple alleles for single copy genes, and RNA editing. These commonly employed "materials and methods" for building and maintaining an effective immune system that might have been part of that ancestral system appear now as a fragmented and likely incomplete set, likely due to the rapid evolutionary change (or loss) of host genes that are under pressure to keep pace with pathogen diversity.

Predicting mutually exclusive spliced exons based on exon length, splice site and reading frame conservation, and exon sequence homology

BACKGROUND

Alternative splicing of pre-mature RNA is an important process eukaryotes utilize to increase their repertoire of different protein products. Several types of different alternative splice forms exist including exon skipping, differential splicing of exons at their 3'- or 5'-end, intron retention, and mutually exclusive splicing. The latter term is used for clusters of internal exons that are spliced in a mutually exclusive manner.

RESULTS

We have implemented an extension to the WebScipio software to search for mutually exclusive exons. Here, the search is based on the precondition that mutually exclusive exons encode regions of the same structural part of the protein product. This precondition provides restrictions to the search for candidate exons concerning their length, splice site conservation and reading frame preservation, and overall homology. Mutually exclusive exons that are not homologous and not of about the same length will not be found. Using the new algorithm, mutually exclusive exons in several example genes, a dynein heavy chain, a muscle myosin heavy chain, and Dscam were correctly identified. In addition, the algorithm was applied to the whole Drosophila melanogaster X chromosome and the results were compared to the Flybase annotation and an ab initio prediction. Clusters of mutually exclusive exons might be subsequent to each other and might encode dozens of exons.

CONCLUSIONS

This is the first implementation of an automatic search for mutually exclusive exons in eukaryotes. Exons are predicted and reconstructed in the same run providing the complete gene structure for the protein query of interest. WebScipio offers high quality gene structure figures with the clusters of mutually exclusive exons colour-coded, and several analysis tools for further manual inspection. The genome scale analysis of all genes of the Drosophila melanogaster X chromosome showed that WebScipio is able to find all but two of the 28 annotated mutually exclusive spliced exons and predicts 39 new candidate exons. Thus, WebScipio should be able to identify mutually exclusive spliced exons in any query sequence from any species with a very high probability. WebScipio is freely available to academics at http://www.webscipio.org.

Penaeus monodon Dscam (PmDscam) has a highly diverse cytoplasmic tail and is the first membrane-bound shrimp Dscam to be reported

Down syndrome cell adhesion molecule (Dscam) seems likely to play a key role in the "alternative adaptive immunity" that has been reported in invertebrates. Dscam consists of a cytoplasmic tail that is involved in signal transduction and a hypervariable extracellular region that might use a pathogen recognition mechanism similar to that used by the vertebrate antibodies. In our previous paper, we isolated a unique tail-less form of Dscam from Litopenaeus vannamei. In this study, we report the first membrane-bound form of shrimp Dscam: PmDscam was isolated from Penaeus monodon, and it occurred in both membrane-bound and tail-less forms. Phylogenetic analysis showed that while the crustacean Dscams from shrimp and water flea did not share a single subclade, they were distinct from the invertebrate Dscam-like molecules and from the insecta Dscams. In the extracellular region, the variable regions of PmDscam were located in N-terminal Ig2, N-terminal Ig3 and the entire Ig7 domain. The PmDscam extracellular variants and transmembrane domain variants were produced by mutually exclusive alternative splicing events. The cytoplasmic tail variants were produced by exon inclusion/exclusion. Based on the genomic organization of Daphnia Dscam's cytoplasmic tail, we propose a model of how the shrimp Dscam genomic locus might use Type III polyadenylation to generate both the tail-less and membrane-bound forms.

Competing RNA secondary structures are required for mutually exclusive splicing of the Dscam exon 6 cluster

Alternative splicing of eukaryotic pre-mRNAs is an important mechanism for generating proteome diversity and regulating gene expression. The Drosophila melanogaster Down Syndrome Cell Adhesion Molecule (Dscam) gene is an extreme example of mutually exclusive splicing. Dscam contains 95 alternatively spliced exons that potentially encode 38,016 distinct mRNA and protein isoforms. We previously identified two sets of conserved sequence elements, the docking site and selector sequences in the Dscam exon 6 cluster, which contains 48 mutually exclusive exons. These elements were proposed to engage in competing RNA secondary structures required for mutually exclusive splicing, though this model has not yet been experimentally tested. Here we describe a new system that allowed us to demonstrate that the docking site and selector sequences are indeed required for exon 6 mutually exclusive splicing and that the strength of these RNA structures determines the frequency of exon 6 inclusion. We also show that the function of the docking site has been conserved for ~500 million years of evolution. This work demonstrates that conserved intronic sequences play a functional role in mutually exclusive splicing of the Dscam exon 6 cluster.

RNA secondary structure in mutually exclusive splicing

Mutually exclusive splicing is a regulated means to generate protein diversity, but the underlying mechanisms are poorly understood. Here comparative genome analysis revealed the built-in intronic elements for controlling mutually exclusive splicing of the 14-3-3ξ pre-mRNA. These elements are clade specific but are evolutionarily conserved at the secondary structure level. Combined evidence revealed the triple functions of these inter-intronic RNA pairings in synergistically ensuring the selection of only one of multiple exons, through activation of the proximal variable exon outside the loop by the approximation of cis elements, and simultaneous repression of the exon within the loop, in combination with the physical competition of RNA pairing. Additionally, under this model, we also deciphered a similar structural code in exon clusters 4 and 9 of Dscam (38,016 isoforms) and Mhc (480 isoforms). Our findings suggest a broadly applicable mechanism to ensure mutually exclusive splicing.

Titin diversity--alternative splicing gone wild

Titin is an extremely large protein found in highest concentrations in heart and skeletal muscle. The single mammalian gene is expressed in multiple isoforms as a result of alternative splicing. Although titin isoform expression is controlled developmentally and in a tissue specific manner, the vast number of potential splicing pathways far exceeds those described in any other alternatively spliced gene. Over 1 million human splice pathways for a single individual can be potentially derived from the PEVK region alone. A new splicing pattern for the human cardiac N2BA isoform type has been found in which the PEVK region includes only the N2B type exons. The alterations in splicing and titin isoform expression in human heart disease provide impetus for future detailed study of the splicing mechanisms for this giant protein.