Strand-Specific RNA-Seq Provides Greater Resolution of Transcriptome Profiling

CRISPR Screening Uncovers a Long-Range Enhancer for ONECUT1 in Pancreatic Differentiation and Links a Diabetes Risk Variant

Functional enhancer annotation is a valuable first step for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants for investigation. However, unbiased enhancer discovery in physiologically relevant contexts remains a major challenge. To discover regulatory elements pertinent to diabetes, we conducted a CRISPR interference screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers uncovered, we focused on a long-range enhancer ∼664 kb from the ONECUT1 promoter, since coding mutations in ONECUT1 cause pancreatic hypoplasia and neonatal diabetes. Homozygous enhancer deletion in hPSCs was associated with a near-complete loss of ONECUT1 gene expression and compromised pancreatic differentiation. This enhancer contains a confidently fine-mapped type 2 diabetes associated variant (rs528350911) which disrupts a GATA motif. Introduction of the risk variant into hPSCs revealed substantially reduced binding of key pancreatic transcription factors (GATA4, GATA6 and FOXA2) on the edited allele, accompanied by a slight reduction of ONECUT1 transcription, supporting a causal role for this risk variant in metabolic disease. This work expands our knowledge about transcriptional regulation in pancreatic development through the characterization of a long-range enhancer and highlights the utility of enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Effect of Curcumin on Hepatic mRNA and lncRNA Co-Expression in Heat-Stressed Laying Hens

Heat stress is an important factor affecting poultry production; birds have a range of inflammatory reactions under high-temperature environments. Curcumin has anti-inflammatory and antioxidant effects. The purpose of this experiment was to investigate the effect of dietary curcumin supplementation on the liver transcriptome of laying hens under heat stress conditions. In the animal experiment, a total of 240 Hy-Line brown hens aged 280 days were divided randomly into four different experimental diets with four replicates, and each replicate consisted of 15 hens during a 42-D experiment. The ambient temperature was adjusted to 34 ± 2 °C for 8 h per day, transiting to a range of 22 °C to 28 °C for the remaining 16 h. In the previous study of our lab, it was found that supplemental 150 mg/kg curcumin can improve production performance, antioxidant enzyme activity, and immune function in laying hens under heat stress. To further investigate the regulatory mechanism of curcumin on heat stress-related genes, in total, six samples of three liver tissues from each of 0 mg/kg and 150 mg/kg curcumin test groups were collected for RNA-seq analysis. In the transcriptome analysis, we reported for the first time that the genes related to heat stress of mRNA, such as HSPA8, HSPH1, HSPA2, and DNAJA4, were co-expressed with lncRNA such as XLOC010450, XLOC037987, XLOC053511, XLOC061207, and XLOC100318, and all of these genes are shown to be down-regulated. These findings provide a scientific basis for the possible benefits of dietary curcumin addition in heat-stressed laying hens.

Merging short and stranded long reads improves transcript assembly

Long-read RNA sequencing has arisen as a counterpart to short-read sequencing, with the potential to capture full-length isoforms, albeit at the cost of lower depth. Yet this potential is not fully realized due to inherent limitations of current long-read assembly methods and underdeveloped approaches to integrate short-read data. Here, we critically compare the existing methods and develop a new integrative approach to characterize a particularly challenging pool of low-abundance long noncoding RNA (lncRNA) transcripts from short- and long-read sequencing in two distinct cell lines. Our analysis reveals severe limitations in each of the sequencing platforms. For short-read assemblies, coverage declines at transcript termini resulting in ambiguous ends, and uneven low coverage results in segmentation of a single transcript into multiple transcripts. Conversely, long-read sequencing libraries lack depth and strand-of-origin information in cDNA-based methods, culminating in erroneous assembly and quantitation of transcripts. We also discover a cDNA synthesis artifact in long-read datasets that markedly impacts the identity and quantitation of assembled transcripts. Towards remediating these problems, we develop a computational pipeline to "strand" long-read cDNA libraries that rectifies inaccurate mapping and assembly of long-read transcripts. Leveraging the strengths of each platform and our computational stranding, we also present and benchmark a hybrid assembly approach that drastically increases the sensitivity and accuracy of full-length transcript assembly on the correct strand and improves detection of biological features of the transcriptome. When applied to a challenging set of under-annotated and cell-type variable lncRNA, our method resolves the segmentation problem of short-read sequencing and the depth problem of long-read sequencing, resulting in the assembly of coherent transcripts with precise 5' and 3' ends. Our workflow can be applied to existing datasets for superior demarcation of transcript ends and refined isoform structure, which can enable better differential gene expression analyses and molecular manipulations of transcripts.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex

BACKGROUND

Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques.

RESULTS

We found a discrete functional association, in which ^H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while ^H3K4me3HP loss played a vital role in synaptic activity. Moreover, ^H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and ^H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by ^H3K4me3HP and ^H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes.

CONCLUSIONS

Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.

LncPLAAT3-AS Regulates PLAAT3-Mediated Adipocyte Differentiation and Lipogenesis in Pigs through miR-503-5p

Animal fat deposition has a significant impact on meat flavor and texture. However, the molecular mechanisms of fat deposition are not well understood. LncPLAAT3-AS is a naturally occurring transcript that is abundant in porcine adipose tissue. Here, we focus on the regulatory role of lncPLAAT3-AS in promoting preadipocyte proliferation and adipocyte differentiation. By overexpressing or repressing lncPLAAT3 expression, we found that lncPLAAT3-AS promoted the transcription of its host gene PLAAT3, a regulator of adipocyte differentiation. In addition, we predicted the region of lncPLAAT3-AS that binds to miR-503-5p and showed by dual luciferase assay that lncPLAAT3-AS acts as a sponge to absorb miR-503-5p. Interestingly, miR-503-5p also targets and represses PLAAT3 expression and helps regulate porcine preadipocyte proliferation and differentiation. Taken together, these results show that lncPLAAT3-AS upregulates PLAAT3 expression by absorbing miR-503-5p, suggesting a potential regulatory mechanism based on competing endogenous RNAs. Finally, we explored lncPLAAT3-AS and PLAAT3 expression in adipose tissue and found that both molecules were expressed at significantly higher levels in fatty pig breeds compared to lean pig breeds. In summary, we identified the mechanism by which lncPLAAT3-AS regulates porcine preadipocyte proliferation and differentiation, contributing to our understanding of the molecular mechanisms of lipid deposition in pigs.

how_are_we_stranded_here: quick determination of RNA-Seq strandedness

Background

Quality control checks are the first step in RNA-Sequencing analysis, which enable the identification of common issues that occur in the sequenced reads. Checks for sequence quality, contamination, and complexity are commonplace, and allow users to implement steps downstream which can account for these issues. Strand-specificity of reads is frequently overlooked and is often unavailable even in published data, yet when unknown or incorrectly specified can have detrimental effects on the reproducibility and accuracy of downstream analyses.

