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

Ectopic and transient expression of VvDIR4 gene in Arabidopsis and grapes enhances resistance to anthracnose via affecting hormone signaling pathways and lignin production

BACKGROUND

DIR (Dirigent) proteins play important roles in the biosynthesis of lignin and lignans and are involved in various processes such as plant growth, development, and stress responses. However, there is less information about VvDIR proteins in grapevine (Vitis vinifera L).

RESULTS

In this study, we used bioinformatics methods to identify members of the DIR gene family in grapevine and identified 18 VvDIR genes in grapevine. These genes were classified into 5 subfamilies based on phylogenetic analysis. In promoter analysis, various plant hormones, stress, and light-responsive cis-elements were detected. Expression profiling of all genes following Colletotrichum gloeosporioides infection and phytohormones (salicylic acid (SA) and jasmonic acid (JA)) application suggested significant upregulation of 17 and 6 VvDIR genes, respectively. Further, we overexpressed the VvDIR4 gene in Arabidopsis thaliana and grapes for functional analysis. Ectopic expression of VvDIR4 in A. thaliana and transient expression in grapes increased resistance against C. gloeosporioides and C. higginsianum, respectively. Phenotypic observations showed small disease lesions in transgenic plants. Further, the expression patterns of genes having presumed roles in SA and JA signaling pathways were also influenced. Lignin contents were measured before and after C. higginsianum infection; the transgenic A. thaliana lines showed higher lignin content than wild-type, and a significant increase was observed after C. higginsianum infection.

CONCLUSIONS

Based on the findings, we surmise that VvDIR4 is involved in hormonal and lignin synthesis pathways which regulate resistance against anthracnose. Our study provides novel insights into the function of VvDIR genes and new candidate genes for grapevine disease resistance breeding programs.

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.

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.

Genome-wide identification and characterization of LcCCR13 reveals its potential role in lignin biosynthesis in Liriodendron chinense

Introduction

Wood formation is closely related to lignin biosynthesis. Cinnamoyl-CoA reductase (CCR) catalyzes the conversion of cinnamoyl-CoA to cinnamaldehydes, which is the initiation of the lignin biosynthesis pathway and a crucial point in the manipulation of associated traits. Liriodendron chinense is an economically significant timber tree. Nevertheless, the underlying mechanism of wood formation in it remains unknown; even the number of LcCCR family members in this species is unclear.

Materials and Results

This study aimed to perform a genome-wide identification of genes(s) involved in lignin biosynthesis in L. chinense via RT-qPCR assays and functional verification. Altogether, 13 LcCCR genes were identified that were divided into four major groups based on structural and phylogenetic features. The gene structures and motif compositions were strongly conserved between members of the same groups. Subsequently, the expression patterns analysis based on RNA-seq data indicated that LcCCR5/7/10/12/13 had high expression in the developing xylem at the stem (DXS). Furthermore, the RT-qPCR assays showed that LcCCR13 had the highest expression in the stem as compared to other tissues. Moreover, the overexpression of the LcCCR13 in transgenic tobacco plants caused an improvement in the CCR activity and lignin content, indicating that it plays a key role in lignin biosynthesis in the stems.

Discussion

Our research lays a foundation for deeper investigation of the lignin synthesis and uncovers the genetic basis of wood formation in L. chinense.

Dirigent protein subfamily function and structure in terrestrial plant phenol metabolism

