Flavonols improve thermotolerance in tomato pollen during germination and tube elongation by maintaining ROS homeostasis.

Elevated temperatures impair pollen performance and reproductive success, resulting in lower crop yields. The Solanum lycopersicum anthocyanin reduced ( are ) mutant has a FLAVANONE 3 HYDROXYLASE ( F3H ) gene mutation resulting in impaired synthesis of flavonol antioxidants. The are mutant has reduced pollen performance and seed set relative to the VF36 parental line, which is accentuated at elevated temperatures. Transformation of are with the wild-type F3H gene, or chemical complementation with flavonols, prevented temperature-dependent ROS accumulation in pollen and reversed are's reduced viability, germination, and tube elongation to VF36 levels. VF36 transformed with an F3H overexpression construct prevented temperature driven ROS increases and impaired pollen performance, revealing thermotolerance results from elevated flavonol synthesis. Although stigmas of are had reduced flavonols and elevated ROS, the growth of are pollen tubes were similarly impaired in both are and VF36 pistils. RNA-Seq was performed at optimal and stress temperatures in are , VF36, and the VF36 F3H overexpression line at multiple timepoints across pollen tube elongation. Differentially expressed gene numbers increased with duration of elevated temperature in all genotypes, with the largest number in are . These findings suggest potential agricultural interventions to combat the negative effects of heat-induced ROS in pollen that leads to reproductive failure. One sentence summary: Flavonol antioxidants reduce the negative impacts of elevated temperatures on pollen performance by reducing levels of heat induced reactive oxygen species and modulation of heat-induced changes in the pollen transcriptome.

Genomic targets of HOP2 are enriched for features found at recombination hotspots.

HOP2 is a conserved protein that plays a positive role in homologous chromosome pairing and a separable role in preventing illegitimate connections between nonhomologous chromosome regions during meiosis. We employed ChIP-seq to discover that Arabidopsis HOP2 binds along the length of all chromosomes, except for centromeric and nucleolar organizer regions, and no binding sites were detected in the organelle genomes. A large number of reads were assigned to the HOP2 locus itself, yet TAIL-PCR and SNP analysis of the aligned sequences indicate that many of these reads originate from the transforming T-DNA, supporting the role of HOP2 in preventing nonhomologous exchanges. The 292 ChIP-seq peaks are largely found in promoter regions and downstream from genes, paralleling the distribution of recombination hotspots, and motif analysis revealed that there are several conserved sequences that are also enriched at crossover sites. We conducted coimmunoprecipitation of HOP2 followed by LC-MS/MS and found enrichment for several proteins, including some histone variants and modifications that are also known to be associated with recombination hotspots. We propose that HOP2 may be directed to chromatin motifs near double strand breaks, where homology checks are proposed to occur.

BioViz Connect: Web Application Linking CyVerse Cloud Resources to Genomic Visualization in the Integrated Genome Browser.

Genomics researchers do better work when they can interactively explore and visualize data. Due to the vast size of experimental datasets, researchers are increasingly using powerful, cloud-based systems to process and analyze data. These remote systems, called science gateways, offer user-friendly, Web-based access to high performance computing and storage resources, but typically lack interactive visualization capability. In this paper, we present BioViz Connect, a middleware Web application that links CyVerse science gateway resources to the Integrated Genome Browser (IGB), a highly interactive native application implemented in Java that runs on the user's personal computer. Using BioViz Connect, users can 1) stream data from the CyVerse data store into IGB for visualization, 2) improve the IGB user experience for themselves and others by adding IGB specific metadata to CyVerse data files, including genome version and track appearance, and 3) run compute-intensive visual analytics functions on CyVerse infrastructure to create new datasets for visualization in IGB or other applications. To demonstrate how BioViz Connect facilitates interactive data visualization, we describe an example RNA-Seq data analysis investigating how heat and desiccation stresses affect gene expression in the model plant Arabidopsis thaliana. The RNA-Seq use case illustrates how interactive visualization with IGB can help a user identify problematic experimental samples, sanity-check results using a positive control, and create new data files for interactive visualization in IGB (or other tools) using a Docker image deployed to CyVerse via the Terrain API. Lastly, we discuss limitations of the technologies used and suggest opportunities for future work. BioViz Connect is available from https://bioviz.org.

