Genome and transcriptome-wide analyses of cellulose synthase gene superfamily in soybean

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.

Advances in genomics for diversity studies and trait improvement in temperate fruit and nut crops under changing climatic scenarios

Genetic improvement of temperate fruit and nut crops through conventional breeding methods is not sufficient alone due to its extreme time-consuming, cost-intensive, and hard-to-handle approach. Again, few other constraints that are associated with these species, viz., their long juvenile period, high heterozygosity, sterility, presence of sexual incompatibility, polyploidy, etc., make their selection and improvement process more complicated. Therefore, to promote precise and accurate selection of plants based on their genotypes, supplement of advanced biotechnological tools, viz., molecular marker approaches along with traditional breeding methods, is highly required in these species. Different markers, especially the molecular ones, enable direct selection of genomic regions governing the trait of interest such as high quality, yield, and resistance to abiotic and biotic stresses instead of the trait itself, thus saving the overall time and space and helping screen fruit quality and other related desired traits at early stages. The availability of molecular markers like SNP (single-nucleotide polymorphism), DArT (Diversity Arrays Technology) markers, and dense molecular genetic maps in crop plants, including fruit and nut crops, led to a revelation of facts from genetic markers, thus assisting in precise line selection. This review highlighted several aspects of the molecular marker approach that opens up tremendous possibilities to reveal valuable information about genetic diversity and phylogeny to boost the efficacy of selection in temperate fruit crops through genome sequencing and thus cultivar improvement with respect to adaptability and biotic and abiotic stress resistance in temperate fruit and nut species.

Identification and bioinformatic analysis of the CaCesA/Csls family members and the expression of the CaCslD1 in the flower buds of CMS/Rf system in pepper

The cellulose synthase gene superfamily contains cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, which synthesize cellulose and hemicellulose in plant cell walls and play a crucial role in plant growth and development. However, the CesA/Csl gene family has not been reported in pepper. Therefore, the genome-wide research of the CaCesA/CaCsl gene family was conducted in pepper. In this study, a total of 39 CaCesA/CaCsls genes (10 CesAs genes and 29 Csls genes) were identified in pepper and unevenly distributed on 11 chromosomes. These CaCesA/Csls were divided into seven subfamilies (CesAs, CslAs, CslBs, CslCs, CslDs, CslEs, CslGs), and most of CaCesA/Csls genes are closely related to AtCesA/Csls genes. The cis-acting elements in the promoters of CaCesA/Csls genes are mainly related to hormone response and stress response. There are ten collinear gene pairs between the CesA/Csls gene family of pepper and Arabidopsis, and four fragment duplication gene pairs of the CaCesA/Csls genes were discovered. RNA-seq analysis shows that the majority of CaCesA/Csls are expressed in a variety of plant tissues, indicating that most CaCesA/Csls gene expression patterns are not organ-specific, and CaCslD1/D4 have the highest expression in anthers, followed by petal, ovary, and F9. RNA-seq analysis shows that most CaCesA/Csls are responsive to five hormones (IAA, GA₃, ABA, SA, and MeJA). The tissue-specific expression analysis of the CaCslD1 gene shows that the CaCslD1 gene is expressed specifically in flowers. In the flower buds IV of cytoplasmic male sterility (CMS) and its restoration of fertility (Rf) system, CaCslD1 reach the highest expression respectively. However, the relative expression level of CaCslD1 in the fertile accessions is extremely significantly higher than in the sterile accessions. This study shows an overall understanding of the CaCesA/Csls gene family and provides a new insight for understanding the function of CaCslD1 in pollen development and exploring the fertility restoration of CMS in pepper.

Genome wide association study to detect genetic regions related to isoflavone content in a mutant soybean population derived from radiation breeding

Isoflavones are major secondary metabolites that are exclusively produced by legumes, including soybean. Soy isoflavones play important roles in human health as well as in the plant defense system. The isoflavone content is influenced by minor-effect quantitative trait loci, which interact with polygenetic and environmental factors. It has been difficult to clarify the regulation of isoflavone biosynthesis because of its complex heritability and the influence of external factors. Here, using a genotype-by-sequencing-based genome-wide association mapping study, 189 mutant soybean genotypes (the mutant diversity pool, MDP) were genotyped on the basis of 25,646 high-quality single nucleotide polymorphisms (SNPs) with minor allele frequency of >0.01 except for missing data. All the accessions were phenotyped by determining the contents of 12 isoflavones in the soybean seeds in two consecutive years (2020 and 2021). Then, quantitative trait nucleotides (QTNs) related to isoflavone contents were identified and validated using multi-locus GWAS models. A total of 112 and 46 QTNs related to isoflavone contents were detected by multiple MLM-based models in 2020 and 2021, respectively. Of these, 12 and 5 QTNs were related to more than two types of isoflavones in 2020 and 2021, respectively. Forty-four QTNs were detected within the 441-Kb physical interval surrounding Gm05:38940662. Of them, four QTNs (Gm05:38936166, Gm05:38936167, Gm05:38940662, and Gm05:38940717) were located at Glyma.05g206900 and Glyma.05g207000, which encode glutathione S-transferase THETA 1 (GmGSTT1), as determined from previous quantitative trait loci annotations and the literature. We detected substantial differences in the transcript levels of GmGSTT1 and two other core genes (IFS1 and IFS2) in the isoflavone biosynthetic pathway between the original cultivar and its mutant. The results of this study provide new information about the factors affecting isoflavone contents in soybean seeds and will be useful for breeding soybean lines with high and stable concentrations of isoflavones.

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.

