Genome-wide characterization of the cellulose synthase gene superfamily in Pyrus bretschneideri and reveal its potential role in stone cell formation

Genome-wide analysis of the cellulose toolbox of Primulina eburnea, a calcium-rich vegetable

BACKGROUND

Human-guided crop domestication has lasted for more than 10,000 years. In terms of the domestication and breeding of vegetables, cellulose content in edible tissues is one of the most important traits. Primulina eburnea is a recently developed calcium-rich vegetable with a high soluble and bioavailable calcium content in its leaves. However, the high cellulose content in the leaves hampers the taste, and no research has been reported on the genetic basis of cellulose biosynthesis in this calcium-rich vegetable.

RESULTS

We identified 36 cellulose biosynthesis-involved genes belonging to eight gene families in the P. eburnea genome. The cellulose accumulated decreasingly throughout leaf development. Nineteen genes were considered core genes in cellulose biosynthesis, which were highly expressed in buds but lowly expressed in mature leaves. In the nitrogen fertilization experiment, exogenous nitrogen decreased the cellulose content in the buds. The expressing pattern of 14 genes were consistent with phenotypic variation in the nitrogen fertilization experiment, and thus they were proposed as cellulose toolbox genes.

CONCLUSIONS

The present study provides a strong basis for the subsequent functional research of cellulose biosynthesis-involved genes in P. eburnea, and provides a reference for breeding and/or engineering this calcium-rich vegetable with decreased leaf cellulose content to improve the taste.

Genome-wide characterization of the cellulose synthase gene family in Ziziphus jujuba reveals its function in cellulose biosynthesis during fruit development

The cellulose synthase (Ces/Csl) is a key enzyme in plant cellulose synthesis. Jujube fruits are rich in cellulose. 29 ZjCesA/Csl genes were identified in jujube genome and showed tissue-specific expression. 13 genes highly expressed in jujube fruit exhibited obviously sequential expressions during the fruit development, indicating that they might play distinct roles during the process. Meanwhile, the correlation analysis showed the expressions of ZjCesA1 and ZjCslA1 were significant positive related to the cellulose synthase activities. Furthermore, transient overexpressions of ZjCesA1 or ZjCslA1 in jujube fruits significantly increased cellulose synthase activities and contents, whereas silencing of ZjCesA1 or ZjCslA1 in jujube seedlings obviously reduced cellulose levels. Moreover, the Y2H assays verified that ZjCesA1 and ZjCslA1 may participate in cellulose synthesis by forming protein complexes. The study not only reveals the bioinformatics characteristics and functions of cellulose synthase genes in jujube, but also provides clues for studying cellulose synthesis in other fruits.

Comparative genomic analysis of the COBRA genes in six Rosaceae species and expression analysis in Chinese white pear (Pyrus bretschneideri)

COBRA-Like (COBL) genes encode a glycosylphosphatidylinositol (GPI) anchoring protein unique to plants. In current study, 87 COBRA genes were identified in 6 Rosaceae species, including Pyrus bretschneideri (16 genes), Malus domestica (22 genes), Fragaria vesca (13 genes), Prunus mume (11 genes), Rubus occidentalis (13 genes) and Prunus avium (12 genes). We revealed the evolution of the COBRA gene in six Rosaceae species by phylogeny, gene structure, conservative sequence, hydrophobicity analysis, gene replication events and sliding window analysis. In addition, based on the analysis of expression patterns in pear fruit combined with bioinformatics, we identified PbCOBL12 and PbCOBL13 as potential genes regulating secondary cell wall (SCW) formation during pear stone cell development. This study aimed to understand the evolutionary relationship of the COBRA gene in Rosaceae species, clarify the potential function of COBRA in pear fruit development, and provide essential theoretical basis and gene resources for improving pear fruit quality through genetical modification mechanism.

Species-Specific Gene Expansion of the Cellulose synthase Gene Superfamily in the Orchidaceae Family and Functional Divergence of Mannan Synthesis-Related Genes in Dendrobium officinale

