Molecular Characterization, Evolution, and Expression Profiling of the Dirigent (DIR) Family Genes in Chinese White Pear (Pyrus bretschneideri)

Genome-Wide Identification, Evolution, and Expression Analysis of the DIR Gene Family in Schima superba

Schima superba, commonly known as the Chinese guger tree, is highly adaptable and tolerant of poor soil conditions. It is one of the primary species forming the evergreen broad-leaved forests in southern China. Dirigent proteins (DIRs) play crucial roles in the synthesis of plant lignin and lignans, secondary metabolism, and response to adversity stress. However, research on the DIR gene family in S. superba is currently limited. This study identified 24 SsDIR genes, categorizing them into three subfamilies. These genes are unevenly distributed across 13 chromosomes, with 83% being intronless. Collinearity analysis indicated that tandem duplication played a more significant role in the expansion of the gene family compared to segmental duplication. Additionally, we analyzed the expression patterns of SsDIRs in different tissues of S. superba. The SsDIR genes exhibited distinct expression patterns across various tissues, with most being specifically expressed in the roots. Further screening identified SsDIR genes that may regulate drought stress, with many showing differential expression under drought stress conditions. In the promoter regions of SsDIRs, various cis-regulatory elements involved in developmental regulation, hormone response, and stress response were identified, which may be closely related to their diverse regulatory functions. This study will contribute to the further functional identification of SsDIR genes, providing insights into the biosynthetic pathways of lignin and lignans and the mechanisms of plant stress resistance.

Functional and evolutionary comparative analysis of the DIR gene family in Nicotiana tabacum L. and Solanum tuberosum L

BACKGROUND

The dirigent (DIR) genes encode proteins that act as crucial regulators of plant lignin biosynthesis. In Solanaceae species, members of the DIR gene family are intricately related to plant growth and development, playing a key role in responding to various biotic and abiotic stresses. It will be of great application significance to analyze the DIR gene family and expression profile under various pathogen stresses in Solanaceae species.

RESULTS

A total of 57 tobacco NtDIRs and 33 potato StDIRs were identified based on their respective genome sequences. Phylogenetic analysis of DIR genes in tobacco, potato, eggplant and Arabidopsis thaliana revealed three distinct subgroups (DIR-a, DIR-b/d and DIR-e). Gene structure and conserved motif analysis showed that a high degree of conservation in both exon/intron organization and protein motifs among tobacco and potato DIR genes, especially within members of the same subfamily. Total 8 pairs of tandem duplication genes (3 pairs in tobacco, 5 pairs in potato) and 13 pairs of segmental duplication genes (6 pairs in tobacco, 7 pairs in potato) were identified based on the analysis of gene duplication events. Cis-regulatory elements of the DIR promoters participated in hormone response, stress responses, circadian control, endosperm expression, and meristem expression. Transcriptomic data analysis under biotic stress revealed diverse response patterns among DIR gene family members to pathogens, indicating their functional divergence. After 96 h post-inoculation with Ralstonia solanacearum L. (Ras), tobacco seedlings exhibited typical symptoms of tobacco bacterial wilt. The qRT-PCR analysis of 11 selected NtDIR genes displayed differential expression pattern in response to the bacterial pathogen Ras infection. Using line 392278 of potato as material, typical symptoms of potato late blight manifested on the seedling leaves under Phytophthora infestans infection. The qRT-PCR analysis of 5 selected StDIR genes showed up-regulation in response to pathogen infection. Notably, three clustered genes (NtDIR2, NtDIR4, StDIR3) exhibited a robust response to pathogen infection, highlighting their essential roles in disease resistance.

CONCLUSION

The genome-wide identification, evolutionary analysis, and expression profiling of DIR genes in response to various pathogen infection in tobacco and potato have provided valuable insights into the roles of these genes under various stress conditions. Our results could provide a basis for further functional analysis of the DIR gene family under pathogen infection conditions.

Genome-wide identification and characterization of DIRIGENT gene family (CcDIR) in pigeonpea (Cajanus cajan L.) provide insights on their spatial expression pattern and relevance to stress response

This study is a thorough characterization of pigeonpea dirigent gene (CcDIR) family, an important component of the lignin biosynthesis pathway. Genome-wide analysis identified 25 CcDIR genes followed by a range of analytical approaches employed to unravel their structural and functional characteristics. Structural examination revealed a classic single exon and no intron arrangement in CcDIRs contributing to our understanding on evolutionary dynamics. Phylogenetic analysis elucidated evolutionary relationships among CcDIR genes with six DIR sub-families, while motif distribution analysis displayed and highlighted ten conserved protein motifs in CcDIRs. Promoter analyses of all the dirigent genes detected 18 stress responsive cis-acting elements offering insights into transcriptional regulation. While spatial expression analyses across six plant tissues showed preferential expression of CcDIR genes, exposure to salt (CcDIR2 and CcDIR9) and herbivory (CcDIR1, CcDIR2, CcDIR3 and CcDIR11), demonstrated potential roles of specific DIRs in plant defense. Interestingly, increased gene expression during herbivory, also correlated with increased lignin content authenticating the specific response. Furthermore, exogenous application of stress hormones, SA and MeJA on leaves significantly induced the expression of CcDIRs that responded to herbivory. Taken together, these findings contribute to a comprehensive understanding of CcDIR genes impacting development and stress response in the important legume pigeonpea.

Comparative genomics of N-acetyl-5-methoxytryptamine members in four Prunus species with insights into bud dormancy and abiotic stress responses in Prunus avium

KEY MESSAGE