Aquatic plant transition to land, and subsequent terrestrial plant species diversification, was accompanied by the emergence and massive elaboration of plant phenol chemo-diversity. Concomitantly, dirigent protein (DP) and dirigent-like protein subfamilies, derived from large multigene families, emerged and became extensively diversified. DP biochemical functions as gateway entry points into new and diverse plant phenol skeletal types then markedly expanded. DPs have at least eight non-uniformly distributed subfamilies, with different DP subfamily members of known biochemical/physiological function now implicated as gateway entries to lignan, lignin, aromatic diterpenoid, pterocarpan and isoflavene pathways. While some other DP subfamily members have jacalin domains, both these and indeed the majority of DPs throughout the plant kingdom await discovery of their biochemical roles. Methods and approaches were developed to discover DP biochemical function as gateway entry points to distinct plant phenol skeletal types in land plants. Various DP 3D X-ray structural determinations enabled structure-based comparative sequence analysis and modeling to understand similarities and differences among the different DP subfamilies. We consider that the core DP β-barrel fold and associated characteristics are likely common to all DPs, with several residues conserved and nearly invariant. There is also considerable variation in residue composition and topography of the putative substrate binding pockets, as well as substantial differences in several loops, such as the β1-β2 loop. All DPs likely bind and stabilize quinone methide intermediates, while guiding distinctive regio- and/or stereo-chemical entry into Nature's chemo-diverse land plant phenol metabolic classes.

The barley DIR gene family: An expanded gene family that is involved in stress responses

Gene family expansion plays a central role in adaptive divergence and, ultimately, speciation is influenced by phenotypic diversity in different environments. Barley (Hordeum vulgare) is the fourth most important cereal crop in the world and is used for brewing purposes, animal feed, and human food. Systematic characterization of expanded gene families is instrumental in the research of the evolutionary history of barley and understanding of the molecular function of their gene products. A total of 31,750 conserved orthologous groups (OGs) were identified using eight genomes/subgenomes, of which 1,113 and 6,739 were rapidly expanded and contracted OGs in barley, respectively. Five expanded OGs containing 20 barley dirigent genes (HvDIRs) were identified. HvDIRs from the same OG were phylogenetically clustered with similar gene structure and domain organization. In particular, 7 and 5 HvDIRs from OG0000960 and OG0001516, respectively, contributed greatly to the expansion of the DIR-c subfamily. Tandem duplication was the driving force for the expansion of the barley DIR gene family. Nucleotide diversity and haplotype network analysis revealed that the expanded HvDIRs experienced severe bottleneck events during barley domestication, and can thus be considered as potential domestication-related candidate genes. The expression profile and co-expression network analysis revealed the critical roles of the expanded HvDIRs in various biological processes, especially in stress responses. HvDIR18, HvDIR19, and HvDIR63 could serve as excellent candidates for further functional genomics studies to improve the production of barley products. Our study revealed that the HvDIR family was significantly expanded in barley and might be involved in different developmental processes and stress responses. Thus, besides providing a framework for future functional genomics and metabolomics studies, this study also identified HvDIRs as candidates for use in improving barley crop resistance to biotic and abiotic stresses.

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.

Mungbean DIRIGENT Gene Subfamilies and Their Expression Profiles Under Salt and Drought Stresses

DIRIGENT (DIR) genes are key players in environmental stress responses that have been identified in many vascular plant species. However, few studies have examined the VrDIR genes in mungbean. In this study, we characterized 37 VrDIR genes in mungbean using a genome-wide identification method. VrDIRs were distributed on seven of the 11 mungbean chromosomes, and chromosome three contained the most VrDIR genes, with seven members. Thirty-two of the 37 VrDIRs contained a typical DIR gene structure, with one exon; the conserved DIR domain (i.e., Pfam domain) occupied most of the protein in 33 of the 37 VrDIRs. The gene structures of VrDIR genes were analyzed, and a total of 19 distinct motifs were detected. VrDIR genes were classified into five groups based on their phylogenetic relationships, and 13 duplicated gene pairs were identified. In addition, a total of 92 cis-acting elements were detected in all 37 VrDIR promoter regions, and VrDIR genes contained different numbers and types of cis-acting elements. As a result, VrDIR genes showed distinct expression patterns in different tissues and in response to salt and drought stress.

Conservation and divergence of the TaSOS1 gene family in salt stress response in wheat (Triticum aestivum L.)