Plant Cellulose synthase genes constitute a supergene family that includes the Cellulose synthase (CesA) family and nine Cellulose synthase-like (Csl) families, the members of which are widely involved in the biosynthesis of cellulose and hemicellulose. However, little is known about the Cellulose synthase superfamily in the family Orchidaceae, one of the largest families of angiosperms. In the present study, we identified and systematically analyzed the CesA/Csl family members in three fully sequenced Orchidaceae species, i.e., Dendrobium officinale, Phalaenopsis equestris, and Apostasia shenzhenica. A total of 125 Cellulose synthase superfamily genes were identified in the three orchid species and classified into one CesA family and six Csl families: CslA, CslC, CslD, CslE, CslG, and CslH according to phylogenetic analysis involving nine representative plant species. We found species-specific expansion of certain gene families, such as the CslAs in D. officinale (19 members). The CesA/Csl families exhibited sequence divergence and conservation in terms of gene structure, phylogeny, and deduced protein sequence, indicating multiple origins via different evolutionary processes. The distribution of the DofCesA/DofCsl genes was investigated, and 14 tandemly duplicated genes were detected, implying that the expansion of DofCesA/DofCsl genes may have originated via gene duplication. Furthermore, the expression profiles of the DofCesA/DofCsl genes were investigated using transcriptome sequencing and quantitative Real-time PCR (qRT-PCR) analysis, which revealed functional divergence in different tissues and during different developmental stages of D. officinale. Three DofCesAs were highly expressed in the flower, whereas DofCslD and DofCslC family genes exhibited low expression levels in all tissues and at all developmental stages. The 19 DofCslAs were differentially expressed in the D. officinale stems at different developmental stages, among which six DofCslAs were expressed at low levels or not at all. Notably, two DofCslAs (DofCslA14 and DofCslA15) showed significantly high expression in the stems of D. officinale, indicating a vital role in mannan synthesis. These results indicate the functional redundancy and specialization of DofCslAs with respect to polysaccharide accumulation. In conclusion, our results provide insights into the evolution, structure, and expression patterns of CesA/Csl genes and provide a foundation for further gene functional analysis in Orchidaceae and other plant species.

MYB transcription factor family in sweet cherry (Prunus avium L.): genome-wide investigation, evolution, structure, characterization and expression patterns

BACK GROUND

MYB Transcription factors (TFs) are most imperative and largest gene family in plants, which participate in development, metabolism, defense, differentiation and stress response. The MYB TFs has been studied in various plant species. However, comprehensive studies of MYB gene family in the sweet cherry (Prunus avium L.) are still unknown.

RESULTS

In the current study, a total of 69 MYB genes were investigated from sweet cherry genome and classified into 28 subfamilies (C1-C28 based on phylogenetic and structural analysis). Microcollinearity analysis revealed that dispersed duplication (DSD) events might play an important role in the MYB genes family expansion. Chromosomal localization, the synonymous (Ks) and nonsynonymous (Ka) analysis, molecular characteristics (pI, weight and length of amino acids) and subcellular localization were accomplished using several bioinformatics tools. Furthermore, the members of distinct subfamilies have diverse cis-acting regions, conserved motifs, and intron-exon architectures, indicating functional heterogeneity in the MYB family. Moreover, the transcriptomic data exposed that MYB genes might play vital role in bud dormancy. The quantitative real-time qRT-PCR was carried out and the expression pattern indicated that MYB genes significantly expressed in floral bud as compared to flower and fruit.

CONCLUSION

Our comprehensive findings provide supportive insights into the evolutions, expansion complexity and functionality of PavMYB genes. These PavMYB genes should be further investigated as they seem to be brilliant candidates for dormancy manipulation in sweet cherry.

Integrative Analysis of the GRAS Genes From Chinese White Pear (Pyrus bretschneideri): A Critical Role in Leaf Regeneration