This study provides novel insights into the evolution, diversification, and functions of melatonin biosynthesis genes in Prunus species, highlighting their potential role in regulating bud dormancy and abiotic stresses. The biosynthesis of melatonin (MEL) in plants is primarily governed by enzymatic reactions involving key enzymes such as serotonin N-acetyltransferase (SNAT), tryptamine 5-hydroxylase (T5H), N-acetylserotonin methyltransferase (ASMT) and tryptophan decarboxylase (TDC). In this study, we analyzed Melatonin genes in four Prunus species such as Prunus avium (Pavi), Prunus pusilliflora (Ppus), Prunus serulata (Pser), and Prunus persica (Pper) based on comparative genomics approach. Among the four Prunus species, a total of 29 TDCs, 998 T5Hs, 16 SNATs, and 115 ASMTs within the genome of four Prunus genomes. A thorough investigation of melatonin-related genes was carried out using systematic biological methods and comparative genomics. Through phylogenetic analysis, orthologous clusters, Go enrichment, syntenic relationship, and gene duplication analysis, we discovered both similarities and variations in Melatonin genes among these Prunus species. Additionally, our study revealed the existence of unique subgroup members in the Melatonin genes of these species, which were distinct from those found in Arabidopsis genes. Furthermore, the transcriptomic expression analysis revealed the potential significance of melatonin genes in bud dormancy regulation and abiotic stresses. Our extensive results offer valuable perspectives on the evolutionary patterns, intricate expansion, and functions of PavMEL genes. Given their promising attributes, PavTDCs, PavT5H, PavNAT, and three PavASMT genes warrant in-depth exploration as prime candidates for manipulating dormancy in sweet cherry. This was done to lay the foundation for future explorations into the structural and functional aspects of these factors in Prunus species. This study offers significant insights into the functions of ASMT, SNAT, T5H, and TDC genes and sheds light on their roles in Prunus avium. Moreover, it established a robust foundation for further exploration functional characterization of melatonin genes in fruit species.

Functional analysis of a dirigent protein AtsDIR23 in Acorustatarinowii

Acorus tatarinowii (A. tatarinowii) is a medicinal plant of the Araceae family. Currently, pharmacology focuses on the study of volatile oils, but there are few reports of another important secondary metabolite, lignan. Dirigent protein is thought to play an important role in plant secondary metabolism and responds to a variety of biotic and abiotic stresses. However, the DIR gene family of A. tatarinowii has not been systematically analyzed, and it is unknown whether it affects lignan synthesis. In this study, a total of 27 AtsDIRs were identified by comprehensive analysis of the genome of the medicinal plant A. tatarinowii, and the candidate gene AtsDIR23 that may be involved in lignan synthesis was screened through bioinformatics and transcriptome analysis. It is worth noting that AtsDIR23 is significantly expressed in rhizomes and is a member of the DIR-a subfamily. Subsequently, subcellular localization revealed that AtsDIR23 was localized in chloroplasts. The functional verification of AtsDIR23 b y the transient transformation of A. tatarinowii and the stable transformation of Arabidopsis thaliana showed that the content of lignans in overexpressed plants increased. Co-expression analysis screening revealed the MYB transcription factor (AtsMYB91) that is highly correlated with AtsDIR23 expression, while yeast one-hybrid assays and double luciferase experiments showed that AtsMYB91 negatively regulated the expression of AtsDIR23 b y binding to the AtsDIR23 promoter. In conclusion, AtsDIR23 can promote the accumulation of lignans, which provides a reference for further research on the regulation of lignans by DIR genes.

Genome-wide identification, characterization, evolution and expression analysis of the DIR gene family in potato (Solanum tuberosum)

The dirigent (DIR) gene is a key player in environmental stress response and has been identified in many multidimensional tube plant species. However, there are few studies on the StDIR gene in potato. In this study, we used genome-wide identification to identify 31 StDIR genes in potato. Among the 12 potato chromosomes, the StDIR gene was distributed on 11 chromosomes, among which the third chromosome did not have a family member, while the tenth chromosome had the most members with 11 members. 22 of the 31 StDIRs had a classical DIR gene structure, with one exon and no intron. The conserved DIR domain accounts for most of the proteins in the 27 StDIRs. The structure of the StDIR gene was analyzed and ten different motifs were detected. The StDIR gene was divided into three groups according to its phylogenetic relationship, and 22 duplicate genes were identified. In addition, four kinds of cis-acting elements were detected in all 31 StDIR promoter regions, most of which were associated with biotic and abiotic stress. The findings demonstrated that the StDIR gene exhibited specific responses to cold stress, salt stress, ABA, and drought stress. This study provides new candidate genes for improving potato's resistance to stress.

Identification and functional characterization of the dirigent gene family in Phryma leptostachya and the contribution of PlDIR1 in lignan biosynthesis

BACKGROUND

Furofuran lignans, the main insecticidal ingredient in Phryma leptostachya, exhibit excellent controlling efficacy against a variety of pests. During the biosynthesis of furofuran lignans, Dirigent proteins (DIRs) are thought to be dominant in the stereoselective coupling of coniferyl alcohol to form ( ±)-pinoresinol. There are DIR family members in almost every vascular plant, but members of DIRs in P. leptostachya are unknown. To identify the PlDIR genes and elucidate their functions in lignan biosynthesis, this study performed transcriptome-wide analysis and characterized the catalytic activity of the PlDIR1 protein.

RESULTS

Fifteen full-length unique PlDIR genes were identified in P. leptostachya. A phylogenetic analysis of the PlDIRs classified them into four subfamilies (DIR-a, DIR-b/d, DIR-e, and DIR-g), and 12 conserved motifs were found among them. In tissue-specific expression analysis, except for PlDIR7, which displayed the highest transcript abundance in seeds, the other PlDIRs showed preferential expression in roots, leaves, and stems. Furthermore, the treatments with signaling molecules demonstrated that PlDIRs could be significantly induced by methyl jasmonate (MeJA), salicylic acid (SA), and ethylene (ETH), both in the roots and leaves of P. leptostachya. In examining the tertiary structure of the protein and the critical amino acids, it was found that PlDIR1, one of the DIR-a subfamily members, might be involved in the region- and stereo-selectivity of the phenoxy radical. Accordingly, LC-MS/MS analysis demonstrated the catalytic activity of recombinant PlDIR1 protein from Escherichia coli to direct coniferyl alcohol coupling into ( +)-pinoresinol. The active sites and hydrogen bonds of the interaction between PlDIR1 and bis-quinone methide (bisQM), the intermediate in ( +)-pinoresinol formation, were analyzed by molecular docking. As a result, 18 active sites and 4 hydrogen bonds (Asp-42, Ala-113, Leu-138, Arg-143) were discovered in the PlDIR1-bisQM complex. Moreover, correlation analysis indicated that the expression profile of PlDIR1 was closely connected with lignan accumulations after SA treatment.