Results

To address these issues, we developed how_are_we_stranded_here, a Python library that helps to quickly infer strandedness of paired-end RNA-Sequencing data. Testing on both simulated and real RNA-Sequencing reads showed that it correctly measures strandedness, and measures outside the normal range may indicate sample contamination.

Conclusions

how_are_we_stranded_here is fast and user friendly, making it easy to implement in quality control pipelines prior to analysing RNA-Sequencing data. how_are_we_stranded_here is freely available at https://github.com/betsig/how_are_we_stranded_here.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04572-7.

Rapid Transient Transcriptional Adaptation to Hypergravity in Jurkat T Cells Revealed by Comparative Analysis of Microarray and RNA-Seq Data

Cellular responses to micro- and hypergravity are rapid and complex and appear within the first few seconds of exposure. Transcriptomic analyses are a valuable tool to analyze these genome-wide cellular alterations. For a better understanding of the cellular dynamics upon altered gravity exposure, it is important to compare different time points. However, since most of the experiments are designed as endpoint measurements, the combination of cross-experiment meta-studies is inevitable. Microarray and RNA-Seq analyses are two of the main methods to study transcriptomics. In the field of altered gravity research, both methods are frequently used. However, the generation of these data sets is difficult and time-consuming and therefore the number of available data sets in this research field is limited. In this study, we investigated the comparability of microarray and RNA-Seq data and applied the results to a comparison of the transcriptomics dynamics between the hypergravity conditions during two real flight platforms and a centrifuge experiment to identify temporal adaptation processes. We performed a comparative study on an Affymetrix HTA2.0 microarray and a paired-end RNA-Seq data set originating from the same Jurkat T cell RNA samples from a short-term hypergravity experiment. The overall agreeability was high, with better sensitivity of the RNA-Seq analysis. The microarray data set showed weaknesses on the level of single upregulated genes, likely due to its normalization approach. On an aggregated level of biotypes, chromosomal distribution, and gene sets, both technologies performed equally well. The microarray showed better performance on the detection of altered gravity-related splicing events. We found that all initially altered transcripts fully adapted after 15 min to hypergravity and concluded that the altered gene expression response to hypergravity is transient and fully reversible. Based on the combined multiple-platform meta-analysis, we could demonstrate rapid transcriptional adaptation to hypergravity, the differential expression of the ATPase subunits ATP6V1A and ATP6V1D, and the cluster of differentiation (CD) molecules CD1E, CD2AP, CD46, CD47, CD53, CD69, CD96, CD164, and CD226 in hypergravity. We could experimentally demonstrate that it is possible to develop methodological evidence for the meta-analysis of individual data.

Non-coding Natural Antisense Transcripts: Analysis and Application

Non-coding natural antisense transcripts (ncNATs) are regulatory RNA sequences that are transcribed in the opposite direction to protein-coding or non-coding transcripts. These transcripts are implicated in a broad variety of biological and pathological processes, including tumorigenesis and oncogenic progression. With this complex field still in its infancy, annotations, expression profiling and functional characterisations of ncNATs are far less comprehensive than those for protein-coding genes, pointing out substantial gaps in the analysis and characterisation of these regulatory transcripts. In this review, we discuss ncNATs from an analysis perspective, in particular regarding the use of high-throughput sequencing strategies, such as RNA-sequencing, and summarize the unique challenges of investigating the antisense transcriptome. Finally, we elaborate on their potential as biomarkers and future targets for treatment, focusing on cancer.

Long Non-coding RNAs: Mechanisms, Experimental, and Computational Approaches in Identification, Characterization, and Their Biomarker Potential in Cancer

Long non-coding RNAs are diverse class of non-coding RNA molecules >200 base pairs of length having various functions like gene regulation, dosage compensation, epigenetic regulation. Dysregulation and genomic variations of several lncRNAs have been implicated in several diseases. Their tissue and developmental specific expression are contributing factors for them to be viable indicators of physiological states of the cells. Here we present an comprehensive review the molecular mechanisms and functions, state of the art experimental and computational pipelines and challenges involved in the identification and functional annotation of lncRNAs and their prospects as biomarkers. We also illustrate the application of co-expression networks on the TCGA-LIHC dataset for putative functional predictions of lncRNAs having a therapeutic potential in Hepatocellular carcinoma (HCC).

Identification and Analysis of Long Non-coding RNAs in Leuciscus waleckii Adapted to Highly Alkaline Conditions

Leuciscus waleckii is a freshwater fish that is known to inhabit the Dali Nor Lake, Inner Mongolia, China. The water in this lake has an HCO3 -/CO3 2- concentration of 54 mM (pH 9.6) and a salinity of 0.6‰. The physiological mechanisms that allow this fish to tolerate these saline/alkaline conditions have yet to be elucidated. Transcriptional component analysis has shown that the expression levels of a large number of genes involved in the pathways responsible for osmo-ionoregulation and arachidonic acid metabolism pathway expression change significantly (p < 0.05) during the regulation of acid-base balance under high alkaline stress. In this study, we investigated the role of long non-coding RNAs (lncRNAs) during adaptation to high alkaline conditions. Fish were challenged to an NaHCO3-adjusted alkalinity of 0 mM, 30 mM (pH 9.44 ± 0.08), and 50 mM (pH 9.55 ± 0.06) for 20 days in the laboratory. Gill and kidney tissues were then collected for high-throughput sequencing assays. A total of 159 million clean reads were obtained by high-throughput sequencing, and 41,248 lncRNA transcripts were identified. Of these, the mean number of exons and the mean length of the lncRNA transcripts were 4.8 and 2,079 bp, respectively. Based on the analysis of differential lncRNA transcript expression, a total of 5,244 and 6,571 lncRNA transcripts were found to be differentially expressed in the gills and kidneys, respectively. Results derived from Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the coding genes were correlated with the lncRNA expression profiles. GO analysis showed that many lncRNAs were enriched in the following processes: "transporter activity," "response to stimulus," and "binding." KEGG analysis further revealed that metabolic pathways were significantly enriched. A random selection of 16 lncRNA transcripts was tested by RT-qPCR; these results were consistent with our sequencing results. We found that a large number of genes, with the same expression profiles as those with differentially expressed lncRNAs, were associated with the regulation of acid-base balance, ion transport, and the excretion of ammonia and nitrogen. Collectively, our data indicate that lncRNA-regulated gene expression plays an important role in the process of adaptation to high alkaline conditions in L. waleckii.

A Scalable Strand-Specific Protocol Enabling Full-Length Total RNA Sequencing From Single Cells

