AT-AC pre-mRNA splicing mechanisms and conservation of minor introns in voltage-gated ion channel genes

Deciphering the impact of coding and non-coding SCN1A gene variants on RNA splicing

Variants that disrupt normal pre-mRNA splicing are increasingly being recognized as a major cause of monogenic disorders. The SCN1A gene, a key epilepsy gene that is linked to various epilepsy phenotypes, is no exception. Approximately 10% of all reported variants in the SCN1A gene are designated as splicing variants, with many located outside of the canonical donor and acceptor splice sites, and most have not been functionally investigated. However, given its restricted expression pattern, functional analysis of splicing variants in the SCN1A gene could not be routinely performed. In this study, we conducted a comprehensive analysis of all reported SCN1A variants and their potential to impact SCN1A splicing and conclude that splicing variants are substantially misannotated and under-represented. We created a splicing reporter system consisting of 18 splicing vectors covering all 26 protein-coding exons with different genomic contexts and several promoters of varying strengths in order to reproduce the wild-type splicing pattern of the SCN1A gene, revealing cis-regulatory elements essential for proper recognition of SCN1A exons. Functional analysis of 95 SCN1A variants was carried out, including all 68 intronic variants reported in the literature, located outside of the splice sites canonical dinucleotides; 21 exonic variants of different classes (synonymous, missense, nonsense and in-frame deletion) and six variants observed in patients with epilepsy. Interestingly, almost 20% of tested intronic variants had no influence on SCN1A splicing, despite being reported as causative in the literature. Moreover, we confirmed that the majority of predicted exonic variants affect splicing unravelling their true molecular mechanism. We used functional data to perform genotype-phenotype correlation, revealing distinct distribution patterns for missense and splice-affecting 'missense' variants and observed no difference in the phenotype severity of variants leading to in-frame and out-of-frame isoforms, indicating that the Nav1.1 protein is highly intolerant to structural variations. Our work demonstrates the importance of functional analysis in proper variant annotation and provides a tool for high-throughput delineation of splice-affecting variants in SCN1A in a whole-gene manner.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

HPV16 and HPV18 Genome Structure, Expression, and Post-Transcriptional Regulation

Human papillomaviruses (HPV) are a group of small non-enveloped DNA viruses whose infection causes benign tumors or cancers. HPV16 and HPV18, the two most common high-risk HPVs, are responsible for ~70% of all HPV-related cervical cancers and head and neck cancers. The expression of the HPV genome is highly dependent on cell differentiation and is strictly regulated at the transcriptional and post-transcriptional levels. Both HPV early and late transcripts differentially expressed in the infected cells are intron-containing bicistronic or polycistronic RNAs bearing more than one open reading frame (ORF), because of usage of alternative viral promoters and two alternative viral RNA polyadenylation signals. Papillomaviruses proficiently engage alternative RNA splicing to express individual ORFs from the bicistronic or polycistronic RNA transcripts. In this review, we discuss the genome structures and the updated transcription maps of HPV16 and HPV18, and the latest research advances in understanding RNA cis-elements, intron branch point sequences, and RNA-binding proteins in the regulation of viral RNA processing. Moreover, we briefly discuss the epigenetic modifications, including DNA methylation and possible APOBEC-mediated genome editing in HPV infections and carcinogenesis.

Non-canonical splice junction processing increases the diversity of RBFOX2 splicing isoforms

The underlying mechanisms of splicing regulation through non-canonical splice junction processing remain largely unknown. Here, we identified two RBFOX2 splicing isoforms by alternative 3' splice site selection of exon 9; the non-canonical splice junction processed RBFOX2 transcript (RBFOX2-N.C.) was expressed by the selection of the 3' splice GG acceptor sequence. The cytoplasmic localization of RBFOX2-C., a canonical splice junction-processed RBFOX2 transcript, was different from that of RBFOX2-N.C., which showed nuclear localization. In addition, we confirmed that RBFOX2-C. showed a significantly stronger localization into stress granules than RBFOX2-N.C. upon sodium arsenite treatment. Next, we investigated the importance of non-canonical 3' splice GG sequence selection of specific cis-regulatory elements using minigene constructs of the RBFOX2 gene. We found that the non-canonical 3' splice GG sequence and suboptimal branch point site adjacent region were critical for RBFOX2-N.C. expression through a non-canonical 3' splice selection. Our results suggest a regulatory mechanism for the non-canonical 3' splice selection in the RBFOX2 gene, providing a basis for studies related to the regulation of alternative pre-mRNA splicing through non-canonical splice junction processing.

Inhibition of minor intron splicing reduces Na+ and Ca2+ channel expression and function in cardiomyocytes

Eukaryotic genomes contain a tiny subset of ‘minor class’ introns with unique sequence elements that require their own splicing machinery. These minor introns are present in certain gene families with specific functions, such as voltage-gated Na⁺ and voltage-gated Ca²⁺ channels. Removal of minor introns by the minor spliceosome has been proposed as a post-transcriptional regulatory layer, which remains unexplored in the heart. Here, we investigate whether the minor spliceosome regulates electrophysiological properties of cardiomyocytes by knocking down the essential minor spliceosome small nuclear snRNA component U6atac in neonatal rat ventricular myocytes. Loss of U6atac led to robust minor intron retention within Scn5a and Cacna1c, resulting in reduced protein levels of Na_v1.5 and Ca_v1.2 channels. Functional consequences were studied through patch-clamp analysis, and revealed reduced Na⁺ and L-type Ca²⁺ currents after loss of U6atac. In conclusion, minor intron splicing modulates voltage-dependent ion channel expression and function in cardiomyocytes. This may be of particular relevance in situations in which minor splicing activity changes, such as in genetic diseases affecting minor spliceosome components, or in acquired diseases in which minor spliceosome components are dysregulated, such as heart failure.

Summary: Knockdown of minor spliceosome component U6atac in cardiomyocytes reveals that expression of the Na⁺ channel Scn5a and the L-type Ca²⁺ channel Cacna1c critically depend on minor intron splicing.

Minor Intron Splicing from Basic Science to Disease

Pre-mRNA splicing is an essential step in gene expression and is catalyzed by two machineries in eukaryotes: the major (U2 type) and minor (U12 type) spliceosomes. While the majority of introns in humans are U2 type, less than 0.4% are U12 type, also known as minor introns (mi-INTs), and require a specialized spliceosome composed of U11, U12, U4atac, U5, and U6atac snRNPs. The high evolutionary conservation and apparent splicing inefficiency of U12 introns have set them apart from their major counterparts and led to speculations on the purpose for their existence. However, recent studies challenged the simple concept of mi-INTs splicing inefficiency due to low abundance of their spliceosome and confirmed their regulatory role in alternative splicing, significantly impacting the expression of their host genes. Additionally, a growing list of minor spliceosome-associated diseases with tissue-specific pathologies affirmed the importance of minor splicing as a key regulatory pathway, which when deregulated could lead to tissue-specific pathologies due to specific alterations in the expression of some minor-intron-containing genes. Consequently, uncovering how mi-INTs splicing is regulated in a tissue-specific manner would allow for better understanding of disease pathogenesis and pave the way for novel therapies, which we highlight in this review.

X-box binding protein 1-mediated COL4A1s secretion regulates communication between vascular smooth muscle and stem/progenitor cells