Salinity is one of the most important problems that adversely affect crops growth, productivity and quality worldwide. Salt Overly Sensitive 1 (SOS1) gene family plays vital roles in plant response to salt stress. Herein, we report the identification of the SOS family in wheat and the exploration of the expression profiles of SOSs under salt stress. Complete genome sequences of T. aestivum were downloaded from Ensembl plant database. Conservation and divergence of TaSOS1 family were conducted by using phylogenetic tree, gene structure and synteny distribution analysis. Expression profiles of TaSOS1s were obtained based on transcriptome and qRT-PCR analysis. Totally, 119 TaSOS1 proteins in wheat were identified at the genome-wide level and classified into three groups. Six motifs were conserved in TaSOS1 gene family. Moreover, 25 TaSOS1 genes had three copies distributing in three sub-genomes (A, B and D). A total of 32, 28 and 29 TaSOS1 genes were located on the sub-genomes A, B and D, respectively. Moreover, there were 19, 12, 6, 7, 28, 5 and 12 genes located on the three homologous of chromosomes 1, 2, 3, 4, 5, 6 and 7, respectively. Two genes were mapped to unattributed scaffolds. The duplication events analysis indicated that tandem repeats contributed to the expansion of the SOS1 family in wheat. Collinearity analysis demonstrated that segmental duplications play an important role in the expansion of SOS1 members. Chromosome 7, 5, 3, and 2 showed collinear relationship. Tissue specific expression pattern analysis revealed that 41 TaSOS1 genes expressed in various tissues, such as root, shoot, leaf, spike and grain. Transcriptomic analysis revealed that 28 and 26 genes were up- and down-regulated under salinity stress, respectively, of which 18 genes were further confirmed by RT-qPCR. The plants with high expression level of these genes displayed higher tolerance to salinity stress, stronger root system, higher Fv/Fm value and water potential. The results could be helpful for further elucidating the molecular mechanism of TaSOS1 related to salt tolerance in wheat and provide a toolkit for improving the salinity tolerance of wheat.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01009-y.

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.

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.

Functional characterization of soybean (Glycine max) DIRIGENT genes reveals an important role of GmDIR27 in the regulation of pod dehiscence

DIRIGENT (DIR) genes play important roles in regulating plant growth and development and have been studied in many plant species. However, information on DIR genes in soybean is limited. Here, we identified and characterized 54 GmDIRs and studied the characteristics of GmDIRs. Most of the GmDIRs contained a classical gene structure, one exon; 26 conserved motifs were found among these GmDIRs. The GmDIRs were grouped into four subfamilies, DIR-a, DIR-b, DIR-e and DIR-f, based on a phylogenetic analysis, and 24 duplicated gene pairs were identified. Differences in the cis-acting elements in the GmDIR promoter regions might result in distinct expression patterns of GmDIRs in different tissues. In addition, GmDIR27 had a close relationship with the pod dehiscence gene GmPdh1, and overexpression of GmDIR27 increased pod dehiscence by affecting several pod dehiscence-related gene expressions. Generally, our results provide essential information that aids future efforts to functionally characterize soybean GmDIR genes.

Genome-wide identification and functional characterization of the Camelina sativa WRKY gene family in response to abiotic stress

Background

WRKY transcription factors are a superfamily of regulators involved in diverse biological processes and stress responses in plants. However, there is limited knowledge about the WRKY family in camelina (Camelina sativa), an important Brassicaceae oil crop with strong tolerance for various stresses. Here, a genome-wide characterization of WRKY proteins is performed to examine their gene structures, phylogenetics, expression, conserved motif organizations, and functional annotation to identify candidate WRKYs that mediate stress resistance regulation in camelinas.

Results

A total of 242 CsWRKY proteins encoded by 224 gene loci distributed unevenly over the chromosomes were identified, and they were classified into three groups by phylogenetic analysis according to their WRKY domains and zinc finger motifs. The 15 CsWRKY gene loci generated 33 spliced variants. Orthologous WRKY gene pairs were identified, with 173 pairs in the C. sativa and Arabidopsis genomes as well as 282 pairs in the C. sativa and B. napus genomes, respectively. A total of 137 segmental duplication events were observed, but there was no tandem duplication in the camelina genome. Ten major conserved motifs were examined, with WRKYGQK being the most conserved, and several variants were present in many CsWRKYs. Expression analysis revealed that 50% more CsWRKY genes were expressed constitutively, and a set of them displayed tissue-specific expression. Notably, 11 CsWRKY genes exhibited significant expression changes in seedlings under cold, salt, and drought stresses, showing a preferentially inducible expression pattern in response to the stress.

