Multi-responses of O-methyltransferase genes to salt stress and fiber development of Gossypium species

Diversity and regioselectivity of O-methyltransferases catalyzing the formation of O-methylated flavonoids

Flavonoids and their methylated derivatives have immense market potential in the food and biomedical industries due to their multiple beneficial effects, such as antimicrobial, anti-inflammatory, and anticancer activities. The biological synthesis of flavonoids and their derivatives is often accomplished via the use of genetically modified microorganisms to ensure large-scale production. Therefore, it is pivotal to understand the properties of O-methyltransferases (OMTs) that mediate the methylation of flavonoids. However, the properties of these OMTs are governed by their: sources, substrate specificity, amino acid residues in the active sites, and the intricate mechanism. In order to obtain a clue for the selection of suitable OMTs for the biosynthesis of a target methylated flavonoid, we made a comprehensive review of the currently reported results, with a particular focus on their comparative regioselectivity for different flavonoid substrates. Additionally, the possible mechanisms for the diversity of this class of enzymes were explored using molecular simulation technology. Finally, major gaps in our understanding and areas for future studies were discussed. The findings of this study may be useful in selecting genes that encode OMTs and designing enzyme-based processes for synthesizing O-methylated flavonoids.

Genetic diversity and candidate genes for transient waterlogging tolerance in mungbean at the germination and seedling stages

Mungbean [Vigna radiata var. radiata (L.) Wilczek] production in Asia is detrimentally affected by transient soil waterlogging caused by unseasonal and increasingly frequent extreme precipitation events. While mungbean exhibits sensitivity to waterlogging, there has been insufficient exploration of germplasm for waterlogging tolerance, as well as limited investigation into the genetic basis for tolerance to identify valuable loci. This research investigated the diversity of transient waterlogging tolerance in a mini-core germplasm collection of mungbean and identified candidate genes for adaptive traits of interest using genome-wide association studies (GWAS) at two critical stages of growth: germination and seedling stage (i.e., once the first trifoliate leaf had fully-expanded). In a temperature-controlled glasshouse, 292 genotypes were screened for tolerance after (i) 4 days of waterlogging followed by 7 days of recovery at the germination stage and (ii) 8 days of waterlogging followed by 7 days of recovery at the seedling stage. Tolerance was measured against drained controls. GWAS was conducted using 3,522 high-quality DArTseq-derived SNPs, revealing five significant associations with five phenotypic traits indicating improved tolerance. Waterlogging tolerance was positively correlated with the formation of adventitious roots and higher dry masses. FGGY carbohydrate kinase domain-containing protein was identified as a candidate gene for adventitious rooting and mRNA-uncharacterized LOC111241851, Caffeoyl-CoA O-methyltransferase At4g26220 and MORC family CW-type zinc finger protein 3 and zinc finger protein 2B genes for shoot, root, and total dry matter production. Moderate to high broad-sense heritability was exhibited for all phenotypic traits, including seed emergence (81%), adventitious rooting (56%), shoot dry mass (81%), root dry mass (79%) and SPAD chlorophyll content (70%). The heritability estimates, marker-trait associations, and identification of sources of waterlogging tolerant germplasm from this study demonstrate high potential for marker-assisted selection of tolerance traits to accelerate breeding of climate-resilient mungbean varieties.

Hub Genes in Stable QTLs Orchestrate the Accumulation of Cottonseed Oil in Upland Cotton via Catalyzing Key Steps of Lipid-Related Pathways

Upland cotton is the fifth-largest oil crop in the world, with an average supply of nearly 20% of vegetable oil production. Cottonseed oil is also an ideal alternative raw material to be efficiently converted into biodiesel. However, the improvement in kernel oil content (KOC) of cottonseed has not received sufficient attention from researchers for a long time, due to the fact that the main product of cotton planting is fiber. Previous studies have tagged QTLs and identified individual candidate genes that regulate KOC of cottonseed. The regulatory mechanism of oil metabolism and accumulation of cottonseed are still elusive. In the current study, two high-density genetic maps (HDGMs), which were constructed based on a recombinant inbred line (RIL) population consisting of 231 individuals, were used to identify KOC QTLs. A total of forty-three stable QTLs were detected via these two HDGM strategies. Bioinformatic analysis of all the genes harbored in the marker intervals of the stable QTLs revealed that a total of fifty-one genes were involved in the pathways related to lipid biosynthesis. Functional analysis via coexpression network and RNA-seq revealed that the hub genes in the co-expression network that also catalyze the key steps of fatty acid synthesis, lipid metabolism and oil body formation pathways (ACX4, LACS4, KCR1, and SQD1) could jointly orchestrate oil accumulation in cottonseed. This study will strengthen our understanding of oil metabolism and accumulation in cottonseed and contribute to KOC improvement in cottonseed in the future, enhancing the security and stability of worldwide food supply.

