Differential gene expression of cotton cultivar CCRI24 during somatic embryogenesis

RNA-Seq with a novel glabrous-ZM24fl reveals some key lncRNAs and the associated targets in fiber initiation of cotton

BACKGROUND

Cotton fiber is an important natural resource for textile industry and an excellent model for cell biology study. Application of glabrous mutant cotton and high-throughput sequencing facilitates the identification of key genes and pathways for fiber development and cell differentiation and elongation. LncRNA is a type of ncRNA with more than 200 nt in length and functions in the ways of chromatin modification, transcriptional and post-transcriptional modification, and so on. However, the detailed lncRNA and associated mechanisms for fiber initiation are still unclear in cotton.

RESULTS

In this study, we used a novel glabrous mutant ZM24fl, which is endowed with higher somatic embryogenesis, and functions as an ideal receptor for cotton genetic transformation. Combined with the high-throughput sequencing, fatty acid pathway and some transcription factors such as MYB, ERF and bHLH families were identified the important roles in fiber initiation; furthermore, 3,288 lncRNAs were identified, and some differentially expressed lncRNAs were also analyzed. From the comparisons of ZM24_0 DPA vs ZM24_-2 DPA and fl_0 DPA vs ZM24_0 DPA, one common lncRNA MSTRG 2723.1 was found that function upstream of fatty acid metabolism, MBY25-mediating pathway, and pectin metabolism to regulate fiber initiation. In addition, other lncRNAs MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1 were also showed potential important roles in fiber development; and the co-expression analysis between lncRNAs and targets showed the distinct models of different lncRNAs and complicated interaction between lncRNAs in fiber development of cotton.

CONCLUSIONS

From the above results, a key lncRNA MSTRG 2723.1 was identified that might mediate some key genes transcription of fatty acid metabolism, MYB25-mediating pathway, and pectin metabolism to regulate fiber initiation of ZM24 cultivar. Co-expression analysis implied that some other important lncRNAs (e.g., MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1) were also showed the different regulatory model and interaction between them, which proposes some valuable clues for the lncRNAs associated mechanisms in fiber development.

Auxin and cytokinin mediated regulation involved in vitro organogenesis of papaya

In vitro organogenesis is a multistep process which is largely controlled by the balance between auxin and cytokinin. Previous studies revealed a complex network regulating in vitro organogenesis in Arabidopsis thaliana; however, our knowledge of the molecular mechanisms underlying de novo shoot formation in papaya (Carica papaya) remains limited. Here, we optimized multiple factors to achieve an efficient and reproducible protocol for the induction of papaya callus formation and shoot regeneration. Subsequently, we analyzed the dynamic transcriptome profiles of samples undergoing this process, identified 5381, 642, 4047, and 2386 differentially expressed genes (DEGs), including 447, 66, 350, and 263 encoding transcription factors (TFs), in four stage comparisons. The DEGs were mainly involved in phytohormone modulation and transduction processes, particularly for auxin and cytokinin. Of these, 21 and 7 candidate genes involved in the auxin and cytokinin pathways, respectively, had distinct expression patterns throughout in vitro organogenesis. Furthermore, we found two genes encoding key TFs, CpLBD19 and CpESR1, were sharply induced on callus induction medium and shoot induction medium, indicating these two TFs may serve as proxies for callus induction and shoot formation in papaya. We therefore report a regulatory network of auxin and cytokinin signaling in papaya according to the one previously modeled for Arabidopsis. Our comprehensive analyses provide insight into the early molecular regulation of callus initiation and shoot formation in papaya, and are useful for the further identification of the regulators governing in vitro organogenesis.

Extensive intraspecific gene order and gene structural variations in upland cotton cultivars