CONCLUSIONS

The results of this study will provide useful clues for uncovering P. leptostachya's lignan biosynthesis pathway as well as facilitate further studies on the DIR family.

Insights into the functional characterization of DIR proteins through genome-wide in silico and evolutionary studies: a systematic review

Dirigent proteins (DIRs) are a new class of proteins that were identified during the 8-8' lignan biosynthetic pathway and involves the formation of ( +) or ( -)-pinoresinol through stereoselective coupling from E-coniferyl alcohol. These proteins are known to play a vital role in the development and stress response in plants. Various studies have reported the functional and structural characterization of dirigent gene family in different plants using in silico approaches. Here, we have summarized the importance of dirigent proteins in plants and their role in plant stress tolerance by analyzing the genome-wide analysis including gene structure, mapping of chromosomes, phylogenetic evolution, conserved motifs, gene structure, and gene duplications in important plants. Overall, this review would help to compare and clarify the molecular and evolutionary characteristics of dirigent gene family in different plants.

Characterization, expression, and functional analysis of the pathogenesis-related gene PtDIR11 in transgenic poplar

Lignins and lignans are important for plant resistance to pathogens. Dirigent (DIR) proteins control the regio- and stereo-selectivity of coniferyl alcohol in lignan and lignin biosynthesis. DIR genes have been implicated in defense-related responses in several plant species, but their role in poplar immunity is unclear. We cloned PtDIR11 from Populus trichocarpa; we found that overexpression of PtDIR11 in poplar improved the lignan biosynthesis and enhanced the resistance of poplar to Septotis populiperda. PtDIR11 has a typical DIR domain; it belongs to the DIR-b/d family and is expressed in the cell membrane. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis showed that PtDIR11 expression was highest in stems, followed by leaves and roots. Furthermore, PtDIR11 expression was induced by S. populiperda, salicylic acid (SA), jasmonate (JA), and ethylene (ET) stresses. The recombinant PtDIR11 protein inhibited the growth of S. populiperda in vitro. Overexpressing (OE) PtDIR11 in "Nanlin 895" poplar enhanced growth. The OE lines exhibited minimal changes in lignin content, but their total lignan and flavonoid contents were significantly greater than in the wild-type (WT) lines. Overexpression of PtDIR11 affected multiple biological pathways of poplar, such as phenylpropanoid biosynthesis. The methanol extracts of OE-PtDIR11 lines showed greater anti-S. populiperda activity than did lignin extracts from the WT lines. Furthermore, OE-PtDIR11 lines upregulated genes that were related to phenylpropanoid biosynthesis and genes associated with the JA and ET signal transduction pathways; it downregulated genes that were related to SA signal transduction compared with the WT line under S. populiperda stress. Therefore, the OE transgenic plants analysis revealed that PtDIR11 is a good candidate gene for breeding of disease resistant poplar.

Overexpression of Sugarcane ScDIR Genes Enhances Drought Tolerance in Nicotiana benthamiana

Dirigent proteins (DIRs) are known to function in lignin biogenesis and to be involved in stress resistance in plants. However, the sugarcane DIRs have not been functionally characterized. In this study, we investigated the DIR-protein-encoding genes in Saccharum spp. (ScDIR) by screening collections of sugarcane databases, monitoring the responses of these genes to drought stress by real-time quantitative PCR, and identifying their heterologous expression in tobacco. Of the 64 ScDIRs identified, four belonging to the DIR-b/d (ScDIR5 and ScDIR11) and DIR-c (ScDIR7 and ScDIR40) subfamilies showed a significant transcriptional response when subjected to drought stress. ScDIR5, ScDIR7, and ScDIR11 are localized in the cell membrane, whereas ScDIR40 is found in the cell wall. The overexpression of these ScDIR genes in tobacco generally increased the drought tolerance of the transgenic lines, with ScDIR7 conferring the highest degree of drought tolerance. The characterization of the physiological and biochemical indicators (superoxide dismutase, catalase, malondialdehyde, and H2O2) confirmed that the ScDIR-overexpressing lines outperformed the wild type. These results demonstrated that specific ScDIRs in sugarcane respond and contribute to tolerance of drought stress, shedding light on potential means of improving drought tolerance in this crop.

Whole-genome analysis of CGS, SAHH, SAMS gene families in five Rosaceae species and their expression analysis in Pyrus bretschneideri

Cystathionine γ-synthase (CGS), S-adenosyl-L-homocysteine hydrolase (SAHH), and S-adenosy-L-methionine synthetase (SAMS) play an important role in the regulation of plant growth, development, and secondary metabolism. In this study, a total of 6 CGS, 6 SAHH, and 28 SAMS genes were identified from five Rosaceae species (Pyrus bretschneideri, Prunus persica, Prunus mume, Fragaria vesca, and Malus domestica). The evolutionary relationship and microsynteny analysis in five Rosaceae species revealed that duplicated regions were conserved between three gene families (CGS, SAHH, SAMS). Moreover, the chromosomal locations, gene structures, conserved motifs, cis-elements, physicochemical properties, and Ka/Ks analysis were performed by using numerous bioinformatics tools. The expression of different organs showed that the CGS, SAHH and SAMS genes of pear have relatively high expression patterns in flowers and stems, except for PbCGS1. RNA-seq and qRT-PCR combined analysis showed that PbSAMS1 may be involved in the regulation of pear stone cell development. In summary, this study provides the basic information of CGS, SAHH and SAMS genes in five Rosaceae species, further revealing the expression patterns in the pear fruit, which provides the theoretical basis for the regulation of pear stone cells.

