Normalization of Single-Cell RNA-Seq Data

Normalization is an important step in the analysis of single-cell RNA-seq data. While no single method outperforms all others in all datasets, the choice of normalization can have profound impact on the results. Data-driven metrics can be used to rank normalization methods and select the best performers. Here, we show how to use R/Bioconductor to calculate normalization factors, apply them to compute normalized data, and compare several normalization approaches. Finally, we briefly show how to perform downstream analysis steps on the normalized data.

E.coli Nissle increases transcription of flagella assembly and formate hydrogenlyase genes in response to colitis

Escherichia coli Nissle (EcN), a probiotic bacterium, has been employed in treating inflammatory bowel disease, but the nature of its therapeutic effect is not fully understood. Intestinal inflammation alters the environment, exposing the microbial population to new stresses and eliciting transcriptional responses. We administered EcN to germ-free mice and then compared its transcriptional response between DSS-treated and untreated conditions using RNA-seq analysis to identify 187 differentially expressed genes (119 upregulated, 68 downregulated) and verifying a subset with qRT-PCR. The upregulated genes included many involved in flagella biosynthesis and motility, as well as several members of the formate hydrogenlyase complex. Despite prior evidence that these pathways are both transcriptionally regulated by nitric oxide, in vitro tests did not establish that nitric oxide exposure alone elicited the transcriptional response. The results provide new information on the transcriptional response of EcN to inflammation and establish a basis for further investigation of its anti-inflammatory activity.

Transcriptional Profiling During Neural Conversion

The processes that underlie neuronal conversion ultimately involve a reorganization of transcriptional networks to establish a neuronal cell fate. As such, transcriptional profiling is a key component toward understanding this process. In this chapter, we will discuss methods of elucidating transcriptional networks during neuronal reprogramming and considerations that should be incorporated in experimental design.

Functional Characterization of the EMBRYONIC FLOWER 2 Gene Involved in Flowering in Ginkgo biloba

Ginkgo biloba has edible, medicinal, and ornamental value. However, the long juvenile phase prevents the development of the G. biloba industry, and there are few reports on the identification and functional analysis of genes regulating the flowering time of G. biloba. EMBRYONIC FLOWER 2 (EMF), an important protein in flower development, functions to promote vegetative growth and repress flowering. In this study, a novel EMF gene (GbEMF2) was cloned and characterized from G. biloba. GbEMF2 contains a 2,193 bp open reading frame (ORF) encoding 730 amino acids. GbEMF2 harbors conserved VEFS-Box domain by the plant EMF protein. The phylogenic analysis showed that GbEMF2 originated from a polycomb-group (Pc-G) protein ancestor and was a member of the EMF2 protein. The quantitative real-time PCR (qRT-PCR) analysis revealed that GbEMF2 was expressed in all detected organs, and it showed a significantly higher level in ovulating strobilus and microstrobilus than in other organs. Compared with emf2 mutant plants, overexpression of GbEMF2 driven by the CaMV 35S promoter in emf2 mutant Arabidopsis plants delayed flowering but earlier than wild-type (WT) plants. This result indicated that GbEMF2 repressed flowering in G. biloba. Moreover, the RNA-seq analysis of GbEMF2 transgenic Arabidopsis plants (GbEMF2-OE/emf2), WT plants, and emf2 mutants screened out 227 differentially expressed genes (DEGs). Among these DEGs, FLC, MAF5, and MAF5-1 genes were related to flower organ development and regulated by GbEMF2. In addition, some genes participating in sugar metabolism, such as Alpha-amylase 1 (AMY1), BAM1, and Sucrose synthase 3 (SUS3) genes, were also controlled by GbEMF2. Overall, our results suggested that GbEMF2 negatively regulates flowering development in G. biloba. This finding provided a foundation and target gene for shortening the Ginkgo juvenile period by genetic engineering technology.

Characteristics of neural growth and cryopreservation of the dorsal root ganglion using three-dimensional collagen hydrogel culture versus conventional culture