RNA sequencing (RNAseq) has been widely used to generate bulk gene expression measurements collected from pools of cells. Only relatively recently have single-cell RNAseq (scRNAseq) methods provided opportunities for gene expression analyses at the single-cell level, allowing researchers to study heterogeneous mixtures of cells at unprecedented resolution. Tumors tend to be composed of heterogeneous cellular mixtures and are frequently the subjects of such analyses. Extensive method developments have led to several protocols for scRNAseq but, owing to the small amounts of RNA in single cells, technical constraints have required compromises. For example, the majority of scRNAseq methods are limited to sequencing only the 3' or 5' termini of transcripts. Other protocols that facilitate full-length transcript profiling tend to capture only polyadenylated mRNAs and are generally limited to processing only 96 cells at a time. Here, we address these limitations and present a novel protocol that allows for the high-throughput sequencing of full-length, total RNA at single-cell resolution. We demonstrate that our method produced strand-specific sequencing data for both polyadenylated and non-polyadenylated transcripts, enabled the profiling of transcript regions beyond only transcript termini, and yielded data rich enough to allow identification of cell types from heterogeneous biological samples.

Single-cell RNA sequencing in human lung cancer: Applications, challenges, and pathway towards personalized therapy

Lung cancer is one of the most prevalent human cancers, and single-cell RNA sequencing (scRNA-seq) has been widely used to study human lung cancer at the cellular, genetic, and molecular level. Even though there are published reviews, which summarized the applications of scRNA-seq in human cancers like breast cancer, there is lack of a comprehensive review, which could effectively highlight the broad use of scRNA-seq in studying lung cancer. This review, therefore, was aimed to summarize the various applications of scRNA-seq in human lung cancer research based on the findings from different published in vitro, in vivo, and clinical studies. The review would first briefly outline the concept and principle of scRNA-seq, followed by the discussion on the applications of scRNA-seq in studying human lung cancer. Finally, the challenges faced when using scRNA-seq to study human lung cancer would be discussed, and the potential applications and challenges of scRNA-seq to facilitate the development of personalized cancer therapy in the future would be explored.

Current advances of epigenetics in periodontology from ENCODE project: a review and future perspectives

BACKGROUND

The Encyclopedia of DNA Elements (ENCODE) project has advanced our knowledge of the functional elements in the genome and epigenome. The aim of this article was to provide the comprehension about current research trends from ENCODE project and establish the link between epigenetics and periodontal diseases based on epigenome studies and seek the future direction.

MAIN BODY

Global epigenome research projects have emphasized the importance of epigenetic research for understanding human health and disease, and current international consortia show an improved interest in the importance of oral health with systemic health. The epigenetic studies in dental field have been mainly conducted in periodontology and have focused on DNA methylation analysis. Advances in sequencing technology have broadened the target for epigenetic studies from specific genes to genome-wide analyses.

CONCLUSIONS

In line with global research trends, further extended and advanced epigenetic studies would provide crucial information for the realization of comprehensive dental medicine and expand the scope of ongoing large-scale research projects.

Hopeful monsters: unintended sequencing of famously malformed mite mitochondrial tRNAs reveals widespread expression and processing of sense-antisense pairs

Although tRNA structure is one of the most conserved and recognizable shapes in molecular biology, aberrant tRNAs are frequently found in the mitochondrial genomes of metazoans. The extremely degenerate structures of several mitochondrial tRNAs (mt-tRNAs) have led to doubts about their expression and function. Mites from the arachnid superorder Acariformes are predicted to have some of the shortest mt-tRNAs, with a complete loss of cloverleaf-like shape. While performing mitochondrial isolations and recently developed tRNA-seq methods in plant tissue, we inadvertently sequenced the mt-tRNAs from a common plant pest, the acariform mite Tetranychus urticae, to a high enough coverage to detect all previously annotated T. urticae tRNA regions. The results not only confirm expression, CCA-tailing and post-transcriptional base modification of these highly divergent tRNAs, but also revealed paired sense and antisense expression of multiple T. urticae mt-tRNAs. Mirrored expression of mt-tRNA genes has been hypothesized but not previously demonstrated to be common in any system. We discuss the functional roles that these divergent tRNAs could have as both decoding molecules in translation and processing signals in transcript maturation pathways, as well as how sense–antisense pairs add another dimension to the bizarre tRNA biology of mitochondrial genomes.

Tight association of genome rearrangements with gene expression in conifer plastomes

BACKGROUND

Our understanding of plastid transcriptomes is limited to a few model plants whose plastid genomes (plastomes) have a highly conserved gene order. Consequently, little is known about how gene expression changes in response to genomic rearrangements in plastids. This is particularly important in the highly rearranged conifer plastomes.

RESULTS

We sequenced and reported the plastomes and plastid transcriptomes of six conifer species, representing all six extant families. Strand-specific RNAseq data show a nearly full transcription of both plastomic strands and detect C-to-U RNA-editing sites at both sense and antisense transcripts. We demonstrate that the expression of plastid coding genes is strongly functionally dependent among conifer species. However, the strength of this association declines as the number of plastomic rearrangements increases. This finding indicates that plastomic rearrangement influences gene expression.

CONCLUSIONS

Our data provide the first line of evidence that plastomic rearrangements not only complicate the plastomic architecture but also drive the dynamics of plastid transcriptomes in conifers.

Endogenous membrane stress induces T6SS activity in Pseudomonas aeruginosa

In this study, we applied a CRISPR interference approach to knock down expression of essential genes involved in cell envelope biogenesis and determined if this can trigger Pseudomonas aeruginosa T6SS apparatus assembly. We found that disruption of envelope biogenesis can be sensed by the bacteria via a sensory pathway that involves phosphorylation leading to profound T6SS dynamic activation. Our data provide further evidence that membrane disruption is a key signal that drives T6SS organelle assembly. This process likely mimics natural activation pathways that allow P. aeruginosa to detect the attack by other aggressive bacteria species and thus sheds further light on the potential key signals that control T6SS organelle biogenesis.

The type 6 secretion system (T6SS) is a dynamic organelle encoded by many gram-negative bacteria that can be used to kill competing bacterial prey species in densely occupied niches. Some predatory species, such as Vibrio cholerae, use their T6SS in an untargeted fashion while in contrast, Pseudomonas aeruginosa assembles and fires its T6SS apparatus only after detecting initial attacks by other bacterial prey cells; this targeted attack strategy has been termed the T6SS tit-for-tat response. Molecules that interact with the P. aeruginosa outer membrane such as polymyxin B can also trigger assembly of T6SS organelles via a signal transduction pathway that involves protein phosphorylation. Recent work suggests that a phospholipase T6SS effector (TseL) of V. cholerae can induce T6SS dynamic activity in P. aeruginosa when delivered to or expressed in the periplasmic space of this organism. Here, we report that inhibiting expression of essential genes involved in outer membrane biogenesis can also trigger T6SS activation in P. aeruginosa. Specifically, we developed a CRISPR interference (CRISPRi) system to knock down expression of bamA, tolB, and lptD and found that these knockdowns activated T6SS activity. This increase in T6SS activity was dependent on the same signal transduction pathway that was previously shown to be required for the tit-for-tat response. We conclude that outer membrane perturbation can be sensed by P. aeruginosa to activate the T6SS even when the disruption is generated by aberrant cell envelope biogenesis.

