How the Ethylene Biosynthesis Pathway of Semi-Halophytes Is Modified with Prolonged Salinity Stress Occurrence?

The mechanism of ethylene (ET)-regulated salinity stress response remains largely unexplained, especially for semi-halophytes and halophytes. Here, we present the results of the multifaceted analysis of the model semi-halophyte Mesembryanthemum crystallinum L. (common ice plant) ET biosynthesis pathway key components' response to prolonged (14 days) salinity stress. Transcriptomic analysis revealed that the expression of 3280 ice plant genes was altered during 14-day long salinity (0.4 M NaCl) stress. A thorough analysis of differentially expressed genes (DEGs) showed that the expression of genes involved in ET biosynthesis and perception (ET receptors), the abscisic acid (ABA) catabolic process, and photosynthetic apparatus was significantly modified with prolonged stressor presence. To some point this result was supported with the expression analysis of the transcript amount (qPCR) of key ET biosynthesis pathway genes, namely ACS6 (1-aminocyclopropane-1-carboxylate synthase) and ACO1 (1-aminocyclopropane-1-carboxylate oxidase) orthologs. However, the pronounced circadian rhythm observed in the expression of both genes in unaffected (control) plants was distorted and an evident downregulation of both orthologs' was induced with prolonged salinity stress. The UPLC-MS analysis of the ET biosynthesis pathway rate-limiting semi-product, namely of 1-aminocyclopropane-1-carboxylic acid (ACC) content, confirmed the results assessed with molecular tools. The circadian rhythm of the ACC production of NaCl-treated semi-halophytes remained largely unaffected by the prolonged salinity stress episode. We speculate that the obtained results represent an image of the steady state established over the past 14 days, while during the first hours of the salinity stress response, the view could be completely different.

A Comparative Transcriptome Analysis Unveils the Mechanisms of Response in Feather Degradation by Pseudomonas aeruginosa Gxun-7

Microbial degradation of feathers offers potential for bioremediation, yet the microbial response mechanisms warrant additional investigation. In prior work, Pseudomonas aeruginosa Gxun-7, which demonstrated robust degradation of feathers at elevated concentrations, was isolated. However, the molecular mechanism of this degradation remains only partially understood. To investigate this, we used RNA sequencing (RNA-seq) to examine the genes that were expressed differentially in P. aeruginosa Gxun-7 when exposed to 25 g/L of feather substrate. The RNA-seq analysis identified 5571 differentially expressed genes; of these, 795 were upregulated and 603 were downregulated. Upregulated genes primarily participated in proteolysis, amino acid, and pyruvate metabolism. Genes encoding proteases, as well as those involved in sulfur metabolism, phenazine synthesis, and type VI secretion systems, were notably elevated, highlighting their crucial function in feather decomposition. Integration of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) taxonomies, combined with a review of the literature, led us to propose that metabolic feather degradation involves environmental activation, reducing agent secretion, protease release, peptide/amino acid uptake, and metabolic processes. Sulfite has emerged as a critical activator of keratinase catalysis, while cysteine serves as a regulatory mediator. qRT-PCR assay results for 11 selected gene subset corroborated the RNA-seq findings. This study enhances our understanding of the transcriptomic responses of P. aeruginosa Gxun-7 to feather degradation and offers insights into potential degradation mechanisms, thereby aiding in the formulation of effective feather waste management strategies in poultry farming.

Comparative Transcriptomics Data Profiling Reveals E2F Targets as an Important Biological Pathway Overexpressed in Intellectual Disability Disorder

Intellectual disability (ID) is an early childhood neurodevelopmental disorder that is characterized by impaired intellectual functioning and adaptive behavior. It is one of the major concerns in the field of neurodevelopmental disorders across the globe. Diversified approaches have been put forward to overcome this problem. Among all these approaches, high throughput transcriptomic analysis has taken an important dimension. The identification of genes causing ID rapidly increased over the past 3 to 5 years owing to the use of sophisticated high throughput sequencing platforms. Early monitoring and preventions are much important for such disorder as their progression occurs during fetal development. This study is an attempt to identify differentially expressed genes (DEGs) and upregulated biological processes involved in development of ID patients through comparative analysis of available transcriptomics data. A total of 7 transcriptomic studies were retrieved from National Center for Biotechnology Information (NCBI) and were subjected to quality check and trimming prior to alignment. The normalization and differential expression analysis were carried out using DESeq2 and EdgeR packages of Rstudio to identify DEGs in ID. In progression of the study, functional enrichment analysis of the results obtained from both DESeq2 and EdgeR was done using gene set enrichment analysis (GSEA) tool to identify major upregulated biological processes involved in ID. Our findings concluded that monitoring the level of E2F targets, estrogen, and genes related to oxidative phosphorylation, DNA repair, and glycolysis during the developmental stage of an individual can help in the early detection of ID disorder.

The elicitor VP2 from Verticillium dahliae triggers defence response in cotton

Verticillium dahliae is a widespread and destructive soilborne vascular pathogenic fungus that causes serious diseases in dicot plants. Here, comparative transcriptome analysis showed that the number of genes upregulated in defoliating pathotype V991 was significantly higher than in the non-defoliating pathotype 1cd3-2 during the early response of cotton. Combined with analysis of the secretome during the V991-cotton interaction, an elicitor VP2 was identified, which was highly upregulated at the early stage of V991 invasion, but was barely expressed during the 1cd3-2-cotton interaction. Full-length VP2 could induce cell death in several plant species, and which was dependent on NbBAK1 but not on NbSOBIR1 in N. benthamiana. Knock-out of VP2 attenuated the pathogenicity of V991. Furthermore, overexpression of VP2 in cotton enhanced resistance to V. dahliae without causing abnormal plant growth and development. Several genes involved in JA, SA and lignin synthesis were significantly upregulated in VP2-overexpressing cotton. The contents of JA, SA, and lignin were also significantly higher than in the wild-type control. In summary, the identified elicitor VP2, recognized by the receptor in the plant membrane, triggers the cotton immune response and enhances disease resistance.

Comparative Transcriptome Analysis of High- and Low-Growth Genotypes of Eucalyptus urophylla in Response to Long-Term Nitrogen Deficiency

