Genome-wide identification and characterization of COMT gene family during the development of blueberry fruit

Genome-Wide Identification of the Sulfate Transporters Gene Family in Blueberry (Vaccinium spp.) and Its Response to Ericoid Mycorrhizal Fungi

Sulfur metabolism plays a major role in plant growth and development, environmental adaptation, and material synthesis, and the sulfate transporters are the beginning of sulfur metabolism. We identified 37 potential VcSULTR genes in the blueberry genome, encoding peptides with 534 to 766 amino acids. The genes were grouped into four subfamilies in an evolutionary analysis. The 37 putative VcSULTR proteins ranged in size from 60.03 to 83.87 kDa. These proteins were predicted to be hydrophobic and mostly localize to the plasma membrane. The VcSULTR genes were distributed on 30 chromosomes; VcSULTR3;5b and VcSULTR3;5c were the only tandemly repeated genes. The VcSULTR promoters contained cis-acting elements related to the fungal symbiosis and stress responses. The transcript levels of the VcSULTRs differed among blueberry organs and changed in response to ericoid mycorrhizal fungi and sulfate treatments. A subcellular localization analysis showed that VcSULTR2;1c localized to, and functioned in, the plasma membrane and chloroplast. The virus-induced gene knock-down of VcSULTR2;1c resulted in a significantly decreased endogenous sulfate content, and an up-regulation of genes encoding key enzymes in sulfur metabolism (VcATPS2 and VcSiR1). These findings enhance our understanding of mycorrhizal-fungi-mediated sulfate transport in blueberry, and lay the foundation for further research on blueberry-mycorrhizal symbiosis.

Evolution and Analysis of Caffeic Acid Transferase (COMT) in Seed Plants

Caffeic acid transferase (COMT) is a key enzyme in the lignin and melatonin synthesis pathways and plays an important role in plant growth and development. All seed plants have two characteristics: they have vascular tissues, phloem, and xylem, and they can produce and reproduce seeds. In order to understand the distribution and evolution of COMTs in seed plants, we performed physicochemical property analysis, subcellular localization, phylogenetic analysis, conserved motif analysis, and protein interaction network analysis of 44 COMT homologs from 26 seed plants through in silico. The results showed that in seed plants, the structure of COMT genes tends to be stable in different plant taxa, while the relationship between the chromosomal positions of different COMT genes in the same plant was more intricate. The conserved distribution of COMT in seed plants reflected its highly specialized function.

Caffeic Acid O-Methyltransferase Gene Family in Mango (Mangifera indica L.) with Transcriptional Analysis under Biotic and Abiotic Stresses and the Role of MiCOMT1 in Salt Tolerance

Caffeic acid O-methyltransferase (COMT) participates in various physiological activities in plants, such as positive responses to abiotic stresses and the signal transduction of phytohormones. In this study, 18 COMT genes were identified in the chromosome-level reference genome of mango, named MiCOMTs. A phylogenetic tree containing nine groups (I-IX) was constructed based on the amino acid sequences of the 71 COMT proteins from seven species. The phylogenetic tree indicated that the members of the MiCOMTs could be divided into four groups. Quantitative real-time PCR showed that all MiCOMT genes have particularly high expression levels during flowering. The expression levels of MiCOMTs were different under abiotic and biotic stresses, including salt and stimulated drought stresses, ABA and SA treatment, as well as Xanthomonas campestris pv. mangiferaeindicae and Colletotrichum gloeosporioides infection, respectively. Among them, the expression level of MiCOMT1 was significantly up-regulated at 6-72 h after salt and stimulated drought stresses. The results of gene function analysis via the transient overexpression of the MiCOMT1 gene in Nicotiana benthamiana showed that the MiCOMT1 gene can promote the accumulation of ABA and MeJA, and improve the salt tolerance of mango. These results are beneficial to future researchers aiming to understand the biological functions and molecular mechanisms of MiCOMT genes.

Characterization of sugarcane mutants developed through gamma irradiations for their lignin content and caffeic acid-O-methyl transferase (COMT) gene mutations

PURPOSE

Bagasse, the residue left after extracting juice from sugarcane stalks, is rich in lignocellulosic biomass. The lignin present in this plant biomass is the key factor that hinders the efficient extraction of ethanol from the bagasse. In the current study, γ-irradiated sugarcane mutants were evaluated for variation in lignin content and its corresponding caffeic acid-O-methyl transferase (COMT) gene.

