TSIS:an R package to infer alternative splicing isoform switches for time-series data

An alternative splicing isoform switch is where a pair of transcript isoforms reverse their relative expression abundances in response to external or internal stimuli. Although computational methods are available to study differential alternative splicing, few tools for detection of isoform switches exist and these are based on pairwise comparisons. Here, we provide the TSIS R package, which is the first tool for detecting significant transcript isoform switches in time-series data. The main steps of TSIS are to search for the isoform switch points in the time-series, characterize the switches and filter the results with user input parameters. All the functions are integrated into a Shiny App for ease of implementation of the analysis

Alternative splicing of barley clock genes in response to low temperature:evidence for alternative splicing conservation

Alternative splicing (AS) is a regulated mechanism that generates multiple transcripts from individual genes. It is widespread in eukaryotic genomes and provides an effective way to control gene expression. At low temperatures, AS regulates Arabidopsis clock genes through dynamic changes in the levels of productive mRNAs. We examined AS in barley clock genes to assess whether temperature-dependent AS responses also occur in a monocotyledonous crop species. We identify changes in AS of various barley core clock genes including the barley orthologues of Arabidopsis AtLHY and AtPRR7 which showed the most pronounced AS changes in response to low temperature. The AS events modulate the levels of functional and translatable mRNAs, and potentially protein levels, upon transition to cold. There is some conservation of AS events and/or splicing behaviour of clock genes between Arabidopsis and barley. In addition, novel temperature-dependent AS of the core clock gene HvPPD-H1 (a major determinant of photoperiod response and AtPRR7 orthologue) is conserved in monocots. HvPPD-H1 showed a rapid, temperature-sensitive isoform switch which resulted in changes in abundance of AS variants encoding different protein isoforms. This novel layer of low temperature control of clock gene expression, observed in two very different species, will help our understanding of plant adaptation to different environments and ultimately offer a new range of targets for plant improvement

High-Resolution RT-PCR Analysis of Alternative Barley Transcripts

Assembly of the barley genome and extensive use of RNA-seq has resulted in an abundance of gene expression data and the recognition of wide-scale production of alternatively spliced transcripts. Here, we describe in detail a high-resolution reverse transcription-PCR based panel (HR RT-PCR) that confirms the accuracy of alternatively spliced transcripts from RNA-seq and allows quantification of changes in the proportion of splice isoforms between different experimental conditions, time points, tissues, genotypes, ecotypes, and treatments. By validating a selection of barley genes, use of the panel gives confidence or otherwise to the genome-wide global changes in alternatively spliced transcripts reported by RNA-seq. This simple assay can readily be applied to perform detailed transcript isoform analysis for any gene in any species.</p

Methods to study splicing from high-throughput RNA Sequencing data

The development of novel high-throughput sequencing (HTS) methods for RNA (RNA-Seq) has provided a very powerful mean to study splicing under multiple conditions at unprecedented depth. However, the complexity of the information to be analyzed has turned this into a challenging task. In the last few years, a plethora of tools have been developed, allowing researchers to process RNA-Seq data to study the expression of isoforms and splicing events, and their relative changes under different conditions. We provide an overview of the methods available to study splicing from short RNA-Seq data. We group the methods according to the different questions they address: 1) Assignment of the sequencing reads to their likely gene of origin. This is addressed by methods that map reads to the genome and/or to the available gene annotations. 2) Recovering the sequence of splicing events and isoforms. This is addressed by transcript reconstruction and de novo assembly methods. 3) Quantification of events and isoforms. Either after reconstructing transcripts or using an annotation, many methods estimate the expression level or the relative usage of isoforms and/or events. 4) Providing an isoform or event view of differential splicing or expression. These include methods that compare relative event/isoform abundance or isoform expression across two or more conditions. 5) Visualizing splicing regulation. Various tools facilitate the visualization of the RNA-Seq data in the context of alternative splicing. In this review, we do not describe the specific mathematical models behind each method. Our aim is rather to provide an overview that could serve as an entry point for users who need to decide on a suitable tool for a specific analysis. We also attempt to propose a classification of the tools according to the operations they do, to facilitate the comparison and choice of methods.Comment: 31 pages, 1 figure, 9 tables. Small corrections adde

Improved DNase-seq protocol facilitates high resolution mapping of DNase I hypersensitive sites in roots in Arabidopsis thaliana

BACKGROUND: Identifying cis-regulatory elements is critical in understanding the direct and indirect regulatory mechanisms of gene expression. Current approaches include DNase-seq, a technique that combines sensitivity to the nonspecific endonuclease DNase I with high throughput sequencing to identify regions of regulatory DNA on a genome-wide scale. While this method was originally developed for human cell lines, later adaptations made the processing of plant tissues possible. Challenges still remain in processing recalcitrant tissues that have low DNA content. RESULTS: By removing steps requiring the use of gel agarose plugs in DNase-seq, we were able to significantly reduce the time required to perform the protocol by at least 2 days, while also making possible the processing of difficult plant tissues. We refer to this simplified protocol as DNase I SIM (for simplified in-nucleus method). We were able to successfully create DNase-seq libraries for both leaf and root tissues in Arabidopsis using DNase I SIM. CONCLUSION: This protocol simplifies and facilitates generation of DNase-seq libraries from plant tissues for high resolution mapping of DNase I hypersensitive sites. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13007-015-0087-1) contains supplementary material, which is available to authorized users

BaRTv1.0:an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq

