Function and Characteristic Analysis of Candidate PEAR Proteins in Populus yunnanensis.

PEAR proteins are a type of plant-specific DNA binding with one finger (Dof) transcription factors that play a key role in the regulation of plant growth, especially during phloem cell growth and seed germination in Arabidopsis. However, the identification, characteristics and function of PEAR proteins, particularly in woody plants, need to be further studied. In the present study, 43 candidate PEAR proteins harboring the conserved Zf-Dof domain were obtained in Populus yunnanensis. Based on phylogenetic and structural analysis, 10 representative PEAR candidates were selected, belonging to different phylogenetic groups. The functions of PEAR proteins in the stress response, signal transduction, and growth regulation of stem cambium and roots undergoing vigorous cell division in Arabidopsis were revealed based on their expression patterns as characterized by qRT-PCR analysis, in accordance with the results of cis-element analysis. In vitro experiments showed that the interaction of transcription factor (E2F) and cyclin indirectly reflects the growth regulation function of PEAR through light signaling and cell-cycle regulation. Therefore, our results provide new insight into the identity of PEAR proteins and their function in stress resistance and vigorous cell division regulation of tissues in P. yunnanensis, which may serve as a basis for further investigation of the functions and characteristics of PEAR proteins in other plants.

Genomic characterization of the NAC transcription factors, directed at understanding their functions involved in endocarp lignification of iron walnut (Juglans sigillata Dode).

The NAC (NAM, ATAF1/2, and CUC2) transcription factors (TF), one of the largest plant-specific gene families, play important roles in the regulation of plant growth and development, stress response and disease resistance. In particular, several NAC TFs have been identified as master regulators of secondary cell wall (SCW) biosynthesis. Iron walnut (Juglans sigillata Dode), an economically important nut and oilseed tree, has been widely planted in the southwest China. The thick and high lignified shell derived endocarp tissues, however, brings troubles in processing processes of products in industry. It is indispensable to dissect the molecular mechanism of thick endocarp formation for further genetic improvement of iron walnut. In the present study, based on genome reference of iron walnut, 117 NAC genes, in total, were identified and characterized in silico, which involves only computational analysis to provide insight into gene function and regulation. We found that the amino acids encoded by these NAC genes varied from 103 to 1,264 in length, and conserved motif numbers ranged from 2 to 10. The JsiNAC genes were unevenly distributed across the genome of 16 chromosomes, and 96 of these genes were identified as segmental duplication genes. Furthermore, 117 JsiNAC genes were divided into 14 subfamilies (A-N) according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana and common walnut (Juglans regia). Furthermore, tissue-specific expression pattern analysis demonstrated that a majority of NAC genes were constitutively expressed in five different tissues (bud, root, fruit, endocarp, and stem xylem), while a total of 19 genes were specifically expressed in endocarp, and most of them also showed high and specific expression levels in the middle and late stages during iron walnut endocarp development. Our result provided a new insight into the gene structure and function of JsiNACs in iron walnut, and identified key candidate JsiNAC genes involved in endocarp development, probably providing mechanistic insight into shell thickness formation across nut species.

Global Transcriptomic Analyses Provide New Insight into the Molecular Mechanisms of Endocarp Formation and Development in Iron Walnut (Juglans sigillata Dode).

Iron walnut (Juglans sigillata Dode) is a native species in southwestern China that exhibits variation in both fruit morphology and shell thickness. However, the underlying molecular processes controlling hardened endocarp development in walnut has not yet been reported. Here, we generated transcriptional profiles of iron walnut endocarp at three developmental stages using "Dapao", the most common commercial variety. Using pairwise comparisons between these three stages, a total of 8555 non-redundant differentially expressed genes (DEGs) were identified, and more than one-half of the total DEGs exhibited significant differential expression in stage I as compared with stage II or stage III, suggesting that the first stage may ultimately determine the final characteristics of the mature walnut shell. Furthermore, in the clustering analysis of the above DEGs, 3682, 2349, and 2388 genes exhibited the highest expression in stages I, II, and III, respectively. GO enrichment analysis demonstrated that the major transcriptional variation among the three developmental stages was caused by differences in cell growth, plant hormones, metabolic process, and phenylpropanoid metabolism. Namely, using the tissue-specific expression analysis and a gene co-expression network, we identified MADS-box transcription factor JsiFBP2 and bHLH transcription factor JsibHLH94 as candidate regulators of endocarp formation in the early stage, and JsiNAC56 and JsiMYB78 might play key roles in regulating the lignification process of endocarp in the late stage. This study provides useful information for further research to dissect the molecular mechanisms governing the shell formation and development of iron walnut.