Nutrients play important roles in the growth and development of most plant species. However, in perennial trees, the function of nutrients in different genotypes is poorly understood. Three different nutrient levels (low, sufficient, and high nutrient levels) were applied to two contrasting Eucalyptus urophylla cultivars (a high-growth cultivar ZQUA44 and a low-growth cultivar ZQUB15), and growth and expression levels were analyzed. Although the growth traits of both genotypes under nutrient starvation treatment were much lower than under abundant nutrients, tree height, crown width, and biomass of different ZQUA44 tissues were much higher than those of ZQUB15 at all three nutrient levels. Differentially expressed genes (DEGs) clustered into six subclusters based on their expression patterns, and functional annotation showed that the DEGs involved in glutathione metabolism and flavonoid biosynthesis may be responsible for nutrient starvation across different genotypes, while the DEGs involved in carotenoid biosynthesis and starch and sucrose metabolism may have a range of functions in different genotypes. The DEGs encoding the MYB-related family may be responsible for nutrient deficiency in all genotypes, while B3 may have different functions in different genotypes. Our results demonstrate that different genotypes may form different pathways to coordinate plant survival when they face abiotic stresses.

Delivery and Transcriptome Assessment of an In Vitro Three-Dimensional Proximal Tubule Model Established by Human Kidney 2 Cells in Clinical Gelatin Sponges

The high prevalence of kidney diseases and the low identification rate of drug nephrotoxicity in preclinical studies reinforce the need for representative yet feasible renal models. Although in vitro cell-based models utilizing renal proximal tubules are widely used for kidney research, many proximal tubule cell (PTC) lines have been indicated to be less sensitive to nephrotoxins, mainly due to altered expression of transporters under a two-dimensional culture (2D) environment. Here, we selected HK-2 cells to establish a simplified three-dimensional (3D) model using gelatin sponges as scaffolds. In addition to cell viability and morphology, we conducted a comprehensive transcriptome comparison and correlation analysis of 2D and 3D cultured HK-2 cells to native human PTCs. Our 3D model displayed stable and long-term growth with a tubule-like morphology and demonstrated a more comparable gene expression profile to native human PTCs compared to the 2D model. Many missing or low expressions of major genes involved in PTC transport and metabolic processes were restored, which is crucial for successful nephrotoxicity prediction. Consequently, we established a cost-effective yet more representative model for in vivo PTC studies and presented a comprehensive transcriptome analysis for the systematic characterization of PTC lines.

Antennae-enriched expression of candidate odorant degrading enzyme genes in the turnip aphid, Lipaphis erysimi

Aphids heavily rely on their olfactory system for foraging behavior. Odorant-degrading enzymes (ODEs) are essential in preserving the olfactory acuity of aphids by removing redundant odorants in the antennae. Certain enzymes within this group stand out as being enriched and/or biased expressed in the antennae, such as carboxylesterases (CXEs), cytochrome P450 (CYPs), glutathione S-transferases (GSTs), and UDP-glycosyltransferases (UGTs). Here, we performed a comparative transcriptome analysis of antennae and body tissue to isolate the antennal ODE genes of turnip aphid Lipaphis erysimi. A dataset of one CXE, seven CYPs, two GSTs, and five UGTs enriched in the antennae was identified and subjected to sequence analysis. Furthermore, qRT-PCR analyses showed that 13 ODE genes (LeCXE6, LeCYP4c1, LeCYP6a2, LeCYP6a13, LeCYP6a14.2, LeCYP6k1, LeCYP18a1, LeGST1, LeUGT1-7, LeUGT2B7, LeUGT2B13, LeUGT2C1.1, and LeUGT2C1.2) were specifically or significantly elevated in antennal tissues. Among these antennae-enriched ODEs, LeCYP4c1, LeCYP6a2, LeCYP6a13, LeCYP6a14.2, LeCYP18a1, LeUGT2B7, and LeUGT2B13 were found to exhibit significantly higher expression levels in alate aphids compared to apterous and nymph aphids, suggesting their putative role in detecting new host plant location. The results presented in this study highlight the identification and expression of ODE genes in L. erysimi, paving the path to investigate their functional role in odorant degradation during the olfactory processes.

Comparative Transcriptome Analysis of Two Chrysomycin-Producing Wild-Type and Mutant Strains of Streptomyces sp. 891

Chrysomycin A (CA), a promising antibiotic agent, usually coexists with two analog chrysomycins B (CB) and C (CC) produced by several wild-type (WT) Streptomyces strains. With the aim to increase CA production, UV mutagenesis-based breeding had been employed on a marine-derived strain Streptomyces sp. 891 in our previous study and afforded an improved strain 891-B6 with enhanced CA yield. By comparative transcriptome analysis, significant differences in chrysomycin BGC-related gene expression between the WT strain 891 and the mutant strain 891-B6 were unveiled in the current study. Among 25 up-regulated genes in mutant 891-B6, chryA, chryB, chryC, chryF, chryG, chryK, chryP, and chryQ, responsible for the biosynthesis of benzonaphthopyranone aglycone, and chryD, chryE, and chryU in charge of production of its deoxyglycoside, were characterized. Furthermore, the expression of genes chryOII, chryOIII, and chryOIV responsible for the formation of 8-vinyl in CA from 8-ethyl in CB were greatly enhanced in strain 891-B6. These findings provide molecular mechanisms for increased yield of CA and decreased yield of CB for mutant 891-B6, which has potential application in industrial CA production.

Comparative Transcriptome Analysis of Two Sugarcane Cultivars in Response to Paclobutrazol Treatment

Sugarcane is an important crop across the globe, and the rapid multiplication of excellent cultivars is an important object of the sugarcane industry. As one of the plant growth regulators, paclobutrazol (PBZ) has been frequently used in the tissue culture of sugarcane seedlings. However, little is known about the molecular mechanisms of response to PBZ in this crop. Here, we performed a comparative transcriptome analysis between sensitive (LC05-136) and non-sensitive (GGZ001) sugarcane cultivars treated by PBZ at three time points (0 d, 10 d, and 30 d) using RNA sequencing (RNA-Seq). The results showed that approximately 70.36 Mb of clean data for each sample were generated and assembled into 239,212 unigenes. A total of 6108 and 4404 differentially expressed genes (DEGs) were identified within the sensitive and non-sensitive sugarcane cultivars, respectively. Among them, DEGs in LC05-136 were most significantly enriched in the photosynthesis and valine, leucine and isoleucine degradation pathways, while in GGZ001, DEGs associated with ion channels and plant-pathogen interaction were mainly observed. Notably, many interesting genes, including those encoding putative regulators, key components of photosynthesis, amino acids degradation and glutamatergic synapse, were identified, revealing their importance in the response of sugarcane to PBZ. Furthermore, the expressions of sixteen selected DEGs were tested by quantitative reverse transcription PCR (RT-qPCR), confirming the reliability of the RNA-seq data used in this study. These results provide valuable information regarding the transcriptome changes in sugarcane treated by PBZ and provide an insight into understanding the molecular mechanisms underlying the resistance to PBZ in sugarcane.