Integrated Genome Browser App Store.

SUMMARY: Rapid progress in genome science requires equally rapid visualization software development so that researchers can better explore and understand novel datasets. To make developing new visualizations faster and easier, we previously re-factored the Integrated Genome Browser (IGB), a desktop Java application with dozens of features, into a pluggable application framework that can accept new functionality as plug-ins, called IGB Apps. However, developers lacked a centralized location for sharing Apps, making it hard to connect with potential users. To fill this gap, we created an App Store for IGB, a user-friendly Web site for developers to release and document Apps, and for users to find them. AVAILABILITY AND IMPLEMENTATION: The IGB App Store is available from https://bioviz.org.

AINTEGUMENTA and AINTEGUMENTA-LIKE6 directly regulate floral homeotic, growth, and vascular development genes in young Arabidopsis flowers.

Arabidopsis flower primordia give rise to organ primordia in stereotypical positions within four concentric whorls. Floral organ primordia in each whorl undergo distinct developmental programs to become one of four organ types (sepals, petals, stamens, and carpels). The Arabidopsis transcription factors AINTEGUMENTA (ANT) and AINTEGUMENTA-LIKE6 (AIL6) are required for correct positioning of floral organ initiation, contribute to the specification of floral organ identity, and regulate the growth and morphogenesis of developing floral organs. To gain insight into the molecular mechanisms by which ANT and AIL6 contribute to floral organogenesis, we identified the genome-wide binding sites of both ANT and AIL6 in stage 3 flower primordia, the developmental stage at which sepal primordia become visible and class B and C floral homeotic genes are first expressed. AIL6 binds to a subset of ANT sites, suggesting that AIL6 regulates some but not all of the same target genes as ANT. ANT- and AIL6-binding sites are associated with genes involved in many biological processes related to meristem and flower organ development. Comparison of genes associated with both ANT and AIL6 ChIP-Seq peaks and those differentially expressed after perturbation of ANT and/or AIL6 activity identified likely direct targets of ANT and AIL6 regulation. These include class B and C floral homeotic genes, growth regulatory genes, and genes involved in vascular development.

Anno genominis XX: 20 years of Arabidopsis genomics.

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

The Arabidopsis transcription factor AINTEGUMENTA orchestrates patterning genes and auxin signaling in the establishment of floral growth and form.

Understanding how flowers form is an important problem in plant biology, as human food supply depends on flower and seed production. Flower development also provides an excellent model for understanding how cell division, expansion and differentiation are coordinated during organogenesis. In the model plant Arabidopsis thaliana, floral organogenesis requires AINTEGUMENTA (ANT) and AINTEGUMENTA-LIKE 6 (AIL6)/PLETHORA 3 (PLT3), two members of the Arabidopsis AINTEGUMENTA-LIKE/PLETHORA (AIL/PLT) transcription factor family. Together, ANT and AIL6/PLT3 regulate aspects of floral organogenesis, including floral organ initiation, growth, identity specification and patterning. Previously, we used RNA-Seq to identify thousands of genes with disrupted expression in ant ail6 mutant flowers, indicating that ANT and AIL6/PLT3 influence a vast transcriptional network. The immediate downstream targets of ANT and AIL6/PLT3 in flowers are unknown, however. To identify direct targets of ANT regulation, we performed an RNA-Seq time-course experiment in which we induced ANT activity in transgenic plants bearing an ANT-glucocorticoid receptor fusion construct. In addition, we performed a ChIP-Seq experiment that identified ANT binding sites in developing flowers. These experiments identified 200 potential ANT target genes based on their proximity to ANT binding sites and differential expression in response to ANT. These 200 candidate target genes were involved in functions such as polarity specification, floral organ development, meristem development and auxin signaling. In addition, we identified several genes associated with lateral organ growth that may mediate the role of ANT in organ size control. These results reveal new features of the ANT transcriptional network by linking ANT to previously unknown regulatory targets.