MATERIALS AND METHODS

The acetyl bromide method was used to estimate lignin content in sugarcane mutants. PCR-based cloning of the COMT gene was performed in low lignin mutants as well as control plants in E. coli (strain DH5α) to understand the mechanism of variation at the molecular level. The Sanger sequencing for cloned gene was performed to check variation in gene sequence.

RESULTS

In comparison to the control (21.5%), the mutant plants' lignin content ranged from 13 to 28%. The Sanger sequencing revealed approximately the same length of the gene from mutants as well as a control plant. In comparison to the reference gene, the mutated gene showed SNPs and indels in different regions, which may have an impact on lignin content.

CONCLUSIONS

Therefore, γ-irradiated mutagenesis is an acceptable approach to develop novel mutants of sugarcane with low lignin content to enhance bioethanol production from waste material using bioprocess technology.

Genome-wide identification of GH9 gene family and the assessment of its role during fruit abscission zone formation in Vaccinium ashei

KEY MESSAGE

The genomic location and stage-specific expression pattern of GH9 genes reveal their critical roles during fruit abscission zone formation in Vaccinium ashei. Glycosyl hydrolase family 9 (GH9) cellulases play a crucial role in both cellulose synthesis and hydrolysis during plant growth and development. Despite this importance, there is currently no study on the involvement of GH9-encoding genes, specifically VaGH9s, in abscission zone formation of rabbiteye blueberries (Vaccinium ashei). In this study, we identified a total of 61 VaGH9s in the genome, which can be classified into 3 subclasses based on conserved motifs and domains, gene structures, and phylogenetic analyses. Our synteny analysis revealed that VaGH9s are more closely related to the GH9s of Populus L. than to those of Arabidopsis, Vitis vinifera, and Citrus sinensis. In silico structural analysis predicted that most of VaGH9s are hydrophilic, and localized in cell membrane and/or cell wall, and the variable sets of cis-acting regulatory elements and functional diversity with four categories of stress response, hormone regulation, growth and development, and transcription factor-related elements are present in the promoter sequence of VaGH9s genes. Transcriptomic analysis showed that there were 22 differentially expressed VaGH9s in fruit abscission zone tissue at the veraison stage, and the expression of VaGH9B2 and VaGH9C10 was continuously increased during fruit maturation, which were in parallel with the increasing levels of cellulase activity and oxidative stress indicators, suggesting that they are involved in the separation stage of fruit abscission in Vaccinium ashei. Our work identified 22 VaGH9s potentially involved in different stages of fruit abscission and would aid further investigation into the molecular regulation of abscission in rabbiteye blueberries fruit.

RNA-Seq and Iso-Seq Reveal the Important Role of COMT and CCoAOMT Genes in Accumulation of Scopoletin in Noni (Morinda citrifolia)

Scopoletin, the main component of clinical drugs and the functional component of health products, is highly abundant in noni fruit (Morinda citrifolia). Multiple enzyme genes regulate scopoletin accumulation. In the present study, differentially expressed genes of noni were analyzed by RNA sequencing (RNA-Seq) and the full-length genes by isoform-sequencing (Iso-Seq) to find the critical genes in the scopoletin accumulation mechanism pathway. A total of 32,682 full-length nonchimeric reads (FLNC) were obtained, out of which 16,620 non-redundant transcripts were validated. Based on KEGG (Kyoto Encyclopedia of Genes and Genomes) annotation and differential expression analysis, two differentially expressed genes, caffeic acid 3-O-methyltransferase (COMT) and caffeoyl-CoA O-methyltransferase (CCoAOMT), were found in the scopoletin accumulation pathway of noni. Real-time quantitative polymerase chain reaction (q-PCR), phylogenetic tree analysis, gene expression analysis, and the change in scopoletin content confirmed that these two proteins are important in this pathway. Based on these results, the current study supposed that COMT and CCoAOMT play a significant role in the accumulation of scopoletin in noni fruit, and COMT (gene number: gene 7446, gene 8422, and gene 6794) and CCoAOMT (gene number: gene 12,084) were more significant. These results provide the importance of COMT and CCoAOMT and a basis for further understanding the accumulation mechanism of scopoletin in noni.

