Genome-wide identification and characterization of laccase gene family in Citrus sinensis

Genome-wide analysis of the common bean (Phaseolus vulgaris) laccase gene family and its functions in response to abiotic stress

BACKGROUND

Laccase (LAC) gene family plays a pivotal role in plant lignin biosynthesis and adaptation to various stresses. Limited research has been conducted on laccase genes in common beans.

RESULTS

29 LAC gene family members were identified within the common bean genome, distributed unevenly in 9 chromosomes. These members were divided into 6 distinct subclades by phylogenetic analysis. Further phylogenetic analyses and synteny analyses indicated that considerable gene duplication and loss presented throughout the evolution of the laccase gene family. Purified selection was shown to be the major evolutionary force through Ka / Ks. Transcriptional changes of PvLAC genes under low temperature and salt stress were observed, emphasizing the regulatory function of these genes in such conditions. Regulation by abscisic acid and gibberellins appears to be the case for PvLAC3, PvLAC4, PvLAC7, PvLAC13, PvLAC14, PvLAC18, PvLAC23, and PvLAC26, as indicated by hormone induction experiments. Additionally, the regulation of PvLAC3, PvLAC4, PvLAC7, and PvLAC14 in response to nicosulfuron and low-temperature stress were identified by virus-induced gene silence, which demonstrated inhibition on growth and development in common beans.

CONCLUSIONS

The research provides valuable genetic resources for improving the resistance of common beans to abiotic stresses and enhance the understanding of the functional roles of the LAC gene family.

Exploring Citrus sinensis Phytochemicals as Potential Inhibitors for Breast Cancer Genes BRCA1 and BRCA2 Using Pharmacophore Modeling, Molecular Docking, MD Simulations, and DFT Analysis

BACKGROUND

Structure-activity relationship (SAR) is considered to be an effective in silico approach when discovering potential antagonists for breast cancer due to gene mutation. Major challenges are faced by conventional SAR in predicting novel antagonists due to the discovery of diverse antagonistic compounds. Methodologyand Results: In predicting breast cancer antagonists, a multistep screening of phytochemicals isolated from the seeds of the Citrus sinensis plant was applied using feasible complementary methodologies. A three-dimensional quantitative structure-activity relationship (3D-QSAR) model was developed through the Flare project, in which conformational analysis, pharmacophore generation, and compound alignment were done. Ten hit compounds were obtained through the development of the 3D-QSAR model. For exploring the mechanism of action of active compounds against cocrystal inhibitors, molecular docking analysis was done through Molegro software (MVD) to identify lead compounds. Three new proteins, namely, 1T15, 3EU7, and 1T29, displayed the best Moldock scores. The quality of the docking study was assessed by a molecular dynamics simulation. Based on binding affinities to the receptor in the docking studies, three lead compounds (stigmasterol P8, epoxybergamottin P28, and nobiletin P29) were obtained, and they passed through absorption, distribution, metabolism, and excretion (ADME) studies via the SwissADME online service, which proved that P28 and P29 were the most active allosteric inhibitors with the lowest toxicity level against breast cancer. Then, density functional theory (DFT) studies were performed to measure the active compound's reactivity, hardness, and softness with the help of Gaussian 09 software.

CONCLUSIONS

This multistep screening of phytochemicals revealed high-reliability antagonists of breast cancer by 3D-QSAR using flare, docking analysis, and DFT studies. The present study helps in providing a proper guideline for the development of novel inhibitors of BRCA1 and BRCA2.

Genome-wide identification and expression patterns of the laccase gene family in response to kiwifruit bacterial canker infection

BACKGROUND

Kiwifruit bacterial canker, caused by Pseudomonas syringae pv. actinidiae (Psa), is a destructive disease worldwide. Resistance genes that respond to Psa infection urgently need to be identified for controlling this disease. Laccase is mainly involved in the synthesis of lignin in the plant cell wall and plays a prominent role in plant growth and resistance to pathogen infection. However, the role of laccase in kiwifruit has not been reported, and whether laccase is pivotal in the response to Psa infection remains unclear.

RESULTS

We conducted a bioinformatics analysis to identify 55 laccase genes (AcLAC1-AcLAC55) in the kiwifruit genome. These genes were classified into five cluster groups (I-V) based on phylogenetic analysis, with cluster groups I and II having the highest number of members. Analysis of the exon-intron structure revealed that the number of exons varied from 1 to 8, with an average of 5 introns. Our evolutionary analysis indicated that fragment duplication played a key role in the expansion of kiwifruit laccase genes. Furthermore, evolutionary pressure analysis suggested that AcLAC genes were under purifying selection. We also performed a cis-acting element analysis and found that AcLAC genes contained multiple hormone (337) and stress signal (36) elements in their promoter regions. Additionally, we investigated the expression pattern of laccase genes in kiwifruit stems and leaves infected with Psa. Our findings revealed that laccase gene expression levels in the stems were higher than those in the leaves 5 days after inoculation with Psa. Notably, AcLAC2, AcLAC4, AcLAC17, AcLAC18, AcLAC26, and AcLAC42 showed significantly higher expression levels (p < 0.001) compared to the non-inoculated control (0 d), suggesting their potential role in resisting Psa infection. Moreover, our prediction indicated that 21 kiwifruit laccase genes are regulated by miRNA397, they could potentially act as negative regulators of lignin biosynthesis.