Comparative Gene Expression Analysis Reveals Similarities and Differences of Chronic Myeloid Leukemia Phases

Chronic myeloid leukemia (CML) is a slowly progressing blood cancer that primarily affects elderly people. Without successful treatment, CML progressively develops from the chronic phase through the accelerated phase to the blast crisis, and ultimately to death. Nowadays, the availability of targeted tyrosine kinase inhibitor (TKI) therapies has led to long-term disease control for the vast majority of patients. Nevertheless, there are still patients that do not respond well enough to TKI therapies and available targeted therapies are also less efficient for patients in accelerated phase or blast crises. Thus, a more detailed characterization of molecular alterations that distinguish the different CML phases is still very important. We performed an in-depth bioinformatics analysis of publicly available gene expression profiles of the three CML phases. Pairwise comparisons revealed many differentially expressed genes that formed a characteristic gene expression signature, which clearly distinguished the three CML phases. Signaling pathway expression patterns were very similar between the three phases but differed strongly in the number of affected genes, which increased with the phase. Still, significant alterations of MAPK, VEGF, PI3K-Akt, adherens junction and cytokine receptor interaction signaling distinguished specific phases. Our study also suggests that one can consider the phase-wise CML development as a three rather than a two-step process. This is in accordance with the phase-specific expression behavior of 24 potential major regulators that we predicted by a network-based approach. Several of these genes are known to be involved in the accumulation of additional mutations, alterations of immune responses, deregulation of signaling pathways or may have an impact on treatment response and survival. Importantly, some of these genes have already been reported in relation to CML (e.g., AURKB, AZU1, HLA-B, HLA-DMB, PF4) and others have been found to play important roles in different leukemias (e.g., CDCA3, RPL18A, PRG3, TLX3). In addition, increased expression of BCL2 in the accelerated and blast phase indicates that venetoclax could be a potential treatment option. Moreover, a characteristic signaling pathway signature with increased expression of cytokine and ECM receptor interaction pathway genes distinguished imatinib-resistant patients from each individual CML phase. Overall, our comparative analysis contributes to an in-depth molecular characterization of similarities and differences of the CML phases and provides hints for the identification of patients that may not profit from an imatinib therapy, which could support the development of additional treatment strategies.

Comparative transcriptome analysis between patient and endometrial cancer cell lines to determine common signaling pathways and markers linked to cancer progression

The rising incidence and mortality of endometrial cancer (EC) in the United States calls for an improved understanding of the disease's progression. Current methodologies for diagnosis and treatment rely on the use of cell lines as models for tumor biology. However, due to inherent heterogeneity and differential growing environments between cell lines and tumors, these comparative studies have found little parallels in molecular signatures. As a consequence, the development and discovery of preclinical models and reliable drug targets are delayed. In this study, we established transcriptome parallels between cell lines and tumors from The Cancer Genome Atlas (TCGA) with the use of optimized normalization methods. We identified genes and signaling pathways associated with regulating the transformation and progression of EC. Specifically, the LXR/RXR activation, neuroprotective role for THOP1 in Alzheimer's disease, and glutamate receptor signaling pathways were observed to be mostly downregulated in advanced cancer stage. While some of these highlighted markers and signaling pathways are commonly found in the central nervous system (CNS), our results suggest a novel function of these genes in the periphery. Finally, our study underscores the value of implementing appropriate normalization methods in comparative studies to improve the identification of accurate and reliable markers.

Comparative Transcriptome and Weighted Gene Co-expression Network Analysis Identify Key Transcription Factors of Rosa chinensis 'Old Blush' After Exposure to a Gradual Drought Stress Followed by Recovery

Rose is one of the most fundamental ornamental crops, but its yield and quality are highly limited by drought. The key transcription factors (TFs) and co-expression networks during rose’s response to drought stress and recovery after drought stress are still limited. In this study, the transcriptomes of leaves of 2-year-old cutting seedlings of Rosa chinensis ‘Old Blush’ from three continuous droughted stages (30, 60, 90 days after full watering) and rewatering were analyzed using RNA sequencing. Weighted gene co-expression network analysis (WGCNA) was used to construct a co-expression network, which was associated with the physiological traits of drought response to discovering the hub TFs involved in drought response. More than 45 million high-quality clean reads were generated from the sample and used for comparison with the rose reference genome. A total of 46433 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that drought stress caused significant changes in signal transduction, plant hormones including ABA, auxin, brassinosteroid (BR), cytokinin, ethylene (ET), jasmonic acid (JA) and salicylic acid (SA), primary and secondary metabolism, and a certain degree of recovery after rewatering. Gene co-expression analysis identified 18 modules, in which four modules showed a high degree of correlation with physiological traits. In addition, 42 TFs including members of NACs, WRKYs, MYBs, AP2/ERFs, ARFs, and bHLHs with high connectivity in navajowhite1 and blue modules were screened. This study provides the transcriptome sequencing report of R. chinensis ‘Old Blush’ during drought stress and rewatering process. The study also identifies the response of candidate TFs to drought stress, providing guidelines for improving the drought tolerance of the rose through molecular breeding in the future.

Use of Comparative Transcriptomics Combined With Physiological Analyses to Identify Key Factors Underlying Cadmium Accumulation in Brassica juncea L

The contamination of soils with cadmium (Cd) has become a serious environmental issue that needs to be addressed. Elucidating the mechanisms underlying Cd accumulation may facilitate the development of plants that accumulate both high and low amounts of Cd. In this study, a combination of phenotypic, physiological, and comparative transcriptomic analyses was performed to investigate the effects of different Cd concentrations (0, 5, 10, 30, 50 mg/kg) on Brassica juncea L. Our results suggest that B. juncea L. seedlings had a degree of tolerance to the 5 mg/kg Cd treatment, whereas higher Cd stress (10-50 mg/kg) could suppress the growth of B. juncea L. seedlings. The contents of soluble protein, as well as MDA (malondialdehyde), were increased, but the activities of CAT (catalase) enzymes and the contents of soluble sugar and chlorophyll were decreased, when B. juncea L. was under 30 and 50 mg/kg Cd treatment. Comparative transcriptomic analysis indicated that XTH18 (xyloglucan endotransglucosylase/hydrolase enzymes), XTH22, and XTH23 were down-regulated, but PME17 (pectin methylesterases) and PME14 were up-regulated, which might contribute to cell wall integrity maintenance. Moreover, the down-regulation of HMA3 (heavy metal ATPase 3) and up-regulation of Nramp3 (natural resistance associated macrophage proteins 3), HMA2 (heavy metal ATPase 2), and Nramp1 (natural resistance associated macrophage proteins 1) might also play roles in reducing Cd toxicity in roots. Taken together, the results of our study may help to elucidate the mechanisms underlying the response of B. juncea L. to various concentrations of Cd.

