ASpli: integrative analysis of splicing landscapes through RNA-Seq assays.

MOTIVATION: Genome-wide analysis of alternative splicing has been a very active field of research since the early days of next generation sequencing technologies. Since then, ever-growing data availability and the development of increasingly sophisticated analysis methods have uncovered the complexity of the general splicing repertoire. A large number of splicing analysis methodologies exist, each of them presenting its own strengths and weaknesses. For instance, methods exclusively relying on junction information do not take advantage of the large majority of reads produced in an RNA-seq assay, isoform reconstruction methods might not detect novel intron retention events, some solutions can only handle canonical splicing events, and many existing methods can only perform pairwise comparisons. RESULTS: In this contribution, we present ASpli, a computational suite implemented in R statistical language, that allows the identification of changes in both, annotated and novel alternative-splicing events and can deal with simple, multi-factor or paired experimental designs. Our integrative computational workflow, that considers the same GLM model applied to different sets of reads and junctions, allows computation of complementary splicing signals. Analyzing simulated and real data, we found that the consolidation of these signals resulted in a robust proxy of the occurrence of splicing alterations. While the analysis of junctions allowed us to uncover annotated as well as non-annotated events, read coverage signals notably increased recall capabilities at a very competitive performance when compared against other state-of-the-art splicing analysis algorithms. AVAILABILITY AND IMPLEMENTATION: ASpli is freely available from the Bioconductor project site https://doi.org/doi:10.18129/B9.bioc.ASpli. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Global transcriptome analysis reveals circadian control of splicing events in Arabidopsis thaliana.

The circadian clock of Arabidopsis thaliana controls many physiological and molecular processes, allowing plants to anticipate daily changes in their environment. However, developing a detailed understanding of how oscillations in mRNA levels are connected to oscillations in co/post-transcriptional processes, such as splicing, has remained a challenge. Here we applied a combined approach using deep transcriptome sequencing and bioinformatics tools to identify novel circadian-regulated genes and splicing events. Using a stringent approach, we identified 300 intron retention, eight exon skipping, 79 alternative 3' splice site usage, 48 alternative 5' splice site usage, and 350 multiple (more than one event type) annotated events under circadian regulation. We also found seven and 721 novel alternative exonic and intronic events. Depletion of the circadian-regulated splicing factor AtSPF30 homologue resulted in the disruption of a subset of clock-controlled splicing events. Altogether, our global circadian RNA-seq coupled with an in silico, event-centred, splicing analysis tool offers a new approach for studying the interplay between the circadian clock and the splicing machinery at a global scale. The identification of many circadian-regulated splicing events broadens our current understanding of the level of control that the circadian clock has over this co/post-transcriptional regulatory layer.

Improvement of alfalfa forage quality and management through the down-regulation of MsFTa1.

Alfalfa (Medicago sativa L.) is one of the most important forage crops worldwide. As a perennial, alfalfa is cut several times each year. Farmers face a dilemma: if cut earlier, forage nutritive value is much higher but regrowth is affected and the longevity of the stand is severely compromised. On the other hand, if alfalfa is cut later at full flower, stands persist longer and more biomass may be harvested, but the nutritive value diminishes. Alfalfa is a strict long-day plant. We reasoned that by manipulating the response to photoperiod, we could delay flowering to improve forage quality and widen each harvesting window, facilitating management. With this aim, we functionally characterized the FLOWERING LOCUS T family of genes, represented by five members: MsFTa1, MsFTa2, MsFTb1, MsFTb2 and MsFTc. The expression of MsFTa1 correlated with photoperiodic flowering and its down-regulation led to severe delayed flowering. Altogether, with late flowering, low expression of MsFTa1 led to changes in plant architecture resulting in increased leaf to stem biomass ratios and forage digestibility. By manipulating photoperiodic flowering, we were able to improve the quality of alfalfa forage and management, which may allow farmers to cut alfalfa of high nutritive value without compromising stand persistence.