CONCLUSIONS

These results are valuable for further analysis of the resistance function and molecular mechanism of laccases in kiwifruit.

Contribution of the LAC Genes in Fruit Quality Attributes of the Fruit-Bearing Plants: A Comprehensive Review

Laccase genes produce laccase enzymes that play a crucial role in the production of lignin and oxidation reactions within plants. Lignin is a complex polymer that provides structure and toughness to the cell walls of numerous fruit plants. The LAC genes that encode laccase enzymes play vital roles in plant physiology, including the synthesis of pigments like PA that contribute to the colors of fruits, and in defending against pathogens and environmental stresses. They are crucial for fruit development, ripening, structural maintenance in plants, and adaptation to various environmental factors. As such, these genes and enzymes are essential for plant growth and development, as well as for various biotechnological applications in environmental remediation and industrial processes. This review article emphasizes the significance of genes encoding laccase enzymes during fruit growth, specifically pertaining to the strengthening of the endocarp through lignification. This process is crucial for ensuring fruit defense and optimizing seed scattering. The information gathered in this article will aid breeders in producing future fruit-bearing plants that are resistant to disease, cost-effective, and nutrient-rich.

Comprehensive analysis of the laccase gene family in tea plant highlights its roles in development and stress responses

BACKGROUND

Laccase (LAC) is the pivotal enzyme responsible for the polymerization of monolignols and stress responses in plants. However, the roles of LAC genes in plant development and tolerance to diverse stresses are still largely unknown, especially in tea plant (Camellia sinensis), one of the most economically important crops worldwide.

RESULTS

In total, 51 CsLAC genes were identified, they were unevenly distributed on different chromosomes and classified into six groups based on phylogenetic analysis. The CsLAC gene family had diverse intron-exon patterns and a highly conserved motif distribution. Cis-acting elements in the promoter demonstrated that promoter regions of CsLACs encode various elements associated with light, phytohormones, development and stresses. Collinearity analysis identified some orthologous gene pairs in C. sinensis and many paralogous gene pairs among C. sinensis, Arabidopsis and Populus. Tissue-specific expression profiles revealed that the majority of CsLACs had high expression in roots and stems and some members had specific expression patterns in other tissues, and the expression patterns of six genes by qRT‒PCR were highly consistent with the transcriptome data. Most CsLACs showed significant variation in their expression level under abiotic (cold and drought) and biotic (insect and fungus) stresses via transcriptome data. Among them, CsLAC3 was localized in the plasma membrane and its expression level increased significantly at 13 d under gray blight treatment. We found that 12 CsLACs were predicted to be targets of cs-miR397a, and most CsLACs showed opposite expression patterns compared to cs-miR397a under gray blight infection. Additionally, 18 highly polymorphic SSR markers were developed, these markers can be widely used for diverse genetic studies of tea plants.

CONCLUSIONS

This study provides a comprehensive understanding of the classification, evolution, structure, tissue-specific profiles, and (a)biotic stress responses of CsLAC genes. It also provides valuable genetic resources for functional characterization towards enhancing tea plant tolerance to multiple (a)biotic stresses.

Genome-wide analysis of the laccase (LAC) gene family in Aeluropus littoralis: A focus on identification, evolution and expression patterns in response to abiotic stresses and ABA treatment

Laccases are plant enzymes with essential functions during growth and development. These monophenoloxidases are involved in lignin polymerization, and their expression respond to environmental stress. However, studies of laccases in some plants and fungi have highlighted that many structural and functional aspects of these genes are still unknown. Here, the laccase gene family in Aeluropus littoralis (AlLAC) is described based on sequence structure and expression patterns under abiotic stresses and ABA treatment. Fifteen non-redundant AlLACs were identified from the A. littoralis genome, which showed differences in physicochemical characteristics and gene structure. Based on phylogenetic analysis, AlLACs and their orthologues were classified into five groups. A close evolutionary relationship was observed between LAC gene family members in rice and A. littoralis. According to the interaction network, AlLACs interact more with proteins involved in biological processes such as iron incorporation into the metallo-sulfur cluster, lignin catabolism, regulation of the symbiotic process and plant-type primary cell wall biogenesis. Gene expression analysis of selected AlLACs using real-time RT (reverse transcription)-PCR revealed that AlLACs are induced in response to abiotic stresses such as cold, salt, and osmotic stress, as well as ABA treatment. Moreover, AlLACs showed differential expression patterns in shoot and root tissues. Our findings indicate that AlLACs are preferentially involved in the late response of A. littoralis to abiotic stress.

Characterization of plant laccase genes and their functions

