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

Genome-wide identification of the HSP70 genes in Pacific oyster Magallana gigas and their response to heat stress

Heat shock protein 70 (HSP70), the most prominent and well-characterized stress protein in animals, plays an important role in assisting animals in responding to various adverse conditions. In the present study, a total of 113 HSP70 gene family members were identified in the updated genome of Magallana gigas (designated MgHSP70) (previously known as Crassostrea gigas). There were 75, 12, 11, and 8 HSP70s located in the cytoplasm, nucleus, mitochondria, and endoplasmic reticulum, respectively, and 7 HSP70s were located in both the nucleus and cytoplasm. Among 113 MgHSP70 genes, 107 were unevenly distributed in 8 chromosomes of M. gigas with the greatest number in chromosome 07 (61 genes, 57.01%). The MgHSP70 gene family members were mainly assigned into five clusters, among which the HSPa12 subfamily underwent lineage-specific expansion, consisting of 89 members. A total of 68 MgHSP70 genes (60.18%) were tandemly duplicated and formed 30 gene pairs, among which 14 gene pairs were under strong positive selection. In general, the expression of MgHSP70s was tissue-specific, with the highest expression in labial palp and gill and the lowest expression in adductor muscle and hemocytes. There were 35, 31, and 47 significantly upregulated genes at 6, 12, and 24 h after heat shock treatment (28 °C), respectively. The expression patterns of different tandemly duplicated genes exhibited distinct characteristics after shock treatment, indicating that these genes may have different functions. Nevertheless, genes within the same tandemly duplicated group exhibit similar expression patterns. Most of the tandemly duplicated HSP70 gene pairs showed the highest expression levels at 24 h. This study provides a comprehensive description of the MgHSP70 gene family in M. gigas and offers valuable insights into the functions of HSP70 in the mollusc adaptation of oysters to environmental stress.

Characterization of the DUF868 gene family in Nicotiana and functional analysis of NtDUF868-E5 involved in pigment metabolism

Domains of unknown function (DUF) proteins represent a large group of uncharacterized protein families. The DUF868 gene family in Nicotiana has not yet been described. In the present study, we identified 12, 11, and 25 DUF868 family members in the genome of Nicotiana sylvestris, N. tomentosiformis, and N. tabacum, respectively. Based on phylogenetic analysis, these were categorized into five groups (A-E). Within each group, the gene structures, motifs, and tertiary structures showed high similarity. NtDUF868 family expansion during evolution was mainly driven by segmental duplication events. MicroRNA (miRNA) target site prediction identified 12 miRNA members that target 16 NtDUF868 family genes. The promoters of these genes contain cis-regulatory elements responsive to light, phytohormones, and abiotic stresses. Expression profiling revealed their tissue- and stage-specific expression patterns. RNA-sequencing and quantitative reverse transcription PCR revealed that the NtDUF868 family genes are potentially involved in the response to abiotic and biotic stresses, particularly drought and hormone stresses, and in the resistance to black shank and bacterial wilt. We generated transformed plants using NtDUF868-E5 overexpression and gene-editing vectors. NtDUF868-E5 overexpression resulted in enhanced tobacco plant growth and development, leading to increased leaf photosynthetic capacity and higher chlorophyll and carotenoid contents. This study provided a comprehensive genome-wide analysis of the DUF868 gene family, shedding light on their potential roles in plant growth and stress responses.

Genome-wide association studies for a comprehensive understanding of the genetic architecture of culm strength and yield traits in rice

Lodging resistance in rice is a complex trait determined by culm morphological and culm physical strength traits, and these traits are a major determinant of yield. We made a detailed analysis of various component traits with the aim of deriving optimized parameters for measuring culm strength. Genotyping by sequencing (GBS)-based genome-wide association study (GWAS) was employed among 181 genotypes for dissecting the genetic control of culm strength traits. The VanRaden kinship algorithm using 6,822 filtered single-nucleotide polymorphisms (SNPs) revealed the presence of two sub-groups within the association panel with kinship values concentrated at<0.5 level, indicating greater diversity among the genotypes. A wide range of phenotypic variation and high heritability for culm strength and yield traits were observed over two seasons, as reflected in best linear unbiased prediction (BLUP) estimates. The multi-locus model for GWAS resulted in the identification of 15 highly significant associations (p< 0.0001) for culm strength traits. Two novel major effect marker-trait associations (MTAs) for section modulus and bending stress were identified on chromosomes 2 and 12 with a phenotypic variance of 21.87% and 10.14%, respectively. Other MTAs were also noted in the vicinity of previously reported putative candidate genes for lodging resistance, providing an opportunity for further research on the biochemical basis of culm strength. The quantitative trait locus (QTL) hotspot identified on chromosome 12 with the synergistic association for culm strength trait (section modulus, bending stress, and internode breaking weight) and grain number can be considered a novel genomic region that can serve a dual purpose of enhancing culm strength and grain yield. Elite donors in the indica background with beneficial alleles of the identified major QTLs could be a valuable resource with greater significance in practical plant breeding programs focusing on improving lodging resistance in rice.

Genome-wide molecular evolution analysis of the GRF and GIF gene families in Plantae (Archaeplastida)

BACKGROUND

Plant growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) interact with each other and collectively have important regulatory roles in plant growth, development, and stress responses. Therefore, it is of great significance to explore the systematic evolution of GRF and GIF gene families. However, our knowledge and understanding of the role of GRF and GIF genes during plant evolution has been fragmentary.

RESULTS

In this study, a large number of genomic and transcriptomic datasets of algae, mosses, ferns, gymnosperms and angiosperms were used to systematically analyze the evolution of GRF and GIF genes during the evolution of plants. The results showed that GRF gene first appeared in the charophyte Klebsormidium nitens, whereas the GIF genes originated relatively early, and these two gene families were mainly expanded by segmental duplication events after plant terrestrialization. During the process of evolution, the protein sequences and functions of GRF and GIF family genes are relatively conservative. As cooperative partner, GRF and GIF genes contain the similar types of cis-acting elements in their promoter regions, which enables them to have similar transcriptional response patterns, and both show higher levels of expression in reproductive organs and tissues and organs with strong capacity for cell division. Based on protein-protein interaction analysis and verification, we found that the GRF-GIF protein partnership began to be established in pteridophytes and is highly conserved across different terrestrial plants.

CONCLUSIONS

These results provide a foundation for further exploration of the molecular evolution and biological functions of GRF and GIF genes.

Cell- and development-specific degradation controls the levels of mixed-linkage glucan in sorghum leaves

Mixed-linkage glucan (MLG) is a component of the cell wall (CW) of grasses and is composed of glucose monomers linked by β-1,3 and β-1,4 bonds. MLG is believed to have several biological functions, such as the mobilizable storage of carbohydrates and structural support of the CW. The extracellular levels of MLG are largely controlled by rates of synthesis mediated by cellulose synthase-like (CSL) enzymes, and turnover by lichenases. Economically important crops like sorghum accumulate MLG to variable levels during development. While in sorghum, like other grasses, there is one major MLG synthase (CSLF6), the identity of lichenases is yet unknown. To fill this gap, we identified three sorghum lichenases (SbLCH1-3) and characterized them in leaves in relation to the expression of SbCSLF6, and the abundance of MLG and starch. We established that SbLCH1-3 are secreted to the apoplast, consistent with a role of degrading MLG extracellularly. Furthermore, while SbCSLF6 expression was associated with cell development, the SbLCH genes exhibited distinct development-, cell-type-specific and diel-regulated expression. Therefore, our study identifies three functional sorghum MLG lichenases and highlights that MLG accumulation in sorghum leaves is likely controlled by the activity of lichenases that tune MLG levels, possibly to suit distinct cell and developmental needs in planta. These findings have important implications for improving the growth, yield, and composition of sorghum as a feedstock.

