Large introns in relation to alternative splicing and gene evolution: a case study of Drosophila bruno-3

Background: Alternative splicing (AS) of maturing mRNA can generate structurally and functionally distinct transcripts from the same gene. Recent bioinformatic analyses of available genome databases inferred a positive correlation between intron length and AS. To study the interplay between intron length and AS empirically and in more detail, we analyzed the diversity of alternatively spliced transcripts (ASTs) in the Drosophila RNA-binding Bruno-3 (Bru-3) gene. This gene was known to encode thirteen exons separated by introns of diverse sizes, ranging from 71 to 41,973 nucleotides in D. melanogaster. Although Bru-3's structure is expected to be conducive to AS, only two ASTs of this gene were previously described. Results: Cloning of RT-PCR products of the entire ORF from four species representing three diverged Drosophila lineages provided an evolutionary perspective, high sensitivity, and long-range contiguity of splice choices currently unattainable by high-throughput methods. Consequently, we identified three new exons, a new exon fragment and thirty-three previously unknown ASTs of Bru-3. All exon-skipping events in the gene were mapped to the exons surrounded by introns of at least 800 nucleotides, whereas exons split by introns of less than 250 nucleotides were always spliced contiguously in mRNA. Cases of exon loss and creation during Bru-3 evolution in Drosophila were also localized within large introns. Notably, we identified a true de novo exon gain: exon 8 was created along the lineage of the obscura group from intronic sequence between cryptic splice sites conserved among all Drosophila species surveyed. Exon 8 was included in mature mRNA by the species representing all the major branches of the obscura group. To our knowledge, the origin of exon 8 is the first documented case of exonization of intronic sequence outside vertebrates. Conclusion: We found that large introns can promote AS via exon-skipping and exon turnover during evolution likely due to frequent errors in their removal from maturing mRNA. Large introns could be a reservoir of genetic diversity, because they have a greater number of mutable sites than short introns. Taken together, gene structure can constrain and/or promote gene evolution

Risks to Birds Traded for African Traditional Medicine: A Quantitative Assessment

Few regional or continent-wide assessments of bird use for traditional medicine have been attempted anywhere in the world. Africa has the highest known diversity of bird species used for this purpose. This study assesses the vulnerability of 354 bird species used for traditional medicine in 25 African countries, from 205 genera, 70 families, and 25 orders. The orders most represented were Passeriformes (107 species), Falconiformes (45 species), and Coraciiformes (24 species), and the families Accipitridae (37 species), Ardeidae (15 species), and Bucerotidae (12 species). The Barn owl (Tyto alba) was the most widely sold species (seven countries). The similarity of avifaunal orders traded is high (analogous to ‘‘morphospecies’’, and using Sørensen’s index), which suggests opportunities for a common understanding of cultural factors driving demand. The highest similarity was between bird orders sold in markets of Benin vs. Burkina Faso (90%), but even bird orders sold in two geographically separated countries (Benin vs. South Africa and Nigeria vs. South Africa) were 87% and 81% similar, respectively. Rabinowitz’s ‘‘7 forms of rarity’’ model, used to group species according to commonness or rarity, indicated that 24% of traded bird species are very common, locally abundant in several habitats, and occur over a large geographical area, but 10% are rare, occur in low numbers in specific habitats, and over a small geographical area. The order with the highest proportion of rare species was the Musophagiformes. An analysis of species mass (as a proxy for size) indicated that large and/or conspicuous species tend to be targeted by harvesters for the traditional medicine trade. Furthermore, based on cluster analyses for species groups of similar risk, vultures, hornbills, and other large avifauna, such as bustards, are most threatened by selective harvesting and should be prioritised for conservation action.University of the Witwatersrand SPARC Prestigious and URC Postdoctoral Fellowships; National Research Foundatio

RsaI repetitive DNA in Buffalo Bubalus bubalis representing retrotransposons, conserved in bovids, are part of the functional genes

<p>Abstract</p> <p>Background</p> <p>Repetitive sequences are the major components of the eukaryotic genomes. Association of these repeats with transcribing sequences and their regulation in buffalo <it>Bubalus bubalis </it>has remained largely unresolved.</p> <p>Results</p> <p>We cloned and sequenced <it>RsaI </it>repeat fragments pDp1, pDp2, pDp3, pDp4 of 1331, 651, 603 and 339 base pairs, respectively from the buffalo, <it>Bubalus bubalis</it>. Upon characterization, these fragments were found to represent retrotransposons and part of some functional genes. The resultant clones showed cross hybridization only with buffalo, cattle, goat and sheep genomic DNA. Real Time PCR, detected ~2 × 10⁴copies of pDp1, ~ 3000 copies of pDp2 and pDp3 and ~ 1000 of pDp4 in buffalo, cattle, goat and sheep genomes, respectively. <it>RsaI </it>repeats are transcriptionally active in somatic tissues and spermatozoa. Accordingly, pDp1 showed maximum expression in lung, pDp2 and pDp3 both in Kidney, and pDp4 in ovary. Fluorescence <it>in situ </it>hybridization showed repeats to be distributed all across the chromosomes.</p> <p>Conclusions</p> <p>The data suggest that <it>RsaI </it>repeats have been incorporated into the exonic regions of various transcribing genes, possibly contributing towards the architecture and evolution of the buffalo and related genomes. Prospects of our present work in the context of comparative and functional genomics are highlighted.</p