SPF45-related splicing factor for phytochrome signaling promotes photomorphogenesis by regulating pre-mRNA splicing in Arabidopsis.

Light signals regulate plant growth and development by controlling a plethora of gene expression changes. Posttranscriptional regulation, especially pre-mRNA processing, is a key modulator of gene expression; however, the molecular mechanisms linking pre-mRNA processing and light signaling are not well understood. Here we report a protein related to the human splicing factor 45 (SPF45) named splicing factor for phytochrome signaling (SFPS), which directly interacts with the photoreceptor phytochrome B (phyB). In response to light, SFPS-RFP (red fluorescent protein) colocalizes with phyB-GFP in photobodies. sfps loss-of-function plants are hyposensitive to red, far-red, and blue light, and flower precociously. SFPS colocalizes with U2 small nuclear ribonucleoprotein-associated factors including U2AF65B, U2A', and U2AF35A in nuclear speckles, suggesting SFPS might be involved in the 3' splice site determination. SFPS regulates pre-mRNA splicing of a large number of genes, of which many are involved in regulating light signaling, photosynthesis, and the circadian clock under both dark and light conditions. In vivo RNA immunoprecipitation (RIP) assays revealed that SFPS associates with EARLY FLOWERING 3 (ELF3) mRNA, a critical link between light signaling and the circadian clock. Moreover, PHYTOCHROME INTERACTING FACTORS (PIFs) transcription factor genes act downstream of SFPS, as the quadruple pif mutant pifq suppresses defects of sfps mutants. Taken together, these data strongly suggest SFPS modulates light-regulated developmental processes by controlling pre-mRNA splicing of light signaling and circadian clock genes.

Combinatorial control of adhesion of Brucella abortus 2308 to host cells by transcriptional rewiring of the trimeric autotransporter btaE gene.

Regulatory network plasticity is a key attribute underlying changes in bacterial gene expression and a source of phenotypic diversity to interact with the surrounding environment. Here, we sought to study the transcriptional circuit of HutC, a regulator of both metabolic and virulence genes of the facultative intracellular pathogen Brucella. Using in silico and biochemical approaches, we identified a novel functional HutC-binding site upstream of btaE, a trimeric-autotransporter adhesin involved in the attachment of Brucella to host extracellular matrix components. Moreover, we identified two additional regulators, one of which, MdrA, acts in concert with HutC to exert a combinatorial control of both btaE promoter activity and attachment of Brucella to HeLa cells. Analysis of btaE promoter sequences of different species indicated that this HutC-binding site was generated de novo by a single point mutation in a virulent Brucella strain, indicative of a transcriptional rewiring event. In addition to major domain organization differences existing between BtaE proteins within the genus Brucella, our analyses revealed that sequences upstream of btaE display high variability probably associated to intrinsic promoter structural features, which may serve as a substrate for reciprocal selection during co-evolution between this pathogen and its mammalian host.

The Dengue Virus NS5 Protein Intrudes in the Cellular Spliceosome and Modulates Splicing.

Dengue virus NS5 protein plays multiple functions in the cytoplasm of infected cells, enabling viral RNA replication and counteracting host antiviral responses. Here, we demonstrate a novel function of NS5 in the nucleus where it interferes with cellular splicing. Using global proteomic analysis of infected cells together with functional studies, we found that NS5 binds spliceosome complexes and modulates endogenous splicing as well as minigene-derived alternative splicing patterns. In particular, we show that NS5 alone, or in the context of viral infection, interacts with core components of the U5 snRNP particle, CD2BP2 and DDX23, alters the inclusion/exclusion ratio of alternative splicing events, and changes mRNA isoform abundance of known antiviral factors. Interestingly, a genome wide transcriptome analysis, using recently developed bioinformatics tools, revealed an increase of intron retention upon dengue virus infection, and viral replication was improved by silencing specific U5 components. Different mechanistic studies indicate that binding of NS5 to the spliceosome reduces the efficiency of pre-mRNA processing, independently of NS5 enzymatic activities. We propose that NS5 binding to U5 snRNP proteins hijacks the splicing machinery resulting in a less restrictive environment for viral replication.