In vertebrates, most somatosensory pathways begin with the activation of dorsal root ganglion (DRG) neurons. The development of an appropriate DRG culture method is a prerequisite for establishing in vitro peripheral nerve disease models and for screening therapeutic drugs. In this study, we compared the changes in morphology, molecular biology, and transcriptomics of chicken embryo DRG cultured on tissue culture plates (T-DRG) versus three-dimensional collagen hydrogels (C-DRG). Our results showed that after 7 days of culture, the transcriptomics of T-DRG and C-DRG were quite different. The upregulated genes in C-DRG were mainly related to neurogenesis, axon guidance, and synaptic plasticity, whereas the downregulated genes in C-DRG were mainly related to cell proliferation and cell division. In addition, the genes related to cycles/pathways such as the synaptic vesicle cycle, cyclic adenosine monophosphate signaling pathway, and calcium signaling pathway were activated, while those related to cell-cycle pathways were downregulated. Furthermore, neurogenesis- and myelination-related genes were highly expressed in C-DRG, while epithelial–mesenchymal transition-, apoptosis-, and cell division-related genes were suppressed. Morphological results indicated that the numbers of branches, junctions, and end-point voxels per C-DRG were significantly greater than those per T-DRG. Furthermore, cells were scattered in T-DRG and more concentrated in C-DRG, with a higher ratio of 5-ethynyl-2′-deoxyuridine (EdU)-positive cells in T-DRG compared with C-DRG. C-DRG also had higher S100 calcium-binding protein B (S100B) and lower α-smooth muscle actin (α-SMA) expression than T-DRG, and contained fewer terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)-positive cells after 48 hours of serum starvation. After cryopreservation, C-DRG maintained more intact morphological characteristics, and had higher viability and less TUNEL-positive cells than T-DRG. Furthermore, newly formed nerve bundles were able to grow along the existing Schwann cells in C-DRG. These results suggest that C-DRG may be a promising in vitro culture model, with better nerve growth and anti-apoptotic ability, quiescent Schwann cells, and higher viability. Results from this study provide a reference for the construction, storage, and transportation of tissue-engineered nerves. The study was approved by the Ethics Committee of Aier School of Ophthalmology, Central South University, China (approval No. 2020-IRB16), on March 15, 2020.

Transcriptional changes of mouse ovary during follicle initial or cyclic recruitment mediated by extra hormone treatment

AIMS

Folliculogenesis contains gonadotropin-independent and -dependent stage. Disruption in any of this process would induce failure in retrieving capable oocytes during clinical treatment. However, there is still limited understanding of the molecular components specifically regulating this process.

MATERIAL AND METHODS

Ovaries of P3, P20 and exogenous gonadotropin-treated P22 mice were sampled and underwent RNA-seq to investigate the transcriptome variance during mouse folliculogenesis.

KEY FINDINGS

In our dataset, 1883 and 626 DEGs were captured for each stage respectively, which were further clustered into eight expression patterns. Pathway enrichment analysis identified distinct biological processes enriched in two stages, with the most prominent being the pathways related to metabolism, gene expression, cell cycle, immune system and DNA methylation. Transcriptional regulator inference yielded eight master transcription factors (i.e. Runx1, Stat3, Sox3, Pou5f1, Gata4, Foxl2, Cebpb, and Esr1) driving folliculogenesis.

SIGNIFICANCE

Our study revealed the temporal transcriptional reprogramming and gene expression dynamics during folliculogenesis mediated by extra hormone treatment, which could provide novel insights to controlled ovarian stimulation in future infertility treatment.

Nucleoredoxin Gene TaNRX1 Positively Regulates Drought Tolerance in Transgenic Wheat (Triticum aestivum L.)

