Search

Introduction: Low temperature inhibits the growth of most microorganisms. However, some microbes can grow well in a low temperature, even a freezing temperature. Methods: In this study, the mechanisms conferring cold resistance in the cryophylactic yeast Metschnikowia (M.) pulcherrima MS612, an isolate of the epidermis of ice grapes, were investigated based on comparative transcriptome analysis. Results: A total of 6018 genes and 374 differentially expressed genes (> 2-fold, p < 0.05) were identified using RNA-Seq. The differentially expressed genes were mainly involved in carbohydrate and energy metabolism, transport mechanisms, antifreeze protection, lipid synthesis, and signal transduction. M. pulcherrima MS612 maintained normal growth at low temperature (5°C) by enhancing energy metabolism, sterol synthesis, metal ion homeostasis, amino acid and MDR transport, while increased synthesis of glycerol and proline transport to improve its resistance to the freezing temperature (-5°C). Furthermore, cAMPPKA and ERAD signaling pathways contribute to resist the low temperature and the freezing temperature, respectively. Conclusion: This study provides new insights into cold resistance in cryophylactic microorganisms for maneuvering various metabolism to resist different cold environment. [ABSTRACT FROM AUTHOR]

Abstract Panax japonicus, a traditional medicinal plant from the Araliaceae family, produces bioactive triterpenes with health benefits. In Traditional Chinese Medicine, its roots have been used, but the chemical basis of its medicinal use is unclear, particularly regarding the metabolism and regulation of triterpene saponin biosynthesis. This study employed an integrative approach using Ultra Performance Liquid Chromatography (UPLC) and transcriptome analysis. Our UPLC analysis showed that ginsenoside Ro and chikusetsusaponin IVa were mainly detected in P. japonicus root. Subsequently, a comparative transcriptome analysis of four P. japonicus tissues (roots, leaves, flowers and fruits) was conducted using Illumina sequencing. As a result, 90,985 unigenes were functionally annotated from a total of 211,650 assembled non-redundant transcripts. Among these, 42,829 unigenes were annotated in NR database. Tissue-specific gene analysis revealed that roots had the highest number of specifically expressed unigenes (11,832). The majority of these unigenes were associated with metabolic processes. Additionally, tissue expression patterns analysis for three common transcription factor families indicated that WRKY family genes showed a significantly root-specific expression pattern, potentially playing a role in triterpene saponin biosynthesis. Notably, we investigated the expression profiles of genes related to the biosynthesis of triterpene saponins and found that four genes, ACCT, HMGS, HMGR and SS, encoding key enzymes in triterpene saponins biosynthesis pathway, were primarily expressed in the root. Overall, our study provides a set of P. japonicus tissue transcriptome data, which will aid in the discovery of triterpene saponin biosynthetic genes and offers valuable genetic information for this medical plant.

BackgroundCross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes.ResultsHere, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues.ConclusionsIn summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.

We generated 73 transcriptomic data of water buffalo, which were integrated with publicly available data in this species, yielding a large dataset of 355 samples representing 20 major tissue categories. We established a multi-tissue gene expression atlas of water buffalo. Furthermore, by comparing them with 4866 cattle transcriptomic data from the cattle genotype–tissue expression atlas (CattleGTEx), we found that the transcriptomes of the two species exhibited conservation in their overall gene expression patterns, tissue-specific gene expression and house-keeping gene expression. We further identified conserved and divergent expression genes between the two species, with the largest number of differentially expressed genes found in the skin, which may be related to structural and functional differences in the skin of the two species. This work provides a source of functional annotation of the buffalo genome and lays the foundations for future genetic and evolutionary studies in water buffalo. [ABSTRACT FROM AUTHOR]

Abstract Background Cross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes. Results Here, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues. Conclusions In summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.