An 'eFP-Seq Browser' for visualizing and exploring RNA sequencing data.

Improvements in next-generation sequencing technologies have resulted in dramatically reduced sequencing costs. This has led to an explosion of '-seq'-based methods, of which RNA sequencing (RNA-seq) for generating transcriptomic data is the most popular. By analysing global patterns of gene expression in organs/tissues/cells of interest or in response to chemical or environmental perturbations, researchers can better understand an organism's biology. Tools designed to work with large RNA-seq data sets enable analyses and visualizations to help generate hypotheses about a gene's function. We present here a user-friendly RNA-seq data exploration tool, called the 'eFP-Seq Browser', that shows the read map coverage of a gene of interest in each of the samples along with 'electronic fluorescent pictographic' (eFP) images that serve as visual representations of expression levels. The tool also summarizes the details of each RNA-seq experiment, providing links to archival databases and publications. It automatically computes the reads per kilobase per million reads mapped expression-level summaries and point biserial correlation scores to sort the samples based on a gene's expression level or by how dissimilar the read map profile is from a gene splice variant, to quickly identify samples with the strongest expression level or where alternative splicing might be occurring. Links to the Integrated Genome Browser desktop visualization tool allow researchers to visualize and explore the details of RNA-seq alignments summarized in eFP-Seq Browser as coverage graphs. We present four cases of use of the eFP-Seq Browser for ABI3, SR34, SR45a and U2AF65B, where we examine expression levels and identify alternative splicing. The URL for the browser is https://bar.utoronto.ca/eFP-Seq_Browser/. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. Tool is at https://bar.utoronto.ca/eFP-Seq_Browser/; RNA-seq data at https://s3.amazonaws.com/iplant-cdn/iplant/home/araport/rnaseq_bam/ and https://s3.amazonaws.com/iplant-cdn/iplant/home/araport/rnaseq_bam/Klepikova/. Code is available at https://github.com/BioAnalyticResource/eFP-Seq-Browser.

Many rice genes are differentially spliced between roots and shoots but cytokinin has minimal effect on splicing.

Alternatively spliced genes produce multiple spliced isoforms, called transcript variants. In differential alternative splicing, transcript variant abundance differs across sample types. Differential alternative splicing is common in animal systems and influences cellular development in many processes, but its extent and significance is not as well known in plants. To investigate differential alternative splicing in plants, we examined RNA-Seq data from rice seedlings. The data included three biological replicates per sample type, approximately 30 million sequence alignments per replicate, and four sample types: roots and shoots treated with exogenous cytokinin delivered hydroponically or a mock treatment. Cytokinin treatment triggered expression changes in thousands of genes but had negligible effect on splicing patterns. However, many genes were differentially spliced between mock-treated roots and shoots, indicating that our methods were sufficiently sensitive to detect differential splicing between data sets. Quantitative fragment analysis of reverse transcriptase-PCR products made from newly prepared rice samples confirmed 9 of 10 differential splicing events between rice roots and shoots. Differential alternative splicing typically changed the relative abundance of splice variants that co-occurred in a data set. Analysis of a similar (but less deeply sequenced) RNA-Seq data set from Arabidopsis showed the same pattern. In both the Arabidopsis and rice RNA-Seq data sets, most genes annotated as alternatively spliced had small minor variant frequencies. Of splicing choices with abundant support for minor forms, most alternative splicing events were located within the protein-coding sequence and maintained the annotated reading frame. A tool for visualizing protein annotations in the context of genomic sequence (ProtAnnot) together with a genome browser (Integrated Genome Browser) were used to visualize and assess effects of differential splicing on gene function. In general, differentially spliced regions coincided with conserved protein domains, indicating that differential alternative splicing is likely to affect protein function between root and shoot tissue in rice.

Coordination of Chloroplast Development through the Action of the GNC and GLK Transcription Factor Families.