GRAS is a transcription regulator factor, which plays an important role in plant growth and development. Previous analyses found that several GRAS functions have been identified, such as axillary bud meristem formation, radial root elongation, gibberellin signaling, light signaling, and abiotic stress. The GRAS family has been comprehensively evaluated in several species. However, little finding is on the GRAS transcription factors (TFs) in Chinese white pear. In this study, 99 PbGRAS were systemically characterized and renamed PbGRAS1 to PbGRAS99 according to their chromosomal localizations. Phylogenetic analysis and structural features revealed that could be classified into eight subfamilies (LISCL, Ls, SHR, HAM, SCL, PAT, SCR, and DELLA). Further analysis of introns/exons and conserved motifs revealed that they are diverse and functionally differentiated in number and structure. Synteny analysis among Pyrus bretschenedri, Prunus mume, Prunus avium, Fragaria vesca, and Prunus persica showed that GRAS duplicated regions were more conserved. Dispersed duplication events are the most common mechanism and may play a crucial role in the expansion of the GRAS gene family. In addition, cis-acting elements of the PbGRAS gene were found in promoter regions associated with hormone and environmental stress responses. Notably, the expression pattern detected by qRT-PCR indicated that PbGRAS genes were differentially expressed under gibberellin (GA), abscisic acid (ABA), and auxin (IAA) conditions, which are responsive to abiotic stress. PbGRAS89 and PbGRAS99 were highly expressed at different stages of hormone treatment and may play important role in leaf development. Therefore, we selected PbGRAS89 and PbGRAS99 to clone and construct pCAMBIA1301-PbGRAS89, 99 and transferred them into Arabidopsis thaliana. Finally, we observed and compared the changes of overexpressed plants and wild-type plants during regeneration. This method was used to analyze their roles in leaf regeneration of Chinese white pear. In addition, we also constructed pCAMBIA1305-PbGRAS89, 99, and transferred them into onion cells to determine the subcellular localization. Subcellular localization experiments showed that PbGRAS89 and PbGRAS99 were localized in the nucleus. In summary, the results of this study indicate that PbGRAS89 and PbGRAS99 are mainly responsible for leaf regeneration of Chinese white pear, which plays a positive role in callus formation and provides rich resources for studying GRAS gene functions.

Identification and Comprehensive Genome-Wide Analysis of Glutathione S-Transferase Gene Family in Sweet Cherry (Prunus avium) and Their Expression Profiling Reveals a Likely Role in Anthocyanin Accumulation

Glutathione S-transferases (GSTs) in plants are multipurpose enzymes that are involved in growth and development and anthocyanins transportation. However, members of the GST gene family were not identified in sweet cherry (Prunus avium). To identify the GST genes in sweet cherry, a genome-wide analysis was conducted. In this study, we identified 67 GST genes in P. avium genome and nomenclature according to chromosomal distribution. Phylogenetic tree analysis revealed that PavGST genes were classified into seven chief subfamily: TCHQD, Theta, Phi, Zeta, Lambda, DHAR, and Tau. The majority of the PavGST genes had a relatively well-maintained exon-intron and motif arrangement within the same group, according to gene structure and motif analyses. Gene structure (introns-exons) and conserved motif analysis revealed that the majority of the PavGST genes showed a relatively well-maintained motif and exons-introns configuration within the same group. The chromosomal localization, GO enrichment annotation, subcellular localization, syntenic relationship, Ka/Ks analysis, and molecular characteristics were accomplished using various bioinformatics tools. Mode of gene duplication showed that dispersed duplication might play a key role in the expansion of PavGST gene family. Promoter regions of PavGST genes contain numerous cis-regulatory components, which are involved in multiple stress responses, such as abiotic stress and phytohormones responsive factors. Furthermore, the expression profile of sweet cherry PavGSTs showed significant results under LED treatment. Our findings provide the groundwork for future research into induced LED anthocyanin and antioxidants deposition in sweet cherries.

Comprehensive Comparative Analysis of the GATA Transcription Factors in Four Rosaceae Species and Phytohormonal Response in Chinese Pear (Pyrus bretschneideri) Fruit