Genome-wide analysis of the heat shock transcription factor family reveals saline-alkali stress responses in Xanthoceras sorbifolium

The heat shock transcription factor (HSF) family is involved in regulating growth, development, and abiotic stress. The characteristics and biological functions of HSF family member in X. sorbifolium, an important oil and ornamental plant, have never been reported. In this study, 21 XsHSF genes were identified from the genome of X. sorbifolium and named XsHSF1-XsHSF21 based on their chromosomal positions. Those genes were divided into three groups, A, B, and C, containing 12, one, and eight genes, respectively. Among them, 20 XsHSF genes are located on 11 chromosomes. Protein structure analysis suggested that XsHSF proteins were conserved, displaying typical DNA binding domains (DBD) and oligomerization domains (OD). Moreover, HSF proteins within the same group contain specific motifs, such as motif 5 in the HSFC group. All XsHSF genes have one intron in the CDS region, except XsHSF1 which has two introns. Promoter analysis revealed that in addition to defense and stress responsiveness elements, some promoters also contained a MYB binding site and elements involved in multiple hormones responsiveness and anaerobic induction. Duplication analysis revealed that XsHSF1 and XsHSF4 genes were segmentally duplicated while XsHSF2, XsHSF9, and XsHSF13 genes might have arisen from transposition. Expression pattern analysis of leaves and roots following salt-alkali treatment using qRT-PCR indicated that five XsHSF genes were upregulated and one XsHSF gene was downregulated in leaves upon NaCl treatment suggesting these genes may play important roles in salt response. Additionally, the expression levels of most XsHSFs were decreased in leaves and roots following alkali-induced stress, indicating that those XsHSFs may function as negative regulators in alkali tolerance. MicroRNA target site prediction indicated that 16 of the XsHSF genes may be regulated by multiple microRNAs, for example XsHSF2 might be regulated by miR156, miR394, miR395, miR408, miR7129, and miR854. And miR164 may effect the mRNA levels of XsHSF3 and XsHSF17, XsHSF9 gene may be regulated by miR172. The expression trends of miR172 and miR164 in leaves and roots on salt treatments were opposite to the expression trend of XsHSF9 and XsHSF3 genes, respectively. Promoter analysis showed that XsHSFs might be involved in light and hormone responses, plant development, as well as abiotic stress responses. Our results thus provide an overview of the HSF family in X. sorbifolium and lay a foundation for future functional studies to reveal its roles in saline-alkali response.

Genome-Wide Identification and Functional Characterization of FAR1-RELATED SEQUENCE (FRS) Family Members in Potato (Solanum tuberosum)

FAR1-RELATED SEQUENCE (FRS) transcription factors are generated by transposases and play vital roles in plant growth and development, light signaling transduction, phytohormone response, and stress resistance. FRSs have been described in various plant species. However, FRS family members and their functions remain poorly understood in vegetative crops such as potato (Solanum tuberosum, St). In the present study, 20 putative StFRS proteins were identified in potato via genome-wide analysis. They were non-randomly localized to eight chromosomes and phylogenetic analysis classified them into six subgroups along with FRS proteins from Arabidopsis and tomato. Conserved protein motif, protein domain, and gene structure analyses supported the evolutionary relationships among the FRS proteins. Analysis of the cis-acting elements in the promoters and the expression profiles of StFRSs in various plant tissues and under different stress treatments revealed the spatiotemporal expression patterns and the potential roles of StFRSs in phytohormonal and stress responses. StFRSs were differentially expressed in the cultivar "Xisen 6", which is exposed to a variety of stresses. Hence, these genes may be critical in regulating abiotic stress. Elucidating the StFRS functions will lay theoretical and empirical foundations for the molecular breeding of potato varieties with high light use efficiency and stress resistance.

Genome-wide identification, characterization and transcriptional profile of the SWEET gene family in Dendrobium officinale

BACKGROUND

Dendrobium officinale Kimura et Migo (D. officinale) is a well-known traditional Chinese medicine with high content polysaccharides in stems. The SWEET (Sugars Will Eventually be Exported Transporters) family is a novel class of sugar transporters mediating sugar translocation among adjacent cells of plants. The expression patterns of SWEETs and whether they are associated with stress response in D. officinale remains uncovered.

RESULTS

Here, 25 SWEET genes were screened out from D. officinale genome, most of which typically contained seven transmembrane domains (TMs) and harbored two conserved MtN3/saliva domains. Using multi-omics data and bioinformatic approaches, the evolutionary relationship, conserved motifs, chromosomal location, expression patterns, correlationship and interaction network were further analyzed. DoSWEETs were intensively located in nine chromosomes. Phylogenetic analysis revealed that DoSWEETs were divided into four clades, and conserved motif 3 specifically existed in DoSWEETs from clade II. Different tissue-specific expression patterns of DoSWEETs suggested the division of their roles in sugar transport. In particular, DoSWEET5b, 5c, and 7d displayed relatively high expression levels in stems. DoSWEET2b and 16 were significantly regulated under cold, drought, and MeJA treatment, which were further verified using RT-qPCR. Correlation analysis and interaction network prediction discovered the internal relationship of DoSWEET family.

CONCLUSIONS

Taken together, the identification and analysis of the 25 DoSWEETs in this study provide basic information for further functional verification in D. officinale.

Genome-Wide Analysis and Expression Profiling of the Glutathione Peroxidase-like Enzyme Gene Family in Solanum tuberosum

Glutathione peroxidase-like enzyme is an important enzymatic antioxidant in plants. It is involved in scavenging reactive oxygen species, which can effectively prevent oxidative damage and improve resistance. GPXL has been studied in many plants but has not been reported in potatoes, the world's fourth-largest food crop. This study identified eight StGPXL genes in potatoes for the first time through genome-wide bioinformatics analysis and further studied the expression patterns of these genes using qRT-PCR. The results showed that the expression of StGPXL1 was significantly upregulated under high-temperature stress, indicating its involvement in potato defense against high-temperature stress, while the expression levels of StGPXL4 and StGPXL5 were significantly downregulated. The expression of StGPXL1, StGPXL2, StGPXL3, and StGPXL6 was significantly upregulated under drought stress, indicating their involvement in potato defense against drought stress. After MeJA hormone treatment, the expression level of StGPXL6 was significantly upregulated, indicating its involvement in the chemical defense mechanism of potatoes. The expression of all StGPXL genes is inhibited under biotic stress, which indicates that GPXL is a multifunctional gene family, which may endow plants with resistance to various stresses. This study will help deepen the understanding of the function of the potato GPXL gene family, provide comprehensive information for the further analysis of the molecular function of the potato GPXL gene family as well as a theoretical basis for potato molecular breeding.

Characterization of HSP70 family in watermelon (Citrullus lanatus): identification, structure, evolution, and potential function in response to ABA, cold and drought stress

Watermelon (Citrullus lanatus) as a crop with important economic value, is widely cultivated around the world. The heat shock protein 70 (HSP70) family in plant is indispensable under stress conditions. However, no comprehensive analysis of watermelon HSP70 family is reported to date. In this study, 12 ClHSP70 genes were identified from watermelon, which were unevenly located in 7 out of 11 chromosomes and divided into three subfamilies. ClHSP70 proteins were predicted to be localized primarily in cytoplasm, chloroplast, and endoplasmic reticulum. Two pairs of segmental repeats and 1 pair of tandem repeats existed in ClHSP70 genes, and ClHSP70s underwent strong purification selection. There were many abscisic acid (ABA) and abiotic stress response elements in ClHSP70 promoters. Additionally, the transcriptional levels of ClHSP70s in roots, stems, true leaves, and cotyledons were also analyzed. Some of ClHSP70 genes were also strongly induced by ABA. Furthermore, ClHSP70s also had different degrees of response to drought and cold stress. The above data indicate that ClHSP70s may be participated in growth and development, signal transduction and abiotic stress response, laying a foundation for further analysis of the function of ClHSP70s in biological processes.