Multiple cotton genomes (diploid and tetraploid) have been assembled. However, genomic variations between cultivars of allotetraploid upland cotton (Gossypium hirsutum L.), the most widely planted cotton species in the world, remain unexplored. Here, we use single-molecule long read and Hi-C sequencing technologies to assemble genomes of the two upland cotton cultivars TM-1 and zhongmiansuo24 (ZM24). Comparisons among TM-1 and ZM24 assemblies and the genomes of the diploid ancestors reveal a large amount of genetic variations. Among them, the top three longest structural variations are located on chromosome A08 of the tetraploid upland cotton, which account for ~30% total length of this chromosome. Haplotype analyses of the mapping population derived from these two cultivars and the germplasm panel show suppressed recombination rates in this region. This study provides additional genomic resources for the community, and the identified genetic variations, especially the reduced meiotic recombination on chromosome A08, will help future breeding.

De novo assembly and comparative analysis of the transcriptome of embryogenic callus formation in bread wheat (Triticum aestivum L.)

BACKGROUND

During asexual reproduction the embryogenic callus can differentiate into a new plantlet, offering great potential for fostering in vitro culture efficiency in plants. The immature embryos (IMEs) of wheat (Triticum aestivum L.) are more easily able to generate embryogenic callus than mature embryos (MEs). To understand the molecular process of embryogenic callus formation in wheat, de novo transcriptome sequencing was used to generate transcriptome sequences from calli derived from IMEs and MEs after 3d, 6d, or 15d of culture (DC).

RESULTS

In total, 155 million high quality paired-end reads were obtained from the 6 cDNA libraries. Our de novo assembly generated 142,221 unigenes, of which 59,976 (42.17%) were annotated with a significant Blastx against nr, Pfam, Swissprot, KOG, KEGG, GO and COG/KOG databases. Comparative transcriptome analysis indicated that a total of 5194 differentially expressed genes (DEGs) were identified in the comparisons of IME vs. ME at the three stages, including 3181, 2085 and 1468 DEGs at 3, 6 and 15 DC, respectively. Of them, 283 overlapped in all the three comparisons. Furthermore, 4731 DEGs were identified in the comparisons between stages in IMEs and MEs. Functional analysis revealed that 271transcription factor (TF) genes (10 overlapped in all 3 comparisons of IME vs. ME) and 346 somatic embryogenesis related genes (SSEGs; 35 overlapped in all 3 comparisons of IME vs. ME) were differentially expressed in at least one comparison of IME vs. ME. In addition, of the 283 overlapped DEGs in the 3 comparisons of IME vs. ME, excluding the SSEGs and TFs, 39 possessed a higher rate of involvement in biological processes relating to response to stimuli, in multi-organism processes, reproductive processes and reproduction. Furthermore, 7 were simultaneously differentially expressed in the 2 comparisons between the stages in IMEs, but not MEs, suggesting that they may be related to embryogenic callus formation. The expression levels of genes, which were validated by qRT-PCR, showed a high correlation with the RNA-seq value.

CONCLUSIONS

This study provides new insights into the role of the transcriptome in embryogenic callus formation in wheat, and will serve as a valuable resource for further studies addressing embryogenic callus formation in plants.

Recent Advances on Genetic and Physiological Bases of In Vitro Somatic Embryo Formation

Somatic embryogenesis involves a broad repertoire of genes, and complex expression patterns controlled by a concerted gene regulatory network. The present work describes this regulatory network focusing on the main aspects involved, with the aim of providing a deeper insight into understanding the total reprogramming of cells into a new organism through a somatic way. To the aim, the chromatin remodeling necessary to totipotent stem cell establishment is described, as the activity of numerous transcription factors necessary to cellular totipotency reprogramming. The eliciting effects of various plant growth regulators on the induction of somatic embryogenesis is also described and put in relation with the activity of specific transcription factors. The role of programmed cell death in the process, and the related function of specific hemoglobins as anti-stress and anti-death compounds is also described. The tools for biotechnology coming from this information is highlighted in the concluding remarks.

Transcriptome profiling and digital gene expression by deep sequencing in early somatic embryogenesis of endangered medicinal Eleutherococcus senticosus Maxim

