Expression profiling and integrative analysis of the CESA/CSL superfamily in rice

Updating Insights into the Catalytic Domain Properties of Plant Cellulose synthase (CesA) and Cellulose synthase-like (Csl) Proteins

The wall is the last frontier of a plant cell involved in modulating growth, development and defense against biotic stresses. Cellulose and additional polysaccharides of plant cell walls are the most abundant biopolymers on earth, having increased in economic value and thereby attracted significant interest in biotechnology. Cellulose biosynthesis constitutes a highly complicated process relying on the formation of cellulose synthase complexes. Cellulose synthase (CesA) and Cellulose synthase-like (Csl) genes encode enzymes that synthesize cellulose and most hemicellulosic polysaccharides. Arabidopsis and rice are invaluable genetic models and reliable representatives of land plants to comprehend cell wall synthesis. During the past two decades, enormous research progress has been made to understand the mechanisms of cellulose synthesis and construction of the plant cell wall. A plethora of cesa and csl mutants have been characterized, providing functional insights into individual protein isoforms. Recent structural studies have uncovered the mode of CesA assembly and the dynamics of cellulose production. Genetics and structural biology have generated new knowledge and have accelerated the pace of discovery in this field, ultimately opening perspectives towards cellulose synthesis manipulation. This review provides an overview of the major breakthroughs gathering previous and recent genetic and structural advancements, focusing on the function of CesA and Csl catalytic domain in plants.

Genome-wide identification of candidate aquaporins involved in water accumulation of pomegranate outer seed coat

Aquaporins (AQPs) are a class of highly conserved integral membrane proteins that facilitate the uptake and transport of water and other small molecules across cell membranes. However, little is known about AQP genes in pomegranate (Punica granatum L.) and their potential role in water accumulation of the outer seed coat. We identified 38 PgrAQP genes in the pomegranate genome and divided them into five subfamilies based on a comparative analysis. Purifying selection played a role in the evolution of PgrAQP genes and a whole-genome duplication event in Myrtales may have contributed to the expansion of PgrTIP, PgrSIP, and PgrXIP genes. Transcriptome data analysis revealed that the PgrAQP genes exhibited different tissue-specific expression patterns. Among them, the transcript abundance of PgrPIPs were significantly higher than that of other subfamilies. The mRNA transcription levels of PgrPIP1.3, PgrPIP2.8, and PgrSIP1.2 showed a significant linear relationship with water accumulation in seed coats, indicating that PgrPIP1.3/PgrPIP2.8 located in the plasma membrane and PgrSIP1.2 proteins located on the tonoplast may be involved in water accumulation and contribute to the cell expansion of the outer seed coat, which then develops into juicy edible flesh. Overall, our results provided not only information on the characteristics and evolution of PgrAQPs, but also insights on the genetic improvement of outer seed coats.

De novo transcriptome analysis of white teak (Gmelina arborea Roxb) wood reveals critical genes involved in xylem development and secondary metabolism

BACKGROUND

Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees.

RESULTS

For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR.

CONCLUSIONS

We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.

The OsmiRNA166b-OsHox32 pair regulates mechanical strength of rice plants by modulating cell wall biosynthesis

The plant cell wall provides mechanical strength to support plant growth and development and to determine plant architecture. Cellulose and mixed-linkage glucan (MLG) present in primary cell wall, whereas cellulose, lignin and hemicellulose exist in secondary cell wall. Biosynthesis of the cell wall biopolymers needs the coordinated transcriptional regulation of all the biosynthetic genes. The module of OsmiR166b-OsHox32 regulates expression levels of the genes related to biosynthesis of MLG, cellulose and lignin. Transgenic plants knocking down miR166b (STTM166b) by short tandem target mimic (STTM) technology or overexpressing OsHox32 (OEHox32) showed drooping leaves and brittle culms. Due to accumulation of less lignin and cellulose, the cell wall thickness of STTM166b and OEHox32 plants was reduced when compared to that of wild-type plants. Overexpression of miR166b (OE166b) in rice plants or knocking down of OsHox32 by RNA interference (RNAiHox32) led to increased thickness of cell walls and enhanced mechanical strength of culms. Molecular analyses showed that OsmiR166b-OsHox32 pair regulates cell wall-related gene expression. OsHox32 binds to the promoters of OsCAD2 and OsCESA7 to suppress the expression levels of these two genes. The suppression of OsCAD2 is synergistic when OsHox32 is co-expressed with OSH15 (Oryza sativa homeobox 15). OsHox32 interacts with OSH15, and the START domain of OsHox32, harbouring the miR166b cleavage site, is required for the interaction of these two proteins. Our results demonstrate that OsmiR166b-OsHox32 pair plays important roles not only in plant growth and development but also in plant architecture by regulating the cell wall-related gene expression.

Genome-Wide Characterization and Expression Analysis of CAMTA Gene Family Under Salt Stress in Cucurbita moschata and Cucurbita maxima

Cucurbita Linn. vegetables have a long history of cultivation and have been cultivated all over the world. With the increasing area of saline–alkali soil, Cucurbita Linn. is affected by salt stress, and calmodulin-binding transcription activator (CAMTA) is known for its important biological functions. Although the CAMTA gene family has been identified in several species, there is no comprehensive analysis on Cucurbita species. In this study, we analyzed the genome of Cucurbita maxima and Cucurbita moschata. Five C. moschata calmodulin-binding transcription activators (CmoCAMTAs) and six C. maxima calmodulin-binding transcription activators (CmaCAMTAs) were identified, and they were divided into three subfamilies (Subfamilies I, II, and III) based on the sequence identity of amino acids. CAMTAs from the same subfamily usually have similar exon–intron distribution and conserved domains (CG-1, TIG, IQ, and Ank_2). Chromosome localization analysis showed that CmoCAMTAs and CmaCAMTAs were unevenly distributed across four and five out of 21 chromosomes, respectively. There were a total of three duplicate gene pairs, and all of which had experienced segmental duplication events. The transcriptional profiles of CmoCAMTAs and CmaCAMTAs in roots, stems, leaves, and fruits showed that these CAMTAs have tissue specificity. Cis-acting elements analysis showed that most of CmoCAMTAs and CmaCAMTAs responded to salt stress. By analyzing the transcriptional profiles of CmoCAMTAs and CmaCAMTAs under salt stress, it was shown that both C. moschata and C. maxima shared similarities against salt tolerance and that it is likely to contribute to the development of these species. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) further demonstrated the key role of CmoCAMTAs and CmaCAMTAs under salt stress. This study provided a theoretical basis for studying the function and mechanism of CAMTAs in Cucurbita Linn.

Genome-Wide Analysis of the G2-Like Transcription Factor Genes and Their Expression in Different Senescence Stages of Tobacco (Nicotiana tabacum L.)

The Golden2-like (GLK) transcription factors play important roles in regulating chloroplast growth, development, and senescence in plants. In this study, a total of 89 NtGLK genes (NtGLK1-NtGLK89) were identified in the tobacco genome and were classified into 10 subfamilies with variable numbers of exons and similar structural organizations based on the gene structure and protein motif analyses. Twelve segmental duplication pairs of NtGLK genes were identified in the genome. These NtGLK genes contain two conserved helix regions related to the HLH structure, and the sequences of the first helix region are less conserved than that of the second helix motif. Cis-regulatory elements of the NtGLK promoters were widely involved in light responsiveness, hormone treatment, and physiological stress. Moreover, a total of 206 GLK genes from tomato, tobacco, maize, rice, and Arabidopsis were retrieved and clustered into eight subgroups. Our gene expression analysis indicated that NtGLK genes showed differential expression patterns in tobacco leaves at five senescence stages. The expression levels of six NtGLK genes in group C were reduced, coinciding precisely with the increment of the degree of senescence, which might be associated with the function of leaf senescence of tobacco. Our results have revealed valuable information for further functional characterization of the GLK gene family in tobacco.

Genome-Wide Identification, Characterization and Expression Profiling of myosin Family Genes in Sebastes schlegelii

Myosins are important eukaryotic motor proteins that bind actin and utilize the energy of ATP hydrolysis to perform a broad range of functions such as muscle contraction, cell migration, cytokinesis, and intracellular trafficking. However, the characterization and function of myosin is poorly studied in teleost fish. In this study, we identified 60 myosin family genes in a marine teleost, black rockfish (Sebastes schlegelii), and further characterized their expression patterns. myosin showed divergent expression patterns in adult tissues, indicating they are involved in different types and compositions of muscle fibers. Among 12 subfamilies, S. schlegelii myo2 subfamily was significantly expanded, which was driven by tandem duplication events. The up-regulation of five representative genes of myo2 in the skeletal muscle during fast-growth stages of juvenile and adult S. schlegelii revealed their active role in skeletal muscle fiber synthesis. Moreover, the expression regulation of myosin during the process of myoblast differentiation in vitro suggested that they contribute to skeletal muscle growth by involvement of both myoblast proliferation and differentiation. Taken together, our work characterized myosin genes systemically and demonstrated their diverse functions in a marine teleost species. This lays foundation for the further studies of muscle growth regulation and molecular mechanisms of indeterminate skeletal muscle growth of large teleost fishes.