Drought is the main abiotic stress factor limiting the growth and yield of wheat (Triticum aestivum L.). Therefore, improving wheat tolerance to drought stress is essential for maintaining yield. Previous studies have reported on the important role of TaNRX1 in conferring drought stress tolerance. Therefore, to elucidate the regulation mechanism by which TaNRX1 confers drought resistance in wheat, we generated TaNRX1 overexpression (OE) and RNA interference (RNAi) wheat lines. The results showed that the tolerance of the OE lines to drought stress were significantly enhanced. The survival rate, leaf chlorophyll, proline, soluble sugar content, and activities of the antioxidant enzymes (catalase, superoxide dismutase, and peroxidase) of the OE lines were higher than those of the wild type (WT); however, the relative electrical conductivity and malondialdehyde, hydrogen peroxide, and superoxide anion levels of the OE lines were lower than those of the WT; the RNAi lines showed the opposite results. RNA-seq results showed that the common differentially expressed genes of TaNRX1 OE and RNAi lines, before and after drought stress, were mainly distributed in the plant-pathogen interaction, plant hormone signal transduction, phenylpropane biosynthesis, starch and sucrose metabolism, and carbon metabolism pathways and were related to the transcription factors, including WRKY, MYB, and bHLH families. This study suggests that TaNRX1 positively regulates drought stress tolerance in wheat.

Coupled Transcriptomics for Differential Expression Analysis and Determination of Transcription Start Sites: Design and Bioinformatics

The term coupled transcriptomics is coined to describe a design of an RNA-seq experiment intended for both differential expression analysis and genome-wide determination of the transcription start sites (TSS). The minimal requirements for the first analysis are two experimental conditions with at least two biological replicates enabling statistical tests. The second analysis involves the bioinformatics comparison of the data generated from a control RNA-seq library with another library enriched in primary transcripts using Terminator^™ 5'-phosphate-dependent exonuclease, in an experiment denominated differential RNA-seq (dRNA-seq). Usually, dRNA-seq is carried out with specific protocols for library construction, different of those used for common differential expression analysis. Our experimental design allows to use the same data for both analyses, reducing the number of libraries to be generated and sequenced. This is a guide for designing a coupled transcriptomics experiment and for the subsequent bioinformatics procedures. The proposed methods can be applied to the detection and study of small RNA genes.

ATP5A1 Participates in Transcriptional and Posttranscriptional Regulation of Cancer-Associated Genes by Modulating Their Expression and Alternative Splicing Profiles in HeLa Cells

Background: Aberrant expression and alternative splicing of oncogenes are the driving events in tumor initiation and development. But how these events are regulated in cancer cells is largely unknown. Functions of ATP5A1, an important mitochondrial ATP synthase gene, in transcriptional and posttranscriptional regulation were explored in this study. Methods: ATP5A1 was overexpressed using plasmid-transformed HeLa cells, and its influence on cell apoptosis and proliferation is evaluated. Transcriptome sequencing was then performed using RNA-seq to study the changes in gene expression and regulation of alternative splicing events. Validation of the implicated genes was achieved using RT-qPCR analysis. Results: It was found that ATP5A1 could significantly promote cellular apoptosis, but it had no influence on cell proliferation. ATP5A1 overexpression significantly increased the expression levels of genes associated with the innate immune response, angiogenesis, and collagen catabolic processes. This included enrichment of MMP2 and MMP19. It was also found that ATP5A1 could interfere with the alternative splicing of hundreds of genes associated with glucose homeostasis, HIF-1 signaling activation, and several pathways associated with cancers. Eight ATP5A1-regulated differentially expressed genes and 3 genes altered by splicing were selected and validated using RT-qPCR analysis. Conclusions: In summary, we illustrate the regulatory functions of ATP5A1 on the transcriptome of HeLa cells by exploring its influence on gene expression and alternative splicing. The results suggest that ATP5A1 may play an important regulatory role in cervical cancer cells by regulating expression and alternative splicing of cancer-associated genes. This study provides novel insights into the current understanding of the mechanisms of ATP5A1 on carcinogenesis and cancer progression.

Genome-Wide Identification and Expression Profiling of Germin-Like Proteins Reveal Their Role in Regulating Abiotic Stress Response in Potato

