An optimized protocol for generation and analysis of Ion Proton sequencing reads for RNA-Seq

OsMADS6-OsMADS32 and REP1 control palea cellular heterogeneity and morphogenesis in rice

Precise regulation of cell proliferation and differentiation is vital for organ morphology. Rice palea, serving as sepal, comprises two distinct regions: the marginal region (MRP) and body of palea (BOP), housing heterogeneous cell populations, which makes it an ideal system for studying organ morphogenesis. We report that the transcription factor (TF) REP1 promotes epidermal cell proliferation and differentiation in the BOP, resulting in hard silicified protrusion cells, by regulating the cyclin-dependent kinase gene, OsCDKB1;1. Conversely, TFs OsMADS6 and OsMADS32 are expressed exclusively in the MRP, where they limit cell division rates by inhibiting OsCDKB2;1 expression and promote endoreduplication, yielding elongated epidermal cells. Furthermore, reciprocal inhibition between the OsMADS6-OsMADS32 complex and REP1 fine-tunes the balance between cell division and differentiation during palea morphogenesis. We further show the functional conservation of these organ identity genes in heterogeneous cell growth in Arabidopsis, emphasizing a critical framework for controlling cellular heterogeneity in organ morphogenesis.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

Transcriptomics-Based Approach Identifies Spinosad-Associated Targets in the Colorado Potato Beetle, Leptinotarsa decemlineata

Simple Summary

The Colorado potato beetle Leptinotarsa decemlineata is a potato pest that can cause substantial damages to potato crops worldwide. Multiple approaches have been leveraged to control this pest including the use of a variety of insecticides. Resistance to different insecticides aimed at controlling this insect has been reported and much work has been conducted in recent years to elucidate the underlying molecular changes associated with insecticide resistance in L. decemlineata. However, information is sparse regarding the molecular impact associated with spinosad treatment in this insect pest. The current study thus explores transcriptional changes associated with spinosad response in L. decemlineata exposed to this compound using high-throughput sequencing. Results presented show multiple transcripts of interest that exhibit differential expression in spinosad-treated L. decemlineata and provide a preliminary footprint of transcripts affected by this insecticide in this potato pest. Select targets identified in this signature should be further explored in follow-up studies to better characterize their contribution, if any, in the process of spinosad resistance.

Abstract

The Colorado potato beetle Leptinotarsa decemlineata is an insect pest that threatens potato crops globally. The primary method to control its damage on potato plants is the use of insecticides, including imidacloprid, chlorantraniliprole and spinosad. However, insecticide resistance has been frequently observed in Colorado potato beetles. The molecular targets and the basis of resistance to imidacloprid and chlorantraniliprole have both been previously quantified. This work was undertaken with the overarching goal of better characterizing the molecular changes associated with spinosad exposure in this insect pest. Next-generation sequencing was conducted to identify transcripts that were differentially expressed between Colorado potato beetles exposed to spinosad versus control insects. Results showed several transcripts that exhibit different expression levels between the two conditions, including ones coding for venom carboxylesterase-6, chitinase 10, juvenile hormone esterase and multidrug resistance-associated protein 4. In addition, several microRNAs, such as miR-12-3p and miR-750-3p, were also modulated in the investigated conditions. Overall, this work reveals a molecular footprint underlying spinosad response in Colorado potato beetles and provides novel leads that could be targeted as part of RNAi-based approaches to control this insect pest.

Exosomal regulation of lymphocyte homing to the gut

Exosomes secreted from T cells have been shown to affect dendritic cells, cancer cells, and other T cells. However, little is known about how T-cell exosomes (T exosomes) modulate endothelial cell functions in the context of tissue-specific homing. Here, we study the roles of T exosomes in the regulation of gut-specific T-cell homing. The gut-tropic T cells induced by retinoic acid secrete the exosomes that upregulate integrin α4β7 binding to the MAdCAM-1 expressed on high endothelial venules in the gut. T exosomes were preferentially distributed to the villi of the small intestine in an α4β7-dependent manner. Exosomes from gut-tropic T cells suppressed the expression of MAdCAM-1 in the small intestine, thereby inhibiting T-cell homing to the gut. Moreover, microRNA (miRNA) profiling analysis has shown that exosomes from gut-tropic T cells were enriched with miRNAs targeting NKX2.3, a transcription factor critical to MAdCAM-1 expression. Taken together, our study proposes that α4β7-expressing T exosomes distribute themselves to the small intestine and modify the expression of microenvironmental tissues such that any subsequent lymphocyte homing is precluded. This may represent a novel mechanism by which excessive lymphocyte homing to the intestinal tissues is downsized.