IntroductionFat-tail in sheep is considered as an important energy reservoir to provide energy as a survival buffer during harsh challenges. However, fat-tail is losing its importance in modern sheep industry systems and thin-tailed breeds are more desirable. Using comparative transcriptome analysis to compare fat-tail tissue between fat- and thin-tailed sheep breeds provides a valuable approach to study the complex genetic factors associated with fat-tail development. However, transcriptomic studies often suffer from issues with reproducibility, which can be improved by integrating multiple studies based on a meta-analysis.MethodsHence, for the first time, an RNA-Seq meta-analysis on sheep fat-tail transcriptomes was performed using six publicly available datasets.Results and discussionA total of 500 genes (221 up-regulated, 279 down-regulated) were identified as differentially expressed genes (DEGs). A jackknife sensitivity analysis confirmed the robustness of the DEGs. Moreover, QTL and functional enrichment analysis reinforced the importance of the DEGs in the underlying molecular mechanisms of fat deposition. Protein-protein interactions (PPIs) network analysis revealed the functional interactions among the DEGs and the subsequent sub-network analysis led to identify six functional sub-networks. According to the results of the network analysis, down-regulated DEGs in green and pink sub-networks (like collagen subunits IV, V, and VI, integrins 1 and 2, SCD, SCD5, ELOVL6, ACLY, SLC27A2, and LPIN1) may impair lipolysis or fatty acid oxidation and cause fat accumulation in tail. On the other hand, up-regulated DEGs, especially those are presented in green and pink sub-networks (like IL6, RBP4, LEPR, PAI-1, EPHX1, HSD11B1, and FMO2), might contribute to a network controlling fat accumulation in the tail of sheep breed through mediating adipogenesis and fatty acid biosynthesis. Our results highlighted a set of known and novel genes/pathways associated with fat-tail development, which could improve the understanding of molecular mechanisms behind fat deposition in sheep fat-tail.

Bacterial canker of tomato is a systemic disease caused by Clavibacter michiganensis (Cm), which poses a grave threat to tomato production worldwide. Towards the identification of genes underlying resistance to Cm infection, the transcriptome of the resistant inbred backcross line IBL2353 carrying the Rcm2.0 locus derived from Solanum habrochaites LA407 and the susceptible Solanum lycopersicum line Ohio88119 was comparatively analyzed after Cm inoculation, and the analysis focused on the genes with different expression patterns between resistant and susceptible lines. Gene ontology (GO) analysis revealed that top terms of differentially expressed genes comprised ubiquitin protein ligases, transcription factors, and receptor kinases. Then we screened out some genes which are potentially associated with the defense response against Cm infection in IBL2353 including the wall-associated receptor kinase-like 20 (WAKL20), and virus-induced gene silencing showed it contributes resistance to Cm infection. In addition to Cm-induced genes related to resistance, the expression of eight homologs from six susceptibility (S) gene families was analyzed. These putative resistance and susceptibility genes are valuable resources for molecular resistance breeding and contribute to the development of new control methods in tomato. [ABSTRACT FROM AUTHOR]

High-throughput sequencing technology has been widely used and continuously optimized. As one kind of high-throughput sequencing technology, RNA sequencing (RNA-seq) can effectively detect complete transcriptome information, which brings prospects for plant researches. Conifers are a family of Coniferopsida, mainly including Cupressaceae, Taxodiaceae, Pinaceae and Taxaceae trees. Presently, the whole genome sequencing (WGS) of many crops and broadleaf trees have been completed, while that of conifer species are limited due to huge data volume and complicated repeat sequences. Therefore, RNA-seq is one of the most important and practical methods for exploring conifers genes. This review introduced the analysis processes of RNA-seq, and highlighted the applications and prospects of RNA-seq in conifers molecular genetics, which might provide references for future genetic improvement and molecular breeding in conifer species. Key message: The review presented RNA-seq applications in conifers, including comparative transcriptome, exploring functional genes, developing molecular markers, and detecting miRNAs. Moreover, the future perspectives were discussed. [ABSTRACT FROM AUTHOR]

Abstract Background The pine wood nematode (PWN), Bursaphelenchus xylophilus, is a devastating pathogen of many Pinus species in China. The aim of this study was to understand the interactive molecular mechanism of PWN and its host by comparing differentially expressed genes and candidate effectors from three transcriptomes of B. xylophilus at different infection stages. Results In total, 62, 69 and 46 candidate effectors were identified in three transcriptomes (2.5 h postinfection, 6, 12 and 24 h postinoculation and 6 and 15 d postinfection, respectively). In addition to uncharacterized pioneers, other candidate effectors were involved in the degradation of host tissues, suppression of host defenses, targeting plant signaling pathways, feeding and detoxification, which helped B. xylophilus survive successfully in the host. Seven candidate effectors were identified in both our study and the B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was identified in all three transcriptomes. These common candidate effectors were upregulated at infection stages, and one of them suppressed pathogen-associated molecular pattern (PAMP) PsXEG1-triggered cell death in Nicotiana benthamiana. Conclusions The results indicated that B. xylophilus secreted various candidate effectors, and some of them continued to function throughout all infection stages. These various candidate effectors were important to B. xylophilus infection and survival, and they functioned in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, promoting detoxification and playing virulence functions). The present results provide valuable resources for in-depth research on the pathogenesis of B. xylophilus from the perspective of effectors.