Short-Term Magnesium Deficiency Triggers Nutrient Retranslocation in Arabidopsis thaliana

Magnesium (Mg) is essential for many biological processes in plant cells, and its deficiency causes yield reduction in crop systems. Low Mg status reportedly affects photosynthesis, sucrose partitioning and biomass allocation. However, earlier physiological responses to Mg deficiency are scarcely described. Here, we report that Mg deficiency in Arabidopsis thaliana first modified the mineral profile in mature leaves within 1 or 2 days, then affected sucrose partitioning after 4 days, and net photosynthesis and biomass production after 6 days. The short-term Mg deficiency reduced the contents of phosphorus (P), potassium, manganese, zinc and molybdenum in mature but not in expanding (young) leaves. While P content decreased in mature leaves, P transport from roots to mature leaves was not affected, indicating that Mg deficiency triggered retranslocation of the mineral nutrients from mature leaves. A global transcriptome analysis revealed that Mg deficiency triggered the expression of genes involved in defence response in young leaves.

Omics for Bioprospecting and Drug Discovery from Bacteria and Microalgae

"Omics" represent a combinatorial approach to high-throughput analysis of biological entities for various purposes. It broadly encompasses genomics, transcriptomics, proteomics, lipidomics, and metabolomics. Bacteria and microalgae exhibit a wide range of genetic, biochemical and concomitantly, physiological variations owing to their exposure to biotic and abiotic dynamics in their ecosystem conditions. Consequently, optimal conditions for adequate growth and production of useful bacterial or microalgal metabolites are critically unpredictable. Traditional methods employ microbe isolation and 'blind'-culture optimization with numerous chemical analyses making the bioprospecting process laborious, strenuous, and costly. Advances in the next generation sequencing (NGS) technologies have offered a platform for the pan-genomic analysis of microbes from community and strain downstream to the gene level. Changing conditions in nature or laboratory accompany epigenetic modulation, variation in gene expression, and subsequent biochemical profiles defining an organism's inherent metabolic repertoire. Proteome and metabolome analysis could further our understanding of the molecular and biochemical attributes of the microbes under research. This review provides an overview of recent studies that have employed omics as a robust, broad-spectrum approach for screening bacteria and microalgae to exploit their potential as sources of drug leads by focusing on their genomes, secondary metabolite biosynthetic pathway genes, transcriptomes, and metabolomes. We also highlight how recent studies have combined molecular biology with analytical chemistry methods, which further underscore the need for advances in bioinformatics and chemoinformatics as vital instruments in the discovery of novel bacterial and microalgal strains as well as new drug leads.

MoAIMS: efficient software for detection of enriched regions of MeRIP-Seq

BACKGROUND

Methylated RNA immunoprecipitation sequencing (MeRIP-Seq) is a popular sequencing method for studying RNA modifications and, in particular, for N6-methyladenosine (m6A), the most abundant RNA methylation modification found in various species. The detection of enriched regions is a main challenge of MeRIP-Seq analysis, however current tools either require a long time or do not fully utilize features of RNA sequencing such as strand information which could cause ambiguous calling. On the other hand, with more attention on the treatment experiments of MeRIP-Seq, biologists need intuitive evaluation on the treatment effect from comparison. Therefore, efficient and user-friendly software that can solve these tasks must be developed.

RESULTS

We developed a software named "model-based analysis and inference of MeRIP-Seq (MoAIMS)" to detect enriched regions of MeRIP-Seq and infer signal proportion based on a mixture negative-binomial model. MoAIMS is designed for transcriptome immunoprecipitation sequencing experiments; therefore, it is compatible with different RNA sequencing protocols. MoAIMS offers excellent processing speed and competitive performance when compared with other tools. When MoAIMS is applied to studies of m6A, the detected enriched regions contain known biological features of m6A. Furthermore, signal proportion inferred from MoAIMS for m6A treatment datasets (perturbation of m6A methyltransferases) showed a decreasing trend that is consistent with experimental observations, suggesting that the signal proportion can be used as an intuitive indicator of treatment effect.

CONCLUSIONS

MoAIMS is efficient and easy-to-use software implemented in R. MoAIMS can not only detect enriched regions of MeRIP-Seq efficiently but also provide intuitive evaluation on treatment effect for MeRIP-Seq treatment datasets.

Identifying inaccuracies in gene expression estimates from unstranded RNA-seq data

RNA-seq methods are widely utilized for transcriptomic profiling of biological samples. However, there are known caveats of this technology which can skew the gene expression estimates. Specifically, if the library preparation protocol does not retain RNA strand information then some genes can be erroneously quantitated. Although strand-specific protocols have been established, a significant portion of RNA-seq data is generated in non-strand-specific manner. We used a comprehensive stranded RNA-seq dataset of 15 blood cell types to identify genes for which expression would be erroneously estimated if strand information was not available. We found that about 10% of all genes and 2.5% of protein coding genes have a two-fold or higher difference in estimated expression when strand information of the reads was ignored. We used parameters of read alignments of these genes to construct a machine learning model that can identify which genes in an unstranded dataset might have incorrect expression estimates and which ones do not. We also show that differential expression analysis of genes with biased expression estimates in unstranded read data can be recovered by limiting the reads considered to those which span exonic boundaries. The resulting approach is implemented as a package available at https://github.com/mikpom/uslcount.

Isolating Male Meiocytes from Maize and Wheat for "-Omics" Analyses

Genome-wide gene expression studies have become a routine approach due to the advances in sequencing technologies, their ease of use, and increasing affordability. Simultaneous investigation of small RNA expression adds further valuable information but is not adopted as widely yet. Both RNA-seq and small RNA-seq benefit from the use of specific cell types. Here, we describe a protocol for the isolation of male meiotic cells from maize or wheat plants, along with the application of downstream RNA sequencing, extendable to other -omics approaches.

Comparative evaluation of RNA-Seq library preparation methods for strand-specificity and low input

Library preparation is a key step in sequencing. For RNA sequencing there are advantages to both strand specificity and working with minute starting material, yet until recently there was no kit available enabling both. The Illumina TruSeq stranded mRNA Sample Preparation kit (TruSeq) requires abundant starting material while the Takara Bio SMART-Seq v4 Ultra Low Input RNA kit (V4) sacrifices strand specificity. The SMARTer Stranded Total RNA-Seq Kit v2 - Pico Input Mammalian (Pico) by Takara Bio claims to overcome these limitations. Comparative evaluation of these kits is important for selecting the appropriate protocol. We compared the three kits in a realistic differential expression analysis. We prepared and sequenced samples from two experimental conditions of biological interest with each of the three kits. We report differences between the kits at the level of differential gene expression; for example, the Pico kit results in 55% fewer differentially expressed genes than TruSeq. Nevertheless, the agreement of the observed enriched pathways suggests that comparable functional results can be obtained. In summary we conclude that the Pico kit sufficiently reproduces the results of the other kits at the level of pathway analysis while providing a combination of options that is not available in the other kits.