The spliceosome assembly factor GEMIN2 attenuates the effects of temperature on alternative splicing and circadian rhythms.

The mechanisms by which poikilothermic organisms ensure that biological processes are robust to temperature changes are largely unknown. Temperature compensation, the ability of circadian rhythms to maintain a relatively constant period over the broad range of temperatures resulting from seasonal fluctuations in environmental conditions, is a defining property of circadian networks. Temperature affects the alternative splicing (AS) of several clock genes in fungi, plants, and flies, but the splicing factors that modulate these effects to ensure clock accuracy throughout the year remain to be identified. Here we show that GEMIN2, a spliceosomal small nuclear ribonucleoprotein assembly factor conserved from yeast to humans, modulates low temperature effects on a large subset of pre-mRNA splicing events. In particular, GEMIN2 controls the AS of several clock genes and attenuates the effects of temperature on the circadian period in Arabidopsis thaliana. We conclude that GEMIN2 is a key component of a posttranscriptional regulatory mechanism that ensures the appropriate acclimation of plants to daily and seasonal changes in temperature conditions.

Role for LSM genes in the regulation of circadian rhythms.

Growing evidence suggests that core spliceosomal components differentially affect RNA processing of specific genes; however, whether changes in the levels or activities of these factors control specific signaling pathways is largely unknown. Here we show that some SM-like (LSM) genes, which encode core components of the spliceosomal U6 small nuclear ribonucleoprotein complex, regulate circadian rhythms in plants and mammals. We found that the circadian clock regulates the expression of LSM5 in Arabidopsis plants and several LSM genes in mouse suprachiasmatic nucleus. Further, mutations in LSM5 or LSM4 in Arabidopsis, or down-regulation of LSM3, LSM5, or LSM7 expression in human cells, lengthens the circadian period. Although we identified changes in the expression and alternative splicing of some core clock genes in Arabidopsis lsm5 mutants, the precise molecular mechanism causing period lengthening remains to be identified. Genome-wide expression analysis of either a weak lsm5 or a strong lsm4 mutant allele in Arabidopsis revealed larger effects on alternative splicing than on constitutive splicing. Remarkably, large splicing defects were not observed in most of the introns evaluated using RNA-seq in the strong lsm4 mutant allele used in this study. These findings support the idea that some LSM genes play both regulatory and constitutive roles in RNA processing, contributing to the fine-tuning of specific signaling pathways.

LNK genes integrate light and clock signaling networks at the core of the Arabidopsis oscillator.

Light signaling pathways and the circadian clock interact to help organisms synchronize physiological and developmental processes with periodic environmental cycles. The plant photoreceptors responsible for clock resetting have been characterized, but signaling components that link the photoreceptors to the clock remain to be identified. Here we describe a family of night light-inducible and clock-regulated genes (LNK) that play a key role linking light regulation of gene expression to the control of daily and seasonal rhythms in Arabidopsis thaliana. A genomewide transcriptome analysis revealed that most light-induced genes respond more strongly to light during the subjective day, which is consistent with the diurnal nature of most physiological processes in plants. However, a handful of genes, including the homologous genes LNK1 and LNK2, are more strongly induced by light in the middle of the night, when the clock is most responsive to this signal. Further analysis revealed that the morning phased LNK1 and LNK2 genes control circadian rhythms, photomorphogenic responses, and photoperiodic dependent flowering, most likely by regulating a subset of clock and flowering time genes in the afternoon. LNK1 and LNK2 themselves are directly repressed by members of the TIMING OF CAB1 EXPRESSION/PSEUDO RESPONSE REGULATOR family of core-clock genes in the afternoon and early night. Thus, LNK1 and LNK2 integrate early light signals with temporal information provided by core oscillator components to control the expression of afternoon genes, allowing plants to keep track of seasonal changes in day length.