Recent researches reveal that during oocyte maturation, species‐specific molecular profile exists and has important functional roles. However, molecular differences between pig (a larger animal model for human reproduction) and human mature oocytes remain unknown. Here, by comparative transcriptome analyses of single‐cell RNA‐seq data, we aimed to identify the common and unique maternal factors and associated biological processes between in vivo and in vitro matured pig oocytes, and between in vitro matured human and pig oocytes. Annotated protein‐coding mRNAs were identified in pig in vivo (11,147) and in vitro (11,997), and human in vitro (14,491) MII oocytes, respectively. For in vivo and in vitro derived pig MII oocytes, 10,551 annotated maternal mRNAs were common, mainly enriched in signalling pathways such as cell cycle, oocyte meiosis, microtubule cytoskeleton, MAPK, RNA processing/binding. Besides, in vivo (596) and in vitro (1446) pig MII‐specific mRNAs and their involved signalling pathways (in vivo: Bmp, calcium‐mediated signalling, PI3K‐Akt; in vitro: growth factor activity, JAK‐STAT, cytokine‐cytokine receptor interaction, calcium signalling pathway) were also found. As for in vitro derived human and pig MII oocytes, 10,285 annotated mRNAs were common, enriched in a variety of signalling pathways (cell cycle, oocyte meiosis, microtubule, AMPK, RNA splicing, protein serine/threonine kinase activity, etc). In vitro MII‐specific mRNAs were found for humans (4206) and pigs (1712), which were also enriched in species‐specific signalling pathways (humans: golgi‐related terms, transcription repressor and hormone activity; pigs: ATP biosynthetic process, G protein‐coupled peptide receptor activity, animoacyl‐tRNA biosynthesis), respectively. These findings improve our understanding on oocyte maturation, and also the limitations of pig model for human oocyte maturation and fertilization. [ABSTRACT FROM AUTHOR]

Bacterial canker of tomato is a systemic disease caused by Clavibacter michiganensis (Cm), which poses a grave threat to tomato production worldwide. Towards the identification of genes underlying resistance to Cm infection, the transcriptome of the resistant inbred backcross line IBL2353 carrying the Rcm2.0 locus derived from Solanum habrochaites LA407 and the susceptible Solanum lycopersicum line Ohio88119 was comparatively analyzed after Cm inoculation, and the analysis focused on the genes with different expression patterns between resistant and susceptible lines. Gene ontology (GO) analysis revealed that top terms of differentially expressed genes comprised ubiquitin protein ligases, transcription factors, and receptor kinases. Then we screened out some genes which are potentially associated with the defense response against Cm infection in IBL2353 including the wall-associated receptor kinase-like 20 (WAKL20), and virus-induced gene silencing showed it contributes resistance to Cm infection. In addition to Cm-induced genes related to resistance, the expression of eight homologs from six susceptibility (S) gene families was analyzed. These putative resistance and susceptibility genes are valuable resources for molecular resistance breeding and contribute to the development of new control methods in tomato.

Polysaccharides separated from Lentinula edodes are well known for their medicinal properties. However, the precise molecular mechanisms of polysaccharide biosynthesis in L. edodes remain unclear. In this study, the fruiting bodies of L. edodes in four developmental stages with significant differences in polysaccharide yield were collected, and the characteristics of polysaccharides were studied. De novo sequencing and comparative transcriptomic analysis were performed by using high-throughput Illumina RNA-sequencing. KS1P30, KS2P30, KS3P30, and KS4P30 were obtained from the four developmental stages, respectively, by hot water extraction and 30% ethanol precipitation. These four polysaccharides had good immune activity in vitro; all of them were β-glucopyranose with a high molecular weight. Glucose was the main monosaccharide component of these polysaccharides. High-quality clean reads (57.88, 53.17, 53.28, and 47.56 million for different growth stages) and mapping ratios ranging from 84.75 to 90.11% were obtained. In total, 11,493 (96.56%) unigenes and 18,924 (97.46%) transcripts were successfully annotated in five public databases. The biosynthetic pathway and related genes of LEFP30 were mined. The molecular mechanism of LEFP30 yield change in the different developmental stages was predicted. The results provide some insights into the possible mechanisms involved in the biosynthetic pathway of this kind of polysaccharide in L. edodes fruiting bodies. They also indicate that candidate genes can be used as important resources for biotechnology and molecular breeding to regulate L. edodes fruiting body polysaccharide biosynthesis.