The reference genome of Miscanthus floridulus illuminates the evolution of Saccharinae

Miscanthus, a member of the Saccharinae subtribe that includes sorghum and sugarcane, has been widely studied as a feedstock for cellulosic biofuel production. Here, we report the sequencing and assembly of the Miscanthus floridulus genome by the integration of PacBio sequencing and Hi-C mapping, resulting in a chromosome-scale, high-quality reference genome of the genus Miscanthus. Comparisons among Saccharinae genomes suggest that Sorghum split first from the common ancestor of Saccharum and Miscanthus, which subsequently diverged from each other, with two successive whole-genome duplication events occurring independently in the Saccharum genus and one whole-genome duplication occurring in the Miscanthus genus. Fusion of two chromosomes occurred during rediploidization in M. floridulus and no significant subgenome dominance was observed. A survey of cellulose synthases (CesA) in M. floridulus revealed quite high expression of most CesA genes in growing stems, which is in agreement with the high cellulose content of this species. Resequencing and comparisons of 75 Miscanthus accessions suggest that M. lutarioriparius is genetically close to M. sacchariflorus and that M. floridulus is more distantly related to other species and is more genetically diverse. This study provides a valuable genomic resource for molecular breeding and improvement of Miscanthus and Saccharinae crops.

Genome-wide analysis and functional characterization of the Dof transcription factor family in rice (Oryza sativa L.)

Exhaustive searches of the rice genome have revealed 30 different potential OsDof (Oryza sativa DNA binding with One Finger) genes. Their subcellular localization, phylogenetic relationship, conserved motifs identification, chromosomal allocation, expression patterns, and interaction networks were analyzed. The Dof (DNA binding with One Finger) family of transcription factors represents a particular class of plant-specific transcriptional regulators, contain a highly conserved region of 50-52 amino acids (Dof domain) and involved in various plant developmental processes and response to various environmental stresses. Few (Oryza sativa) OsDof genes have been demonstrated previously for their biological functions but there is no comprehensive study on most of the Dof genes of rice. In the current study, exhaustive searches of the rice genome revealed 30 different potential OsDof genes, and then their subcellular localization, phylogenetic relationship, conserved motifs identification, chromosomal allocation, expression patterns, and interaction networks were analyzed. Phylogenetic analysis of Dof proteins in rice showed that they are distributed in 4 groups. By genome-wide observation of gene expression profiles, we found that OsDof genes showed significant variances in expression levels in different tissues across multiple developmental stages. Protein-protein correlation network analysis, shows a statically significant overlap of some OsDofs, suggesting their similar functions and a high degree of co-expression. The Dof family transcription factors have been reported for their involvement in the regulation of various gene expression processes in rice but still, most of the Dof genes are not characterized for their specific physiological functions. This study revealed useful information and clues about predicting the potential roles of OsDofs in rice by combining their genome-wide characterization, expression profiling, protein-protein interactions, and for further studies to develop high-quality rice varieties.

Genome-Wide Identification of Brassinosteroid Signaling Downstream Genes in Nine Rosaceae Species and Analyses of Their Roles in Stem Growth and Stress Response in Apple

Brassinosteroid signaling downstream genes regulate many important agronomic traits in rice. However, information on such genes is limited in Arabidopsis and Rosaceae species. We identified these genes in Arabidopsis and nine Rosaceae species. They were, respectively, named based on chromosomal locations. Segmental duplication and whole-genome duplication under purifying selection, as determined by Ka/Ks analysis, likely contributed to Rosaceae gene expansion. Apple (Malus domestica), Arabidopsis, and rice genes were generally similar, while several Rosaceae genes differed from their rice homologs in various characteristics, such as gene length, subcellular localization, transmembrane topology, conserved domains, secondary structures, and responses to external signals. The brassinosteroid downstream genes in apple were, respectively, induced or repressed by five phytohormones. Furthermore, these apple downstream genes were differentially expressed in different apple grafting combinations (“Nagafu No. 2”/“Malling 9” and “Nagafu No. 2”/“Nagafu No. 2”) and long–short shoot varieties (“Yanfu No. 6” and “Nagafu No. 2”). Responses of the MdBZR genes to diverse stress signals were examined and candidate hub genes were identified. These findings indicated that several brassinosteroid signaling downstream genes in Rosaceae functionally differed from their rice homologs, and certain apple genes may play roles in plant height and stress responses. This study provided valuable information and presented enriched biological theories on brassinosteroid signaling downstream genes in apple. Identification of such genes serve to help expand apple breeding and growth. This study provides useful information for brassinosteroid signaling downstream genes.

Molecular Characteristics, Synthase, and Food Application of Cereal β-Glucan

Cereal β-glucan is a type of valuable dietary fiber that mainly exists in the aleurone, subaleurone, and endosperm of some cereal grains. β-Glucan is acknowledged as a functional food ingredient owing to its multiple health benefits, including the prevention of diabetes, reduction in the incidence of cardiovascular disease, antitumor effects, antioxidant activities, and immunostimulation. It is well documented that cellulose synthase-like CslF/H/J genes encode synthases responsible for β-glucan biosynthesis in cereal grains. Recently, β-glucan has been widely applied as an emulsion stabilizer, thickening agent, fat substitute, and bioactive ingredient in the food industry due to its water solubility, viscosity, gelation property, and health benefits. Therefore, the present paper aims to review the molecular characteristics, synthase gene family, and food application of cereal β-glucan in recent years.

Genome-wide identification and expression analysis of the NF-Y transcription factor family in Populus

In the past few years, many studies have reported that the transcription factor Nuclear Factor Y (NF-Y) gene family plays important roles in embryonic development, photosynthesis, flowering time regulation and stress response, in various plants. Although the NF-Y gene family has been systematically studied in many species, little is known about NF-Y genes in Populus. In this study, the NF-Y gene family in the Populus genome was identified and its structural characteristics were described. Fifty-two NF-Y genes were authenticated in the Populus trichocarpa genome and categorized into three subfamilies (NF-YA/B/C) by phylogenetic analysis. Chromosomal localization of these genes revealed that they were distributed randomly across 17 of the 19 chromosomes. Segmental duplication played a vital role in the amplification of Populus NF-Y gene family. Moreover, microsynteny analysis indicated that, among Populus trichocarpa, Arabidopsis thaliana, Vitis vinifera and Carica papaya, NF-Y duplicated regions were more conserved between Populus trichocarpa and Vitis vinifera. Redundant stress-related cis-elements were also found in the promoters of most 13 NF-YA genes and their expression levels varied widely following drought, salt, ABA and cold treatments. Subcellular localization experiments in tobacco showed that PtNF-YA3 was localized in nucleus and cytomembrane, while PtNF-YA4 was only in the nucleus in tobacco. According to the transcriptional activity experiments, neither of them had transcriptional activity in yeast. In summary, a comprehensive analysis of the Populus NF-Y gene family was performed to establish a theoretical basis for further functional studies on this family.

Curled Flag Leaf 2, Encoding a Cytochrome P450 Protein, Regulated by the Transcription Factor Roc5, Influences Flag Leaf Development in Rice

Moderate curling generally causes upright leaf blades, which favors the establishment of ideal plant architecture and increases the photosynthetic efficiency of the population, both of which are desirable traits for super hybrid rice (Oryza sativa L.). In this study, we identified a novel curled-leaf mutant, curled flag leaf 2 (cfl2), which shows specific curling at the base of the flag leaf owing to abnormal epidermal development, caused by enlarged bulliform cells and increased number of papillae with the disordered distribution. Map-based cloning reveals that CFL2 encodes a cytochrome P450 protein and corresponds to the previously reported OsCYP96B4. CFL2 was expressed in all analyzed tissues with differential abundance and was downregulated in the clf1 mutant [a mutant harbors a mutation in the homeodomain leucine zipper IV (HD-ZIP IV) transcription factor Roc5]. Yeast one-hybrid and transient expression assays confirm that Roc5 could directly bind to the cis-element L1 box in the promoter of CFL2 before activating CFL2 expression. RNA sequencing reveals that genes associated with cellulose biosynthesis and cell wall-related processes were significantly upregulated in the cfl2 mutant. The components of cell wall, such as lignin, cellulose, and some kinds of monosaccharide, were altered dramatically in the cfl2 mutant when compared with wild-type "Jinhui10" (WT). Taken together, CFL2, as a target gene of Roc5, plays an important role in the regulation of flag leaf shape by influencing epidermis and cell wall development.

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.

The TaCslA12 gene expressed in the wheat grain endosperm synthesizes wheat-like mannan when expressed in yeast and Arabidopsis