Fundamental questions regarding how chloroplasts develop from proplastids remain poorly understood despite their central importance to plant life. Two families of nuclear transcription factors, the GATA NITRATE-INDUCIBLE CARBON-METABOLISM-INVOLVED (GNC) and GOLDEN TWO-LIKE (GLK) families, have been implicated in directly and positively regulating chloroplast development. Here, we determined the degree of functional overlap between the two transcription factor families in Arabidopsis (Arabidopsis thaliana), characterizing their ability to regulate chloroplast biogenesis both alone and in concert. We determined the DNA-binding motifs for GNC and GLK2 using protein-binding microarrays; the enrichment of these motifs in transcriptome datasets indicates that GNC and GLK2 are repressors and activators of gene expression, respectively. ChIP-seq analysis of GNC identified PHYTOCHROME INTERACTING FACTOR and brassinosteroid activity genes as targets whose repression by GNC facilitates chloroplast biogenesis. In addition, GNC targets and represses genes involved in ERECTA signaling and thereby facilitates stomatal development. Our results define key regulatory features of the GNC and GLK transcription factor families that contribute to the control of chloroplast biogenesis and photosynthetic activity, including areas of independence and cross talk.

Cell-specific cis-natural antisense transcripts (cis-NATs) in the sperm and the pollen vegetative cells of Arabidopsis thaliana.

Background: cis-NATs (cis-natural antisense transcripts ) are transcribed from opposite strands of adjacent genes and have been shown to regulate gene expression by generating small RNAs from the overlapping region. cis-NATs are important for plant development and resistance to pathogens and stress. Several genome-wide investigations identified a number of cis-NAT pairs, but these investigations predicted cis-NATS using expression data from bulk samples that included lots of cell types. Some cis-NAT pairs identified from those investigations might not be functional, because both transcripts of cis-NAT pairs need to be co-expressed in the same cell. Pollen only contains two cell types, two sperm and one vegetative cell, which makes cell-specific investigation of cis-NATs possible. Methods: We investigated potential protein-coding cis-NATs in pollen and sperm using pollen RNA-seq data and TAIR10 gene models using the Integrated Genome Browser.  We then used sperm microarray data and sRNAs in sperm and pollen to determine possibly functional cis-NATs in the sperm or vegetative cell, respectively. Results: We identified 1471 potential protein-coding cis-NAT pairs, including 131 novel pairs that were not present in TAIR10 gene models. In pollen, 872 possibly functional pairs were identified. 72 and 56 pairs were potentially functional in sperm and vegetative cells, respectively. sRNAs were detected at 794 genes, belonging to 739 pairs. Conclusion: These potential candidates in sperm and the vegetative cell are tools for understanding gene expression mechanisms in pollen.

Genome-wide characterization of differential transcript usage in Arabidopsis thaliana.

Alternative splicing and the usage of alternate transcription start- or stop sites allows a single gene to produce multiple transcript isoforms. Most plant genes express certain isoforms at a significantly higher level than others, but under specific conditions this expression dominance can change, resulting in a different set of dominant isoforms. These events of differential transcript usage (DTU) have been observed for thousands of Arabidopsis thaliana, Zea mays and Vitis vinifera genes, and have been linked to development and stress response. However, neither the characteristics of these genes, nor the implications of DTU on their protein coding sequences or functions, are currently well understood. Here we present a dataset of isoform dominance and DTU for all genes in the AtRTD2 reference transcriptome based on a protocol that was benchmarked on simulated data and validated through comparison with a published reverse transciptase-polymerase chain reaction panel. We report DTU events for 8148 genes across 206 public RNA-Seq samples, and find that protein sequences are affected in 22% of the cases. The observed DTU events show high consistency across replicates, and reveal reproducible patterns in response to treatment and development. We also demonstrate that genes with different evolutionary ages, expression breadths and functions show large differences in the frequency at which they undergo DTU, and in the effect that these events have on their protein sequences. Finally, we showcase how the generated dataset can be used to explore DTU events for genes of interest or to find genes with specific DTU in samples of interest.

Cytokinin induces genome-wide binding of the type-B response regulator ARR10 to regulate growth and development in Arabidopsis.