Genome-wide investigation and expression profiling of LOR gene family in rapeseed under salinity and ABA stress

The Brassica napus (B. napus) LOR (Lurp-One-Related) gene family is a little-known gene family characterized by a conserved LOR domain in the proteins. Limited research in Arabidopsis showed that LOR family members played important roles in Hyaloperonospora parasitica (Hpa) defense. Nevertheless, there is a paucity of research investigating the role of the LOR gene family towards their responses to abiotic stresses and hormone treatments. This study encompassed a comprehensive survey of 56 LOR genes in B. napus, which is a prominent oilseed crop that holds substantial economic significance in China, Europe, and North America. Additionally, the study evaluated the expression profiles of these genes in response to salinity and ABA stress. Phylogenetic analysis showed that 56 BnLORs could be divided into 3 subgroups (8 clades) with uneven distribution on 19 chromosomes. 37 out of 56 BnLOR members have experienced segmental duplication and 5 of them have undergone tandem repeats events with strong evidence of purifying selection. Cis-regulatory elements (CREs) analysis indicated that BnLORs involved in process such as light response, hormone response, low temperature response, heat stress response, and dehydration response. The expression pattern of BnLOR family members revealed tissue specificity. RNA-Seq and qRT-PCR were used to validate BnLOR gene expression under temperature, salinity and ABA stress, revealing that most BnLORs showed inducibility. This study enhanced our comprehension of the B. napus LOR gene family and could provide valuable information for identifying and selecting genes for stress resistant breeding.

Identification of ZmBK2 Gene Variation Involved in Regulating Maize Brittleness

Maize stalk strength is a crucial agronomic trait that affects lodging resistance. We used map-based cloning and allelic tests to identify a maize mutant associated with decreased stalk strength and confirmed that the mutated gene, ZmBK2, is a homolog of Arabidopsis AtCOBL4, which encodes a COBRA-like glycosylphosphatidylinositol (GPI)-anchored protein. The bk2 mutant exhibited lower cellulose content and whole-plant brittleness. Microscopic observations showed that sclerenchymatous cells were reduced in number and had thinner cell walls, suggesting that ZmBK2 affects the development of cell walls. Transcriptome sequencing of differentially expressed genes in the leaves and stalks revealed substantial changes in the genes associated with cell wall development. We constructed a cell wall regulatory network using these differentially expressed genes, which revealed that abnormal cellulose synthesis may be a reason for brittleness. These results reinforce our understanding of cell wall development and provide a foundation for studying the mechanisms underlying maize lodging resistance.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Single-molecular insights into the breakpoint of cellulose nanofibers assembly during saccharification

Plant cellulose microfibrils are increasingly employed to produce functional nanofibers and nanocrystals for biomaterials, but their catalytic formation and conversion mechanisms remain elusive. Here, we characterize length-reduced cellulose nanofibers assembly in situ accounting for the high density of amorphous cellulose regions in the natural rice fragile culm 16 (Osfc16) mutant defective in cellulose biosynthesis using both classic and advanced atomic force microscopy (AFM) techniques equipped with a single-molecular recognition system. By employing individual types of cellulases, we observe efficient enzymatic catalysis modes in the mutant, due to amorphous and inner-broken cellulose chains elevated as breakpoints for initiating and completing cellulose hydrolyses into higher-yield fermentable sugars. Furthermore, effective chemical catalysis mode is examined in vitro for cellulose nanofibers conversion into nanocrystals with reduced dimensions. Our study addresses how plant cellulose substrates are digestible and convertible, revealing a strategy for precise engineering of cellulose substrates toward cost-effective biofuels and high-quality bioproducts.

Genome-wide identification and analysis of the IQM gene family in soybean

IQM, a plant-specific calmodulin-binding protein, plays multiple roles in plant growth and development. Although a comprehensive analysis has been carried out on the IQM family genes in Arabidopsis and rice, the number and functions of IQM genes in other species have not been explored. In this study, we identified 15 members of the soybean (Glycine max) IQM gene family using BLASTP tools. These members were distributed on 12 soybean chromosomes and constitute six pairs caused by fragment duplication events. According to phylogeny, the 15 genes were divided into three subfamilies (I, II, and III), and members of the same subfamily had similar gene and protein structures. Yeast two-hybrid experiments revealed that the IQ motif is critical for the binding of GmIQM proteins to GmCaM, and its function is conserved in soybean, Arabidopsis, and rice. Based on real-time PCR, the soybean IQM genes were strongly induced by PEG and NaCl, suggesting their important biological functions in abiotic stress responses. Overall, this genome-wide analysis of the soybean IQM gene family lays a solid theoretical foundation for further research on the functions of GmIQM genes and could serve as a reference for the improvement and breeding of soybean stress resistance traits.

Identification of Alfalfa SPL gene family and expression analysis under biotic and abiotic stresses

The SQUAMOSA promoter binding-like protein (SPL) is a specific transcription factor that affects plant growth and development. The SPL gene family has been explored in various plants, but information about these genes in alfalfa is limited. This study, based on the whole genome data of alfalfa SPL, the fundamental physicochemical properties, phylogenetic evolution, gene structure, cis-acting elements, and gene expression of members of the MsSPL gene family were analyzed by bioinformatics methods. We identified 82 SPL sequences in the alfalfa, which were annotated into 23 genes, including 7 (30.43%) genes with four alleles, 10 (43.47%) with three, 3 (13.04%) with two, 3 (13.04%) with one allele. These SPL genes were divided into six groups, that are constructed from A. thaliana, M. truncatula and alfalfa. Chromosomal localization of the identified SPL genes showed arbitary distribution. The subcellular localization predictions showed that all MsSPL proteins were located in the nucleus. A total of 71 pairs of duplicated genes were identified, and segmental duplication mainly contributed to the expansion of the MsSPL gene family. Analysis of the Ka/Ks ratios indicated that paralogs of the MsSPL gene family principally underwent purifying selection. Protein-protein interaction analysis of MsSPL proteins were performed to predict their roles in potential regulatory networks. Twelve cis-acting elements including phytohormone and stress elements were detected in the regions of MsSPL genes. We further analyzed that the MsSPLs had apparent responses to abiotic stresses such as drought and salt and the biotic stress of methyl jasmonate. These results provide comprehensive information on the MsSPL gene family in alfalfa and lay a solid foundation for elucidating the biological functions of MsSPLs. This study also provides valuable on the regulation mechanism and function of MsSPLs in response to biotic and abiotic stresses.

OsCSLD1 Mediates NH4+-Dependent Root Hair Growth Suppression and AMT1;2 Expression in Rice (Oryza sativa L.)