Functional Genome Analyses Reveal the Molecular Basis of Oil Accumulation in Developing Seeds of Castor Beans.

Castor (Ricinus communis L.) seeds produce abundant ricinoleic acid during seed maturation, which is important for plant development and human demands. Ricinoleic acid, as a unique hydroxy fatty acid (HFA), possesses a distinct bond structure that could be used as a substitute for fossil fuels. Here, we identified all homologous genes related to glycolysis, hydroxy fatty acid biosynthesis, and triacylglycerol (TAG) accumulation in castor seeds. Furthermore, we investigated their expression patterns globally during five seed development stages. We characterized a total of 66 genes involved in the glycolysis pathway, with the majority exhibiting higher expression levels during the early stage of castor bean seed development. This metabolic process provided abundant acetyl-CoA for fatty acid (FA) biosynthesis. Subsequently, we identified 82 genes involved in the processes of de novo FA biosynthesis and TAG assembly, with the majority exhibiting high expression levels during the middle or late stages. In addition, we examined the expression patterns of the transcription factors involved in carbohydrate and oil metabolism. For instance, RcMYB73 and RcERF72 exhibited high expression levels during the early stage, whereas RcWRI1, RcABI3, and RcbZIP67 showed relatively higher expression levels during the middle and late stages, indicating their crucial roles in seed development and oil accumulation. Our study suggests that the high HFA production in castor seeds is attributed to the interaction of multiple genes from sugar transportation to lipid droplet packaging. Therefore, this research comprehensively characterizes all the genes related to glycolysis, fatty acid biosynthesis, and triacylglycerol (TAG) accumulation in the castor and provides novel insight into exploring the genetic mechanisms underlying seed oil accumulation in the endosperm of castor beans.

Genome-wide characterization and evolutionary analysis of linker histones in castor bean (Ricinus communis).

H1s, or linker histones, are ubiquitous proteins in eukaryotic cells, consisting of a globular GH1 domain flanked by two unstructured tails. Whilst it is known that numerous non-allelic variants exist within the same species, the degree of interspecific and intraspecific variation and divergence of linker histones remain unknown. The conserved basic binding sites in GH1 and evenly distributed strong positive charges on the C-terminal domain (CTD) are key structural characters for linker histones to bind chromatin. Based on these features, we identified five linker histones from 13 GH1-containing proteins in castor bean (Ricinus communis), which were named as RcH1.1, RcH1.2a, RcH1.2b, RcH1.3, and RcH1.4 based on their phylogenetic relationships with the H1s from five other economically important Euphorbiaceae species (Hevea brasiliensis Jatropha curcas, Manihot esculenta Mercurialis annua, and Vernicia fordii) and Arabidopsis thaliana. The expression profiles of RcH1 genes in a variety of tissues and stresses were determined from RNA-seq data. We found three RcH1 genes (RcH1.1, RcH1.2a, and RcH1.3) were broadly expressed in all tissues, suggesting a conserved role in stabilizing and organizing the nuclear DNA. RcH1.2a and RcH1.4 was preferentially expressed in floral tissues, indicating potential involvement in floral development in castor bean. Lack of non-coding region and no expression detected in any tissue tested suggest that RcH1.2b is a pseudogene. RcH1.3 was salt stress inducible, but not induced by cold, heat and drought in our investigation. Structural comparison confirmed that GH1 domain was highly evolutionarily conserved and revealed that N- and C-terminal domains of linker histones are divergent between variants, but highly conserved between species for a given variant. Although the number of H1 genes varies between species, the number of H1 variants is relatively conserved in more closely related species (such as within the same family). Through comparison of nucleotide diversity of linker histone genes and oil-related genes, we found similar mutation rate of these two groups of genes. Using Tajima's D and ML-HKA tests, we found RcH1.1 and RcH1.3 may be under balancing selection.