Germin and germin-like proteins (GLPs) perform a significant role in plants against biotic and abiotic stress. To understand the role of GLPs in potato, a comprehensive genome-wide analysis was performed in the potato genome. This study identified a total of 70 StGLPs genes in the potato genome, distributed among 11 chromosomes. Phylogenetic analysis exhibited that StGLPs were categorized into six groups with high bootstrap values. StGLPs gene structure and motifs analysis showed a relatively well-maintained intron-exon and motif formation within the cognate group. Additionally, several cis-elements in the promoter regions of GLPs were hormones, and stress-responsive and different families of miRNAs target StGLPs. Gene duplication under selection pressure also exhibited positive and purifying selections in StGLPs. In our results, the StGLP5 gene showed the highest expression in response to salt stress among all expressed StGLPs. Totally 19 StGLPs genes were expressed in response to heat stress. Moreover, three genes, StGLP30, StGLP17, and StGLP14, exhibited a relatively higher expression level in the potato after heat treatment. In total, 22 genes expressed in response to abscisic acid (ABA) treatment indicated that ABA performed an essential role in the plant defense or tolerance mechanism to environmental stress. RNA-Seq data validated by RT-qPCR also confirm that the StGLP5 gene showed maximum expression among selected genes under salt stress. Concisely, our results provide a platform for further functional exploration of the StGLPs against salt and heat stress conditions.

Neoantigen prediction in human breast cancer using RNA sequencing data

Neoantigens have attracted attention as biomarkers or therapeutic targets. However, accurate prediction of neoantigens is still challenging, especially in terms of its accuracy and cost. Variant detection using RNA sequencing (RNA-seq) data has been reported to be a low-accuracy but cost-effective tool, but the feasibility of RNA-seq data for neoantigen prediction has not been fully examined. In the present study, we used whole-exome sequencing (WES) and RNA-seq data of tumor and matched normal samples from six breast cancer patients to evaluate the utility of RNA-seq data instead of WES data in variant calling to detect neoantigen candidates. Somatic variants were called in three protocols using: (i) tumor and normal WES data (DNA method, Dm); (ii) tumor and normal RNA-seq data (RNA method, Rm); and (iii) combination of tumor RNA-seq and normal WES data (Combination method, Cm). We found that the Rm had both high false-positive and high false-negative rates because this method depended greatly on the expression status of normal transcripts. When we compared the results of Dm with those of Cm, only 14% of the neoantigen candidates detected in Dm were identified in Cm, but the majority of the missed candidates lacked coverage or variant allele reads in the tumor RNA. In contrast, about 70% of the neoepitope candidates with higher expression and rich mutant transcripts could be detected in Cm. Our results showed that Cm could be an efficient and a cost-effective approach to predict highly expressed neoantigens in tumor samples.

Usage of Metatranscriptomics to Understand Oral Disease

Metatranscriptomics is a method used to comprehensively capture bacterial activity within microbiota at the transcription level. It has become an alternative to the 16S rDNA sequencing, which uses only the 16S rRNA gene for predicting bacterial composition. By conducting metatranscriptomics, investigators can obtain substantial information about what types of genes are transcribed at the time of sampling and which bacterial taxa are responsible for their transcription. Here, I describe a protocol for metatranscriptomics for oral microbiota by using high-throughput sequencing technology. A remarkable feature of this protocol is that it uses the level of rRNA expression as the internal control for measuring transcriptional activity of each bacterial taxon. The normalized mRNA level is given by the mRNA/rRNA ratio, which indicates the extent of transcriptional activity.

Physiological and Transcripts Analyses Reveal the Mechanism by Which Melatonin Alleviates Heat Stress in Chrysanthemum Seedlings