Comparative Transcriptome Analysis of Salt Stress-Induced Leaf Senescence in Medicago truncatula

Leaves are the most critical portion of forage crops such as alfalfa (Medicago sativa). Leaf senescence caused by environmental stresses significantly impacts the biomass and quality of forages. To understand the molecular mechanisms and identify the key regulator of the salt stress-induced leaf senescence process, we conducted a simple and effective salt stress-induced leaf senescence assay in Medicago truncatula, which was followed by RNA-Seq analysis coupled with physiological and biochemical characterization. By comparing the observed expression data with that derived from dark-induced leaf senescence at different time points, we identified 3,001, 3,787, and 4,419 senescence-associated genes (SAGs) for salt stress-induced leaf senescence on day 2, 4, and 6, respectively. There were 1546 SAGs shared by dark and salt stress treatment across the three time points. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that the 1546 SAGs were mainly related to protein and amino acids metabolism, photosynthesis, chlorophyll metabolism, and hormone signaling during leaf senescence. Strikingly, many different transcription factors (TFs) families out of the 1546 SAGs, including NAC, bHLH, MYB, and ERF, were associated with salt stress-induced leaf senescence processes. Using the transient expression system in Nicotiana benthamiana, we verified that three functional NAC TF genes from the 1546 SAGs were related to leaf senescence. These results clarify SAGs under salt stress in M. truncatula and provide new insights and additional genetic resources for further forage crop breeding.

Comparative Transcriptome Analysis of Toxic and Non-Toxic Nassarius Communities and Identification of Genes Involved in TTX-Adaptation

Nassarius has caused serious people poisoning and death incident as a popular food due to tetrodotoxin (TTX) accumulation in their body. Understanding the genetic basis of tetrodotoxin (TTX) transformation and resistance in animals could lead to significant insights into adaptive evolution to toxins and toxin poisoning cures in medicine. Here we performed comparative transcriptome analysis for toxic and non-toxic communities in Nassarius succinctus and Nassarius variciferus to reveal their genetic expression and mutation patterns. For both species, the cellular and metabolic process, and binding and catalytic activity accounted for the top classification categories, and the toxic communities generally produced more up-regulated genes than non-toxic communities. Most unigenes and different expression genes were related to disease, e.g., heat shock protein and tissue factor pathway inhibitors, which involve detoxification and coagulation. In mutation levels, the sodium channel gene of N. succinctus had one amino acid mutation "L", which is different from that of other animals. In conclusion, the comparative transcriptome analysis of different species and populations provided an important genetic basis for adaptive evolution to toxins, health and toxin poisoning cure research for TTX in marine gastropoda mollusk. Future studies will focus on the action mechanism of the important functional gene for TTX accumulation and resistance.

Comparative Transcriptome Analysis of the Regenerating Zebrafish Telencephalon Unravels a Resource With Key Pathways During Two Early Stages and Activation of Wnt/β-Catenin Signaling at the Early Wound Healing Stage

Owing to its pronounced regenerative capacity in many tissues and organs, the zebrafish brain represents an ideal platform to understand the endogenous regeneration mechanisms that restore tissue integrity and function upon injury or disease. Although radial glial and neuronal cell populations have been characterized with respect to specific marker genes, comprehensive transcriptomic profiling of the regenerating telencephalon has not been conducted so far. Here, by processing the lesioned and unlesioned hemispheres of the telencephalon separately, we reveal the differentially expressed genes (DEGs) at the early wound healing and early proliferative stages of regeneration, i.e., 20 h post-lesion (hpl) and 3 days post-lesion (dpl), respectively. At 20 hpl, we detect a far higher number of DEGs in the lesioned hemisphere than in the unlesioned half and only 7% of all DEGs in both halves. However, this difference disappears at 3 dpl, where the lesioned and unlesioned hemispheres share 40% of all DEGs. By performing an extensive comparison of the gene expression profiles in these stages, we unravel that the lesioned hemispheres at 20 hpl and 3 dpl exhibit distinct transcriptional profiles. We further unveil a prominent activation of Wnt/β-catenin signaling at 20 hpl, returning to control level in the lesioned site at 3 dpl. Wnt/β-catenin signaling indeed appears to control a large number of genes associated primarily with the p53, apoptosis, forkhead box O (FoxO), mitogen-activated protein kinase (MAPK), and mammalian target of rapamycin (mTOR) signaling pathways specifically at 20 hpl. Based on these results, we propose that the lesioned and unlesioned hemispheres react to injury dynamically during telencephalon regeneration and that the activation of Wnt/β-catenin signaling at the early wound healing stage plays a key role in the regulation of cellular and molecular events.

Global Reprogramming of Gene Transcription in Trichoderma reesei by Overexpressing an Artificial Transcription Factor for Improved Cellulase Production and Identification of Ypr1 as an Associated Regulator

Synthetic biology studies on filamentous fungi are providing unprecedented opportunities for optimizing this important category of microbial cell factory. Artificial transcription factor can be designed and used to offer novel modes of regulation on gene transcription network. Trichoderma reesei is commonly used for cellulase production. In our previous studies, a plasmid library harboring genes encoding artificial zinc finger proteins (AZFPs) was constructed for engineering T. reesei, and the mutant strains with improved cellulase production were selected. However, the underlying mechanism by which AZFP function remain unclear. In this study, a T. reesei Rut-C30 mutant strain T. reesei U5 bearing an AZFP named as AZFP-U5 was focused, which secretes high level protein and shows significantly improved cellulase and xylanase production comparing with its parental strain. In addition, enhanced sugar release was achieved from lignocellulosic biomass using the crude cellulase from T. reesei U5. Comparative transcriptome analysis was further performed, which showed reprogramming of global gene transcription and elevated transcription of genes encoding glycoside hydrolases by overexpressing AZFP-U5. Furthermore, 15 candidate regulatory genes which showed remarkable higher transcription levels by AZFP-U5 insertion were overexpressed in T. reesei Rut-C30 to examine their effects on cellulase biosynthesis. Among these genes, TrC30_93861 (ypr1) and TrC30_74374 showed stimulating effects on filter paper activity (FPase), but deletion of these two genes did not affect cellulase activity. In addition, increased yellow pigment production in T. reesei Rut-C30 by overexpression of gene ypr1 was observed, and changes of cellulase gene transcription were revealed in the ypr1 deletion mutant, suggesting possible interaction between pigment production and cellulase gene transcription. The results in this study revealed novel aspects in regulation of cellulase gene expression by the artificial regulators. In addition, the candidate genes and processes identified in the transcriptome data can be further explored for synthetic biology design and metabolic engineering of T. reesei to enhance cellulase production.