Analysis of pooled genome sequences from Djallonke and Sahelian sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment

BACKGROUND

The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually.

RESULTS

Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis.

CONCLUSIONS

Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.

Transcriptome profile of Corynebacterium pseudotuberculosis in response to iron limitation

BACKGROUND

Iron is an essential micronutrient for the growth and development of virtually all living organisms, playing a pivotal role in the proliferative capability of many bacterial pathogens. The impact that the bioavailability of iron has on the transcriptional response of bacterial species in the CMNR group has been widely reported for some members of the group, but it hasn't yet been as deeply explored in Corynebacterium pseudotuberculosis. Here we describe for the first time a comprehensive RNA-seq whole transcriptome analysis of the T1 wild-type and the Cp13 mutant strains of C. pseudotuberculosis under iron restriction. The Cp13 mutant strain was generated by transposition mutagenesis of the ciuA gene, which encodes a surface siderophore-binding protein involved in the acquisition of iron. Iron-regulated acquisition systems are crucial for the pathogenesis of bacteria and are relevant targets to the design of new effective therapeutic approaches.

RESULTS

Transcriptome analyses showed differential expression in 77 genes within the wild-type parental T1 strain and 59 genes in Cp13 mutant under iron restriction. Twenty-five of these genes had similar expression patterns in both strains, including up-regulated genes homologous to the hemin uptake hmu locus and two distinct operons encoding proteins structurally like hemin and Hb-binding surface proteins of C. diphtheriae, which were remarkably expressed at higher levels in the Cp13 mutant than in the T1 wild-type strain. These hemin transport protein genes were found to be located within genomic islands associated with known virulent factors. Down-regulated genes encoding iron and heme-containing components of the respiratory chain (including ctaCEF and qcrCAB genes) and up-regulated known iron/DtxR-regulated transcription factors, namely ripA and hrrA, were also identified differentially expressed in both strains under iron restriction.

CONCLUSION

Based on our results, it can be deduced that the transcriptional response of C. pseudotuberculosis under iron restriction involves the control of intracellular utilization of iron and the up-regulation of hemin acquisition systems. These findings provide a comprehensive analysis of the transcriptional response of C. pseudotuberculosis, adding important understanding of the gene regulatory adaptation of this pathogen and revealing target genes that can aid the development of effective therapeutic strategies against this important pathogen.

A Novel lncRNA Regulates the Toll-Like Receptor Signaling Pathway and Related Immune Function by Stabilizing FOS mRNA as a Competitive Endogenous RNA

Long non-coding RNAs (lncRNAs) have recently emerged as new regulatory molecules with diverse functions in regulating gene expression and significant roles in the immune response. However, the function of many unknown lncRNAs is still unclear. By studying the regulatory effect of daidzein (DA) on immunity, we identified a novel lncRNA with an immune regulatory function: lncRNA- XLOC_098131. In vivo, DA treatment upregulated the expression of lncRNA- XLOC_098131, FOS, and JUN in chickens and affected the expression of activator protein 1 (AP-1) to regulate MAPK signaling, Toll-like receptor signaling, and related mRNA expression. It also enhanced macrophage activity and increased the numbers of blood neutrophils and mononuclear cells, which can improve the body's ability to respond to stress and bacterial and viral infections. Furthermore, DA treatment also reduced B lymphocyte apoptosis and promoted the differentiation of B lymphocytes into plasma cells, which in turn resulted in the production of more immunoglobulins and the promotion of antigen presentation. In vitro, using HEK293FT cells, we demonstrated that mir-548s could bind to and decrease the expression of both FOS and lncRNA- XLOC_098131. LncRNA- XLOC_098131 served as a competitive endogenous RNA to stabilize FOS by competitively binding to miR-548s and thereby reducing its inhibitory effect of FOS expression. Therefore, we concluded that the novel lncRNA XLOC_098131 acts as a key regulatory molecule that can regulate the Toll-like receptor signaling pathway and related immune function by serving as a competitive endogenous RNA to stabilize FOS mRNA expression.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