Elucidating the editome: bioinformatics approaches for RNA editing detection

RNA editing is a widespread co/posttranscriptional mechanism affecting primary RNAs by specific nucleotide modifications, which plays relevant roles in molecular processes including regulation of gene expression and/or the processing of noncoding RNAs. In recent years, the detection of editing sites has been improved through the availability of high-throughput RNA sequencing (RNA-Seq) technologies. Accurate bioinformatics pipelines are essential for the analysis of next-generation sequencing (NGS) data to ensure the correct identification of edited sites. Several pipelines, using various read mappers and variant callers with a wide range of adjustable parameters, are available for the detection of RNA editing events. In this review, we discuss some of the most recent and popular tools and provide guidelines for RNA-Seq data generation and analysis for the detection of RNA editing in massive transcriptome data. Using simulated and real data sets, we provide an overview of their behavior, emphasizing the fact that the RNA editing detection in NGS data sets remains a challenging task.

Production of 14α-hydroxysteroids by a recombinant Saccharomyces cerevisiae biocatalyst expressing of a fungal steroid 14α-hydroxylation system

The 14α-hydroxysteroids have specific anti-gonadotropic and carcinolytic biological activities and can be produced by microbial biotransformation. The steroid 11β-/14α-hydroxylase P-450_lun from Cochliobolus lunatus is the only fungal cytochrome P450 enzyme identified to date with steroid C14 hydroxylation ability. Previous work has mainly revealed the 11β-hydroxylation activity of the P-450_lun towards cortexolone (RSS) substrate; however, the potential steroid 14α-hydroxylation activity of this enzyme, especially for androstenedione (AD) substrate, has not yet conducted in-depth testing. In this work, we further tested the steroid 14α-hydroxylation activity of the P-450_lun towards RSS and AD in the Saccharomyces cerevisiae system. We demonstrated that P-450_lun functions as the specific 14α-hydroxylase towards the AD substrate (regiospecificity > 99%); however, it showed a poor C14-hydroxylation regiospecificity (around 40%) for the RSS substrate. In addition, through transcriptome analysis combined with gene functional characterizations, we also identified and cloned the gene for the P-450_lun-associated redox partner CPR_lun. Finally, through codon optimization, knockout of genes for the side reactions related enzymes GCY1 and YPR1, and increasing copies of the P-450_lun and CPR_lun, we developed a recombinant S. cerevisiae biocatalyst based on the C. lunatus steroid 14α-hydroxylation system to produce 14α-hydroxysteroids. Initial production of 14α-OH-AD (150 mg/L day productivity, 99% regioisomeric purity, and 60% w/w yield) and 14α-OH-RSS (64 mg/L day productivity, 40% regioisomeric purity, and 26% w/w yield) were separately achieved in shake flasks; these results represent the highest level of 14α-hydroxysteroid production in the current yeast system.

Comparison of library construction kits for mRNA sequencing in the Illumina platform

BACKGROUND

The emergence of next-generation sequencing (NGS) technologies has made a tremendous contribution to the deciphering and significance of transcriptome analysis in biological fields. Since the advent of NGS technology in 2007, Illumina, Inc. has provided one of the most widely used sequencing platforms for NGS analysis.

OBJECTIVE

Although reagents and protocols provided by Illumina are adequately performed in transcriptome sequencing, recently, alternative reagents and protocols which are relatively cost effective are accessible. However, the kits derived from various manufacturers have advantages and disadvantages when researchers carry out the transcriptome library construction.

METHODS

We compared them using a variety of protocols to produce Illumina-compatible libraries based on transcriptome. Three different mRNA sequencing kits were selected for this study: TruSeq^® RNA Sample Preparation V2 (Illumina, Inc., USA), Universal Plus mRNA-Seq (NuGEN, Ltd., UK), and NEBNext^® Ultra™ Directional RNA Library Prep Kit for Illumina^® (New England BioLabs, Ltd., USA). We compared them focusing on cost, experimental time, and data output.

RESULTS

The quality and quantity of sequencing data obtained through the NGS technique were strongly influenced by the type of the sequencing library kits. It suggests that for transcriptome studies, researchers should select a suitable library construction kit according to the goal and resources of experiments.

CONCLUSION

The present work will help researchers to choose the right sequencing library construction kit for transcriptome analyses.

Long non-coding and coding RNA profiling using strand-specific RNA-seq in human hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) represents one of the most common heritable heart diseases. However, the signalling pathways and regulatory networks underlying the pathogenesis of HCM remain largely unknown. Here, we present a strand-specific RNA-seq dataset for both coding and lncRNA profiling in myocardial tissues from 28 HCM patients and 9 healthy donors. This dataset constitutes a valuable resource for the community to examine the dysregulated coding and lncRNA genes in HCM versus normal conditions.

Mapping the dsRNA World

Long double-stranded RNAs (dsRNAs) are abundantly expressed in animals, in which they frequently occur in introns and 3' untranslated regions of mRNAs. Functions of long, cellular dsRNAs are poorly understood, although deficiencies in adenosine deaminases that act on RNA, or ADARs, promote their recognition as viral dsRNA and an aberrant immune response. Diverse dsRNA-binding proteins bind cellular dsRNAs, hinting at additional roles. Understanding these roles is facilitated by mapping the genomic locations that express dsRNA in various tissues and organisms. ADAR editing provides a signature of dsRNA structure in cellular transcripts. In this review, we detail approaches to map ADAR editing sites and dsRNAs genome-wide, with particular focus on high-throughput sequencing methods and considerations for their successful application to the detection of editing sites and dsRNAs.

RNA-Seq profiling in peripheral blood mononuclear cells of amyotrophic lateral sclerosis patients and controls

Coding and long non-coding RNA (lncRNA) metabolism is now revealing its crucial role in Amyotrophic Lateral Sclerosis (ALS) pathogenesis. In this work, we present a dataset obtained via Illumina RNA-seq analysis on Peripheral Blood Mononuclear Cells (PBMCs) from sporadic and mutated ALS patients (mutations in FUS, TARDBP, SOD1 and VCP genes) and healthy controls. This dataset allows the whole-transcriptome characterization of PBMCs content, both in terms of coding and non-coding RNAs, in order to compare the disease state to the healthy controls, both for sporadic patients and for mutated patients. Our dataset is a starting point for the omni-comprehensive analysis of coding and lncRNAs, from an easy to withdraw, manage and store tissue that shows to be a suitable model for RNA profiling in ALS.

