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.

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.

Two related Arabidopsis transcription factors may regulate flower development through effects on the cell wall polysaccharide pectin and act in plant defense pa

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. 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. 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.

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. 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. 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. 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.

INTRODUCTION: Microscopic colitis is a relatively common cause of chronic diarrhea and may be linked to luminal factors. Given the essential role of the microbiome in human gut health, analysis of microbiome changes associated with microscopic colitis could provide insights into the development of the disease. METHODS: We enrolled patients who underwent colonoscopy for diarrhea. An experienced pathologist classified patients as having microscopic colitis (n = 52) or controls (n = 153). Research biopsies were taken from the ascending (ASC) and descending (DES) colon, and the microbiome was characterized with Illumina sequencing. We analyzed the associations between microscopic colitis and microbiome with a series of increasingly complex models adjusted for a range of demographic and health factors. RESULTS: We found that alpha diversity was significantly lower in cases with microscopic colitis compared with that in controls in the DES colon microbiome. In the DES colon, a series of models that adjusted for an increasing number of covariates found taxa significantly associated with microscopic colitis, including Proteobacteria that was enriched in cases and Collinsella that was enriched in controls. While the alpha diversity and taxa were not significantly associated with microscopic colitis in the ASC colon microbiome, the inference P values based on ASC and DES microbiomes were highly correlated. DISCUSSION: Our study demonstrates an altered microbiome in cases with microscopic colitis compared with that in controls. Because both the cases and controls experienced diarrhea, we have identified candidate taxa that could be mechanistically responsible for the development of microscopic colitis independent of changes to the microbial community caused by diarrhea.