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

Exclusive Leptoproduction of rho^0 Mesons from Hydrogen at Intermediate Virtual Photon Energies

Measurements of the cross section for exclusive virtual-photoproduction of rho^0 mesons from hydrogen are reported. The data were collected by the HERMES experiment using 27.5 GeV positrons incident on a hydrogen gas target in the HERA storage ring. The invariant mass W of the photon-nucleon system ranges from 4.0 to 6.0 GeV, while the negative squared four-momentum Q^2 of the virtual photon varies from 0.7 to 5.0 GeV^2. The present data together with most of the previous data at W > 4 GeV are well described by a model that infers the W-dependence of the cross section from the dependence on the Bjorken scaling variable x of the unpolarized structure function for deep-inelastic scattering. In addition, a model calculation based on Off-Forward Parton Distributions gives a fairly good account of the longitudinal component of the rho^0 production cross section for Q^2 > 2 GeV^2.Comment: 10 pages, 6 embedded figures, LaTeX for SVJour(epj) document class. Revisions: curves added to Fig. 1, several clarifications added to tex

Using C. elegans to decipher the cellular and molecular mechanisms underlying neurodevelopmental disorders

Prova tipográfica (uncorrected proof)Neurodevelopmental disorders such as epilepsy, intellectual disability (ID), and autism spectrum disorders (ASDs) occur in over 2 % of the population, as the result of genetic mutations, environmental factors, or combination of both. In the last years, use of large-scale genomic techniques allowed important advances in the identification of genes/loci associated with these disorders. Nevertheless, following association of novel genes with a given disease, interpretation of findings is often difficult due to lack of information on gene function and effect of a given mutation in the corresponding protein. This brings the need to validate genetic associations from a functional perspective in model systems in a relatively fast but effective manner. In this context, the small nematode, Caenorhabditis elegans, presents a good compromise between the simplicity of cell models and the complexity of rodent nervous systems. In this article, we review the features that make C. elegans a good model for the study of neurodevelopmental diseases. We discuss its nervous system architecture and function as well as the molecular basis of behaviors that seem important in the context of different neurodevelopmental disorders. We review methodologies used to assess memory, learning, and social behavior as well as susceptibility to seizures in this organism. We will also discuss technological progresses applied in C. elegans neurobiology research, such as use of microfluidics and optogenetic tools. Finally, we will present some interesting examples of the functional analysis of genes associated with human neurodevelopmental disorders and how we can move from genes to therapies using this simple model organism.The authors would like to acknowledge Fundação para a Ciência e Tecnologia (FCT) (PTDC/SAU-GMG/112577/2009). AJR and CB are recipients of FCT fellowships: SFRH/BPD/33611/2009 and SFRH/BPD/74452/2010, respectively

Expedited batch processing and analysis of transposon insertions

<p>Abstract</p> <p>Background</p> <p>With advances in sequencing technology, greater and greater amounts of eukaryotic genome data are becoming available. Often, large portions of these genomes consist of transposable elements, frequently accounting for 50% or more in vertebrates. Each transposable element family may have thousands or tens of thousands of individual copies within a given genome, and therefore it can take an exorbitant amount of time and effort to process data in a meaningful fashion.</p> <p>Findings</p> <p>In order to combat this problem, we developed a set of bioinformatics techniques and programs to streamline the analysis. This includes a unique Perl script which automates the process of taking BLAST, Repeatmasker and similar data to extract and manipulate the hit sequences from the genome. This script, called Process_hits uses an object-oriented methodology to compile all hit locations from a given file for processing, organize this data into useable categories, and output it in multiple formats.</p> <p>Conclusions</p> <p>The program proved capable of handling large amounts of transposon data in an efficient fashion. It is equipped with a number of useful sub-functions, each of which is contained within its own sub-module to allow for greater expandability and as a foundation for future program design.</p

Reliability analysis of the Ahringer Caenorhabditis elegans RNAi feeding library: a guide for genome-wide screens