Hidden Gems in the Transcriptome Maps of Competent Streptococci

Natural transformation is regarded as an important mechanism in bacteria that allows for adaptation to different environmental stressors by ensuring genome plasticity. Since the discovery of this phenomenon in Streptococcus pneumoniae, remarkable progress has been made in the understanding of the molecular mechanisms and pathways coordinating this process. Recently, the advent of high-throughput sequencing allows the posing of questions that address the system at a larger scale but also allow for the creation of high-resolution maps of transcription. Thus, while much is already known about genetic competence in streptococci, recent studies continue to reveal intricate novel regulation pathways and components. In this perspective article, we highlight the use of transcriptional profiling and mapping as a valuable resource in the identification and characterization of “hidden gems” pertinent to the natural transformation system. Such strategies have recently been employed in a variety of different species. In S. mutans, for example, genome editing combined with the power of promoter mapping and RNA-Seq allowed for the identification of a link between the ComCDE and the ComRS systems, a ComR positive feedback loop mediated by SigX, and the XrpA peptide, encoded within sigX, which inhibits competence. In S. pneumoniae, a novel member of the competence regulon termed BriC was found to be directly under control of ComE and to promote biofilm formation and nasopharyngeal colonization but not competence. Together these new technologies enable us to discover new links and to revisit old pathways in the compelling study of natural genetic transformation.

Human Olfactory Receptors: Novel Cellular Functions Outside of the Nose

Olfactory receptors (ORs) are not exclusively expressed in the olfactory sensory neurons; they are also observed outside of the olfactory system in all other human tissues tested to date, including the testis, lung, intestine, skin, heart, and blood. Within these tissues, certain ORs have been determined to be exclusively expressed in only one tissue, whereas other ORs are more widely distributed in many different tissues throughout the human body. For most of the ectopically expressed ORs, limited data are available for their functional roles. They have been shown to be involved in the modulation of cell-cell recognition, migration, proliferation, the apoptotic cycle, exocytosis, and pathfinding processes. Additionally, there is a growing body of evidence that they have the potential to serve as diagnostic and therapeutic tools, as ORs are highly expressed in different cancer tissues. Interestingly, in addition to the canonical signaling pathways activated by ORs in olfactory sensory neurons, alternative pathways have been demonstrated in nonolfactory tissues. In this review, the existing data concerning the expression, as well as the physiological and pathophysiological functions, of ORs outside of the nose are highlighted to provide insights into future lines of research.

Crosstalk in competing endogenous RNA networks reveals new circular RNAs involved in the pathogenesis of early HIV infection

Background

The events in early HIV infection (EHI) are important determinants of disease severity and progression rate to AIDS, but the mechanisms of pathogenesis in EHI have not been fully understood. Circular RNAs (circRNAs) have been verified as “microRNA sponges” that regulate gene expression through competing endogenous RNA (ceRNA) networks, but circRNA expression profiles and their contribution to EHI pathogenesis are still unclear.

Methods

Two different libraries were constructed with RNA from human peripheral blood mononuclear cells from 3 HARRT-naive EHI patients and 3 healthy controls (HCs). The complete transcriptomes were sequenced with RNA sequencing (RNA-Seq) and miRNA sequencing (miRNA-Seq). The differentially expressed (DE) RNAs were validated with RT-qPCR. The circRNA profile and circRNA-associated-ceRNA network in EHI were analyzed with the integrated data of RNA-Seq and miRNA-Seq. Gene ontology (GO) analysis was used to annotate the circRNAs involved in the circRNA-associated-ceRNA networks.

Results

A total of 1365 circRNAs, 30 miRNAs, and 2049 mRNAs were differentially expressed between HARRT-naive EHI patients and HCs. A ceRNA network was constructed with 516 DE circRNAs and 903 DE mRNAs that shared miR response elements with 21 DE miRNAs. GO analysis demonstrated the multiple roles of the circRNAs enriched in EHI with circRNA-associated-ceRNA networks, such as immune response, inflammatory response and defense responses to virus, 67 circRNAs were revealed to be potentially involved in HIV-1 replication through regulating the expression of CCNK, CDKN1A and IL-15.

Conclusions

This study, for the first time, revealed a large circRNA profile and complex pathogenesis roles of circRNAs in EHI. A group of enriched circRNAs and associated circRNA-associated-ceRNA networks might contribute to HIV replication regulation and provide novel potential targets for both the pathogenesis of EHI and antiviral therapy.

Electronic supplementary material

The online version of this article (10.1186/s12967-018-1706-1) contains supplementary material, which is available to authorized users.

Comparative evaluation of cDNA library construction approaches for RNA-Seq analysis from low RNA-content human specimens

With the emergence of RNA sequencing technologies, metatranscriptomic studies are rapidly gaining attention as they simultaneously provide insight into gene expression profiles and therefore disease association pathways of microbial pathogens and their hosts. This approach, therefore, holds promise for applicability in infectious disease diagnostics. A challenge of this approach in the clinical setting is the low amount and quality of RNA, especially microbial RNA in most clinically-infected specimens. Here, we compared two commercially available stranded cDNA library preparation kits, the NuGEN Ovation SoLo RNA-Seq System and the Illumina TruSeq Stranded Total RNA, using RNA extracted from synovial and sonicate fluids from a subject with periprosthetic joint infection. The Ovation SoLo RNA-Seq System provided more useful transcriptomic data for the infecting bacterium, whereas the TruSeq Stranded Total RNA kit provided more useful human transcriptomic data.

Widespread antisense transcription of Populus genome under drought

Antisense transcription is widespread in many genomes and plays important regulatory roles in gene expression. The objective of our study was to investigate the extent and functional relevance of antisense transcription in forest trees. We employed Populus, a model tree species, to probe the antisense transcriptional response of tree genome under drought, through stranded RNA-seq analysis. We detected nearly 48% of annotated Populus gene loci with antisense transcripts and 44% of them with co-transcription from both DNA strands. Global distribution of reads pattern across annotated gene regions uncovered that antisense transcription was enriched in untranslated regions while sense reads were predominantly mapped in coding exons. We further detected 1185 drought-responsive sense and antisense gene loci and identified a strong positive correlation between the expression of antisense and sense transcripts. Additionally, we assessed the antisense expression in introns and found a strong correlation between intronic expression and exonic expression, confirming antisense transcription of introns contributes to transcriptional activity of Populus genome under drought. Finally, we functionally characterized drought-responsive sense-antisense transcript pairs through gene ontology analysis and discovered that functional groups including transcription factors and histones were concordantly regulated at both sense and antisense transcriptional level. Overall, our study demonstrated the extensive occurrence of antisense transcripts of Populus genes under drought and provided insights into genome structure, regulation pattern and functional significance of drought-responsive antisense genes in forest trees. Datasets generated in this study serve as a foundation for future genetic analysis to improve our understanding of gene regulation by antisense transcription.