Heat stress limits the growth and development of chrysanthemum seedlings. Although melatonin (MT) has been linked to the heat stress response in plants, research on the underlying molecular mechanisms is scarce. In this study, the regulatory networks of MT on heat stress in chrysanthemum seedlings were explored. Physiological measurements suggested that MT not only reduced malondialdehyde accumulation, hydrogen peroxide content, and superoxide anion free radical generation rate, but also significantly promoted osmotic regulation substance synthesis (proline and soluble protein), antioxidant accumulation (GSH and AsA), and the antioxidant enzyme activities (SOD, POD, CAT, and APX) in chrysanthemum leaves under heat stress. Furthermore, MT increased the fresh weight, dry weight, chlorophyll content, photosynthesis rate, and gas exchange indexes. Further, RNA-seq results revealed 33,497 and 36,740 differentially expressed genes in the S/Con and SMT/ConMT comparisons, respectively. The differences in the comparisons revealed that MT regulated heat shock transcription factors (HSFs) and heat shock proteins (HSPs), and the genes involved in Ca²⁺ signal transduction (CNGCs and CAM/CMLs), starch and sucrose metabolism (EDGL, BGLU, SuS, and SPS), hormone (PP2Cs, AUX/IAAs, EBFs, and MYC2), chlorophyll metabolism (HEMA and PORA), flavonoid biosynthesis (CHS, DFR, and FNS), and carotenoid biosynthesis (DXPS, GGDP, and PSY). MT effectively improved chrysanthemum seedling heat-resistance. Our study, thus, provides novel evidence of a gene network regulated by MT under heat stress.

Root Transcriptome and Metabolome Profiling Reveal Key Phytohormone-Related Genes and Pathways Involved Clubroot Resistance in Brassica rapa L

Plasmodiophora brassicae, an obligate biotrophic pathogen-causing clubroot disease, can seriously affect Brassica crops worldwide, especially Chinese cabbage. Understanding the transcriptome and metabolome profiling changes during the infection of P. brassicae will provide key insights in understanding the defense mechanism in Brassica crops. In this study, we estimated the phytohormones using targeted metabolome assays and transcriptomic changes using RNA sequencing (RNA-seq) in the roots of resistant (BrT24) and susceptible (Y510-9) plants at 0, 3, 9, and 20 days after inoculation (DAI) with P. brassicae. Differentially expressed genes (DEGs) in resistant vs. susceptible lines across different time points were identified. The weighted gene co-expression network analysis of the DEGs revealed six pathways including "Plant-pathogen interaction" and "Plant hormone signal transduction" and 15 hub genes including pathogenic type III effector avirulence factor gene (RIN4) and auxin-responsive protein (IAA16) to be involved in plants immune response. Inhibition of Indoleacetic acid, cytokinin, jasmonate acid, and salicylic acid contents and changes in related gene expression in R-line may play important roles in regulation of clubroot resistance (CR). Based on the combined metabolome profiling and hormone-related transcriptomic responses, we propose a general model of hormone-mediated defense mechanism. This study definitely enhances our current understanding and paves the way for improving CR in Brassica rapa.

Computational Identification of ceRNA and Reconstruction of ceRNA Regulatory Network Based on RNA-seq and Small RNA-seq Data in Plants

Competing endogenous RNAs (ceRNAs) are transcripts with the ability to competitively titrate microRNAs (miRNAs) against miRNA repressing target genes to post-transcriptionally regulate the expression of corresponding miRNAs. It is a newly discovered gene regulation pattern between longer RNA and miRNA molecules. Recent research has gradually revealed the functional significance of ceRNAs in regulating normal development and stress response processes in plants and animals, as well as in cancer genesis and metastasis. Therefore, ceRNA identification is an important and necessary step to deepen our understanding of the regulation mechanisms of various biological processes. Here, we provide a pipeline used to computationally identify plant ceRNAs and reconstruct ceRNA regulatory networks based on RNA-seq and small RNA-seq data.

Effects of the Thyroid Endocrine System on Gonadal Sex Ratios and Sex-Related Gene Expression in the Pufferfish Takifugu rubripes