Identification of Rhoptry Trafficking Determinants and Evidence for a Novel Sorting Mechanism in the Malaria Parasite Plasmodium falciparum

The rhoptry of the malaria parasite Plasmodium falciparum is an unusual secretory organelle that is thought to be related to secretory lysosomes in higher eukaryotes. Rhoptries contain an extensive collection of proteins that participate in host cell invasion and in the formation of the parasitophorous vacuole, but little is known about sorting signals required for rhoptry protein targeting. Using green fluorescent protein chimeras and in vitro pull-down assays, we performed an analysis of the signals required for trafficking of the rhoptry protein RAP1. We provide evidence that RAP1 is escorted to the rhoptry via an interaction with the glycosylphosphatidyl inositol-anchored rhoptry protein RAMA. Once within the rhoptry, RAP1 contains distinct signals for localisation within a sub-compartment of the organelle and subsequent transfer to the parasitophorous vacuole after invasion. This is the first detailed description of rhoptry trafficking signals in Plasmodium

Alu Exonization Events Reveal Features Required for Precise Recognition of Exons by the Splicing Machinery

Despite decades of research, the question of how the mRNA splicing machinery precisely identifies short exonic islands within the vast intronic oceans remains to a large extent obscure. In this study, we analyzed Alu exonization events, aiming to understand the requirements for correct selection of exons. Comparison of exonizing Alus to their non-exonizing counterparts is informative because Alus in these two groups have retained high sequence similarity but are perceived differently by the splicing machinery. We identified and characterized numerous features used by the splicing machinery to discriminate between Alu exons and their non-exonizing counterparts. Of these, the most novel is secondary structure: Alu exons in general and their 5′ splice sites (5′ss) in particular are characterized by decreased stability of local secondary structures with respect to their non-exonizing counterparts. We detected numerous further differences between Alu exons and their non-exonizing counterparts, among others in terms of exon–intron architecture and strength of splicing signals, enhancers, and silencers. Support vector machine analysis revealed that these features allow a high level of discrimination (AUC = 0.91) between exonizing and non-exonizing Alus. Moreover, the computationally derived probabilities of exonization significantly correlated with the biological inclusion level of the Alu exons, and the model could also be extended to general datasets of constitutive and alternative exons. This indicates that the features detected and explored in this study provide the basis not only for precise exon selection but also for the fine-tuned regulation thereof, manifested in cases of alternative splicing

The Origin Recognition Complex Interacts with a Subset of Metabolic Genes Tightly Linked to Origins of Replication

The origin recognition complex (ORC) marks chromosomal sites as replication origins and is essential for replication initiation. In yeast, ORC also binds to DNA elements called silencers, where its primary function is to recruit silent information regulator (SIR) proteins to establish transcriptional silencing. Indeed, silencers function poorly as chromosomal origins. Several genetic, molecular, and biochemical studies of HMR-E have led to a model proposing that when ORC becomes limiting in the cell (such as in the orc2-1 mutant) only sites that bind ORC tightly (such as HMR-E) remain fully occupied by ORC, while lower affinity sites, including many origins, lose ORC occupancy. Since HMR-E possessed a unique non-replication function, we reasoned that other tight sites might reveal novel functions for ORC on chromosomes. Therefore, we comprehensively determined ORC “affinity” genome-wide by performing an ORC ChIP–on–chip in ORC2 and orc2-1 strains. Here we describe a novel group of orc2-1–resistant ORC–interacting chromosomal sites (ORF–ORC sites) that did not function as replication origins or silencers. Instead, ORF–ORC sites were comprised of protein-coding regions of highly transcribed metabolic genes. In contrast to the ORC–silencer paradigm, transcriptional activation promoted ORC association with these genes. Remarkably, ORF–ORC genes were enriched in proximity to origins of replication and, in several instances, were transcriptionally regulated by these origins. Taken together, these results suggest a surprising connection among ORC, replication origins, and cellular metabolism

