- Open Access
Comparative transcriptome study provides insights into acquisition of embryogenic ability in upland cotton during somatic embryogenesis
- SUN Ruibin †1,
- TIAN Ruiping †1,
- MA Dan 1,
- WANG Shaohui 1 and
- LIU Chuanliang 1, 2Email authorView ORCID ID profile
© The Author(s) 2018
- Received: 14 May 2018
- Accepted: 20 August 2018
- Published: 31 August 2018
The conversion from non-embryogenic callus (NEC) to embryogenic callus (EC) is the key bottleneck step in regeneration of upland cotton (Gossypium hirsutum), and hinders the transgenic breeding of upland cotton. To investigate molecular mechanisms underlying acquisition of embryogenic potential during this process, comparation analysis of transcriptome dynamics between two upland cotton cultivars with different somatic embryogenesis abilities was conducted.
Differentially expressed genes involved in the transformation from NEC to EC were detected in the two different cultivars. Principal component analysis based on DEGs showed that the NEC tissues of the two cultivars were highly heterogeneous, whereas the derived EC tissues were similar, which suggested the homogeneousness of EC between different lines. In the highly embryogenic cultivar CCRI 24, more of these genes were down-regulated, whereas, in the recalcitrant cultivar CCRI 12, more were up-regulated. Bioinformatics analysis on these DEGs showed that the vast majority of differentially expressed genes were enriched in metabolism and secondary metabolites biosynthesis pathways. Flavonoid biosynthesis and phenylpropanoid biosynthesis pathways were enriched in both cultivars, and the associated genes were down-regulated more in CCRI 24 than in CCRI 12. We deduced that vigorous secondary metabolism in CCRI 12 may hinder primary metabolism, resulting in tardiness of cell differentiation. Interestingly, genes involved in the plant hormone signal transduction pathway were enriched in the recalcitrant cultivar CCRI 12, but not in CCRI 24, suggesting more radical regulation of hormone signal transduction in the recalcitrant cultivar. Signal transduction rather than biosynthesis of plant hormones is more likely to be the determining factor triggering NEC to EC transition in recalcitrant cotton lines. Transcription factor encoding genes showed differential regulation between two cultivars.
Our study provides valuable information about the molecular mechanism of conversion from NEC to EC in cotton and allows for identification of novel genes involved. By comparing transcriptome changes in transformation from NEC to EC between the two cultivars, we identified 46 transcripts that may contribute to initiating embryogenic shift.
- Upland cotton
- Non-embryogenic callus (NEC)
- Embryogenic callus (EC)
- Somatic embryogenesis
Upland cotton (Gossypium hirsutum) is an important economic crop, grown in many countries, as the major natural fiber, an important raw material for textile, and nutrient-rich edible oil that can be extracted from cotton seeds (Sunilkumar et al. 2006; Shi et al. 2006). For many years, many efforts have been made to improve cotton yield, fiber quality, stress tolerance, and other traits. In recent decades, genetic modification technologies have been widely applied in the breeding of many plants, including cotton. Incorporating exogenous genes into explants by Agrobacterium-mediation or gene gun bombardment, and regeneration of transformed cells through somatic embryogenesis (SE) is the main and the most effective method used to produce transgenic cotton (Liu et al. 2014; Zhang et al. 2000; Leelavathi et al. 2004). Genetic modification technologies have great potential to improve desirable cotton traits, including yield, fiber quality, and pest, disease, drought, and herbicide resistance. However, cotton genetic modification is severely constrained by the poor regenerative capacity of cotton plants (Zhang et al. 2011; Shang et al. 2009). Many factors, such as the explants, culture medium components, and light and temperature during culture, affect the regeneration of cotton cells (Wu et al. 2004; Chi et al. 2004). Nevertheless, the inherent genetic background is the essential limiting factor (Trolinder and Chen 1989; Dong 1991; Deo et al. 2010). Although lots of effort have been made to optimize the somatic regeneration culture system, only very few upland cotton cultivars can be regenerated from Agrobacterium-mediated or biolistic genetic transformed cells (Leelavathi et al. 2004; Shang et al. 2009; Zhang et al. 2000; Rajeswari et al. 2010; Xie et al. 2007; Zhang et al. 2009; Sun et al. 2009; Dong 1991; Liu et al. 2004; Han et al. 2009).
Plant regeneration through SE is a long and complicated process. First, the transformed cells dedifferentiate into calli on the culture medium. After several weeks, an embryogenic callus (EC) is acquired and then differentiates into an embryo, which ultimately grows into a plant. The process of converting a non-embryogenic callus (NEC) into an EC is the key bottleneck step (Zhang et al. 2011; Shang et al. 2009; Wu et al. 2004; Rajeswari et al. 2010). SE is a complex dynamic process that involves many intracellular metabolic changes and is influenced by various external environment factors. The underlying temporal-specific expression of genes plays an important role in this process. Analyzing the transcriptional activity of genes during the transition from NEC to EC may help to reveal the molecular mechanisms involved in the acquirement of embryogenic potential, which, in turn, may provide ideas to facilitate the induction of EC to further promote regeneration of a wider range of cotton cultivars for genetic modification. Differential-display of reverse transcription polymerase chain reaction (PCR), cDNA-amplified fragment length polymorphisms, suppression subtractive hybridization, and microarrays have been used to study gene expression during cotton SE, and a number of specifically or differentially expressed genes (DEGs) associated with cotton SE have been identified (Wu et al. 2009; Zeng et al. 2006; Zhu et al. 2008; Leng et al. 2007), but the data were insufficient or incomplete. Next-generation sequencing combined with bioinformatics analyses have proved to be powerful tools for more comprehensive genome and transcriptome analysis. Using RNA sequencing-based digital gene expression technology, Yang et al. (2012) reported the first genome-wide gene expression profiling of different SE stages (from hypocotyl to cotyledon embryos) in cotton. A total of 3 329 genes linked to various pathways were found to be differential expression during the transition from NEC to EC, and 137 genes were specifically switched on or off in this process. Xu et al. (2013) performed RNA sequencing to study the transcriptome changes in two upland cotton lines with different somatic embryogenic capacities in response to SE induction. Focusing on transcripts associated with phytohormone metabolism, seven auxin-related transcripts and nine cytokinin-related transcripts were specifically found in embryonic calli of a highly embryogenic line. All of these existing transcriptome studies on cotton SE by RNA sequencing are based on a de novo assembled transcriptome, which may contain many mistakes as the high complexity of plant genome. Now that genome sequencing of upland cotton has been completed (Li et al. 2015), transcriptome sequences can be obtained more easily and accurately, making it conducive for transcriptome analysis. Here we conducted a comparative analysis on transcriptomes dynamics during NEC to EC conversion process between a recalcitrant embryogenic cotton cultivar and a highly embryogenic cotton cultivar by RNA sequencing. Different from previous studies, we mapped sequencing reads to the published genome as reference to assemble transcriptome. Besides, we also obtain EC from a recalcitrant cultivar, that the transcriptome difference between NEC and EC was compared, rather than just compare the difference of NEC before inductive treatment with after treatment.
Plant materials and induction of callus
CCRI 12 and CCRI 24, two commercial planting varieties of upland cotton in China, were used in our study. CCRI 12 is an excellent cultivar that has superb disease resistance and high yield, and produces high quality fiber. Because of its superior comprehensive performance, CCRI 12 has large acreage, and is an important material in cotton breeding. However, CCRI 12 shows recalcitrant embryogenesis from explant in genetic transformation, which hinders its application in transgenic breeding. CCRI 24 is the other major cultivar with high embryogenesis efficiency in China, and it has been developed into a model cultivar for cotton genetic transformation. We clarified the embryogenic differences between CCRI 12 and CCRI 24 by a regeneration test as described below.
Seeds were delinted using concentrated sulfuric acid, then sterilized in 3% hydrogen peroxide solution for 12 h, and rinsed with sterile water for five times before sowing on solid medium containing Murashige & Skoog (MS) (Murashige and Skoog 1962) macroelements salts, sucrose, and agar. After culturing at 28 °C under a 16 h light/8 h dark photoperiod, the hypocotyls of 7-day-old seedlings were cut into 1 cm long segments, then placed on calli induction medium containing MS, 0.1 mg·L− 1 3-indole acetic acid (IAA), 0.1 mg·L− 1 kinetin (KT), 0.1 mg·L− 1 2,4-dichlorophenoxyacetic acid (2,4-D), 30 g·L− 1 glucose, and 2.5 g·L− 1 Gelrite, adjusted to pH 6.0. NEC of each cultivar were sampled from three different seedlings as biological replicates after 40-day culture under the same conditions as seedling growth, the number of NEC was counted before sampling. Then, the calli were transferred to subculture medium, in which all the exogenous hormone levels were reduced to 1/10 of those in the calli induction medium. EC was sampled from three different seedlings as biological replicates when the EC of each cultivar formed finally, after subculturing several times. The number of EC was counted before sampling. All the calli tissues were frozen in liquid nitrogen and stored at − 80 °C for RNA extraction. The NEC and EC of CCRI 12 and CCRI 24 represented four different sample types - named, respectively, by AN (NEC of CCRI 12), AE (EC of CCRI 12), BN (NEC of CCRI 24), and BE (EC of CCRI 24), with each sample type containing three samples from different individuals as biological replicates. These 12 samples were used for subsequent RNA library construction and sequencing.
We quantified the embryogenic ability by computing the regeneration rate of NEC (the ratio of NEC to explant), and the regeneration rate of EC (the ratio of EC to NEC).
Total RNA was extracted using a TIANGEN RNAprep Pure Plant Kit, following the manufacturer’s protocol (Tiangen n.d.). The quality of total RNA was tested using a NanoDrop 2000 instrument (Thermo Scientific) and agarose gel electrophoresis. The RNA samples were used for library construction and sequencing.
Library construction and RNA sequencing
We provided 20 μg of total RNA to the Beijing Genomics Institute (BGI) where the libraries were constructed and sequenced on an Illumina HiSeq 2000 platform. Briefly, mRNA was enriched using oligo-dT magnetic beads and fragmented into short fragments. cDNA was synthesized using the mRNA fragments as templates primed by random hexamers. Then, the cDNA fragments were purified for end repairing and single nucleotide A (adenine) addition, then connected with adapters. Fragments of about 250 base pair (bp) were selected by agarose gel electrophoresis, and purified as templates for PCR amplification. After PCR, qualification and quantification of the libraries were conducted using an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System. The qualified libraries were sequenced using an Illumina HiSeq™ 2000 sequencer.
Genome mapping and bioinformatics analysis of (DEGs)
The raw sequence data were filtered to remove adapter sequences, reads with > 10% of unknown bases, or > 50% of low quality bases. The filtered clean data were used for further analysis. The BWA (Li and Durbin 2009) and Bowtie2 (Langmead and Salzberg 2012) software were used to map the clean reads to the published Gossypium hirsutum TM-1 genome (Li et al. 2015) and gene reference sequences, respectively. The mapping information was used to detect specific and common genes between the different samples. The expression levels of the identified genes and isoforms were quantified using the RSEM (RNA sequencing by Expectation Maximization) (Li and Dewey 2011) software package. Fragments per kilobase of transcript per million mapped reads (FPKM) values were calculated to evaluate expression levels of genes or transcripts. Based on the results, correlation and principal component analyses were conducted pair-wise between samples. The NOIseq software package (Tarazona et al. 2012) was used to screen differential expression analysis based on the noise distribution model, genes with fold change ≥2, and divergence probability ≥0.8 were recognized as DEGs. Genes with similar expression patterns usually have some functional correlation. A hierarchical clustering analysis was performed of DEGs among the sample pairs using Cluster and Java Treeview software (Saldanha 2004). The Gene Ontology (GO) (Ashburner et al. 2000) database and canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto 2000) and plant gene datasets selected were used for GO functional and pathway enrichment analyses of the DEGs.
Quantitative real-time PCR (qRT-PCR) was carried out to verify the expression profile results calculated from the RNA sequencing data. We selected 13 genes that were differentially expressed during SE for qRT-PCR validation. The PCR primers were designed using the NCBI online Primer-BLAST. RNA samples were treated to remove genomic DNA and reversely transcribed using a TransGen AH341 kit. The obtained cDNA templates were diluted five times prior to amplification. The qRT-PCR reaction system contained 1 μL cDNA templates, 0.5 μL each of 10 μmol·L-1 forward and reverse primers, 2× TransStart Top Green qPCR Super Mix, add ddH2O up to 20 μL. All the reactions had three technical replicates and were run on a Roche Light Cycler 480 instrument. The amplification program was as follows: pre-incubation at 95 °C for 5 min, followed by 40 cycles of denaturation at 95 °C for 10 s, annealing at 60 °C for 10 s and extension at 72 °C for 15 s. A melting curve was constructed to detect the specificity of the products. GhHis3 (GenBank: AF024716) was used as the reference gene. The relative expression values were computed by the 2−ΔCT method.
Obtaining and characterizing the NEC and EC tissues of CCRI 12 and CCRI 24
Regeneration rates in different stages of SE in CCRI 12 and CCRI 24
Average frequency /%
Summary of the RNA sequencing data
After filtering the raw data, we obtained a total of 598.2 million clean reads of 90 bp in length for further analysis. Each library produced at least 4.37 G bases clean data, which were sufficient for gene expression profiling. Samples which showed bad correlations with other duplicates were removed. An average of 79.5% clean reads were mapped to the reference genome, and an average of 71.3% clean reads were mapped within known genes. Integrating sequencing data of twelve samples, a total of 60 908 genes (79% of the annotated genes in the reference genome) were detected. 80.7% of these genes were common in the four libraries. 1 121, 850, 803, 1 136 genes were specifically expressed in AE, AN, BE, and BN, respectively.
qRT-PCR validation of gene expression from the RNA sequencing data
Principal component analysis (PCA) shows homogeneousness of EC between different cultivars
Global analysis of differential gene expression during the transformation from NEC to EC
TF encoding genes expression analysis during the NEC to EC transformation
Summary of differentially expressed transcription factor (TF) encoding genes compared between two cultivars and different callus states
Sample type pair
Differentially expressed TF encoding genes
Up-regulated TF encoding genes
Down-regulated TF encoding genes
AN vs AE
BN vs BE
Although CCRI 12 and CCRI 24 have different SE capabilities, we obtained embryogenic calli from both cultivars, suggesting that both of them have genetic potentials for SE and embryonic callus formation competence. Under the same culture conditions, different SE efficiencies of CCRI 12 and CCRI 24 are likely due to differences in gene expression (Xu et al. 2013). Based on the transcriptome profile obtained by RNA sequencing, we compared gene expression patterns between NEC and EC of CCRI 12 and CCRI 24 to dissect the underlying molecular mechanism of transformation from NEC to EC, which is the key restriction in SE. In the two cultivars, most of the DEGs were implicated and enriched significantly in metabolism and secondary metabolites biosynthesis, especially flavonoid biosynthesis and phenylpropanoid biosynthesis pathways which are closely related. Flavonoids are synthesized from phenylpropanoid derivatives by condensation with malonyl-CoA (De Vries 2000; Vogt 2010; Saito et al. 2013). Genes encoding key enzymes of phenylpropanoid biosynthesis, such as dihydroflavonol-4-reductase, chalcone synthase, chalcone isomerase, flavonoid 3′-hydroxylase, and flavonoid 3′, 5′-hydroxylase, were down-regulated from NEC to EC. Similarly, genes encoding key enzymes in flavonoid biosynthesis, such as phenylalanine ammonialyase, cinnamate-4-hydroxylase, and 4-coumarate-CoA ligase, were also down-regulated from NEC to EC. These results are consistent with previous studies on metabolites, which reported that the content of flavonoids and phenolic compounds in EC were lower than NEC (Zhang et al. 1992; Li et al. 1994; Wu et al. 1999). Flavonoids are a large class of plant secondary metabolites with many functions, including pigmentation, antioxidant activity, and response to biological and abiotic stresses (De Vries 2000; Saito et al. 2013). Down-regulation of genes encoding key enzymse in flavonoid biosynthesis and phenylpropanoid biosynthesis may indicate a decreased secondary metabolism in EC, while major primary metabolic enhancement may favor cell differentiation (Zhang et al. 1992). More of these involved genes were down-regulated in CCRI 24 (120 in 137) comparing with CCRI 12 (104 in 133) during the transformation from NEC to EC. For example, all the seven putative 2-oxoglutarate and Fe (II)-dependent oxygenase superfamily protein encoding genes in CCRI 24 were down-regulated, whereas three of them showed highly up-regulated (log2Ratio values were 2, 4, and 9, respectively) in CCRI 12. The NAD (P)-binding Rossmann-fold superfamily protein encoding genes were up-regulated significantly in CCRI 12, while down-regulated in CCRI 24. Our results showed that flavonoid and phenylpropanoid biosynthesis processes were down-regulated more greatly in CCRI 24 than CCRI 12. More vigorous secondary metabolism in CCRI 12 may hinder the main primary metabolism required, and result in tardiness of cell differentiation and gain of embryonic callus.
Plant hormones were found to participate in the regulation of SE (Fehér et al. 2003; Jiménez 2005). Exogenous plant hormones such as IAA, 2,4-D, and KT are commonly added to medium to induce and regulate SE (Liu et al. 2014). In the present study, from NEC to EC, the plant hormone signal transduction pathway was enriched in the recalcitrant cultivar CCRI 12, but not in CCRI 24. While comparing NEC with EC, 148 plant hormone signal transduction-related genes were differentially expressed in CCRI 12, containing 77 up-regulated genes and 71 down-regulated genes. One hundred nine plant hormone signal transduction-related genes were differentially expressed in CCRI 24, containing 34 up-regulated genes and 75 down-regulated genes. Concerning hormone signal transduction-related genes, CCRI 12 showed significant more DEGs up-regulated than CCRI 24. Besides, all up-regulated DEGs showed significant higher ratio of up-regulation in CCRI 12 than CCRI 24. All of these indicate a more dramatical regulation in CCRI 12, implying that hormone signal transduction may be more crucial for the recalcitrant cultivar CCRI 12 to obtain embryonic competence than CCRI 24. This seems to contradict a previous study, which found that signal transduction of plant hormones was more likely to be the determining factor for recalcitrant cotton lines to obtain EC than biosynthesis (Xu et al. 2013).
Principal component analysis showed that NEC tissues of CCRI 12 and CCRI 24 were highly heterogeneous, whereas the EC tissues of CCRI 12 and CCRI 24 were similar. It may mean that the EC tissues derived from different (even any) genotypes have similar molecular status, despite the difference of NEC. Therefore, we hypothesized that the differentially expressed transcripts detected in AE vs BE comparison were attributed to the genetic background differences between cultivars, and were unrelated to SE. While the differentially expressed transcripts detected in AN vs BN comparison were more likely to be SE-related. Of the 490 DEGs detected in the AN vs BN comparison, 175 DEGs showed differential expressions during the transformation from NEC to EC in both CCRI 12 and CCRI 24 (except for eight genes detected in AE vs BE comparison). All of these DEGs showed same expression pattern between CCRI 12 and CCRI 24 during the transformation, with 155 DEGs up-regulated and 20 DEGs down-regulated. All the 155 up-regulated genes had higher expression levels in BN than in AN, including five genes that were not detected in AN. The expression of these genes was positively correlated with the formation of EC, indicating that the higher expression levels of these genes in BN may account for the more effective formation of EC. Of the 20 down-regulated DEGs, 13 genes had higher expression levels in AN than BN, and other 7 genes had higher expression levels in BN than AN. The expression of these genes was negatively correlated with the formation of EC.
SE marks the end of the explant gene expression pattern and the beginning of the embryonic gene expression pattern (Han et al. 2009). We assumed that some transcripts may be critical or favorable for SE, deficiency of these transcripts in CCRI 12 may hinder the initiation of transformation from NEC to EC. By comparing AN with BN, we attempted to identify genes that contributed to initiation of the embryogenic transformation. Thirty transcripts were detected in BN but not in AN. Among them, one gene encoded each of the following proteins, GSTU43, disease resistance-responsive family protein, a phosphatidyltransferase, phloem protein 2-A1, nucleotide binding protein, DnaJ subfamily C member 2, a heat shock protein, WD-repeat protein, cysteine-rich protein kinase, seed storage 2S albumin superfamily protein, SNARE-like superfamily protein, and ELF4-like protein, 5 genes encoded uncharacterized proteins, 7 encoded hypothetical proteins (five were L484_002552), and 6 had no annotations. Most of these transcripts were up-regulated in CCRI 12, implying they may be crucial for effective embryogenesis, and accumulation of these transcripts may promote the initiation of embryogenic shift. Among these transcripts, 7 transcripts, including 5 for hypothetical protein L484_002552, one for uncharacterized protein TCM_032689 and one for phosphatidyltransferase, showed opposite expression patterns, up-regulated in CCRI 12 and down-regulated in CCRI 24. We considered that these transcripts may be significantly up-regulated to activate the reprogramming for embryonic transformation and subsequently down-regulated, which implied that these transcripts may play important roles specifically in the initiation of transformation from NEC to EC. Ten up-regulated transcripts in both CCRI 12 and CCRI 24 and 6 up-regulated transcripts in CCRI 12 may be required throughout the NEC to EC transformation to maintain the transformation process. Six transcripts were only detected in CCRI 24, including four transcripts encoding cysteine-rich protein kinase, seed storage 2S albumin superfamily protein, SNARE-like superfamily protein, and ELF4-like protein respectively, and two encoding uncharacterized proteins. The transcript that encoded seed storage 2S albumin superfamily protein was uniquely detected in BN, which is interesting because this transcript showed down-regulated in BE comparing with BN, implying its specific role in initiation of transformation. This transcript may be significantly up-regulated in BN to activate the initiation of embryonic transformation, then down-regulated through the whole transformation process, as indicated by its absence in AN and AE.
In contrast, some transcripts may be repressive for SE, accumulation of these transcripts in CCRI 12 may also hinder the NEC initiate the transformation to EC. Sixteen transcripts were detected in AN but not in BN, and most of them were absent in both BN and BE. Thirteen of the transcripts existed only in CCRI 12, including one transcript that uniquely existed in AN. Among these thirteen transcripts, three were down-regulated in AN compared with AE, one encoded calcium-binding EF-hand family protein, one encoded fiber protein Fb17, and one encoded multi-vesicular body protein. It is possible that these three transcripts may account for the recalcitrance of CCRI 12, especially the fiber protein Fb17 coding one which was detected only in AN. The down-regulation of these genes may be conducive to the removal of suppression and the subsequent initiation of embryonic transformation.
In current study, during the transformation from NEC to EC, we sampled only at the initial and end points. Because many biological processes are gradual, and the transition from NEC to EC spans a long time, the two end-timepoint samplings may have missed gene expression information in the intermediate processes. For example, some transcripts may appear in a special middle period, but not at either end. Sampling at more timepoints during the process will provide more information about the gene expression dynamics that takes place from NEC to EC. More precise transcriptome profile will provide insight into more actual cellular regulation during the transformation process.
Comparative analysis of transcriptome profiles derived from RNA sequencing indicate that EC tissues are highly homogenous, despite between genotypes with different SE abilities. The transformation from NEC to EC involves extensive metabolism changes. In SE-recalcitrant cultivars, excessive secondary metabolism, especially flavonoid and phenylpropanoid biosynthesis, may hinder cell differentiation and delay the emergence of embryogenic attributes. Plant hormones play important roles in EC induction, and we found that the recalcitrant genotype likely had poorer ability to perceive the exogenous hormones and signaling. Take published upland cotton genomes as references, the present transcriptome analysis study provided new insight into molecular mechanism of the intractable conversion of NEC to EC in cotton, and may help for related functional genes mining.
We gratefully acknowledge Li WJ for her encourage and support.
This work was supported by National Science and Technology Major Project (2016ZX08010004), China.
Availability of data and materials
Plant materials are available from the corresponding author on reasonable request. The datasets supporting the results of this article are included within the article and its additional files.
Tian RP designed and carried out the experiment, Sun RB analyzed the data and wrote draft manuscript. Ma D and Liu CL supervised the research and gave final approval of the version to be published. All the authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ashburner M, Ball CA, Blake JA, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9. https://doi.org/10.1038/75556 View ArticlePubMedPubMed CentralGoogle Scholar
- Chi JN, Li XH, Wang XF, et al. The factors affecting somatic embryogenesis and plant regeneration in cotton. Acta Gossypii Sinica. 2004;16(1):55–61.Google Scholar
- De Vries GE. Flavonoid biosynthetic pathway. Trends Plant Sci. 2000;5(1):7. https://doi.org/10.1016/S1360-1385(99)01519-8 Google Scholar
- Deo PC, Tyagi AP, Taylor M, et al. Factors affecting somatic embryogenesis and transformation in modern plant breeding. S Pac J Natl Appl Sci. 2010;28(1):27–40. https://doi.org/10.1071/SP10002 View ArticleGoogle Scholar
- Dong HZ. Cotton somatic embryogenesis of different geotypes. J Laiyang Agric Coll. 1991;8(2):97–101.Google Scholar
- Fehér A, Pasternak TP, Dudits D. Transition of somatic plant cells to an embryogenic state. Plant Cell Tissue Organ Cult. 2003;74(3):201–28. https://doi.org/10.1023/A:1024033216561 View ArticleGoogle Scholar
- Han GY, Wang XF, Zhang GY, Ma ZY. Somatic embryogenesis and plant regeneration of recalcitrant cottons (Gossypium hirsutum). Afr J Biotechnol. 2009;8(3):432–7.Google Scholar
- Jiménez VM. Involvement of plant hormones and plant growth regulators on in vitro somatic embryogenesis. Plant Growth Regul. 2005;47(2):91–110. https://doi.org/10.1007/s10725-005-3478-x View ArticleGoogle Scholar
- Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. https://doi.org/10.1093/nar/28.1.27
- Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923 View ArticlePubMedPubMed CentralGoogle Scholar
- Leelavathi S, Sunnichan VG, Kumria R, et al. A simple and rapid Agrobacterium-mediated transformation protocol for cotton (Gossypium hirsutum L.): embryogenic calli as a source to generate large numbers of transgenic plants. Plant Cell Rep. 2004;22(7):465–70. https://doi.org/10.1007/s00299-003-0710-x View ArticlePubMedGoogle Scholar
- Leng CX, Li FG, Chen GY, Liu CL. cDNA-AFLP analysis of somatic embryogenesis at early stage in TM-1(Gossypium hirsutum L.). Acta Bot Boreal-Occident Sin. 2007;27(2):233–7.Google Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323. https://doi.org/10.1186/1471-2105-12-323 View ArticlePubMedPubMed CentralGoogle Scholar
- Li FG, Fan GY, Lu CR, et al. Genome sequence of cultivated upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat Biotechnol. 2015;33(5):524–30. https://doi.org/10.1038/nbt.3208 View ArticlePubMedGoogle Scholar
- Li FG, Li XL, Li FL. Changes of major biochemical metabolites in somatic embryogenesis of cotton. Acta Agric Univ Henanensis. 1994;28(3):313–6. https://doi.org/10.16445/j.cnki.1000-2340.1994.04.020 Google Scholar
- Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60. https://doi.org/10.1093/bioinformatics/btp324 View ArticlePubMedPubMed CentralGoogle Scholar
- Liu CL, Tian RP, Kong DP, et al. Establishment and application of efficient transformation system for cotton. Sci Agric Sin. 2014;47(21):4183–97. https://doi.org/10.3864/j.issn.0578-1752.2014.21.005 Google Scholar
- Liu CL, Wu ZX, Zhang CJ, et al. Study on large-scale and high efficient transformation system mediated by Agrobacterium tumefaciens of cotton. Acta Bot Boreal-Occident Sin. 2004;24(5):768–75.Google Scholar
- Murashige T, Skoog F. A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol Plantarum. 1962;15(3):473–97. https://doi.org/10.1111/j.1399-3054.1962.tb08052.x View ArticleGoogle Scholar
- Rajeswari S, Muthuramu S, Chandirakala R, et al. Callus induction, somatic embryogenesis and plant regeneration in cotton (Gossypium hirsutum L.). Electron J Plant Breed. 2010;1(4):1186–90.Google Scholar
- TIANGEN. RNAprep Pure Kit (For Plant) Handbook n.d. http://www.tiangen.com/asset/imsupload/up0390628001433129719.pdf. Accessed 5 Oct 2016.
- Saito K, Yonekura-Sakakibara K, Nakabayashi R, et al. The flavonoid biosynthetic pathway in Arabidopsis: structural and genetic diversity. Plant Physiol Biochem. 2013;72:21–34. https://doi.org/10.1016/j.plaphy.2013.02.001 View ArticlePubMedGoogle Scholar
- Saldanha AJ. Java Treeview—extensible visualization of microarray data. Bioinformatics. 2004;20(17):3246–8. https://doi.org/10.1093/bioinformatics/bth349 View ArticlePubMedGoogle Scholar
- Shang HH, Liu CL, Zhang CJ, et al. Progress in mechanisms of cotton somatic embryogenesis. Acta Bot Boreal-Occident Sin. 2009;29(3):637–42.Google Scholar
- Shi YH, Zhu SW, Mao XZ, et al. Transcriptome profiling, molecular biological, and physiological studies reveal a major role for ethylene in cotton fiber cell elongation. The Plant Cell Online. 2006;18(3):651–64. https://doi.org/10.1105/tpc.105.040303 View ArticleGoogle Scholar
- Sun JY, Li WM, Zhang HS, et al. Somatic embryogenesis and plant regeneration in glandless upland cotton (Gossypium hirsutum L.). Front Agric China. 2009;3(3):279–83. https://doi.org/10.1007/s11703-009-0060-8 View ArticleGoogle Scholar
- Sunilkumar G, Campbell LM, Puckhaber L, et al. Engineering cottonseed for use in human nutrition by tissue-specific reduction of toxic gossypol. Proc Natl Acad Sci. 2006;103(48):18054–9. https://doi.org/10.1073/pnas.0605389103 View ArticlePubMedGoogle Scholar
- Tarazona S, García F, Ferrer A, et al. NOIseq: a RNA-seq differential expression method robust for sequencing depth biases. EMBnetJournal. 2012;17(B):18–9. https://doi.org/10.14806/ej.17.B.265.Google Scholar
- Trolinder NL, Chen X. Genotype specificity of the somatic embryogenesis response in cotton. Plant Cell Rep. 1989;8(3):133–6. https://doi.org/10.1007/BF00716824 View ArticlePubMedGoogle Scholar
- Vogt T. Phenylpropanoid biosynthesis. Mol Plant. 2010;3(1):2–20. https://doi.org/10.1093/mp/ssp106 View ArticlePubMedGoogle Scholar
- Wu JH, Chen ZX, Li SJ, et al. Primary study on metabolic variation of callus in the procedure of cotton regeneration. J Biol. 1999;16(2):25–6.Google Scholar
- Wu JH, Zhang XL, Nie YC, et al. Factors affecting somatic embryogenesis and plant regeneration from a range of recalcitrant genotypes of Chinese cottons (Gossypium hirsutum L.). in vitro Cell Dev Biol Plant. 2004;40(4):371–5. https://doi.org/10.1079/IVP2004535
- Wu XM, Li FG, Zhang CJ, et al. Differential gene expression of cotton cultivar CCRI24 during somatic embryogenesis. J Plant Physiol. 2009;166(12):1275–83. https://doi.org/10.1016/j.jplph.2009.01.012 View ArticlePubMedGoogle Scholar
- Xie DY, Jin SX, Guo XP, Zhang XL. Somatic embryogenesis and plant regeneration in cotton cultivars from yellow and Yangtze river planting areas. Acta Agron Sin. 2007;33(3):394–400.Google Scholar
- Xu ZZ, Zhang CJ, Zhang XY, et al. Transcriptome profiling reveals auxin and cytokinin regulating somatic embryogenesis in different sister lines of cotton cultivar CCRI24. J Integr Plant Biol. 2013;55(7):631–42. https://doi.org/10.1111/jipb.12073 View ArticlePubMedGoogle Scholar
- Yang XY, Zhang XL, Yuan DJ, et al. Transcript profiling reveals complex auxin signalling pathway and transcription regulation involved in dedifferentiation and redifferentiation during somatic embryogenesis in cotton. BMC Plant Biol. 2012;12(1):110. https://doi.org/10.1186/1471-2229-12-110 View ArticlePubMedPubMed CentralGoogle Scholar
- Zeng FC, Zhang XL, Zhu LF, et al. Isolation and characterization of genes associated to cotton somatic embryogenesis by suppression subtractive hybridization and macroarray. Plant Mol Biol. 2006;60(2):167–83. https://doi.org/10.1007/s11103-005-3381-x View ArticlePubMedGoogle Scholar
- Zhang BH, Liu F, Yao CB. Plant regeneration via somatic embryogenesis in cotton. Plant Cell Tissue Organ Cult. 2000;60(2):89–94. https://doi.org/10.1023/A:1006488119200 View ArticleGoogle Scholar
- Zhang BH, Wang QL, Fang L, et al. Highly efficient plant regeneration through somatic embryogenesis in 20 elite commercial cotton (Gossypium hirsutum L.) cultivars. Plant OMICS. 2009;2(6):259–68.Google Scholar
- Zhang C, Yu S, Fan S, et al. Inheritance of somatic embryogenesis using leaf petioles as explants in upland cotton. Euphytica. 2011;181(1):55–63. https://doi.org/10.1007/s10681-011-0380-7 View ArticleGoogle Scholar
- Zhang XL, Sun JZ, Liu JL. A comparative study on the biochemical metabolites between the embryogenic and non-embryogenic calli from the variety “Coker 201” of Gossypium hirsutum L. Acta Agron Sin. 1992;18(3):176–82.Google Scholar
- Zhu HG, Tu LL, Jin SX, et al. Analysis of genes differentially expressed during initial cellular dedifferentiation in cotton. Chin Sci Bull. 2008;53(23):3666–76. https://doi.org/10.1007/s11434-008-0468-1 View ArticleGoogle Scholar