Polysaccharides separated from Lentinula edodes are well known for their medicinal properties. However, the precise molecular mechanisms of polysaccharide biosynthesis in L. edodes remain unclear. In this study, the fruiting bodies of L. edodes in four developmental stages with significant differences in polysaccharide yield were collected, and the characteristics of polysaccharides were studied. De novo sequencing and comparative transcriptomic analysis were performed by using high-throughput Illumina RNA-sequencing. KS1P30, KS2P30, KS3P30, and KS4P30 were obtained from the four developmental stages, respectively, by hot water extraction and 30% ethanol precipitation. These four polysaccharides had good immune activity in vitro ; all of them were β-glucopyranose with a high molecular weight. Glucose was the main monosaccharide component of these polysaccharides. High-quality clean reads (57.88, 53.17, 53.28, and 47.56 million for different growth stages) and mapping ratios ranging from 84.75 to 90.11% were obtained. In total, 11,493 (96.56%) unigenes and 18,924 (97.46%) transcripts were successfully annotated in five public databases. The biosynthetic pathway and related genes of LEFP30 were mined. The molecular mechanism of LEFP30 yield change in the different developmental stages was predicted. The results provide some insights into the possible mechanisms involved in the biosynthetic pathway of this kind of polysaccharide in L. edodes fruiting bodies. They also indicate that candidate genes can be used as important resources for biotechnology and molecular breeding to regulate L. edodes fruiting body polysaccharide biosynthesis. [ABSTRACT FROM AUTHOR]

Abstract Background Hericium erinaceus, also known as lion’s mane mushroom, is a widely distributed edible and medicinal fungus in Asian countries. H. erinaceus harbors diverse bioactive metabolites with anticancer, immunomodulating, anti-inflammatory, antimicrobial, antihypertensive, antidiabetic and neuroprotective properties. Although the chemical synthesis processes of these bioactive metabolites are known, the biosynthetic processes remain unknown. Results In this study, we obtained the transcriptomes of six H. erinaceus strains using next-generation RNA sequencing and investigated the characteristics of the transcriptomes and biosynthesis of bioactive compounds, especially polysaccharides. The transcriptomes ranged in size from 46.58 to 58.14 Mb, with the number of unigenes ranging from 20,902 to 37,259 across the six H. erinaceus strains. Approximately 60% of the unigenes were successfully annotated by comparing sequences against different databases, including the nonredundant (NR), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), clusters of orthologous groups for eukaryotic complete genomes (KOG) and Swiss-Prot databases. Most of the transcripts were putatively involved in signal transduction, carbohydrate metabolism, translation, transport and catabolism, and amino acid metabolism. Genes involved in polysaccharide biosynthesis were identified, and these genes encoded phosphoglucomutase (PGM), glucose phosphate isomerase (PGI), UDP-glucose pyrophosphorylase (UGP), glycoside hydrolase family proteins, glycosyltransferase family proteins and other proteins. Moreover, the putative pathway for the intracellular polysaccharide biosynthesis of H. erinaceus was analyzed. Additionally, the open reading frames (ORFs) and simple sequence repeats (SSRs) were predicted from the transcriptome data of the six strains. Conclusions Overall, the present study may facilitate the discovery of polysaccharide biosynthesis processes in H. erinaceus and provide useful information for exploring the secondary metabolites in other members of the Basidiomycetes genus.

Background: The pine wood nematode (PWN), Bursaphelenchus xylophilus, is a devastating pathogen of many Pinus species in China. The aim of this study was to understand the interactive molecular mechanism of PWN and its host by comparing differentially expressed genes and candidate effectors from three transcriptomes of B. xylophilus at different infection stages. Results: In total, 62, 69 and 46 candidate effectors were identified in three transcriptomes (2.5 h postinfection, 6, 12 and 24 h postinoculation and 6 and 15 d postinfection, respectively). In addition to uncharacterized pioneers, other candidate effectors were involved in the degradation of host tissues, suppression of host defenses, targeting plant signaling pathways, feeding and detoxification, which helped B. xylophilus survive successfully in the host. Seven candidate effectors were identified in both our study and the B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was identified in all three transcriptomes. These common candidate effectors were upregulated at infection stages, and one of them suppressed pathogen-associated molecular pattern (PAMP) PsXEG1-triggered cell death in Nicotiana benthamiana. Conclusions: The results indicated that B. xylophilus secreted various candidate effectors, and some of them continued to function throughout all infection stages. These various candidate effectors were important to B. xylophilus infection and survival, and they functioned in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, promoting detoxification and playing virulence functions). The present results provide valuable resources for in-depth research on the pathogenesis of B. xylophilus from the perspective of effectors. [ABSTRACT FROM AUTHOR]