Identification and Functional Analysis of the Caffeic Acid O-Methyltransferase (COMT) Gene Family in Rice (Oryza sativa L.)

Caffeic acid O-methyltransferase (COMT) is one of the core enzymes involved in lignin synthesis. However, there is no systematic study on the rice COMT gene family. We identified 33 COMT genes containing the methyltransferase-2 domain in the rice genome using bioinformatic methods and divided them into Group I (a and b) and Group II. Motifs, conserved domains, gene structure and SNPs density are related to the classification of OsCOMTs. The tandem phenomenon plays a key role in the expansion of OsCOMTs. The expression levels of fourteen and thirteen OsCOMTs increased or decreased under salt stress and drought stress, respectively. OsCOMTs showed higher expression levels in the stem. The lignin content of rice was measured in five stages; combined with the expression analysis of OsCOMTs and multiple sequence alignment, we found that OsCOMT8, OsCOMT9 and OsCOMT15 play a key role in the synthesis of lignin. Targeted miRNAs and gene ontology annotation revealed that OsCOMTs were involved in abiotic stress responses. Our study contributes to the analysis of the biological function of OsCOMTs, which may provide information for future rice breeding and editing of the rice genome.

There and back again; historical perspective and future directions for Vaccinium breeding and research studies

The genus Vaccinium L. (Ericaceae) contains a wide diversity of culturally and economically important berry crop species. Consumer demand and scientific research in blueberry (Vaccinium spp.) and cranberry (Vaccinium macrocarpon) have increased worldwide over the crops' relatively short domestication history (~100 years). Other species, including bilberry (Vaccinium myrtillus), lingonberry (Vaccinium vitis-idaea), and ohelo berry (Vaccinium reticulatum) are largely still harvested from the wild but with crop improvement efforts underway. Here, we present a review article on these Vaccinium berry crops on topics that span taxonomy to genetics and genomics to breeding. We highlight the accomplishments made thus far for each of these crops, along their journey from the wild, and propose research areas and questions that will require investments by the community over the coming decades to guide future crop improvement efforts. New tools and resources are needed to underpin the development of superior cultivars that are not only more resilient to various environmental stresses and higher yielding, but also produce fruit that continue to meet a variety of consumer preferences, including fruit quality and health related traits.

Genome-Wide Comparative Analysis of the R2R3-MYB Gene Family in Five Solanaceae Species and Identification of Members Regulating Carotenoid Biosynthesis in Wolfberry

The R2R3-MYB is a large gene family involved in various plant functions, including carotenoid biosynthesis. However, this gene family lacks a comprehensive analysis in wolfberry (Lycium barbarum L.) and other Solanaceae species. The recent sequencing of the wolfberry genome provides an opportunity for investigating the organization and evolutionary characteristics of R2R3-MYB genes in wolfberry and other Solanaceae species. A total of 610 R2R3-MYB genes were identified in five Solanaceae species, including 137 in wolfberry. The LbaR2R3-MYB genes were grouped into 31 subgroups based on phylogenetic analysis, conserved gene structures, and motif composition. Five groups only of Solanaceae R2R3-MYB genes were functionally divergent during evolution. Dispersed and whole duplication events are critical for expanding the R2R3-MYB gene family. There were 287 orthologous gene pairs between wolfberry and the other four selected Solanaceae species. RNA-seq analysis identified the expression level of LbaR2R3-MYB differential gene expression (DEGs) and carotenoid biosynthesis genes (CBGs) in fruit development stages. The highly expressed LbaR2R3-MYB genes are co-expressed with CBGs during fruit development. A quantitative Real-Time (qRT)-PCR verified seven selected candidate genes. Thus, Lba11g0183 and Lba02g01219 are candidate genes regulating carotenoid biosynthesis in wolfberry. This study elucidates the evolution and function of R2R3-MYB genes in wolfberry and the four Solanaceae species.

Systematic analysis reveals O-methyltransferase gene family members involved in flavonoid biosynthesis in grape