Mannan is a class of cell wall polysaccharides widespread in the plant kingdom. Mannan structure and properties vary according to species and organ. The cell walls of cereal grains have been extensively studied due to their role in cereal processing and to their beneficial effect on human health as dietary fiber. Recently, we showed that mannan in wheat (Triticum aestivum) grain endosperm has a linear structure of β-1,4-linked mannose residues. The aim of this work was to study the biosynthesis and function of wheat grain mannan. We showed that mannan is deposited in the endosperm early during grain development, and we identified candidate mannan biosynthetic genes expressed in the endosperm. The functional study in wheat was unsuccessful therefore our best candidate genes were expressed in heterologous systems. The endosperm-specificTaCslA12 gene expressed in Pichia pastoris and in an Arabidopsis thaliana mutant depleted in glucomannan led to the production of wheat-like linear mannan lacking glucose residues and with moderate acetylation. Therefore, this gene encodes a mannan synthase and is likely responsible for the synthesis of wheat endosperm mannan.

Genome-wide member identification, phylogeny and expression analysis of PEBP gene family in wheat and its progenitors

The phosphatidylethanolamine binding protein (PEBP) family comprises ancient proteins found throughout the biosphere that play an important role in plant growth and development, flowering, seed development and dormancy. However, not all PEBP genes have been identified or analyzed in common wheat (Triticum aestivum L.) and its progenitors. In this study, we identified the PEBP genes in common wheat, Triticum dicoccoides, Triticum urartu and Aegilops tauschii by searching whole genome sequences, and characterized these genes by phylogenetic and transcriptome analyses. A total of 76, 38, 16 and 22 PEBP genes were identified in common wheat, T. dicoccoides, T. urartu and Ae. tauschii, respectively. Phylogenetic analysis classified the PEBP genes into four subfamilies (PEBP-like, MFT-like, TFL-like and FT-like); the PEBP-like subfamily was identified as a new subfamily with genes in this subfamily were conserved in plants. Group 2, 3 and 5 chromosomes of common wheat and its progenitors contained more PEBP genes than other chromosomes. The PEBP genes were conserved in wheat during evolution, and tandem duplication played a more important role in the amplification of PEBP genes than segmental duplication. Furthermore, transcriptome analysis revealed that PEBP genes showed tissue/organ-specific expression profiles and some PEBP genes were induced to express by biotic stresses. Quantitative real-time PCR (qRT-PCR) analysis revealed that seven randomly selected PEBP genes expressed differently during seed germination under cold, drought, flood, heat and salt stress treatments, and five of these genes (TaPEBP1, TaPEBP5, TaPEBP9, TaPEBP66 and TaPEBP69) showed significantly higher expression under different stress treatments, indicating that these genes play important roles during seed germination under stress conditions.

Genome-Wide Identification and Expression Analysis of the Class III Peroxidase Gene Family in Potato (Solanum tuberosum L.)

Class III peroxidases (PRXs) are plant-specific enzymes and play important roles in plant growth, development and stress response. In this study, a total of 102 non-redundant PRX gene members (StPRXs) were identified in potato (Solanum tuberosum L.). They were divided into 9 subfamilies based on phylogenetic analysis. The members of each subfamily were found to contain similar organizations of the exon/intron structures and protein motifs. The StPRX genes were not equally distributed among chromosomes. There were 57 gene pairs of segmental duplication and 26 gene pairs of tandem duplication. Expression pattern analysis based on the RNA-seq data of potato from public databases indicated that StPRX genes were expressed differently in various tissues and responded specifically to heat, salt and drought stresses. Most of the StPRX genes were expressed at significantly higher levels in root than in other tissues. In addition, real-time quantitative PCR (qRT-PCR) analysis for 7 selected StPRX genes indicated that these genes displayed various expression levels under abiotic stresses. Our results provide valuable information for better understanding the evolution of StPRX gene family in potato and lay the vital foundation for further exploration of PRX gene function in plants.

New insights into the evolution and functional divergence of the CIPK gene family in Saccharum

Background

Calcineurin B-like protein (CBL)-interacting protein kinases (CIPKs) are the primary components of calcium sensors, and play crucial roles in plant developmental processes, hormone signaling transduction, and in the response to exogenous stresses.

Results

In this study, 48 CIPK genes (SsCIPKs) were identified from the genome of Saccharum spontaneum. Phylogenetic reconstruction suggested that the SsCIPK gene family may have undergone six gene duplication events from the last common ancestor (LCA) of SsCIPKs. Whole-genome duplications (WGDs) served as the driving force for the amplification of SsCIPKs. The Nonsynonymous to synonymous substitution ratio (Ka/Ks) analysis showed that the duplicated genes were possibly under strong purifying selection pressure. The divergence time of these duplicated genes had an average duplication time of approximately 35.66 Mya, suggesting that these duplication events occurred after the divergence of the monocots and eudicots (165 Mya). The evolution of gene structure analysis showed that the SsCIPK family genes may involve intron losses. Ten ScCIPK genes were amplified from sugarcane (Saccharum spp. hybrids). The results of real-time quantitative polymerase chain reaction (qRT-PCR) demonstrated that these ten ScCIPK genes had different expression patterns under abscisic acid (ABA), polyethylene glycol (PEG), and sodium chloride (NaCl) stresses. Prokaryotic expression implied that the recombinant proteins of ScCIPK3, − 15 and − 17 could only slightly enhance growth under salinity stress conditions, but the ScCIPK21 did not. Transient N. benthamiana plants overexpressing ScCIPKs demonstrated that the ScCIPK genes were involved in responding to external stressors through the ethylene synthesis pathway as well as to bacterial infections.

Conclusions

In generally, a comprehensive genome-wide analysis of evolutionary relationship, gene structure, motif composition, and gene duplications of SsCIPK family genes were performed in S. spontaneum. The functional study of expression patterns in sugarcane and allogenic expressions in E. coli and N. benthamiana showed that ScCIPKs played various roles in response to different stresses. Thus, these results improve our understanding of the evolution of the CIPK gene family in sugarcane as well as provide a basis for in-depth functional studies of CIPK genes in sugarcane.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07264-9.

Genome-wide characterization and expression analyses of the auxin/indole-3-acetic acid (Aux/IAA) gene family in apple (Malus domestica)

Auxin is a necessary phytohormone for fruit development, accompanying the whole process of fruit growth and development. The Aux/IAA gene family is one of the early auxin-responsive gene families. At present, there were few reports involved in Aux/IAA genes in the fruit, especially in apple. In our study, we identified 42 MdAux/IAAs, phylogenetic analysis showed that Aux/IAA proteins from apple, tomato, and strawberry were clustered into 5 groups, 42 MdAux/IAAs randomly distributed on 13 chromosomes. Additionally, a comprehensive analysis of Aux/IAA gene family was completed, including gene structures, conserved motifs, phylogenetic analysis, chromosome mapping, orthologous identification, selection pressure analyses, synteny analysis, and protein interaction. We also tested the expression of MdAux/IAAs in different tissues and fruit development stages using quantitative reverse transcription-polymerase chain reaction (qRT-PCR), we found that the most members of Aux/IAA showed higher expression in seeds compared within stem and leaves, indicating they may play a role in regulating fruit development. We also declared that the expression of Aux/IAA gene was not consistent in the pericarp and seeds at the same developmental stage, 3 MdAux/IAAs of the pericarp were upregulated over 20-fold at 90 d and 5 MdAux/IAAs of the seeds were upregulated over 40-fold at 90 d. It was MdAux/IAA23 that showed extreme up-regulated expression in both pericarp and seeds. This study proved that the Aux/IAA gene families may perform a different function in apple fruit development and ripening, more importantly, it provided a foundation for further exploring the biological function of these MdAux/IAAs.

Genome-Wide Identification and Coexpression Network Analysis of DNA Methylation Pathway Genes and Their Differentiated Functions in Ginkgo biloba L

DNA methylation plays a vital role in diverse biological processes. DNA methyltransferases (DNMTs) genes and RNA-directed DNA methylation (RdDM)-related genes are key genes responsible for establishing and maintaining genome DNA methylation in plants. In the present study, we systematically identified nine GbDNMTs in Ginkgo biloba, including the three common families of GbMET1a/1b, GbCMT2, and GbDRMa/b/2a/2b/2c, and a fourth family—GbDNMT3—which is absent in most angiosperms. We also identified twenty RdDM-related genes, including four GbDCLs, six GbAGOs, and ten GbRDRs. Expression analysis of the genes showed the different patterns of individual genes, and 15 of 29 genes displayed expression change under five types of abiotic stress. Gene coexpression analysis and weighted gene co-expression network analysis (WGCNA) using 126 public transcriptomic datasets revealed that these genes were clustered into two groups. In group I, genes covered members from all six families which were preferentially expressed in the ovulate strobile and fruit. A gene ontology (GO) enrichment analysis of WGCNA modules indicated that group I genes were most correlated with the biological process of cell proliferation. Group II only consisted of RdDM-related genes, including GbDRMs, GbAGOs, and GbRDRs, but no GbDCLs, and these genes were specifically expressed in the cambium, suggesting that they may function in a dicer-like (DCL)-independent RdDM pathway in specific tissues. The gene module related to group II was most enriched in signal transduction, cell communication, and the response to the stimulus. These results demonstrate that gene family members could be conserved or diverged across species, and multi-member families in the same pathway may cluster into different modules to function differentially. The study provides insight into the DNA methylation genes and their possible functions in G. biloba, laying a foundation for the further study of DNA methylation in gymnosperms.