Somatic embryogenesis (SE) has been studied as a model system to understand molecular events in physiology, biochemistry, and cytology during plant embryo development. In particular, it is exceedingly difficult to access the morphological and early regulatory events in zygotic embryos. To understand the molecular mechanisms regulating early SE in Eleutherococcus senticosus Maxim., we used high-throughput RNA-Seq technology to investigate its transcriptome. We obtained 58,327,688 reads, which were assembled into 75,803 unique unigenes. To better understand their functions, the unigenes were annotated using the Clusters of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes databases. Digital gene expression libraries revealed differences in gene expression profiles at different developmental stages (embryogenic callus, yellow embryogenic callus, global embryo). We obtained a sequencing depth of >5.6 million tags per sample and identified many differentially expressed genes at various stages of SE. The initiation of SE affected gene expression in many KEGG pathways, but predominantly that in metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction. This information on the changes in the multiple pathways related to SE induction in E. senticosus Maxim. embryogenic tissue will contribute to a more comprehensive understanding of the mechanisms involved in early SE. Additionally, the differentially expressed genes may act as molecular markers and could play very important roles in the early stage of SE. The results are a comprehensive molecular biology resource for investigating SE of E. senticosus Maxim.

Extensive modulation of the transcription factor transcriptome during somatic embryogenesis in Arabidopsis thaliana

Molecular mechanisms controlling plant totipotency are largely unknown and studies on somatic embryogenesis (SE), the process through which already differentiated cells reverse their developmental program and become embryogenic, provide a unique means for deciphering molecular mechanisms controlling developmental plasticity of somatic cells. Among various factors essential for embryogenic transition of somatic cells transcription factors (TFs), crucial regulators of genetic programs, are believed to play a central role. Herein, we used quantitative real-time polymerase chain reaction (qRT-PCR) to identify TF genes affected during SE induced by in vitro culture in Arabidopsis thaliana. Expression profiles of 1,880 TFs were evaluated in the highly embryogenic Col-0 accession and the non-embryogenic tanmei/emb2757 mutant. Our study revealed 729 TFs whose expression changes during the 10-days incubation period of SE; 141 TFs displayed distinct differences in expression patterns in embryogenic versus non-embryogenic cultures. The embryo-induction stage of SE occurring during the first 5 days of culture was associated with a robust and dramatic change of the TF transcriptome characterized by the drastic up-regulation of the expression of a great majority (over 80%) of the TFs active during embryogenic culture. In contrast to SE induction, the advanced stage of embryo formation showed attenuation and stabilization of transcript levels of many TFs. In total, 519 of the SE-modulated TFs were functionally annotated and transcripts related with plant development, phytohormones and stress responses were found to be most abundant. The involvement of selected TFs in SE was verified using T-DNA insertion lines and a significantly reduced embryogenic response was found for the majority of them. This study provides comprehensive data focused on the expression of TF genes during SE and suggests directions for further research on functional genomics of SE.

Transcriptome profiling reveals auxin and cytokinin regulating somatic embryogenesis in different sister lines of cotton cultivar CCRI24

To get a broader view on the molecular mechanisms underlying somatic embryogenesis (SE) in cotton (Gossypium hirsutum L.), global analysis of cotton transcriptome dynamics during SE in different sister lines was performed using RNA-Seq. A total of 204 349 unigenes were detected by de novo assembly of the 214 977 462 Illumina reads. The quantitative reverse transcription-polymerase chain reaction (qRT-PCR) measurements were positively correlated with the RNA-Seq results for almost all the tested genes (R(2)  = 0.841, correlation was significant at the 0.01 level). Different phytohormone (auxin and cytokinin) concentration ratios in medium and the endogenous content changes of these two phytohormones at two stages in different sister lines suggested the roles of auxin and cytokinin during cotton SE. On the basis of global gene regulation of phytohormone-related genes, numerous genes from all the differentially expressed transcripts were involved in auxin and cytokinin biosynthesis and signal transduction pathways. Analyses of differentially expressed genes that were involved in these pathways revealed the substantial changes in gene type and abundance between two sister lines. Isolation, cloning and silencing/overexpressing the genes that revealed remarkable up- or down-expression during cotton SE were important. Furthermore, auxin and cytokinin play a primary role in SE, but potential cross-talk with each other or other factors remains unclear.