Genome wide comparative analysis of the effects of PRMT5 and PRMT4/CARM1 arginine methyltransferases on the Arabidopsis thaliana transcriptome.

BACKGROUND: Methylation at arginine residues (R) is an important post-translational modification that regulates a myriad of essential cellular processes in eukaryotes, such as transcriptional regulation, RNA processing, signal transduction and DNA repair. Arginine methylation is catalyzed by a family of enzymes known as protein arginine methyltransferases (PRMTs). PRMTs are classified as Type I or Type II, depending on the position of the methyl group on the guanidine of the methylated arginine. Previous reports have linked symmetric R methylation to transcriptional repression, while asymmetric R methylation is generally associated with transcriptional activation. However, global studies supporting this conclusion are not available. RESULTS: Here we compared side by side the physiological and molecular roles of the best characterized plant PRMTs, the Type II PRMT5 and the Type I PRMT4, also known as CARM1 in mammals. We found that prmt5 and prmt4a;4b mutants showed similar alterations in flowering time, photomorphogenic responses and salt stress tolerance, while only prmt5 mutants exhibited alterations in circadian rhythms. An RNA-seq analysis revealed that expression and splicing of many differentially regulated genes was similarly enhanced or repressed by PRMT5 and PRMT4s. Furthermore, PRMT5 and PRMT4s co-regulated the expression and splicing of key regulatory genes associated with transcription, RNA processing, responses to light, flowering, and abiotic stress tolerance, being candidates to mediate the physiological alterations observed in the mutants. CONCLUSIONS: Our global analysis indicates that two of the most important Type I and Type II arginine methyltransferases, PRTM4 and PRMT5, have mostly overlapping as well as specific, but not opposite, roles in the global regulation of gene expression in plants.

Genome-wide annotation of the soybean WRKY family and functional characterization of genes involved in response to Phakopsora pachyrhizi infection.

BACKGROUND: Many previous studies have shown that soybean WRKY transcription factors are involved in the plant response to biotic and abiotic stresses. Phakopsora pachyrhizi is the causal agent of Asian Soybean Rust, one of the most important soybean diseases. There are evidences that WRKYs are involved in the resistance of some soybean genotypes against that fungus. The number of WRKY genes already annotated in soybean genome was underrepresented. In the present study, a genome-wide annotation of the soybean WRKY family was carried out and members involved in the response to P. pachyrhizi were identified. RESULTS: As a result of a soybean genomic databases search, 182 WRKY-encoding genes were annotated and 33 putative pseudogenes identified. Genes involved in the response to P. pachyrhizi infection were identified using superSAGE, RNA-Seq of microdissected lesions and microarray experiments. Seventy-five genes were differentially expressed during fungal infection. The expression of eight WRKY genes was validated by RT-qPCR. The expression of these genes in a resistant genotype was earlier and/or stronger compared with a susceptible genotype in response to P. pachyrhizi infection. Soybean somatic embryos were transformed in order to overexpress or silence WRKY genes. Embryos overexpressing a WRKY gene were obtained, but they were unable to convert into plants. When infected with P. pachyrhizi, the leaves of the silenced transgenic line showed a higher number of lesions than the wild-type plants. CONCLUSIONS: The present study reports a genome-wide annotation of soybean WRKY family. The participation of some members in response to P. pachyrhizi infection was demonstrated. The results contribute to the elucidation of gene function and suggest the manipulation of WRKYs as a strategy to increase fungal resistance in soybean plants.

Intron detention tightly regulates the stemness/differentiation switch in the adult neurogenic niche.