Dietary genistein supplementation for breeders and their offspring improves the growth performance and immune function of broilers

Genistein (GEN) is mainly extracted from soy plants and has potential functions as an antioxidant and in promoting immune function and growth. This study evaluated the effects of feeding breeders and their offspring dietary GEN on the immune function and growth performance of broiler chicks. Breeders were assigned to a control diet or GEN diet (control diet +400 mg/kg GEN), and their offspring were fed a control diet or GEN diet (control diet +40 mg/kg GEN). GEN treatment increased the body weight gain, tibial length, tibial width and slaughter performance of broilers and decreased the feed conversion ratio. The treatment also affected skeletal muscle myosin assembly and growth and increased growth hormone levels and IGF-I and IGFBP1 expression. Following GEN treatment, antigen processing and presentation, macrophage activation, B lymphocyte, NK cell and helper T cell proliferation, and CD4+ T lymphocyte differentiation all increased significantly. Increases were also observed in IgM and IgG concentrations, antibody titers, and antioxidant capacity. In addition, GEN treatment activated the Toll-like receptor signaling pathway and MAPK cascade signaling pathway. In summary, dietary GEN supplementation for breeders and their offspring can improve the growth performance and immune function of broiler chicks.

Transcriptomics-Related Mechanisms of Supplementing Laying Broiler Breeder Hens with Dietary Daidzein to Improve the Immune Function and Growth Performance of Offspring

Daidzein (DA) is an isoflavone that is primarily extracted from soy plants. This study evaluated the effects of supplementing laying broiler breeder hens with dietary DA on the immune function and growth performance of their offspring and the underlying mechanism. A total of 720 breeders were divided into three treatment groups that were fed either a control diet (CON), a DA-low-supplemented diet (DLS, CON+20 mg/kg DA), or a DA-high-supplemented diet (DHS, CON+100 mg/kg DA) for 8 weeks, and eggs were collected for hatching during the final week. The broiler offspring received a basal diet for 42 days, and blood, livers, and immune organs were collected at 21 and 42 days of age. DLS treatment promoted embryonic development and increased growth hormone levels, body weight, feed intake, and carcass traits on days 21 and 42 of broilers. Additionally, the IgA and IgG concentrations, antibody titers, and antioxidant capacity of broilers were increased at 21 days of age, and B lymphocyte differentiation was increased at 42 days. Besides, DLS treatment upregulated the expression of genes related to embryonic and muscle development in offspring and regulated mitogen-activated protein kinase (MAPK), transforming growth factor beta (TGF-β), nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), and Toll-like receptor signaling. DHS treatment decreased the percentage of abdominal fat in the broilers at 42 days, but it did not significantly affect embryonic development, growth performance, or IgA and IgG concentrations. In summary, providing dietary DA supplementation at 20 mg/kg to broiler breeders can improve their immune function and growth performance.

Target discovery screens using pooled shRNA libraries and next-generation sequencing: A model workflow and analytical algorithm

In the search for novel therapeutic targets, RNA interference screening has become a valuable tool. High-throughput technologies are now broadly accessible but their assay development from baseline remains resource-intensive and challenging. Focusing on this assay development process, we here describe a target discovery screen using pooled shRNA libraries and next-generation sequencing (NGS) deconvolution in a cell line model of Ewing sarcoma. In a strategy designed for comparative and synthetic lethal studies, we screened for targets specific to the A673 Ewing sarcoma cell line. Methods, results and pitfalls are described for the entire multi-step screening procedure, from lentiviral shRNA delivery to bioinformatics analysis, illustrating a complete model workflow. We demonstrate that successful studies are feasible from the first assay performance and independent of specialized screening units. Furthermore, we show that a resource-saving screen depth of 100-fold average shRNA representation can suffice to generate reproducible target hits despite heterogeneity in the derived datasets. Because statistical analysis methods are debatable for such datasets, we created ProFED, an analysis package designed to facilitate descriptive data analysis and hit calling using an aim-oriented profile filtering approach. In its versatile design, this open-source online tool provides fast and easy analysis of shRNA and other count-based datasets to complement other analytical algorithms.