Characteristics and Expression Analyses of Trehalose-6-Phosphate Synthase Family in Prunus mume Reveal Genes Involved in Trehalose Biosynthesis and Drought Response

Trehalose and its key synthase (trehalose-6-phosphate synthase, TPS) can improve the drought tolerance of plants. However, little is known about the roles of trehalose and the TPS family in Prunus mume response to drought. In our study, we discovered that the trehalose content in leaf, root, and stem tissues significantly increased in P. mume in response to drought. Therefore, the characteristics and functions of the TPS family are worth investigating in P. mume. We identified nine TPS family members in P. mume, which were divided into two sub-families and characterized by gene structure, promoter elements, protein conserved domains, and protein motifs. We found that the Hydrolase_3 domain and several motifs were highly conserved in Group II instead of Group I. The distinctions between the two groups may result from selective constraints, which we estimated by the d_N/d_S (ω) ratio. The ω values of all the PmTPS family gene pairs were evaluated as less than 1, indicating that purity selection facilitated their divergence. A phylogenetic tree was constructed using 92 TPSs from 10 Rosaceae species, which were further divided into five clusters. Based on evolutionary analyses, the five clusters of TPS family proteins mainly underwent varied purity selection. The expression patterns of PmTPSs under drought suggested that the TPS family played an important role in the drought tolerance of P. mume. Combining the expression patterns of PmTPSs and the trehalose content changes in leaf, stem, and root tissues under normal conditions and drought stress, we found that the PmTPS2 and PmTPS6 mainly function in the trehalose biosynthesis in P. mume. Our findings not only provide valuable information about the functions of trehalose and TPSs in the drought response of P. mume, but they also contribute to the future drought breeding of P. mume.

Genome-wide identification of the maize 2OGD superfamily genes and their response to Fusarium verticillioides and Fusarium graminearum

In maize, eat rot and stalk rot caused by Fusarium verticillioides and Fusarium graminearum lead to contamination of moldy grains to produce mycotoxins. Identification of resistance genes against these pathogens for maize breeding is an effective way for disease control. Several 2-oxoglutarate-dependent dioxygenase (2OGD) proteins have been found to confer resistance to different pathogens in diverse plant species. However, little is known about the 2OGD superfamily in maize. Here, we identified 103 putative 2OGD genes in maize from a genome-wide analysis, and divided them into three classes - DOXA, DOXB, and DOXC. We further comprehensively investigated their gene structure, chromosome distribution, phylogenetic tree, gene-function enrichment, and expression profiles among different tissues. The genes encoding three 2OGD proteins, ACO, F3H, and NCS involved in ethylene biosynthesis, flavonoids biosynthesis, and alkaloids biosynthesis pathways, respectively, were identified to be induced by F. verticillioides and F. graminearum. The promoters of the three genes contain the binding sites for the transcription factor ZmDOF and ZmHSF, which are also induced by the two pathogens. The results imply that the three 2OGDs and the two transcription factors might be involved in the resistance to the two pathogens. This study provided a comprehensive understanding of the 2OGD superfamily in maize and laid the foundation for the further functional analysis of their roles in maize resistance to eat rot and stalk rot.

Plant cell wall integrity maintenance in model plants and crop species-relevant cell wall components and underlying guiding principles

The walls surrounding the cells of all land-based plants provide mechanical support essential for growth and development as well as protection from adverse environmental conditions like biotic and abiotic stress. Composition and structure of plant cell walls can differ markedly between cell types, developmental stages and species. This implies that wall composition and structure are actively modified during biological processes and in response to specific functional requirements. Despite extensive research in the area, our understanding of the regulatory processes controlling active and adaptive modifications of cell wall composition and structure is still limited. One of these regulatory processes is the cell wall integrity maintenance mechanism, which monitors and maintains the functional integrity of the plant cell wall during development and interaction with environment. It is an important element in plant pathogen interaction and cell wall plasticity, which seems at least partially responsible for the limited success that targeted manipulation of cell wall metabolism has achieved so far. Here, we provide an overview of the cell wall polysaccharides forming the bulk of plant cell walls in both monocotyledonous and dicotyledonous plants and the effects their impairment can have. We summarize our current knowledge regarding the cell wall integrity maintenance mechanism and discuss that it could be responsible for several of the mutant phenotypes observed.

Genome-wide investigation of WRKY transcription factors in sweet osmanthus and their potential regulation of aroma synthesis

WRKY transcription factors, one of the largest transcription factor families, play important roles in regulating the synthesis of secondary metabolites. In sweet osmanthus (Osmanthus fragrans), the monoterpenes have been demonstrated as the most important volatile compounds, and the W-box, which is the cognate binding site of WRKY transcription factors, could be identified in most of the terpene-synthesis-related genes' promoters. However, the role of the WRKY family in terpene synthesis in sweet osmanthus has rarely been examined. In this study, 154 WRKY genes with conserved WRKY domain were identified and classified into three groups. The group II was further divided into five subgroups, and almost all members of IId contained a plant zinc cluster domain. Eight OfWRKYs (OfWRKY7/19/36/38/42/84/95/139) were screened from 20 OfWRKYs for their flower-specific expression patterns in different tissues. Simultaneously, the expression patterns of OfWRKYs and emission patterns of volatile compounds during the flowering process were determined and gas chromatography-mass spectrometry results showed that monoterpenes, such as linalool and ocimene, accounted for the highest proportion, contributing to the floral scent of sweet osmanthus in two cultivars. In addition, correlation analysis revealed the expression patterns of OfWRKYs (OfWRKY7/19/36/139) were each correlated with distinct monoterpenes (linalool, linalool derivatives, ocimene and ocimene derivatives). Subcellular localization analysis showed that p35S::GFP-OfWRKY7/38/95/139 were localized in the nucleus and OfWRKY139 had very strong transactivation activity. Collectively, the results indicated potential roles of OfWRKY139 and OfWRKYs with plant zinc cluster domain in regulating synthesis of aromatic compounds in sweet osmanthus, laying the foundation for use of OfWRKYs to improve the aroma of ornamental plants.

Genome-wide identification and analysis of the Q-type C2H2 gene family in potato (Solanum tuberosum L.)

Plant Q-type C2H2 zinc finger proteins play an important role in plant tolerance to abiotic stresses. Although the Q-type C2H2 gene family has been identified in many plants, little is known about it in potato (Solanum tuberosum). In the present study, a total of 79 Q-type C2H2 proteins in potato (StZFPs) were identified and their distribution on chromosomes, gene structure, and conserved motifs was assessed. According to their protein structural and phylogenetic features, these 79 StZFPs were classified into 12 distinct subclasses. Collinearity analysis showed that tandem and segmental duplication events played a crucial role in expansion of the StZFP gene family. Synteny analysis indicated that 11 and 21 StZFP genes were orthologous to Arabidopsis and wheat (Triticum aestivum), respectively. RNA-seq data were used to analyze the tissue-specific expression and abiotic stress responses of the StZFP genes. Furthermore, we analyzed the expression of StZFP genes in drought-sensitive and drought-tolerant potato cultivars under drought stress. Subsequently, we used qPCR (Quantitative real-time-PCR) to calculate the relative expression of candidate genes in potato plantlets treated with NaCl (100 mM) and PEG 6000 (10% w/v) for 24 h. Such candidate genes could provide valuable information for abiotic stress resistance research in potato.

Characterization of the Tomato (Solanum lycopersicum) Pectin Methylesterases: Evolution, Activity of Isoforms and Expression During Fruit Ripening

Pectin methylesterase (PME, EC 3.1.1.11) is a hydrolytic enzyme of pectin that plays multiple roles in different plant development processes and responses to biotic stress. To characterize the molecular evolution and functional divergence of the PME gene family, a genome-wide analysis of the PME gene family in the tomato was performed. In total, 57 non-redundant PME genes were identified, and these PME genes were divided into five groups based on their phylogeneny; their classification was supported by similar gene structures and domain distributions. The PME genes were found to be unevenly distributed among 12 chromosomes of the tomato. In addition, 11 segmental duplication and 11 tandem duplication events occurred in these PME genes, implying that both contributed to the expansion of the tomato PME gene family. Non-synonymous/synonymous mutation ratio analysis revealed that positive selection played a key role in the functional divergence of PME genes. Interspecific collinear analysis indicated a large divergence in the PME gene family after the divergence of monocot and dicot plants in ancient times. Gene expression pattern analysis suggested that PMEs plays roles in the different parts of the tomato plant, including the fruit. Three newly identified candidate genes (Solyc03g083360, Solyc07g071600, and Solyc12g098340) may have functions during fruit ripening. Immunoassays suggested that the tomato isoform PE1 and PE2 may change pectin structure at cell junctions, which could be associated with fruit softening. In addition, our analysis indicate that two undescribed PE isoforms might be active in leaves and fruits. This study increases our understanding of the PME gene family in the tomato and may facilitate further functional analyses to elucidate PME function, especially during fruit ripening.