Laccase is a copper-containing polyphenol oxidase found in different organisms. The multigene family that encodes laccases is widely distributed in plant genomes. Plant laccases oxidize monolignols to produce lignin which is important for plant growth and stress responses. Industrial applications of fungal and bacterial laccases are extensively explored and addressed. Recently many studies have focused on the significance of plant laccase, particularly in crop yield, and its functions in different environmental conditions. This review summarizes the transcriptional and posttranscriptional regulation of plant laccase genes and their functions in plant growth and development. It especially describes the responses of laccase genes to various stresses and their contributions to plant biotic and abiotic stress resistance. In-depth explanations and scientific advances will serve as foundations for research into plant laccase genes' function, mechanism, and possible applications.

Genome-wide analysis of the laccase gene family in wheat and relationship with arbuscular mycorrhizal colonization

We identified 156 laccase genes belonging to 11 subfamilies in the wheat genome, and the natural variation of laccase genes significantly affected the development of wheat-arbuscular mycorrhizal symbiosis. Laccases (LACs) have a variety of functions in plant lignification, cell elongation and stress responses. This study aimed to reveal the phylogeny, chromosomal spatial distribution, coexpression and evolution of LAC genes in the wheat genome and to investigate the possible roles of LAC genes during arbuscular mycorrhizal (AM) symbiosis. The genomic characteristics of LAC genes were analyzed by using bioinformatics analysis methods, and the polymorphisms of LAC genes were analyzed by using a diverse wheat panel composed of 289 wheat cultivars. We identified 156 LAC genes belonging to 11 subfamilies in the wheat genome, and segmental duplication dominated the amplification of the LAC gene family in the wheat genome. LACs are dominantly located in the R2 region of wheat chromosomes. Some LACs are collinear with the characterized LACs in Arabidopsis thaliana or rice. A number of genes encoding transcription factors, kinases, and phosphatases were coexpressed with LAC genes in wheat. TaLACs may be potential targets for some miRNAs. Most TaLACs are mainly expressed in the roots and stems of plants. The expression of TaLACs could be regulated by the inoculation of Fusarium graminearum or AM fungi. The polymorphisms of TaLACs mainly accumulate by random drift instead of by selection. Through candidate gene association analysis, we found that the natural variations in TaLACs significantly affected root colonization by AM fungi. The present study provides useful information for further study of the biological functions of LAC genes in wheat, especially the roles of LAC genes during the development of AM symbiosis.

The Citrus Laccase Gene CsLAC18 Contributes to Cold Tolerance

Plant laccases, as multicopper oxidases, play an important role in monolignol polymerization, and participate in the resistance response of plants to multiple biotic/abiotic stresses. However, little is currently known about the role of laccases in the cold stress response of plants. In this study, the laccase activity and lignin content of C. sinensis leaves increased after the low-temperature treatment, and cold treatment induced the differential regulation of 21 CsLACs, with 15 genes being upregulated and 6 genes being downregulated. Exceptionally, the relative expression level of CsLAC18 increased 130.17-fold after a 48-h treatment. The full-length coding sequence of CsLAC18 consists of 1743 nucleotides and encodes a protein of 580 amino acids, and is predominantly expressed in leaves and fruits. CsLAC18 was phylogenetically related to AtLAC17, and was localized in the cell membrane. Overexpression of CsLAC18 conferred enhanced cold tolerance on transgenic tobacco; however, virus-induced gene silencing (VIGS)-mediated suppression of CsLAC18 in Poncirus trifoliata significantly impaired resistance to cold stress. As a whole, our findings revealed that CsLAC18 positively regulates a plant's response to cold stress, providing a potential target for molecular breeding or gene editing.

Genome-Wide Identification of Switchgrass Laccases Involved in Lignin Biosynthesis and Heavy-Metal Responses

Plant laccase genes belong to a multigene family, play key roles in lignin polymerization, and participate in the resistance of plants to biotic and abiotic stresses. Switchgrass is an important resource for forage and bioenergy production, yet information about the switchgrass laccase gene family is scarce. Using bioinformatic approaches, a genome-wide analysis of the laccase multigene family in switchgrass was carried out in this study. In total, 49 laccase genes (PvLac1 to PvLac49) were identified; these can be divided into five subclades, and 20 of them were identified as targets of miR397. The tandem and segmental duplication of laccase genes on Chr05 and Chr08 contributed to the expansion of the laccase family. The laccase proteins shared conserved signature sequences but displayed relatively low sequence similarity, indicating the potential functional diversity of switchgrass laccases. Switchgrass laccases exhibited distinct tissue/organ expression patterns, revealing that some laccases might be involved in the lignification process during stem development. All five of the laccase isoforms selected from different subclades responded to heavy metal. The immediate response of lignin-related laccases, as well as the delayed response of low-abundance laccases, to heavy-metal treatment shed light on the multiple roles of laccase isoforms in response to heavy-metal stress.

Involvement of Laccase2 in Cuticle Sclerotization of the Whitefly, Bemisia tabaci Middle East–Asia Minor 1