Separation and Purification of Two Saponins from Paris polyphylla var. yunnanensis by a Macroporous Resin.

An effective method for separating and purifying critical saponins (polyphyllin II and polyphyllin VII) from a Paris polyphylla var. yunnanensis extract was developed in this study which was environmentally friendly and economical. Static adsorption kinetics, thermodynamics, and the dynamic adsorption-desorption of macroporous resins were investigated, and then the conditions of purification and separation were optimized by fitting with an adsorption thermodynamics equation and a kinetic equation. Effective NKA-9 resin from seven macroporous resins was screened out to separate and purify the two saponins. The static adsorption and dynamic adsorption were chemical and physical adsorption dual-processes on the NKA-9 resin. Under the optimum parameters, the contents of polyphyllin II and polyphyllin VII in the product were 17.3-fold and 28.6-fold those in plant extracts, respectively. The total yields of the two saponins were 93.16%. This research thus provides a theoretical foundation for the large-scale industrial production of the natural drugs polyphyllin II and polyphyllin VII.

Function and Evolution of C1-2i Subclass of C2H2-Type Zinc Finger Transcription Factors in POPLAR.

C2H2 zinc finger (C2H2-ZF) transcription factors participate in various aspects of normal plant growth regulation and stress responses. C1-2i C2H2-ZFs are a special subclass of conserved proteins that contain two ZnF-C2H2 domains. Some C1-2i C2H2-ZFs in Arabidopsis (ZAT) are involved in stress resistance and other functions. However, there is limited information on C1-2i C2H2-ZFs in Populus trichocarpa (PtriZATs). To analyze the function and evolution of C1-2i C2H2-ZFs, eleven PtriZATs were identified in P. trichocarpa, which can be classified into two subgroups. The protein structure, conserved ZnF-C2H2 domains and QALGGH motifs, showed high conservation during the evolution of PtriZATs in P. trichocarpa. The spacing between two ZnF-C2H2 domains, chromosomal locations and cis-elements implied the original proteins and function of PtriZATs. Furthermore, the gene expression of different tissues and stress treatment showed the functional differentiation of PtriZATs subgroups and their stress response function. The analysis of C1-2i C2H2-ZFs in different Populus species and plants implied their evolution and differentiation, especially in terms of stress resistance. Cis-elements and expression pattern analysis of interaction proteins implied the function of PtriZATs through binding with stress-related genes, which are involved in gene regulation by via epigenetic modification through histone regulation, DNA methylation, ubiquitination, etc. Our results for the origin and evolution of PtriZATs will contribute to understanding the functional differentiation of C1-2i C2H2-ZFs in P. trichocarpa. The interaction and expression results will lay a foundation for the further functional investigation of their roles and biological processes in Populus.

Transcriptomic Analyses Shed Light on Critical Genes Associated with Bibenzyl Biosynthesis in Dendrobium officinale.

The Dendrobium plants (members of the Orchidaceae family) are used as traditional Chinese medicinal herbs. Bibenzyl, one of the active compounds in Dendrobium officinale, occurs in low amounts among different tissues. However, market demands require a higher content of thes compounds to meet the threshold for drug production. There is, therefore, an immediate need to dissect the physiological and molecular mechanisms underlying how bibenzyl compounds are biosynthesized in D. officinale tissues. In this study, the accumulation of erianin and gigantol in tissues were studied as representative compounds of bibenzyl. Exogenous application of Methyl-Jasmonate (MeJA) promotes the biosynthesis of bibenzyl compounds; therefore, transcriptomic analyses were conducted between D. officinale-treated root tissues and a control. Our results show that the root tissues contained the highest content of bibenzyl (erianin and gigantol). We identified 1342 differentially expressed genes (DEGs) with 912 up-regulated and 430 down-regulated genes in our transcriptome dataset. Most of the identified DEGs are functionally involved in the JA signaling pathway and the biosynthesis of secondary metabolites. We also identified two candidate cytochrome P450 genes and nine other enzymatic genes functionally involved in bibenzyl biosynthesis. Our study provides insights on the identification of critical genes associated with bibenzyl biosynthesis and accumulation in Dendrobium plants, paving the way for future research on dissecting the physiological and molecular mechanisms of bibenzyl synthesis in plants as well as guide genetic engineering for the improvement of Dendrobium varieties through increasing bibenzyl content for drug production and industrialization.