Vascular smooth muscle cells (VSMCs) contribute to the deposition of extracellular matrix proteins (ECMs), including Type IV collagen, in the vessel wall. ECMs coordinate communication among different cell types, but mechanisms underlying this communication remain unclear. Our previous studies have demonstrated that X-box binding protein 1 (XBP1) is activated and contributes to VSMC phenotypic transition in response to vascular injury. In this study, we investigated the participation of XBP1 in the communication between VSMCs and vascular progenitor cells (VPCs). Immunofluorescence and immunohistology staining revealed that Xbp1 gene was essential for type IV collagen alpha 1 (COL4A1) expression during mouse embryonic development and vessel wall ECM deposition and stem cell antigen 1-positive (Sca1⁺)-VPC recruitment in response to vascular injury. The Western blot analysis elucidated an Xbp1 gene dose-dependent effect on COL4A1 expression and that the spliced XBP1 protein (XBP1s) increased protease-mediated COL4A1 degradation as revealed by Zymography. RT-PCR analysis revealed that XBP1s in VSMCs not only upregulated COL4A1/2 transcription but also induced the occurrence of a novel transcript variant, soluble type IV collagen alpha 1 (COL4A1s), in which the front part of exon 4 is joined with the rear part of exon 42. Chromatin-immunoprecipitation, DNA/protein pulldown and in vitro transcription demonstrated that XBP1s binds to exon 4 and exon 42, directing the transcription from exon 4 to exon 42. This leads to transcription complex bypassing the internal sequences, producing a shortened COL4A1s protein that increased Sca1⁺-VPC migration. Taken together, these results suggest that activated VSMCs may recruit Sca1⁺-VPCs via XBP1s-mediated COL4A1s secretion, leading to vascular injury repair or neointima formation.

Minor snRNA gene delivery improves the loss of proprioceptive synapses on SMA motor neurons

Spinal muscular atrophy (SMA) is an inherited neuromuscular disorder caused by reduced expression of the survival motor neuron (SMN) protein. SMN has key functions in multiple RNA pathways, including the biogenesis of small nuclear ribonucleoproteins that are essential components of both major (U2-dependent) and minor (U12-dependent) spliceosomes. Here we investigated the specific contribution of U12 splicing dysfunction to SMA pathology through selective restoration of this RNA pathway in mouse models of varying phenotypic severity. We show that virus-mediated delivery of minor snRNA genes specifically improves select U12 splicing defects induced by SMN deficiency in cultured mammalian cells, as well as in the spinal cord and dorsal root ganglia of SMA mice without increasing SMN expression. This approach resulted in a moderate amelioration of several parameters of the disease phenotype in SMA mice, including survival, weight gain, and motor function. Importantly, minor snRNA gene delivery improved aberrant splicing of the U12 intron-containing gene Stasimon and rescued the severe loss of proprioceptive sensory synapses on SMA motor neurons, which are early signatures of motor circuit dysfunction in mouse models. Taken together, these findings establish the direct contribution of U12 splicing dysfunction to synaptic deafferentation and motor circuit pathology in SMA.

An Integrated Model of Minor Intron Emergence and Conservation

Minor introns constitute <0.5% of the introns in the human genome and have remained an enigma since their discovery. These introns are removed by a distinct splicing complex, the minor spliceosome. Both are ancient, tracing back to the last eukaryotic common ancestor (LECA), which is reflected by minor intron enrichment in specific gene families, such as the mitogen activated-protein kinase kinases, voltage-gated sodium and calcium ion channels, and E2F transcription factors. Most minor introns occur as single introns in genes with predominantly major introns. Due to this organization, minor intron-containing gene (MIG) expression requires the coordinated action of two spliceosomes, which increases the probability of missplicing. Thus, one would expect loss of minor introns via purifying selection. This has resulted in complete minor intron loss in at least nine eukaryotic lineages. However, minor introns are highly conserved in land plants and metazoans, where their importance is underscored by embryonic lethality when the minor spliceosome is inactivated. Conditional inactivation of the minor spliceosome has shown that rapidly dividing progenitor cells are highly sensitive to minor spliceosome loss. Indeed, we found that MIGs were significantly enriched in a screen for genes essential for survival in 341 cycling cell lines. Here, we propose that minor introns inserted randomly into genes in LECA or earlier and were subsequently conserved in genes crucial for cycling cell survival. We hypothesize that the essentiality of MIGs allowed minor introns to endure through the unicellularity of early eukaryotic evolution. Moreover, we identified 59 MIGs that emerged after LECA, and that many of these are essential for cycling cell survival, reinforcing our essentiality model for MIG conservation. This suggests that minor intron emergence is dynamic across eukaryotic evolution, and that minor introns should not be viewed as molecular fossils. We also posit that minor intron splicing was co-opted in multicellular evolution as a regulatory switch for en masse control of MIG expression and the biological processes they regulate. Specifically, this mode of regulation could control cell proliferation and thus body size, an idea supported by domestication syndrome, wherein MIGs are enriched in common candidate animal domestication genes.

Alternative splicing of P2RX7 pre-messenger RNA in health and diseases: Myth or reality?

Alternative splicing (AS) tremendously increases the use of genetic information by generating protein isoforms that differ in protein-protein interactions, catalytic activity and/or subcellular localization. This review is not dedicated to AS in general, but rather we focus our attention on AS of P2RX7 pre-mRNA. Whereas P2RX7 mRNA is expressed by virtually all eukaryotic mammalian cells, the expression of this channel receptor is restrained to certain cells. When expressed at the cell membrane, P2RX7 controls downstream events including release of inflammatory molecules, phagocytosis, cell proliferation and death and metabolic events. Therefore, P2RX7 is an important actor of health and diseases. In this review, we summarize the general mechanisms leading to AS. Further, we recapitulate our current knowledge concerning the functional regions in P2RX7, identified at the genetic or exonic levels, and how AS may affect the expression of these regions. Finally, the potential of P2RX7 splice variants to control the fate of cancer cells is discussed.

The main WAP isoform usually found in camel milk arises from the usage of an improbable intron cryptic splice site in the precursor to mRNA in which a GC-AG intron occurs

Background

Whey acidic protein (WAP) is a major protein identified in the milk of several mammalian species with cysteine-rich domains known as four-disulfide cores (4-DSC). The organization of the eutherian WAP genes is highly conserved through evolution. It has been proposed that WAP could play an important role in regulating the proliferation of mammary epithelial cells. A bacteriostatic activity was also reported. Conversely to the other mammalian species expressing WAP in their milk, camel WAP contains 4 additional amino acid residues at the beginning of the second 4-DSC domain, introducing a phosphorylation site. The aim of this study was to elucidate the origin of this specificity, which possibly impacts its physiological functions.

Results

Using LC-ESI-MS, we identified in Camelus bactrianus from Kazakhstan a phosphorylated whey protein, exhibiting a molecular mass (12,596 Da), 32 Da higher than the original WAP (12,564 Da) and co-eluting with WAP. cDNA sequencing revealed a transition G/A, which modifies an amino acid residue of the mature protein (V12 M), accounting for the mass difference observed between WAP genetic variants. We also report the existence of two splicing variants of camel WAP precursors to mRNA, arising from an alternative usage of the canonical splice site recognized as such in the other mammalian species. However, the major camel WAP isoform results from the usage of an unlikely intron cryptic splice site, extending camel exon 3 upstream by 12-nucleotides encoding the 4 additional amino acid residues (VSSP) in which a potentially phosphorylable Serine residue occurs. Combining protein and cDNA sequences with genome data available (NCBI database), we report another feature of the camel WAP gene which displays a very rare GC-AG type intron. This result was confirmed by sequencing a genomic DNA fragment encompassing exon 3 to exon 4, suggesting for the GC donor site a compensatory effect in terms of consensus at the acceptor exon position.