Background: Histone deacetylase (HDAC)-1, a Class-I HDAC family member, forms three types of complexes, the nucleosome remodeling deacetylase, Sin3, and CoREST complexes with the specific corepressor components chromodomain-helicase-DNA-binding protein 3 (Mi2/CHD-3), Sin3, and REST corepressor 1 (RCOR1), respectively, in humans. Objective: To elucidate the functional relationships among the three transcriptional corepressors during embryogenesis. Methods: The activities of HDA-1, LET-418, SIN-3, and SPR-1, the homologs of HDAC-1, Mi2, Sin3, and RCOR1 in Caenorhabditis elegans during embryogenesis were investigated through measurement of relative mRNA expression levels and embryonic lethality given either gene knockdown or deletion. Additionally, the terminal phenotypes of each knockdown and mutant embryo were observed using a differential-interference contrast microscope. Finally, the functional relationships among the three corepressors were examined through genetic interactions and transcriptome analyses. Results: Here, we report that each of the corepressors LET-418, SIN-3, and SPR-1 are expressed and have essential roles in C. elegans embryonic development. Our terminal phenotype observations of single mutants further implied that LET-418, SIN-3, and SPR-1 play similar roles in promoting advancement to the middle and late embryonic stages. Combined analysis of genetic interactions and gene ontology of these corepressors indicate a prominent overlapping role among SIN-3, SPR-1, and LET-418 and between SIN-3 and SPR-1. Conclusion: Our findings suggest that the class-I HDAC-1 corepressors LET-418, SIN-3, and SPR-1 may cooperatively regulate the expression levels of some genes during C. elegans embryogenesis or may have some similar roles but functioning independently within a specific cell. [ABSTRACT FROM AUTHOR]

The perennial woody plants of citrus are one of the most important fruit crops in the world and largely depends on arbuscular mycorrhizal symbiosis (AMS) to obtain essential nutrients from soil. However, the molecular aspects of AMS in citrus and perennial woody plants in general have largely been understudied. We used RNA-sequencing to identify differentially expressed genes in roots of Poncirus trifoliata upon mycorrhization by the AM fungus Glomus versiforme and evaluated their conservation by comparative transcriptome analyses with four herbaceous model plants. We identified 282 differentially expressed genes in P. trifoliata, including orthologs of 21 genes with characterized roles in AMS and 83 genes that are considered to be conserved in AM-host plants. Comparative transcriptome analysis revealed a 'core set' of 156 genes from P. trifoliata whose orthologous genes from at least three of the five species also exhibited similar transcriptional changes during AMS. Functional analysis of one of these conserved AM-induced genes, a 3-keto-acyl-ACP reductase (FatG) involved in fatty acid biosynthesis, confirmed its involvement in AMS in Medicago truncatula. Our results identify a core transcriptional program for AMS that is largely conserved between P. trifoliata and other plants. The comparative transcriptomics approach adds to previous phylogenomics studies to identify conserved genes required for AMS. [ABSTRACT FROM AUTHOR]

Background: Rust fungi constitute the largest group of plant fungal pathogens. However, a paucity of data, including genomic sequences, transcriptome sequences, and associated molecular markers, hinders the development of inhibitory compounds and prevents their analysis from an evolutionary perspective. Gymnosporangium yamadae and G. asiaticum are two closely related rust fungal species, which are ecologically and economically important pathogens that cause apple rust and pear rust, respectively, proved to be devastating to orchards. In this study, we investigated the transcriptomes of these two Gymnosporangium species during the telial stage of their lifecycles. The aim of this study was to understand the evolutionary patterns of these two related fungi and to identify genes that developed by selection. Results: The transcriptomes of G. yamadae and G. asiaticum were generated from a mixture of RNA from three biological replicates of each species. We obtained 49,318 and 54,742 transcripts, with N50 values of 1957 and 1664, for G. yamadae and G. asiaticum, respectively. We also identified a repertoire of candidate effectors and other gene families associated with pathogenicity. A total of 4947 pairs of putative orthologues between the two species were identified. Estimation of the non-synonymous/synonymous substitution rate ratios for these orthologues identified 116 pairs with Ka/Ks values greater than1 that are under positive selection and 170 pairs with Ka/Ks values of 1 that are under neutral selection, whereas the remaining 4661 genes are subjected to purifying selection. We estimate that the divergence time between the two species is approximately 5.2 Mya. Conclusion: This study constitutes a de novo assembly and comparative analysis between the transcriptomes of the two rust species G. yamadae and G. asiaticum. The results identified several orthologous genes, and many expressed genes were identified by annotation. Our analysis of Ka/Ks ratios identified orthologous genes subjected to positive or purifying selection. An evolutionary analysis of these two species provided a relatively precise divergence time. Overall, the information obtained in this study increases the genetic resources available for research on the genetic diversity of the Gymnosporangium genus. [ABSTRACT FROM AUTHOR]