The adult mammalian brain retains some capacity to replenish neurons and glia, holding promise for brain regeneration. Thus, understanding the mechanisms controlling adult neural stem cell (NSC) differentiation is crucial. Paradoxically, adult NSCs in the subependymal zone transcribe genes associated with both multipotency maintenance and neural differentiation, but the mechanism that prevents conflicts in fate decisions due to these opposing transcriptional programmes is unknown. Here we describe intron detention as such control mechanism. In NSCs, while multiple mRNAs from stemness genes are spliced and exported to the cytoplasm, transcripts from differentiation genes remain unspliced and detained in the nucleus, and the opposite is true under neural differentiation conditions. We also show that m6A methylation is the mechanism that releases intron detention and triggers nuclear export, enabling rapid and synchronized responses. m6A RNA methylation operates as an on/off switch for transcripts with antagonistic functions, tightly controlling the timing of NSCs commitment to differentiation.

Genomic deletions explain the generation of alternative BRAF isoforms conferring resistance to MAPK inhibitors in melanoma.

Resistance to MAPK inhibitors (MAPKi), the main cause of relapse in BRAF-mutant melanoma, is associated with the production of alternative BRAF mRNA isoforms (altBRAFs) in up to 30% of patients receiving BRAF inhibitor monotherapy. These altBRAFs have been described as being generated by alternative pre-mRNA splicing, and splicing modulation has been proposed as a therapeutic strategy to overcome resistance. In contrast, we report that altBRAFs are generated through genomic deletions. Using different in vitro models of altBRAF-mediated melanoma resistance, we demonstrate the production of altBRAFs exclusively from the BRAF V600E allele, correlating with corresponding genomic deletions. Genomic deletions are also detected in tumor samples from melanoma and breast cancer patients expressing altBRAFs. Along with the identification of altBRAFs in BRAF wild-type and in MAPKi-naive melanoma samples, our results represent a major shift in our understanding of mechanisms leading to the generation of BRAF transcripts variants associated with resistance in melanoma.

MsTFL1A delays flowering and regulates shoot architecture and root development in Medicago sativa.

KEY MESSAGE: MsTFL1A is an important gene involved in flowering repression in alfalfa (Medicago sativa) which conditions not only above-ground plant shoot architecture but also root development and growth. Delayed flowering is an important trait for forage species, as it allows harvesting of high-quality forage for a longer time before nutritional values decline due to plant architecture changes related to flowering onset. Despite the relevance of delayed flowering, this trait has not yet been thoroughly exploited in alfalfa. This is mainly due to its complex genetics, sensitivity to inbreeding and to the fact that delayed flowering would be only advantageous if it allowed increased forage quality without compromising seed production. To develop new delayed-flowering varieties, we have characterized the three TERMINAL FLOWERING 1 (TFL1) family of genes in alfalfa: MsTFL1A, MsTFL1B and MsTFL1C. Constitutive expression of MsTFL1A in Arabidopsis caused late flowering and changes in inflorescence architecture, indicating that MsTFL1A is the ortholog of Arabidopsis TFL1. Overexpression of MsTFL1A in alfalfa consistently led to delayed flowering in both controlled and natural field conditions, coupled to an increase in leaf/stem ratio, a common indicator of forage quality. Additionally, overexpression of MsTFL1A reduced root development, reinforcing the role of MsTFL1A not only as a flowering repressor but also as a regulator of root development.We conclude that the precise manipulation of MsTFL1A gene expression may represent a powerful tool to improve alfalfa forage quality.

Diagnosis of mitochondrial disorders applying massive pyrosequencing.

Mitochondrial disorders are a frequent cause of neurological disability affecting children and adults. Traditionally, molecular diagnosis of mitochondrial diseases was mostly accomplished by the use of Sanger sequencing and PCR-RFLP. However, there are particular drawbacks associated with the use of these methods. Recent multidisciplinary advances have led to new sequencing methods that may overcome these limitations. Our goal was to explore the use of a next generation sequencing platform in the molecular diagnosis of mitochondrial diseases reporting our findings in adult patients that present with a clinical-pathological diagnosis of a mitochondrial encephalomyopathy. Complete genomic sequences of mitochondrial DNA were obtained by 454 massive pyrosequencing from blood samples. The analysis of these sequences allowed us to identify two diagnostic pathogenic mutations and 74 homoplasmic polymorphisms, useful for obtaining high-resolution mitochondrial haplogroups. In summary, molecular diagnosis of mitochondrial disorders could be efficiently done from readily accessible samples, such as blood, with the use of a new sequencing platform.