Conclusions

Combining proteomic and molecular biology approaches we report: the characterization of a new genetic variant of camel WAP, the usage of an unlikely intron cryptic splice site, and the occurrence of an extremely rare GC-AG type of intron.

Late diagnosis of Alstrom syndrome in a Yemenite-Jewish child

BACKGROUND

We describe the ophthalmologic, clinical, and genetic findings in a patient of Yemenite-Jewish origin diagnosed with Alstrom syndrome due to a novel splice-site mutation 10 years after a clinical misdiagnosis of Leber congenital amaurosis.

METHODS

Ophthalmological evaluations included visual acuity, cycloplegic refraction, slit-lamp, and optical coherent tomography. Genetic analyses included whole exome sequencing followed by bioinformatics analysis and segregation analysis. An in vitro splicing assay was used to evaluate the effect of the identified mutation on splicing. Taqman assay was used to determine the need for population screening for the identified mutation.

RESULTS

Ophthalmologic findings at age 6 were impaired vision, nystagmus, and hyperopia. At age 16 years, the patient presented with obesity, hypothyroidism, and elevated transaminase levels in addition to reduced vision, wandering nystagmus, disc pallor, and degenerative retinal changes. Targeted genetic analysis of ALMS1 revealed a homozygous transversion, c.11544 + 3A>T, suggesting a novel splicing mutation, with elimination of the donor splice site and insertion of 73 nucleotides at the end of exon 16. These changes were validated by Sanger sequencing and co-segregation on family members.

CONCLUSIONS

Ophthalmologists should be alert to the differential diagnosis of inherited retinal degeneration in young patients who present with impaired vision, especially if systemic symptoms are mild and there is no known family history. In the present case, targeted genetic analysis of a child with a syndromic cone-rod dystrophy yielded a novel splicing mutation in ALMS1 causing Alstrom syndrome.

Detection of Viral RNA Splicing in Diagnostic Virology

This chapter is the first one to introduce the detection of viral RNA splicing as a new tool for clinical diagnosis of virus infections. These include various infections caused by influenza viruses, human immunodeficiency viruses (HIV), human T-cell leukemia viruses (HTLV), Torque teno viruses (TTV), parvoviruses, adenoviruses, hepatitis B virus, polyomaviruses, herpesviruses, and papillomaviruses. Detection of viral RNA splicing for active viral gene expression in a clinical sample is a nucleic acid-based detection. The interpretation of the detected viral RNA splicing results is straightforward without concern for carry-over DNA contamination, because the spliced RNA is smaller than its corresponding DNA template.

Although many methods can be used, a simple method to detect viral RNA splicing is reverse transcription-polymerase chain reaction (RT-PCR). In principle, the detection of spliced RNA transcripts by RT-PCR depends on amplicon selection and primer design. The most common approach is the amplification over the intron regions by a set of primers in flanking exons. A larger product than the predicted size of smaller, spliced RNA is in general an unspliced RNA or contaminating viral genomic DNA. A spliced mRNA always gives a smaller RT-PCR product than its unspliced RNA due to removal of intron sequences by RNA splicing. The contaminating viral DNA can be determined by a minus RT amplification (PCR). Alternatively, specific amplification of a spliced RNA can be obtained by using an exon-exon junction primer because the sequence at exon-exon junction is not present in the unspliced RNA nor in viral genomic DNA.

Short and narrow flag leaf1, a GATA zinc finger domain-containing protein, regulates flag leaf size in rice (Oryza sativa)

BACKGROUND

The flag leaf of rice (Oryza sativa L.) is an important determinant of plant type characteristics and grain yield. Identification of flag leaf mutants of rice is crucial to elucidate the molecular mechanism of flag-leaf development, and for exploitation of rice germplasm resources.

RESULTS

In this study, we describe a mutant designated short and narrow flag leaf 1 (snfl1). Histological analysis showed that the length of epidermal cells and number of longitudinal veins were decreased in the flag leaf of the snfl1 mutant. Map-based cloning indicated that a member of the GATA family of transcription factors is a candidate gene for SNFL1. A single-nucleotide transition at the last base in the single intron of snfl1 led to variation in alternative splicing and early termination of translation. Complemented transgenic plants harbouring the candidate SNFL1 gene rescued the snfl1 mutant. Analysis of RT-PCR and the SNFL1 promoter by means of a GUS fusion expression assay showed that abundance of SNFL1 transcripts was higher in the culm, leaf sheath, and root. Expression of the SNFL1-GFP fusion protein in rice protoplasts showed that SNFL1 was localized in nucleus.

CONCLUSIONS

We conclude that SNFL1 is an important regulator of leaf development, the identification of which might have important implications for future research on GATA transcription factors.

Modeling the structural implications of an alternatively spliced Exoc3l2, a paralog of the tunneling nanotube-forming M-Sec

The exocyst is a molecular tether that retains secretory vesicles at the plasma membrane prior to SNARE-mediated docking and fusion. However, individual exocyst complex components (EXOCs) may also function independently of exocyst assembly. Alternative splice variants of EXOC mRNA and paralogs of EXOC genes have been described and several have been attributed functions that may be independent of the exocyst complex. Here we describe a novel splice variant of murine Exoc3l2, which we term Exoc3l2a. We discuss possible functional implications of the resulting domain excision from this isoform of EXOC3L2 based on structural similarities with its paralog M-Sec (EXOC3L3), which is implicated in tunneling nanotube formation. The identification of this Exoc3l2 splice variant expands the potential for subunit diversity within the exocyst and for alternative functionality of this component independently of the exocyst.

The matrices and constraints of GT/AG splice sites of more than 1000 species/lineages

To provide a resource for the splice sites (SS) of different species, we calculated the matrices of nucleotide compositions of about 38 million splice sites from >1000 species/lineages. The matrices are enriched of aGGTAAGT (5'SS) or (Y)₆N(C/t)AG(g/a)t (3'SS) overall; however, they are quite diverse among hundreds of species. The diverse matrices remain prominent even under sequence selection pressures, suggesting the existence of diverse constraints as well as U snRNAs and other spliceosomal factors and/or their interactions with the splice sites. Using an algorithm to measure and compare the splice site constraints across all species, we demonstrate their distinct differences quantitatively. As an example of the resource's application to answering specific questions, we confirm that high constraints of particular positions are significantly associated with transcriptome-wide, increased occurrences of alternative splicing when uncommon nucleotides are present. More interestingly, the abundance of alternative splicing in 16 species correlates with the average constraint index of splice sites in a bell curve. This resource will allow users to assess specific sequences/splice sites against the consensus of every Ensembl-annotated species, and to explore the evolutionary changes or relationship to alternative splicing and transcriptome diversity. Web-search or update features are also included.

The emerging role of minor intron splicing in neurological disorders

Pre-mRNA splicing is an essential step in eukaryotic gene expression. Mutations in cis-acting sequence elements within pre-mRNA molecules or trans-acting factors involved in pre-mRNA processing have both been linked to splicing dysfunction that give rise to a large number of human diseases. These mutations typically affect the major splicing pathway, which excises more than 99% of all introns in humans. However, approximately 700-800 human introns feature divergent intron consensus sequences at their 5' and 3' ends and are recognized by a separate pre-mRNA processing machinery denoted as the minor spliceosome. This spliceosome has been studied less than its major counterpart, but has received increasing attention during the last few years as a novel pathomechanistic player on the stage in neurodevelopmental and neurodegenerative diseases. Here, we review the current knowledge on minor spliceosome function and discuss its potential pathomechanistic role and impact in neurodegeneration.