Diversity of Cellulase-Producing Filamentous Fungi From Tibet and Transcriptomic Analysis of a Superior Cellulase Producer Trichoderma harzianum LZ117

Filamentous fungi are widely used for producing cellulolytic enzymes to degrade lignocellulosic biomass. Microbial resources from Tibet have received great attention due to the unique geographic and climatic conditions in the Qinghai-Tibet Plateau. However, studies on cellulase producing fungal strains originated from Tibet remain very limited, and so far no studies have been focused on regulation of cellulase production of the specific strains thereof. Here, filamentous fungal strains were isolated from soil, plant, and other environments in Tibet, and cellulase-producing strains were further investigated. A total of 88 filamentous fungal strains were identified, and screening of cellulase-producing fungi revealed that 16 strains affiliated with the genera Penicillium, Trichoderma, Aspergillus, and Talaromyces exhibited varying cellulolytic activities. Among these strains, T. harzianum isolate LZ117 is the most potent producer. Comparative transcriptome analysis using T. harzianum LZ117 and the control strain T. harzianum K223452 cultured on cellulose indicated an intensive modulation of gene transcription related to protein synthesis and quality control. Furthermore, transcription of xyr1 which encodes the global transcriptional activator for cellulase expression was significantly up-regulated. Transcription of cre1 and other predicted repressors controlling cellulase gene expression was decreased in T. harzianum LZ117, which may contribute to enhancing formation of primary cellulases. To our knowledge, this is the first report that the transcription landscape at the early enzyme production stage of T. harzianum was comprehensively described, and detailed analysis on modulation of transporters, regulatory proteins as well as protein synthesis and processing was presented. Our study contributes to increasing the catalog of publicly available transcriptome data from T. harzianum, and provides useful clues for unraveling the biotechnological potential of this species for lignocellulosic biorefinery.

Improving Azo Dye Decolorization Performance and Halotolerance of Pichia occidentalis A2 by Static Magnetic Field and Possible Mechanisms Through Comparative Transcriptome Analysis

A halotolerant yeast, Pichia occidentalis A2, was recently isolated that can decolorize various azo dyes. The azo dye decolorization performance of this strain was characterized, including the degradation pathway and detoxification effects of this yeast. Additionally, the effect of static magnetic field (SMF) on this decolorization process was investigated. Activities of key enzymes were analyzed to estimate the change of metabolic activity. Furthermore, possible mechanisms were analyzed through detecting differentially expressed genes between yeast A2 in the absence and presence of SMF. The results indicated that yeast A2 displayed the optimal decolorization performance when the concentrations (in g/L) of glucose, (NH4)2SO4, yeast extract, and NaCl were 4.0, 1.0, 0.1, and ≤30.0, respectively. Meanwhile, the optimal rotation speed, temperature, and pH were 160 rpm, 30°C, and 5.0, respectively. Acid Red B was decolorized and detoxified by yeast A2 through successive steps, including cleavage of the naphthalene-amidine bond, reductive deamination, oxidative deamination/desulfurization, open-loop of hydroxy-substituted naphthalene, and tricarboxylic acid cycle. The dye decolorization efficiency and halotolerance of yeast A2 were enhanced by 206.3 mT SMF. The activities of manganese peroxidase, and laccase were elevated 1.37- and 1.16-fold by 206.3 mT SMF, but lignin peroxidase activity showed little change. It was suggested from the transcriptome sequence that the enhanced halotolerance might be related to the upregulated genes encoding the enzymes or functional proteins related to intracellular synthesis and accumulation of glycerol.

Analysis of evolutionary relationships provides new clues to the origins of weedy rice

Weedy rice (WR) (Oryza sativa f. spontanea) is considered to be a pest in modern rice production systems because it competes for resources, has poor yield characteristics, and subsequently has a negative effect on rice grain yield. The evolutionary relationships among WR, landrace rice (LR), improved rice (IR) cultivars, and wild rice are largely unknown. In this study, we conducted a population genetic analysis based on neutral markers and gene haplotypes in 524 rice accessions and a comparative transcriptomic analysis using 15 representative samples. The results showed that WR populations have the highest level of genetic diversity (H e = 0.8386) and can be divided into two groups (japonica-type and indica-type). The japonica-type WR accessions from Heilongjiang province (HLJ), Jilin province (JL), Liaoning province (LN), and NX provinces clustered with the landraces grown in these same provinces. The indica-types from Jiangsu province (JS) also clustered with the indica-type landraces from JS province. Comparative transcriptome analysis of WR' IR and LR from HLJ, JL, and LN provinces showed that the WR still clustered with the LR, and that the IR lines comprise a single population. Thirty-two differentially expressed genes were shared by the IR and LR groups as well as between the IR and WR groups. Using Gene ontology (GO) analysis, we identified 19 shared GO terms in the IR and LR groups as well as between the IR and WR groups. Our results suggest that WR populations in China have diverse origins, and comparative transcriptome analysis of different types of rice from HLJ, JL, and LN provinces suggests that IR populations have become a end point in the evolution of WR, which provides a new perspective for the study of WR origins and lays a solid foundation for rice breeding.

C3orf70 Is Involved in Neural and Neurobehavioral Development