A xylan glucuronosyltransferase gene exhibits pleiotropic effects on cellular composition and leaf development in rice

Leaf chlorophyll content is an important physiological indicator of plant growth, metabolism and nutritional status, and it is highly correlated with leaf nitrogen content and photosynthesis. In this study, we report the cloning and identification of a xylan glucuronosyltransferase gene (OsGUX1) that affects relative chlorophyll content in rice leaf. Using a set of chromosomal segment substitution lines derived from a cross of wild rice accession ACC10 and indica variety Zhenshan 97 (ZS97), we identified numerous quantitative trait loci for relative chlorophyll content. One major locus of them for relative chlorophyll content was mapped to a 10.3-kb region that contains OsGUX1. The allele OsGUX1^AC from ACC10 significantly decreases nitrogen content and chlorophyll content of leaf compared with OsGUX1^ZS from ZS97. The overexpression of OsGUX1 reduced chlorophyll content, and the suppression of this gene increased chlorophyll content of rice leaf. OsGUX1 is located in Golgi apparatus, and highly expressed in seedling leaf and the tissues in which primary cell wall synthesis occurring. Our experimental data indicate that OsGUX1 is responsible for addition of glucuronic acid residues onto xylan and participates in accumulation of cellulose and hemicellulose in the cell wall deposition, thus thickening the primary cell wall of mesophyll cells, which might lead to reduced chlorophyll content in rice leaf. These findings provide insights into the association of cell wall components with leaf nitrogen content in rice.

Genome-wide identification and expression analysis of the PHD-finger gene family in Solanum tuberosum

Plant homeodomain (PHD) proteins are prevalent in eukaryotes and play important roles in plant growth, development and abiotic stress response. In this study, the comprehensive study of the PHD family (StPHD) was performed in potato (Solanum tuberosum L.). Seventy-two PHD genes (named StPHD1-72) were identified and grouped into 10 subfamilies based on phylogenetic analysis. Similar structure organizations were found within each subfamily according to the exon/intron structures and protein motif analysis. These genes were unequally scattered on the chromosomes of potato, with 9 pairs of segmental duplicated genes and 6 pairs of tandem duplicated genes showing that both segmental duplicated and tandem duplicated events contributed to the expansion of the potato PHD family. The gene ontology (GO) analysis suggests that StPHD mainly functioned at the intracellular level and was involved in various binding, metabolic and regulation processes. The analysis of expression patterns of StPHD genes showed that these genes were differentially expressed in 10 different tissues and responded specifically to heat, salt and drought stress based on the FPKM (Fragments per kilobase of transcript per million mapped reads) values of the RNA-seq data. Furthermore, the real-time quantitative PCR for 12 selected StPHD genes revealed the various levels of gene expression corresponding to abiotic stress. Our results provide useful information for a better understanding of PHD genes and provide the foundation for additional functional exploration of the potato PHD gene family.

Improving Lodging Resistance: Using Wheat and Rice as Classical Examples

One of the most chronic constraints to crop production is the grain yield reduction near the crop harvest stage by lodging worldwide. This is more prevalent in cereal crops, particularly in wheat and rice. Major factors associated with lodging involve morphological and anatomical traits along with the chemical composition of the stem. These traits have built up the remarkable relationship in wheat and rice genotypes either prone to lodging or displaying lodging resistance. In this review, we have made a comparison of our conceptual perceptions with foregoing published reports and proposed the fundamental controlling techniques that could be practiced to control the devastating effects of lodging stress. The management of lodging stress is, however, reliant on chemical, agronomical, and genetic factors that are reducing the risk of lodging threat in wheat and rice. But, still, there are many questions remain to be answered to elucidate the complex lodging phenomenon, so agronomists, breeders, physiologists, and molecular biologists require further investigation to address this challenging problem.

Genome-wide identification, characterization and expression analysis of the HD-Zip gene family in the stem development of the woody plant Prunus mume

The homeodomain-leucine zipper (HD-Zip) gene family, a group of plant-specific transcriptional factors (TFs), participates in regulating growth, development, and environmental responses. However, the characteristics and biological functions of HD-Zip genes in Prunus mume, which blooms in late winter or early spring, have not been reported. In this study, 32 HD-Zip genes, named PmHB1-PmHB32 based on their chromosomal positions, were identified in the genome of P. mume. These genes are distributed among seven chromosomes and are phylogenetically clustered into four major groups. Gene structure and motif composition were mostly conserved in each group. The Ka/Ks ratios showed that purifying selection has played a leading role in the long-term evolution of the genes, which maintained the function of this family. MicroRNA target site prediction indicated that the genes of the HD-Zip III subfamily may be regulated by miR165/166. Expression pattern analysis showed that the 32 genes were differentially expressed across five different tissues (leaf, flower bud, stem, fruit, and root) and at different stages of stem and leaf-bud development, suggesting that 10 of the genes may play important roles in stem development. Protein-protein interaction predictions showed that the subfamily III genes may regulate vascular development and shoot apical meristem (SAM) maintenance. Promoter analysis showed that the HD-Zip III genes might be involved in responses to light, hormones, and abiotic stressors and stem development. Taken together, our results provide an overview of the HD-Zip family in P. mume and lay the foundation for the molecular breeding of woody ornamental plants.

Genome-Wide Identification, Expression Pattern Analysis and Evolution of the Ces/Csl Gene Superfamily in Pineapple (Ananas comosus)

The cellulose synthase (Ces) and cellulose synthase-like (Csl) gene families belonging to the cellulose synthase gene superfamily, are responsible for the biosynthesis of cellulose and hemicellulose of the plant cell wall, and play critical roles in plant development, growth and evolution. However, the Ces/Csl gene family remains to be characterized in pineapple, a highly valued and delicious tropical fruit. Here, we carried out genome-wide study and identified a total of seven Ces genes and 25 Csl genes in pineapple. Genomic features and phylogeny analysis of Ces/Csl genes were carried out, including phylogenetic tree, chromosomal locations, gene structures, and conserved motifs identification. In addition, we identified 32 pineapple AcoCes/Csl genes with 31 Arabidopsis AtCes/Csl genes as orthologs by the syntenic and phylogenetic approaches. Furthermore, a RNA-seq investigation exhibited the expression profile of several AcoCes/Csl genes in various tissues and multiple developmental stages. Collectively, we provided comprehensive information of the evolution and function of pineapple Ces/Csl gene superfamily, which would be useful for screening out and characterization of the putative genes responsible for tissue development in pineapple. The present study laid the foundation for future functional characterization of Ces/Csl genes in pineapple.

Identification of TPX2 Gene Family in Upland Cotton and Its Functional Analysis in Cotton Fiber Development

Microtubules (MTs) are of importance to fiber development. The Xklp2 (TPX2) proteins as a class of microtubule-associated proteins (MAPs) play a key role in plant growth and development by regulating the dynamic changes of microtubules (MTs). However, the mechanism underlying this is unknown. The interactions between TPX2 proteins and tubulin protein, which are the main structural components, have not been studied in fiber development of upland cotton. Therefore, a genome-wide analysis of the TPX2 family was firstly performed in Gossypiumhirsutum L. This study identified 41 GhTPX2 sequences in the assembled G. hirsutum genome by a series of bioinformatic methods. Generally, this gene family is phylogenetically grouped into six subfamilies, and 41 G. hirsutum TPX2 genes (GhTPX2s) are distributed across 21 chromosomes. A heatmap of the TPX2 gene family showed that homologous GhTPX2 genes, GhWDLA2/7 and GhWDLA4/9, have large differences in expression levels between two upland cotton recombinant inbred lines (69307 and 69362) that are different in fiber quality at 15 and 20 days post anthesis. The relative data indicate that these four genes are down-regulated under oryzalin, which causes microtubule depolymerization, as determined via qRT-PCR. A subcellular localization experiment suggested that GhWDLA2 and GhWDLA7 are localized to the microtubule cytoskeleton, and GhWDLA4 and GhWDLA9 are only localized to the nucleus. However, only GhWDLA7 between GhWDLA2 and GhWDLA7 interacted with GhTUA2 in the yeast two-hybrid assay. These results lay the foundation for further function study of the TPX2 gene family.

Selection of Suitable Reference Genes for RT-qPCR Gene Expression Analysis in Siberian Wild Rye (Elymus sibiricus) under Different Experimental Conditions