Minor spliceosome and disease

The U12-dependent (minor) spliceosome excises a rare group of introns that are characterized by a highly conserved 5' splice site and branch point sequence. Several new congenital or somatic diseases have recently been associated with mutations in components of the minor spliceosome. A common theme in these diseases is the detection of elevated levels of transcripts containing U12-type introns, of which a subset is associated with other splicing defects. Here we review the present understanding of minor spliceosome diseases, particularly those associated with the specific components of the minor spliceosome. We also present a model for interpreting the molecular-level consequences of the different diseases.

Consideration of non-canonical splice sites improves gene prediction on the Arabidopsis thaliana Niederzenz-1 genome sequence

Objective

The Arabidopsis thaliana Niederzenz-1 genome sequence was recently published with an ab initio gene prediction. In depth analysis of the predicted gene set revealed some errors involving genes with non-canonical splice sites in their introns. Since non-canonical splice sites are difficult to predict ab initio, we checked for options to improve the annotation by transferring annotation information from the recently released Columbia-0 reference genome sequence annotation Araport11.

Results

Incorporation of hints generated from Araport11 enabled the precise prediction of non-canonical splice sites. Manual inspection of RNA-Seq read mapping and RT-PCR were applied to validate the structural annotations of non-canonical splice sites. Predictions of untranslated regions were also updated by harnessing the potential of Araport11’s information, which was generated by using high coverage RNA-Seq data. The improved gene set of the Nd-1 genome assembly (GeneSet_Nd-1_v1.1) was evaluated via comparison to the initial gene prediction (GeneSet_Nd-1_v1.0) as well as against Araport11 for the Col-0 reference genome sequence. GeneSet_Nd-1_v1.1 contains previously missed non-canonical splice sites in 1256 genes. Reciprocal best hits for 24,527 (89.4%) of all nuclear Col-0 genes against the GeneSet_Nd-1_v1.1 indicate a high gene prediction quality.

Electronic supplementary material

The online version of this article (10.1186/s13104-017-2985-y) contains supplementary material, which is available to authorized users.

RNA-binding proteins in neurodegeneration: mechanisms in aggregate

Neurodegeneration is a leading cause of death in the developed world and a natural, albeit unfortunate, consequence of longer-lived populations. Despite great demand for therapeutic intervention, it is often the case that these diseases are insufficiently understood at the basic molecular level. What little is known has prompted much hopeful speculation about a generalized mechanistic thread that ties these disparate conditions together at the subcellular level and can be exploited for broad curative benefit. In this review, we discuss a prominent theory supported by genetic and pathological changes in an array of neurodegenerative diseases: that neurons are particularly vulnerable to disruption of RNA-binding protein dosage and dynamics. Here we synthesize the progress made at the clinical, genetic, and biophysical levels and conclude that this perspective offers the most parsimonious explanation for these mysterious diseases. Where appropriate, we highlight the reciprocal benefits of cross-disciplinary collaboration between disease specialists and RNA biologists as we envision a future in which neurodegeneration declines and our understanding of the broad importance of RNA processing deepens.

Splicing features in the expression of the complementary-sense genes of Beet curly top Iran virus

Beet curly top Iran virus (BCTIV) is a distinct geminivirus which has been reported from sugar-beet-growing farms in Iran. In this study, the role of the splicing in expression of complementary-sense genes of BCTIV was studied. Total RNA was extracted from BCTIV-infected tissue, and the predicted intron position of complementary-sense mRNA transcripts was amplified by RT-PCR followed by cloning of the amplicons. Sequence confirmed that both spliced and unspliced mRNAs are synthesized by the same transcription unit. Sequence comparison showed that a 155-nt segment (intron) corresponding to nucleotides 1890-2044 of the viral genome has been removed from the latter transcript and therefore fusion of the C1:C2 genes resulted creation of a continuous reading frame for potential production of intact replication initiator protein (Rep). BCTIV intron comprises of most consensus splicing signals required for splicing in eukaryotes and several plant viruses including mastre- and capulaviruses.

Sequence variation between 462 human individuals fine-tunes functional sites of RNA processing

Recent advances in the cost-efficiency of sequencing technologies enabled the combined DNA- and RNA-sequencing of human individuals at the population-scale, making genome-wide investigations of the inter-individual genetic impact on gene expression viable. Employing mRNA-sequencing data from the Geuvadis Project and genome sequencing data from the 1000 Genomes Project we show that the computational analysis of DNA sequences around splice sites and poly-A signals is able to explain several observations in the phenotype data. In contrast to widespread assessments of statistically significant associations between DNA polymorphisms and quantitative traits, we developed a computational tool to pinpoint the molecular mechanisms by which genetic markers drive variation in RNA-processing, cataloguing and classifying alleles that change the affinity of core RNA elements to their recognizing factors. The in silico models we employ further suggest RNA editing can moonlight as a splicing-modulator, albeit less frequently than genomic sequence diversity. Beyond existing annotations, we demonstrate that the ultra-high resolution of RNA-Seq combined from 462 individuals also provides evidence for thousands of bona fide novel elements of RNA processing-alternative splice sites, introns, and cleavage sites-which are often rare and lowly expressed but in other characteristics similar to their annotated counterparts.

Physiological and Pathophysiological Insights of Nav1.4 and Nav1.5 Comparison

Mutations in Nav1.4 and Nav1.5 α-subunits have been associated with muscular and cardiac channelopathies, respectively. Despite intense research on the structure and function of these channels, a lot of information is still missing to delineate the various physiological and pathophysiological processes underlying their activity at the molecular level. Nav1.4 and Nav1.5 sequences are similar, suggesting structural and functional homologies between the two orthologous channels. This also suggests that any characteristics described for one channel subunit may shed light on the properties of the counterpart channel subunit. In this review article, after a brief clinical description of the muscular and cardiac channelopathies related to Nav1.4 and Nav1.5 mutations, respectively, we compare the knowledge accumulated in different aspects of the expression and function of Nav1.4 and Nav1.5 α-subunits: the regulation of the two encoding genes (SCN4A and SCN5A), the associated/regulatory proteins and at last, the functional effect of the same missense mutations detected in Nav1.4 and Nav1.5. First, it appears that more is known on Nav1.5 expression and accessory proteins. Because of the high homologies of Nav1.5 binding sites and equivalent Nav1.4 sites, Nav1.5-related results may guide future investigations on Nav1.4. Second, the analysis of the same missense mutations in Nav1.4 and Nav1.5 revealed intriguing similarities regarding their effects on membrane excitability and alteration in channel biophysics. We believe that such comparison may bring new cues to the physiopathology of cardiac and muscular diseases.

Genetic characterization of interleukins (IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18) with relevant biological roles in lagomorphs