Cuticle sclerotization is critical for insect survival. Laccase2 (Lac2) is a phenol oxidase that plays a key role in cuticle formation and pigmentation in a variety of insects. However, the function of Lac2 in whitefly, Bemisia tabaci, remains unclear. In this study, we identified a BtLac2 gene in B. tabaci MEAM1 and found that BtLac2 was expressed in all stages. It was highly expressed in the egg stage, followed by nymph and adult. Moreover, the expression of BtLac2 was higher in the cuticle than in other tissues. Knockdown of BtLac2 in nymphs produced thinner and fragile cuticles, which significantly increased the mortality rate, extended the development duration of nymphs, and decreased the emergence rate of adults. This result demonstrates that BtLac2 plays an important role in the cuticle hardening of B. tabaci and suggests a potential management strategy using RNAi to knock down BtLac2 expression.

Integrated Analysis of Physiological, mRNA Sequencing, and miRNA Sequencing Data Reveals a Specific Mechanism for the Response to Continuous Cropping Obstacles in Pogostemon cablin Roots

Pogostemon cablin (patchouli) is a commercially important medicinal and industrial crop grown worldwide for its medicinal and aromatic properties. Patchoulol and pogostone, derived from the essential oil of patchouli, are considered valuable components in the cosmetic and pharmaceutical industries. Due to its high application value in the clinic and industry, the demand for patchouli is constantly growing. Unfortunately, patchouli cultivation has suffered due to severe continuous cropping obstacles, resulting in a significant decline in yield and quality. Moreover, the physiological and transcriptional changes in patchouli in response to continuous cropping obstacles remain unclear. This has greatly restricted the development of the patchouli industry. To explore the mechanism underlying the rapid response of patchouli roots to continuous cropping stress, integrated analysis of the transcriptome and miRNA profiles of patchouli roots under continuous and noncontinuous cropping conditions in different growth periods was conducted using RNA sequencing (RNA-seq) and miRNA-seq and complemented with physiological data. The physiological and biochemical results showed that continuous cropping significantly inhibited root growth, decreased root activity, and increased the activity of antioxidant enzymes (superoxide dismutase, peroxidase, and catalase) and the levels of osmoregulators (malondialdehyde, soluble protein, soluble sugar, and proline). Subsequently, we found 4,238, 3,494, and 7,290 upregulated and 4,176, 3,202, and 8,599 downregulated differentially expressed genes (DEGs) in the three growth periods of continuously cropped patchouli, many of which were associated with primary carbon and nitrogen metabolism, defense responses, secondary metabolite biosynthesis, and transcription factors. Based on miRNA-seq, 927 known miRNAs and 130 novel miRNAs were identified, among which 67 differentially expressed miRNAs (DEMIs) belonging to 24 miRNA families were induced or repressed by continuous cropping. By combining transcriptome and miRNA profiling, we obtained 47 miRNA-target gene pairs, consisting of 18 DEMIs and 43 DEGs, that likely play important roles in the continuous cropping response of patchouli. The information provided in this study will contribute to clarifying the intricate mechanism underlying the patchouli response to continuous cropping obstacles. In addition, the candidate miRNAs and genes can provide a new strategy for breeding continuous cropping-tolerant patchouli.

Genome-wide identification and characterization of laccase family members in eggplant (Solanum melongena L.)

Laccase, as a copper-containing polyphenol oxidase, primarily functions in the process of lignin, anthocyanin biosynthesis, and various abiotic/biotic stresses. In this study, forty-eight laccase members were identified in the eggplant genome. Only forty-two laccase genes from eggplant (SmLACs) were anchored unevenly in 12 chromosomes, the other six SmLACs were mapped on unanchored scaffolds. Phylogenetic analysis indicated that only twenty-five SmLACs were divided into six different groups on the basis of groups reported in Arabidopsis. Gene structure analysis revealed that the number of exons ranged from one to 13. Motif analysis revealed that SmLACs included six conserved motifs. In aspects of gene duplication analysis, twenty-one SmLACs were collinear with LAC genes from Arabidopsis, tomato or rice. Cis-regulatory elements analysis indicated many SmLACs may be involved in eggplant morphogenesis, flavonoid biosynthesis, diverse stresses and growth/development processes. Expression analysis further confirmed that a few SmLACs may function in vegetative and reproductive organs at different developmental stages and also in response to one or multiple stresses. This study would help to further understand and enrich the physiological function of the SmLAC gene family in eggplant, and may provide high-quality genetic resources for eggplant genetics and breeding.

Genome-Wide Comprehensive Analysis of PtLACs: Prediction and Verification of the Functional Divergence of Tandem-Duplicated Genes

Laccases (EC 1.10.3.2) have been widely considered to participate in the metabolic processes of lignin synthesis, osmotic stress response, and flavonoid oxidation in higher plants. The research into Populus trichocarpa laccase focused on the synthesis of lignin in the past few years. In this study, for the first time, a comprehensive analysis of 53 laccase copies in the P. trichocarpa genome was conducted. Positive selection analysis using the branch-site model indicated that LAC genes in terrestrial plants have undergone selective pressure for adaptive evolution. On the basis of the phylogenetic relationship, we reconstructed the evolutionary process of terrestrial plant laccase and found that this gene family began to expand during the evolution of angiosperms. Tandem duplication is the main form of expansion of the PtLAC gene family. The analysis of the sequence characteristics, gene structure, expression pattern, and gene synonymous mutation rate of PtLACs provided a theoretical basis for the functional divergence of tandem duplicated genes. The synonymous mutation rate was used to quantify the divergence time of 11 tandem duplicated gene clusters. Cluster 2, with the earliest divergence time and lower share of sequence similarity, and cluster 5, with the latest divergence time and higher share of similarity, were selected in this study to explore the functional divergence of tandem-duplicated gene clusters. Tobacco subcellular localization and Arabidopsis transgenes verified the functional differentiation of PtLAC genes in cluster 2 and the functional non-differentiation of PtLAC genes in cluster 5. The results of this study provide a reference for the functional differentiation of tandem-duplicated PtLAC.