<p>Abstract</p> <p>Background</p> <p>The Ahringer <it>C. elegans </it>RNAi feeding library prepared by cloning genomic DNA fragments has been widely used in genome-wide analysis of gene function. However, the library has not been thoroughly validated by direct sequencing, and there are potential errors, including: 1) mis-annotation (the clone with the retired gene name should be remapped to the actual target gene); 2) nonspecific PCR amplification; 3) cross-RNAi; 4) mis-operation such as sample loading error, <it>etc</it>.</p> <p>Results</p> <p>Here we performed a reliability analysis on the Ahringer <it>C. elegans </it>RNAi feeding library, which contains 16,256 bacterial strains, using a bioinformatics approach. Results demonstrated that most (98.3%) of the bacterial strains in the library are reliable. However, we also found that 2,851 (17.54%) bacterial strains need to be re-annotated even they are reliable. Most of these bacterial strains are the clones having the retired gene names. Besides, 28 strains are grouped into unreliable category and 226 strains are marginal because of probably expressing unrelated double-stranded RNAs (dsRNAs). The accuracy of the prediction was further confirmed by direct sequencing analysis of 496 bacterial strains. Finally, a freely accessible database named CelRNAi (<url>http://biocompute.bmi.ac.cn/CelRNAi/</url>) was developed as a valuable complement resource for the feeding RNAi library by providing the predicted information on all bacterial strains. Moreover, submission of the direct sequencing result or any other annotations for the bacterial strains to the database are allowed and will be integrated into the CelRNAi database to improve the accuracy of the library. In addition, we provide five candidate primer sets for each of the unreliable and marginal bacterial strains for users to construct an alternative vector for their own RNAi studies.</p> <p>Conclusions</p> <p>Because of the potential unreliability of the Ahringer <it>C. elegans </it>RNAi feeding library, we strongly suggest the user examine the reliability information of the bacterial strains in the CelRNAi database before performing RNAi experiments, as well as the post-RNAi experiment analysis.</p

Linkage mapping bovine EST-based SNP

BACKGROUND: Existing linkage maps of the bovine genome primarily contain anonymous microsatellite markers. These maps have proved valuable for mapping quantitative trait loci (QTL) to broad regions of the genome, but more closely spaced markers are needed to fine-map QTL, and markers associated with genes and annotated sequence are needed to identify genes and sequence variation that may explain QTL. RESULTS: Bovine expressed sequence tag (EST) and bacterial artificial chromosome (BAC)sequence data were used to develop 918 single nucleotide polymorphism (SNP) markers to map genes on the bovine linkage map. DNA of sires from the MARC reference population was used to detect SNPs, and progeny and mates of heterozygous sires were genotyped. Chromosome assignments for 861 SNPs were determined by twopoint analysis, and positions for 735 SNPs were established by multipoint analyses. Linkage maps of bovine autosomes with these SNPs represent 4585 markers in 2475 positions spanning 3058 cM . Markers include 3612 microsatellites, 913 SNPs and 60 other markers. Mean separation between marker positions is 1.2 cM. New SNP markers appear in 511 positions, with mean separation of 4.7 cM. Multi-allelic markers, mostly microsatellites, had a mean (maximum) of 216 (366) informative meioses, and a mean 3-lod confidence interval of 3.6 cM Bi-allelic markers, including SNP and other marker types, had a mean (maximum) of 55 (191) informative meioses, and were placed within a mean 8.5 cM 3-lod confidence interval. Homologous human sequences were identified for 1159 markers, including 582 newly developed and mapped SNP. CONCLUSION: Addition of these EST- and BAC-based SNPs to the bovine linkage map not only increases marker density, but provides connections to gene-rich physical maps, including annotated human sequence. The map provides a resource for fine-mapping quantitative trait loci and identification of positional candidate genes, and can be integrated with other data to guide and refine assembly of bovine genome sequence. Even after the bovine genome is completely sequenced, the map will continue to be a useful tool to link observable phenotypes and animal genotypes to underlying genes and molecular mechanisms influencing economically important beef and dairy traits

Molecular diversity of phospholipase D in angiosperms

BACKGROUND: The phospholipase D (PLD) family has been identified in plants by recent molecular studies, fostered by the emerging importance of plant PLDs in stress physiology and signal transduction. However, the presence of multiple isoforms limits the power of conventional biochemical and pharmacological approaches, and calls for a wider application of genetic methodology. RESULTS: Taking advantage of sequence data available in public databases, we attempted to provide a prerequisite for such an approach. We made a complete inventory of the Arabidopsis thaliana PLD family, which was found to comprise 12 distinct genes. The current nomenclature of Arabidopsis PLDs was refined and expanded to include five newly described genes. To assess the degree of plant PLD diversity beyond Arabidopsis we explored data from rice (including the genome draft by Monsanto) as well as cDNA and EST sequences from several other plants. Our analysis revealed two major PLD subfamilies in plants. The first, designated C2-PLD, is characterised by presence of the C2 domain and comprises previously known plant PLDs as well as new isoforms with possibly unusual features-catalytically inactive or independent on Ca(2+). The second subfamily (denoted PXPH-PLD) is novel in plants but is related to animal and fungal enzymes possessing the PX and PH domains. CONCLUSIONS: The evolutionary dynamics, and inter-specific diversity, of plant PLDs inferred from our phylogenetic analysis, call for more plant species to be employed in PLD research. This will enable us to obtain generally valid conclusions

Linking microarray reporters with protein functions

<p>Abstract</p> <p>Background</p> <p>The analysis of microarray experiments requires accurate and up-to-date functional annotation of the microarray reporters to optimize the interpretation of the biological processes involved. Pathway visualization tools are used to connect gene expression data with existing biological pathways by using specific database identifiers that link reporters with elements in the pathways.</p> <p>Results</p> <p>This paper proposes a novel method that aims to improve microarray reporter annotation by BLASTing the original reporter sequences against a species-specific EMBL subset, that was derived from and crosslinked back to the highly curated UniProt database. The resulting alignments were filtered using high quality alignment criteria and further compared with the outcome of a more traditional approach, where reporter sequences were BLASTed against EnsEMBL followed by locating the corresponding protein (UniProt) entry for the high quality hits. Combining the results of both methods resulted in successful annotation of > 58% of all reporter sequences with UniProt IDs on two commercial array platforms, increasing the amount of Incyte reporters that could be coupled to Gene Ontology terms from 32.7% to 58.3% and to a local GenMAPP pathway from 9.6% to 16.7%. For Agilent, 35.3% of the total reporters are now linked towards GO nodes and 7.1% on local pathways.</p> <p>Conclusion</p> <p>Our methods increased the annotation quality of microarray reporter sequences and allowed us to visualize more reporters using pathway visualization tools. Even in cases where the original reporter annotation showed the correct description the new identifiers often allowed improved pathway and Gene Ontology linking. These methods are freely available at http://www.bigcat.unimaas.nl/public/publications/Gaj_Annotation/.</p

Preservation of Genes Involved in Sterol Metabolism in Cholesterol Auxotrophs: Facts and Hypotheses

Background: It is known that primary sequences of enzymes involved in sterol biosynthesis are well conserved in organisms that produce sterols de novo. However, we provide evidence for a preservation of the corresponding genes in two animals unable to synthesize cholesterol (auxotrophs): Drosophila melanogaster and Caenorhabditis elegans. Principal Findings: We have been able to detect bona fide orthologs of several ERG genes in both organisms using a series of complementary approaches. We have detected strong sequence divergence between the orthologs of the nematode and of the fruitfly; they are also very divergent with respect to the orthologs in organisms able to synthesize sterols de novo (prototrophs). Interestingly, the orthologs in both the nematode and the fruitfly are still under selective pressure. It is possible that these genes, which are not involved in cholesterol synthesis anymore, have been recruited to perform different new functions. We propose a more parsimonious way to explain their accelerated evolution and subsequent stabilization. The products of ERG genes in prototrophs might be involved in several biological roles, in addition to sterol synthesis. In the case of the nematode and the fruitfly, the relevant genes would have lost their ancestral function in cholesterogenesis but would have retained the other function(s), which keep them under pressure. Conclusions: By exploiting microarray data we have noticed a strong expressional correlation between the orthologs of ERG24 and ERG25 in D. melanogaster and genes encoding factors involved in intracellular protein trafficking and folding an

Effect of valdecoxib pretreatment on pain and secondary hyperalgesia: a randomized controlled trial in healthy volunteers [ISRCTN05282752, NCT00260325]

BACKGROUND: Induction of the COX-2 isoenzyme appears to play a major role in the genesis of central sensitization after nociceptive stimulation. This study aimed to investigate the efficacy of a single, oral dose of the specific COX-2 inhibitor-valdecoxib in attenuating the central sensitization – induced secondary hyperalgesia in a heat/capsaicin pain model in healthy volunteers. METHODS: The study was a randomized, double blind, placebo controlled, crossover, single dose efficacy trial using 20 healthy volunteers. Two hours following placebo or 40 mg, PO valdecoxib, participants underwent skin sensitization with heat/capsaicin, as well as supra-threshold pain and re-kindling measurements according to an established, validated pain model. Subjects rated pain intensity and unpleasantness on a visual analog scale and the area of secondary hyperalgesia was serially mapped. RESULTS: The area of secondary hyperalgesia produced after 40 mg of valdecoxib was no different than that after placebo. Furthermore, there were no significantly relevant differences when volunteers were treated with valdecoxib or placebo in relation to either cold- or hot pain threshold or the intensity of pain after supra-threshold, thermal pain stimulation. CONCLUSION: We demonstrated that a single, oral dose of valdecoxib when does not attenuate secondary hyperalgesia induced by heat/capsaicin in a cutaneous sensitization pain model in healthy volunteers