Conclusions

The present article describes a detailed analysis of the CsWRKY gene family and its expression profiles in 12 tissues and under several stress conditions. Segmental duplication is the major force underlying the broad expansion of this gene family, and a strong purifying pressure occurred for CsWRKY proteins during their evolution. CsWRKY proteins play important roles in plant development, with differential functions in different tissues. Exceptionally, eleven CsWRKYs, particularly five alternative spliced isoforms, were found to be the possible key players in mediating plant responses to various stresses. Overall, our results provide a foundation for understanding the roles of CsWRKYs and the precise mechanism through which CsWRKYs regulate high stress resistance as well as the development of stress tolerance cultivars among Cruciferae crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07189-3.

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.

Genome-wide identification and characterization of SRO gene family in wheat: Molecular evolution and expression profiles during different stresses

SRO (SIMILAR TO RCD ONE), a type of plant-specific small protein family, play important roles in plant growth and development, as well as in response to biotic/abiotic stresses. Although characterization of SROs have been performed in model plants, little is known about their function in wheat, especially under stress conditions. In this study, 30 SRO genes were identified from the wheat genome (TaSRO). They were phylogenetically separated into two groups with distinct structures. The cis-regulatory elements in the promoter region of TaSROs were analyzed and numerous elements functionally associated with stress responding and hormones were interpreted, implying the reason for induction expression patterns of TaSROs during abiotic and biotic stresses in wheat. Whole-genome replication events in the SRO gene family of wheat and seven other species (Arabidopsis thaliana, rice, maize, barley, soybean, upland cotton, and cucumber) were analyzed, resulting in 1, 12, 9, 23, 6, 5, and 3 of gene pairs, respectively. The tissue-specific expression pattern profiling revealed that most TaSROs are highly expressed in one or more tissues and may play an important role in wheat growth and development. In addition, qRT-PCR results further confirmed that these TaSRO genes are involved in wheat stress response. In summary, our study laid a theoretical basis for molecular function deciphering of TaSROs, especially in plant hormones and biotic/abiotic stress responses.

RNA Sequencing-Associated Study Identifies GmDRR1 as Positively Regulating the Establishment of Symbiosis in Soybean

In soybean (Glycine max)-rhizobium interactions, the type III secretion system (T3SS) of rhizobium plays a key role in regulating host specificity. However, the lack of information on the role of T3SS in signaling networks limits our understanding of symbiosis. Here, we conducted an RNA sequencing analysis of three soybean chromosome segment substituted lines, one female parent and two derived lines with different chromosome-substituted segments of wild soybean and opposite nodulation patterns. By analyzing chromosome-linked differentially expressed genes in the substituted segments and quantitative trait loci (QTL)-assisted selection in the substituted-segment region, genes that may respond to type III effectors to mediate plant immunity-related signaling were identified. To narrow down the number of candidate genes, QTL assistant was used to identify the candidate region consistent with the substituted segments. Furthermore, one candidate gene, GmDRR1, was identified in the substituted segment. To investigate the role of GmDRR1 in symbiosis establishment, GmDRR1-overexpression and RNA interference soybean lines were constructed. The nodule number increased in the former compared with wild-type soybean. Additionally, the T3SS-regulated effectors appeared to interact with the GmDDR1 signaling pathway. This finding will allow the detection of T3SS-regulated effectors involved in legume-rhizobium interactions.

Genome-Wide Identification and Analysis of the Growth-Regulating Factor (GRF) Gene Family and GRF-Interacting Factor Family in Triticum aestivum L