Strand-specific RNA-seq analysis of the Acidithiobacillus ferrooxidans transcriptome in response to magnesium stress

Bioleaching is a promising process for 350 million tons Jinchuan low-grade pentlandite. But, Jinchuan pentlandite has lots of magnesium and high concentration of Mg²⁺ is harmful to bioleaching microorganisms. Thus, finding a way to improve the adaption of microorganisms to Mg²⁺ is a key for bioleaching. In the study, we found that oxidizing activity, bioleaching ability and biofilm formation of A.f were inhibited by Mg²⁺ stress. In addition, we analyzed mRNA and small RNA (sRNA) of Acidithiobacillus ferrooxidans (A.f) under Mg²⁺ stress by strand-specific RNA-sequencing (ssRNA-seq). After the bioinformatics process, 2475 coding genes were obtained, and there were 33 differential expression genes (DEGs) in 0.1 M-VS-Con, including 28 down-regulated and 5 up-regulated, whereas 52 DEGs were obtained in 0.5 M-VS-Con, including 28 down-regulated and 24 up-regulated. Gene ontology analysis showed most of DEGs were involved in catalytic activity, metabolic process and single-organism process. Furthermore, we identified 636 sRNA and some differential expression sRNA that may respond to Mg²⁺ stress. Further analysis of DEGs suggested that Mg²⁺ stress reduced biofilm formation perhaps through inhibiting Type IV Pili-related gene expression and inhibited bacterial activity perhaps through affecting carbon fixation. The study provided the foundation to understand the mechanisms of Mg²⁺ resistance in A.f and may be helpful to improve bioleaching ability for pentlandit.

Preparing cDNA Libraries from Lytic Phage-Infected Cells for Whole Transcriptome Analysis by RNA-Seq

Whole genome wide analysis of transcription using RNA-Seq methods is a powerful way to elucidate differential expression of gene features in bacteria across different conditions as well as for discovering previously exotic RNA species. Indeed, RNA sequencing has revolutionized the study of bacterial transcription with the diversity and quantity of small noncoding RNA elements that have been found and its ability to clearly define operons, promoters , and terminators . We discuss our experience with applying RNA sequencing technology to analyzing the lytic cycle, including extraction, processing, and a guide to the customized statistical analysis necessary for analyzing differential host and phage transcription.

Increasing quality, throughput and speed of sample preparation for strand-specific messenger RNA sequencing

BACKGROUND

RNA-Sequencing (RNA-seq) is now commonly used to reveal quantitative spatiotemporal snapshots of the transcriptome, the structures of transcripts (splice variants and fusions) and landscapes of expressed mutations. However, standard approaches for library construction typically require relatively high amounts of input RNA, are labor intensive, and are time consuming. METHODS: Here, we report the outcome of a systematic effort to optimize and streamline steps in strand-specific RNA-seq library construction. RESULTS: This work has resulted in the identification of an optimized messenger RNA isolation protocol, a potent reverse transcriptase for cDNA synthesis, and an efficient chemistry and a simplified formulation of library construction reagents. We also present an optimization of bead-based purification and size selection designed to maximize the recovery of cDNA fragments.

CONCLUSIONS

These developments have allowed us to assemble a rapid high throughput pipeline that produces high quality data from amounts of total RNA as low as 25 ng. While the focus of this study is on RNA-seq sample preparation, some of these developments are also relevant to other next-generation sequencing library types.

The emerging role of circular RNAs in transcriptome regulation

Circular RNAs (circRNAs) are a recently discovered form of RNA that has been found to regulate mammalian transcription. CircRNAs are covalently closed, single-stranded transcripts produced from precursor mRNA. While initially circRNAs were considered to be splicing artefacts, next-generation RNA sequencing of non-polyadenylated transcriptomes has recently shown that the expression of circRNAs is widespread and over 20% of expressed genes in examined cells and tissues can produce these transcripts. Until now thousands of circRNAs have been discovered in organisms ranging from Drosophila melanogaster to Homo sapiens. Functional studies indicate that these transcripts regulate expression of protein-coding linear transcripts and thus comprise an important component of gene expression regulation. Here we provide a comprehensive overview on the biology of circRNAs, including the expression patterns and function. Moreover, we discuss current methodologies for the discovery and validation of circular transcripts. Finally, perspectives on the utilization of circRNA as molecular markers of complex diseases are presented.

Global transcriptome analysis reveals extensive gene remodeling, alternative splicing and differential transcription profiles in non-seed vascular plant Selaginella moellendorffii

Background

Selaginella moellendorffii, a lycophyte, is a model plant to study the early evolution and development of vascular plants. As the first and only sequenced lycophyte to date, the genome of S. moellendorffii revealed many conserved genes and pathways, as well as specialized genes different from flowering plants. Despite the progress made, little is known about long noncoding RNAs (lncRNA) and the alternative splicing (AS) of coding genes in S. moellendorffii. Its coding gene models have not been fully validated with transcriptome data. Furthermore, it remains important to understand whether the regulatory mechanisms similar to flowering plants are used, and how they operate in a non-seed primitive vascular plant.

Results

RNA-sequencing (RNA-seq) was performed for three S. moellendorffii tissues, root, stem, and leaf, by constructing strand-specific RNA-seq libraries from RNA purified using RiboMinus isolation protocol. A total of 176 million reads (44 Gbp) were obtained from three tissue types, and were mapped to S. moellendorffii genome. By comparing with 22,285 existing gene models of S. moellendorffii, we identified 7930 high-confidence novel coding genes (a 35.6% increase), and for the first time reported 4422 lncRNAs in a lycophyte. Further, we refined 2461 (11.0%) of existing gene models, and identified 11,030 AS events (for 5957 coding genes) revealed for the first time for lycophytes. Tissue-specific gene expression with functional implication was analyzed, and 1031, 554, and 269 coding genes, and 174, 39, and 17 lncRNAs were identified in root, stem, and leaf tissues, respectively. The expression of critical genes for vascular development stages, i.e. formation of provascular cells, xylem specification and differentiation, and phloem specification and differentiation, was compared in S. moellendorffii tissues, indicating a less complex regulatory mechanism in lycophytes than in flowering plants. The results were further strengthened by the evolutionary trend of seven transcription factor families related to vascular development, which was observed among four representative species of seed and non-seed vascular plants, and nonvascular land and aquatic plants.

Conclusions

The deep RNA-seq study of S. moellendorffii discovered extensive new gene contents, including novel coding genes, lncRNAs, AS events, and refined gene models. Compared to flowering vascular plants, S. moellendorffii displayed a less complexity in both gene structure, alternative splicing, and regulatory elements of vascular development. The study offered important insight into the evolution of vascular plants, and the regulation mechanism of vascular development in a non-seed plant.