ILs, as essential innate immune modulators, are involved in an array of biological processes. In the European rabbit (Oryctolagus cuniculus) IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18 have been implicated in inflammatory processes and in the immune response against rabbit hemorrhagic disease virus and myxoma virus infections. In this study we characterized these ILs in six Lagomorpha species (European rabbit, pygmy rabbit, two cottontail rabbit species, European brown hare and American pika). Overall, these ILs are conserved between lagomorphs, including in their exon/intron structure. Most differences were observed between leporids and American pika. Indeed, when comparing both, some relevant differences were observed in American pika, such as the location of the stop codon in IL-1α and IL-2, the existence of a different transcript in IL8 and the number of cysteine residues in IL-1β. Changes at N-glycosylation motifs were also detected in IL-1, IL-10, IL-12B and IL-15. IL-1α is the protein that presents the highest evolutionary distances, which is in contrast to IL-12A where the distances between lagomorphs are the lowest. For all these ILs, sequences of human and European rabbit are more closely related than between human and mouse or European rabbit and mouse.

Dinoflagellate Gene Structure and Intron Splice Sites in a Genomic Tandem Array

Dinoflagellates are one of the last major lineages of eukaryotes for which little is known about genome structure and organization. We report here the sequence and gene structure of a clone isolated from a cosmid library which, to our knowledge, represents the largest contiguously sequenced, dinoflagellate genomic, tandem gene array. These data, combined with information from a large transcriptomic library, allowed a high level of confidence of every base pair call. This degree of confidence is not possible with PCR‐based contigs. The sequence contains an intron‐rich set of five highly expressed gene repeats arranged in tandem. One of the tandem repeat gene members contains an intron 26,372 bp long. This study characterizes a splice site consensus sequence for dinoflagellate introns. Two to nine base pairs around the 3′ splice site are repeated by an identical two to nine base pairs around the 5′ splice site. The 5′ and 3′ splice sites are in the same locations within each repeat so that the repeat is found only once in the mature mRNA. This identically repeated intron boundary sequence might be useful in gene modeling and annotation of genomes.

Assays for the identification and prioritization of drug candidates for spinal muscular atrophy

Spinal muscular atrophy (SMA) is an autosomal recessive genetic disorder resulting in degeneration of α-motor neurons of the anterior horn and proximal muscle weakness. It is the leading cause of genetic mortality in children younger than 2 years. It affects ∼1 in 11,000 live births. In 95% of cases, SMA is caused by homozygous deletion of the SMN1 gene. In addition, all patients possess at least one copy of an almost identical gene called SMN2. A single point mutation in exon 7 of the SMN2 gene results in the production of low levels of full-length survival of motor neuron (SMN) protein at amounts insufficient to compensate for the loss of the SMN1 gene. Although no drug treatments are available for SMA, a number of drug discovery and development programs are ongoing, with several currently in clinical trials. This review describes the assays used to identify candidate drugs for SMA that modulate SMN2 gene expression by various means. Specifically, it discusses the use of high-throughput screening to identify candidate molecules from primary screens, as well as the technical aspects of a number of widely used secondary assays to assess SMN messenger ribonucleic acid (mRNA) and protein expression, localization, and function. Finally, it describes the process of iterative drug optimization utilized during preclinical SMA drug development to identify clinical candidates for testing in human clinical trials.

Pectinases from Sphenophorus levis Vaurie, 1978 (Coleoptera: Curculionidae): putative accessory digestive enzymes

The cell wall in plants offers protection against invading organisms and is mainly composed of the polysaccharides pectin, cellulose, and hemicellulose, which can be degraded by plant cell wall degrading enzymes (PCWDEs). Such enzymes are often synthesized by free living microorganisms or endosymbionts that live in the gut of some animals, including certain phytophagous insects. Thus, the ability of an insect to degrade the cell wall was once thought to be related to endosymbiont enzyme activity. However, recent studies have revealed that some phytophagous insects are able to synthesize their own PCWDEs by endogenous genes, although questions regarding the origin of these genes remain unclear. This study describes two pectinases from the sugarcane weevil, Sphenophorus levis Vaurie, 1978 (Sl-pectinases), which is considered one of the most serious agricultural pests in Brazil. Two cDNA sequences identified in a cDNA library of the insect larvae coding for a pectin methylesterase (PME) and an endo-polygalacturonase (endo-PG)-denominated Sl-PME and Sl-endoPG, respectively-were isolated and characterized. The quantitative real-time reverse transcriptase polymerase chain reaction expression profile for both Sl-pectinases showed mRNA production mainly in the insect feeding stages and exclusively in midgut tissue of the larvae. This analysis, together Western blotting data, suggests that Sl-pectinases have a digestive role. Phylogenetic analyses indicate that Sl-PME and Sl-endoPG sequences are closely related to bacteria and fungi, respectively. Moreover, the partial genomic sequences of the pectinases were amplified from insect fat body DNA, which was certified to be free of endosymbiotic DNA. The analysis of genomic sequences revealed the existence of two small introns with 53 and 166 bp in Sl-endoPG, which is similar to the common pattern in fungal introns. In contrast, no intron was identified in the Sl-PME genomic sequence, as generally observed in bacteria. These data support the theory of horizontal gene transfer proposed for the origin of insect pectinases, reinforcing the acquisition of PME genes from bacteria and endo-PG genes from fungi.

VEGF Splicing and the Role of VEGF Splice Variants: From Physiological-Pathological Conditions to Specific Pre-mRNA Splicing

During this past decade, the vascular endothelial growth factor (VEGF) pathway has been extensively studied. VEGF is a paradigm of molecular regulation since its expression is controlled at all possible steps including transcription, mRNA stability, translation, and pre-mRNA splicing. The latter form of molecular regulation is probably the least studied. This field has been neglected; yet different forms of VEGF with different sizes and different physiological properties issued from alternative splicing have been described a long time ago. Recently a new level of complexity was added to the field of splicing of VEGF pre-mRNA. Whereas thousands of publications have described VEGF as a pro-angiogenic factor, an alternative splicing event generates specific anti-angiogenic forms of VEGF that only differ from the others by a modification in the last six amino acids of the protein. According to the scientists who discovered these isoforms, which are indistinguishable from the pro-angiogenic ones with pan VEGF antibodies, some of the literature on VEGF is at least inexact if not completely false. Moreover, the presence of anti-angiogenic forms of VEGF may explain the disappointing efficacy of anti-VEGF therapies on the overall survival of patients with different forms of cancers and with wet age-related macular degeneration. This review focuses on the existence of the different alternative splice variants of VEGF and the molecular mechanisms associated with their expression and function.

Therapeutic strategies based on modified U1 snRNAs and chaperones for Sanfilippo C splicing mutations

BACKGROUND

Mutations affecting RNA splicing represent more than 20% of the mutant alleles in Sanfilippo syndrome type C, a rare lysosomal storage disorder that causes severe neurodegeneration. Many of these mutations are localized in the conserved donor or acceptor splice sites, while few are found in the nearby nucleotides.

METHODS

In this study we tested several therapeutic approaches specifically designed for different splicing mutations depending on how the mutations affect mRNA processing. For three mutations that affect the donor site (c.234 + 1G > A, c.633 + 1G > A and c.1542 + 4dupA), different modified U1 snRNAs recognizing the mutated donor sites, have been developed in an attempt to rescue the normal splicing process. For another mutation that affects an acceptor splice site (c.372-2A > G) and gives rise to a protein lacking four amino acids, a competitive inhibitor of the HGSNAT protein, glucosamine, was tested as a pharmacological chaperone to correct the aberrant folding and to restore the normal trafficking of the protein to the lysosome.

RESULTS