Investigating the Mechanism of Unilateral Cross Incompatibility in Longan (Dimocarpus longan Lour.) Cultivars (Yiduo × Shixia)

Longan (Dimocarpus longan Lour.) is an important subtropical fruit tree in China. Nearly 90% of longan fruit imports from Thailand are from the cultivar Yiduo. However, we have observed that there exists a unilateral cross incompatibility (UCI) when Yiduo is used as a female parent and Shixia (a famous Chinese cultivar) as a male parent. Here, we performed a comparative transcriptome analysis coupled with microscopy of pistils from two reciprocal pollination combinations [Shixia♂ × Yiduo♀(SY) and Yiduo♀ × Shixia♂(YS)] 4, 8, 12, and 24 h after pollination. We also explored endogenous jasmonic acid (JA) and jasmonyl isoleucine (JA-Ile) levels in pistils of the crosses. The microscopic observations showed that the UCI was sporophytic. The endogenous JA and JA-Ile levels were higher in YS than in SY at the studied time points. We found 7,251 differentially expressed genes from the transcriptome analysis. Our results highlighted that genes associated with JA biosynthesis and signaling, pollen tube growth, cell wall modification, starch and sucrose biosynthesis, and protein processing in endoplasmic reticulum pathways were differentially regulated between SY and YS. We discussed transcriptomic changes in the above-mentioned pathways regarding the observed microscopic and/or endogenous hormone levels. This is the first report on the elaboration of transcriptomic changes in longan reciprocal pollination combination showing UCI. The results presented here will enable the longan breeding community to better understand the mechanisms of UCI.

Genome-wide bioinformatics analysis of Cellulose Synthase gene family in common bean (Phaseolus vulgaris L.) and the expression in the pod development

Background

CesA and Csl gene families, which belong to the cellulose synthase gene superfamily, plays an important role in the biosynthesis of the plant cell wall. Although researchers have investigated this gene superfamily in several model plants, to date, no comprehensive analysis has been conducted in the common bean.

Results

In this study, we identified 39 putative cellulose synthase genes from the common bean genome sequence. Then, we performed a bioinformatics analysis of this gene family involving sequence alignment, phylogenetic analysis, gene structure, collinearity analysis and chromosome location. We found all members possess a cellulose_synt domain. Phylogenetic analysis revealed that these cellulose synthase genes may be classified into five subfamilies, and that members in the same subfamily share conserved exon-intron distribution and motif compositions. Abundant and distinct cis-acting elements in the 2 k basepairs upstream regulatory regions indicate that the cellulose synthase gene family may plays a vital role in the growth and development of common bean. Moreover, the 39 cellulose synthase genes are distributed on 10 of the 11 chromosomes. Additionally expression analysis shows that all CesA/Csl genes selected are constitutively expressed in the pod development.

Conclusions

This research reveals both the putative biochemical and physiological functions of cellulose synthase genes in common bean and implies the importance of studying non-model plants to understand the breadth and diversity of cellulose synthase genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01026-0.

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.

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

Members of the cellulose synthase (CesA) and cellulose synthase-like (Csl) families from the cellulose synthase gene superfamily participate in cellulose and hemicellulose synthesis in the plasma membrane. The members of this superfamily are vital for cell wall construction during plant growth and development. However, little is known about their function in pear fruit, a model for Rosaceae species and for fleshy fruit development. In our research, a total of 36 CesA/Csl family members were identified from the pear and were grouped into six subfamilies (CesA, CslB, CslC, CslD, CslE, and CslG) according to phylogenetic relationships. We performed a protein sequence physicochemical analysis, phylogenetic tree construction, a gene structure, a conserved domain, and chromosomal localization analysis. The results indicated that most of the CesA/Csl genes from pear are closely related to genes in Arabidopsis, but these families have unique characteristics in terms of their gene structure, chromosomal localization, phylogeny, and deduced protein sequences, suggesting that they have evolved through different processes. Tissue expression analysis results showed that most of the CesA/Csl genes were constitutively expressed at different levels in different organs. Furthermore, the expression levels of four genes (Pbr032894.2, Pbr016107.1, Pbr00518.1, and Pbr034218.1) tended to first increase and then decrease during fruit development, implying that these four genes may be involved in the development of stone cells of pear fruit. Our results may help elucidate the evolutionary history and functional differences of the CesA/Csl genes in pear and lay a foundation for further investigation of the CesA/Csl genes in pear and other Rosaceae species.