Wuschel overexpression promotes somatic embryogenesis and induces organogenesis in cotton (Gossypium hirsutum L.) tissues cultured in vitro

This work shows that overexpression of the WUS gene from Arabidopsis enhanced the expression of embryogenic competence and triggered organogenesis from some cells of the regenerated embryo-like structures. Agrobacterium-mediated genetic transformation of cotton was described in the late 1980s, but is still time consuming and largely genotype dependant due to poor regeneration. To help solve this bottleneck, we over-expressed the WUSCHEL (WUS) gene, a homeobox transcription factor cloned in Arabidopsis thaliana, known to stimulate organogenesis and/or somatic embryogenesis in Arabidopsis tissues cultured in vitro. The AtWUS gene alone, and AtWUS gene fused to the GFP marker were compared to the GFP gene alone and to an empty construct used as a control. Somatic embryogenesis was improved in WUS expressed calli, as the percentage of explants giving rise to embryogenic tissues was significantly higher (×3) when WUS gene was over-expressed than in the control. An interesting result was that WUS embryogenic lines evolved in green embryo-like structures giving rise to ectopic organogenesis never observed in any of our previous transformation experiments. Using our standard in vitro culture protocol, the overexpression of AtWUS in tissues of a recalcitrant variety did not result in the production of regenerated plants. This achievement will still require the optimization of other non-genetic factors, such as the balance of exogenous phytohormones. However, our results suggest that targeted expression of the WUS gene is a promising strategy to improve gene transfer in recalcitrant cotton cultivars.

A chimeric arabinogalactan protein promotes somatic embryogenesis in cotton cell culture

Arabinogalactan proteins (AGPs) are a family of extracellular plant proteoglycans implicated in many aspects of plant growth and development, including in vitro somatic embryogenesis (SE). We found that specific AGPs were produced by cotton (Gossypium hirsutum) calli undergoing SE and that when these AGPs were isolated and incorporated into tissue culture medium, cotton SE was promoted. When the AGPs were partly or fully deglycosylated, SE-promoting activity was not diminished. Testing of AGPs separated by reverse-phase high-performance liquid chromatography revealed that the SE-promoting activity resided in a hydrophobic fraction. We cloned a full-length complementary DNA (cotton PHYTOCYANIN-LIKE ARABINOGALACTAN-PROTEIN1 [GhPLA1]) that encoded the protein backbone of an AGP in the active fraction. It has a chimeric structure comprising an amino-terminal signal sequence, a phytocyanin-like domain, an AGP-like domain, and a hydrophobic carboxyl-terminal domain. Recombinant production of GhPLA1 in tobacco (Nicotiana tabacum) cells enabled us to purify and analyze a single glycosylated AGP and to demonstrate that this chimeric AGP promotes cotton SE. Furthermore, the nonglycosylated phytocyanin-like domain from GhPLA1, which was bacterially produced, also promoted SE, indicating that the glycosylated AGP domain was unnecessary for in vitro activity.

Transcript profiling reveals complex auxin signalling pathway and transcription regulation involved in dedifferentiation and redifferentiation during somatic embryogenesis in cotton

BACKGROUND

Somatic embryogenesis (SE), by which somatic cells of higher plants can dedifferentiate and reorganize into new plants, is a notable illustration of cell totipotency. However, the precise molecular mechanisms regulating SE remain unclear. To characterize the molecular events of this unique process, transcriptome analysis, in combination with biochemical and histological approaches, were conducted in cotton, a typical plant species in SE. Genome-wide profiling of gene expression allowed the identification of novel molecular markers characteristic of this developmental process.

RESULTS