Elymus sibiricus, which is a perennial and self-pollinated grass, is the typical species of the genus Elymus, which plays an important role in forage production and ecological restoration. No reports have, so far, systematically described the selection of optimal reference genes for reverse transcriptase quantitative real-time polymerase chain reaction (RT-qPCR) analysis in E. sibiricus. The goals of this study were to evaluate the expression stability of 13 candidate reference genes in different experimental conditions, and to determine the appropriate reference genes for gene expression analysis in E. sibiricus. Five methods including Delta Ct (ΔCt), BestKeeper, NormFinder, geNorm, and RefFinder were used to assess the expression stability of 13 potential reference genes. The results of the RefFinder analysis showed that TBP2 and HIS3 were the most stable reference genes in different genotypes. TUA2 and PP2A had the most stable expression in different developmental stages. TBP2 and PP2A were suitable reference genes in different tissues. Under salt stress, ACT2 and TBP2 were identified as the most stable reference genes. ACT2 and TUA2 showed the most stability under heat stress. For cold stress, PP2A and ACT2 presented the highest degree of expression stability. DNAJ and U2AF were considered as the most stable reference genes under osmotic stress. The optimal reference genes were selected to investigate the expression pattern of target gene CSLE6 in different conditions. This study provides suitable reference genes for further gene expression analysis using RT-qPCR in E. sibiricus.

Genome-wide identification, characterization and expression analysis of the non-specific lipid transfer proteins in potato

BACKGROUND

Plant non-specific lipid transfer proteins (nsLTPs) are small, basic proteins that are abundant in higher plants. They have been reported to play an important role in various plant physiological processes, such as lipid transfer, signal transduction, and pathogen defense. To date, a comprehensive analysis of the potato nsLTP gene family is still lacking after the completion of potato (Solanum tuberosum L.) genome sequencing. A genome-wide characterization, classification and expression analysis of the StnsLTP gene family was performed in this study.

RESULTS

In this study, a total of 83 nsLTP genes were identified and categorized into eight types based on Boutrot's method. Multiple characteristics of these genes, including phylogeny, gene structures, conserved motifs, protein domains, chromosome locations, and cis-elements in the promoter sequences, were analyzed. The chromosome distribution and the collinearity analyses suggested that the expansion of the StnsLTP gene family was greatly enhanced by the tandem duplications. Ka/Ks analysis showed that 47 pairs of duplicated genes tended to undergo purifying selection during evolution. Moreover, the expression of StnsLTP genes in various tissues was analyzed by using RNA-seq data and verified by quantitative real-time PCR, revealing that the StnsLTP genes were mainly expressed in younger tissues. These results indicated that StnsLTPs may played significant and functionally varied roles in the development of different tissues.

CONCLUSION

In this study, we comprehensively analyzed nsLTPs in potato, providing valuable information to better understand the functions of StnsLTPs in different tissues and pathways, especially in response to abiotic stress.

Meta-analysis of drought-tolerant genotypes in Oryza sativa: A network-based approach

BACKGROUND

Drought is a severe environmental stress. It is estimated that about 50% of the world rice production is affected mainly by drought. Apart from conventional breeding strategies to develop drought-tolerant crops, innovative computational approaches may provide insights into the underlying molecular mechanisms of stress response and identify drought-responsive markers. Here we propose a network-based computational approach involving a meta-analytic study of seven drought-tolerant rice genotypes under drought stress.

RESULTS

Co-expression networks enable large-scale analysis of gene-pair associations and tightly coupled clusters that may represent coordinated biological processes. Considering differentially expressed genes in the co-expressed modules and supplementing external information such as resistance/tolerance QTLs, transcription factors, network-based topological measures, we identify and prioritize drought-adaptive co-expressed gene modules and potential candidate genes. Using the candidate genes that are well-represented across the datasets as 'seed' genes, two drought-specific protein-protein interaction networks (PPINs) are constructed with up- and down-regulated genes. Cluster analysis of the up-regulated PPIN revealed ABA signalling pathway as a central process in drought response with a probable crosstalk with energy metabolic processes. Tightly coupled gene clusters representing up-regulation of core cellular respiratory processes and enhanced degradation of branched chain amino acids and cell wall metabolism are identified. Cluster analysis of down-regulated PPIN provides a snapshot of major processes associated with photosynthesis, growth, development and protein synthesis, most of which are shut down during drought. Differential regulation of phytohormones, e.g., jasmonic acid, cell wall metabolism, signalling and posttranslational modifications associated with biotic stress are elucidated. Functional characterization of topologically important, drought-responsive uncharacterized genes that may play a role in important processes such as ABA signalling, calcium signalling, photosynthesis and cell wall metabolism is discussed. Further transgenic studies on these genes may help in elucidating their biological role under stress conditions.

CONCLUSION

Currently, a large number of resources for rice functional genomics exist which are mostly underutilized by the scientific community. In this study, a computational approach integrating information from various resources such as gene co-expression networks, protein-protein interactions and pathway-level information is proposed to provide a systems-level view of complex drought-responsive processes across the drought-tolerant genotypes.

Genome-wide identification, classification and expression analysis of the JmjC domain-containing histone demethylase gene family in maize

Background

Histone methylation mainly occurs on the lysine residues and plays a crucial role during flowering and stress responses of plants, through changing the methylation status or ratio of lysine residues. Histone lysine residues of plants can arise in three forms of methylation (single, double and triple) and the corresponding demethylation can also ensue on certain occasions, by which the plants can accommodate the homeostasis of histone methylation by means of lysine methyltransferase and demethylase. The JmjC domain-containing proteins, an important family of histone lysine demethylases, play a vital role in maintaining homeostasis of histone methylation in vivo.

Results

In this study, we have identified 19 JmjC domain-containing histone demethylase (JHDM) proteins in maize. Based on structural characteristics and a comparison of phylogenetic relationships of JHDM gene families from Arabidopsis, rice and maize, all 19 JHDM proteins in maize were categorized into three different subfamilies. Furthermore, chromosome location and schematic structure revealed an unevenly distribution on chromosomes and structure features of ZmJMJ genes in maize, respectively. Eventually, the 19 ZmJMJ genes displayed different expression patterns at diverse developmental stages of maize based on transcriptome analysis. Further, quantitative real-time PCR analysis showed that all 19 ZmJMJ genes were responsive to heat stress treatment, suggesting their potential roles in heat stress response.

Conclusions

Overall, our study will serve to present an important theoretical basis for future functional verification of JHDM genes to further unravel the mechanisms of epigenetic regulation in plants.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5633-1) contains supplementary material, which is available to authorized users.

Genome-wide analysis of the Chinese cabbage IQD gene family and the response of BrIQD5 in drought resistance

KEY MESSAGE

Thirty-five IQD genes were identified and analysed in Chinese cabbage and BrIQD5 transgenic plants enhanced the drought resistance of plants. The IQD (IQ67-domain) family plays an important role in various abiotic stress responses in plant species. However, the roles of IQD genes in the Chinese cabbage response to abiotic stress remain unclear. Here, 35 IQD genes, from BrIQD1 to BrIQD35, were identified in Chinese cabbage (Brassica rapa ssp. pekinensis). Based on the phylogenetic analysis, these genes were clustered into three subfamilies (I-III), and members within the same subfamilies shared conserved exon-intron distribution and motif composition. The 35 BrIQD genes were unevenly distributed on 9 of the 10 chromosomes with 4 segmental duplication events. Ka/Ks ratios showed that the duplicated BrIQDs had mainly experienced strong purifying selection. Quantitative real-time polymerase chain reaction of 35 BrIQDs under PEG6000 indicated that BrIQD5 was significantly induced by PEG6000. To verify BrIQD5 function, BrIQD5 was heterologously overexpressed in tobacco and was silenced in Chinese cabbage. BrIQD5-overexpressed plants showed more tolerance to drought stress than wild-type plants, while BrIQD5-silenced plants in Chinese cabbage showed decreased drought tolerance. Additionally, six BrIQD5 potential interactive proteins were isolated by the yeast two-hybrid assay, including BrCaMa, BrCaMb and four other stress-related proteins. Motif IQ1 of BrIQD5 is important for the interaction with BrCaMa and BrCaMb, and the isoleucine in motif IQ1 is an essential amino acid for calmodulin binding to BrIQD5. The identification and cloning of the new Chinese cabbage drought tolerance genes will promote the drought-resistant breeding of Chinese cabbage and help to better understand the mechanism of IQD involved in the drought tolerance of plants.

Computational approach to understand molecular mechanism involved in BPH resistance in Bt- rice plant

In silico approach was utilised to identify differentially expressed key hub genes during BPH infestation on Bt rice plant, under laboratory conditions. Re-analysis of GSE74745 data with in-house R scripts and STRING database reveals that only 5 key hub genes, namely Os05g0176100, Os06g0683200, Os07g0208500, Os07g0252400 and Os07g0424400, belonging to cellulose synthase family, are differentially expressed and have confidence score ≥0.9 among themselves. Conserve domain analysis of all proteins encoded via these 5 key hub genes reveals that they have a common cellulose synthase domain, in which "Plant-Conserved Region" (PCR) is highly conserved. After binding with other domains of cellulose synthase proteins or other accessory proteins, like sucrose synthase, PCR serves as a metabolic channel to deliver UDP-Glucose, which is the main substrate for cellulose synthesis, into the active site of cellulose synthase and initiate cellulose synthesis. Simulation study of recently solved topological model of PCR [PDB ID: 5JNP] and molecular docking studies of PCR with UDP-glucose reveals that, during BPH infestation, in nearby phloem tissue where BPH suck sap, there is an increase interaction of UDP-glucose with PCR and other accessory proteins which in turn increases both the stability of PCR and the production of cellulose, finally causing callose deposition at that site and hence causing longer nymphal developmental period and lower fertility of BPH infested on Bt rice. In near future, these differentially identified 5 hub genes could be possible targets for controlling BPH infestation in rice plant under field conditions and increasing rice yield globally.