Genome-wide investigation and comparative analysis of MATE gene family in Rosaceae species and their regulatory role in abiotic stress responses in Chinese pear (Pyrus bretschneideri)

The Multidrug and Toxic Compound Extrusion (MATE) protein belongs to a secondary transporter gene family, which plays a primary role in transporting many kinds of substrates such as organic compounds, secondary metabolites, and phytohormones. MATE protein members exist in both prokaryotes and eukaryotes. However, evolution and comprehensive analysis of the MATE genes has not been performed in Rosaceae species. In the present study, a total of 404 MATEs genes were identified from six Rosaceae genomes (Prunus avium, Pyrus bretschneideri, Prunus persica, Fragaria vesca, Prunus mume, and Malus domestica) and classified into eight main subfamilies (I-VII) based on structural and phylogenetic analysis. Microcollinearity analysis showed that whole-genome duplication events might play a vital role in the expansion of the MATE genes family. The Ka/Ks analysis, chromosomal localization, subcellular localization, and molecular characteristics (length, weight, and pI) were performed using various bioinformatics tools. Furthermore, different subfamilies have different introns-exons structures, cis-acting elements, and conserved motifs analysis, indicating functional divergence in the MATE family. Subsequently, RNA-seq analysis and real-time qRT-PCR were conducted during Chinese pear fruit development. Moreover, PbMATE genes were significantly expressed under hormonal treatments of MeJA (methyl jasmonate), SA (salicylic acid), and ABA (abscisic acid). Overall, our results provide helpful insights into the functions, expansion complexity, and evolutions of the MATE genes in Chinese pear and five Rosaceae species.

The whole-genome and expression profile analysis of WRKY and RGAs in Dactylis glomerata showed that DG6C02319.1 and DgWRKYs may cooperate in the immunity against rust

Orchardgrass (Dactylis glomerata) is one of the top four perennial forages worldwide and, despite its large economic advantages, often threatened by various environmental stresses. WRKY transcription factors (TFs) can regulate a variety of plant processes, widely participate in plant responses to biotic and abiotic stresses, and are one of the largest gene families in plants. WRKYs can usually bind W-box elements specifically. In this study, we identified a total of 93 DgWRKY genes and 281 RGAs, including 65, 169 and 47 nucleotide-binding site-leucine-rich repeats (NBS-LRRs), leucine-rich repeats receptor-like protein kinases (LRR-RLKs), and leucine-rich repeats receptor-like proteins (LRR-RLPs), respectively. Through analyzing the expression of DgWRKY genes in orchardgrass under different environmental stresses, it was found that many DgWRKY genes were differentially expressed under heat, drought, submergence, and rust stress. In particular, it was found that the greatest number of genes were differentially expressed under rust infection. Consistently, GO and KEGG enrichment analysis of all genes showed that 78 DgWRKY TFs were identified in the plant–pathogen interaction pathway, with 59 of them differentially expressed. Through cis-acting element prediction, 154 RGAs were found to contain W-box elements. Among them, DG6C02319.1 (a member of the LRR-RLK family) was identified as likely to interact with 14 DGWRKYs. Moreover, their expression levels in susceptible plants after rust inoculation were first up-regulated and then down-regulated, while those in the resistant plants were always up-regulated. In general, DgWRKYs responded to both biotic stress and abiotic stress. DgWRKYs and RGAs may synergistically respond to the response of orchardgrass to rust. This study provides meaningful insight into the molecular mechanisms of WRKY proteins in orchardgrass.

Genome-Wide Characterization of Dirigent Proteins in Populus: Gene Expression Variation and Expression Pattern in Response to Marssonina brunnea and Phytohormones

Marssonina brunnea causes a major disease that limits poplar growth. Lignin and lignan play essential roles in protecting plants from various biological stresses. Dirigent (DIR) proteins are thought to control the stereoselective coupling of coniferyl alcohol in the formation of lignan and lignin. DIR family members have been well studied in several plant species, but no previous detailed genome-wide analysis has been carried out in forest trees, such as poplar. We identified 40 PtDIR genes in Populus trichocarpa and classified them into three subgroups (DIR-a, DIR-b/d, and DIR-e) based on phylogenetic analyses. These genes are distributed on 11 poplar chromosomes, and 80% of PtDIRs (32/40) are intronless. The cis-element analysis inferred that PtDIRs possess many types of biological and abiotic stress-response cis-elements. We also analyzed intra- and inter-specific collinearity, which provided deep insights into the evolutionary characteristics of the poplar DIR genes. Analyses of the protein tertiary structure and critical amino acid residues showed that PtDIR7–10 and PtDIR13–16, which belong to the DIR-a subfamily, might be involved in the regio- and stereo-selectivity of bimolecular phenoxy radical coupling in poplars. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) analysis revealed different expression patterns for the PtDIR genes of P. trichocarpa and the PeDIR genes of ‘Nanlin 895’ in various tissues. Additionally, we analyzed responses of PeDIRs to M. brunnea and different phytohormone treatments (abscisic acid, salicylic acid, methyl jasmonate, and ethylene) in ‘Nanlin 895’. The results showed that at least 18 genes responded strongly to M. brunnea, and these PeDIRs also showed significant responses to phytohormones. These results suggest that DIR genes are involved in the poplar defense response against M. brunnea, and this study will provide fundamental insights for future research on poplar DIR genes.

Genome-wide identification and expression analysis of dirigent-jacalin genes from plant chimeric lectins in Moso bamboo (Phyllostachys edulis)