To examine the effect and mechanism of thyroid hormone on gonadal sex differentiation, Takifugu rubripes larvae were treated with goitrogen (methimazole, MET, 1000 g/g), and thyroxine (T4, 2nM) from 25 to 80 days after hatching (dah). Gonadal histology and sex ratios of fish were then determined at 80 dah. MET treatment induced masculinization, but T4 treatment did not induce feminization in T. rubripes larvae. Transcriptomic analysis of gonads at 80 dah was then conducted. Among the large number of differentially expressed genes between the groups, the expression of foxl2, cyp19a1a, and dmrt1 was altered. The expression of foxl2, cyp19a1a, dmrt1 and gsdf at 25, 40, 55 days after treatment (dat) was further analyzed by qPCR. MET treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 in genetic females (p < 0.05). Additionally, T4 treatment induced an increase in the expression of cyp19a1a in genetic XY gonads only at 25 dat. However, the increase in cyp19a1a expression did not continue to 40 and 55 dat. This may explain why feminization of larvae was not found in the T4-treated group. Thus, the present study provides the first evidence that MET treatment causes masculinization in teleost fish. The effects of MET-induced masculinization in T. rubripes may act primarily via suppression of the expression of foxl2 and cyp19a1a, and stimulation of the expression of dmrt1. Moreover, the effects of higher concentrations of T4 or different concentrations of T3, on sex differentiation require further testing.

Identification of Genes Post-Transcriptionally Regulated from RNA-seq: The Case Study of Liver Hepatocellular Carcinoma

The field of transcriptional regulation generally assumes that changes in transcripts levels reflect changes in transcriptional status of the corresponding gene. While this assumption might hold true for a large population of transcripts, a considerable and still unrecognized fraction of the variation might involve other steps of the RNA lifecycle, that is the processing of the premature RNA, and degradation of the mature RNA. Discrimination between these layers requires complementary experimental techniques, such as RNA metabolic labeling or block of transcription experiments. Nonetheless, the analysis of the premature and mature RNA, derived from intronic and exonic read counts in RNA-seq data, allows distinguishing between transcriptionally and post-transcriptionally regulated genes, although not recognizing the specific step involved in the post-transcriptional response, that is processing, degradation, or a combination of the two. We illustrate how the INSPEcT R/Bioconductor package could be used to infer post-transcriptional regulation in TCGA RNA-seq samples for Hepatocellular Carcinoma.

Gene expression profile of THZ1-treated nasopharyngeal carcinoma cell lines indicates its involvement in the inhibition of the cell cycle

BACKGROUND

The aim of this study was to identify downstream target genes and pathways regulated by THZ1 in nasopharyngeal carcinoma (NPC).

METHODS

The gene expression profile of GSE95750 in two NPC cell lines, untreated group and treated with THZ1 group, was analyzed. Differentially expressed genes (DEGs) were compared using the R-software. Then Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathways (KEGG) was analyzed using Database for Annotation, Visualization, and Integrated Discovery (DAVID). Cytoscape was used for protein-protein interaction (PPI) analysis. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to verified the gene expression.

RESULTS

We identified 25 genes with increased expression and 567 genes with decreased expression in THZ1-treated NPC cells. The top 10 significantly DEGs between untreated group and THZ1 treated group were identified by qRT-PCR and the results were in agreement with RNA-seq. The total 592 DEGs were found enriched in 1,148 GO terms and 38 KEGG pathways. The most important enriched pathways identified were cell cycle related, and several related node genes were identified, such as CDC6, CDC34, CDK7, CDK9, CCNA2, CCNB1, CDT1, KIF11, LIN9, PLK1, and POLR family, which consistent with RNA-seq.

CONCLUSIONS

Our results emphasize the differential genes and pathways occurring in THZ1-treated NPC cells, which increases our understanding of the anti-tumor mechanisms of THZ1.

QuickIsoSeq for Isoform Quantification in Large-Scale RNA Sequencing