Genomic insights into the origin, domestication and genetic basis of agronomic traits of castor bean.

BACKGROUND: Castor bean (Ricinus communis L.) is an important oil crop, which belongs to the Euphorbiaceae family. The seed oil of castor bean is currently the only commercial source of ricinoleic acid that can be used for producing about 2000 industrial products. However, it remains largely unknown regarding the origin, domestication, and the genetic basis of key traits of castor bean. RESULTS: Here we perform a de novo chromosome-level genome assembly of the wild progenitor of castor bean. By resequencing and analyzing 505 worldwide accessions, we reveal that the accessions from East Africa are the extant wild progenitors of castor bean, and the domestication occurs ~ 3200 years ago. We demonstrate that significant genetic differentiation between wild populations in Kenya and Ethiopia is associated with past climate fluctuation in the Turkana depression ~ 7000 years ago. This dramatic change in climate may have caused the genetic bottleneck in wild castor bean populations. By a genome-wide association study, combined with quantitative trait locus analysis, we identify important candidate genes associated with plant architecture and seed size. CONCLUSIONS: This study provides novel insights of domestication and genome evolution of castor bean, which facilitates genomics-based breeding of this important oilseed crop and potentially other tree-like crops in future.

Bulked segregant analysis reveals candidate genes responsible for dwarf formation in woody oilseed crop castor bean.

Plant dwarfism is a desirable agronomic trait in non-timber trees, but little is known about the physiological and molecular mechanism underlying dwarfism in woody plants. Castor bean (Ricinus communis) is a typical woody oilseed crop. We performed cytological observations within xylem, phloem and cambia tissues, revealing that divergent cell growth in all tissues might play a role in the dwarf phenotype in cultivated castor bean. Based on bulked segregant analyses for a F2 population generated from the crossing of a tall and a dwarf accession, we identified two QTLs associated with plant height, covering 325 candidate genes. One of these, Rc5NG4-1 encoding a putative IAA transport protein localized in the tonoplast was functionally characterized. A non-synonymous SNP (altering the amino acid sequence from Y to C at position 218) differentiated the tall and dwarf plants and we confirmed, through heterologous yeast transformation, that the IAA uptake capacities of Rc5NG4-1Y and Rc5NG4-1C were significantly different. This study provides insights into the physiological and molecular mechanisms of dwarfing in woody non-timber economically important plants, with potential to aid in the genetic breeding of castor bean and other related crops.

Changes and Associations of Genomic Transcription and Histone Methylation with Salt Stress in Castor Bean.

Soil salinity is a major source of abiotic plant stress, adversely affecting plant growth, development and productivity. Although the physiological and molecular mechanisms that underlie plant responses to salt stress are becoming increasingly understood, epigenetic modifications, such as histone methylations and their potential regulation of the transcription of masked genes at the genome level in response to salt stress, remain largely unclear. Castor bean, an important nonedible oil crop, has evolved the capacity to grow under salt stress. Here, based on high-throughput RNA-seq and ChIP-seq data, we systematically investigated changes in genomic transcription and histone methylation using typical histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 tri-methylated lysine 27 (H3K27me3) markers in castor bean leaves subjected to salt stress. The results showed that gain or loss of histone methylation was closely associated with activated or repressed gene expression, though variations in both transcriptome and histone methylation modifications were relatively narrow in response to salt stress. Diverse salt responsive genes and switched histone methylation sites were identified in this study. In particular, we found for the first time that the transcription of the key salt-response regulator RADIALIS-LIKE SANT (RSM1), a MYB-related transcription factor involved in ABA(abscisic acid)-mediated salt stress signaling, was potentially regulated by bivalent H3K4me3-H3K27me3 modifications. Combining phenotypic variations with transcriptional and epigenetic changes, we provide a comprehensive profile for understanding histone modification, genomic transcription and their associations in response to salt stress in plants.