Root hairs play crucial roles in the roots, including nutrient uptake, water assimilation, and anchorage with soil, along with supporting rhizospheric microorganisms. In rice, ammonia uptake is mediated by a specialized ammonium transporter (AMT). AMT1;1, AMT1;2, and AMT1;3 have been extensively studied in relation to nitrogen signaling. Cellulose synthase-like D1 (CSLD1) is essential for cell expansion and is highly specific to root hair cells. csld1 mutants showed successful initiation but failed to elongate. However, when nitrogen was depleted, csld1 root hairs resumed elongation. Further experiments revealed that in the presence of ammonium (NH4+), csld1 roots failed to elongate. csld1 elongated normally in the presence of nitrate (NO3−). Expression analysis showed an increase in root hair-specific AMT1;2 expression in csld1. CSLD1 was positively co-expressed with AMT1;2 changing nitrogen concentration in the growth media. CSLD1 showed increased expression in the presence of both ammonium and nitrate. Methylammonium (MeA) treatment of CSLD1 overexpression lines suggests that CSLD1 does not directly participate in nitrogen transport. Further studies on the root hair elongation mutant sndp1 showed that nitrogen assimilation is unlikely to depend on root hair length. Therefore, these results suggest that CSLD1 is closely involved in nitrogen-dependent root hair elongation and regulation of AMT1;2 expression in rice roots.

Genome-wide identification, phylogenetic, and expression analysis under abiotic stress conditions of Whirly (WHY) gene family in Medicago sativa L

The WHY family is a group of plant-specific transcription factors, that can bind to single-stranded DNA molecules and play a variety of functions in plant nuclei and organelles, participating in the regulation of plant leaf senescence. It has been identified and analyzed in many species, however, the systematic identification and analysis of the WHY genes family have not yet been reported in alfalfa (Medicago sativa L.). Therefore, to explore the function of alfalfa the WHY genes, and 10 MsWHY genes were identified and further characterized their evolutionary relationship and expression patterns by analyzing the recently published genome of alfalfa. Comprehensive analysis of the chromosome location, physicochemical properties of the protein, evolutionary relationship, conserved motifs, and responses to abiotic stresses of the WHY gene family in alfalfa using bioinformatics methods. The results showed that 10 MsWHY genes were distributed on 10 chromosomes, and collinearity analysis showed that many MsWHYs might be derived from segmental duplications, and these genes are under purifying selection. Based on phylogenetic analyses, the WHY gene family of alfalfa can be divided into four subfamilies: I-IV subfamily, and approximately all the WHY genes within the same subfamily share similar gene structures. The 10 MsWHY gene family members contained 10 motifs, of which motif 2 and motif 4 are the conserved motifs shared by these genes. Furthermore, the analysis of cis-regulatory elements indicated that regulatory elements related to transcription, cell cycle, development, hormone, and stress response are abundant in the promoter sequence of the MsWHY genes. Real-time quantitative PCR demonstrated that MsWHYs gene expression is induced by drought, salt, and methyl jasmonate. The present study serves as a basic foundation for future functional studies on the alfalfa WHY family.

Genome-wide identification and adaptive evolution of CesA/Csl superfamily among species with different life forms in Orchidaceae

Orchidaceae, with more than 25,000 species, is one of the largest flowering plant families that can successfully colonize wide ecological niches, such as land, trees, or rocks, and its members are divided into epiphytic, terrestrial, and saprophytic types according to their life forms. Cellulose synthase (CesA) and cellulose synthase-like (Csl) genes are key regulators in the synthesis of plant cell wall polysaccharides, which play an important role in the adaptation of orchids to resist abiotic stresses, such as drought and cold. In this study, nine whole-genome sequenced orchid species with three types of life forms were selected; the CesA/Csl gene family was identified; the evolutionary roles and expression patterns of CesA/Csl genes adapted to different life forms and abiotic stresses were investigated. The CesA/Csl genes of nine orchid species were divided into eight subfamilies: CesA and CslA/B/C/D/E/G/H, among which the CslD subfamily had the highest number of genes, followed by CesA, whereas CslB subfamily had the least number of genes. Expansion of the CesA/Csl gene family in orchids mainly occurred in the CslD and CslF subfamilies. Conserved domain analysis revealed that eight subfamilies were conserved with variations in orchids. In total, 17 pairs of CesA/Csl homologous genes underwent positive selection, of which 86%, 14%, and none belonged to the epiphytic, terrestrial, and saprophytic orchids, respectively. The inter-species collinearity analysis showed that the CslD genes expanded in epiphytic orchids. Compared with terrestrial and saprophytic orchids, epiphytic orchids experienced greater strength of positive selection, with expansion events mostly related to the CslD subfamily, which might have resulted in strong adaptability to stress in epiphytes. Experiments on stem expression changes under abiotic stress showed that the CslA might be a key subfamily in response to drought stress for orchids with different life forms, whereas the CslD might be a key subfamily in epiphytic and saprophytic orchids to adapt to freezing stress. This study provides the basic knowledge for the further systematic study of the adaptive evolution of the CesA/Csl superfamily in angiosperms with different life forms, and research on orchid-specific functional genes related to life-history trait evolution.

Rice cellulose synthase-like protein OsCSLD4 coordinates the trade-off between plant growth and defense

Plant cell wall is a complex and changeable structure, which is very important for plant growth and development. It is clear that cell wall polysaccharide synthases have critical functions in rice growth and abiotic stress, yet their role in plant response to pathogen invasion is poorly understood. Here, we describe a dwarf and narrowed leaf in Hejiang 19 (dnl19) mutant in rice, which shows multiple growth defects such as reduced plant height, enlarged lamina joint angle, curled leaf morphology, and a decrease in panicle length and seed setting. MutMap analysis, genetic complementation and gene knockout mutant show that cellulose synthase-like D4 (OsCSLD4) is the causal gene for DNL19. Loss function of OsCSLD4 leads to a constitutive activation of defense response in rice. After inoculation with rice blast and bacterial blight, dnl19 displays an enhanced disease resistance. Widely targeted metabolomics analysis reveals that disruption of OsCSLD4 in dnl19 resulted in significant increase of L-valine, L-asparagine, L-histidine, L-alanine, gentisic acid, but significant decrease of L-aspartic acid, malic acid, 6-phosphogluconic acid, glucose 6-phosphate, galactose 1-phosphate, gluconic acid, D-aspartic acid. Collectively, our data reveals the importance of OsCSLD4 in balancing the trade-off between rice growth and defense.

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.

The class III peroxidase gene family is involved in ascorbic acid induced delay of internal browning in pineapple

Excessive production of reactive oxygen species (ROS) leads to potential toxicity in an organism. Class III peroxidases (PRXs) play an important role in maintaining ROS homeostasis in plants. Internal browning (IB) limits industrial development of pineapple, which is the third most important fruit trade in the world. IB is mainly caused by ROS, and the mechanism underlying IB is still unknown from the perspective of ROS. Here, we soaked pineapples in ascorbic acid after harvest and before storage to decrease excessive ROS and polyphenol oxidase (PPO) activity, ultimately restraining the spread and deterioration of IB. Using phylogenetic analysis; we identified 78 pineapple PRX genes (AcPRXs) and divided them into five subgroups. Gene structure analysis indicated that the exon numbers ranged from 2 to 14, and conserved motif analysis verified that all of the AcPRXs identified here have standard peroxidase domains. Analysis of duplication events suggested that tandem and segmental duplication events may have played equal and important roles in expanding the AcPRX family. Comprehensive transcriptomic analysis uncovered that AcPRXs may play an important role in negatively regulating the occurrence of IB. In summary, we found that ROS scavenging delayed IB occurrence. The results of characterized AcPRX family revealed that AcPRXs family responded to growth and development, and negatively regulated to IB occurrence in storage stage. This research provides potential target genes for future in-depth analysis of the molecular mechanisms underlying IB and contributes to develop IB-resistant pineapple varieties.

Comprehensive identification, evolutionary patterns and the divergent response of PRX genes in Phaseolus vulgaris under biotic and abiotic interactions