De Novo Transcriptome Assembly of Agave H11648 by Illumina Sequencing and Identification of Cellulose Synthase Genes in Agave Species

Agave plants are important crassulacean acid metabolism (CAM) plants with multiple agricultural uses, such as being used in tequila and fiber production. Agave hybrid H11648 ((A. amaniensis Trel. and Nowell × A. angustifolia Haw.) × A. amaniensis) is the main cultivated Agave species for fiber production in large tropical areas around the world. In this study, we conducted a transcriptome analysis of A. H11648. About 49.25 million clean reads were obtained by Illumina paired-end sequencing. De novo assembly produced 148,046 unigenes with more than 40% annotated in public databases, or matched homologs in model plants. More homologous gene pairs were found in Asparagus genome than in Arabidopsis or rice, which indicated a close evolutionary relationship between Asparagus and A. H11648. CAM-related gene families were also characterized as previously reported in A. americana. We further identified 12 cellulose synthase genes (CesA) in Asparagus genome and 38 CesA sequences from A. H11648, A. americana, A. deserti and A. tequilana. The full-length CesA genes were used as references for the cloning and assembly of their homologs in other Agave species. As a result, we obtained CesA1/3/4/5/7 genes with full-length coding region in the four Agave species. Phylogenetic and expression analysis revealed a conserved evolutionary pattern, which could not explain the distinct fiber traits in different Agave species. We inferred that transcriptional regulation might be responsible for Agave fiber development. This study represents the transcriptome of A. H11648, which would expand the number of Agave genes and benefit relevant studies of Agave fiber development.

Cell Wall Enzymes in Zygnema circumcarinatum UTEX 1559 Respond to Osmotic Stress in a Plant-Like Fashion

Previous analysis of charophyte green algal (CGA) genomes and transcriptomes for specific protein families revealed that numerous land plant characteristics had already evolved in CGA. In this study, we have sequenced and assembled the transcriptome of Zygnema circumcarinatum UTEX 1559, and combined its predicted protein sequences with those of 13 additional species [five embryophytes (Emb), eight charophytes (Cha), and two chlorophytes (Chl) as the outgroup] for a comprehensive comparative genomics analysis. In total 25,485 orthologous gene clusters (OGCs, equivalent to protein families) of the 14 species were classified into nine OGC groups. For example, the Cha+Emb group contains 4,174 OGCs found in both Cha and Emb but not Chl species, representing protein families that have evolved in the common ancestor of Cha and Emb. Different OGC groups were subjected to a Gene Ontology (GO) enrichment analysis with the Chl+Cha+Emb group (including 5,031 OGCs found in Chl and Cha and Emb) as the control. Interestingly, nine of the 20 top enriched GO terms in the Cha+Emb group are cell wall-related, such as biological processes involving celluloses, pectins, lignins, and xyloglucans. Furthermore, three glycosyltransferase families (GT2, 8, 43) were selected for in-depth phylogenetic analyses, which confirmed their presence in UTEX 1559. More importantly, of different CGA groups, only Zygnematophyceae has land plant cellulose synthase (CesA) orthologs, while other charophyte CesAs form a CGA-specific CesA-like (Csl) subfamily (likely also carries cellulose synthesis activity). Quantitative real-time-PCR experiments were performed on selected GT family genes in UTEX 1559. After osmotic stress treatment, significantly elevated expression was found for GT2 family genes ZcCesA, ZcCslC and ZcCslA-like (possibly mannan and xyloglucan synthases, respectively), as well as for GT8 family genes (possibly pectin synthases). All these suggest that the UTEX 1559 cell wall polysaccharide synthesis-related genes respond to osmotic stress in a manner that is similar to land plants.

Overproduction of native endo-β-1,4-glucanases leads to largely enhanced biomass saccharification and bioethanol production by specific modification of cellulose features in transgenic rice

BACKGROUND

Genetic modification of plant cell walls has been implemented to reduce lignocellulosic recalcitrance for biofuel production. Plant glycoside hydrolase family 9 (GH9) comprises endo-β-1,4-glucanase in plants. Few studies have examined the roles of GH9 in cell wall modification. In this study, we independently overexpressed two genes from GH9B subclasses (OsGH9B1 and OsGH9B3) and examined cell wall features and biomass saccharification in transgenic rice plants.

RESULTS

Compared with the wild type (WT, Nipponbare), the OsGH9B1 and OsGH9B3 transgenic rice plants, respectively, contained much higher OsGH9B1 and OsGH9B3 protein levels and both proteins were observed in situ with nonspecific distribution in the plant cells. The transgenic lines exhibited significantly increased cellulase activity in vitro than the WT. The OsGH9B1 and OsGH9B3 transgenic plants showed a slight alteration in three wall polymer compositions (cellulose, hemicelluloses, and lignin), in their stem mechanical strength and biomass yield, but were significantly decreased in the cellulose degree of polymerization (DP) and lignocellulose crystalline index (CrI) by 21-22%. Notably, the crude cellulose substrates of the transgenic lines were more efficiently digested by cellobiohydrolase (CBHI) than those of the WT, indicating the significantly increased amounts of reducing ends of β-1,4-glucans in cellulose microfibrils. Finally, the engineered lines generated high sugar yields after mild alkali pretreatments and subsequent enzymatic hydrolysis, resulting in the high bioethanol yields obtained at 22.5% of dry matter.

CONCLUSIONS

Overproduction of OsGH9B1/B3 enzymes should have specific activity in the postmodification of cellulose microfibrils. The increased reducing ends of β-1,4-glucan chains for reduced cellulose DP and CrI positively affected biomass enzymatic saccharification. Our results demonstrate a potential strategy for genetic modification of cellulose microfibrils in bioenergy crops.

Evolution of Cell Wall Polymers in Tip-Growing Land Plant Gametophytes: Composition, Distribution, Functional Aspects and Their Remodeling

During evolution of land plants, the first colonizing species presented leafy-dominant gametophytes, found in non-vascular plants (bryophytes). Today, bryophytes include liverworts, mosses, and hornworts. In the first seedless vascular plants (lycophytes), the sporophytic stage of life started to be predominant. In the seed producing plants, gymnosperms and angiosperms , the gametophytic stage is restricted to reproduction. In mosses and ferns, the haploid spores germinate and form a protonema, which develops into a leafy gametophyte producing rhizoids for anchorage, water and nutrient uptakes. The basal gymnosperms (cycads and Ginkgo) reproduce by zooidogamy. Their pollen grains develop a multi-branched pollen tube that penetrates the nucellus and releases flagellated sperm cells that swim to the egg cell. The pollen grain of other gymnosperms (conifers and gnetophytes) as well as angiosperms germinates and produces a pollen tube that directly delivers the sperm cells to the ovule (siphonogamy). These different gametophytes, which are short or long-lived structures, share a common tip-growing mode of cell expansion. Tip-growth requires a massive cell wall deposition to promote cell elongation, but also a tight spatial and temporal control of the cell wall remodeling in order to modulate the mechanical properties of the cell wall. The growth rate of these cells is very variable depending on the structure and the species, ranging from very slow (protonemata, rhizoids, and some gymnosperm pollen tubes), to a slow to fast-growth in other gymnosperms and angiosperms. In addition, the structural diversity of the female counterparts in angiosperms (dry, semi-dry vs wet stigmas, short vs long, solid vs hollow styles) will impact the speed and efficiency of sperm delivery. As the evolution and diversity of the cell wall polysaccharides accompanied the diversification of cell wall structural proteins and remodeling enzymes, this review focuses on our current knowledge on the biochemistry, the distribution and remodeling of the main cell wall polymers (including cellulose, hemicelluloses, pectins, callose, arabinogalactan-proteins and extensins), during the tip-expansion of gametophytes from bryophytes, pteridophytes (lycophytes and monilophytes), gymnosperms and the monocot and eudicot angiosperms.

Genome-wide identification of SERK genes in apple and analyses of their role in stress responses and growth

BACKGROUND

Somatic embryogenesis receptor-like kinases (SERKs) are leucine-rich repeat receptor-like kinases associated with various signaling pathways. These kinases have a relationship with stress signals, and they are also believed to be important for regulating plant growth. However, information about this protein family in apple is limited.

RESULTS