Background: The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants.Results: A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts - BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427-433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20-28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage.Conclusion: A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.</p

Efficiency of gene silencing in \u3ci\u3eArabidopsis\u3c/i\u3e: direct inverted repeats vs. transitive RNAi vectors

We investigated the efficiency of RNA interference (RNAi) in Arabidopsis using transitive and homologous inverted repeat (hIR) vectors. hIR constructs carry self-complementary intron-spliced fragments of the target gene whereas transitive vectors have the target sequence fragment adjacent to an intron-spliced, inverted repeat of heterologous origin. Both transitive and hIR constructs facilitated specific and heritable silencing in the three genes studied (AP1 , ETTIN and TTG1 ). Both types of vectors produced a phenotypic series that phenocopied reduction of function mutants for the respective target gene. The hIR yielded up to fourfold higher proportions of events with strongly manifested reduction of function phenotypes compared to transitive RNAi. We further investigated the efficiency and potential off-target effects of AP1 silencing by both types of vectors using genome-scale microarrays and quantitative RT-PCR. The depletion of AP1 transcripts coincided with reduction of function phenotypic changes among both hIR and transitive lines and also showed similar expression patterns among differentially regulated genes. We did not detect significant silencing directed against homologous potential off-target genes when constructs were designed with minimal sequence similarity. Both hIR and transitive methods are useful tools in plant biotechnology and genomics. The choice of vector will depend on specific objectives such as cloning throughput, number of events and degree of suppression required

Expression, Splicing, and Evolution of the Myosin Gene Family in Plants

Plants possess two myosin classes, VIII and XI. The myosins XI are implicated in organelle transport, filamentous actin organization, and cell and plant growth. Due to the large size of myosin gene families, knowledge of these molecular motors remains patchy. Using deep transcriptome sequencing and bioinformatics, we systematically investigated myosin genes in two model plants, Arabidopsis (Arabidopsis thaliana) and Brachypodium (Brachypodium distachyon). We improved myosin gene models and found that myosin genes undergo alternative splicing. We experimentally validated the gene models for Arabidopsis myosin XI-K, which plays the principal role in cell interior dynamics, as well as for its Brachypodium ortholog. We showed that the Arabidopsis gene dubbed HDK (for headless derivative of myosin XI-K), which emerged through a partial duplication of the XI-K gene, is developmentally regulated. A gene with similar architecture was also found in Brachypodium. Our analyses revealed two predominant patterns of myosin gene expression, namely pollen/stamen-specific and ubiquitous expression throughout the plant. We also found that several myosins XI can be rhythmically expressed. Phylogenetic reconstructions indicate that the last common ancestor of the angiosperms possessed two myosins VIII and five myosins XI, many of which underwent additional lineage-specific duplications.This is the publisher’s final pdf. The published article is copyrighted by the American Society of Plant Biologists and can be found at: http://www.plantphysiol.org/

Genome-wide SNPs and re-sequencing of growth habit and inflorescence genes in barley: implications for association mapping in germplasm arrays varying in size and structure

<p>Abstract</p> <p>Background</p> <p>Considerations in applying association mapping (AM) to plant breeding are population structure and size: not accounting for structure and/or using small populations can lead to elevated false-positive rates. The principal determinants of population structure in cultivated barley are growth habit and inflorescence type. Both are under complex genetic control: growth habit is controlled by the epistatic interactions of several genes. For inflorescence type, multiple loss-of-function alleles in one gene lead to the same phenotype. We used these two traits as models for assessing the effectiveness of AM. This research was initiated using the CAP Core germplasm array (n = 102) assembled at the start of the Barley Coordinated Agricultural Project (CAP). This array was genotyped with 4,608 SNPs and we re-sequenced genes involved in morphology, growth and development. Larger arrays of breeding germplasm were subsequently genotyped and phenotyped under the auspices of the CAP project. This provided sets of 247 accessions phenotyped for growth habit and 2,473 accessions phenotyped for inflorescence type. Each of the larger populations was genotyped with 3,072 SNPs derived from the original set of 4,608.</p> <p>Results</p> <p>Significant associations with SNPs located in the vicinity of the loci involved in growth habit and inflorescence type were found in the CAP Core. Differentiation of true and spurious associations was not possible without <it>a priori </it>knowledge of the candidate genes, based on re-sequencing. The re-sequencing data were used to define allele types of the determinant genes based on functional polymorphisms. In a second round of association mapping, these synthetic markers based on allele types gave the most significant associations. When the synthetic markers were used as anchor points for analysis of interactions, we detected other known-function genes and candidate loci involved in the control of growth habit and inflorescence type. We then conducted association analyses - with SNP data only - in the larger germplasm arrays. For both vernalization sensitivity and inflorescence type, the most significant associations in the larger data sets were found with SNPs coincident with the synthetic markers used in the CAP Core and with SNPs detected via interaction analysis in the CAP Core.</p> <p>Conclusions</p> <p>Small and highly structured collections of germplasm, such as the CAP Core, are cost-effectively phenotyped and genotyped with high-throughput markers. They are also useful for characterizing allelic diversity at loci in germplasm of interest. Our results suggest that discovery-oriented exercises in AM in such small arrays may generate a large number of false-positives. However, if haplotypes in candidate genes are available, they may be used as anchors in an analysis of interactions to identify other candidate regions harboring genes determining target traits. Using larger germplasm arrays, genome regions where the principal genes determining vernalization sensitivity and row type are located were identified.</p