Analysis of Alternative Splicing During the Combinatorial Response to Simultaneous Copper and Iron Deficiency in Arabidopsis Reveals Differential Events in Genes Involved in Amino Acid Metabolism.

Copper (Cu) and iron (Fe) constitute fundamental nutrients for plant biology but are often limited due to low bioavailability. Unlike responses to single Cu or Fe deprivation, the consequences of simultaneous Cu and Fe deficiency have not yet been fully deciphered. Previously, it was demonstrated that Cu and Fe deficiency applied in combination imposes transcriptome, proteome, and metabolome changes different from those triggered under each deficiency individually. Here, we evaluated the effect of alternative splicing (AS) on the transcriptome of rosette leaves under single and simultaneous Cu and Fe deficiency. Differentially spliced genes (DSGs) and differentially expressed genes (DEGs) coincided in number (2,600 approx.) although the overlapping fraction was minimal (15%). Functional annotation of changes exclusively detected under simultaneous Cu and Fe deficiency revealed that DEGs participated in general stress responses and translation, while DSGs were involved in metabolic reactions, especially amino acid biosynthesis. Interestingly, transcripts encoding central features for tryptophan (Trp) and asparagine (Asn) synthesis - two significantly altered metabolites under simultaneous Cu and Fe deficiency - underwent exclusive intron retention events under the double deficiency. However, transcript and protein amounts for these enzymes did not correlate with Trp and Asn concentration. In consequence, we propose that AS might act as a regulatory mechanism to modify the stability and/or functionality of the enzymes and therefore fine-tune amino acid production during the combinatorial response to simultaneous Cu and Fe deficiency.

Genome sequence of Sphingomonas sp. S17, isolated from an alkaline, hyperarsenic, and hypersaline volcano-associated lake at high altitude in the Argentinean Puna.

The high-altitude Andean lakes (HAAL) in the Argentinean Puna-high Andes region represent an almost unexplored ecosystem exposed to extreme conditions (high UV irradiation, hypersalinity, drastic temperature changes, desiccation, and high pH). Here we present the first genome sequence, a Sphingomonas sp., isolated from this extreme environment.

Rhythmic Behavior Is Controlled by the SRm160 Splicing Factor in Drosophila melanogaster.

Circadian clocks organize the metabolism, physiology, and behavior of organisms throughout the day-night cycle by controlling daily rhythms in gene expression at the transcriptional and post-transcriptional levels. While many transcription factors underlying circadian oscillations are known, the splicing factors that modulate these rhythms remain largely unexplored. A genome-wide assessment of the alterations of gene expression in a null mutant of the alternative splicing regulator SR-related matrix protein of 160 kDa (SRm160) revealed the extent to which alternative splicing impacts on behavior-related genes. We show that SRm160 affects gene expression in pacemaker neurons of the Drosophila brain to ensure proper oscillations of the molecular clock. A reduced level of SRm160 in adult pacemaker neurons impairs circadian rhythms in locomotor behavior, and this phenotype is caused, at least in part, by a marked reduction in period (per) levels. Moreover, rhythmic accumulation of the neuropeptide PIGMENT DISPERSING FACTOR in the dorsal projections of these neurons is abolished after SRm160 depletion. The lack of rhythmicity in SRm160-downregulated flies is reversed by a fully spliced per construct, but not by an extra copy of the endogenous locus, showing that SRm160 positively regulates per levels in a splicing-dependent manner. Our findings highlight the significant effect of alternative splicing on the nervous system and particularly on brain function in an in vivo model.