Neurogenesis is the process by which undifferentiated progenitor cells develop into mature and functional neurons. Defects in neurogenesis are associated with neurodevelopmental and neuropsychiatric disorders; therefore, elucidating the molecular mechanisms underlying neurogenesis can advance our understanding of the pathophysiology of these disorders and facilitate the discovery of novel therapeutic targets. In this study, we performed a comparative transcriptomic analysis to identify common targets of the proneural transcription factors Neurog1/2 and Ascl1 during neurogenesis of human and mouse stem cells. We successfully identified C3orf70 as a novel common target gene of Neurog1/2 and Ascl1 during neurogenesis. Using in situ hybridization, we demonstrated that c3orf70a and c3orf70b, two orthologs of C3orf70, were expressed in the midbrain and hindbrain of zebrafish larvae. We generated c3orf70 knockout zebrafish using CRISPR/Cas9 technology and demonstrated that loss of c3orf70 resulted in significantly decreased expression of the mature neuron markers elavl3 and eno2. We also found that expression of irx3b, a zebrafish ortholog of IRX3 and a midbrain/hindbrain marker, was significantly reduced in c3orf70 knockout zebrafish. Finally, we demonstrated that neurobehaviors related to circadian rhythm and altered light–dark conditions were significantly impaired in c3orf70 knockout zebrafish. These results suggest that C3orf70 is involved in neural and neurobehavioral development and that defects in C3orf70 may be associated with midbrain/hindbrain-related neurodevelopmental and neuropsychiatric disorders.

Corrigendum: The Differences between NAD-ME and NADP-ME Subtypes of C4 Photosynthesis: More than Decarboxylating Enzymes

[This corrects the article DOI: 10.3389/fpls.2016.01525.].

Comparative Transcriptomic Analysis of Streptococcus thermophilus TH1436 and TH1477 Showing Different Capability in the Use of Galactose

Streptococcus thermophilus is a species widely used in the dairy industry for its capability to rapidly ferment lactose and lower the pH. The capability to use galactose produced from lactose hydrolysis is strain dependent and most of commercial S. thermophilus strains are galactose-negative (Gal⁻), although galactose-positive (Gal⁺) would be more technologically advantageous because this feature could provide additional metabolic products and prevent galactose accumulation in foods. In this study, a next generation sequencing transcriptome approach was used to compare for the first time a Gal⁺ and a Gal⁻ strain to characterize their whole metabolism and shed light on their different properties, metabolic performance and gene regulation. Transcriptome analysis revealed that all genes of the gal operon were expressed very differently in Gal⁺ and in the Gal⁻ strains. The expression of several genes involved in mixed acid fermentation, PTS sugars transporter and stress response were found enhanced in Gal⁺. Conversely, genes related to amino acids, proteins metabolism and CRISPR associated proteins were under-expressed. In addition, the strains showed a diverse series of predicted genes controlled by the transcriptional factor catabolite control protein A (CcpA). Overall, transcriptomic analysis suggests that the Gal⁺ strain underwent a metabolic remodeling to cope with the changed environmental conditions.

The Novel Secreted Meloidogyne incognita Effector MiISE6 Targets the Host Nucleus and Facilitates Parasitism in Arabidopsis

Meloidogyne incognita is highly specialized parasite that interacts with host plants using a range of strategies. The effectors are synthesized in the esophageal glands and secreted into plant cells through a needle-like stylet during parasitism. In this study, based on RNA-seq and bioinformatics analysis, we predicted 110 putative Meloidogyne incognita effectors that contain nuclear localization signals (NLSs). Combining the Burkholderia glumae-pEDV based screening system with subcellular localization, from 20 randomly selected NLS effector candidates, we identified an effector MiISE6 that can effectively suppress B. glumae-induced cell death in Nicotiana benthamiana, targets to the nuclei of plant cells, and is highly expressed in early parasitic J2 stage. Sequence analysis showed that MiISE6 is a 157-amino acid peptide, with an OGFr_N domain and two NLS motifs. Hybridization in situ verified that MiISE6 is expressed in the subventral esophageal glands. Yeast invertase secretion assay validated the function of the signal peptide harbored in MiISE6. Transgenic Arabidopsis thaliana plants expressing MiISE6 become more susceptible to M. incognita. Inversely, the host-derived RNAi of MiISE6 of the nematode can decrease its parasitism on host. Based on transcriptome analysis of the MiISE6 transgenic Arabidopsis samples and the wild-type samples, we obtained 852 differentially expressed genes (DEGs). Integrating Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, we found that expression of MiISE6 in Arabidopsis can suppress jasmonate signaling pathway. In addition, the expression of genes related to cell wall modification and the ubiquitination proteasome pathway also have detectable changes in the transgenic plants. Results from the present study suggest that MiISE6 is involved in interaction between nematode-plant, and plays an important role during the early stages of parasitism by interfering multiple signaling pathways of plant. Moreover, we found homologs of MiISE6 in other sedentary nematodes, Meloidogyne hapla and Globodera pallida. Our experimental results provide evidence to decipher the molecular mechanisms underlying the manipulation of host immune defense responses by plant parasitic nematodes, and transcriptome data also provide useful information for further study nematode-plant interactions.

The Mechanism of Starch Over-Accumulation in Chlamydomonas reinhardtii High-Starch Mutants Identified by Comparative Transcriptome Analysis

The focus of this study was the mechanism of starch accumulation in Chlamydomonas reinhardtii high-starch mutants. Three C. reinhardtii mutants showing high-starch content were generated using gamma irradiation. When grown under nitrogen-deficient conditions, these mutants had more than twice as much starch than a wild-type control. The mechanism of starch over-accumulation in these mutants was studied with comparative transcriptome analysis. In all mutants, induction of phosphoglucomutase 1 (PGM1) expression was detected; PGM1 catalyzes the inter-conversion of glucose 1-phosphate and glucose 6-phosphate in both starch biosynthetic and glycolytic pathway. Interestingly, transcript levels of phosphoglucose isomerase 1 (PGI1), fructose 1,6-bisphosphate aldolase 1 and 2 (FBA1 and FBA2) were down-regulated in all mutants; PGI1, FBA1, and FBA2 act on downstream of glucose 6-phosphate conversion in glycolytic pathway. Therefore, down-regulations of PGI1, FBA1, and FBA2 may lead to accumulation of upstream metabolites, notably glucose 6-phosphate, resulting in induction of PGM1 expression through feed-forward regulation and that PGM1 overexpression caused starch over-accumulation in mutants. These results suggest that PGI1, FBA1, FBA2, and PGM1 correlate with each other in terms of coordinated transcriptional regulation and play central roles for starch over-accumulation in C. reinhardtii.

The Differences between NAD-ME and NADP-ME Subtypes of C4 Photosynthesis: More than Decarboxylating Enzymes