Genome Wide Identification of Respiratory Burst Oxidase Homolog (Rboh) Genes in Citrus sinensis and Functional Analysis of CsRbohD in Cold Tolerance

Respiratory burst oxidase homologs (Rbohs) are critical enzymes involved in the generation of reactive oxygen species (ROS) that play an important role in plant growth and development as well as various biotic and abiotic stresses in plants. Thus far, there have been few reports on the characterization of the Rboh gene family in Citrus. In this study, seven Rboh genes (CsRbohA~CsRbohG) were identified in the Citrus sinensis genome. The CsRboh proteins were predicted to localize to the cell membrane. Most CsRbohs contained four conserved domains, an EF-hand domain, and a transmembrane region. Phylogenetic analysis demonstrated that the CsRbohs were divided into five groups, suggesting potential distinct functions and evolution. The expression profiles revealed that these seven CsRboh genes displayed tissue-specific expression patterns, and five CsRboh genes were responsive to cold stress. Fourteen putative cis-acting elements related to stress response, hormone response, and development regulation were present within the promoters of CsRboh genes. The in-silico microRNA target transcript analyses indicated that CsRbohE might be targeted by csi-miR164. Further functional and physiological analyses showed that the knockdown of CsRbohD in trifoliate orange impaired resistance to cold stress. As a whole, our results provide valuable information for further functional studies of the CsRboh genes in response to cold stress.

Genome-wide identification studies - A primer to explore new genes in plant species

Genome data have accumulated rapidly in recent years, doubling roughly after every 6 months due to the influx of next-generation sequencing technologies. A plethora of plant genomes are available in comprehensive public databases. This easy access to data provides an opportunity to explore genome datasets and recruit new genes in various plant species not possible a decade ago. In the past few years, many gene families have been published using these public datasets. These genome-wide studies identify and characterize gene members, gene structures, evolutionary relationships, expression patterns, protein interactions and gene ontologies, and predict putative gene functions using various computational tools. Such studies provide meaningful information and an initial framework for further functional elucidation. This review provides a concise layout of approaches used in these gene family studies and demonstrates an outline for employing various plant genome datasets in future studies.

Transcriptomic Response of Huanglongbing-Infected Citrus sinensis Following Field Application of a Microbial Fermentation Product

Huanglongbing (HLB) is considered the most destructive disease in Citrus production and threatens the future of the industry. Microbial-derived defense elicitors have gained recognition for their role in plant defense priming. This work assessed a 5% (V/V) microbial fermentation application (MFA) and its role in the elicitation of defense responses in HLB-infected Citrus sinensis trees following a foliar application with a pump sprayer. Using a PCR detection method, HLB infection levels were monitored in healthy and infected trees for 20months. Nutrient analysis assessed N, P, K, Ca, Mg, Mn, Zn, Fe, B, and Cu concentrations in the trees. MFA significantly increased Cu concentrations in treated trees and resulted in the stabilization of disease index (DI) in infected trees. Initial real-time qPCR analysis of defense-associated genes showed a significant increase in pathogenesis-related protein 2 (PR2) and phenylalanine ammonia lyase (PAL) gene expression in healthy and HLB-infected trees in response to MFA. Gene expression of PR2 and PAL peaked 6h post-microbial fermentation application during an 8-h sampling period. A transcriptomic assessment using GeneChip microarray of the hour 6 samples revealed differential expression of 565 genes when MFA was applied to healthy trees and 909 genes when applied infected citrus trees when compared to their respective controls. There were 403 uniquely differentially expressed genes in response to MFA following an intersectional analysis of both healthy and infected citrus trees. The transcriptomic analysis revealed that several genes associated with plant development, growth, and defense were upregulated in response to MFA, including multiple PR genes, lignin formation genes, ROS-related genes, hormone synthases, and hormone regulators. This study provides further evidence that MFA may play an important role as a plant elicitor in an integrated pest management strategy in citrus and other agronomically important crops.

Plant Copper Metalloenzymes As Prospects for New Metabolism Involving Aromatic Compounds