Peroxiredoxins (Prxs) are novel cysteine-based peroxidases which are involved in protecting cells from oxidative damage by catalyzing the reduction of different peroxides. The present study addressed, for the first time, genome-wide identification, evolutionary patterns and expression dynamics of Phaseolus vulgaris Prx gene family (PvPrx). Nine Prx proteins were identified in P. vulgaris based on homology searches. The phylogeny analysis of Prxs from seven plant species revealed that Prx proteins can be clustered into four groups (1C-Prx, 2C-Prxs, PrxQ and type II Prxs). Both tandem and segmental duplication contributed to PvPrx gene family expansion. Intragenic reorganizations including gain/loss of exon/intron and insertions/deletions have also contributed to PvPrx gene diversification. The collinearity analysis revealed the presence of some orthologous Prx gene pairs between A. thaliana and P. vulgaris genomes. The Ka/Ks ratio indicated that two of the three PvPrx duplicated gene pairs have undergone a purifying selection. Redundant stress-related cis-acting elements were also found in the promoters of most PvPrx genes. RT q-PCR analysis revealed an upregulation of key PvPrx members in response to symbiosis and different abiotic factors. The upregulation of targeted PvPrx members, particularly in leaves exposed to salinity or drought, was accompanied by an accumulation of hydrogen peroxide (H₂O₂). When exogenously applied, H₂O₂ modulated almost all PvPrx genes, suggesting a potential H₂O₂-scavenging role for these proteins. Collectively, our analysis provided valuable information for further functional analysis of key PvPrx members to improve common bean stress tolerance and/or its symbiotic performance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03246-8.

Genome-Wide Analysis and Expression Profiling of the Phenylalanine Ammonia-Lyase Gene Family in Solanum tuberosum

Phenylalanine ammonia-lyase is one of the most widely studied enzymes in the plant kingdom. It is a crucial pathway from primary metabolism to significant secondary phenylpropanoid metabolism in plants, and plays an essential role in plant growth, development, and stress defense. Although PAL has been studied in many actual plants, only one report has been reported on potato, one of the five primary staple foods in the world. In this study, 14 StPAL genes were identified in potato for the first time using a genome-wide bioinformatics analysis, and the expression patterns of these genes were further investigated using qRT-PCR. The results showed that the expressions of StPAL1, StPAL6, StPAL8, StPAL12, and StPAL13 were significantly up-regulated under drought and high temperature stress, indicating that they may be involved in the stress defense of potato against high temperature and drought. The expressions of StPAL1, StPAL2, and StPAL6 were significantly up-regulated after MeJa hormone treatment, indicating that these genes are involved in potato chemical defense mechanisms. These three stresses significantly inhibited the expression of StPAL7, StPAL10, and StPAL11, again proving that PAL is a multifunctional gene family, which may give plants resistance to multiple and different stresses. In the future, people may improve critical agronomic traits of crops by introducing other PAL genes. This study aims to deepen the understanding of the versatility of the PAL gene family and provide a valuable reference for further genetic improvement of the potato.

Genome-wide identification and expression analysis of the R2R3-MYB gene family in tobacco (Nicotiana tabacum L.)

BACKGROUND

The R2R3-MYB transcription factor is one of the largest gene families in plants and involved in the regulation of plant development, hormone signal transduction, biotic and abiotic stresses. Tobacco is one of the most important model plants. Therefore, it will be of great significance to investigate the R2R3-MYB gene family and their expression patterns under abiotic stress and senescence in tobacco.

RESULTS

A total of 174 R2R3-MYB genes were identified from tobacco (Nicotiana tabacum L.) genome and were divided into 24 subgroups based on phylogenetic analysis. Gene structure (exon/intron) and protein motifs were especially conserved among the NtR2R3-MYB genes, especially members within the same subgroup. The NtR2R3-MYB genes were distributed on 24 tobacco chromosomes. Analysis of gene duplication events obtained 3 pairs of tandem duplication genes and 62 pairs of segmental duplication genes, suggesting that segmental duplications is the major pattern for R2R3-MYB gene family expansion in tobacco. Cis-regulatory elements of the NtR2R3-MYB promoters were involved in cellular development, phytohormones, environmental stress and photoresponsive. Expression profile analysis showed that NtR2R3-MYB genes were widely expressed in different maturity tobacco leaves, and however, the expression patterns of different members appeared to be diverse. The qRT-PCR analysis of 15 NtR2R3-MYBs confirmed their differential expression under different abiotic stresses (cold, salt and drought), and notably, NtMYB46 was significantly up-regulated under three treatments.

CONCLUSIONS

In summary, a genome-wide identification, evolutionary and expression analysis of R2R3-MYB gene family in tobacco were conducted. Our results provided a solid foundation for further biological functional study of NtR2R3-MYB genes in tobacco.

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.

Transcriptome Mining Provides Insights into Cell Wall Metabolism and Fiber Lignification in Agave tequilana Weber

Resilience of growing in arid and semiarid regions and a high capacity of accumulating sugar-rich biomass with low lignin percentages have placed Agave species as an emerging bioenergy crop. Although transcriptome sequencing of fiber-producing agave species has been explored, molecular bases that control wall cell biogenesis and metabolism in agave species are still poorly understood. Here, through RNAseq data mining, we reconstructed the cellulose biosynthesis pathway and the phenylpropanoid route producing lignin monomers in A. tequilana, and evaluated their expression patterns in silico and experimentally. Most of the orthologs retrieved showed differential expression levels when they were analyzed in different tissues with contrasting cellulose and lignin accumulation. Phylogenetic and structural motif analyses of putative CESA and CAD proteins allowed to identify those potentially involved with secondary cell wall formation. RT-qPCR assays revealed enhanced expression levels of AtqCAD5 and AtqCESA7 in parenchyma cells associated with extraxylary fibers, suggesting a mechanism of formation of sclerenchyma fibers in Agave similar to that reported for xylem cells in model eudicots. Overall, our results provide a framework for understanding molecular bases underlying cell wall biogenesis in Agave species studying mechanisms involving in leaf fiber development in monocots.

Rice STOMATAL CYTOKINESIS DEFECTIVE2 regulates cell expansion by affecting vesicular trafficking in rice

Vesicular trafficking plays critical roles in cell expansion in yeast and mammals, but information linking vesicular trafficking and cell expansion in plants is limited. Here, we isolated and characterized a rice (Oryza sativa) mutant, decreased plant height 1-1 (dph1-1), which exhibited a wide spectrum of developmental phenotypes, including reduced plant height and smaller panicles and grains. Cytological analysis revealed that limited cell expansion was responsible for the dph1-1 mutant phenotype compared to the wild-type. Map-based cloning revealed that DPH1 encodes a plant-specific protein, OsSCD2, which is homologous to Arabidopsis (Arabidopsis thaliana) STOMATAL CYTOKINESIS DEFECTIVE2 (SCD2). Subcellular localization revealed that OsSCD2 is associated with clathrin. Confocal microscopy showed that the dph1-1 mutant has defective endocytosis and post-Golgi trafficking. Biochemical and confocal data indicated that OsSCD2 physically interacts with OsSCD1 and that they are associated with intracellular structures that colocalize with microtubules. Furthermore, we found that cellulose synthesis was affected in the dph1-1 mutant, evidenced by reduced cellulose synthase gene accumulation at the transcript and protein levels, most likely resulting from an impaired localization pattern. Our results suggest that OsSCD2 is involved in clathrin-related vesicular trafficking with an important role in maintaining plant growth in rice.

Overexpression of PnMYB2 from Panax notoginseng induces cellulose and lignin biosynthesis during cell wall formation