Growth-regulating factors (GRFs) are unique transcription factors in plants. GRFs can interact with SNH (SYT N-terminal homology) domains in GRF-interacting factor (GIF) proteins via the N-terminal QLQ (Gln, Leu, Gln) domain to form functional complexes and participate in the regulation of downstream gene expression. In this study, we systematically identified the GRF gene family and GIF gene family in wheat and its relatives comprising Triticum urartu, Triticum dicoccoides, and Aegilops tauschii. Thirty GRF gene members are present in wheat, which are distributed on 12 chromosomes and they have 2-5 protein-coding regions. They all contain QLQ and WRC (Trp, Arg, Cys) conserved domains. Wheat possesses only eight members of the GIF gene family, which are distributed on six chromosomes. All wheat GIF (TaGIF) proteins have highly conserved SNH and QG (Gln, Gly) domains. The wheat GRF (TaGRF) gene family has 13 pairs of segmental duplication genes and no tandem duplication genes; the TaGIF gene family has two pairs of segmental duplication genes and no tandem duplication genes. It is speculated that segmental duplication events may be the main reason for the amplification of TaGRF gene family and TaGIF gene family. Based on published transcriptome data and qRT-PCR results of 8 TaGRF genes and 4 TaGIF genes, all of the genes responded strongly to osmotic stress, and the expression levels of TaGRF21 and TaGIF5 were also significantly upregulated under drought and cold stress conditions. The results obtained in this study may facilitate further investigations of the functions of TaGRF genes and TaGIF genes in order to identify candidate genes for use in stress-resistant wheat breeding programs.

Insights into Lignan Composition and Biosynthesis in Stinging Nettle (Urtica dioica L.)

Stinging nettle (Urtica dioica L.) has been used as herbal medicine to treat various ailments since ancient times. The biological activity of nettle is chiefly attributed to a large group of phenylpropanoid dimers, namely lignans. Despite the pharmacological importance of nettle lignans, there are no studies addressing lignan biosynthesis in this plant. We herein identified 14 genes encoding dirigent proteins (UdDIRs) and 3 pinoresinol-lariciresinol reductase genes (UdPLRs) in nettle, which are two gene families known to be associated with lignan biosynthesis. Expression profiling of these genes on different organs/tissues revealed a specific expression pattern. Particularly, UdDIR7, 12 and 13 displayed a remarkable high expression in the top internode, fibre tissues of bottom internodes and roots, respectively. The relatively high expression of UdPLR1 and UdPLR2 in the young internodes, core tissue of bottom internode and roots is consistent with the high accumulation of lariciresinol and secoisolariciresinol in these tissues. Lignan quantification showed a high abundance of pinoresinol in roots and pinoresinol diglucosides in young internodes and leaves. This study sheds light on lignan composition and biosynthesis in nettle, providing a good basis for further functional analysis of DIRs and PLRs and, ultimately, engineering lignan metabolism in planta and in cell cultures.

The Dynamics of the Cell Wall Proteome of Developing Alfalfa Stems

In this study, the cell-wall-enriched subproteomes at three different heights of alfalfa stems were compared. Since these three heights correspond to different states in stem development, a view on the dynamics of the cell wall proteome during cell maturation is obtained. This study of cell wall protein-enriched fractions forms the basis for a description of the development process of the cell wall and the linking cell wall localized proteins with the evolution of cell wall composition and structure. The sequential extraction of cell wall proteins with CaCl2, EGTA, and LiCl-complemented buffers was combined with a gel-based proteome approach and multivariate analysis. Although the highest similarities were observed between the apical and intermediate stem regions, the proteome patterns are characteristic for each region. Proteins that bind carbohydrates and have proteolytic activity, as well as enzymes involved in glycan remobilization, accumulate in the basal stem region. Beta-amylase and ferritin likewise accumulate more in the basal stem segment. Therefore, remobilization of nutrients appears to be an important process in the oldest stem segment. The intermediate and apical regions are sites of cell wall polymer remodeling, as suggested by the high abundance of proteins involved in the remodeling of the cell wall, such as xyloglucan endoglucosylase, beta-galactosidase, or the BURP-domain containing polygalacturonase non-catalytic subunit. However, the most striking change between the different stem parts is the strong accumulation of a DUF642-conserved domain containing protein in the apical region of the stem, which suggests a particular role of this protein during the early development of stem tissues.