Copper is an important transition metal cofactor in plant metabolism, which enables diverse biocatalysis in aerobic environments. Multiple classes of plant metalloenzymes evolved and underwent genetic expansions during the evolution of terrestrial plants and, to date, several representatives of these copper enzyme classes have characterized mechanisms. In this review, we give an updated overview of chemistry, structure, mechanism, function and phylogenetic distribution of plant copper metalloenzymes with an emphasis on biosynthesis of aromatic compounds such as phenylpropanoids (lignin, lignan, flavonoids) and cyclic peptides with macrocyclizations via aromatic amino acids. We also review a recent addition to plant copper enzymology in a copper-dependent peptide cyclase called the BURP domain. Given growing plant genetic resources, a large pool of copper biocatalysts remains to be characterized from plants as plant genomes contain on average more than 70 copper enzyme genes. A major challenge in characterization of copper biocatalysts from plant genomes is the identification of endogenous substrates and catalyzed reactions. We highlight some recent and future trends in filling these knowledge gaps in plant metabolism and the potential for genomic discovery of copper-based enzymology from plants.

Vitis OneGenE: A Causality-Based Approach to Generate Gene Networks in Vitis vinifera Sheds Light on the Laccase and Dirigent Gene Families

The abundance of transcriptomic data and the development of causal inference methods have paved the way for gene network analyses in grapevine. Vitis OneGenE is a transcriptomic data mining tool that finds direct correlations between genes, thus producing association networks. As a proof of concept, the stilbene synthase gene regulatory network obtained with OneGenE has been compared with published co-expression analysis and experimental data, including cistrome data for MYB stilbenoid regulators. As a case study, the two secondary metabolism pathways of stilbenoids and lignin synthesis were explored. Several isoforms of laccase, peroxidase, and dirigent protein genes, putatively involved in the final oxidative oligomerization steps, were identified as specifically belonging to either one of these pathways. Manual curation of the predicted sequences exploiting the last available genome assembly, and the integration of phylogenetic and OneGenE analyses, identified a group of laccases exclusively present in grapevine and related to stilbenoids. Here we show how network analysis by OneGenE can accelerate knowledge discovery by suggesting new candidates for functional characterization and application in breeding programs.

The Laccase Gene Family Mediate Multi-Perspective Trade-Offs during Tea Plant (Camellia sinensis) Development and Defense Processes

Laccase (LAC) plays important roles in different plant development and defense processes. In this study, we identified laccase genes (CsLACs) in Camellia sinensis cv ‘Longjing43′ cultivars, which were classified into six subclades. The expression patterns of CsLACs displayed significant spatiotemporal variations across different tissues and developmental stages. Most members in subclades II, IV and subclade I exhibited contrasting expression patterns during leaf development, consistent with a trade-off model for preferential expression in the early and late developmental stages. The extensive transcriptional changes of CsLACs under different phytohormone and herbivore treatment were observed and compared, with the expression of most genes in subclades I, II and III being downregulated but genes in subclades IV, V and VI being upregulated, suggesting a growth and defense trade-off model between these subclades. Taken together, our research reveal that CsLACs mediate multi-perspective trade-offs during tea plant development and defense processes and are involved in herbivore resistance in tea plants. More in-depth research of CsLACs upstream regulation and downstream targets mediating herbivore defense should be conducted in the future.

The role of JrLACs in the lignification of walnut endocarp

BACKGROUND

The walnut shell, which is composed of a large number of sclereids originating from the lignified parenchyma of the endocarp, plays an important role in fruit development and during harvesting and storage. The physical and chemical properties of walnut shells are closely related to the lignin content. Laccase is the key enzyme responsible for lignin biosynthesis by the polymerization of monolignols and plays crucial roles in secondary cell wall formation in plants. In this study, we screened and identified laccase family genes from the walnut genome and investigated the expression of laccase during endocarp lignification in walnut.

RESULTS

A total of 37 laccase genes were screened from the walnut genome and distributed on nine chromosomes and classified into 6 subfamilies, among which subfamily IV showed distinct expansion. We observed that endocarp lignification started 44 days after flowering (DAF), and at later periods, the lignin content increased rapidly, with growth peaks at 44-50 DAF and 100-115 DAF. The lignification of the endocarp proceeded from the outside to the inside, as demonstrated by section staining in combination with endocarp staining. Furthermore, the changes in the expression of laccase family genes in the endocarp at different developmental stages were studied, and JrLACs showed different expression trends. The expression of nine genes showed significant increase after 44 DAF, and among these, JrLAC12-1, JrLAC12-2 and JrLAC16 showed a significant change in expression at the lignification stage. A study of the expression of JrLACs in different tissues and at various endocarp developmental stages revealed, that most JrLACs were expressed at low levels in mature tissues and at high levels in young tissues, in particular, JrLAC12-1 showed high expression in the young stems. A significant positive correlation was found between the expression of JrLAC12-1 and the variation in the lignin content in the endocarp.

CONCLUSION

Laccase genes play an important role in the lignification of the walnut endocarp, and JrLACs play different roles during fruit development. This study shows that JrLAC12-1 may play a key role in the lignification of endocarp.

Comparative analysis of alfalfa (Medicago sativa L.) seedling transcriptomes reveals genotype-specific drought tolerance mechanisms