Partial correction of c.234 + 1G > A mutation was achieved with a modified U1 snRNA that completely matches the splice donor site suggesting that these molecules may have a therapeutic potential for some splicing mutations. Furthermore, the importance of the splice site sequence context is highlighted as a key factor in the success of this type of therapy. Additionally, glucosamine treatment resulted in an increase in the enzymatic activity, indicating a partial recovery of the correct folding.

CONCLUSIONS

We have assayed two therapeutic strategies for different splicing mutations with promising results for the future applications.

Lack of variation in voltage-gated sodium channels of common bottlenose dolphins (Tursiops truncatus) exposed to neurotoxic algal blooms

In coastal marine ecosystems, neurotoxins produced by harmful algal blooms (HABs) often result in large-scale mortality events of many marine species. Historical and frequent exposure to HABs therefore may provide a strong selective pressure for adaptations that result in toxin resistance. Neurotoxin resistance has independently evolved in a variety of terrestrial and marine species via mutations in genes encoding the toxin binding sites within the voltage-gated sodium channel gene complex. Accordingly, we tested the hypothesis that genetic variation in the putative binding site of brevetoxins in common bottlenose dolphins (Tursiops truncatus) explains differences among individuals or populations in resistance to harmful Karenia brevis blooms in the Gulf of Mexico. We found very little variation in the sodium channel exons encoding the putative brevetoxin binding site among bottlenose dolphins from central-west Florida and the Florida Panhandle. Our study included samples from several bottlenose dolphin mortality events associated with HABs, but we found no association between genetic variation and survival. We observed a significant effect of geographic region on genetic variation for some sodium channel isoforms, but this can be primarily explained by rare private alleles and is more likely a reflection of regional genetic differentiation than the cause of different levels of HAB resistance between regions. In contrast to many other previously studied neurotoxin-resistant species, we conclude that bottlenose dolphins have not evolved resistance to HABs via mutations in genes encoding the brevetoxin binding site on the voltage-gated sodium channels.

The Arabidopsis U11/U12-65K is an indispensible component of minor spliceosome and plays a crucial role in U12 intron splicing and plant development

The U12-dependent introns have been identified in a wide range of eukaryotes and are removed from precursor-mRNAs by U12 intron-specific minor spliceosome. Although several proteins unique to minor spliceosome have been identified, the nature of their effect on U12 intron splicing as well as plant growth and development remain largely unknown. Here, we characterized the functional role of an U12-type spliceosomal protein, U11/U12-65K in Arabidopsis thaliana. The transgenic knockdown plants generated by artificial miRNA-mediated silencing strategy exhibited severe defect in growth and development, such as severely arrested primary inflorescence stems, serrated leaves, and the formation of many rosette leaves after bolting. RNA sequencing and reverse transcription polymerase chain reaction (RT-PCR) analyses revealed that splicing of 198 out of the 234 previously predicted U12 intron-containing genes and 32 previously unidentified U12 introns was impaired in u11/u12-65k mutant. Moreover, the U11/U12-65K mutation affected alternative splicing, as well as U12 intron splicing, of many introns. Microarray analysis revealed that the genes involved in cell wall biogenesis and function, plant development, and metabolic processes are differentially expressed in the mutant plants. U11/U12-65K protein bound specifically to U12 small nuclear RNA (snRNA), which is necessary for branch-point site recognition. Taken together, these results provide clear evidence that U11/U12-65K is an indispensible component of minor spliceosome and involved in U12 intron splicing and alternative splicing of many introns, which is crucial for plant development.

Deducing intracellular distributions of metabolic pathways from genomic data

In the recent years, a large number of genomes from a variety of different organisms have been sequenced. Most of the sequence data has been publicly released and can be assessed by interested users. However, this wealth of information is currently underexploited by scientists not directly involved in genome annotation. This is partially because sequencing, assembly, and automated annotation can be done much faster than the identification, classification, and prediction of the intracellular localization of the gene products. This part of the annotation process still largely relies on manual curation and addition of contextual information. Users of genome databases who are unfamiliar with the types of data available from (whole) genomes might therefore find themselves either overwhelmed by the vast amount and multiple layers of data or dissatisfied with less-than-meaningful analyses of the data.In this chapter we present procedures and approaches to identify and characterize gene models of enzymes involved in metabolic pathways based on their similarity to known sequences. Furthermore we describe how to predict the subcellular location of the proteins using publicly available prediction servers and how to interpret the obtained results. The strategies we describe are generally applicable to organisms with primary plastids such as land plants or green algae. Additionally, we describe strategies suitable for those groups of algae with secondary plastids (for instance diatoms), which are characterized by a different cellular topology and a larger number of intracellular compartments compared to plants.

A second rhodopsin-like protein in Cyanophora paradoxa: gene sequence and protein expression in a cell-free system

Here we report the identification and expression of a second rhodopsin-like protein in the alga Cyanophora paradoxa (Glaucophyta), named Cyanophopsin_2. This new protein was identified due to a serendipity event, since the RACE reaction performed to complete the sequence of Cyanophopsin_1, (the first rhodopsin-like protein of C. paradoxa identified in 2009 by our group), amplified a 619 bp sequence corresponding to a portion of a new gene of the same protein family. The full sequence consists of 1175 bp consisting of 849 bp coding DNA sequence and 4 introns of 326 bp. The protein is characterized by an N-terminal region of 47 amino acids, followed by a region with 7 α-helices of 213 amino acids and a C-terminal region of 22 amino acids. This protein showed high identity with Cyanophopsin_1 and other rhodopsin-like proteins of Archea, Bacteria, Fungi and Algae. Cyanophosin_2 (CpR2) was expressed in a cell-free expression system, and characterized by means of absorption spectroscopy.

Aberrant splicing in neurological diseases

Splicing of precursor messenger RNA (pre-mRNA) removes the intervening sequences (introns) and joins the expressed regions (exons) in the nucleus, before an intron-containing eukaryotic mRNA transcript can be exported and translated into proteins in the cytoplasm. While some sequences are always included or excluded (constitutive splicing), others can be selectively used (alternative splicing) in this process. Particularly by alternative splicing, up to tens of thousands of variant transcripts can be produced from a single gene, which contributes greatly to the proteomic diversity for such complex cellular functions as 'wiring' neurons in the nervous system. Disruption of this process leads to aberrant splicing, which accounts for the defects of up to 50% of mutations that cause certain human genetic diseases. In this review, we describe the different mechanisms of aberrant splicing that cause or have been associated with neurological diseases.

Splicing therapy for neuromuscular disease

Duchenne muscular dystrophy (DMD) and spinal muscular atrophy (SMA) are two of the most common inherited neuromuscular diseases in humans. Both conditions are fatal and no clinically available treatments are able to significantly alter disease course in either case. However, by manipulation of pre-mRNA splicing using antisense oligonucleotides, defective transcripts from the DMD gene and from the SMN2 gene in SMA can be modified to once again produce protein and restore function. A large number of in vitro and in vivo studies have validated the applicability of this approach and an increasing number of preliminary clinical trials have either been completed or are under way. Several different oligonucleotide chemistries can be used for this purpose and various strategies are being developed to facilitate increased delivery efficiency and prolonged therapeutic effect. As these novel therapeutic compounds start to enter the clinical arena, attention must also be drawn to the question of how best to facilitate the clinical development of such personalised genetic therapies and how best to implement their provision.

The significant other: splicing by the minor spliceosome