The plant hormone cytokinin affects a diverse array of growth and development processes and responses to the environment. How a signaling molecule mediates such a diverse array of outputs and how these response pathways are integrated with other inputs remain fundamental questions in plant biology. To this end, we characterized the transcriptional network initiated by the type-B ARABIDOPSIS RESPONSE REGULATORs (ARRs) that mediate the cytokinin primary response, making use of chromatin immunoprecipitation sequencing (ChIP-seq), protein-binding microarrays, and transcriptomic approaches. By ectopic overexpression of ARR10, Arabidopsis lines hypersensitive to cytokinin were generated and used to clarify the role of cytokinin in regulation of various physiological responses. ChIP-seq was used to identify the cytokinin-dependent targets for ARR10, thereby defining a crucial link between the cytokinin primary-response pathway and the transcriptional changes that mediate physiological responses to this phytohormone. Binding of ARR10 was induced by cytokinin with binding sites enriched toward the transcriptional start sites for both induced and repressed genes. Three type-B ARR DNA-binding motifs, determined by use of protein-binding microarrays, were enriched at ARR10 binding sites, confirming their physiological relevance. WUSCHEL was identified as a direct target of ARR10, with its cytokinin-enhanced expression resulting in enhanced shooting in tissue culture. Results from our analyses shed light on the physiological role of the type-B ARRs in regulating the cytokinin response, mechanism of type-B ARR activation, and basis by which cytokinin regulates diverse aspects of growth and development as well as responses to biotic and abiotic factors.

Characterization of the cytokinin-responsive transcriptome in rice.

BACKGROUND: Cytokinin activates transcriptional cascades important for development and the responses to biotic and abiotic stresses. Most of what is known regarding cytokinin-regulated gene expression comes from studies of the dicotyledonous plant Arabidopsis thaliana. To expand the understanding of the cytokinin-regulated transcriptome, we employed RNA-Seq to analyze gene expression in response to cytokinin in roots and shoots of the monocotyledonous plant rice. RESULTS: We identified over 4,600 and approximately 2,400 genes differentially expressed in response to cytokinin in roots and shoots respectively. There were some similarities in the sets of cytokinin-regulated genes identified in rice and Arabidopsis, including an up-regulation of genes that act to reduce cytokinin function. Consistent with this, we found that the preferred DNA-binding motif of a rice type-B response regulator is similar to those from Arabidopsis. Analysis of the genes regulated by cytokinin in rice revealed a large number of transcription factors, receptor-like kinases, and genes involved in protein degradation, as well as genes involved in development and the response to biotic stress. Consistent with the over-representation of genes involved in biotic stress, there is a substantial overlap in the genes regulated by cytokinin and those differentially expressed in response to pathogen infection, suggesting that cytokinin plays an integral role in the transcriptional response to pathogens in rice, including the induction of a large number of WRKY transcription factors. CONCLUSIONS: These results begin to unravel the complex gene regulation after cytokinin perception in a crop of agricultural importance and provide insight into the processes and responses modulated by cytokinin in monocots.

RNA-Seq Links the Transcription Factors AINTEGUMENTA and AINTEGUMENTA-LIKE6 to Cell Wall Remodeling and Plant Defense Pathways.

AINTEGUMENTA (ANT) and AINTEGUMENTA-LIKE6 (AIL6) are two related transcription factors in Arabidopsis (Arabidopsis thaliana) that have partially overlapping roles in several aspects of flower development, including floral organ initiation, identity specification, growth, and patterning. To better understand the biological processes regulated by these two transcription factors, we performed RNA sequencing (RNA-Seq) on ant ail6 double mutants. We identified thousands of genes that are differentially expressed in the double mutant compared with the wild type. Analyses of these genes suggest that ANT and AIL6 regulate floral organ initiation and growth through modifications to the cell wall polysaccharide pectin. We found reduced levels of demethylesterified homogalacturonan and altered patterns of auxin accumulation in early stages of ant ail6 flower development. The RNA-Seq experiment also revealed cross-regulation of AIL gene expression at the transcriptional level. The presence of a number of overrepresented Gene Ontology terms related to plant defense in the set of genes differentially expressed in ant ail6 suggest that ANT and AIL6 also regulate plant defense pathways. Furthermore, we found that ant ail6 plants have elevated levels of two defense hormones: salicylic acid and jasmonic acid, and show increased resistance to the bacterial pathogen Pseudomonas syringae These results suggest that ANT and AIL6 regulate biological pathways that are critical for both development and defense.