Drought is one of the main abiotic factors that affect alfalfa yield. The identification of genes that control this complex trait can provide important insights for alfalfa breeding. However, little is known about how alfalfa responds and adapts to drought stress, particularly in cultivars of differing drought tolerance. In this study, the drought-tolerant cultivar Dryland 'DT' and the drought-sensitive cultivar WL343HQ 'DS' were used to characterize leaf and root physiological responses and transcriptional changes in response to water deficit. Under drought stress, Dryland roots (DTR) showed more differentially expressed genes than WL343HQ roots (DSR), whereas WL343HQ leaves (DSL) showed more differentially expressed genes than Dryland leaves (DTL). Many of these genes were involved in stress-related pathways, carbohydrate metabolism, and lignin and wax biosynthesis, which may have improved the drought tolerance of alfalfa. We also observed that several genes related to ABA metabolism, root elongation, peroxidase activity, cell membrane stability, ubiquitination, and genetic processing responded to drought stress in alfalfa. We highlighted several candidate genes, including sucrose synthase, xylan 1,4-beta-xylosidase, primary-amine oxidase, and alcohol-forming fatty acyl-CoA reductase, for future studies on drought stress resistance in alfalfa and other plant species. In summary, our results reveal the unique drought adaptation and resistance characteristics of two alfalfa genotypes. These findings, which may be valuable for drought resistance breeding, warrant further gene functional analysis to augment currently available information and to clarify the drought stress regulatory mechanisms of alfalfa and other plants.

Genome-wide identification of the restorer-of-fertility-like (RFL) gene family in Brassica napus and expression analysis in Shaan2A cytoplasmic male sterility

Background

Cytoplasmic male sterility (CMS) is very important in hybrid breeding. The restorer-of-fertility (Rf) nuclear genes rescue the sterile phenotype. Most of the Rf genes encode pentatricopeptide repeat (PPR) proteins.

Results

We investigated the restorer-of-fertility-like (RFL) gene family in Brassica napus. A total of 53 BnRFL genes were identified. While most of the BnRFL genes were distributed on 10 of the 19 chromosomes, gene clusters were identified on chromosomes A9 and C8. The number of PPR motifs in the BnRFL proteins varied from 2 to 19, and the majority of BnRFL proteins harbored more than 10 PPR motifs. An interaction network analysis was performed to predict the interacting partners of RFL proteins. Tissue-specific expression and RNA-seq analyses between the restorer line KC01 and the sterile line Shaan2A indicated that BnRFL1, BnRFL5, BnRFL6, BnRFL8, BnRFL11, BnRFL13 and BnRFL42 located in gene clusters on chromosomes A9 and C8 were highly expressed in KC01.

Conclusions

In the present study, identification and gene expression analysis of RFL gene family in the CMS system were conducted, and seven BnRFL genes were identified as candidates for the restorer genes in Shaan2A CMS. Taken together, this method might provide new insight into the study of Rf genes in other CMS systems.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07163-z.

Expression Profile of Laccase Gene Family in White-Rot Basidiomycete Lentinula edodes under Different Environmental Stresses

Laccases belong to ligninolytic enzymes and play important roles in various biological processes of filamentous fungi, including fruiting-body formation and lignin degradation. The process of fruiting-body development in Lentinula edodes is complex and is greatly affected by environmental conditions. In this paper, 14 multicopper oxidase-encoding (laccase) genes were analyzed in the draft genome sequence of L. edodes strain W1-26, followed by a search of multiple stress-related Cis-elements in the promoter region of these laccase genes, and then a transcription profile analysis of 14 laccase genes (Lelcc) under the conditions of different carbon sources, temperatures, and photoperiods. All laccase genes were significantly regulated by varying carbon source materials. The expression of only two laccase genes (Lelcc5 and Lelcc6) was induced by sodium-lignosulphonate and the expression of most laccase genes was specifically upregulated in glucose medium. Under different temperature conditions, the expression levels of most laccase genes decreased at 39 °C and transcription was significantly increased for Lelcc1, Lelcc4, Lelcc5, Lelcc9, Lelcc12, Lelcc13, and Lelcc14 after induction for 24 h at 10 °C, indicating their involvement in primordium differentiation. Tyrosinase, which is involved in melanin synthesis, was clustered with the same group as Lelcc4 and Lelcc7 in all the different photoperiod treatments. Meanwhile, five laccase genes (Lelcc8, Lelcc9, Lelcc12, Lelcc13, and Lelcc14) showed similar expression profiles to that of two blue light receptor genes (LephrA and LephrB) in the 12 h light/12 h dark treatment, suggesting the involvement of laccase genes in the adaptation process of L. edodes to the changing environment and fruiting-body formation. This study contributes to our understanding of the function of the different Lelcc genes and facilitates the screening of key genes from the laccase gene family for further functional research.

Genome-wide analysis, transcription factor network approach and gene expression profile of GH3 genes over early somatic embryogenesis in Coffea spp

Background

Coffee production relies on plantations with varieties from Coffea arabica and Coffea canephora species. The first, the most representative in terms of coffee consumption, is mostly propagated by seeds, which leads to management problems regarding the plantations maintenance, harvest and processing of grains. Therefore, an efficient clonal propagation process is required for this species cultivation, which is possible by reaching a scalable and cost-effective somatic embryogenesis protocol. A key process on somatic embryogenesis induction is the auxin homeostasis performed by Gretchen Hagen 3 (GH3) proteins through amino acid conjugation. In this study, the GH3 family members were identified on C. canephora genome, and by performing analysis related to gene and protein structure and transcriptomic profile on embryogenic tissues, we point a GH3 gene as a potential regulator of auxin homeostasis during early somatic embryogenesis in C. arabica plants.