Panax notoginseng PnMYB2 is a transcriptional activator of primary and secondary cell wall formation by promoting the PCW-specific gene CesA3 and key lignin biosynthetic gene CCoAOMT1, respectively. R2R3-MYB transcription factors play important roles in regulation secondary cell wall (SCW) formation. However, there are few reports on the functions of MYB transcription factors which involved in both primary cell wall (PCW) and SCW formation. Here, we isolated an R2R3-MYB transcription factor, PnMYB2, from Panax notoginseng roots which are widely used in Chinese traditional medicines and contain abundant cellulose and lignin. The expression pattern of PnMYB2 was similar to the accumulation pattern of cellulose and lignin contents in different organs. PnMYB2 localized in the nucleus and may function as a transcriptional activator. Overexpression of PnMYB2 in Arabidopsis thaliana enhanced cellulose and lignin biosynthesis, and remarkably increased thickness of PCW and SCW in the stem of transgenic plants compared with wild-type plants. The expression levels of genes associated with PCW-specific cellulose synthase (CesA) genes and key SCW-specific lignin biosynthetic genes were significantly increased in PnMYB2-overexpressing plants compared to the wild type plants. Furthermore, yeast one-hybrid, dual-luciferase reporter assays and electrophoretic mobility shift assays (EMSA) results verified that PnMYB2 could bind and activate the promoters of AtCesA3 and PnCesA3, which are the PCW-specific cellulose biosynthetic genes, and AtCCoAOMT1 and PnCCoAOMT1, which are the key lignin biosynthetic genes. These results demonstrated the central role of PnMYB2 in PCW-specific cellulose formation and SCW-specific lignin biosynthesis.

Genomic Architecture and Evolution of the Cellulose synthase Gene Superfamily as Revealed by Phylogenomic Analysis

The Cellulose synthase superfamily synthesizes cellulose and different hemicellulosic polysaccharides in plant cell walls. While much has been discovered about the evolution and function of these genes, their genomic architecture and relationship with gene (sub-)functionalization and evolution remains unclear. By using 242 genomes covering plant evolution from green algae to eudicots, we performed a large-scale analysis of synteny, phylogenetic, and functional data of the CesA superfamily. Results revealed considerable gene copy number variation across species and gene families, and also two patterns - singletons vs. tandem arrays - in chromosomic gene arrangement. Synteny analysis revealed exceptional conservation of gene architecture across species, but also lineage-specific patterns across gene (sub-)families. Synteny patterns correlated with gene sub-functionalization into primary and secondary CesAs and distinct CslD functional isoforms. Furthermore, a genomic context shift of a group of cotton secondary CesAs was associated with peculiar properties of cotton fiber synthesis. Finally, phylogenetics suggested that primary CesA sequences appeared before the secondary CesAs, while phylogenomic analyses unveiled the genomic trace of the CslD duplication that initiated the CslF family. Our results describe in detail the genomic architecture of the CesA superfamily in plants, highlighting its crucial relevance for gene diversification and sub-functionalization, and for understanding their evolution.

Characterization of pectin methylesterase gene family and its possible role in juice sac granulation in navel orange (Citrus sinensis Osbeck)

BACKGROUND

Citrus is one of the most important fresh fruit crops worldwide. Juice sac granulation is a physiological disorder, which leads to a reduction in soluble solid concentration, total sugar, and titratable acidity of citrus fruits. Pectin methylesterase (PME) catalyzes the de-methylesterification of homogalacturonans and plays crucial roles in cell wall modification during plant development and fruit ripening. Although PME family has been well investigated in various model plants, little is known regarding the evolutionary property and biological function of PME family genes in citrus.

RESULTS

In this study, 53 non-redundant PME genes were identified from Citrus sinensis genome, and these PME genes were divided into four clades based on the phylogenetic relationship. Subsequently, bioinformatics analyses of gene structure, conserved domain, chromosome localization, gene duplication, and collinearity were performed on CsPME genes, providing important clues for further research on the functions of CsPME genes. The expression profiles of CsPME genes in response to juice sac granulation and low-temperature stress revealed that CsPME genes were involved in the low temperature-induced juice sac granulation in navel orange fruits. Subcellular localization analysis suggested that CsPME genes were localized on the apoplast, endoplasmic reticulum, plasma membrane, and vacuole membrane. Moreover, yeast one-hybrid screening and dual luciferase activity assay revealed that the transcription factor CsRVE1 directly bound to the promoter of CsPME3 and activated its activity.

CONCLUSION

In summary, this study conducts a comprehensive analysis of the PME gene family in citrus, and provides a novel insight into the biological functions and regulation patterns of CsPME genes during juice sac granulation of citrus.

Global N6-Methyladenosine Profiling Revealed the Tissue-Specific Epitranscriptomic Regulation of Rice Responses to Salt Stress

N⁶-methyladenosine (m⁶A) methylation represents a new layer of the epitranscriptomic regulation of plant development and growth. However, the effects of m⁶A on rice responses to environmental stimuli remain unclear. In this study, we performed a methylated-RNA immunoprecipitation sequencing analysis and compared the changes in m⁶A methylation and gene expression in rice under salt stress conditions. Salt stress significantly increased the m⁶A methylation in the shoots (p value < 0.05). Additionally, 2537 and 2304 differential m⁶A sites within 2134 and 1997 genes were identified in the shoots and roots, respectively, under salt stress and control conditions. These differential m⁶A sites were largely regulated in a tissue-specific manner. A unique set of genes encoding transcription factors, antioxidants, and auxin-responsive proteins had increased or decreased m⁶A methylation levels only in the shoots or roots under salt stress, implying m⁶A may mediate salt tolerance by regulating transcription, ROS homeostasis, and auxin signaling in a tissue-specific manner. Integrating analyses of m⁶A modifications and gene expression changes revealed that m⁶A changes regulate the expression of genes controlling plant growth, stress responses, and ion transport under saline conditions. These findings may help clarify the regulatory effects of m⁶A modifications on rice salt tolerance.

Systematic analysis of the non-specific lipid transfer protein gene family in Nicotiana tabacum reveal its potential roles in stress responses

Plant non-specific lipid transfer proteins (nsLTPs) are characterized by an eight-cysteine motif backbone stabilized by four disulfide bonds; these proteins can bind or transfer lipids. NsLTPs play important roles in plant growth and development, and in the responses to abiotic and biotic stresses. In this study, 50, 51, and 100 nsLTPs from Nicotiana sylvestris, N. tomentosiformis, and their descendant N. tabacum, respectively, were identified and classified into six types (I, II, IV, V, VII, and VIII). The phylogeny, gene structures, motifs, tertiary structures, gene duplications and expression patterns were systematically analyzed. The intron/exon patterns and the conserved motifs were highly similar among the same types of nsLTP genes. Purifying selection and segmental duplication dominated the expansion of the nsLTPs family during evolution. Cis-regulatory elements of the NtLTP promoters were involved in light responsiveness, abiotic stress, and phytohormone responsiveness. Expression pattern analysis using RNA-seq and qPCR revealed that NtLTP family genes exhibited tissue-specific expression patterns and they have potential roles in response to abiotic and biotic stresses, especially drought stress, and resistance to black shank and bacterial wilt. Furthermore, overexpression of NtLTPI.38 in tobacco increased drought tolerance by improving the antioxidant defense ability, through reducing O2•- and H2O2 accumulation and increasing the number of lateral roots. These results provide a comprehensive overview of this gene family and provide valuable insights for the functional characterization of nsLTP family genes.

Identification and Fine Mapping of the Recessive Gene BK-5, Which Affects Cell Wall Biosynthesis and Plant Brittleness in Maize