The removal of non-coding sequences, introns, from the mRNA precursors is an essential step in eukaryotic gene expression. U12-type introns are a minor subgroup of introns, distinct from the major or U2-type introns. U12-type introns are present in most eukaryotes but only account for less than 0.5% of all introns in any given genome. They are processed by a specific U12-dependent spliceosome, which is similar to, but distinct from, the major spliceosome. U12-type introns are spliced somewhat less efficiently than the major introns, and it is believed that this limits the expression of the genes containing such introns. Recent findings on the role of U12-dependent splicing in development and human disease have shown that it can also affect multiple cellular processes not directly related to the functions of the host genes of U12-type introns. At the same time, advances in understanding the regulation and phylogenetic distribution of the minor spliceosome are starting to shed light on how the U12-type introns and the minor spliceosome may have evolved. © 2012 John Wiley & Sons, Ltd.

SMN-inducing compounds for the treatment of spinal muscular atrophy

Spinal muscular atrophy (SMA) is a leading genetic cause of infant mortality. A neurodegenerative disease, it is caused by loss of SMN1, although low, but essential, levels of SMN protein are produced by the nearly identical gene SMN2. While no effective treatment or therapy currently exists, a new wave of therapeutics has rapidly progressed from cell-based and preclinical animal models to the point where clinical trials have initiated for SMA-specific compounds. There are several reasons why SMA has moved relatively rapidly towards novel therapeutics, including: SMA is monogenic; the molecular understanding of SMN gene regulation has been building for nearly 20 years; and all SMA patients retain one or more copies of SMN2 that produces low levels of full-length, fully functional SMN protein. This review primarily focuses upon the biology behind the disease and examines SMN1- and SMN2-targeted therapeutics.

Spliceosomal small nuclear ribonucleoprotein biogenesis defects and motor neuron selectivity in spinal muscular atrophy

The SMN protein is essential and participates in the assembly of macromolecular complexes of RNA and protein in all cells. The best-characterized function of SMN is as an assembler of spliceosomal small nuclear ribonucleoproteins (snRNPs). SMN performs this function as part of a complex with several other proteins called Gemins. snRNPs are assembled in the cytoplasm in a stepwise manner and then are imported to the nucleus where they participate globally in the splicing of pre-mRNA. Mutations in the SMN1 gene result in the motor neuron disease, spinal muscular atrophy (SMA). Most of these mutations result in a reduction in the expression levels of the SMN protein, which, in turn, results in a reduction in snRNP assembly capacity. This review highlights current studies that have investigated the mechanism of SMN-dependent snRNP assembly, as well as the downstream effects on pre-mRNA splicing that result from a decrease in SMN. This article is part of a Special Issue entitled "RNA-Binding Proteins".

U12-type spliceosomal introns of Insecta

Most of eukaryotic genes are interrupted by introns that need to be removed from pre-mRNAs before they can perform their function. This is done by complex machinery called spliceosome. Many eukaryotes possess two separate spliceosomal systems that process separate sets of introns. The major (U2) spliceosome removes majority of introns, while minute fraction of intron repertoire is processed by the minor (U12) spliceosome. These two populations of introns are called U2-type and U12-type, respectively. The latter fall into two subtypes based on the terminal dinucleotides. The minor spliceosomal system has been lost independently in some lineages, while in some others few U12-type introns persist. We investigated twenty insect genomes in order to better understand the evolutionary dynamics of U12-type introns. Our work confirms dramatic drop of U12-type introns in Diptera, leaving these genomes just with a handful cases. This is mostly the result of intron deletion, but in a number of dipteral cases, minor type introns were switched to a major type, as well. Insect genes that harbor U12-type introns belong to several functional categories among which proteins binding ions and nucleic acids are enriched and these few categories are also overrepresented among these genes that preserved minor type introns in Diptera.

Selection of reference genes for qPCR in hairy root cultures of peanut

Background

Hairy root cultures produced via Agrobacterium rhizogenes-mediated transformation have emerged as practical biological models to elucidate the biosynthesis of specialized metabolites. To effectively understand the expression patterns of the genes involved in the metabolic pathways of these compounds, reference genes need to be systematically validated under specific experimental conditions as established by the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) guidelines. In the present report we describe the first validation of reference genes for RT-qPCR in hairy root cultures of peanut which produce stilbenoids upon elicitor treatments.

Results

A total of 21 candidate reference genes were evaluated. Nineteen genes were selected based on previous qPCR studies in plants and two were from the T-DNAs transferred from A. rhizogenes. Nucleotide sequences of peanut candidate genes were obtained using their homologous sequences in Arabidopsis. To identify the suitable primers, calibration curves were obtained for each candidate reference gene. After data analysis, 12 candidate genes meeting standard efficiency criteria were selected. The expression stability of these genes was analyzed using geNorm and NormFinder algorithms and a ranking was established based on expression stability of the genes. Candidate reference gene expression was shown to have less variation in methyl jasmonate (MeJA) treated root cultures than those treated with sodium acetate (NaOAc).

Conclusions

This work constitutes the first effort to validate reference genes for RT-qPCR in hairy roots. While these genes were selected under conditions of NaOAc and MeJA treatment, we anticipate these genes to provide good targets for reference genes for hairy roots under a variety of stress conditions. The lead reference genes were a gene encoding for a TATA box binding protein (TBP2) and a gene encoding a ribosomal protein (RPL8C). A commonly used reference gene GAPDH showed low stability of expression suggesting that its use may lead to inaccurate gene expression profiles when used for data normalization in stress-stimulated hairy roots. Likewise the A. rhizogenes transgene rolC showed less expression stability than GAPDH. This study proposes that a minimum of two reference genes should be used for a normalization procedure in gene expression profiling using elicited hairy roots.

A mutation in a rare type of intron in a sodium-channel gene results in aberrant splicing and causes myotonia

Many mutations in the skeletal-muscle sodium-channel gene SCN4A have been associated with myotonia and/or periodic paralysis, but so far all of these mutations are located in exons. We found a patient with myotonia caused by a deletion/insertion located in intron 21 of SCN4A, which is an AT-AC type II intron. This is a rare class of introns that, despite having AT-AC boundaries, are spliced by the major or U2-type spliceosome. The patient's skeletal muscle expressed aberrantly spliced SCN4A mRNA isoforms generated by activation of cryptic splice sites. In addition, genetic suppression experiments using an SCN4A minigene showed that the mutant 5' splice site has impaired binding to the U1 and U6 snRNPs, which are the cognate factors for recognition of U2-type 5' splice sites. One of the aberrantly spliced isoforms encodes a channel with a 35-amino acid insertion in the cytoplasmic loop between domains III and IV of Nav1.4. The mutant channel exhibited a marked disruption of fast inactivation, and a simulation in silico showed that the channel defect is consistent with the patient's myotonic symptoms. This is the first report of a disease-associated mutation in an AT-AC type II intron, and also the first intronic mutation in a voltage-gated ion channel gene showing a gain-of-function defect.

The proteogenomic mapping tool

Background

High-throughput mass spectrometry (MS) proteomics data is increasingly being used to complement traditional structural genome annotation methods. To keep pace with the high speed of experimental data generation and to aid in structural genome annotation, experimentally observed peptides need to be mapped back to their source genome location quickly and exactly. Previously, the tools to do this have been limited to custom scripts designed by individual research groups to analyze their own data, are generally not widely available, and do not scale well with large eukaryotic genomes.

Results

The Proteogenomic Mapping Tool includes a Java implementation of the Aho-Corasick string searching algorithm which takes as input standardized file types and rapidly searches experimentally observed peptides against a given genome translated in all 6 reading frames for exact matches. The Java implementation allows the application to scale well with larger eukaryotic genomes while providing cross-platform functionality.

