The somatic embryogenesis (SE) process of plants, as one of the typical responses to abiotic stresses with hormone, occurs through the dynamic expression of different proteins that constitute a complex regulatory network in biological activities and promotes plant totipotency. Plant SE includes two critical stages: primary embryogenic calli redifferentiation and somatic embryos development initiation, which leads to totipotency. The isobaric labels tandem mass tags (TMT) large-scale and quantitative proteomics technique was used to identify the dynamic protein expression changes in nonembryogenic calli (NEC), primary embryogenic calli (PEC) and globular embryos (GEs) of cotton. A total of 9369 proteins (6730 quantified) were identified; 805, 295 and 1242 differentially accumulated proteins (DAPs) were identified in PEC versus NEC, GEs versus PEC and GEs versus NEC, respectively. Eight hundred and five differentially abundant proteins were identified, 309 of which were upregulated and 496 down regulated in PEC compared with NEC. Of the 295 DAPs identified between GEs and PEC, 174 and 121 proteins were up- and down regulated, respectively. Of 1242 differentially abundant proteins, 584 and 658 proteins were up- and down regulated, respectively, in GEs versus NEC. We have also complemented the authenticity and accuracy of the proteomic analysis. Systematic analysis indicated that peroxidase, photosynthesis, environment stresses response processes, nitrogen metabolism, phytohormone response/signal transduction, transcription/posttranscription and modification were involved in somatic embryogenesis. The results generated in this study demonstrate a proteomic molecular basis and provide a valuable foundation for further investigation of the roles of DAPs in the process of SE transdifferentiation during cotton totipotency.

As a notable illustration of totipotency and plant regeneration, somatic embryogenesis (SE) is the developmental reprogramming of somatic cells toward the embryogenesis pathway, the key step for genetic engineering. Investigations examining the totipotency process are of great fundamental and practical importance in crop biotechnology. However, high-frequency regeneration of cotton via SE has been limited due to genotype-dependent response. The molecular basis deciphering SE genotype recalcitrance remains largely unexplored in cotton. In the current study, to comprehensively investigate the dynamic transcriptional profiling and gene regulatory patterns involved in SE process, a genome-wide RNA sequencing analysis was performed in two cotton genotypes with distinct embryogenic abilities, the highly embryogenic genotype Yuzao 1 (YZ) and the recalcitrant genotype Lumian 1 (LM). Three typical developmental staged cultures of early SE—hypocotyls (HY), nonembryogenic calli (NEC) and primary embryogenic calli (PEC)—were selected to establish the transcriptional profiles. Our data revealed that a total of 62,562 transcripts were present amongst different developmental stages in the two genotypes. Of these, 18,394 and 26,514 differentially expressed genes (DEGs) were identified during callus dedifferentiation (NEC-VS-HY) and embryogenic transdifferentiation (PEC-VS-NEC), respectively in the recalcitrant genotype, 21,842 and 22,343 DEGs in the highly embryogenic genotype. Furthermore, DEGs were clustered into six expression patterns during cotton SE process in the two genotypes. Moreover, functional enrichment analysis revealed that DEGs were significantly enriched in fatty acid, tryptophan and pyruvate metabolism in the highly embryogenic genotype and in DNA conformation change otherwise in the recalcitrant genotype. In addition, critical SE-associated expressed transcription factors, as well as alternative splicing events, were notably and preferentially activated during embryogenic transdifferentiation in the highly embryogenic genotype compared with the recalcitrant genotype. Taken together, by systematically comparing two genotypes with distinct embryogenic abilities, the findings in our study revealed a comprehensive overview of the dynamic gene regulatory patterns and uncharacterized complex regulatory pathways during cotton SE genotype-dependent response. Our work provides insights into the molecular basis and important gene resources for understanding the underlying genotype recalcitrance during SE process and plant regeneration, thereby holding great promise for accelerating the application of biotechnology to cotton for improving its breeding efficiency.

Intercellular material transport and information transmission in plants are carried out through the plasmodesmata (PD). The amount of callose around the PD controls channel permeability. In plants, β-1,3-glucanase can degrade callose and affect plant growth and development. In this study, the gene producing PD-localized β-1,3-glucanase and regulating the leaf trichomes is identified and named PdBG4. Based on functional analysis through a series of genetic manipulation assays, we found that the high expression of PdBG4 was associated with strong PD permeability and short Arabidopsis thaliana leaf trichomes. Conversely, the low expression of PdBG4 correlated with weak PD permeability and long Arabidopsis thaliana leaf trichomes. This study revealed that the PdBG4 gene negatively modulates leaf trichome growth and development by regulating PD permeability.

Auxin has a profound impact on plant physiology and participates in almost all aspects of plant development processes. Auxin exerts profound pleiotropic effects on plant growth and differentiation by regulating the auxin response genes’ expressions. The classical auxin reaction is usually mediated by auxin response factors (ARFs), which bind to the auxin response element (AuxRE) in the promoter region of the target gene. Experiments have generated only a limited number of plant genes with well-characterized functions. It is still unknown how many genes respond to exogenous auxin treatment. An economical and effective method was proposed for the genome-wide discovery of genes responsive to auxin in a model plant, Arabidopsis thaliana (A. thaliana). Our method relies on cis-regulatory-element-based targeted gene finding across different promoters in a genome. We first exploit and analyze auxin-specific cis-regulatory elements for the transcription of the target genes, and then identify putative auxin responsive genes whose promoters contain the elements in the collection of over 25,800 promoters in the A. thaliana genome. Evaluating our result by comparing with a published database and the literature, we found that this method has an accuracy rate of 65.2% (309/474) for predicting candidate genes responsive to auxin. Chromosome distribution and annotation of the putative auxin-responsive genes predicted here were also mined. The results can markedly decrease the number of identified but merely potential auxin target genes and also provide useful clues for improving the annotation of gene that lack functional information.