ProtAnnot: an App for Integrated Genome Browser to display how alternative splicing and transcription affect proteins.

UNLABELLED: One gene can produce multiple transcript variants encoding proteins with different functions. To facilitate visual analysis of transcript variants, we developed ProtAnnot, which shows protein annotations in the context of genomic sequence. ProtAnnot searches InterPro and displays profile matches (protein annotations) alongside gene models, exposing how alternative promoters, splicing and 3' end processing add, remove, or remodel functional motifs. To draw attention to these effects, ProtAnnot color-codes exons by frame and displays a cityscape graphic summarizing exonic sequence at each position. These techniques make visual analysis of alternative transcripts faster and more convenient for biologists. AVAILABILITY AND IMPLEMENTATION: ProtAnnot is a plug-in App for Integrated Genome Browser, an open source desktop genome browser available from http://www.bioviz.org CONTACT: aloraine@uncc.edu.

Integrated genome browser: visual analytics platform for genomics.

MOTIVATION: Genome browsers that support fast navigation through vast datasets and provide interactive visual analytics functions can help scientists achieve deeper insight into biological systems. Toward this end, we developed Integrated Genome Browser (IGB), a highly configurable, interactive and fast open source desktop genome browser. RESULTS: Here we describe multiple updates to IGB, including all-new capabilities to display and interact with data from high-throughput sequencing experiments. To demonstrate, we describe example visualizations and analyses of datasets from RNA-Seq, ChIP-Seq and bisulfite sequencing experiments. Understanding results from genome-scale experiments requires viewing the data in the context of reference genome annotations and other related datasets. To facilitate this, we enhanced IGB's ability to consume data from diverse sources, including Galaxy, Distributed Annotation and IGB-specific Quickload servers. To support future visualization needs as new genome-scale assays enter wide use, we transformed the IGB codebase into a modular, extensible platform for developers to create and deploy all-new visualizations of genomic data. AVAILABILITY AND IMPLEMENTATION: IGB is open source and is freely available from http://bioviz.org/igb CONTACT: aloraine@uncc.edu.

A Comparison of transgenic and wild type soybean seeds: analysis of transcriptome profiles using RNA-Seq.

BACKGROUND: Soybean (Glycine max) has been bred for thousands of years to produce seeds rich in protein for human and animal consumption, making them an appealing bioreactor for producing valuable recombinant proteins at high levels. However, the effects of expressing recombinant protein at high levels on bean physiology are not well understood. To address this, we investigated whether gene expression within transgenic soybean seed tissue is altered when large amounts of recombinant proteins are being produced and stored exclusively in the seeds. We used RNA-Seq to survey gene expression in three transgenic soybean lines expressing recombinant protein at levels representing up to 1.61 % of total protein in seed tissues. The three lines included: ST77, expressing human thyroglobulin protein (hTG), ST111, expressing human myelin basic protein (hMBP), and 764, expressing a mutant, nontoxic form of a staphylococcal subunit vaccine protein (mSEB). All lines selected for analysis were homozygous and contained a single copy of the transgene. METHODS: Each transgenic soybean seed was screened for transgene presence and recombinant protein expression via PCR and western blotting.  Whole seed mRNA was extracted and cDNA libraries constructed for Illumina sequencing.  Following alignment to the soybean reference genome, differential gene expression analysis was conducted using edgeR and cufflinks.  Functional analysis of differentially expressed genes was carried out using the gene ontology analysis tool AgriGO. RESULTS: The transcriptomes of nine seeds from each transgenic line were sequenced and compared with wild type seeds. Native soybean gene expression was significantly altered in line 764 (mSEB) with more than 3000 genes being upregulated or downregulated. ST77 (hTG) and ST111 (hMBP) had significantly less differences with 52 and 307 differentially expressed genes respectively. Gene ontology enrichment analysis found that the upregulated genes in the 764 line were annotated with functions related to endopeptidase inhibitors and protein synthesis, but suppressed expression of genes annotated to the nuclear pore and to protein transport. No significant gene ontology terms were detected in ST77, and only a few genes involved in photosynthesis and thylakoid functions were downregulated in ST111. Despite these differences, transgenic plants and seeds appeared phenotypically similar to non-transgenic controls. There was no correlation between recombinant protein expression level and the quantity of differentially expressed genes detected. CONCLUSIONS: Measurable unscripted gene expression changes were detected in the seed transcriptomes of all three transgenic soybean lines analyzed, with line 764 being substantially altered. Differences detected at the transcript level may be due to T-DNA insert locations, random mutations following transformation or direct effects of the recombinant protein itself, or a combination of these. The physiological consequences of such changes remain unknown.