Twelve apple SERK genes distributed across eight chromosomes were identified. These genes clustered into three distinct groups in a phylogenetic analysis. All of the encoded proteins contained typical SERK domains. The chromosomal locations, gene/protein structures, synteny, promoter sequences, protein-protein interactions, and physicochemical characteristics of MdSERK genes were analyzed. Bioinformatics analyses demonstrated that gene duplications have likely contributed to the expansion and evolution of SERK genes in the apple genome. Six homologs of SERK genes were identified between apple and Arabidopsis. Quantitative real-time PCR analyses revealed that the MdSERK genes showed different expression patterns in various tissues. Eight MdSERK genes were responsive to stress signals, such as methyl jasmonate, salicylic acid, abscisic acid, and salt (NaCl). The application of exogenous brassinosteroid and auxin increased the growth and endogenous hormone contents of Malus hupehensis seedlings. The expression levels of seven MdSERK genes were significantly upregulated by brassinosteroid and auxin. In addition, several MdSERK genes showed higher expression levels in standard trees of 'Nagafu 2' (CF)/CF than in dwarf trees of CF/'Malling 9' (M.9), and in CF than in the spur-type bud mutation "Yanfu 6" (YF).

CONCLUSION

This study represents the first comprehensive investigation of the apple SERK gene family. These data indicate that apple SERKs may function in adaptation to adverse environmental conditions and may also play roles in controlling apple tree growth.

Genome-wide identification and expression analyses of the pectate lyase (PEL) gene family in cotton (Gossypium hirsutum L.)

Background

Pectin is a major component and structural polysaccharide of the primary cell walls and middle lamella of higher plants. Pectate lyase (PEL, EC 4.2.2.2), a cell wall modification enzyme, degrades de-esterified pectin for cell wall loosening, remodeling and rearrangement. Nevertheless, there have been few studies on PEL genes and no comprehensive analysis of the PEL gene family in cotton.

Results

We identified 53, 42 and 83 putative PEL genes in Gossypium raimondii (D5), Gossypium arboreum (A2), and Gossypium hirsutum (AD1), respectively. These PEL genes were classified into five subfamilies (I-V). Members from the same subfamilies showed relatively conserved gene structures, motifs and protein domains. An analysis of gene chromosomal locations and gene duplication revealed that segmental duplication likely contributed to the expansion of the GhPELs. The 2000 bp upstream sequences of all the GhPELs contained auxin response elements. A transcriptomic data analysis showed that 62 GhPELs were expressed in various tissues. Notably, most (29/32) GhPELs of subfamily IV were preferentially expressed in the stamen, and five GhPELs of subfamily V were prominently expressed at the fiber elongation stage. In addition, qRT-PCR analysis revealed the expression characteristics of 24 GhPELs in four pollen developmental stages and significantly different expression of some GhPELs between long- and short-fiber cultivars. Moreover, some members were responsive to IAA treatment. The results indicate that GhPELs play significant and functionally diverse roles in the development of different tissues.

Conclusions

In this study, we comprehensively analyzed PELs in G. hirsutum, providing a foundation to better understand the functions of GhPELs in different tissues and pathways, especially in pollen, fiber and the auxin signaling pathway.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5047-5) contains supplementary material, which is available to authorized users.

Promoters of Arabidopsis Casein kinase I-like 2 and 7 confer specific high-temperature response in anther

KEY MESSAGE: (1) We systematically analyze the promoter activities of AtCKLs in various tissues; (2) AtCKL2 and AtCKL7 were expressed in early developmental anthers under high temperature (HT) conditions; (3) AtMYB24 may function as a positive regulator of AtCKL2 and AtCKL7 expression under HT. High temperature (HT) can seriously impede plant growth and development, causing severe loss of crop production. In Arabidopsis, AtCKL genes show high similarity to GhCKI, a gene reported to disrupt tapetal programmed cell death in cotton. However, most of AtCKL genes are not well characterized. Here, we systematically analyzed the expression patterns of AtCKLs in various tissues. The expression of AtCKL2 and AtCKL7 was induced in early anther development under HT, which is similar to the case of GhCKI. In silico analysis of AtCKL2 and AtCKL7 promoters indicated that four types of transcription factors (TFs) (MADS, NAC, WRKY and R2R3-MYB) might bind to AtCKL2 and AtCKL7 promoters. Furthermore, three MADS, three NAC, one WRKY, and three R2R3-MYB TFs were up-regulated in stage 1-8 anthers and three R2R3-MYB TFs were up-regulated in stage 9-10 anthers under HT, implying the important roles of R2R3-MYB genes in the response of anthers to HT. Among the R2R3-MYB genes, AtMYB24 showed the similar expression as AtCKL2 and AtCKL7 in the anthers under HT. Additionally, yeast one-hybrid and dual-luciferase reporter system assays verified that AtMYB24 could bind to AtCKL2 and AtCKL7 promoters and activate the expression of these two genes. In brief, this study provides the overall expression profiles of AtCKLs, useful information for unraveling the molecular mechanism of AtCKL2 and AtCKL7 gene expression in early anther development under HT, and important clues for elucidating the mechanism of transcriptional regulation of CKI genes in plant anther under HT, which are critical to the reduction of crop yield loss resulting from HT.

Functional Analysis of Cellulose Synthase CesA4 and CesA6 Genes in Switchgrass (Panicum virgatum) by Overexpression and RNAi-Mediated Gene Silencing

Switchgrass (Panicum virgatum L.) is a leading lignocellulosic bioenergy feedstock. Cellulose is a major component of the plant cell walls and the primary substrate for saccharification. Accessibility of cellulose to enzymatic breakdown into fermentable sugars is limited by the presence of lignin in the plant cell wall. In this study, putatively novel switchgrass secondary cell wall cellulose synthase PvCesA4 and primary cell wall PvCesA6 genes were identified and their functional role in cellulose synthesis and cell wall composition was examined by overexpression and knockdown of the individual genes in switchgrass. The endogenous expression of PvCesA4 and PvCesA6 genes varied among including roots, leaves, stem, and reproductive tissues. Increasing or decreasing PvCesA4 and PvCesA6 expression to extreme levels in the transgenic lines resulted in decreased biomass production. PvCesA6-overexpressing lines had reduced lignin content and syringyl/guaiacyl lignin monomer ratio accompanied by increased sugar release efficiency, suggesting an impact of PvCesA6 expression levels on lignin biosynthesis. Cellulose content and cellulose crystallinity were decreased, while xylan content was increased in PvCesA4 and PvCesA6 overexpression or knockdown lines. The increase in xylan content suggests that the amount of non-cellulosic cell wall polysaccharide was modified in these plants. Taken together, the results show that the manipulation of the cellulose synthase genes alters the cell wall composition and availability of cellulose as a bioprocessing substrate.

Genome-Wide Identification and Expression Profiling Analysis of the Xyloglucan Endotransglucosylase/Hydrolase Gene Family in Tobacco (Nicotiana tabacum L.)

Xyloglucan endotransglucosylase/hydrolase genes (XTHs) encode enzymes required for the reconstruction and modification of xyloglucan backbones, which will result in changes of cell wall extensibility during growth. A total of 56 NtXTH genes were identified from common tobacco, and 50 cDNA fragments were verified by PCR amplification. The 56 NtXTH genes could be classified into two subfamilies: Group I/II and Group III according to their phylogenetic relationships. The gene structure, chromosomal localization, conserved protein domains prediction, sub-cellular localization of NtXTH proteins and evolutionary relationships among Nicotiana tabacum, Nicotiana sylvestrisis, Nicotiana tomentosiformis, Arabidopsis, and rice were also analyzed. The NtXTHs expression profiles analyzed by the TobEA database and qRT-PCR revealed that NtXTHs display different expression patterns in different tissues. Notably, the expression patterns of 12 NtXTHs responding to environment stresses, including salinity, alkali, heat, chilling, and plant hormones, including IAA and brassinolide, were characterized. All the results would be useful for the function study of NtXTHs during different growth cycles and stresses.

Exploring the loblolly pine (Pinus taeda L.) genome by BAC sequencing and Cot analysis

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

Genome-wide identification, classification and expression analysis of the serine carboxypeptidase-like protein family in poplar

Previous studies have shown that the serine carboxypeptidase-like (SCPL) proteins in several plants play a key part in plant growth, development and stress responses. However, little is known about the functions of the SCPL genes in poplar. We identified 57 SCPL genes and divided into 3 subfamilies, which were unevenly distributed on 19 poplar chromosomes. Gene structure indicated that SCPL genes contain more introns, and motifs of each subfamily were relatively conserved. There were a total of 14 pairs of paralogs, with 6 pairs of these paralogs generated by segmental duplication and 1 generated by tandem duplication. In microsynteny analysis, large-scale duplication events played a key part in the expansion of Carboxypeptidase III genes. Expression of these genes was higher in mature leaf. Quantitative real-time PCR showed that majority of the SCPL genes were induced by methyl jasmonate (MeJA) treatment. PtSCPL27 and PtSCPL40 were located on the cytomembrane by conducting subcellular localization analysis. Our paper provides a theoretical basis for further functional research of PtSCPL genes and will benefit the molecular breeding for resistance to disease in poplar.