Comparative proteomic and transcriptomic analyses provide new insight into the formation of seed size in castor bean.

BACKGROUND: Little is known about the molecular basis of seed size formation in endospermic seed of dicotyledons. The seed of castor bean (Ricinus communis L.) is considered as a model system in seed biology studies because of its persistent endosperms throughout seed development. RESULTS: We compared the size of endosperm and endospermic cells between ZB107 and ZB306 and found that the larger seed size of ZB107 resulted from a higher cell count in the endosperm, which occupy a significant amount of the total seed volume. In addition, fresh weight, dry weight, and protein content of seeds were remarkably higher in ZB107 than in ZB306. Comparative proteomic and transcriptomic analyses were performed between large-seed ZB107 and small-seed ZB306, using isobaric tags for relative and absolute quantification (iTRAQ) and RNA-seq technologies, respectively. A total of 1416 protein species were identified, of which 173 were determined as differentially abundant protein species (DAPs). Additionally, there were 9545 differentially expressed genes (DEGs) between ZB306 and ZB107. Functional analyses revealed that these DAPs and DEGs were mainly involved in cell division and the metabolism of carbohydrates and proteins. CONCLUSIONS: These findings suggest that both cell number and storage-component accumulation are critical for the formation of seed size, providing new insight into the potential mechanisms behind seed size formation in endospermic seeds.

Application of a high-resolution genetic map for chromosome-scale genome assembly and fine QTLs mapping of seed size and weight traits in castor bean.

Castor bean (Ricinus communis L., Euphorbiaceae) is a critical biodiesel crop and its seed derivatives have important industrial applications. Due to lack of a high-density genetic map, the breeding and genetic improvement of castor bean has been largely restricted. In this study, based on a recombinant inbred line (RIL) population consisting of 200 individuals, we generated 8,896 high-quality genomic SNP markers and constructed a high-resolution genetic map with 10 linkage groups (LGs), spanning 1,852.33 centiMorgan (cM). Based on the genetic map, 996 scaffolds from the draft reference genome were anchored onto 10 pseudo-chromosomes, covering 84.43% of the castor bean genome. Furthermore, the quality of the pseudo-chromosome scale assembly genome was confirmed via genome collinearity analysis within the castor bean genome as well as between castor bean and cassava. Our results provide new evidence that the phylogenetic position of castor bean is relatively solitary from other taxa in the Euphorbiaceae family. Based on the genetic map, we identified 16 QTLs that control seed size and weight (covering 851 candidate genes). The findings will be helpful for further research into potential new mechanisms controlling seed size and weight in castor bean. The genetic map and improved pseudo-chromosome scale genome provide crucial foundations for marker-assisted selection (MAS) of QTL governing important agronomic traits, as well as the accelerated molecular breeding of castor bean in a cost-effective pattern.

Global Gene Expression of Seed Coat Tissues Reveals a Potential Mechanism of Regulating Seed Size Formation in Castor Bean.

The physiological and molecular basis of seed size formation is complex, and the development of seed coat (derived from integument cells) might be a critical factor that determines seed size formation for many endospermic seeds. Castor bean (Ricinus communis L.), a model system of studying seed biology, has large and persistent endosperm with a hard seed coat at maturity. Here, we investigated the potential molecular mechanisms underlying seed size formation in castor bean by comparing the difference between global gene expression within developing seed coat tissues between the large-seed ZB107 and small-seed ZB306. First, we observed the cell size of seed coat and concluded that the large seed coat area of ZB107 resulted from more cell numbers (rather than cell size). Furthermore, we found that the lignin proportion of seed coat was higher in ZB306. An investigation into global gene expression of developing seed coat tissues revealed that 815 genes were up-regulated and 813 were down-regulated in ZB306 relative to ZB107. Interestingly, we found that many genes involved in regulating cell division were up-regulated in ZB107, whereas many genes involved in regulating lignin biosynthesis (including several NAC members, as well as MYB46/83 and MYB58/63) and in mediating programmed cell death (such as CysEP1 and ßVPE) were up-regulated in ZB306. Furthermore, the expression patterns of the genes mentioned above indicated that the lignification of seed coat tissues was enhanced and occurred earlier in the developing seeds of ZB306. Taken together, we tentatively proposed a potential scenario for explaining the molecular mechanisms of seed coat governing seed size formation in castor bean by increasing the cell number and delaying the onset of lignification in seed coat tissues in large-seed ZB107. This study not only presents new information for possible modulation of seed coat related genes to improve castor seed yield, but also provides new insights into understanding the molecular basis of seed size formation in endospermic seeds with hard seed coat.