Dirigent-jacalin (D-J) genes belong to the plant chimeric lectin family, and play vital roles in plant growth and resistance to abiotic and biotic stresses. To explore the functions of the D-J family in the growth and development of Moso bamboo (Phyllostachys edulis), their physicochemical properties, phylogenetic relationships, gene and protein structures, and expression patterns were analyzed in detail. Four putative PeD-J genes were identified in the Moso bamboo genome, and microsynteny and phylogenetic analyses indicated that they represent a new branch in the evolution of plant lectins. PeD-J proteins were found to be composed of a dirigent domain and a jacalin-related lectin domain, each of which contained two different motifs. Multiple sequence alignment and homologous modeling analysis indicated that the three-dimensional structure of the PeD-J proteins was significantly different compared to other plant lectins, primarily due to the tandem dirigent and jacalin domains. We surveyed the upstream putative promoter regions of the PeD-Js and found that they mainly contained cis-acting elements related to hormone and abiotic stress response. An analysis of the expression patterns of root, leaf, rhizome and panicle revealed that four PeD-J genes were highly expressed in the panicle, indicating that they may be required during the formation and development of several different tissue types in Moso bamboo. Moreover, PeD-J genes were shown to be involved in the rapid growth and development of bamboo shoots. Quantitative Real-time PCR (qRT PCR) assays further verified that D-J family genes were responsive to hormones and stresses. The results of this study will help to elucidate the biological functions of PeD-Js during bamboo growth, development and stress response.

Comparative Analysis, Characterization and Evolutionary Study of Dirigent Gene Family in Cucurbitaceae and Expression of Novel Dirigent Peptide against Powdery Mildew Stress

Dirigent (DIR) proteins are induced under various stress conditions and involved in sterio- and regio-selective coupling of monolignol. A striking lack of information about dirigent genes in cucurbitaceae plants underscores the importance of functional characterization. In this study, 112 DIR genes were identified in six species, and 61 genes from major cultivated species were analyzed. DIRs were analyzed using various bioinformatics tools and complemented by expression profiling. Phylogenetic analysis segregated the putative DIRs into six distinctively known subgroups. Chromosomal mapping revealed uneven distribution of genes, whereas synteny analysis exhibited that duplication events occurred during gene evolution. Gene structure analysis suggested the gain of introns during gene diversification. Gene ontology (GO) enrichment analysis indicates the participation of proteins in lignification and pathogen resistance activities. We also determined their organ-specific expression levels in three species revealing preferential expression in root and leaves. Furthermore, the number of CmDIR (CmDIR1, 6, 7 and 12) and ClDIR (ClDIR2, 5, 8, 9 and 17) genes exhibited higher expression in resistant cultivars after powdery mildew (PM) inoculation. In summary, based on the expression and in-silico analysis, we propose a role of DIRs in disease resistance mechanisms.

Evolution, expression and functional analysis of cultivated allotetraploid cotton DIR genes

BACKGROUND

Dirigent (DIR) proteins mediate regioselectivity and stereoselectivity during lignan biosynthesis and are also involved in lignin, gossypol and pterocarpan biosynthesis. This gene family plays a vital role in enhancing stress resistance and in secondary cell-wall development, but systematical understanding is lacking in cotton.

RESULTS

In this study, 107 GbDIRs and 107 GhDIRs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. Most of these genes have a classical gene structure without intron and encode proteins containing a signal peptide. Phylogenetic analysis showed that cotton DIR genes were classified into four distinct subfamilies (a, b/d, e, and f). Of these groups, DIR-a and DIR-e were evolutionarily conserved, and segmental and tandem duplications contributed equally to their formation. In contrast, DIR-b/d mainly expanded by recent tandem duplications, accompanying with a number of gene clusters. With the rapid evolution, DIR-b/d-III was a Gossypium-specific clade involved in atropselective synthesis of gossypol. RNA-seq data highlighted GhDIRs in response to Verticillium dahliae infection and suggested that DIR gene family could confer Verticillium wilt resistance. We also identified candidate DIR genes related to fiber development in G. barbadense and G. hirsutum and revealed their differential expression. To further determine the involvement of DIR genes in fiber development, we overexpressed a fiber length-related gene GbDIR78 in Arabidopsis and validated its function in trichomes and hypocotyls.

CONCLUSIONS

These findings contribute novel insights towards the evolution of DIR gene family and provide valuable information for further understanding the roles of DIR genes in cotton fiber development as well as in stress responses.

In silico identification and validation of miRNA and their DIR specific targets in Oryza sativa Indica under abiotic stress

Several biotic (bacterial and viral pathogenesis) and abiotic stress factors like salt, drought, cold, and extreme temperatures significantly reduce crop productivity and grain quality throughout the world. MicroRNAs (miRNAs) are small (~22 nucleotides) non-coding endogenous RNA molecules which negatively regulate gene expression at the post-transcriptional level either by degrading the target protein-coding mRNA genes or suppressing translation in plants. Dirigent (DIR) gene protein plays a crucial role as they are involved to dictate the stereochemistry of a compound synthesized by other enzymes as well as in lignifications against biotic and abiotic stress. In plants, several miRNAs, as well as their targets, are known to regulate stress response but systematic identification of the same is limited. The present work has been designed for in silico identification of miRNAs against a total of sixty-one DIR genes in Oryza sativa Indica followed by target prediction of identified miRNAs through the computational approach and thereafter validation of potential miRNAs in rice genotypes. We systematically identified 3 miRNA and their respective DIR specific target gene in Oryza sativa Indica. The expression of these three miRNAs and their respective DIR specific targets were validated in rice seedlings subjected to five different abiotic stress conditions (heavy metal, high temperature, low temperature, salinity and drought) by quantitative Real-Time PCR (qRT-PCR). Expression analysis indicated that miRNA under stress conditions regulates the gene expression of the DIR gene in rice. To the best of our knowledge this is this is the first report in any organism showing the expression of ath-miRf10317-akr, and osamiRf10761-akr miRNAs in response to various abiotic stresses.

•

Total 61 DIR proteins were identified & classified into 6 groups based on phylogeny analysis in Oryza sativa Indica.

•

Three miRNAs ath-miRf10317-akr, cre-miR910 and osa-miRf10761-akr were identified via computational approach.

•

These 3 miRNAs in response to abiotic stresses showed inverse expression pattern in the respective target genes.

•

This is the first report on expression of ath-miRf10317-akr, and osa-miRf10761-akr miRNAs in response to abiotic stresses.

Gene structure, evolution and expression analysis of the P-ATPase gene family in Chinese pear (Pyrus bretschneideri)