Evolution of Alternative Splicing Regulation: Changes in Predicted Exonic Splicing Regulators Are Not Associated with Changes in Alternative Splicing Levels in Primates

Alternative splicing is tightly regulated in a spatio-temporal and quantitative manner. This regulation is achieved by a complex interplay between spliceosomal (trans) factors that bind to different sequence (cis) elements. cis-elements reside in both introns and exons and may either enhance or silence splicing. Differential combinations of cis-elements allows for a huge diversity of overall splicing signals, together comprising a complex ‘splicing code’. Many cis-elements have been identified, and their effects on exon inclusion levels demonstrated in reporter systems. However, the impact of interspecific differences in these elements on the evolution of alternative splicing levels has not yet been investigated at genomic level. Here we study the effect of interspecific differences in predicted exonic splicing regulators (ESRs) on exon inclusion levels in human and chimpanzee. For this purpose, we compiled and studied comprehensive datasets of predicted ESRs, identified by several computational and experimental approaches, as well as microarray data for changes in alternative splicing levels between human and chimpanzee. Surprisingly, we found no association between changes in predicted ESRs and changes in alternative splicing levels. This observation holds across different ESR exon positions, exon lengths, and 5′ splice site strengths. We suggest that this lack of association is mainly due to the great importance of context for ESR functionality: many ESR-like motifs in primates may have little or no effect on splicing, and thus interspecific changes at short-time scales may primarily occur in these effectively neutral ESRs. These results underscore the difficulties of using current computational ESR prediction algorithms to identify truly functionally important motifs, and provide a cautionary tale for studies of the effect of SNPs on splicing in human disease

A Critical Analysis of Atoh7 (Math5) mRNA Splicing in the Developing Mouse Retina

The Math5 (Atoh7) gene is transiently expressed during retinogenesis by progenitors exiting mitosis, and is essential for ganglion cell (RGC) development. Math5 contains a single exon, and its 1.7 kb mRNA encodes a 149-aa polypeptide. Mouse Math5 mutants have essentially no RGCs or optic nerves. Given the importance of this gene in retinal development, we thoroughly investigated the possibility of Math5 mRNA splicing by Northern blot, 3′RACE, RNase protection assays, and RT-PCR, using RNAs extracted from embryonic eyes and adult cerebellum, or transcribed in vitro from cDNA clones. Because Math5 mRNA contains an elevated G+C content, we used graded concentrations of betaine, an isostabilizing agent that disrupts secondary structure. Although ∼10% of cerebellar Math5 RNAs are spliced, truncating the polypeptide, our results show few, if any, spliced Math5 transcripts exist in the developing retina (<1%). Rare deleted cDNAs do arise via RT-mediated RNA template switching in vitro, and are selectively amplified during PCR. These data differ starkly from a recent study (Kanadia and Cepko 2010), which concluded that the vast majority of Math5 and other bHLH transcripts are spliced to generate noncoding RNAs. Our findings clarify the architecture of the Math5 gene and its mechanism of action. These results have implications for all members of the bHLH gene family, for any gene that is alternatively spliced, and for the interpretation of all RT-PCR experiments

Interaction between Purkinje Cells and Inhibitory Interneurons May Create Adjustable Output Waveforms to Generate Timed Cerebellar Output

We develop a new model that explains how the cerebellum may generate the timing in classical delay eyeblink conditioning. Recent studies show that both Purkinje cells (PCs) and inhibitory interneurons (INs) have parallel signal processing streams with two time scales: an AMPA receptor-mediated fast process and a metabotropic glutamate receptor (mGluR)-mediated slow process. Moreover, one consistent finding is an increased excitability of PC dendrites (in Larsell's lobule HVI) in animals when they acquire the classical delay eyeblink conditioning naturally, in contrast to in vitro studies, where learning involves long-term depression (LTD). Our model proposes that the delayed response comes from the slow dynamics of mGluR-mediated IP3 activation, and the ensuing calcium concentration change, and not from LTP/LTD. The conditioned stimulus (tone), arriving on the parallel fibers, triggers this slow activation in INs and PC spines. These excitatory (from PC spines) and inhibitory (from INs) signals then interact at the PC dendrites to generate variable waveforms of PC activation. When the unconditioned stimulus (puff), arriving on the climbing fibers, is coupled frequently with this slow activation the waveform is amplified (due to an increased excitability) and leads to a timed pause in the PC population. The disinhibition of deep cerebellar nuclei by this timed pause causes the delayed conditioned response. This suggested PC-IN interaction emphasizes a richer role of the INs in learning and also conforms to the recent evidence that mGluR in the cerebellar cortex may participate in slow motor execution. We show that the suggested mechanism can endow the cerebellar cortex with the versatility to learn almost any temporal pattern, in addition to those that arise in classical conditioning