Results

We have searched within the published C. canephora genome and found 17 GH3 family members. We checked the conserved domains for GH3 proteins and clustered the members in three main groups according to phylogenetic relationships. We identified amino acids sets in four GH3 proteins that are related to acidic amino acid conjugation to auxin, and using a transcription factor (TF) network approach followed by RT-qPCR we analyzed their possible transcriptional regulators and expression profiles in cells with contrasting embryogenic potential in C. arabica. The CaGH3.15 expression pattern is the most correlated with embryogenic potential and with CaBBM, a C. arabica ortholog of a major somatic embryogenesis regulator.

Conclusion

Therefore, one out of the GH3 members may be influencing on coffee somatic embryogenesis by auxin conjugation with acidic amino acids, which leads to the phytohormone degradation. It is an indicative that this gene can serve as a molecular marker for coffee cells with embryogenic potential and needs to be further studied on how much determinant it is for this process. This work, together with future studies, can support the improvement of coffee clonal propagation through in vitro derived somatic embryos.

Analysis of the laccase gene family and miR397-/miR408-mediated posttranscriptional regulation in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most commonly used traditional Chinese medicine materials. It contains important bioactive phenolic compounds, such as salvianolic acids, flavonoids and anthocyanins. Elucidation of phenolic compound biosynthesis and its regulatory mechanism is of great significance for S. miltiorrhiza quality improvement. Laccases (LACs) are multicopper-containing enzymes potentially involved in the polymerization of phenolic compounds. So far, little has been known about LAC genes in S. miltiorrhiza. Through systematic investigation of the whole genome sequence and transcriptomes of S. miltiorrhiza, we identified 65 full-length SmLAC genes (SmLAC1–SmLAC65). Phylogenetic analysis showed that 62 of the identified SmLACs clustered with LACs from Arabidopsis and Populus trichocarpa in seven clades (C1–C7), whereas the other three fell into one S. miltiorrhiza-specific clade (C8). All of the deduced SmLAC proteins contain four conserved signature sequences and three typical Cu-oxidase domains, and gene structures of most LACs from S. miltiorrhiza, Arabidopsis and P. trichocarpa were highly conserved, however SmLACs encoding C8 proteins showed distinct intron-exon structures. It suggests the conservation and diversity of plant LACs in gene structures. The majority of SmLACs exhibited tissue-specific expression patterns, indicates manifold functions of SmLACs played in S. miltiorrhiza. Analysis of high-throughput small RNA sequences and degradome data and experimental validation using the 5′ RACE method showed that 23 SmLACs were targets of Smi-miR397. Among them, three were also targeted by Smi-miR408. It suggests the significance of miR397 and miR408 in posttranscriptional regulation of SmLAC genes. Our results provide a foundation for further demonstrating the functions of SmLACs in the production of bioactive phenolic compounds in S. miltiorrhiza.

Physiological responses and small RNAs changes in maize under nitrogen deficiency and resupply

BACKGROUND

Maize is an important crop in the world, nitrogen stress severely reduces maize yield. Although a large number of studies have identified the expression changes of microRNAs (miRNAs) under N stress in several species, the miRNAs expression patterns of N-deficient plants under N resupply remain unclear.

OBJECTIVE

The primary objective of this study was to identify miRNAs in response to nitrogen stress and understand relevant physiological changes in nitrogen-deficient maize after nitrogen resupply.

METHODS

Physiological parameters were measured to study relevant physiological changes under different nitrogen conditions. Small RNA sequencing and qRT-PCR analysis were performed to understand the response of miRNAs under different nitrogen conditions.

RESULTS

The content of chlorophyll, soluble protein and nitrate nitrogen decreased than CK by 0.52, 0.49 and 0.82 times after N deficiency treatment and increased than ND by 0.52, 1.36 and 0.65 times after N resupply, respectively. Conversely, the activity of superoxide dismutase (SOD) and peroxidase (POD) increased by 0.67 and 1.64 times than CK after N deficiency, respectively, and decreased by 0.09 and 0.35 times than ND after N resupply. A total of 226 known miRNAs were identified by sRNA sequencing; 106 miRNAs were differentially expressed between the control and N-deficient groups, and 103 were differentially expressed between the N-deficient and N-resupply groups (P < 0.05). Real-time quantitative PCR (qPCR) was used to further validate and analyze the expression of the identified miRNAs. A total of 1609 target genes were identified by target prediction, and some differentially expressed miRNAs were predicted to target transcription factors and functional proteins. Gene Ontology (GO) analysis was used to determine the biological function of these targets and revealed that some miRNAs, such as miR169, miR1214, miR2199, miR398, miR408 and miR827 might be involved in nitrogen metabolism regulation.

CONCLUSION

Our study comprehensively provides important information on miRNA functions and molecular mechanisms in response to N stress. These findings may assist to improve nitrogen availability in plants.