Evolution and diversification of the O-methyltransferase (OMT) gene family in Solanaceae

O-methyltransferases (OMTs) are a group of enzymes involved in several fundamental biological processes in plants, including lignin biosynthesis, pigmentation, and aroma production. Despite the intensive investigation of the role of OMTs in plant secondary metabolism, the evolution and diversification of this gene family in Solanaceae remain poorly understood. Here, we conducted a genome-wide survey of OMT genes in six Solanaceae species, reconstructing gene phylogenetic trees, predicting the potential involvement in biological processes, and investigating the exon/intron structure and chromosomal location. We identified 57 caffeoyl-CoA OMTs (CCoAOMTs) and 196 caffeic acid OMTs (COMTs) in the studied species. We observed a conserved gene block on chromosome 2 that consisted of tandem duplicated copies of OMT genes. Our results suggest that the expansion of the OMT gene family in Solanaceae was driven by whole genome duplication, segmental duplication, and tandem duplication, with multiple genes being retained by neofunctionalization and subfunctionalization. This study represents an essential first step in unraveling the evolutionary history of OMTs in Solanaceae. Our findings deepen our understanding of the crucial role of OMTs in several biological processes and highlight their significance as potential biotechnological targets.

QTL Verification and Candidate Gene Screening of Fiber Quality and Lint Percentage in the Secondary Segregating Population of Gossypium hirsutum

Fiber quality traits, especially fiber strength, length, and micronaire (FS, FL, and FM), have been recognized as critical fiber attributes in the textile industry, while the lint percentage (LP) was an important indicator to evaluate the cotton lint yield. So far, the genetic mechanism behind the formation of these traits is still unclear. Quantitative trait loci (QTL) identification and candidate gene validation provide an effective methodology to uncover the genetic and molecular basis of FL, FS, FM, and LP. A previous study identified three important QTL/QTL cluster loci, harboring at least one of the above traits on chromosomes A01, A07, and D12 via a recombinant inbred line (RIL) population derived from a cross of Lumianyan28 (L28) × Xinluzao24 (X24). A secondary segregating population (F2) was developed from a cross between L28 and an RIL, RIL40 (L28 × RIL40). Based on the population, genetic linkage maps of the previous QTL cluster intervals on A01 (6.70-10.15 Mb), A07 (85.48-93.43 Mb), and D12 (0.40-1.43 Mb) were constructed, which span 12.25, 15.90, and 5.56 cM, with 2, 14, and 4 simple sequence repeat (SSR) and insertion/deletion (Indel) markers, respectively. QTLs of FL, FS, FM, and LP on these three intervals were verified by composite interval mapping (CIM) using WinQTL Cartographer 2.5 software via phenotyping of F2 and its derived F2:3 populations. The results validated the previous primary QTL identification of FL, FS, FM, and LP. Analysis of the RNA-seq data of the developing fibers of L28 and RIL40 at 10, 20, and 30 days post anthesis (DPA) identified seven differentially expressed genes (DEGs) as potential candidate genes. qRT-PCR verified that five of them were consistent with the RNA-seq result. These genes may be involved in regulating fiber development, leading to the formation of FL, FS, FM, and LP. This study provides an experimental foundation for further exploration of these functional genes to dissect the genetic mechanism of cotton fiber development.

Identification and Characterization of Phytocyanin Family Genes in Cotton Genomes

Phytocyanins (PCs) are a class of plant-specific blue copper proteins that have been demonstrated to play a role in electron transport and plant development. Through analysis of the copper ligand residues, spectroscopic properties, and domain architecture of the protein, PCs have been grouped into four subfamilies: uclacyanins (UCs), stellacyanins (SCs), plantacyanins (PLCs), and early nodulin-like proteins (ENODLs). The present study aimed to identify and characterise the PCs present in three distinct cotton species (Gossypium hirsutum, Gossyium arboreum, and Gossypium raimondii) through the identification of 98, 63, and 69 genes respectively. We grouped PCs into four clades by using bioinformatics analysis and sequence alignment, which exhibit variations in gene structure and motif distribution. PCs are distributed across all chromosomes in each of the three species, with varying numbers of exons per gene and multiple conserved motifs, and with a minimum of 1 and maximum of 11 exons found on one gene. Transcriptomic data and qRT-PCR analysis revealed that two highly differentiated PC genes were expressed at the fibre initiation stage, while three highly differentiated PCs were expressed at the fibre elongation stage. These findings serve as a foundation for further investigations aimed at understanding the contribution of this gene family in cotton fibre production.