Identification of QTNs and Their Candidate Genes for 100-Seed Weight in Soybean (Glycine max L.) Using Multi-Locus Genome-Wide Association Studies

100-seed weight (100-SW) in soybeans is a yield component trait and controlled by multiple genes with different effects, but limited information is available for its quantitative trait nucleotides (QTNs) and candidate genes. To better understand the genetic architecture underlying the trait and improve the precision of marker-assisted selection, a total of 43,834 single nucleotide polymorphisms (SNPs) in 250 soybean accessions were used to identify significant QTNs for 100-SW in four environments and their BLUP values using six multi-locus and one single-locus genome-wide association study methods. As a result, a total of 218 significant QTNs were detected using multi-locus methods, whereas eight QTNs were identified by a single-locus method. Among 43 QTNs or QTN clusters identified repeatedly across various environments and/or approaches, all of them exhibited significant trait differences between their corresponding alleles, 33 were found in the genomic region of previously reported QTLs, 10 were identified as new QTNs, and three (qHSW-4-1, qcHSW-7-3, and qcHSW-10-4) were detected in all the four environments. The number of seed weight (SW) increasing alleles for each accession ranged from 8 (18.6%) to 36 (83.72%), and three accessions (Yixingwuhuangdou, Nannong 95C-5, and Yafanzaodou) had more than 35 SW increasing alleles. Among 36 homologous seed-weight genes in Arabidopsis underlying the above 43 stable QTNs, more importantly, Glyma05g34120, GmCRY1, and GmCPK11 had known seed-size/weight-related genes in soybean, and Glyma07g07850, Glyma10g03440, and Glyma10g36070 were candidate genes identified in this study. These results provide useful information for genetic foundation, marker-assisted selection, genomic prediction, and functional genomics of 100-SW.

Functional Analysis of the PgCesA3 White Spruce Cellulose Synthase Gene Promoter in Secondary Xylem

Cellulose is an essential structural component of the plant cell wall. Its biosynthesis involves genes encoding cellulose synthase enzymes and a complex transcriptional regulatory network. Three cellulose synthases have been identified in conifers as being potentially involved in secondary cell wall biosynthesis because of their preferential expression in xylem tissues; however, no direct functional association has been made to date. In the present work, we characterized the white spruce [Picea glauca (Moench) Voss] cellulose synthase PgCesA3 gene and 5' regulatory elements. Phylogenetic analysis showed that PgCesA1-3 genes grouped with secondary cell wall-associated Arabidopsis cellulose synthase genes, such as AtCesA8, AtCesA4, and AtCesA7. We produced transgenic spruce expressing the GUS reporter gene driven by the PgCesA3 promoter. We observed blue staining in differentiating xylem cells from stem and roots, and in foliar guard cells indicating that PgCesA3 is clearly involved in secondary cell wall biosynthesis. The promoter region sequence of PgCesA3 contained several putative MYB cis-regulatory elements including AC-I like motifs and secondary wall MYB-responsive element (SMRE); however, it lacked SMRE4, 7 and 8 that correspond to the sequences of AC-I, II, and III. Based on these findings and results of previous transient trans-activation assays that identified interactions between the PgCesA3 promoter and different MYB transcription factors, we performed electrophoretic mobility shift assays with MYB recombinant proteins and cis-regulatory elements present in the PgCesA3 promoter. We found that PgMYB12 bound to a canonical AC-I element identified in the Pinus taeda PAL promoter and two AC-I like elements. We hypothesized that the PgMYB12 could regulate PgCesA3 in roots based on previous expression results. This functional study of PgCesA3 sequences and promoter opens the door for future studies on the interaction between PgMYBs and the PgCesA3 regulatory elements.

Fiber Quality Improvement in Upland Cotton (Gossypium hirsutum L.): Quantitative Trait Loci Mapping and Marker Assisted Selection Application

Genetic improvement in fiber quality is one of the main challenges for cotton breeders. Fiber quality traits are controlled by multiple genes and are classified as complex quantitative traits, with a negative relationship with yield potential, so the genetic gain is low in traditional genetic improvement by phenotypic selection. The availability of Gossypium genomic sequences facilitates the development of high-throughput molecular markers, quantitative trait loci (QTL) fine mapping and gene identification, which helps us to validate candidate genes and to use marker assisted selection (MAS) on fiber quality in breeding programs. Based on developments of high density linkage maps, QTLs fine mapping, marker selection and omics, we have performed trait dissection on fiber quality traits in diverse populations of upland cotton. QTL mapping combined with multi-omics approaches such as, RNA sequencing datasets to identify differentially expressed genes have benefited the improvement of fiber quality. In this review, we discuss the application of molecular markers, QTL mapping and MAS for fiber quality improvement in upland cotton.