A comparison of Illumina and Ion Torrent sequencing platforms in the context of differential gene expression

Background

Though Illumina has largely dominated the RNA-Seq field, the simultaneous availability of Ion Torrent has left scientists wondering which platform is most effective for differential gene expression (DGE) analysis. Previous investigations of this question have typically used reference samples derived from cell lines and brain tissue, and do not involve biological variability. While these comparisons might inform studies of tissue-specific expression, marked by large-scale transcriptional differences, this is not the common use case.

Results

Here we employ a standard treatment/control experimental design, which enables us to evaluate these platforms in the context of the expression differences common in differential gene expression experiments. Specifically, we assessed the hepatic inflammatory response of mice by assaying liver RNA from control and IL-1β treated animals with both the Illumina HiSeq and the Ion Torrent Proton sequencing platforms. We found the greatest difference between the platforms at the level of read alignment, a moderate level of concordance at the level of DGE analysis, and nearly identical results at the level of differentially affected pathways. Interestingly, we also observed a strong interaction between sequencing platform and choice of aligner. By aligning both real and simulated Illumina and Ion Torrent data with the twelve most commonly-cited aligners in the literature, we observed that different aligner and platform combinations were better suited to probing different genomic features; for example, disentangling the source of expression in gene-pseudogene pairs.

Conclusions

Taken together, our results indicate that while Illumina and Ion Torrent have similar capacities to detect changes in biology from a treatment/control experiment, these platforms may be tailored to interrogate different transcriptional phenomena through careful selection of alignment software.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4011-0) contains supplementary material, which is available to authorized users.

RED-ML: a novel, effective RNA editing detection method based on machine learning

With the advancement of second generation sequencing techniques, our ability to detect and quantify RNA editing on a global scale has been vastly improved. As a result, RNA editing is now being studied under a growing number of biological conditions so that its biochemical mechanisms and functional roles can be further understood. However, a major barrier that prevents RNA editing from being a routine RNA-seq analysis, similar to gene expression and splicing analysis, for example, is the lack of user-friendly and effective computational tools. Based on years of experience of analyzing RNA editing using diverse RNA-seq datasets, we have developed a software tool, RED-ML: RNA Editing Detection based on Machine learning (pronounced as "red ML"). The input to RED-ML can be as simple as a single BAM file, while it can also take advantage of matched genomic variant information when available. The output not only contains detected RNA editing sites, but also a confidence score to facilitate downstream filtering. We have carefully designed validation experiments and performed extensive comparison and analysis to show the efficiency and effectiveness of RED-ML under different conditions, and it can accurately detect novel RNA editing sites without relying on curated RNA editing databases. We have also made this tool freely available via GitHub . We have developed a highly accurate, speedy and general-purpose tool for RNA editing detection using RNA-seq data. With the availability of RED-ML, it is now possible to conveniently make RNA editing a routine analysis of RNA-seq. We believe this can greatly benefit the RNA editing research community and has profound impact to accelerate our understanding of this intriguing posttranscriptional modification process.

Optimized approach for Ion Proton RNA sequencing reveals details of RNA splicing and editing features of the transcriptome

RNA-sequencing (RNA-seq) has become the standard method for unbiased analysis of gene expression but also provides access to more complex transcriptome features, including alternative RNA splicing, RNA editing, and even detection of fusion transcripts formed through chromosomal translocations. However, differences in library methods can adversely affect the ability to recover these different types of transcriptome data. For example, some methods have bias for one end of transcripts or rely on low-efficiency steps that limit the complexity of the resulting library, making detection of rare transcripts less likely. We tested several commonly used methods of RNA-seq library preparation and found vast differences in the detection of advanced transcriptome features, such as alternatively spliced isoforms and RNA editing sites. By comparing several different protocols available for the Ion Proton sequencer and by utilizing detailed bioinformatics analysis tools, we were able to develop an optimized random primer based RNA-seq technique that is reliable at uncovering rare transcript isoforms and RNA editing features, as well as fusion reads from oncogenic chromosome rearrangements. The combination of optimized libraries and rapid Ion Proton sequencing provides a powerful platform for the transcriptome analysis of research and clinical samples.