O-methyltransferases (OMTs) are an important group of enzymes involved in the methylation of various secondary metabolites, including flavonoids. However, the features and functions of OMTs have not been comprehensively studied in grape (Vitis vinifera), a rich source of methylated flavonoids. Here, 47 OMT members were identified in grape genome. They were unevenly distributed on grape chromosomes and some genes were tandem duplicated, indicating the role of duplication processes in the expansion of this gene family. Based on the phylogenetic relationship, these OMTs were clustered into CCoAOMT and COMT subclades, which were further supported by the results of conserved motif and gene structure analysis. Correlation analysis revealed that three members (VvCCoAOMT1, VvCCoAOMT4, and VvCOMT1) were potentially involved in the synthesis of most methylated flavonoids in the berry skins. Expression profiling based on RNA-seq data and qRT-PCR experiments indicated that VvCCoAOMT1 and VvCCoAOMT4 had specific and high expression in berry skins, and responded to abscisic acid and high temperature treatments; and that VvCOMT1 expression was significantly induced during berry development and UVC treatment. Cis-regulatory element analysis suggested important roles of OMTs in growth, development, and defense against stresses. We further demonstrated the transcriptional regulation of VvCCoAOMT4 by VvMYBA1, a master regulator of grape berry anthocyanin, and verified the protein localization of VvCCoAOMT4 in membrane and nucleus. These findings facilitate a better understanding of the characteristics of OMT gene family, especially of the potential members involved in the formation of O-methylated flavonoids in grape.

Genome-Wide Identification and Characterization of Caffeic Acid O-Methyltransferase Gene Family in Soybean

Soybean is one of the most important legumes, providing high-quality protein for humans. The caffeic acid O-methyltransferase (COMT) gene has previously been demonstrated to be a critical gene that regulates lignin production in plant cell walls and plays an important function in plant growth and development. However, the COMT gene family has not been studied in soybeans. In this study, 55 COMT family genes in soybean were identified by phylogenetic analysis and divided into two groups, I and II. The analysis of conserved domains showed that all GmCOMTs genes contained Methyltransferase-2 domains. Further prediction of cis-acting elements showed that GmCOMTs genes were associated with growth, light, stress, and hormonal responses. Eventually, based on the genomic data of soybean under different stresses, the results showed that the expression of GmCOMTs genes was different under different stresses, such as salt and drought stress. This study has identified and characterized the COMT gene family in soybean, which provides an important theoretical basis for further research on the biological functions of COMT genes and promotes revealing the role of GmCOMTs genes under stress resistance.

Genome-Wide Identification and Characterization of DnaJ Gene Family in Grape (Vitis vinifera L.)

Grape production in southern China suffers great loss due to various environmental stresses. To understand the mechanism of how the grape plants respond to these stresses is an active area of research in developing cultivation techniques. Plant stress resistance is known to rely on special proteins. Amongst them, DnaJ protein (HSP40) serves as co-chaperones of HSP70, playing crucial roles in various stress response. However, the DnaJ proteins encoded by the DnaJ gene family in Vitis vinifera L. have not been fully described yet. In this study, we identified 78 VvDnaJs in the grape genome that can be classified into three groups—namely, DJA, DJB, and DJC. To reveal the evolutionary and stress response mechanisms for the VvDnaJ gene family, their evolutionary and expression patterns were analyzed using the bioinformatic approach and qRT-PCR. We found that the members in the same group exhibited a similar gene structure and protein domain organization. Gene duplication analysis demonstrated that segmental and tandem duplication may not be the dominant pathway of gene expansion in the VvDnaJ gene family. Codon usage pattern analysis showed that the codon usage pattern of VvDnaJs differs obviously from the monocotyledon counterparts. Tissue-specific analysis revealed that 12 VvDnaJs present a distinct expression profile, implying their distinct roles in various tissues. Cis-acting element analysis showed that almost all VvDnaJs contained the elements responsive to either hormones or stresses. Therefore, the expression levels of VvDnaJs subjected to exogenous hormone applications and stress treatments were determined, and we found that VvDnaJs were sensitive to hormone treatments and shade, salt, and heat stresses, especially VIT_00s0324g00040. The findings of this study could provide comprehensive information for the further investigation on the genetics and protein functions of the DnaJ gene family in grape.

Genome-Wide Identification and Characterization of GhCOMT Gene Family during Fiber Development and Verticillium Wilt Resistance in Cotton