P-ATPase are a large protein family of integral membrane, playing an important role in plant growth, development and stress. P-ATPase genes family have been identified and characterized in several model plants such as cotton, grapes, tobacco, rice, rubber plant and Arabidopsis. However, still lack of comprehensive study of P-ATPase genes in Chinese pear (Pyrus bretschneideri). A systematic analysis was performed and identified 30 P-ATPase genes from the pear genome to evaluate the qualities and diversity of P-ATPase proteins. Phylogenetic analysis was performed using A. thaliana P-ATPase genes as a model, allowing us to categorize into 4 subfamilies (PbHMA, PbECA, PbACA, and PbAHA) and two subfamilies (ALA and P₅) is absent in pear. Even Within the same subclade, P-ATPase genes also shows the similar exon-intron structure and conserved motif structure. Continuing chromosomal localization analysis showed that 23 P-ATPase genes were distributed among 13 chromosome and 7 gene on the scaffold of pear. Promoter regions of P-ATPase genes revealed that several cis-acting elements were involved in plant growth/development, stress responses as well as hormone responses. Additionally, P-ATPase genes were also differentially expressed under hormones treatments of ABA (abscisic acid) and SA (salicylic acid) treatments. Remarkably, the transcriptome data exposed that P-ATPase gene might play an important role in lignin biosynthesis during fruit development. The real time qRT-PCR was performed, and the expression analysis indicated that various P-ATPase genes extremely expressed during different developmental stages of fruit. Our study provides valuable information about the P-ATPase gene family in pear fruit development and lignin polymerization.

Effects of Metaxenia on Stone Cell Formation in Pear (Pyrus bretschneideri) Based on Transcriptomic Analysis and Functional Characterization of the Lignin-Related Gene PbC4H2

The deposition of lignin in flesh parenchyma cells for pear stone cells, and excessive stone cells reduce the taste and quality of the fruit. The effect of metaxenia on the quality of fruit has been heavily studied, but the effect of metaxenia on stone cell formation has not been fully elucidated to date. This study used P. bretschneideri (Chinese white pear) cv. ‘Yali’ (high-stone cell content) and P. pyrifolia (Sand pear) cv. ‘Cuiguan’ (low-stone cell content) as pollination trees to pollinate P. bretschneideri cv. ‘Lianglizaosu’ separately to fill this gap in the literature. The results of quantitative determination, histochemical staining and electron microscopy indicated that the content of stone cells and lignin in YL fruit (‘Yali’ (pollen parent) × ‘Lianglizaosu’ (seed parent)) was significantly higher than that in CL fruit (‘Cuiguan’ (pollen parent) × ‘Lianglizaosu’ (seed parent)). The transcriptome sequencing results that were obtained from the three developmental stages of the two types of hybrid fruits indicated that a large number of differentially expressed genes (DEGs) related to auxin signal transduction (AUX/IAAs and ARFs), lignin biosynthesis, and lignin metabolism regulation (MYBs, LIMs, and KNOXs) between the CL and YL fruits at the early stage of fruit development. Therefore, metaxenia might change the signal transduction process of auxin in pear fruit, thereby regulating the expression of transcription factors (TFs) related to lignin metabolism, and ultimately affecting lignin deposition and stone cell development. In addition, we performed functional verification of a differentially expressed gene, PbC4H2 (cinnamate 4-hydroxylase). Heterologous expression of PbC4H2 in the c4h mutant not only restored its collapsed cell wall, but also significantly increased the lignin content in the inflorescence stem. The results of our research help to elucidate the metaxenia-mediated regulation of pear stone cell development and clarify the function of PbC4H2 in cell wall development and lignin synthesis, which establishes a foundation for subsequent molecular breeding.

Biochemical characterization and expression analysis of lignification in two pear (Pyrus ussuriensis Maxim.) varieties with contrasting stone cell content

As lignified stone cells reduce fruit quality, we investigated lignin deposition, phenolic metabolites, and expression of lignin biosynthetic genes during fruit development to elucidate the molecular mechanism of stone cell lignification using histological, biochemical, and transcriptional data from two Ussurian pear varieties (Jianba and Nanguo) with contrasting stone cell content. Lignin content and distribution coincided with stone cell accumulation. As per LC-MS analysis, Jianba exhibited higher levels of lignin monomers and hydroxycinnamates than Nanguo, consistently with lignin amount in each case. However, flavonoid content was much higher in Nanguo. Transcriptional data showed that most monolignol biosynthesis-related genes were particularly upregulated in Jianba during lignin accumulation; especially CCR and LAC, two monolignol biosynthesis-specific genes, were substantially upregulated in Jianba fruits at critical stages. Therefore, differences in stone cell content between "Jianba" and "Nanguo" may result from differential expression of lignin synthase genes located downstream of the lignin biosynthesis pathway. Taken together, our data may provide a deeper understanding of the molecular mechanism for stone cell lignification in pear fruit.

Identification of Homeobox Genes Associated with Lignification and Their Expression Patterns in Bamboo Shoots

: Homeobox (HB) genes play critical roles in regulating various aspects of plant growth and development. However, little is known about HB genes in bamboo. In this study, a total of 115 HB genes (PeHB001‒PeHB115) were identified from moso bamboo (Phyllostachys edulis) and grouped into 13 distinct classes (BEL, DDT, HD-ZIP I‒IV, KNOX, NDX, PHD, PINTOX, PLINC, SAWADEE, and WOX) based on the conserved domains and phylogenetic analysis. The number of members in the different classes ranged from 2 to 24, and they usually varied in terms of exon‒intron distribution pattern and length. There were 20 conserved motifs found in 115 PeHBs, with motif 1 being the most common. Gene ontology (GO) analysis showed that PeHBs had diverse molecular functions, with 19 PeHBs being annotated as having xylem development, xylem, and phloem pattern formation functions. Co-expression network analysis showed that 10 of the 19 PeHBs had co-expression correlations, and three members of the KNOX class were hub proteins that interacted with other transcription factors (TFs) such as MYB, bHLH, and OVATE, which were associated with lignin synthesis. Yeast two-hybridization results further proved that PeHB037 (BEL class) interacted with PeHB057 (KNOX class). Transcriptome expression profiling indicated that all PeHBs except PeHB017 were expressed in at least one of the seven tissues of moso bamboo, and 90 PeHBs were expressed in all the tissues. The qRT-PCR results of the 19 PeHBs showed that most of them were upregulated in shoots as the height increased. Moreover, a KNOX binding site was found in the promoters of the key genes involved in lignin synthesis such as Pe4CL, PeC3H, PeCCR, and PeCOMT, which had positive expression correlations with five KNOX genes. Similar results were found in winter bamboo shoots with prolonged storage time, which was consistent with the degree of lignification. These results provide basic data on PeHBs in moso bamboo, which will be helpful for future functional research on PeHBs with positive regulatory roles in the process of lignification.