Genetic Diversity Analysis Reveals Genetic Differentiation and Strong Population Structure in Calotropis Plants.

The genus Calotropis (Asclepiadaceae) is comprised of two species, C. gigantea and C. procera, which both show significant economic potential for use of their seed fibers in the textile industry, and of their bioactive compounds as new medicinal resources. The available wild-sourced germplasm contains limited genetic information that restricts further germplasm exploration for the purposes of domestication. We here developed twenty novel EST-SSR markers and applied them to assess genetic diversity, population structure and differentiation within Calotropis. The polymorphic information index of these markers ranged from 0.102 to 0.800; indicating that they are highly informative. Moderate genetic diversity was revealed in both species, with no difference between species in the amount of genetic diversity. Population structure analysis suggested five main genetic groups (K = 5) and relatively high genetic differentiation (FST = 0.528) between the two species. Mantel test analysis showed strong correlation between geographical and genetic distance in C. procera (r = 0.875, p = 0.020) while C. gigantea showed no such correlation (r = 0.390, p = 0.210). This study provides novel insights into the genetic diversity and population structure of Calotropis, which will promote further resource utilization and the development of genetic improvement strategies for Calotropis.

Sequence analysis of the Hsp70 family in moss and evaluation of their functions in abiotic stress responses.

The 70-kD heat shock proteins (Hsp70s) are highly conserved molecular chaperones that play essential roles in cellular processes including abiotic stress responses. Physcomitrella patens serves as a representative of the first terrestrial plants and can recover from serious dehydration. To assess the possible relationship between P. patens Hsp70s and dehydration tolerance, we analyzed the P. patens genome and found at least 21 genes encoding Hsp70s. Gene structure and motif composition were relatively conserved in each subfamily. The intron-exon structure of PpcpHsp70-2 was different from that of other PpcpHsp70s; this gene exhibits several forms of intron retention, indicating that introns may play important roles in regulating gene expression. We observed expansion of Hsp70s in P. patens, which may reflect adaptations related to development and dehydration tolerance, and results mainly from tandem and segmental duplications. Expression profiles of rice, Arabidopsis and P. patens Hsp70 genes revealed that more than half of the Hsp70 genes were responsive to ABA, salt and drought. The presence of overrepresented cis-elements (DOFCOREZM and GCCCORE) among stress-responsive Hsp70s suggests that they share a common regulatory pathway. Moss plants overexpressing PpcpHsp70-2 showed salt and dehydration tolerance, further supporting a role in adaptation to land. This work highlights directions for future functional analyses of Hsp70s.

Roles of Hsp70s in Stress Responses of Microorganisms, Plants, and Animals.

Hsp70s (heat shock protein 70s) are a class of molecular chaperones that are highly conserved and ubiquitous in organisms ranging from microorganisms to plants and humans. Most research on Hsp70s has focused on the mechanisms of their functions as molecular chaperones, but recently, studies on stress responses are coming to the forefront. Hsp70s play key roles in cellular development and protecting living organisms from environmental stresses such as heat, drought, salinity, acidity, and cold. Moreover, functions of human Hsp70s are related to diseases including neurological disorders, cancer, and virus infection. In this review, we provide an overview of the specific roles of Hsp70s in response to stress, particularly abiotic stress, in all living organisms.