Caffeic acid O-methyltransferases (COMTs) play an essential role in lignin synthesis procession, especially in the plant’s phenylalanine metabolic pathway. The content of COMT genes in cotton and the relationship between their expression patterns have not been studied clearly in cotton. In this study, we have identified 190 COMT genes in cotton, which were classified into three groups (I, II and III), and mapped on the cotton chromosomes. In addition, we found that 135 of the 190 COMT genes result from dispersed duplication (DSD) and whole-genome duplication (WGD), indicating that DSD and WGD were the main forces driving COMT gene expansion. The Ka/Ks analysis showed that GhCOMT43 and GhCOMT41 evolved from GaCOMT27 and GrCOMT14 through positive selection. The results of qRT-PCR showed that GhCOMT13, GhCOMT28, GhCOMT39 and GhCOMT55 were related to lignin content during the cotton fiber development. GhCOMT28, GhCOMT39, GhCOMT55, GhCOMT56 and GhCOMT57 responded to Verticillium Wilt (VW) and maybe related to VW resistance through lignin synthesis. Conclusively, this study found that GhCOMTs were highly expressed in the secondary wall thickening stage and VW. These results provide a clue for studying the functions of GhCOMTs in the development of cotton fiber and VW resistance and could lay a foundation for breeding cotton cultivates with higher quantity and high resistance to VW.

Genome-Wide Identification and Analysis of the WRKY Gene Family and Cold Stress Response in Acer truncatum

WRKY transcription factors constitute one of the largest gene families in plants and are involved in many biological processes, including growth and development, physiological metabolism, and the stress response. In earlier studies, the WRKY gene family of proteins has been extensively studied and analyzed in many plant species. However, information on WRKY transcription factors in Acer truncatum has not been reported. In this study, we conducted genome-wide identification and analysis of the WRKY gene family in A. truncatum, 54 WRKY genes were unevenly located on all 13 chromosomes of A. truncatum, the highest number was found in chromosomes 5. Phylogenetic relationships, gene structure, and conserved motif identification were constructed, and the results affirmed 54 AtruWRKY genes were divided into nine subgroup groups. Tissue species analysis of AtruWRKY genes revealed which were differently exhibited upregulation in flower, leaf, root, seed and stem, and the upregulation number were 23, 14, 34, 18, and 8, respectively. In addition, the WRKY genes expression in leaf under cold stress showed that more genes were significantly expressed under 0, 6 and 12 h cold stress. The results of this study provide a new insight the regulatory function of WRKY genes under abiotic and biotic stresses.

Transcriptome analysis and identification of genes associated with floral transition and fruit development in rabbiteye blueberry (Vaccinium ashei)

Flowering and fruit set are important traits affecting fruit quality and yield in rabbiteye blueberry (Vaccinium ashei). Intense efforts have been made to elucidate the influence of vernalization and phytohormones on flowering, but the molecular mechanisms of flowering and fruit set remain unclear. To unravel these mechanisms, we performed transcriptome analysis to explore blueberry transcripts from flowering to early fruit stage. We divided flowering and fruit set into flower bud (S2), initial flower (S3), bloom flower (S4), pad fruit (S5), and cup fruit (S6) based on phenotype and identified 1,344, 69, 658, and 189 unique differentially expressed genes (DEGs) in comparisons of S3/S2, S4/S3, S5/S4, and S6/S5, respectively. There were obviously more DEGs in S3/S2 and S5/S4 than in S4/S3, and S6/S5, suggesting that S3/S2 and S5/S4 represent major transitions from buds to fruit in blueberry. GO and KEGG enrichment analysis indicated these DEGs were mostly enriched in phytohormone biosynthesis and signaling, transporter proteins, photosynthesis, anthocyanins biosynthesis, disease resistance protein and transcription factor categories, in addition, transcript levels of phytohormones and transporters changed greatly throughout the flowering and fruit set process. Gibberellic acid and jasmonic acid mainly acted on the early stage of flowering development like expression of the florigen gene FT, while the expression of auxin response factor genes increased almost throughout the process from bud to fruit development. Transporter proteins were mainly associated with minerals during the early flowering development stage and sugars during the early fruit stage. At the early fruit stage, anthocyanins started to accumulate, and the fruit was susceptible to diseases such as fungal infection. Expression of the transcription factor MYB86 was up-regulated during initial fruit development, which may promote anthocyanin accumulation. These results will aid future studies exploring the molecular mechanism underlying flowering and fruit set of rabbiteye blueberry.