Lignin and Its Pathway-Associated Phytoalexins Modulate Plant Defense against Fungi

Fungi infections cause approximately 60-70% yield loss through diseases such as rice blast, powdery mildew, Fusarium rot, downy mildew, etc. Plants naturally respond to these infections by eliciting an array of protective metabolites to confer physical or chemical protection. Among plant metabolites, lignin, a phenolic compound, thickens the middle lamella and the secondary cell walls of plants to curtail fungi infection. The biosynthesis of monolignols (lignin monomers) is regulated by genes whose transcript abundance significantly improves plant defense against fungi. The catalytic activities of lignin biosynthetic enzymes also contribute to the accumulation of other defense compounds. Recent advances focus on modifying the lignin pathway to enhance plant growth and defense against pathogens. This review presents an overview of monolignol regulatory genes and their contributions to fungi immunity, as reported over the last five years. This review expands the frontiers in lignin pathway engineering to enhance plant defense.

Comparative DNA Methylome of Phytoplasma Associated Retrograde Metamorphosis in Sesame (Sesamum indicum L.)

Phytoplasma-associated diseases such as phyllody and little leaf are critical threats to sesame cultivation worldwide. The mechanism of the dramatic conversion of flowers to leafy structures leading to yield losses and the drastic reduction in leaf size due to Phytoplasma infection remains yet to be identified. Cytosine methylation profiles of healthy and infected sesame plants studied using Whole Genome Bisulfite Sequencing (WGBS) and Quantitative analysis of DNA methylation with the real-time PCR (qAMP) technique revealed altered DNA methylation patterns upon infection. Phyllody was associated with global cytosine hypomethylation, though predominantly in the CHH (where H = A, T or C) context. Interestingly, comparable cytosine methylation levels were observed between healthy and little leaf-affected plant samples in CG, CHG and CHH contexts. Among the different genomic fractions, the highest number of differentially methylated Cytosines was found in the intergenic regions, followed by promoter, exonic and intronic regions in decreasing order. Further, most of the differentially methylated genes were hypomethylated and were mainly associated with development and defense-related processes. Loci for STOREKEEPER protein-like, a DNA-binding protein and PP2-B15, an F-Box protein, responsible for plugging sieve plates to maintain turgor pressure within the sieve tubes were found to be hypomethylated by WGBS, which was confirmed by methylation-dependent restriction digestion and qPCR. Likewise, serine/threonine-protein phosphatase-7 homolog, a positive regulator of cryptochrome signaling involved in hypocotyl and cotyledon growth and probable O-methyltransferase 3 locus were determined to be hypermethylated. Phytoplasma infection-associated global differential methylation as well as the defense and development-related loci reported here for the first time significantly elucidate the mechanism of phytoplasma-associated disease development.

Genome-Wide Identification and Expression Analysis Elucidates the Potential Role of PFK Gene Family in Drought Stress Tolerance and Sugar Metabolism in Cotton

Drought has been identified as a major threat for global crop production worldwide. Phosphofructokinase (PFK) is vital for sugar metabolism. During phosphorylation, plants have two enzymes: ATP-dependent phosphofructokinase (PFK) and pyrophosphate-dependent fructose-6-phosphate phosphotransferase (PFP). Genome-wide identification led to the identification of 80 PFK genes, 26 genes in G. hirsutum and G. barbadense, and 14 genes in G. arboreum and G. raimondii. Phylogenetic, gene structure, and motif analyses showed that PFK genes were grouped into two main categories, namely, PFK and PFP, with 18 and 8 genes in the allotetraploid species and 10 PFK and 4 PFP genes in the diploid species, respectively. Using the RNA-seq expressions of 26 genes from GhPFK, a co-expression network analysis was performed to identify the hub genes. GhPFK04, GhPFK05, GhPFK09, GhPFK11, GhPFK13, GhPFK14, and GhPFK17 in leaves and GhPFK02, GhPFK09, GhPFK11, GhPFK15, GhPFK16, and GhPFK17 in root tissues were found as hub genes. RT-qPCR analysis validated the expressions of identified hub genes. Interestingly, GhPFK11 and GhPFK17 were identified as common hub genes, and these might be the true candidate genes involved in the drought stress tolerance. In the KEGG enrichment analysis, amino acids such as L-valine, L-histidine, L-glutamine, L-serine, L-homoserine, L-methionine, L-cysteine, and gluconic acid were significantly upregulated, whereas sugars, mainly fructose-1-phosphate, D-mannitol, D-sorbitol, dulcitol, and lactose, were significantly downregulated during drought stress. Genome-wide analysis paves the way for a deeper understanding of the PFK genes and establishes the groundwork for future research into PFK’s role in enhancing drought stress tolerance and sugar metabolism in cotton.