Background: Somatic embryogenesis is a notable illustration of cell totipotency, by which somatic cells undergo dedifferentiation and then differentiate into somatic embryos. Our previous work demonstrated that pretreatment of immature zygotic embryos with 0.5 M sucrose solution for 72 h efficiently induced somatic embryo initiation in camphor tree. To better understand the molecular basis of somatic embryogenesis induced by osmotic stress, de novo transcriptome sequencing of three tissues of camphor tree (immature zygotic embryos, sucrose-pretreated immature zygotic embryos, and somatic embryos induced from sucrose-pretreated zygotic embryos) were conducted using Illumina Hiseq 2000 platform. Results: A total of 30.70 G high quality clean reads were obtained from cDNA libraries of the three samples. The overall de novo assembly of cDNA sequence data generated 205592 transcripts, with an average length of 998 bp. 114229 unigenes (55.56 % of all transcripts) with an average length of 680 bp were annotated with gene descriptions, gene ontology terms or metabolic pathways based on Blastx search against Nr, Nt, Swissprot, GO, COG/KOG, and KEGG databases. CEGMA software identified 237 out of 248 ultra-conserved core proteins as 'complete' in the transcriptome assembly, showing a completeness of 95.6 %. A total of 897 genes previously annotated to be potentially involved in somatic embryogenesis were identified. Comparative transcriptome analysis showed that a total of 3335 genes were differentially expressed in the three samples. The differentially expressed genes were divided into six groups based on K-means clustering. Expression level analysis of 52 somatic embryogenesis-related genes indicated a high correlation between RNA-seq and qRT-PCR data. Gene enrichment analysis showed significantly differential expression of genes responding to stress and stimulus. Conclusions: The present work reported a de novo transcriptome assembly and global analysis focused on gene expression changes during initiation and formation of somatic embryos in camphor tree. Differential expression of somatic embryogenesis-related genes indicates that sucrose induced somatic embryogenesis may share or partly share the mechanisms of somatic embryogenesis induced by plant hormones. This study provides comprehensive transcript information and gene expression data for camphor tree. It could also serve as an important platform resource for further functional studies in plant embryogenesis. [ABSTRACT FROM AUTHOR]

The Asian yellow pond turtle (Mauremys mutica) displays temperature-dependent sex determination (TSD), in which incubation temperature during embryonic development determines the sexual fate of the individual. However, the mechanism of the sex determination/differentiation of Mauremys mutica remains a mystery. Here, we first analyzed the temperature-specific gonadal transcriptomes of Mauremys mutica prior to gonad formation and gonads during the thermosensitive period. We uncovered a list of candidates that respond to temperature stimuli enriched in several categories, such as heat shock protein family members dnajb6a , dnaja4 , hspa8 and hsp90aa1 , temperature sensor genes mmp17 and mmp28 , and putative novel temperature-responsive genes tmco6 , gria3 and eif3f. Notably, striking differences were identified in the expression profiles of genes underlying sexual development, such as tex15 , insr , igf1r , cirbp , esr1 , dmrt2 and Serpinh1. Moreover, we analyzed the similarity and divergence of the timecourse of gene expression among Mauremys mutica and two other reported TSD turtles (Trachemys scripta and Chrysemys picta). The shared genes revealed the common gonad-specific regulatory mechanisms existing in these three TSD turtles that initiate their sexual development. Therefore, our findings could provide basic data to elucidate the mechanisms of sex determination/differentiation of M. mutica , even contributing to further understanding of these mechanisms in other TSD turtles. [Display omitted] • A list of candidates genes that respond to temperature stimuli were identified in Mauremys mutica. • Striking differences were identified in the expression profiles of genes underlying sexual development of Mauremys mutica. • Eight shared genes associate with common gonad-specific regulatory mechanisms were identified in the three TSD turtles. [ABSTRACT FROM AUTHOR]