As an adaptation to changing climatic conditions that caused high rates of photorespiration, C4 plants have evolved to display higher photosynthetic efficiency than C3 plants under elevated temperature, high light intensities, and drought. The C4 plants independently evolved more than 60 times in 19 families of angiosperms to establish similar but not uniform C4 mechanisms to concentrate CO2 around the carboxylating enzyme Rubisco (ribulose bisphosphate carboxylase oxygenase). C4 photosynthesis is divided into at least two basic biochemical subtypes based on the primary decarboxylating enzymes, NAD-dependent malic enzyme (NAD-ME) and NADP-dependent malic enzyme (NADP-ME). The multiple polygenetic origins of these subtypes raise questions about the association of C4 variation between biochemical subtypes and diverse lineages. This review addresses the differences in evolutionary scenario, leaf anatomy, and especially C4 metabolic flow, C4 transporters, and cell-specific function deduced from recently reported cell-specific transcriptomic, proteomic, and metabolic analyses of NAD-ME and NADP-ME subtypes. Current omic analysis has revealed the extent to which component abundances differ between the two biochemical subtypes, leading to a better understanding of C4 photosynthetic mechanisms in NAD-ME and NADP-ME subtypes.

Comparative Transcriptome Analysis Identifies CCDC80 as a Novel Gene Associated with Pulmonary Arterial Hypertension

Pulmonary arterial hypertension (PAH) is a heterogeneous disorder associated with a progressive increase in pulmonary artery resistance and pressure. Although various therapies have been developed, the 5-year survival rate of PAH patients remains low. There is thus an important need to identify novel genes that are commonly dysregulated in PAH of various etiologies and could be used as biomarkers and/or therapeutic targets. In this study, we performed comparative transcriptome analysis of five mammalian PAH datasets downloaded from a public database. We identified 228 differentially expressed genes (DEGs) from a rat PAH model caused by inhibition of vascular endothelial growth factor receptor under hypoxic conditions, 379 DEGs from a mouse PAH model associated with systemic sclerosis, 850 DEGs from a mouse PAH model associated with schistosomiasis, 1598 DEGs from one cohort of human PAH patients, and 4260 DEGs from a second cohort of human PAH patients. Gene-by-gene comparison identified four genes that were differentially upregulated or downregulated in parallel in all five sets of DEGs. Expression of coiled-coil domain containing 80 (CCDC80) and anterior gradient two genes was significantly increased in the five datasets, whereas expression of SMAD family member six and granzyme A was significantly decreased. Weighted gene co-expression network analysis revealed a connection between CCDC80 and collagen type I alpha 1 (COL1A1) expression. To validate the function of CCDC80 in vivo, we knocked out ccdc80 in zebrafish using the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system. In vivo imaging of zebrafish expressing a fluorescent protein in endothelial cells showed that ccdc80 deletion significantly increased the diameter of the ventral artery, a vessel supplying blood to the gills. We also demonstrated that expression of col1a1 and endothelin-1 mRNA was significantly decreased in the ccdc80-knockout zebrafish. Finally, we examined Ccdc80 immunoreactivity in a rat PAHmodel and found increased expression in the hypertrophied media and adventitia of the pre-acinar pulmonary arteries (PAs) and in the thickened intima, media, and adventitia of the obstructed intra-acinar PAs. These results suggest that increased expression of CCDC80 may be involved in the pathogenesis of PAH, potentially by modulating the expression of endothelin-1 and COL1A1.

EP300 Protects from Light-Induced Retinopathy in Zebrafish

Exposure of rhodopsin to bright white light can induce photoreceptor cell damage and degeneration. However, a comprehensive understanding of the mechanisms underlying light-induced retinopathy remains elusive. In this study, we performed comparative transcriptome analysis of three rodent models of light-induced retinopathy, and we identified 37 genes that are dysregulated in all three models. Gene ontology analysis revealed that this gene set is significantly associated with a cytokine signaling axis composed of signal transducer and activator of transcription 1 and 3 (STAT1/3), interleukin 6 signal transducer (IL6ST), and oncostatin M receptor (OSMR). Furthermore, the analysis suggested that the histone acetyltransferase EP300 may be a key upstream regulator of the STAT1/3–IL6ST/OSMR axis. To examine the role of EP300 directly, we developed a larval zebrafish model of light-induced retinopathy. Using this model, we demonstrated that pharmacological inhibition of EP300 significantly increased retinal cell apoptosis, decreased photoreceptor cell outer segments, and increased proliferation of putative Müller cells upon exposure to intense light. These results suggest that EP300 may protect photoreceptor cells from light-induced damage and that activation of EP300 may be a novel therapeutic approach for the treatment of retinal degenerative diseases.

E2F4 Promotes Neuronal Regeneration and Functional Recovery after Spinal Cord Injury in Zebrafish

Mammals exhibit poor recovery after spinal cord injury (SCI), whereas non-mammalian vertebrates exhibit significant spontaneous recovery after SCI. The mechanisms underlying this difference have not been fully elucidated; therefore, the purpose of this study was to investigate these mechanisms. Using comparative transcriptome analysis, we demonstrated that genes related to cell cycle were significantly enriched in the genes specifically dysregulated in zebrafish SCI. Most of the cell cycle-related genes dysregulated in zebrafish SCI were down-regulated, possibly through activation of e2f4. Using a larval zebrafish model of SCI, we demonstrated that the recovery of locomotive function and neuronal regeneration after SCI were significantly inhibited in zebrafish treated with an E2F4 inhibitor. These results suggest that activation of e2f4 after SCI may be responsible, at least in part, for the significant recovery in zebrafish. This provides novel insight into the lack of recovery after SCI in mammals and informs potential therapeutic strategies.

Arabidopsis type B cytokinin response regulators ARR1, ARR10, and ARR12 negatively regulate plant responses to drought

In this study, we used a loss-of-function approach to elucidate the functions of three Arabidopsis type B response regulators (ARRs)--namely ARR1, ARR10, and ARR12--in regulating the Arabidopsis plant responses to drought. The arr1,10,12 triple mutant showed a significant increase in drought tolerance versus WT plants, as indicated by its higher relative water content and survival rate on drying soil. This enhanced drought tolerance of arr1,10,12 plants can be attributed to enhanced cell membrane integrity, increased anthocyanin biosynthesis, abscisic acid (ABA) hypersensitivity, and reduced stomatal aperture, but not to altered stomatal density. Further drought-tolerance tests of lower-order double and single mutants indicated that ARR1, ARR10, and ARR12 negatively and redundantly control plant responses to drought, with ARR1 appearing to bear the most critical function among the three proteins. In agreement with these findings, a comparative genome-wide analysis of the leaves of arr1,10,12 and WT plants under both normal and dehydration conditions suggested a cytokinin (CK) signaling-mediated network controlling plant adaptation to drought via many dehydration/drought- and/or ABA-responsive genes that can provide osmotic adjustment and protection to cellular and membrane structures. Expression of all three ARR genes was repressed by dehydration and ABA treatments, inferring that plants down-regulate these genes as an adaptive mechanism to survive drought. Collectively, our results demonstrate that repression of CK response, and thus CK signaling, is one of the strategies plants use to cope with water deficit, providing novel insight for the design of drought-tolerant plants by genetic engineering.

