Identification of Conserved Splicing Motifs in Mutually Exclusive Exons of 15 Insect Species

Background: During alternative splicing, the inclusion of an exon in the final mRNA molecule is determined by nuclear proteins that bind cis-regulatory sequences in a target pre-mRNA molecule. A recent study suggested that the regulatory codes of individual RNA-binding proteins may be nearly immutable between very diverse species such as mammals and insects. The model system Drosophila melanogaster therefore presents an excellent opportunity for the study of alternative splicing due to the availability of quality EST annotations in FlyBase. Methods: In this paper, we describe an in silico analysis pipeline to extract putative exonic splicing regulatory sequences from a multiple alignment of 15 species of insects. Our method, ESTs-to-ESRs (E2E), uses graph analysis of EST splicing graphs to identify mutually exclusive (ME) exons and combines phylogenetic measures, a sliding window approach along the multiple alignment and the Welch’s t statistic to extract conserved ESR motifs. Results: The most frequent 100% conserved word of length 5 bp in different insect exons was “ATGGA”. We identified 799 statistically significant “spike” hexamers, 218 motifs with either a left or right FDR corrected spike magnitude p-value < 0.05 and 83 with both left and right uncorrected p < 0.01. 11 genes were identified with highly significant motifs in one ME exon but not in the other, suggesting regulation of ME exon splicing through these highly conserved hexamers. The majority of these genes have been shown to have regulated spatiotemporal expression. 10 elements were found to match three mammalian splicing regulator databases. A putative ESR motif, GATGCAG, was identified in the ME-13b but not in the ME-13a of Drosophila N-Cadherin, a gene that has been shown to have a distinct spatiotemporal expression pattern of spliced isoforms in a recent study. Conclusions: Analysis of phylogenetic relationships and variability of sequence conservation as implemented in the E2E spikes method may lead to improved identification of ESRs. We found that approximately half of the putative ESRs in common between insects and mammals have a high statistical support (p < 0.01). Several Drosophila genes with spatiotemporal expression patterns were identified to contain putative ESRs located in one exon of the ME exon pairs but not in the other.

Simulating Mechanisms for Dispersal, Production and Stranding of Small Forage Fish in Temporary Wetland Habitats

Movement strategies of small forage fish (<8 cm total length) between temporary and permanent wetland habitats affect their overall population growth and biomass concentrations, i.e., availability to predators. These fish are often the key energy link between primary producers and top predators, such as wading birds, which require high concentrations of stranded fish in accessible depths. Expansion and contraction of seasonal wetlands induce a sequential alternation between rapid biomass growth and concentration, creating the conditions for local stranding of small fish as they move in response to varying water levels. To better understand how landscape topography, hydrology, and fish behavior interact to create high densities of stranded fish, we first simulated population dynamics of small fish, within a dynamic food web, with different traits for movement strategy and growth rate, across an artificial, spatially explicit, heterogeneous, two-dimensional marsh slough landscape, using hydrologic variability as the driver for movement. Model output showed that fish with the highest tendency to invade newly flooded marsh areas built up the largest populations over long time periods with stable hydrologic patterns. A higher probability to become stranded had negative effects on long-term population size, and offset the contribution of that species to stranded biomass. The model was next applied to the topography of a 10 km × 10 km area of Everglades landscape. The details of the topography were highly important in channeling fish movements and creating spatiotemporal patterns of fish movement and stranding. This output provides data that can be compared in the future with observed locations of fish biomass concentrations, or such surrogates as phosphorus ‘hotspots’ in the marsh.