The cellulose of the plant cell wall indirectly affects the cell shape and straw stiffness of the plant. Here, the novel brittleness mutant brittle stalk-5 (bk-5) of the maize inbred line RP125 was characterized. We found that the mutant displayed brittleness of the stalk and even the whole plant, and that the brittleness phenotype existed during the whole growth period from germination to senescence. The compressive strength was reduced, the cell wall was thinner, and the cellulose content was decreased compared to that of the wild type. Genetic analysis and map-based cloning indicated that bk-5 was controlled by a single recessive nuclear gene and that it was located in a 90.2-Kb region on chromosome 3 that covers three open reading frames (ORFs). Sequence analysis revealed a single non-synonymous missense mutation, T-to-A, in the last exon of Zm00001d043477 (B73: version 4, named BK-5) that caused the 951th amino acid to go from leucine to histidine. BK-5 encodes a cellulose synthase catalytic subunit (CesA), which is involved with cellulose synthesis. We found that BK-5 was constitutively expressed in all tissues of the germinating stage and silking stage, and highly expressed in the leaf, auricula, and root of the silking stage and the 2-cm root and bud of the germinating stage. We found that BK-5 mainly localized to the Golgi apparatus, suggesting that the protein might move to the plasma membrane with the aid of Golgi in maize. According to RNA-seq data, bk-5 had more downregulated genes than upregulated genes, and many of the downregulated genes were enzymes and transcription factors related to cellulose, hemicellulose, and lignin biosynthesis of the secondary cell wall. The other differentially expressed genes were related to metabolic and cellular processes, and were significantly enriched in hormone signal transduction, starch and sucrose metabolism, and the plant-pathogen interaction pathway. Taken together, we propose that the mutation of gene BK-5 causes the brittle stalk phenotype and provides important insights into the regulatory mechanism of cellulose biosynthesis and cell wall development in maize.

Genome-wide expression analysis of phospholipase A1 (PLA1) gene family suggests phospholipase A1-32 gene responding to abiotic stresses in cotton

Cotton is the most important crop for the production of natural fibres used in the textile industry. High salinity, drought, cold and high temperature represent serious abiotic stresses, which seriously threaten cotton production. Phospholipase A_S has an irreplaceable role in lipid signal transmission, growth and development and stress events. Phospholipase A can be divided into three families: PLA1, PLA2 and pPLA. Among them, the PLA1 family is rarely studied in plants. In order to study the potential functions of the PLA1 family in cotton, the bioinformatics analysis of the PLA1 family was correlated with cotton adversity, and tissue-specific analysis was performed. Explore the structure-function relationship of PLA1 members. It is found that the expression of GbPLA1-32 gene is affected by a variety of environmental stimuli, indicating that it plays a very important role in stress and hormone response, and closely associates the cotton adversity with this family. Through further functional verification, we found that virus-induced GbPLA1-32 gene silencing (VIGS) caused Gossypium barbadense to be sensitive to salt stress. This research provides an important basis for further research on the molecular mechanism of cotton resistance to abiotic stress.

Deep transcriptomic study reveals the role of cell wall biosynthesis and organization networks in the developing shell of peanut pod

BACKGROUND

Peanut (Arachis hypogaea L.) belongs to an exceptional group of legume plants, wherein the flowers are produced aerially, but the pods develop under the ground. In such a unique environment, the pod's outer shell plays a vital role as a barrier against mechanical damage and soilborne pathogens. Recent studies have reported the uniqueness and importance of gene expression patterns that accompany peanut pods' biogenesis. These studies focused on biogenesis and pod development during the early stages, but the late developmental stages and disease resistance aspects still have gaps. To extend this information, we analyzed the transcriptome generated from four pod developmental stages of two genotypes, Hanoch (Virginia-type) and IGC53 (Peruvian-type), which differs significantly in their pod shell characteristics and pathogen resistance.

RESULTS

The transcriptome study revealed a significant reprogramming of the number and nature of differentially expressed (DE) genes during shell development. Generally, the numbers of DE genes were higher in IGC53 than in Hanoch, and the R5-R6 transition was the most dynamic in terms of transcriptomic changes. Genes related to cell wall biosynthesis, modification and transcription factors (TFs) dominated these changes therefore, we focused on their differential, temporal and spatial expression patterns. Analysis of the cellulose synthase superfamily identified specific Cellulose synthase (CesAs) and Cellulose synthase-like (Csl) genes and their coordinated interplay with other cell wall-related genes during the peanut shell development was demonstrated. TFs were also identified as being involved in the shell development process, and their pattern of expression differed in the two peanut genotypes. The shell component analysis showed that overall crude fiber, cellulose, lignin, hemicelluloses and dry matter increased with shell development, whereas K, N, protein, and ash content decreased. Genotype IGC53 contained a higher level of crude fiber, cellulose, NDF, ADF, K, ash, and dry matter percentage, while Hanoch had higher protein and nitrogen content.

CONCLUSIONS

The comparative transcriptome analysis identified differentially expressed genes, enriched processes, and molecular processes like cell wall biosynthesis/modifications, carbohydrate metabolic process, signaling, transcription factors, transport, stress, and lignin biosynthesis during the peanut shell development between two contrasting genotypes. TFs and other genes like chitinases were also enriched in peanut shells known for pathogen resistance against soilborne major pathogens causing pod wart disease and pod damages. This study will shed new light on the biological processes involved with underground pod development in an important legume crop.

Pumpkin (Cucurbita moschata) HSP20 Gene Family Identification and Expression Under Heat Stress

Pumpkin (Cucurbita moschata) is an important cucurbit vegetable crop that has strong resistance to abiotic stress. While heat shock protein 20 (HSP20) has been implicated in vegetable response to heat stress, little is known regarding activity of HSP20 family proteins in C. moschata. Here, we performed a comprehensive genome-wide analysis to identify and characterize the functional dynamics of the Cucurbita moschata HSP20 (CmoHSP20) gene family. A total of 33 HSP20 genes distributed across 13 chromosomes were identified from the pumpkin genome. Our phylogenetic analysis determined that the CmoHSP20 proteins fell into nine distinct subfamilies, a division supported by the conserved motif composition and gene structure analyses. Segmental duplication events were shown to play a key role in expansion of the CmoHSP20 gene family. Synteny analysis revealed that 19 and 18 CmoHSP20 genes were collinear with those in the cucumber and melon genomes, respectively. Furthermore, the expression levels of pumpkin HSP20 genes were differentially induced by heat stress. The transcript level of CmoHSP20-16, 24 and 25 were down-regulated by heat stress, while CmoHSP20-7, 13, 18, 22, 26 and 32 were up-regulated by heat stress, which could be used as heat tolerance candidate genes. Overall, these findings contribute to our understanding of vegetable HSP20 family genes and provide valuable information that can be used to breed heat stress resistance in cucurbit vegetable crops.

Genome-wide identification, genomic organization, and expression profiling of the CONSTANS-like (COL) gene family in petunia under multiple stresses

Background

CONSTANS-like (CO-like, COL) are putative zinc-finger transcription factors known to play vital role in various plant biological processes such as control of flowering time, regulation of plant growth and development and responses to stresses. However, no systematic analysis of COL family gene regarding the plant development and stress response has been previously performed in any solanaceous crop. In the present study, a comprehensive genome-wide analysis of COL family genes in petunia has been conducted to figure out their roles in development of organs and stress response.

Results

A total of 33 COL genes, 15 PaCOL genes in P. axillaris and 18 PiCOL genes in P. inflata, were identified in petunia. Subsequently, a genome-wide systematic analysis was performed in 15 PaCOL genes. Considering the domain composition and sequence similarity the 15 PaCOL and 18 PiCOL genes were phylogenetically classified into three groups those are conserved among the flowering plants. Moreover, all of the 15 PaCOL proteins were localized in nucleus. Furthermore, differential expression patterns of PaCOL genes were observed at different developmental stages of petunia. Additionally, transcript expression of 15 PaCOL genes under various abiotic and phytohormone treatments showed their response against stresses. Moreover, several cis-elements related to stress, light-responsive, hormone signaling were also detected in different PaCOL genes.