De novo assembly and comparative analysis of root transcriptomes from different varieties of Panax ginseng C. A. Meyer grown in different environments

Panax ginseng C. A. Meyer is an important traditional herb in eastern Asia. It contains ginsenosides, which are primary bioactive compounds with medicinal properties. Although ginseng has been cultivated since at least the Ming dynasty to increase production, cultivated ginseng has lower quantities of ginsenosides and lower disease resistance than ginseng grown under natural conditions. We extracted root RNA from six varieties of fifth-year P. ginseng cultivars representing four different growth conditions, and performed Illumina paired-end sequencing. In total, 163,165,706 raw reads were obtained and used to generate a de novo transcriptome that consisted of 151,763 contigs (76,336 unigenes), of which 100,648 contigs (66.3%) were successfully annotated. Differential expression analysis revealed that most differentially expressed genes (DEGs) were upregulated (246 out of 258, 95.3%) in ginseng grown under natural conditions compared with that grown under artificial conditions. These DEGs were enriched in gene ontology (GO) terms including response to stimuli and localization. In particular, some key ginsenoside biosynthesis-related genes, including HMG-CoA synthase (HMGS), mevalonate kinase (MVK), and squalene epoxidase (SE), were upregulated in wild-grown ginseng. Moreover, a high proportion of disease resistance-related genes were upregulated in wild-grown ginseng. This study is the first transcriptome analysis to compare wild-grown and cultivated ginseng, and identifies genes that may produce higher ginsenoside content and better disease resistance in the wild; these genes may have the potential to improve cultivated ginseng grown in artificial environments.

Regulatory RNAs in photosynthetic cyanobacteria

Regulatory RNAs play versatile roles in bacteria in the coordination of gene expression during various physiological processes, especially during stress adaptation. Photosynthetic bacteria use sunlight as their major energy source. Therefore, they are particularly vulnerable to the damaging effects of excess light or UV irradiation. In addition, like all bacteria, photosynthetic bacteria must adapt to limiting nutrient concentrations and abiotic and biotic stress factors. Transcriptome analyses have identified hundreds of potential regulatory small RNAs (sRNAs) in model cyanobacteria such as Synechocystis sp. PCC 6803 or Anabaena sp. PCC 7120, and in environmentally relevant genera such as Trichodesmium, Synechococcus and Prochlorococcus. Some sRNAs have been shown to actually contain μORFs and encode short proteins. Examples include the 40-amino-acid product of the sml0013 gene, which encodes the NdhP subunit of the NDH1 complex. In contrast, the functional characterization of the non-coding sRNA PsrR1 revealed that the 131 nt long sRNA controls photosynthetic functions by targeting multiple mRNAs, providing a paradigm for sRNA functions in photosynthetic bacteria. We suggest that actuatons comprise a new class of genetic elements in which an sRNA gene is inserted upstream of a coding region to modify or enable transcription of that region.

Hundreds of potentially regulatory small RNAs (sRNAs) have been identified in model cyanobacteria and, despite recent significant progress, their functional characterization is substantial work and continues to provide surprising insights of general interest.

Comparative analysis of the primary transcriptome of Synechocystis sp. PCC 6803

RNA-seq and especially differential RNA-seq-type transcriptomic analyses (dRNA-seq) are powerful analytical tools, as they not only provide insights into gene expression changes but also provide detailed information about all promoters active at a given moment, effectively giving a deep insight into the transcriptional landscape. Synechocystis sp. PCC 6803 (Synechocystis 6803) is a unicellular model cyanobacterium that is widely used in research fields from ecology, photophysiology to systems biology, modelling and biotechnology. Here, we analysed the response of the Synechocystis 6803 primary transcriptome to different, environmentally relevant stimuli. We established genome-wide maps of the transcriptional start sites active under 10 different conditions relevant for photosynthetic growth and identified 4,091 transcriptional units, which provide information about operons, 5′ and 3′ untranslated regions (UTRs). Based on a unique expression factor, we describe regulons and relevant promoter sequences at single-nucleotide resolution. Finally, we report several sRNAs with an intriguing expression pattern and therefore likely function, specific for carbon depletion (CsiR1), nitrogen depletion (NsiR4), phosphate depletion (PsiR1), iron stress (IsaR1) or photosynthesis (PsrR1). This dataset is accompanied by comprehensive information providing extensive visualization and data access to allow an easy-to-use approach for the design of experiments, the incorporation into modelling studies of the regulatory system and for comparative analyses.

Differentially expressed genes distributed over chromosomes and implicated in certain biological processes for site insertion genetically modified rice Kemingdao

Release of genetically modified (GM) plants has sparked off intensive debates worldwide partly because of concerns about potential adverse unintended effects of GM plants to the agro system and the safety of foods. In this study, with the aim of revealing the molecular basis for unintended effects of a single site insertion GM Kemingdao (KMD) rice transformed with a synthetic cry1Ab gene, and bridging unintended effects of KMD rice through clues of differentially expressed genes, comparative transcriptome analyses were performed for GM KMD rice and its parent rice of Xiushui11 (XS11). The results showed that 680 differentially expressed transcripts were identified from 30-day old seedlings of GM KMD rice. The absolute majority of these changed expression transcripts dispersed and located over all rice chromosomes, and existed physical distance on chromosome from the insertion site, while only two transcripts were found to be differentially expressed within the 21 genes located within 100 kb up and down-stream of the insertion site. Pathway and biology function analyses further revealed that differentially expressed transcripts of KMD rice were involved in certain biological processes, and mainly implicated in two types of pathways. One type was pathways implicated in plant stress/defense responses, which were considerably in coordination with the reported unintended effects of KMD rice, which were more susceptible to rice diseases compared to its parent rice XS11; the other type was pathways associated with amino acids metabolism. With this clue, new unintended effects for changes in amino acids synthesis of KMD rice leaves were successfully revealed. Such that an actual case was firstly provided for identification of unintended effects in GM plants by comparative transciptome analysis.