Conclusions

The Proteogenomic Mapping Tool provides a standalone application for mapping peptides back to their source genome on a number of operating system platforms with standard desktop computer hardware and executes very rapidly for a variety of datasets. Allowing the selection of different genetic codes for different organisms allows researchers to easily customize the tool to their own research interests and is recommended for anyone working to structurally annotate genomes using MS derived proteomics data.

SMN in spinal muscular atrophy and snRNP biogenesis

Ribonucleoprotein (RNP) complexes function in nearly every facet of cellular activity. The spliceosome is an essential RNP that accurately identifies introns and catalytically removes the intervening sequences, providing exquisite control of spatial, temporal, and developmental gene expressions. U-snRNPs are the building blocks for the spliceosome. A significant amount of insight into the molecular assembly of these essential particles has recently come from a seemingly unexpected area of research: neurodegeneration. Survival motor neuron (SMN) performs an essential role in the maturation of snRNPs, while the homozygous loss of SMN1 results in the development of spinal muscular atrophy (SMA), a devastating neurodegenerative disease. In this review, the function of SMN is examined within the context of snRNP biogenesis and evidence is examined which suggests that the SMN functional defects in snRNP biogenesis may account for the motor neuron pathology observed in SMA.

Multiple neurofilament subunits are present in lamprey CNS

In mammals, there are three neurofilament (NF) subunits (NF-L, NF-M, and NF-H), but it was thought that only a single NF, NF180, exists in lamprey. However, NF180 lacked the ability to self-assemble, suggesting that like mammalian NFs, lamprey NFs are heteropolymers, and that additional NF subunits may exist. The present study provides evidence for the existence of a lamprey NF-L homolog (L-NFL). Genes encoding two new NF-M isoforms (NF132 and NF95) also have been isolated and characterized. With NF180, this makes three NF-M-like isoforms. In situ hybridization showed that all three newly cloned NFs are expressed in spinal cord neurons and in spinal-projecting neurons of the brainstem. Like NF180, there were no KSP multiphosphorylation repeat motifs in the tail regions of NF132 or NF95. NF95 was highly identical to homologous parts of NF180, sharing 2 common pieces of DNA with it. Northern blots suggested that NF95 may be expressed at very low levels in older larvae. The presence of L-NFL in lamprey CNS may support the hypothesis that as in mammals, NFs in lamprey are obligate heteropolymers, in which NF-L is a required subunit.

A comprehensive survey of human polymorphisms at conserved splice dinucleotides and its evolutionary relationship with alternative splicing

BACKGROUND

Alternative splicing (AS) is a key molecular process that endows biological functions with diversity and complexity. Generally, functional redundancy leads to the generation of new functions through relaxation of selective pressure in evolution, as exemplified by duplicated genes. It is also known that alternatively spliced exons (ASEs) are subject to relaxed selective pressure. Within consensus sequences at the splice junctions, the most conserved sites are dinucleotides at both ends of introns (splice dinucleotides). However, a small number of single nucleotide polymorphisms (SNPs) occur at splice dinucleotides. An intriguing question relating to the evolution of AS diversity is whether mutations at splice dinucleotides are maintained as polymorphisms and produce diversity in splice patterns within the human population. We therefore surveyed validated SNPs in the database dbSNP located at splice dinucleotides of all human genes that are defined by the H-Invitational Database.

RESULTS

We found 212 validated SNPs at splice dinucleotides (sdSNPs); these were confirmed to be consistent with the GT-AG rule at either allele. Moreover, 53 of them were observed to neighbor ASEs (AE dinucleotides). No significant differences were observed between sdSNPs at AE dinucleotides and those at constitutive exons (CE dinucleotides) in SNP properties including average heterozygosity, SNP density, ratio of predicted alleles consistent with the GT-AG rule, and scores of splice sites formed with the predicted allele. We also found that the proportion of non-conserved exons was higher for exons with sdSNPs than for other exons.

CONCLUSIONS

sdSNPs are found at CE dinucleotides in addition to those at AE dinucleotides, suggesting two possibilities. First, sdSNPs at CE dinucleotides may be robust against sdSNPs because of unknown mechanisms. Second, similar to sdSNPs at AE dinucleotides, those at CE dinucleotides cause differences in AS patterns because of the arbitrariness in the classification of exons into alternative and constitutive type that varies according to the dataset. Taking into account the absence of differences in sdSNP properties between those at AE and CE dinucleotides, the increased proportion of non-conserved exons found in exons flanked by sdSNPs suggests the hypothesis that sdSNPs are maintained at the splice dinucleotides of newly generated exons at which negative selection pressure is relaxed.

A rhodopsin-like protein in Cyanophora paradoxa: gene sequence and protein immunolocalization

Here, we report the DNA sequence of the rhodopsin gene in the alga Cyanophora paradoxa (Glaucophyta). The primers were designed according to the conserved regions of prokaryotic and eukaryotic rhodopsin-like proteins deposited in the GenBank. The sequence consists of 1,272 bp comprised of 5 introns. The correspondent protein, named Cyanophopsin, showed high identity to rhodopsin-like proteins of Archea, Bacteria, Fungi, and Algae. At the N-terminal, the protein is characterized by a region with no transmembrane alpha-helices (80 aa), followed by a region with 7alpha-helices (219 aa) and a shorter 35-aa C-terminal region. The DNA sequence of the N-terminal region was expressed in E. coli and the recombinant purified peptide was used as antigen in hens to obtain polyclonal antibodies. Indirect immunofluorescence in C. paradoxa cells showed a marked labeling of the muroplast (aka cyanelle) membrane.

Evolutionary dynamics of U12-type spliceosomal introns

Background

Many multicellular eukaryotes have two types of spliceosomes for the removal of introns from messenger RNA precursors. The major (U2) spliceosome processes the vast majority of introns, referred to as U2-type introns, while the minor (U12) spliceosome removes a small fraction (less than 0.5%) of introns, referred to as U12-type introns. U12-type introns have distinct sequence elements and usually occur together in genes with U2-type introns. A phylogenetic distribution of U12-type introns shows that the minor splicing pathway appeared very early in eukaryotic evolution and has been lost repeatedly.

Results

We have investigated the evolution of U12-type introns among eighteen metazoan genomes by analyzing orthologous U12-type intron clusters. Examination of gain, loss, and type switching shows that intron type is remarkably conserved among vertebrates. Among 180 intron clusters, only eight show intron loss in any vertebrate species and only five show conversion between the U12 and the U2-type. Although there are only nineteen U12-type introns in Drosophila melanogaster, we found one case of U2 to U12-type conversion, apparently mediated by the activation of cryptic U12 splice sites early in the dipteran lineage. Overall, loss of U12-type introns is more common than conversion to U2-type and the U12 to U2 conversion occurs more frequently among introns of the GT-AG subtype than among introns of the AT-AC subtype. We also found support for natural U12-type introns with non-canonical terminal dinucleotides (CT-AC, GG-AG, and GA-AG) that have not been previously reported.

Conclusions

Although complete loss of the U12-type spliceosome has occurred repeatedly, U12 introns are extremely stable in some taxa, including eutheria. Loss of U12 introns or the genes containing them is more common than conversion to the U2-type. The degeneracy of U12-type terminal dinucleotides among natural U12-type introns is higher than previously thought.