Conclusion

The phylogenetic clustering, organ specific expression pattern and stress responsive expression profile of conserved petunia COL genes indicating their involvement in plant growth and development and stress response mechanism. This work provide a significant foundation for understanding the biological roles of petunia COL genes in plant growth, development and in stress response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08019-w.

Genome-Wide Identification and Characterization of Polygalacturonase Gene Family in Maize (Zea mays L.)

Polygalacturonase (PG, EC 3.2.1.15) is a crucial enzyme for pectin degradation and is involved in various developmental processes such as fruit ripening, pollen development, cell expansion, and organ abscission. However, information on the PG gene family in the maize (Zea mays L.) genome and the specific members involved in maize anther development are still lacking. In this study, we identified 55 PG family genes from the maize genome and further characterized their evolutionary relationship and expression patterns. Phylogenetic analysis revealed that ZmPGs are grouped into six Clades, and gene structures of the same Clade are highly conserved, suggesting their functional conservation. The ZmPGs are randomly distributed across maize chromosomes, and collinearity analysis showed that many ZmPGs might be derived from tandem duplications and segmental duplications, and these genes are under purifying selection. Furthermore, gene expression analysis provided insights into possible functional divergence among ZmPGs. Based on the RNA-seq data analysis, we found that many ZmPGs are expressed in various tissues while 18 ZmPGs are highly expressed in maize anther, and their detailed expression profiles in different anther developmental stages were further investigated by using RT-qPCR analysis. These results provide valuable information for further functional characterization and application of the ZmPGs in maize.

Genome-Wide Analysis and Expression Profiling of the Phospholipase D Gene Family in Solanum tuberosum

Phospholipase D (PLD) is the most important phospholipid hydrolase in plants, which can hydrolyze phospholipids into phosphatidic acid (PA) and choline. When plants encounter low temperature, drought and high salt stress, phospholipase D and its products play an important role in regulating plant growth and development and coping with stress. In this study, 16 members of StPLD gene family were identified in potato genome, which were distributed in α, β, δ, and ζ subfamilies, and their expression patterns under salt, high temperature, drought, and ABA stress were detected by qRT-PCR method. Gene expression analysis showed that the expression of StPLD genes in potato was upregulated and downregulated to varying degrees under the four stresses, indicating that the PLD gene family is involved in the interaction of potato plant hormones and abiotic stress signals. Chromosome distribution showed that StPLD gene was unevenly distributed on 8 chromosomes, and only one pair of tandem repeat genes was found. All StPLD promoters contain hormone and stress-related cis-regulatory elements to respond to different stresses. Structural analysis showed that StPLD genes in the same subgroup had a similar exon-intron structure. Our study provides a valuable reference for further research of the function and structure of PLD gene.

Changes in secondary metabolites and fiber quality of cotton (Gossypium hirsutum) seed under consecutive water stress and in silico analysis of cellulose synthase and xyloglucan endotransglucosylase

Global warming has led to severe drought conditions. The selection of plant varieties that can withstand drought and produce increased yields are of utmost importance. In the current study, secondary metabolites, seed trait and fiber characteristic of cottonseeds (Gossypium hirsutum) exposed to double and third water stress exposure was investigated. Total phenol and tannin content in W₁S₃₃ increased significantly after third water stress exposure. Accumulation of wax was enhanced in seeds of W₃S₃₃ and W₃S₃₄ that were subjected to third water stress. Fiber quality parameters decreased when cottonseeds were rainfed. High irrigation resulted in fragile and delicate fiber. Seeds grown under 66% FC irrigation saved water and produced seeds that had the potential of producing high quality fibers. In silico analysis was performed on cellulose synthase A (CesA) and xyloglucan endotransglycosylase (XET) enzymes present in Gossypium hirsutum. The intracellular locations of the CesA and XET1 enzymes are the plasma membrane and cell wall, respectively. Proline is conserved in the C-terminal of the CesA enzyme and plays an important role in enzyme functionality. This study provides a better understanding as to the mechanisms by which the plant can tolerate and combat water stress conditions as well as reduce water consumption. In order to grow cotton seeds with desirable morphometric characteristics and optimal fibers under water stress exposure and in dry areas, it is better to use seeds that are irrigated under optimal irrigation conditions, ie 66% FC.

Genome-Wide Association Study in Rice Revealed a Novel Gene in Determining Plant Height and Stem Development, by Encoding a WRKY Transcription Factor

Semi-dwarfism is a main agronomic trait in crop breeding. In this study, we performed genome-wide association study (GWAS) and identified a new quantitative trait nucleotide (QTN) for rice shoot length. The peak QTN (C/T) was located in the first coding region of a group III WRKY transcription factor OsWRKY21 (LOC_Os01g60640). Interestingly, further haplotype analysis showed that C/T difference only existed in the indica group but not in the japonica group, resulting in significant differences in plant height among the different indica rice varieties. OsWRKY21 was expressed in embryo, radicle, shoots, leaves, and stems. Most notably, overexpressing OsWRKY21 resulted in the semi-dwarf phenotype, early heading date and short internodes compared to the wild type, while the knockout mutant plants by CRISPR/Cas9 technology yielded the opposite. The overexpressing lines exhibited the decreased length of the cells near sclerenchyma epidermis, accompanied with the lower levels of indole-3-acetic acid (IAA) and gibberellin 3 (GA₃), but increased levels of the abscisic acid (ABA) and salicylic acid (SA) in the internodes at heading stage. Moreover, the semi-dwarf phenotype could be fully rescued by exogenous GA₃ application at seedling stage. The RNA-seq and qRT-PCR analysis confirmed the differential expression levels of genes in development and the stress responses in rice, including GA metabolism (GA20ox2, GA2ox6, and YABY1) and cell wall biosynthesis (CesA4, 7, and 9) and regulation (MYB103L). These data suggest the essential role of OsWRKY21 in regulation of internode elongation and plant height in rice.

Genome-wide identification of alcohol dehydrogenase (ADH) gene family under waterlogging stress in wheat (Triticum aestivum)

Background

Alcohol dehydrogenase (ADH) plays an important role in plant survival under anaerobic conditions. Although some research about ADH in many plants have been carried out, the bioinformatics analysis of the ADH gene family from Triticum aestivum and their response to abiotic stress is unclear.

Methods

A total of 22 ADH genes were identified from the wheat genome, and these genes could be divided into two subfamilies (subfamily I and subfamily II). All TaADH genes belonged to the Medium-chain ADH subfamily. Sequence alignment analysis showed that all TaADH proteins contained a conservative GroES-like domain and Zinc-binding domain. A total of 64 duplicated gene pairs were found, and the Ka/Ks value of these gene pairs was less than 1, which indicated that these genes were relatively conservative and did not change greatly in the process of duplication.

Results

The organizational analysis showed that nine TaADH genes were highly expressed in all organs, and the rest of TaADH genes had tissue specificity. Cis-acting element analysis showed that almost all of the TaADH genes contained an anaerobic response element. The expression levels of ADH gene in waterlogging tolerant and waterlogging sensitive wheat seeds were analyzed by quantitative real-time PCR (qRT-PCR). This showed that some key ADH genes were significantly responsive to waterlogging stress at the seed germination stage, and the response of waterlogging tolerant and waterlogging sensitive wheat seeds to waterlogging stress was regulated by different ADH genes. The results may be helpful to further study the function of TaADH genes and to determine the candidate gene for wheat stress resistance breeding.

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.