Analysis of pollen-specific alternative splicing in Arabidopsis thaliana via semi-quantitative PCR.

Alternative splicing enables a single gene to produce multiple mRNA isoforms by varying splice site selection. In animals, alternative splicing of mRNA isoforms between cell types is widespread and supports cellular differentiation. In plants, at least 20% of multi-exon genes are alternatively spliced, but the extent and significance of tissue-specific splicing is less well understood, partly because it is difficult to isolate cells of a single type. Pollen is a useful model system to study tissue-specific splicing in higher plants because pollen grains contain only two cell types and can be collected in large amounts without damaging cells. Previously, we identified pollen-specific splicing patterns by comparing RNA-Seq data from Arabidopsis pollen and leaves. Here, we used semi-quantitative PCR to validate pollen-specific splicing patterns among genes where RNA-Seq data analysis indicated splicing was most different between pollen and leaves. PCR testing confirmed eight of nine alternative splicing patterns, and results from the ninth were inconclusive. In four genes, alternative transcriptional start sites coincided with alternative splicing. This study highlights the value of the low-cost PCR assay as a method of validating RNA-Seq results.

RNA-Seq analysis and annotation of a draft blueberry genome assembly identifies candidate genes involved in fruit ripening, biosynthesis of bioactive compounds, and stage-specific alternative splicing.

BACKGROUND: Blueberries are a rich source of antioxidants and other beneficial compounds that can protect against disease. Identifying genes involved in synthesis of bioactive compounds could enable the breeding of berry varieties with enhanced health benefits. RESULTS: Toward this end, we annotated a previously sequenced draft blueberry genome assembly using RNA-Seq data from five stages of berry fruit development and ripening. Genome-guided assembly of RNA-Seq read alignments combined with output from ab initio gene finders produced around 60,000 gene models, of which more than half were similar to proteins from other species, typically the grape Vitis vinifera. Comparison of gene models to the PlantCyc database of metabolic pathway enzymes identified candidate genes involved in synthesis of bioactive compounds, including bixin, an apocarotenoid with potential disease-fighting properties, and defense-related cyanogenic glycosides, which are toxic. Cyanogenic glycoside (CG) biosynthetic enzymes were highly expressed in green fruit, and a candidate CG detoxification enzyme was up-regulated during fruit ripening. Candidate genes for ethylene, anthocyanin, and 400 other biosynthetic pathways were also identified. Homology-based annotation using Blast2GO and InterPro assigned Gene Ontology terms to around 15,000 genes. RNA-Seq expression profiling showed that blueberry growth, maturation, and ripening involve dynamic gene expression changes, including coordinated up- and down-regulation of metabolic pathway enzymes and transcriptional regulators. Analysis of RNA-seq alignments identified developmentally regulated alternative splicing, promoter use, and 3' end formation. CONCLUSIONS: We report genome sequence, gene models, functional annotations, and RNA-Seq expression data that provide an important new resource enabling high throughput studies in blueberry.