The GATA gene family is one of the most important transcription factors (TFs). It extensively exists in plants, contributes to diverse biological processes such as the development process, and responds to environmental stress. Although the GATA gene family has been comprehensively and systematically studied in many species, less is known about GATA genes in Chinese pears (Pyrus bretschneideri). In the current study, the GATA gene family in the four Rosaceae genomes was identified, its structural characteristics identified, and a comparative analysis of its properties was carried out. Ninety-two encoded GATA proteins were authenticated in the four Rosaceae genomes (Pyrus bretschneideri, Prunus avium, Prunus mume, and Prunus persica) and categorized into four subfamilies (Ⅰ-Ⅳ) according to phylogeny. The majority of GATA genes contained one to two introns and conserved motif composition analysis revealed their functional divergence. Whole-genome duplications (WGDs) and dispersed duplication (DSD) played a key role in the expansion of the GATA gene family. The microarray indicated that, among P. bretschneideri, P. avium, P. mume and P. persica, GATA duplicated regions were more conserved between Pyrus bretschneideri and Prunus persica with 32 orthologous genes pairs. The physicochemical parameters, duplication patterns, non-synonymous (ka), and synonymous mutation rate (ks) and GO annotation ontology were performed using different bioinformatics tools. cis-elements respond to various phytohormones, abiotic/biotic stress, and light-responsive were found in the promoter regions of GATA genes which were induced via stimuli. Furthermore, subcellular localization of the PbGATA22 gene product was investigated, showing that it was present in the nucleus of tobacco (Nicotiana tabacum) epidermal cells. Finally, in silico analysis was performed on various organs (bud, leaf, stem, ovary, petal, and sepal) and different developmental stages of fruit. Subsequently, the expression profiles of PbGATA genes were extensively expressed under exogenous hormonal treatments of SA (salicylic acid), MeJA (methyl jasmonate), and ABA (abscisic acid) indicating that play important role in hormone signaling pathways. A comprehensive analysis of GATA transcription factors was performed through systematic biological approaches and comparative genomics to establish a theoretical base for further structural and functional investigations in Rosaceae species.

The cellulose synthase (CesA) gene family in four Gossypium species: phylogenetics, sequence variation and gene expression in relation to fiber quality in Upland cotton

Cellulose synthases (CesAs) are multi-subunit enzymes found on the plasma membrane of plant cells and play a pivotal role in cellulose production. The cotton fiber is mainly composed of cellulose, and the genetic relationships between CesA genes and cotton fiber yield and quality are not fully understood. Through a phylogenetic analysis, the CesA gene family in diploid Gossypium arboreum and Gossypium raimondii, as well as tetraploid Gossypium hirsutum ('TM-1') and Gossypium barbadense ('Hai-7124' and '3-79'), was divided into 6 groups and 15 sub-groups, with each group containing two to five homologous genes. Most CesA genes in the four species are highly collinear. Among the five cotton genomes, 440 and 1929 single nucleotide polymorphisms (SNPs) in the CesA gene family were identified in exons and introns, respectively, including 174 SNPs resulting in amino acid changes. In total, 484 homeologous SNPs between the A and D genomes were identified in diploids, while 142 SNPs were detected between the two tetraploids, with 32 and 82 SNPs existing within G. hirsutum and G. barbadense, respectively. Additionally, 74 quantitative trait loci near 18 GhCesA genes were associated with fiber quality. One to four GhCesA genes were differentially expressed (DE) in ovules at 0 and 3 days post anthesis (DPA) between two backcross inbred lines having different fiber lengths, but no DE genes were identified between these lines in developing fibers at 10 DPA. Twenty-seven SNPs in above DE CesA genes were detected among seven cotton lines, including one SNP in Ghi_A08G03061 that was detected in four G. hirsutum genotypes. This study provides the first comprehensive characterization of the cotton CesA gene family, which may play important roles in determining cotton fiber quality.

Genome-Wide Identification of Banana Csl Gene Family and Their Different Responses to Low Temperature between Chilling-Sensitive and Tolerant Cultivars

The cell wall plays an important role in responses to various stresses. The cellulose synthase-like gene (Csl) family has been reported to be involved in the biosynthesis of the hemicellulose backbone. However, little information is available on their involvement in plant tolerance to low-temperature (LT) stress. In this study, a total of 42 Csls were identified in Musa acuminata and clustered into six subfamilies (CslA, CslC, CslD, CslE, CslG, and CslH) according to phylogenetic relationships. The genomic features of MaCsl genes were characterized to identify gene structures, conserved motifs and the distribution among chromosomes. A phylogenetic tree was constructed to show the diversity in these genes. Different changes in hemicellulose content between chilling-tolerant and chilling-sensitive banana cultivars under LT were observed, suggesting that certain types of hemicellulose are involved in LT stress tolerance in banana. Thus, the expression patterns of MaCsl genes in both cultivars after LT treatment were investigated by RNA sequencing (RNA-Seq) technique followed by quantitative real-time PCR (qPCR) validation. The results indicated that MaCslA4/12, MaCslD4 and MaCslE2 are promising candidates determining the chilling tolerance of banana. Our results provide the first genome-wide characterization of the MaCsls in banana, and open the door for further functional studies.