The Resurgence of Dirigent Story: Time for a Bacterial Chapter

For several decades, dirigent (DIR) domain-containing proteins have been assumed to be green lineage-specific, responsible for the defence response and lignan/lignin biosynthesis. Despite their high potential in terms of biotechnology and chemistry, to date there have been very few well-studied plant DIRs. However, recent achievements in sequencing technologies have allowed for discovery of DIR genes in bacteria. This prospective study suggests expansion of the focus of research to consider the existence of bacterial DIRs. It also considers the outlook for understanding DIR functioning with respect to the fields of green lineage evolution, organic synthesis, and biotechnology.

In Silico Genome-Wide Analysis of the Pear (Pyrus bretschneideri) KNOX Family and the Functional Characterization of PbKNOX1, an Arabidopsis BREVIPEDICELLUS Orthologue Gene, Involved in Cell Wall and Lignin Biosynthesis

Stone cells are a characteristic trait of pear fruit, but the contents and sizes of stone cells negatively correlate with fruit texture and flavor. Secondary cell wall thickening and lignification have been established as key steps of stone cell development. KNOTTED-LIKE HOMEOBOX (KNOX) proteins play important roles in plant cell growth and development, including cell wall formation and lignification. Although the characteristics and biological functions of KNOX proteins have been investigated in other plants, this gene family has not been functionally characterized in pear. Eighteen PbKNOX genes were identified in the present study, and all of the identified family members contained the KNOX I and/or KNOX II domains. Based on the phylogenetic tree and chromosomal localization, the 18 PbKNOX genes were divided into five subfamilies [SHOOT MERISTEMLESS (STM)-like, BREVIPEDICELLUS (BP)-like, KNOTTED ARABIDOPSIS THALIANA 2/6 (KNAT2/6)-like, KNAT7-like, and KNAT3-5-like] and were distributed among 10 chromosomes. In addition, we identified 9, 11, and 11 KNOX genes in the genomes of grape, mei, and strawberry, respectively, and the greatest number of collinear KNOX gene pairs formed between pears and peaches. Analyses of the spatiotemporal expression patterns showed that the tissue specificity of PbKNOX gene expression was not very significant and that the level of the PbKNOX1 transcript showed an opposite trend to the levels of stone cells and lignin accumulation. Furthermore, PbKNOX1 has high sequence identity and similarity with Arabidopsis BP. Compared with wild-type Arabidopsis, plants overexpressing PbKNOX1 not only showed an approximately 19% decrease in the secondary cell wall thickness of vessel cells but also exhibited an approximately 13% reduction in the lignin content of inflorescence stems. Moreover, the expression of several genes involved in lignin biosynthesis was downregulated in transgenic lines. Based on our results, PbKNOX1/BP participates in cell wall-thickening and lignin biosynthesis and represses the transcription of key structural genes involved in lignin synthesis, providing genetic evidence for the roles of KNOX in cell wall thickening and lignin biosynthesis in pear.

In Silico Genome-Wide Analysis of Respiratory Burst Oxidase Homolog (RBOH) Family Genes in Five Fruit-Producing Trees, and Potential Functional Analysis on Lignification of Stone Cells in Chinese White Pear

: The accumulation of lignin in fruit has a significant negative impact on the quality of fruit-producing trees, and in particular the lignin formation stimulates the development of stone cells in pear fruit. Reactive oxygen species (ROS) are essential for lignin polymerization. However, knowledge of the RBOH family, a key enzyme in ROS metabolism, remains unknown in most fruit trees. In this study, a total of 40 RBOHs were identified from five fruit-producing trees (Pyrusbretschneideri, Prunuspersica, Citrussinensis, Vitisvinifera, and Prunusmume), and 10 of these sequences came from Pyrusbretschneideri. Multiple sequence alignments revealed that all 10 PbRBOHs contained the NADPH_Ox domain and the six alpha-helical transmembrane domains (TM-I to TM-VI). Chromosome localization and interspecies phylogenetic tree analysis showed that 10 PbRBOHs irregularly distributed on 8 chromosomes and 3 PbRBOHs (PbRBOHA, PbRBOHB, and PbRBOHD) are closely related to known lignification-related RBOHs. Furthermore, hormone response pattern analysis showed that the transcription of PbRBOHs is regulated by SA, ABA and MeJA. Reverse transcription-quantitative real-time polymerase chain reaction (qRT-PCR) and transcriptome sequencing analysis showed that PbRBOHA, PbRBOHB, and PbRBOHD accumulated high transcript abundance in pear fruit, and the transcriptional trends of PbRBOHA and PbRBOHD was consistent with the change of stone cell content during fruit development. In addition, subcellular localization revealed that PbRBOHA and PbRBOHD are distributed on the plasma membrane. Combining the changes of apoplastic superoxide (O2.-) content and spatio-temporal expression analysis, these results indicate that PbRBOHA and PbRBOHD, which are candidate genes, may play an important role in ROS metabolism during the lignification of pear stone cells. This study not only provided insight into the molecular characteristics of the RBOH family in fruit-producing trees, but also lays the foundation for studying the role of ROS in plant lignification.