Genome-wide characterization of the cellulose synthase gene superfamily in Solanum lycopersicum

The cellulose synthase gene superfamily, which includes the cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, plays a vital role in the biosynthesis of cellulose and hemicellulose in plants. However, these genes have not been extensively studied in tomato (Solanum lycopersicum), a model for Solanaceae plants and for fleshy fruit development. Here, we identified and systematically analyzed 38 CesA/Csl family members that contained cellulose synthase domain regions, and categorized their encoded proteins into 6 subfamilies (CesA, CslA, CslB, CslD, CslE, and CslG) based on phylogenetic analysis. Most CesA/Csl genes from tomato are closely related to those from Arabidopsis, but the families have distinct features regarding gene structure, chromosome distribution and localization, phylogeny, and deduced protein sequence, indicating that they arose via different evolutionary process. Furthermore, expression analysis of CesA/Csl genes in different tissues at various developmental stages showed that most CesAs were constitutively expressed with differential expression levels in various organs; three CslD genes were expressed specifically in flowers, and four CesA and five Csl putative genes were preferentially expressed in fruits. Our results provide insight into the general characteristics of the CesA/Csl genes in tomato, and lay the foundation for further functional studies of CesA/Csl genes in tomato and other Solanaceae species.

Characterization of Cellulose Synthase A (CESA) Gene Family in Eudicots

Cellulose synthase A (CESA) is a key enzyme involved in the complex process of plant cell wall biosynthesis, and it remains a productive subject for research. We employed systems biology approaches to explore structural diversity of eudicot CESAs by exon-intron organization, mode of duplication, synteny, and splice site analyses. Using a combined phylogenetics and comparative genomics approach coupled with co-expression networks we reconciled the evolution of cellulose synthase gene family in eudicots and found that the basic forms of CESA proteins are retained in angiosperms. Duplications have played an important role in expansion of CESA gene family members in eudicots. Co-expression networks showed that primary and secondary cell wall modules are duplicated in eudicots. We also identified 230 simple sequence repeat markers in 103 eudicot CESAs. The 13 identified conserved motifs in eudicots will provide a basis for gene identification and functional characterization in other plants. Furthermore, we characterized (in silico) eudicot CESAs against senescence and found that expression levels of CESAs decreased during leaf senescence.

Transcription factors WRKY11 and WRKY17 are involved in abiotic stress responses in Arabidopsis

Plant WRKY transcription factors play a vital role in abiotic stress tolerance and regulation of plant defense responses. This study examined AtWRKY11 and AtWRKY17 expression under ABA, salt, and osmotic stress at different developmental stages in Arabidopsis. We used reverse transcriptase PCR, quantitative real-time PCR, and promoter:GUS lines to analyze expression. Both genes were upregulated in response to abiotic stress. Next, we applied the same stressors to seedlings of T-DNA insertion wrky11 and 17 knock-out mutants (single and double). Under stress, the mutants exhibited slower germination and compromised root growth compared with the wild type. In most cases, double-mutant seedlings were more affected than single mutants. These results suggest that wrky11 and wrky17 are not strictly limited to plant defense responses but are also involved in conferring stress tolerance.

Systems Identification and Characterization of Cell Wall Reassembly and Degradation Related Genes in Glycine max (L.) Merill, a Bioenergy Legume

Soybean is a promising biomass resource for generation of second-generation biofuels. Despite the utility of soybean cellulosic biomass and post-processing residues in biofuel generation, there is no comprehensive information available on cell wall loosening and degradation related gene families. In order to achieve enhanced lignocellulosic biomass with softened cell walls and reduced recalcitrance, it is important to identify genes involved in cell wall polymer loosening and degrading. Comprehensive genome-wide analysis of gene families involved in cell wall modifications is an efficient stratagem to find new candidate genes for soybean breeding for expanding biofuel industry. We report the identification of 505 genes distributed among 12 gene families related to cell wall loosening and degradation. 1262 tandem duplication events contributed towards expansion and diversification of studied gene families. We identified 687 Simple Sequence Repeat markers and 5 miRNA families distributed on 316 and 10 genes, respectively. Publically available microarray datasets were used to explore expression potential of identified genes in soybean plant developmental stages, 68 anatomical parts, abiotic and biotic stresses. Co-expression networks revealed transcriptional coordination of different gene families involved in cell wall loosening and degradation process.