RNA-sequencing (RNA-seq) is a powerful technology for transcriptome profiling. While most RNA-seq projects focus on gene-level quantification and analysis, there is growing evidence that most mammalian genes are alternatively spliced to generate different isoforms that can be subsequently translated to protein molecules with diverse or even opposing biological functions. Quantifying the expression levels of these isoforms is key to understanding the genes biological functions in healthy tissues and the progression of diseases. Among open source tools developed for isoform quantification, Salmon, Kallisto, and RSEM are recommended based upon previous systematic evaluation of these tools using both experimental and simulated RNA-seq datasets. However, isoform quantification in practical RNA-seq data analysis needs to deal with many QC issues, such as the abundance of rRNAs in mRNA-seq, the efficiency of globin RNA depletion in whole blood samples, and potential sample swapping. To overcome these practical challenges, QuickIsoSeq was developed for large-scale RNA-seq isoform quantification along with QC. In this chapter, we describe the pipeline and detailed the steps required to deploy and use it to analyze RNA-seq datasets in practice. The QuickIsoSeq package can be downloaded from https://github.com/shanrongzhao/QuickIsoSeq.

Identification of Intronic Lariat-Derived Circular RNAs in Arabidopsis by RNA Deep Sequencing

Lariat RNAs are well-known by-products of pre-mRNA splicing in eukaryotes, which are produced by the excised introns when the 5' splice site (5' ss) joins with the branchpoint (BP) during splicing. In general, most of lariat RNAs are usually linearized by RNA debranching enzyme 1 (DBR1), followed by degradation for intron turnover. However, with the high-throughput RNA sequencing technology and bioinformatics methods, increasing evidences have shown that many lariat RNAs can stably accumulate under physiological conditions in both animals and plants. Here, we describe a large-scale analysis to systematically identify the lariat RNAs (i.e., intronic circular RNAs) in Arabidopsis by utilizing the RNA-sequencing data.

Architecture and Dynamics of the Wounding-Induced Gene Regulatory Network During the Oolong Tea Manufacturing Process (Camellia sinensis)

Understanding extensive transcriptional reprogramming events mediated by wounding during the oolong tea manufacturing process is essential for improving oolong tea quality. To improve our comprehension of the architecture of the wounding-induced gene regulatory network, we systematically analyzed the high-resolution transcriptomic and metabolomic data from wounding-treated (after turnover stage) tea leaves at 11 time points over a 220-min period. The results indicated that wounding activates a burst of transcriptional activity within 10 min and that the temporal expression patterns over time could be partitioned into 18 specific clusters with distinct biological processes. The transcription factor (TF) activity linked to the TF binding motif participated in specific biological processes within different clusters. A chronological model of the wounding-induced gene regulatory network provides insight into the dynamic transcriptional regulation event after wounding treatment (the turnover stage). Time series data of wounding-induced volatiles reveal the scientific significance of resting for a while after wounding treatment during the actual manufacturing process of oolong tea. Integrating information-rich expression data with information on volatiles allowed us to identify many high-confidence TFs participating in aroma formation regulation after wounding treatment by using weighted gene co-expression network analysis (WGCNA). Collectively, our research revealed the complexity of the wounding-induced gene regulatory network and described wounding-mediated dynamic transcriptional reprogramming events, serving as a valuable theoretical basis for the quality formation of oolong tea during the post-harvest manufacturing process.

Processing and Analysis of RNA-seq Data from Public Resources

Advances in next generation sequencing (NGS) technologies resulted in a broad array of large-scale gene expression studies and an unprecedented volume of whole messenger RNA (mRNA) sequencing data, or the transcriptome (also known as RNA sequencing, or RNA-seq). These include the Genotype Tissue Expression project (GTEx) and The Cancer Genome Atlas (TCGA), among others. Here we cover some of the commonly used datasets, provide an overview on how to begin the analysis pipeline, and how to explore and interpret the data provided by these publicly available resources.

Transcriptome analysis of the circadian clock gene BMAL1 deletion with opposite carcinogenic effects