Comprehensive genome-wide analysis of the pear (Pyrus bretschneideri) laccase gene (PbLAC) family and functional identification of PbLAC1 involved in lignin biosynthesis

The content and size of stone cell clusters affects the quality of pear fruit, and monolignol polymerization and deposition in the cell walls constitute a required step for stone cell formation. Laccase (LAC) is the key enzyme responsible for the polymerization of monolignols. However, there are no reports on the LAC family in pear (Pyrus bretschneideri), and the identity of the members responsible for lignin synthesis has not been clarified. Here, 41 LACs were identified in the whole genome of pear. All Pyrus bretschneideri LACs (PbLACs) were distributed on 13 chromosomes and divided into four phylogenetic groups (I-IV). In addition, 16 segmental duplication events were found, implying that segmental duplication was a primary reason for the expansion of the PbLAC family. LACs from the genomes of three Rosaceae species (Prunus mummer, Prunus persica, and Fragaria vesca) were also identified, and an interspecies collinearity analysis was performed. The phylogenetic analysis, sequence alignments and spatiotemporal expression pattern analysis suggested that PbLAC1, 5, 6, 29, 36 and 38 were likely associated with lignin synthesis and stone cell formation in fruit. The two target genes of Pyr-miR1890 (a microRNA identified from pear fruit that is associated with lignin and stone cell accumulation), PbLAC1 and PbLAC14, were selected for genetic transformation. Interfamily transfer of PbLAC1 into Arabidopsis resulted in a significant increase (approximately 17%) in the lignin content and thicker cell walls in interfascicular fibre and xylem cells, which demonstrated that PbLAC1 is involved in lignin biosynthesis and cell wall development. However, the lignin content and cell wall thickness were not changed significantly in the PbLAC14-overexpressing transgenic Arabidopsis plants. This study revealed the function of PbLAC1 in lignin synthesis and provides important insights into the characteristics and evolution of the PbLAC family.

Family-1 UDP glycosyltransferases in pear (Pyrus bretschneideri): Molecular identification, phylogenomic characterization and expression profiling during stone cell formation

Stone cells are a characteristic trait of pear fruits, and excessive stone cell formation has a significant negative impact on the texture and flavour of the pulp. Lignin is one of the main components of stone cells. Family-1 uridine diphosphate-glycosyltransferases (UGTs) are responsible for the glycosylation modification of monolignols. However, information remains limited regarding the relationship between UGTs and stone cell formation. To address this problem, we identified 139 UGTs from the pear genome, which were distributed in 15 phylogenetic groups (A-M, O, and P). We also performed a collinearity analysis of UGTs among four Rosaceae plants (pear, peach, mei, and strawberry). Phylogenetic analysis suggested that 13 PbUGTs might be related to the glycosylation of monolignols. Analysis of expression patterns demonstrated that most putative monolignol glycosylation-related PbUGTs not only showed high expression levels in flowers and buds but were also induced by exogenous ABA, SA, and MeJA. In addition, the transcript level of Pbr005014.1 (named PbUGT72AJ2) was consistent with the changing trend of lignin content in pear fruit, and the transcript level was also higher in 'Dangshan Su' pear with higher lignin and stone cell contents. Subcellular localization results showed that PbUGT72AJ2 was located mainly in the cytomembrane and cytoplasm. Based on our study, PbUGT72AJ2 is considered to be a monolignol glycosylation-related UGT. Our results provide an important source for the identification of UGTs and a foundation for the future understanding and manipulation of lignin metabolism and stone cell formation in pear fruit.

Genome-Wide Identification and Characterization of DIR Genes in Medicago truncatula

Dirigent proteins (DIRs) are critically involved in the formation of lignans, a diverse and widely distributed class of secondary plant metabolites exhibiting interesting pharmacological activities and implicated in natural plant defense. However, no detailed information is available about DIR gene family in Medicago truncatula. In this study, a total of 45 DIR genes were identified in M. truncatula. DIR proteins have variability in sequence. Most MtDIR genes have no intron. All MtDIR proteins contain single dirigent domain. A large number of MtDIR genes were expanded via gene duplication, and 37 MtDIR genes were duplicated in tandem. Digital expression data showed that 40% MtDIR genes had a higher expression level in the root. Analysis of RNA-seq and microarray data indicated that more than 30% MtDIR genes were responsive to biotic and/or abiotic treatments. This study will facilitate further studies on DIR family and provide useful clues for functional validation of DIR genes in higher plants.

Molecular identification, phylogenomic characterization and expression patterns analysis of the LIM (LIN-11, Isl1 and MEC-3 domains) gene family in pear (Pyrus bretschneideri) reveal its potential role in lignin metabolism

Lignin is the main component of stone cells, which are a key factor in determining pear quality. Therefore, modification of lignin biosynthesis has important implications for regulating stone cell formation. LIMs are involved in plant development, stress response and metabolism. However, there is still a lack of knowledge about the pear LIM family and lignin-related LIMs. To address this problem, we identified 14 LIMs from the pear genome and named them. Phylogenomic and feature domain analysis showed that they were divided into CRP- and DA&DAR-LIM groups and five subclades. LIMs from the genomes of four rosids (Prunus mummer, Prunus persica, Fragaria vesca and Vitis vinifera) were also identified, and microsynteny analysis revealed the most orthologous gene pairs in the cross of pear/grape and pear/mei. The transcript levels of PbLIMs were significantly affected by SA, ABA and MeJA. Spatio-temporal expression analysis showed that PbLIMs of the δLIM2 subfamily were highly expressed in the flowers. Changes in the expression levels of PbWLIM1a and PbWLIM1b during fruit development was consistent with the changes in lignin content. Combining phylogenetic analyses, protein three-dimensional structure determination and sequence alignment analyses, these two genes were suggested as lignin-related PbLIMs. Subcellular localization results showed that PbWLIM1a and PbWLIM1b were located mainly in the chloroplast. This study lays the foundation for revealing the mechanism of LIM-mediated lignin metabolism to regulate stone cell formation.