Electronic supplementary material

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

InFusion: Advancing Discovery of Fusion Genes and Chimeric Transcripts from Deep RNA-Sequencing Data

Analysis of fusion transcripts has become increasingly important due to their link with cancer development. Since high-throughput sequencing approaches survey fusion events exhaustively, several computational methods for the detection of gene fusions from RNA-seq data have been developed. This kind of analysis, however, is complicated by native trans-splicing events, the splicing-induced complexity of the transcriptome and biases and artefacts introduced in experiments and data analysis. There are a number of tools available for the detection of fusions from RNA-seq data; however, certain differences in specificity and sensitivity between commonly used approaches have been found. The ability to detect gene fusions of different types, including isoform fusions and fusions involving non-coding regions, has not been thoroughly studied yet. Here, we propose a novel computational toolkit called InFusion for fusion gene detection from RNA-seq data. InFusion introduces several unique features, such as discovery of fusions involving intergenic regions, and detection of anti-sense transcription in chimeric RNAs based on strand-specificity. Our approach demonstrates superior detection accuracy on simulated data and several public RNA-seq datasets. This improved performance was also evident when evaluating data from RNA deep-sequencing of two well-established prostate cancer cell lines. InFusion identified 26 novel fusion events that were validated in vitro, including alternatively spliced gene fusion isoforms and chimeric transcripts that include intergenic regions. The toolkit is freely available to download from http:/bitbucket.org/kokonech/infusion.

Transcriptional regulation of long-term potentiation

Long-term potentiation (LTP), the persistent strengthening of synapses following high levels of stimulation, is a form of synaptic plasticity that has been studied extensively as a possible mechanism for learning and memory formation. The strengthening of the synapse that occurs during LTP requires cascades of complex molecular processes and the coordinated remodeling of pre-synaptic and post-synaptic neurons. Despite over four decades of research, our understanding of the transcriptional mechanisms and molecular processes underlying LTP remains incomplete. Identification of all the proteins and non-coding RNA transcripts expressed during LTP may provide greater insight into the molecular mechanisms involved in learning and memory formation.

Characterization of the Olfactory Receptors Expressed in Human Spermatozoa

The detection of external cues is fundamental for human spermatozoa to locate the oocyte in the female reproductive tract. This task requires a specific chemoreceptor repertoire that is expressed on the surface of human spermatozoa, which is not fully identified to date. Olfactory receptors (ORs) are candidate molecules and have been attributed to be involved in sperm chemotaxis and chemokinesis, indicating an important role in mammalian spermatozoa. An increasing importance has been suggested for spermatozoal RNA, which led us to investigate the expression of all 387 OR genes. This study provides the first comprehensive analysis of OR transcripts in human spermatozoa of several individuals by RNA-Seq. We detected 91 different transcripts in the spermatozoa samples that could be aligned to annotated OR genes. Using stranded mRNA-Seq, we detected a class of these putative OR transcripts in an antisense orientation, indicating a different function, rather than coding for a functional OR protein. Nevertheless, we were able to detect OR proteins in various compartments of human spermatozoa, indicating distinct functions in human sperm. A panel of various OR ligands induced Ca(2+) signals in human spermatozoa, which could be inhibited by mibefradil. This study indicates that a variety of ORs are expressed at the mRNA and protein level in human spermatozoa.

RNA-Seq as a Tool to Study the Tumor Microenvironment

The transcriptome is composed of different types of RNA molecules including mRNAs, tRNAs, rRNAs, and other noncoding RNAs that are found inside a cell at a given time. Analyzing transcriptome patterns can shed light on the functional state of the cell as well as on the dynamics of cellular behavior associated with genomic and environmental changes. Likewise, transcriptome analysis has been a major help in solving biological issues and understanding the molecular basis of many diseases including human cancers. Specifically, since targeted and whole genome sequencing studies are becoming more common in identifying the driving factors of cancer, a comprehensive and high-resolution analysis of the transcriptome, as provided by RNA-Sequencing (RNA-Seq), plays a key role in investigating the functional relevance of the identified genomic aberrations. Here, we describe experimental procedures of RNA-Seq and downstream data processing and analysis, with a focus on the identification of abnormally expressed transcripts and genes.

Alzheimer's as a Systems-Level Disease Involving the Interplay of Multiple Cellular Networks

Alzheimer's disease (AD), and many neurodegenerative disorders, are multifactorial in nature. They involve a combination of genomic, epigenomic, interactomic and environmental factors. Progress is being made, and these complex diseases are beginning to be understood as having their origin in altered states of biological networks at the cellular level. In the case of AD, genomic susceptibility and mechanisms leading to (or accompanying) the impairment of the central Amyloid Precursor Protein (APP) processing and tau networks are widely accepted as major contributors to the diseased state. The derangement of these networks may result in both the gain and loss of functions, increased generation of toxic species (e.g., toxic soluble oligomers and aggregates) and imbalances, whose effects can propagate to supra-cellular levels. Although well sustained by empirical data and widely accepted, this global perspective often overlooks the essential roles played by the main counteracting homeostatic networks (e.g., protein quality control/proteostasis, unfolded protein response, protein folding chaperone networks, disaggregases, ER-associated degradation/ubiquitin proteasome system, endolysosomal network, autophagy, and other stress-protective and clearance networks), whose relevance to AD is just beginning to be fully realized. In this chapter, an integrative perspective is presented. Alzheimer's disease is characterized to be a result of: (a) intrinsic genomic/epigenomic susceptibility and, (b) a continued dynamic interplay between the deranged networks and the central homeostatic networks of nerve cells. This interplay of networks will underlie both the onset and rate of progression of the disease in each individual. Integrative Systems Biology approaches are required to effect its elucidation. Comprehensive Systems Biology experiments at different 'omics levels in simple model organisms, engineered to recapitulate the basic features of AD may illuminate the onset and sequence of events underlying AD. Indeed, studies of models of AD in simple organisms, differentiated cells in culture and rodents are beginning to offer hope that the onset and progression of AD, if detected at an early stage, may be stopped, delayed, or even reversed, by activating or modulating networks involved in proteostasis and the clearance of toxic species. In practice, the incorporation of next-generation neuroimaging, high-throughput and computational approaches are opening the way towards early diagnosis well before irreversible cell death. Thus, the presence or co-occurrence of: (a) accumulation of toxic Aβ oligomers and tau species; (b) altered splicing and transcriptome patterns; (c) impaired redox, proteostatic, and metabolic networks together with, (d) compromised homeostatic capacities may constitute relevant 'AD hallmarks at the cellular level' towards reliable and early diagnosis. From here, preventive lifestyle changes and tailored therapies may be investigated, such as combined strategies aimed at both lowering the production of toxic species and potentiating homeostatic responses, in order to prevent or delay the onset, and arrest, alleviate, or even reverse the progression of the disease.