RNA-Seq was used to identify 5,076 differentially expressed genes during cotton SE. Expression profile and functional assignments of these genes indicated significant transcriptional complexity during this process, associated with morphological, histological changes and endogenous indole-3-acetic acid (IAA) alteration. Bioinformatics analysis showed that the genes were enriched for basic processes such as metabolic pathways and biosynthesis of secondary metabolites. Unigenes were abundant for the functions of protein binding and hydrolase activity. Transcription factor-encoding genes were found to be differentially regulated during SE. The complex pathways of auxin abundance, transport and response with differentially regulated genes revealed that the auxin-related transcripts belonged to IAA biosynthesis, indole-3-butyric acid (IBA) metabolism, IAA conjugate metabolism, auxin transport, auxin-responsive protein/indoleacetic acid-induced protein (Aux/IAA), auxin response factor (ARF), small auxin-up RNA (SAUR), Aux/IAA degradation, and other auxin-related proteins, which allow an intricate system of auxin utilization to achieve multiple purposes in SE. Quantitative real-time PCR (qRT-PCR) was performed on selected genes with different expression patterns and functional assignments were made to demonstrate the utility of RNA-Seq for gene expression profiles during cotton SE.

CONCLUSION

We report here the first comprehensive analysis of transcriptome dynamics that may serve as a gene expression profile blueprint in cotton SE. Our main goal was to adapt the RNA-Seq technology to this notable development process and to analyse the gene expression profile. Complex auxin signalling pathway and transcription regulation were highlighted. Together with biochemical and histological approaches, this study provides comprehensive gene expression data sets for cotton SE that serve as an important platform resource for further functional studies in plant embryogenesis.

Identification of genes up-regulated during somatic embryogenesis of cucumber

Somatic embryogenesis is a method of plant regeneration, but it can also be used as a model to study plant development. A normalized library of cDNA fragments representing genes up-regulated after the induction of somatic embryogenesis in cucumber suspension cultures was constructed using the suppression subtractive hybridization technique. Candidate cDNA fragments (119) were classified according to their similarity to genes encoding known proteins and the presence of potential functional domains. Of the translation products with homology to known proteins, about 23% were possibly involved in metabolism, 13% represented proteins with a probable role in cellular communication and signal transduction, about 12% were likely to participate in protein synthesis, while around 10% were potential transcription factors. The genes corresponding to four of the cDNAs were subsequently analyzed in more detail: CsSEF2, CsSEM1 and CsSESTK1 encoding putative transcription factors or co-activators, and CsSECAD1 encoding cinnamyl alcohol dehydrogenase. Full-length cDNAs were isolated and analyzed. RT-PCR confirmed the up-regulation of these genes after the induction of somatic embryogenesis and showed the presence of their transcripts in other tissues. The in situ localization of transcripts of the CsSEF2 and CsSEM1 genes demonstrated that signalling in somatic embryo tissues involving these factors is concentrated in the cotyledon primordia and roots.

Organogenic nodule formation in hop: a tool to study morphogenesis in plants with biotechnological and medicinal applications

The usage of Humulus lupulus for brewing increased the demand for high-quality plant material. Simultaneously, hop has been used in traditional medicine and recently recognized with anticancer and anti-infective properties. Tissue culture techniques have been reported for a wide range of species, and open the prospect for propagation of disease-free, genetically uniform and massive amounts of plants in vitro. Moreover, the development of large-scale culture methods using bioreactors enables the industrial production of secondary metabolites. Reliable and efficient tissue culture protocol for shoot regeneration through organogenic nodule formation was established for hop. The present review describes the histological, and biochemical changes occurring during this morphogenic process, together with an analysis of transcriptional and metabolic profiles. We also discuss the existence of common molecular factors among three different morphogenic processes: organogenic nodules and somatic embryogenesis, which strictly speaking depend exclusively on intrinsic developmental reprogramming, and legume nitrogen-fixing root nodules, which arises in response to symbiosis. The review of the key factors that participate in hop nodule organogenesis and the comparison with other morphogenic processes may have merit as a study presenting recent advances in complex molecular networks occurring during morphogenesis and together, these provide a rich framework for biotechnology applications.