We have previously reported that the deletion of BMAL1 gene has opposite effects in respect to its contribution to the pathways that are effective in the multistage carcinogenesis process. BMAL1 deletion sensitized nearly normal breast epithelial (MCF10A) and invasive breast cancer cells (MDA-MB-231) to cisplatin- and doxorubicin-induced apoptosis, while this deletion also aggravated the invasive potential of MDA-MB-231 cells. However, the mechanistic relationship of the seemingly opposite contribution of BMAL1 deletion to carcinogenesis process is not known at genome-wide level. In this study, an RNA-seq approach was taken to uncover the differentially expressed genes (DEGs) and pathways after treating BMAL1 knockout (KO) or wild-type (WT) MDA-MB-231 cells with cisplatin and doxorubicin to initiate apoptosis. Gene set enrichment analysis with the DEGs demonstrated that enrichment in multiple genes/pathways contributes to sensitization to cisplatin- or doxorubicin-induced apoptosis in BMAL1-dependent manner. Additionally, our DEG analysis suggested that non-coding transcript RNA (such as lncRNA and processed pseudogenes) may have role in cisplatin- or doxorubicin-induced apoptosis. Protein-protein interaction network obtained from common DEGs in cisplatin and doxorubicin treatments revealed that GSK3β, NACC1, and EGFR are the principal genes regulating the response of the KO cells. Moreover, the analysis of DEGs among untreated BMAL1 KO and WT cells revealed that epithelial-mesenchymal transition genes are up-regulated in KO cells. As a negative control, we have also analyzed the DEGs following treatment with an endoplasmic reticulum (ER) stress-inducing agent, tunicamycin, which was affected by BMAL1 deletion minimally. Collectively, the present study suggests that BMAL1 regulates many genes/pathways of which the alteration in BMAL1 KO cells may shed light on pleotropic phenotype observed.

Transcriptome Analysis of Resistance to Fusarium Wilt in Mung Bean (Vigna radiata L.)

Fusarium wilt is a destructive soil-borne disease that threatens the production of mung bean. Mung bean lines Zheng8-4 and Zheng8-20 show high resistance and high susceptibility to Fusarium wilt, respectively. Transcriptome analysis was carried out to identify candidate genes involved in Fusarium wilt resistance using Zheng8-4 and Zheng8-20 at 0, 0.5, 1, 2, and 4 days post inoculation (dpi). Differential expression analysis showed that 3,254 genes responded to pathogen infection and were differentially expressed in the resistant and susceptible lines. Weighted gene co-expression network analysis (WGCNA) was also performed to identify five modules highly correlated with Fusarium wilt resistance, in which 453 differentially expressed genes (DEGs) were considered likely to be involved in Fusarium wilt resistance. Among these DEGs, we found 24 genes encoding resistance (R) proteins, 22 encoding protein kinases, 20 belonging to transcription factor families, 34 encoding proteins with oxidoreductase activity, 17 involved in stimulation/stress responses, and 54 annotated to pathogen resistance-related pathways. Finally, 27 annotated genes were further selected as candidate genes of Fusarium wilt resistance in mung bean. This study identifies novel potential resistance-related genes against Fusarium wilt and provides a theoretical basis for further investigation of Fusarium wilt resistance in mung bean breeding.

RNA Editing in Human and Mouse Tissues

Adenosine-to-inosine (A-to-I) RNA editing is a fundamental posttranscriptional mechanism that greatly diversifies the transcriptome in many living organisms, including mammals. Multiple studies have demonstrated the importance of this process not just in normal development and physiology but also in various human diseases. Importantly, the precise editing level of a site may have downstream consequences on cellular behavior. Hence, the editing levels should be quantified as accurately as possible. In this chapter, we describe how to examine RNA editing in human and mouse tissues. The rapid development of next-generation sequencing technologies is affording us an unprecedented ability to accurately measure the editing levels of numerous sites simultaneously. Our experimental workflow includes the harvesting of high-quality RNA samples and the construction of different high-throughput sequencing libraries. We also delineate the computational steps needed to analyze the sequencing data from an Illumina platform.