The Pivotal Role of Major Chromosomes of Sub-Genomes A and D in Fiber Quality Traits of Cotton

Lack of precise information about the candidate genes involved in a complex quantitative trait is a major obstacle in the cotton fiber quality improvement, and thus, overall genetic gain in conventional phenotypic selection is low. Recent molecular interventions and advancements in genome sequencing have led to the development of high-throughput molecular markers, quantitative trait locus (QTL) fine mapping, and single nucleotide polymorphisms (SNPs). These advanced tools have resolved the existing bottlenecks in trait-specific breeding. This review demonstrates the significance of chromosomes 3, 7, 9, 11, and 12 of sub-genomes A and D carrying candidate genes for fiber quality. However, chromosome 7 carrying SNPs for stable and potent QTLs related to fiber quality provides great insights for fiber quality-targeted research. This information can be validated by marker-assisted selection (MAS) and transgene in Arabidopsis and subsequently in cotton.

Structure basis of the caffeic acid O-methyltransferase from Ligusiticum chuanxiong to understand its selective mechanism

Caffeic acid O-methyltransferase from Ligusticum chuanxiong (LcCOMT) showed strict regiospecificity despite a relative degree of preference. Compared with caffeic acid, methyl caffeate was the preferential substrate by its low Km and high Kcat. In this study, we obtained the SAM binary (1.80 Å) and SAH binary (1.95 Å) complex LcCOMT crystal structures, and established the ternary complex structure with methyl caffeate by molecular docking. The active site of LcCOMT included phenolic substrate pocket, SAM/SAH ligand pocket and conserved catalytic residues as well. The regiospecificity of LcCOMT that permitted only 3-hydroxyl group to be methylated arise from the interactions between the active site and the phenyl ring. However, the propanoid tail governed the relative preference of LcCOMT. The ester group in methyl caffeate stabilized the anionic intermediate caused by His268-Asp269 pair, whereas caffeic acid was unable to stabilize the anionic intermediate due to the adjacent carboxylate anion in the propanoid tail. Ser183 residue formed an additional hydrogen bond with SAH and its role was identified by S183A mutation. Ile318 residue might be a potential site for determination of substrate preference, and its mutation led to the change of tertiary conformation. The results supported the selective mechanism of LcCOMT.

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.

Study of Triticum aestivum Resistome in Response to Wheat dwarf India Virus Infection

Susceptible and resistant germplasm respond differently to pathogenic attack, including virus infections. We compared the transcriptome changes between a resistant wheat cultivar, Sonalika, and a susceptible cultivar, WL711, to understand this process in wheat against wheat dwarf India virus (WDIV) infection. A total of 2760 and 1853 genes were differentially expressed in virus-infected and mock-inoculated Sonalika, respectively, compared to WL711. The overrepresentation of genes involved in signaling, hormone metabolism, enzymes, secondary metabolites, proteolysis, and transcription factors was documented, including the overexpression of multiple PR proteins. We hypothesize that the virus resistance in Sonalika is likely due to strong intracellular surveillance via the action of multiple PR proteins (PR1, RAR1, and RPM1) and ChiB. Other genes such as PIP1, LIP1, DnaJ, defensins, oxalate oxidase, ankyrin repeat protein, serine-threonine kinase, SR proteins, beta-1,3-glucanases, and O-methyltransferases had a significant differential expression and play roles in stress tolerance, may also be contributing towards the virus resistance in Sonalika. In addition, we identified putative genes with unknown functions, which are only expressed in response to WDIV infection in Sonalika. The role of these genes could be further validated and utilized in engineering resistance in wheat and other crops.