In-silico analysis of cis-acting regulatory elements of pathogenesis-related proteins of Arabidopsis thaliana and Oryza sativa

Global identification of LIM genes in response to different heat stress regimes in Lactuca sativa

BACKGROUND

LIM (Lineage-11 (LIN-11), Insulin-1 (ISL-1), and Mechanotransduction-3 (MEC-3)) genes belong to a family that hold ubiquitous properties contributing to organ, seed, and pollen development as well as developmental and cellular responses to biotic and abiotic stresses. Lettuce (Lactuca sativa) is a highly consumed vegetable crop susceptible heat stress. High temperatures limit lettuce's overall yield, quality and marketability. Lettuce LIM genes have not been identified and their role in response to high temperatures is not known. Aiming to identify potential new targets for thermoresilience, we searched for LIM genes in lettuce and compared them with orthologous of several dicotyledons and monocotyledons plant species.

RESULTS

We identified fourteen lettuce LIM genes distributed into eight different subgroups using a genome-wide analysis strategy. Three belonging to DAR (DA means "large" in Chinese) class I, two DAR class II, one in the WLIM1, two in the WLIM2, one in the PLIM1, two in the PLIM2 class, one ßLIM and two δLIMs. No DAR-like were identified in any of the species analyzed including lettuce. Interestingly, unlike other gene families in lettuce which underwent large genome tandem duplications, LIM genes did not increase in number compared to other plant species. The response to heat stress induced a dynamic transcriptional response on LsLIM genes. All heat stress regimes, including night stress, day stress and day and night stress were largely responsible for changes in LIM transcriptional expression.

CONCLUSIONS

Our global analysis at the genome level provides a detailed identification of LIM genes in lettuce and other dicotyledonous and monocotyledonous plant species. Gene structure, physical and chemical properties as well as chromosomal location and Cis-regulatory element analysis together with our gene expression analysis under different temperature regimes identified LsWLIM1, LsWLIM2b, LsDAR3 and LsDAR5 as candidate genes that could be used by breeding programs aiming to produce lettuce varieties able to withstand high temperatures.

Genome-wide identification, and gene expression analysis of NBS-LRR domain containing R genes in Chenopodium quinoa for unveiling the dynamic contribution in plant immunity against Cercospora cf. chenopodii

The plant R genes encode the NLR proteins comprising nucleotide-binding sites (NBS) and variable-length C-terminal leucine-rich repeat domains. The proteins act as intracellular immune receptors and recognize effector proteins of phytopathogens, which convene virulence. Among stresses, diseases contribute majorly to yield loss in crop plants, and R genes confer disease resistance against phytopathogens. We investigated the NLRome of Chenopodium quinoa for intraspecific diversity, characterization, and contribution to immune response regulation against phytopathogens. One eighty-three NBS proteins were identified and grouped into four distinct classes. Exon-intron organization displayed discrimination in gene structure patterns among NLR proteins. Thirty-eight NBS proteins revealed ontology with defense response, ADP binding, and inter alia cellular components. These proteins had shown functional homology with disease-resistance proteins involved in the plant-pathogen interaction pathway. Likewise, expression analysis demonstrated that NLRs encoding genes showed differential expression patterns. However, most genes displayed high expression levels in plant defense response with varying magnitude compared to ADP binding and cellular components. Twenty-four NBS genes were selected based on Heatmap analysis for quantitative polymerase chain reaction under Cercospora disease stress, and their progressive expression pattern provides insights into their functional role under stress conditions. The protein-protein interaction analysis revealed functional enrichment of NLR proteins in regulating hypersensitive, immune, and stress responses. This study, the first to identify and characterize NBS genes in C. quinoa, reveals their contribution to disease response and divulges their dynamic involvement in inducing plant immunity against phytopathogens.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01475-0.

Genome-Wide Analysis of Serine Carboxypeptidase-like Genes in Soybean and Their Roles in Stress Resistance

The serine carboxypeptidase-like (SCPL) gene family plays a crucial role in the regulation of plant growth, development, and stress response through activities such as acyltransferases in plant secondary metabolism pathways. Although SCPL genes have been identified in various plant species, their specific functions and characteristics in soybean (Glycine max) have not yet been studied. We identified and characterized 73 SCPL genes, grouped into three subgroups based on gene structure and phylogenetic relationships. These genes are distributed unevenly across 20 soybean chromosomes and show varied codon usage patterns influenced by both mutation and selection pressures. Gene ontology (GO) enrichment suggests these genes are involved in plant cell wall regulation and stress responses. Expression analysis in various tissues and under stress conditions, including the presence of numerous stress-related cis-acting elements, indicated that these genes have varied expression patterns. This suggests that they play specialized roles such as modulating plant defense mechanisms against nematode infections, enhancing tolerance to drought and high salinity, and responding to cold stress, thereby helping soybean adapt to environmental stresses. Moreover, the expression of specific GmSCPLs was significantly affected following exposure to nematode infection, drought, high salt (NaCl), and cold stresses. Our findings underscore the potential of SCPL genes in enhancing stress resistance in soybean, providing a valuable resource for future genetic improvement and breeding strategies.

Involvement of the GH38 Family Exoglycosidase α-Mannosidase in Strawberry Fruit Ripening

Exoglycosidase enzymes hydrolyze the N-glycosylations of cell wall enzymes, releasing N-glycans that act as signal molecules and promote fruit ripening. Vesicular exoglycosidase α-mannosidase enzymes of the GH38 family (EC 3.2.1.24; α-man) hydrolyze N-glycans in non-reduced termini. Strawberry fruit (Fragaria × ananassa) is characterized by rapid softening as a result of cell wall modifications during the fruit ripening process. Enzymes acting on cell wall polysaccharides explain the changes in fruit firmness, but α-man has not yet been described in F. × ananassa, meaning that the indirect effects of N-glycan removal on its fruit ripening process are unknown. The present study identified 10 GH38 α-man sequences in the F. × ananassa genome with characteristic conserved domains and key residues. A phylogenetic tree built with the neighbor-joining method and three groups of α-man established, of which group I was classified into three subgroups and group III contained only Poaceae spp. sequences. The real-time qPCR results demonstrated that FaMAN genes decreased during fruit ripening, a trend mirrored by the total enzyme activity from the white to ripe stages. The analysis of the promoter regions of these FaMAN genes was enriched with ripening and phytohormone response elements, and contained cis-regulatory elements related to stress responses to low temperature, drought, defense, and salt stress. This study discusses the relevance of α-man in fruit ripening and how it can be a useful target to prolong fruit shelf life.

Genome-wide identification and characterization of DTX family genes highlighting their locations, functions, and regulatory factors in banana (Musa acuminata)

The detoxification efflux carriers (DTX) are a significant group of multidrug efflux transporter family members that play diverse functions in all kingdoms of living organisms. However, genome-wide identification and characterization of DTX family transporters have not yet been performed in banana, despite its importance as an economic fruit plant. Therefore, a detailed genome-wide analysis of DTX family transporters in banana (Musa acuminata) was conducted using integrated bioinformatics and systems biology approaches. In this study, a total of 37 DTX transporters were identified in the banana genome and divided into four groups (I, II, III, and IV) based on phylogenetic analysis. The gene structures, as well as their proteins' domains and motifs, were found to be significantly conserved. Gene ontology (GO) annotation revealed that the predicted DTX genes might play a vital role in protecting cells and membrane-bound organelles through detoxification mechanisms and the removal of drug molecules from banana cells. Gene regulatory analyses identified key transcription factors (TFs), cis-acting elements, and post-transcriptional regulators (miRNAs) of DTX genes, suggesting their potential roles in banana. Furthermore, the changes in gene expression levels due to pathogenic infections and non-living factor indicate that banana DTX genes play a role in responses to both biotic and abiotic stresses. The results of this study could serve as valuable tools to improve banana quality by protecting them from a range of environmental stresses.

How to awaken a sleeping giant: antagonistic expression of Flowering locus T homologs and elements of the age-related pathway are associated with the flowering transition in Agave tequilana

KEY MESSAGE

Antagonistic expression of Flowering locus T proteins and the ageing pathway via miRNAs and sugar metabolism regulate the initiation of flowering in A. tequilana. Flowering in commercial plantations of Agave tequilana signals that plants are ready to harvest for tequila production. However, time of flowering is often unpredictable and a detailed understanding of the process would be beneficial in the field, for breeding and for the development of future research. This report describes the functional analysis of A. tequilana FLOWERING LOCUS T (FT) genes by heterologous expression in A. thaliana and in situ hybridization in agave plants. The gene structures of the Agave tequilana FT family are also described and putative regulatory promoter elements were identified. Most Agave species have monocarpic, perennial life cycles that can last over 25 years during which plants do not respond to the normal environmental signals which induce flowering, suggesting that the ageing pathway as described in Arabidopsis may play an important role in determining flowering time in these species. Elements of this pathway were analyzed and in silico data is presented that supports the regulation of SQUAMOSA PROMOTER BINDING LIKE proteins (SPL), APETALA2 (AP2) proteins and members of Plant Glycoside Hydrolase Family 32 (PGHF32) by interactions with miRNAs 156, 172 and 164 during the initiation of flowering in A. tequilana.

The impact of multiple abiotic stresses on ns-LTP2.8 gene transcript and ns-LTP2.8 protein accumulation in germinating barley (Hordeum vulgare L.) embryos

Abiotic stresses occur more often in combination than alone under regular field conditions limiting in more severe way crop production. Stress recognition in plants primarily occurs in the plasma membrane, modification of which is necessary to maintain homeostasis in response to it. It is known that lipid transport proteins (ns-LTPs) participate in modification of the lipidome of cell membranes. Representative of this group, ns-LTP2.8, may be involved in the reaction to abiotic stress of germinating barley plants by mediating the intracellular transport of hydrophobic particles, such as lipids, helping to maintain homeostasis. The ns-LTP2.8 protein was selected for analysis due to its ability to transport not only linear hydrophobic molecules but also compounds with a more complex spatial structure. Moreover, ns-LTP2.8 has been qualified as a member of pathogenesis-related proteins, which makes it particularly important in relation to its high allergenic potential. This paper demonstrates for the first time the influence of various abiotic stresses acting separately as well as in their combinations on the change in the ns-LTP2.8 transcript, ns-LTP2.8 protein and total soluble protein content in the embryonal axes of germinating spring barley genotypes with different ns-LTP2.8 allelic forms and stress tolerance. Tissue localization of ns-LTP2.8 transcript as well as ns-LTP2.8 protein were also examined. Although the impact of abiotic stresses on the regulation of gene transcription and translation processes remains not fully recognized, in this work we managed to demonstrate different impact on applied stresses on the fundamental cellular processes in very little studied tissue of the embryonal axis of barley.

Genome-wide identification and in silico characterization of major RNAi gene families in date palm (Phoenix dactylifera)

BACKGROUND

Dates contain various minerals that are essential for good health. The major RNA interference (RNAi) gene families play a vital role in plant growth and development by controlling the expression of protein-coding genes against different biotic and abiotic stresses. However, these gene families for date palm are not yet studied. Therefore, this study has explored major RNAi genes and their characteristics in date palm.

RESULTS

We have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome by using AtRNAi genes as query sequences in BLASTp search. Domain analysis of predicted RNAi genes has revealed the Helicase_C, Dicer_dimer, PAZ, RNase III, and Piwi domains that are associated with the gene silencing mechanisms. Most PdRNAi proteins have been found in the nucleus and cytosol associated with the gene silencing actions. The gene ontology (GO) enrichment analysis has revealed some important GO terms including RNA interference, dsRNA fragmentation, and ribonuclease_III activity that are related to the protein-coding gene silencing mechanisms. Gene regulatory network (GRN) analysis has identified PAZ and SNF2 as the transcriptional regulators of PdRNAi genes. Top-ranked 10 microRNAs including Pda-miR156b, Pda-miR396a, Pda-miR166a, Pda-miR167d, and Pda-miR529a have been identified as the key post-transcriptional regulators of PdRNAi genes that are associated with different biotic/abiotic stresses. The cis-acting regulatory element analysis of PdRNAi genes has detected some vital cis-acting elements including ABRE, MBS, MYB, MYC, Box-4, G-box, I-box, and STRE that are linked with different abiotic stresses.

CONCLUSION

The results of this study might be valuable resources for the improvement of different characteristics in date palm by further studies in wet-lab.

Genome-Wide Comparative Analysis of Five Amaranthaceae Species Reveals a Large Amount of Repeat Content

Amaranthus is a genus of C4 dicotyledonous herbaceous plant species that are widely distributed in Asia, Africa, Australia, and Europe and are used as grain, vegetables, forages, and ornamental plants. Amaranth species have gained significant attention nowadays as potential sources of nutritious food and industrial products. In this study, we performed a comparative genome analysis of five amaranth species, namely, Amaranthus hypochondriacus, Amaranthus tuberculatus, Amaranthus hybridus, Amaranthus palmeri, and Amaranthus cruentus. The estimated repeat content ranged from 54.49% to 63.26% and was not correlated with the genome sizes. Out of the predicted repeat classes, the majority of repetitive sequences were Long Terminal Repeat (LTR) elements, which account for about 13.91% to 24.89% of all amaranth genomes. Phylogenetic analysis based on 406 single-copy orthologous genes revealed that A. hypochondriacus is most closely linked to A. hybridus and distantly related to A. cruentus. However, dioecious amaranth species, such as A. tuberculatus and A. palmeri, which belong to the subgenera Amaranthus Acnida, have formed their distinct clade. The comparative analysis of genomic data of amaranth species will be useful to identify and characterize agronomically important genes and their mechanisms of action. This will facilitate genomics-based, evolutionary studies, and breeding strategies to design faster, more precise, and predictable crop improvement programs.

Genome-Wide Comprehensive Identification and In Silico Characterization of Lectin Receptor-Like Kinase Gene Family in Barley (Hordeum vulgare L.)

Lectin receptor-like kinases (LecRLKs) are a significant subgroup of the receptor-like kinases (RLKs) protein family. They play crucial roles in plant growth, development, immune responses, signal transduction, and stress tolerance. However, the genome-wide identification and characterization of LecRLK genes and their regulatory elements have not been explored in a major cereal crop, barley (Hordeum vulgare L.). Therefore, in this study, integrated bioinformatics tools were used to identify and characterize the LecRLK gene family in barley. Based on the phylogenetic tree and domain organization, a total of 113 LecRLK genes were identified in the barley genome (referred to as HvlecRLK) corresponding to the LecRLK genes of Arabidopsis thaliana. These putative HvlecRLK genes were classified into three groups: 62 G-type LecRLKs, 1 C-type LecRLK, and 50 L-type LecRLKs. They were unevenly distributed across eight chromosomes, including one unknown chromosome, and were predominantly located in the plasma membrane (G-type HvlecRLK (96.8%), C-type HvlecRLK (100%), and L-type HvlecRLK (98%)). An analysis of motif composition and exon-intron configuration revealed remarkable homogeneity with the members of AtlecRLK. Notably, most of the HvlecRLKs (27 G-type, 43 L-type) have no intron, suggesting their rapid functionality. The Ka/Ks and syntenic analysis demonstrated that HvlecRLK gene pairs evolved through purifying selection and gene duplication was the major factor for the expansion of the HvlecRLK gene family. Exploration of gene ontology (GO) enrichment indicated that the identified HvlecRLK genes are associated with various cellular processes, metabolic pathways, defense mechanisms, kinase activity, catalytic activity, ion binding, and other essential pathways. The regulatory network analysis identified 29 transcription factor families (TFFs), with seven major TFFs including bZIP, C2H2, ERF, MIKC_MADS, MYB, NAC, and WRKY participating in the regulation of HvlecRLK gene functions. Most notably, eight TFFs were found to be linked to the promoter region of both L-type HvleckRLK64 and HvleckRLK86. The promoter cis-acting regulatory element (CARE) analysis of barley identified a total of 75 CARE motifs responsive to light responsiveness (LR), tissue-specific (TS), hormone responsiveness (HR), and stress responsiveness (SR). The maximum number of CAREs was identified in HvleckRLK11 (25 for LR), HvleckRLK69 (17 for TS), and HvleckRLK80 (12 for HR). Additionally, HvleckRLK14, HvleckRLK16, HvleckRLK33, HvleckRLK50, HvleckRLK52, HvleckRLK56, and HvleckRLK110 were predicted to exhibit higher responses in stress conditions. In addition, 46 putative miRNAs were predicted to target 81 HvlecRLK genes and HvlecRLK13 was the most targeted gene by 8 different miRNAs. Protein-protein interaction analysis demonstrated higher functional similarities of 63 HvlecRLKs with 7 Arabidopsis STRING proteins. Our overall findings provide valuable information on the LecRLK gene family which might pave the way to advanced research on the functional mechanism of the candidate genes as well as to develop new barley cultivars in breeding programs.

Molecular Identification and Characterization of Hevein Antimicrobial Peptide Genes in Two-Row and Six-Row Cultivars of Barley (Hordeum vulgare L.)

Heveins are one of the most important groups of plant antimicrobial peptides. So far, various roles in plant growth and development and in response to biotic and abiotic stresses have reported for heveins. The present study aimed to identify and characterize the hevein genes in two-row and six-row cultivars of barley. In total, thirteen hevein genes were identified in the genome of two-row and six-row cultivars of barley. The identified heveins were identical in two-row and six-row cultivars of barley and showed a high similarity with heveins from other plant species. The hevein coding sequences produced open reading frames (ORFs) ranged from 342 to 1002 bp. Most of the identified hevein genes were intronless, and the others had only one intron. The hevein ORFs produced proteins ranged from 113 to 333 amino acids. Search for conserved functional domains showed CBD and LYZ domains in barley heveins. All barley heveins comprised extracellular signal peptides ranged from 19 to 35 amino acids. The phylogenetic analysis divided barley heveins into two groups. The promoter analysis showed regulatory elements with different frequencies between two-row and six-row cultivars. These cis-acting elements included elements related to growth and development, hormone response, and environmental stresses. The expression analysis showed high expression level of heveins in root and reproductive organs of both two-row and six-row cultivars. The expression analysis also showed that barley heveins is induced by both biotic and abiotic stresses. The results of antimicrobial activity prediction showed the highest antimicrobial activity in CBD domain of barley heveins. The findings of the current study can improve our knowledge about the role of hevein genes in plant and can be used for future studies.

Cacao pod transcriptome profiling of seven genotypes identifies features associated with post-penetration resistance to Phytophthora palmivora

The oomycete Phytophthora palmivora infects the fruit of cacao trees (Theobroma cacao) causing black pod rot and reducing yields. Cacao genotypes vary in their resistance levels to P. palmivora, yet our understanding of how cacao fruit respond to the pathogen at the molecular level during disease establishment is limited. To address this issue, disease development and RNA-Seq studies were conducted on pods of seven cacao genotypes (ICS1, WFT, Gu133, Spa9, CCN51, Sca6 and Pound7) to better understand their reactions to the post-penetration stage of P. palmivora infection. The pod tissue-P. palmivora pathogen assay resulted in the genotypes being classified as susceptible (ICS1, WFT, Gu133 and Spa9) or resistant (CCN51, Sca6 and Pound7). The number of differentially expressed genes (DEGs) ranged from 1625 to 6957 depending on genotype. A custom gene correlation approach identified 34 correlation groups. De novo motif analysis was conducted on upstream promoter sequences of differentially expressed genes, identifying 76 novel motifs, 31 of which were over-represented in the upstream sequences of correlation groups and associated with gene ontology terms related to oxidative stress response, defense against fungal pathogens, general metabolism and cell function. Genes in one correlation group (Group 6) were strongly induced in all genotypes and enriched in genes annotated with defense-responsive terms. Expression pattern profiling revealed that genes in Group 6 were induced to higher levels in the resistant genotypes. An additional analysis allowed the identification of 17 candidate cis-regulatory modules likely to be involved in cacao defense against P. palmivora. This study is a comprehensive exploration of the cacao pod transcriptional response to P. palmivora spread after infection. We identified cacao genes, promoter motifs, and promoter motif combinations associated with post-penetration resistance to P. palmivora in cacao pods and provide this information as a resource to support future and ongoing efforts to breed P. palmivora-resistant cacao.

Cloning of maize chitinase 1 gene and its expression in genetically transformed rice to confer resistance against rice blast caused by Pyricularia oryzae

Fungal pathogens are one of the major reasons for biotic stress on rice (Oryza sativa L.), causing severe productivity losses every year. Breeding for host resistance is a mainstay of rice disease management, but conventional development of commercial resistant varieties is often slow. In contrast, the development of disease resistance by targeted genome manipulation has the potential to deliver resistant varieties more rapidly. The present study reports the first cloning of a synthetic maize chitinase 1 gene and its insertion in rice cv. (Basmati 385) via Agrobacterium-mediated transformation to confer resistance to the rice blast pathogen, Pyricularia oryzae. Several factors for transformation were optimized; we found that 4-week-old calli and an infection time of 15 minutes with Agrobacterium before colonization on co-cultivation media were the best-suited conditions. Moreover, 300 μM of acetosyringone in co-cultivation media for two days was exceptional in achieving the highest callus transformation frequency. Transgenic lines were analyzed using molecular and functional techniques. Successful integration of the gene into rice lines was confirmed by polymerase chain reaction with primer sets specific to chitinase and hpt genes. Furthermore, real-time PCR analysis of transformants indicated a strong association between transgene expression and elevated levels of resistance to rice blast. Functional validation of the integrated gene was performed by a detached leaf bioassay, which validated the efficacy of chitinase-mediated resistance in all transgenic Basmati 385 plants with variable levels of enhanced resistance against the P. oryzae. We concluded that overexpression of the maize chitinase 1 gene in Basmati 385 improved resistance against the pathogen. These findings will add new options to resistant germplasm resources for disease resistance breeding. The maize chitinase 1 gene demonstrated potential for genetic improvement of rice varieties against biotic stresses in future transformation programs.

Analysis of Germin-like protein genes family in Vitis vinifera (VvGLPs) using various in silico approaches

Abstract Germin-like proteins (GLPs) play an important role against various stresses. Vitis vinifera L. genome contains 7 GLPs; many of them are functionally unexplored. However, the computational analysis may provide important new insight into their function. Currently, physicochemical properties, subcellular localization, domain architectures, 3D structures, N-glycosylation & phosphorylation sites, and phylogeney of the VvGLPs were investigated using the latest computational tools. Their functions were predicted using the Search tool for the retrieval of interacting genes/proteins (STRING) and Blast2Go servers. Most of the VvGLPs were extracellular (43%) in nature but also showed periplasmic (29%), plasma membrane (14%), and mitochondrial- or chloroplast-specific (14%) expression. The functional analysis predicted unique enzymatic activities for these proteins including terpene synthase, isoprenoid synthase, lipoxygenase, phosphate permease, receptor kinase, and hydrolases generally mediated by Mn+ cation. VvGLPs showed similarity in the overall structure, shape, and position of the cupin domain. Functionally, VvGLPs control and regulate the production of secondary metabolites to cope with various stresses. Phylogenetically VvGLP1, -3, -4, -5, and VvGLP7 showed greater similarity due to duplication while VvGLP2 and VvGLP6 revealed a distant relationship. Promoter analysis revealed the presence of diverse cis-regulatory elements among which CAAT box, MYB, MYC, unnamed-4 were common to all of them. The analysis will help to utilize VvGLPs and their promoters in future food programs by developing resistant cultivars against various biotic (Erysiphe necator and in Powdery Mildew etc.) and abiotic (Salt, drought, heat, dehydration, etc.) stresses.

Nodule-specific Cu+ -chaperone NCC1 is required for symbiotic nitrogen fixation in Medicago truncatula root nodules

Cu+ -chaperones are a diverse group of proteins that allocate Cu+ ions to specific copper proteins, creating different copper pools targeted to specific physiological processes. Symbiotic nitrogen fixation carried out in legume root nodules indirectly requires relatively large amounts of copper, for example for energy delivery via respiration, for which targeted copper deliver systems would be required. MtNCC1 is a nodule-specific Cu+ -chaperone encoded in the Medicago truncatula genome, with a N-terminus Atx1-like domain that can bind Cu+ with picomolar affinities. MtNCC1 is able to interact with nodule-specific Cu+ -importer MtCOPT1. MtNCC1 is expressed primarily from the late infection zone to the early fixation zone and is located in the cytosol, associated with plasma and symbiosome membranes, and within nuclei. Consistent with its key role in nitrogen fixation, ncc1 mutants have a severe reduction in nitrogenase activity and a 50% reduction in copper-dependent cytochrome c oxidase activity. A subset of the copper proteome is also affected in the ncc1 mutant nodules. Many of these proteins can be pulled down when using a Cu+ -loaded N-terminal MtNCC1 moiety as a bait, indicating a role in nodule copper homeostasis and in copper-dependent physiological processes. Overall, these data suggest a pleiotropic role of MtNCC1 in copper delivery for symbiotic nitrogen fixation.

Genome-wide identification and characterization of ABC transporter superfamily in the legume Cajanus cajan

Plant ATP-binding cassette (ABC) protein family is the largest multifunctional highly conserved protein superfamily that transports diverse substrates across biological membranes by the hydrolysis of ATP and is also the part of the several other biological processes like cellular detoxification, growth and development, stress biology, and signaling processes. In the agriculturally important legume crop Cajanus cajan, a genome-wide identification and characterization of the ABC gene family was carried out. A total of 159 ABC genes were identified that belong to eight canonical classes CcABCA to CcABCG and CcABCI based on the phylogenetic analysis. The number of genes was highest in CcABCG followed by CcABCC and CcABCB class. A total of 85 CcABC genes were found on 11 chromosomes and 74 were found on scaffold. Tandem duplication was the major driver of CcABC gene family expansion. The dN/dS ratio revealed the purifying selection. The phylogenetic analysis revealed class-specific eight superclades which reflect their functional importance. The largest clade was found to be CcABCG which reflects their functional significance. CcABC proteins were mainly basic in nature and found to be localized in the plasma membrane. The secondary structure prediction revealed the dominance of α-helix. The canonical transmembrane and nucleotide binding domain, signature motif LSSGQ, Walker A, Walker B region, and Q loop were also identified. A class-specific exon-intron pattern was also observed. In addition to core elements, different cis-acting regulatory elements like stress, hormone, and cellular responsive were also identified. Expression profiling of CcABC genes at various developmental stages of different anatomical tissues was performed and it was noticed that CcABCF3, CcABCF4, CcABCF5, CcABCG66, and CcABCI3 had the highest expression. The results of the current study endow us with the further functional analysis of Cajanus ABC in the future.

Genome-Wide Identification and Characterization of Major RNAi Genes Highlighting Their Associated Factors in Cowpea (Vigna unguiculata (L.) Walp.)

In different regions of the world, cowpea (Vigna unguiculata (L.) Walp.) is an important vegetable and an excellent source of protein. It lessens the malnutrition of the underprivileged in developing nations and has some positive effects on health, such as a reduction in the prevalence of cancer and cardiovascular disease. However, occasionally, certain biotic and abiotic stresses caused a sharp fall in cowpea yield. Major RNA interference (RNAi) genes like Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are essential for the synthesis of their associated factors like domain, small RNAs (sRNAs), transcription factors, micro-RNAs, and cis-acting factors that shield plants from biotic and abiotic stresses. In this study, applying BLASTP search and phylogenetic tree analysis with reference to the Arabidopsis RNAi (AtRNAi) genes, we discovered 28 VuRNAi genes, including 7 VuDCL, 14 VuAGO, and 7 VuRDR genes in cowpea. We looked at the domains, motifs, gene structures, chromosomal locations, subcellular locations, gene ontology (GO) terms, and regulatory factors (transcription factors, micro-RNAs, and cis-acting elements (CAEs)) to characterize the VuRNAi genes and proteins in cowpea in response to stresses. Predicted VuDCL1, VuDCL2(a, b), VuAGO7, VuAGO10, and VuRDR6 genes might have an impact on cowpea growth, development of the vegetative and flowering stages, and antiviral defense. The VuRNAi gene regulatory features miR395 and miR396 might contribute to grain quality improvement, immunity boosting, and pathogen infection resistance under salinity and drought conditions. Predicted CAEs from the VuRNAi genes might play a role in plant growth and development, improving grain quality and production and protecting plants from biotic and abiotic stresses. Therefore, our study provides crucial information about the functional roles of VuRNAi genes and their associated components, which would aid in the development of future cowpeas that are more resilient to biotic and abiotic stress. The manuscript is available as a preprint at this link: doi:10.1101/2023.02.15.528631v1.

Genome-wide identification, classification, and characterization of lectin gene superfamily in sweet orange (Citrus sinensis L.)

Lectins are sugar-binding proteins found abundantly in plants. Lectin superfamily members have diverse roles, including plant growth, development, cellular processes, stress responses, and defense against microbes. However, the genome-wide identification and functional analysis of lectin genes in sweet orange (Citrus sinensis L.) remain unexplored. Therefore, we used integrated bioinformatics approaches (IBA) for in-depth genome-wide identification, characterization, and regulatory factor analysis of sweet orange lectin genes. Through genome-wide comparative analysis, we identified a total of 141 lectin genes distributed across 10 distinct gene families such as 68 CsB-Lectin, 13 CsLysin Motif (LysM), 4 CsChitin-Bind1, 1 CsLec-C, 3 CsGal-B, 1 CsCalreticulin, 3 CsJacalin, 13 CsPhloem, 11 CsGal-Lec, and 24 CsLectinlegB.This classification relied on characteristic domain and phylogenetic analysis, showing significant homology with Arabidopsis thaliana's lectin gene families. A thorough analysis unveiled common similarities within specific groups and notable variations across different protein groups. Gene Ontology (GO) enrichment analysis highlighted the predicted genes' roles in diverse cellular components, metabolic processes, and stress-related regulation. Additionally, network analysis of lectin genes with transcription factors (TFs) identified pivotal regulators like ERF, MYB, NAC, WRKY, bHLH, bZIP, and TCP. The cis-acting regulatory elements (CAREs) found in sweet orange lectin genes showed their roles in crucial pathways, including light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and more. These findings will aid in the in-depth molecular examination of these potential genes and their regulatory elements, contributing to targeted enhancements of sweet orange species in breeding programs.

Exploring the genetic landscape of nitrogen uptake in durum wheat: genome-wide characterization and expression profiling of NPF and NRT2 gene families

Nitrate uptake by plants primarily relies on two gene families: Nitrate transporter 1/peptide transporter (NPF) and Nitrate transporter 2 (NRT2). Here, we extensively characterized the NPF and NRT2 families in the durum wheat genome, revealing 211 NPF and 20 NRT2 genes. The two families share many Cis Regulatory Elements (CREs) and Transcription Factor binding sites, highlighting a partially overlapping regulatory system and suggesting a coordinated response for nitrate transport and utilization. Analyzing RNA-seq data from 9 tissues and 20 cultivars, we explored expression profiles and co-expression relationships of both gene families. We observed a strong correlation between nucleotide variation and gene expression within the NRT2 gene family, implicating a shared selection mechanism operating on both coding and regulatory regions. Furthermore, NPF genes showed highly tissue-specific expression profiles, while NRT2s were mainly divided in two co-expression modules, one expressed in roots (NAR2/NRT3 dependent) and the other induced in anthers and/ovaries during maturation. Our evidences confirmed that the majority of these genes were retained after small-scale duplication events, suggesting a neo- or sub-functionalization of many NPFs and NRT2s. Altogether, these findings indicate that the expansion of these gene families in durum wheat could provide valuable genetic variability useful to identify NUE-related and candidate genes for future breeding programs in the context of low-impact and sustainable agriculture.

Isolation and Characterization of Erianthus arundinaceus Phosphate Transporter 1 (PHT1) Gene Promoter and 5' Deletion Analysis of Transcriptional Regulation Regions under Phosphate Stress in Transgenic Tobacco

Phosphorus deficiency highly interferes with plant growth and development. Plants respond to persistent P deficiency by coordinating the expression of genes involved in the alleviation of stress. Promoters of phosphate transporter genes are a great choice for the development of genetically modified plants with enhanced phosphate uptake abilities, which improve crop yields in phosphate-deficient soils. In our previous study, the sugarcane phosphate transporter PHT1;2 gene showed a significantly high expression under salinity stress. In this study, the Erianthus arundinaceus EaPHT1;2 gene was isolated and characterized using various in silico tools. The deduced 542 amino acid residues have 10 transmembrane domains, with a molecular weight and isoelectric point of 58.9 kDa and 9.80, respectively. They displayed 71-96% similarity with Arabidopsis thaliana, Zea mays, and the Saccharum hybrid. To elucidate the function of the 5' regulatory region, the 1.1 kb promoter was isolated and validated in tobacco transgenics under Pi stress. The EaPHT1;2 promoter activity was detected using a β-glucuronidase (GUS) assay. The EaPHT1;2 promoter showed 3- to 4.2-fold higher expression than the most widely used CaMV35S promoter. The 5' deletion analysis with and without 5' UTRs revealed a small-sized 374 bp fragment with the highest promoter activity among 5' truncated fragments, which was 2.7 and 4.2 times higher than the well-used CaMV35S promoter under normal and Pi deprivation conditions, respectively. The strong and short promoter of EaPHT1;2 with 374 bp showed significant expression in low-Pi-stress conditions and it could be a valuable source for the development of stress-tolerant transgenic crops.

Genome-wide identification of glutamate receptor-like gene family in soybean

Glutamate receptor-like genes (GLRs) are essential in the growth and development of plants and many physiological and biochemical processes; however, related information in soybean is lacking. In this study, 105 GLRs, including 67 Glycine soja and 38 Glycine max GLRs, were identified and divided into two clades (Clades II and III) according to their phylogenetic relationships. GLR members in the same branch had a relatively conservative motif composition and genetic structure. Furthermore, the soybean GLR family mainly experienced purification selection during evolution. Cis-acting element analysis, gene ontology, and Kyoto Encyclopedia of Genes and Genomic annotations indicated the complexity of the gene regulation and functional diversity of the soybean GLR. Moreover, transcriptome data analysis showed that these GLRs had different expression profiles in different tissues, and Clade III members had higher and more common expression patterns. Additionally, the expression profiles under jasmonic acid treatment and salt stress indicate that the GLR participated in the jasmonic acid signaling pathway and plays a role in salt treatment. This study provides information for a comprehensive understanding of the soybean GLR family and a reference for further functional research and genetic improvement.

Comprehensive genomic characterisation of the NAC transcription factor family and its response to drought stress in Eucommia ulmoides

The NAC transcription factor family enhances plant adaptation to environmental challenges by participating in signalling pathways triggered by abiotic stressors and hormonal cues. We identified 69 NAC genes in the Eucommia ulmoides genome and renamed them according to their chromosomal distribution. These EuNAC proteins were clustered into 13 sub-families and distributed on 16 chromosomes and 2 scaffolds. The gene structures suggested that the number of exons varied from two to eight among these EuNACs, with a multitude of them containing three exons. Duplicated events resulted in a large gene family; 12 and four pairs of EuNACs were the result of segmental and tandem duplicates, respectively. The drought-stress response pattern of 12 putative EuNACs was observed under drought treatment, revealing that these EuNACs could play crucial roles in mitigating the effects of drought stress responses and serve as promising candidate genes for genetic engineering aimed at enhancing the drought stress tolerance of E. ulmoides. This study provides insight into the evolution, diversity, and characterisation of NAC genes in E. ulmoides and will be helpful for future characterisation of putative EuNACs associated with water deficit.

A genome‑wide approach to the systematic and comprehensive analysis of LIM gene family in sorghum (Sorghum bicolor L.)

The LIM domain-containing proteins are dominantly found in plants and play a significant role in various biological processes such as gene transcription as well as actin cytoskeletal organization. Nevertheless, genome-wide identification as well as functional analysis of the LIM gene family have not yet been reported in the economically important plant sorghum (Sorghum bicolor L.). Therefore, we conducted an in silico identification and characterization of LIM genes in S. bicolor genome using integrated bioinformatics approaches. Based on phylogenetic tree analysis and conserved domain, we identified five LIM genes in S. bicolor (SbLIM) genome corresponding to Arabidopsis LIM (AtLIM) genes. The conserved domain, motif as well as gene structure analyses of the SbLIM gene family showed the similarity within the SbLIM and AtLIM members. The gene ontology (GO) enrichment study revealed that the candidate LIM genes are directly involved in cytoskeletal organization and various other important biological as well as molecular pathways. Some important families of regulating transcription factors such as ERF, MYB, WRKY, NAC, bZIP, C2H2, Dof, and G2-like were detected by analyzing their interaction network with identified SbLIM genes. The cis-acting regulatory elements related to predicted SbLIM genes were identified as responsive to light, hormones, stress, and other functions. The present study will provide valuable useful information about LIM genes in sorghum which would pave the way for the future study of functional pathways of candidate SbLIM genes as well as their regulatory factors in wet-lab experiments.

An insight into Pisum sativum HSF gene family-Genome-wide identification, phylogenetic, expression, and analysis of transactivation potential of pea heat shock transcription factor

Field pea (Pisum sativum L, 2n = 14) is a popular temperate legume with high economic value. Heat shock factors (HSFs) are the core element in the regulatory mechanism of heat stress responses. HSFs in pea (P. sativum) have not been characterized and their role remains unclear in different abiotic stresses. To address this knowledge gap, the current study aimed to characterize the HSF gene family in pea. We identified 38 PsHsf members in P. sativum, which are distributed on the seven chromosomes, and based on phylogenetic analysis, we classified them into three representative classes i.e. A, B, and C. Conserved motif and gene structure analysis confirmed a high degree of similarity among the members of the same class. Additionally, identified cis-acting regulatory elements (CAREs) related to abiotic responses, development, growth, and hormone signaling provides crucial insights into the regulatory mechanisms of PsHsfs. Our research revealed instances of gene duplication in PsHsf gene family, suggesting that this mechanism could be driving the expansion of the PsHsf gene family. Moreover, Expression analysis of PsHsfs exhibited upregulation under heat stress (HS), salt stress (SS), and drought stress (DS) showing their phenomenal role in stress conditions. PsHsfs protein interaction network suggested their involvement in stress-responsive mechanisms. Further transactivation potential was checked for spliced variant of PsHsfA2a (PsHsfA2aI, PsHsfA2aII, and PsHsfA2aIII), PsHsfA3, PsHsfA6b, PsHsfA9, PsHsfB1a, and PsHsfB2a. Overall, these findings provide valuable insight into the evolutionary relationship of PsHsf gene family and their role in abiotic stress responses.

Genome-wide identification and characterization of glutathione S-transferase gene family in quinoa (Chenopodium quinoa Willd.)

The present investigation was envisaged for large scale in-silico genome wide identification and characterization of glutathione S-transferases (GSTs) in Chenopodium quinoa. In this study, a total of 120 GST genes (CqGSTs) were identified and divided into 11 classes of which tau and phi were highest in numbers. The average protein length of protein was found to be 279.06 with their corresponding average molecular weight of 31,819.4 kDa. The subcellular localization analysis results showed that proteins were centrally localized in the cytoplasm followed by chloroplast, mitochondria and plastids. Structural analysis revealed the presence of 2 -14 exons in CqGST genes. Most of the proteins possessed two exon one intron organization. MEME analysis identified 15 significantly conserved motifs with a width of 6-50 amino acids. Motifs 1, 3, 2, 5, 6, 8, 9 and 13 were found specifically in tau class family; motifs 3, 4, 5, 6, 7 and 9 were found in phi class gene family, while motifs 3, 4, 13 and 14 were found in metaxin class. Multiple sequence alignment revealed highly conserved N-terminus with active site serine (Ser; S) or cysteine (Cys; C) residue for the activation of GSH binding and GST catalytic activity. The gene loci were found to be unevenly distributed across 18 different chromosomes with a maximum of 17 genes located on chromosome number 7. Dominance of alpha helix was followed by coil, extended strand and beta turns. Gene duplication analysis revealed that segmental duplication and purifying type selection were highest in number and found to be main source of expansion of GST gene family. Cis acting regulatory elements analysis showed the presence of 21 different elements involved in stress, hormone and light response and cellular development. The evolutionary relationship of CqGST proteins carried out using maximum likelihood method revealed that all the tau and phi class GSTs were closely associated with those of G. max, O. sativa and A. thaliana. Molecular docking of GST molecules with the fungicide metalaxyl showed that the CqGSTF1 had the lowest binding energy. The comprehensive study of CqGST gene family in quinoa provides groundwork for further functional analysis of CqGST genes in the species at molecular level and has potential applications in plant breeding.

Characterization of myo-inositol oxygenase from rice (OsMIOX): influence of salinity stress in different indica rice cultivars

Myo-inositol oxygenase (MIOX), the only catabolic enzyme of the inositol pathway, catalyzes conversion of myo-inositol to D-GlcA (glucuronic acid). The present study encompasses bioinformatic analysis of MIOX gene across phylogenetically related plant lineages and representative animal groups. Comparative motif analysis of the MIOX gene(s) across various plant groups suggested existence of abiotic- stress related cis-acting elements such as, DRE, MYB, MYC, STRE, MeJa among others. A detailed analysis revealed a single isoform of MIOX gene, located in chromosome 6 of indica rice (Oryza sativa) with an open reading frame of 938 bp coding for 308 amino acids producing a protein of ~ 35 kD. Secondary structure prediction of the protein gave the predicted number of 144 alpha helices and 154 random coils. The three-dimensional structure suggested it to be a monomeric protein with a single domain. Bacterial overexpression of the protein, purification and enzyme assay showed optimal catalytic activity at pH 7.5-8 at an optimal temperature of 37 °C with Michaelis constant of 40.92 mM. The range of Km was determined as 22.74-28.7 mM and the range of Vmax was calculated as 3.51-3.6 µM/min, respectively. Four salt-tolerant and salt-sensitive rice cultivars displayed differential gene expression of OsMIOX at different time points in different tissues under salinity and drought stress as observed from qRT-PCR data, microarray results and protein expression profile in immunoblot analysis. Gel volumetric analysis confirmed a very high expression of MIOX in roots and leaves on 7th day following germination. Microarray data showed high expression of MIOX at all developmental stages including seedling growth and reproduction. These data suggest that OsMIOX might have a role to play in rice abiotic stress responses mediated through the myo-inositol oxidation pathway.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01340-6.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory element analyses in sunflower (Helianthus annuus)

RNA interference (RNAi) regulates a variety of eukaryotic gene expressions that are engaged in response to stress, growth, and the conservation of genomic stability during developmental phases. It is also intimately connected to the post-transcriptional gene silencing (PTGS) process and chromatin modification levels. The entire process of RNA interference (RNAi) pathway gene families mediates RNA silencing. The main factors of RNA silencing are the Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) gene families. To the best of our knowledge, genome-wide identification of RNAi gene families like DCL, AGO, and RDR in sunflower (Helianthus annuus) has not yet been studied despite being discovered in some species. So, the goal of this study is to find the RNAi gene families like DCL, AGO, and RDR in sunflower based on bioinformatics approaches. Therefore, we accomplished an inclusive in silico investigation for genome-wide identification of RNAi pathway gene families DCL, AGO, and RDR through bioinformatics approaches such as (sequence homogeneity, phylogenetic relationship, gene structure, chromosomal localization, PPIs, GO, sub-cellular localization). In this study, we have identified five DCL (HaDCLs), fifteen AGO (HaAGOs), and ten RDR (HaRDRs) in the sunflower genome database corresponding to the RNAi genes of model plant Arabidopsis thaliana based on genome-wide analysis and a phylogenetic method. The analysis of the gene structure that contains exon-intron numbers, conserved domain, and motif composition analyses for all HaDCL, HaAGO, and HaRDR gene families indicated almost homogeneity among the same gene family. The protein-protein interaction (PPI) network analysis illustrated that there exists interconnection among identified three gene families. The analysis of the Gene Ontology (GO) enrichment showed that the detected genes directly contribute to the RNA gene-silencing and were involved in crucial pathways. It was observed that the cis-acting regulatory components connected to the identified genes were shown to be responsive to hormone, light, stress, and other functions. That was found in HaDCL, HaAGO, and HaRDR genes associated with the development and growth of plants. Finally, we are able to provide some essential information about the components of sunflower RNA silencing through our genome-wide comparison and integrated bioinformatics analysis, which open the door for further research into the functional mechanisms of the identified genes and their regulatory elements.

Diversity and expression analysis of ZIP transporters and associated metabolites under zinc and iron stress in Capsicum

The members of ZRT, IRT-like protein (ZIP) family are involved in the uptake and transportation of several metal ions. Here, we report a comprehensive identification of ZIP transporter genes from Capsicum annuum, C. chinense, and C. baccatum, and their expression analysis under Zn and Fe stress. Changes in root morphology and differential accumulation of several metabolites from sugars, amino acids, carboxylic acids, and fatty acids in root and leaf tissues of plants in the absence of Zn and Fe were observed. Further, metabolites such as L-aspartic acid, 2-ketoglutaric acids, β-L-fucopyranose, quininic acid, chlorogenic acid, and aucubin were significantly upregulated in root and leaf tissues under Zn/Fe deprived conditions. qRT-PCR analysis of 17 CaZIPs in different tissues revealed tissue-specific expression of CaZIP1-2, CaZIP4-8, CaZIP13, and CaZIP16-17 under normal conditions. However, the absence of Zn and Fe significantly induced the expression of CaZIP4-5, CaZIP7-9, and CaZIP14 genes in root and leaf tissues. Additionally, in the absence of Fe, upregulation of CaZIP4-5 and CaZIP8 and increased uptake of mineral elements Cu, Zn, Mg, P, and S were observed in roots, suggesting their potential role in metal-ion uptake in Capsicum. The identified genes provide the basis for future studies of mineral uptake and their biofortification to increase the nutritional values in Capsicum.

Cloning of Three Cytokinin Oxidase/Dehydrogenase Genes in Bambusa oldhamii

Cytokinin oxidase/dehydrogenase (CKX) catalyzes the irreversible breakdown of active cytokinins, which are a class of plant hormones that regulate cell division. According to conserved sequences of CKX genes from monocotyledons, PCR primers were designed to synthesize a probe for screening a bamboo genomic library. Cloned results of three genes encoding cytokinin oxidase were named as follows: BoCKX1, BoCKX2, and BoCKX3. In comparing the exon-intron structures among the above three genes, there are three exons and two introns in BoCKX1 and BoCKX3 genes, whereas BoCKX2 contains four exons and three introns. The amino acid sequence of BoCKX2 protein shares 78% and 79% identity with BoCKX1 and BoCKX3 proteins, respectively. BoCKX1 and BoCKX3 genes are particularly closely related given that the amino acid and nucleotide sequence identities are more than 90%. These three BoCKX proteins carried putative signal peptide sequences typical of secretion pathway, and a GHS-motif was found at N-terminal flavin adenine dinucleotide (FAD) binding domain, suggesting that BoCKX proteins might covalently conjugate with an FAD cofactor through a predicted histidine residue.

Genome-Wide Characterization and Gene Expression Analyses of Malate Dehydrogenase (MDH) Genes in Low-Phosphorus Stress Tolerance of Chinese Fir (Cunninghamia lanceolata)

Malate dehydrogenase (MDH) genes play vital roles in developmental control and environmental stress tolerance in sessile plants by modulating the organic acid-malic acid level. However, MDH genes have not yet been characterized in gymnosperm, and their roles in nutrient deficiency are largely unexplored. In this study, 12 MDH genes were identified in Chinese fir (Cunninghamia lanceolata), namely, ClMDH-1, -2, -3, …, and -12. Chinese fir is one of the most abundant commercial timber trees in China, and low phosphorus has limited its growth and production due to the acidic soil of southern China. According to the phylogenetic analysis, MDH genes were classified into five groups, and Group 2 genes (ClMDH-7, -8, -9, and 10) were only found to be present in Chinese fir but not in Arabidopsis thaliana and Populus trichocarpa. In particular, the Group 2 MDHs also had specific functional domains-Ldh_1_N (malidase NAD-binding functional domain) and Ldh_1_C (malate enzyme C-terminal functional domain)-indicating a specific function of ClMDHs in the accumulation of malate. All ClMDH genes contained the conserved MDH gene characteristic functional domains Ldh_1_N and Ldh_1_C, and all ClMDH proteins exhibited similar structures. Twelve ClMDH genes were identified from eight chromosomes, involving fifteen ClMDH homologous gene pairs, each with a Ka/Ks ratio of <1. The analysis of cis-elements, protein interactions, and transcription factor interactions of MDHs showed that the ClMDH gene might play a role in plant growth and development, and in response to stress mechanisms. The results of transcriptome data and qRT-PCR validation based on low-phosphorus stress showed that ClMDH1, ClMDH6, ClMDH7, ClMDH2, ClMDH4, ClMDH5, ClMDH10 and ClMDH11 were upregulated under low-phosphorus stress and played a role in the response of fir to low-phosphorus stress. In conclusion, these findings lay a foundation for further improving the genetic mechanism of the ClMDH gene family in response to low-phosphorus stress, exploring the potential function of this gene, promoting the improvement of fir genetics and breeding, and improving production efficiency.

Genome-wide identification and analysis of the evolution and expression pattern of the HVA22 gene family in three wild species of tomatoes

Wild tomato germplasm is a valuable resource for improving biotic and abiotic stresses in tomato breeding. The HVA22 is widely present in eukaryotes and involved in growth and development as well as stress response, such as cold, salt, drought, and biotic stress. In the present study, we identified 45 HVA22 genes in three wild species of tomatoes. The phylogenetic relationships, gene localization to chromosomes, gene structure, gene collinearity, protein interactions, and cis-acting element prediction of all 45 HVA22 genes (14 in Solanum pennellii, 15 in S. pimpinellifolium, and 16 in S. lycopersicoides) were analyzed. The phylogenetic analysis showed that the all HVA22 proteins from the family Solanaceae were divided into three branches. The identified 45 HVA22 genes were grouped into four subfamilies, which displayed similar number of exons and expanded in a fragmentary replication manner. The distribution of HVA22 genes on the chromosomes of the three wild tomato species was also highly similar. RNA-seq and qRT-PCR revealed that HVA22 genes were expressed in different tissues and induced by drought, salt, and phytohormone treatments. These results might be useful for explaining the evolution, expression patterns, and functional divergence of HVA22 genes in Lycopersicon.

Genome-wide identification and characterization of parthenocarpic fruit set-related gene homologs in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.), a major horticultural crop, in the family Cucurbitaceae is grown and consumed globally. Parthenocarpy is an ideal trait for many fruit and vegetables which produces seedless fruit desired by consumers. The seedlessness occurs when fruit develops without fertilization which can be either natural or induced. So far, a limited number of genes regulating parthenocarpic fruit set have been reported in several fruit or vegetable crops, most of which are involved in hormone biosynthesis or signalling. Although parthenocarpic cucumber has been widely used in commercial production for a long time; its genetic basis is not well understood. In this study, we retrieved thirty five parthenocarpy fruit-set related genes (PRGs) from bibliomic data in various plants. Thirty-five PRG homologs were identified in the cucumber genome via homology-based search. An in silico analysis was performed on phylogenetic tree, exon-intron structure, cis-regulatory elements in the promoter region, and conserved domains of their deduced proteins, which provided insights into the genetic make-up of parthenocarpy-related genes in cucumber. Simple sequence repeat (SSR) sequences were mined in these PRGs, and 31 SSR markers were designed. SSR genotyping identified three SSRs in two polymorphic genes. Quantitative real-time PCR of selected genes was conducted in five cucumber lines with varying degrees of parthenocarpic fruit set capacities, which revealed possible association of their expression with parthenocarpy. The results revealed that homologs CsWD40 and CsPIN-4 could be considered potential genes for determination of parthenocarpy as these genes showed parental polymorphism and differential gene expression in case of parthenocarpic and non-parthenocarpic parents.

Overexpression of OsGF14C enhances salinity tolerance but reduces blast resistance in rice

High-salinity and blast disease are two major stresses that cause dramatic yield loss in rice production. GF14 (14-3-3) genes have been reported to play important role in biotic and abiotic stresses in plants. However, the roles of OsGF14C remain unknown. To understand the functions and regulatory mechanisms of OsGF14C in regulating salinity tolerance and blast resistance in rice, we have conducted OsGF14C-overexpressing transgenic experiments in the present study. Our results showed that overexpression of OsGF14C enhanced salinity tolerance but reduced blast resistance in rice. The enhanced salinity tolerance is related to the reduction of methylglyoxal and Na⁺ uptake instead of exclusion or compartmentation and the negative role of OsGF14C in blast resistance is associated with the suppression of OsGF14E, OsGF14F and PR genes. Our results together with the results from the previous studies suggest that the lipoxygenase gene LOX2 which is regulated by OsGF14C may play roles in coordinating salinity tolerance and blast resistance in rice. The current study for the first time revealed the possible roles of OsGF14C in regulating salinity tolerance and blast resistance in rice, and laid down a foundation for further functional study and crosstalk regulation between salinity and blast resistance in rice.

NF-YB family transcription factors in Arabidopsis: Structure, phylogeny, and expression analysis in biotic and abiotic stresses

Nuclear factor-Y (NF-Y) transcription factors (TFs) are conserved heterotrimeric complexes present and widespread across eukaryotes. Three main subunits make up the structural and functional aspect of the NF-Y TFs: NF-YA, NF-YB and NF-YC, which bind to the conserved CCAAT- box of the promoter region of specific genes, while also interacting with each other, thereby forming myriad combinations. The NF-YBs are expressed differentially in various tissues and plant development stages, likely impacting many of the cellular processes constitutively and under stress conditions. In this study, ten members of NF-YB family from Arabidopsis thaliana were identified and expression profiles were mined from microarray data under different biotic and abiotic conditions, revealing key insights into the involvement of this class of proteins in the cellular and biological processes in Arabidopsis. Analysis of cis-acting regulatory elements (CAREs) indicated the presence of abiotic and biotic stress-related transcription factor binding sites (TFBs), shedding light on the multifaceted roles of these TFs. Microarray data analysis inferred distinct patterns of expression in various tissues under differing treatments such as drought, cold and heat stress as well as bacterial, fungal, and viral stress, indicating their likelihood of having an expansive range of regulatory functions under native and stressed conditions; while quantitative real-time PCR (qRT-PCR) based expression analysis revealed that these TFs get real-time-modulated in a stress dependent manner. This study, overall, provides an understanding of the AtNF-YB family of TFs in their regulation and participation in various morphogenetic and defense- related pathways and can provide insights for development of transgenic plants for trait dependent studies.

Integration of meta-analysis, machine learning and systems biology approach for investigating the transcriptomic response to drought stress in Populus species

In Populus, drought is a major problem affecting plant growth and development which can be closely reflected by corresponding transcriptomic changes. Nevertheless, how these changes in Populus are not fully understood. Here, we first used meta-analysis and machine learning methods to identify water stress-responsive genes and then performed a systematic approach to discover important gene networks. Our analysis revealed that large transcriptional variations occur during drought stress. These changes were more associated with the response to stress, cellular catabolic process, metabolic pathways, and hormone-related genes. The differential gene coexpression analysis highlighted two acetyltransferase NATA1-like and putative cytochrome P450 genes that have a special contribution in response to drought stress. In particular, the findings showed that MYBs and MAPKs have a prominent role in the drought stress response that could be considered to improve the drought tolerance of Populus. We also suggest ARF2-like and PYL4-like genes as potential markers for use in breeding programs. This study provides a better understanding of how Populus responses to drought that could be useful for improving tolerance to stress in Populus.

The NAC transcription factor family in Eucommia ulmoides: Genome-wide identification, characterization, and network analysis in relation to the rubber biosynthetic genes

The NAC transcription factor family is a large plant gene family, participating in plant growth and development, secondary metabolite synthesis, biotic and abiotic stresses responses, and hormone signaling. Eucommia ulmoides is a widely planted economic tree species in China that can produce trans-polyisoprene: Eucommia rubber (Eu-rubber). However, genome-wide identification of the NAC gene family has not been reported in E. ulmoides. In this study, 71 NAC proteins were identified based on genomic database of E. ulmoides. Phylogenetic analysis showed that the EuNAC proteins were distributed in 17 subgroups based on homology with NAC proteins in Arabidopsis, including the E. ulmoides-specific subgroup Eu_NAC. Gene structure analysis suggested that the number of exons varied from 1 to 7, and multitudinous EuNAC genes contained two or three exons. Chromosomal location analysis revealed that the EuNAC genes were unevenly distributed on 16 chromosomes. Three pairs of genes of tandem duplicates genes and 12 segmental duplications were detected, which indicated that segmental duplications may provide the primary driving force of expansion of EuNAC. Prediction of cis-regulatory elements indicated that the EuNAC genes were involved in development, light response, stress response and hormone response. For the gene expression analysis, the expression levels of EuNAC genes in various tissues were quite different. To explore the effect of EuNAC genes on Eu-rubber biosynthesis, a co-expression regulatory network between Eu-rubber biosynthesis genes and EuNAC genes was constructed, which indicated that six EuNAC genes may play an important role in the regulation of Eu-rubber biosynthesis. In addition, this six EuNAC genes expression profiles in E. ulmoides different tissues were consistent with the trend in Eu-rubber content. Quantitative real-time PCR analysis showed that EuNAC genes were responsive to different hormone treatment. These results will provide a useful reference for further studies addressing the functional characteristics of the NAC genes and its potential role in Eu-rubber biosynthesis.

Identification, characterization of Apyrase (APY) gene family in rice (Oryza sativa) and analysis of the expression pattern under various stress conditions

Apyrase (APY) is a nucleoside triphosphate (NTP) diphosphohydrolase (NTPDase) which is a member of the superfamily of guanosine diphosphatase 1 (GDA1)-cluster of differentiation 39 (CD39) nucleoside phosphatase. Under various circumstances like stress, cell growth, the extracellular adenosine triphosphate (eATP) level increases, causing a detrimental influence on cells such as cell growth retardation, ROS production, NO burst, and apoptosis. Apyrase hydrolyses eATP accumulated in the extracellular membrane during stress, wounds, into adenosine diphosphate (ADP) and adenosine monophosphate (AMP) and regulates the stress-responsive pathway in plants. This study was designed for the identification, characterization, and for analysis of APY gene expression in Oryza sativa. This investigation discovered nine APYs in rice, including both endo- and ecto-apyrase. According to duplication event analysis, in the evolution of OsAPYs, a significant role is performed by segmental duplication. Their role in stress control, hormonal responsiveness, and the development of cells is supported by the corresponding cis-elements present in their promoter regions. According to expression profiling by RNA-seq data, the genes were expressed in various tissues. Upon exposure to a variety of biotic as well as abiotic stimuli, including anoxia, drought, submergence, alkali, heat, dehydration, salt, and cold, they showed a differential expression pattern. The expression analysis from the RT-qPCR data also showed expression under various abiotic stress conditions, comprising cold, salinity, cadmium, drought, submergence, and especially heat stress. This finding will pave the way for future in-vivo analysis, unveil the molecular mechanisms of APY genes in stress response, and contribute to the development of stress-tolerant rice varieties.

Studies on the Requirement of Transthyretin Protein (BxTTR-52) for the Suppression of Host Innate Immunity in Bursaphelenchus xylophilus

The pinewood nematode, Bursaphelenchus xylophilus, has been determined as one of the world's top ten plant-parasitic nematodes. It causes pine wilt, a progressive disease that affects the economy and ecologically sustainable development in East Asia. B. xylophilus secretes pathogenic proteins into host plant tissues to promote infection. However, little is known about the interaction between B. xylophilus and pines. Previous studies reported transthyretin proteins in some species and their strong correlation with immune evasion, which has also been poorly studied in B. xylophilus. In this study, we cloned and functionally validated the B. xylophilus pathogenic protein BxTTR-52, containing a transthyretin domain. An in situ hybridization assay demonstrated that BxTTR-52 was expressed mainly in the esophageal glands of B. xylophilus. Confocal microscopy revealed that BxTTR-52-RFP localized to the nucleus, cytoplasm, and plasma membrane. BxTTR-52 recombinant proteins produced by Escherichia coli could be suppressed by hydrogen peroxide and antioxidant enzymes in pines. Moreover, silencing BxTTR-52 significantly attenuated the morbidity of Pinus thunbergii infected with B. xylophilus. It also suppressed the expression of pathogenesis-related genes in P. thunbergii. These results suggest that BxTTR-52 suppresses the plant immune response in the host pines and might contribute to the pathogenicity of B. xylophilus in the early infection stages.

Stem and leaf rust-induced miRNAome in bread wheat near-isogenic lines and their comparative analysis

Wheat rusts remain a major threat to global wheat production and food security. The R-gene-mediated resistance has been employed as an efficient approach to develop rust-resistant varieties. However, evolution of new fungal races and infection strategies put forward the urgency of unravelling novel molecular players, including non-coding RNAs for plant response. This study identified microRNAs associated with Sr36 and Lr45 disease resistance genes in response to stem and leaf rust, respectively. Here, small RNA sequencing was performed on susceptible and resistant wheat near-isogenic lines inoculated with stem and leaf rust pathotypes. microRNA mining in stem rust-inoculated cultivars revealed a total of distinct 26 known and 7 novel miRNAs, and leaf rust libraries culminated with 22 known and 4 novel miRNAs. The comparative analysis between two disease sets provides a better understanding of altered miRNA profiles associated with respective R-genes and infections. Temporal differential expression pattern of miRNAs pinpoints their role during the progress of infection. Differential expression pattern of miRNAs among various treatments as well as time-course expression of miRNAs revealed stem and leaf rust-responsive miRNAs and their possible role in balancing disease resistance/susceptibility. Disclosure of guide strand, passenger strand and a variant of novel-Tae-miR02 from different subgenome origins might serve as a potential link between stem and leaf rust defence mechanisms downstream to respective R-genes. The outcome from the analysis of microRNA dynamics among two rust diseases and further characterization of identified microRNAs can contribute to significant novel insights on wheat-rust interactions and rust management. KEY POINTS: • Identification and comparative analysis of stem and leaf rust-responsive miRNAs. • Chromosomal location and functional prediction of miRNAs. • Time-course expression analysis of pathogen-responsive miRNAs.

Genome-wide identification, comprehensive characterization of transcription factors, cis-regulatory elements, protein homology, and protein interaction network of DREB gene family in Solanum lycopersicum

Tomato is a drought-sensitive crop which has high susceptibility to adverse climatic changes. Dehydration-responsive element-binding (DREB) are significant plant transcription factors that have a vital role in regulating plant abiotic stress tolerance by networking with DRE/CRT cis-regulatory elements in response to stresses. In this study, bioinformatics analysis was performed to conduct the genome-wide identification and characterization of DREB genes and promoter elements in Solanum lycopersicum. In genome-wide coverage, 58 SlDREB genes were discovered on 12 chromosomes that justified the criteria of the presence of AP2 domain as conserved motifs. Intron-exon organization and motif analysis showed consistency with phylogenetic analysis and confirmed the absence of the A3 class, thus dividing the SlDREB genes into five categories. Gene expansion was observed through tandem duplication and segmental duplication gene events in SlDREB genes. Ka/Ks values were calculated in ortholog pairs that indicated divergence time and occurrence of purification selection during the evolutionary period. Synteny analysis demonstrated that 32 out of 58 and 47 out of 58 SlDREB genes were orthologs to Arabidopsis and Solanum tuberosum, respectively. Subcellular localization predicted that SlDREB genes were present in the nucleus and performed primary functions in DNA binding to regulate the transcriptional processes according to gene ontology. Cis-acting regulatory element analysis revealed the presence of 103 motifs in 2.5-kbp upstream promoter sequences of 58 SlDREB genes. Five representative SlDREB proteins were selected from the resultant DREB subgroups for 3D protein modeling through the Phyre2 server. All models confirmed about 90% residues in the favorable region through Ramachandran plot analysis. Moreover, active catalytic sites and occurrence in disorder regions indicated the structural and functional flexibility of SlDREB proteins. Protein association networks through STRING software suggested the potential interactors that belong to different gene families and are involved in regulating similar functional and biological processes. Transcriptome data analysis has revealed that the SlDREB gene family is engaged in defense response against drought and heat stress conditions in tomato. Overall, this comprehensive research reveals the identification and characterization of SlDREB genes that provide potential knowledge for improving abiotic stress tolerance in tomato.

Detection and validation of cis-regulatory motifs in osmotic stress-inducible synthetic gene switches via computational and experimental approaches

Synthetic cis -regulatory modules can improve our understanding of gene regulatory networks. We applied an ensemble approach for de novo cis motif discovery among the promoters of 181 drought inducible differentially expressed soybean (Glycine max L.) genes. A total of 43 cis motifs were identified in promoter regions of all gene sets using the binding site estimation suite of tools (BEST). Comparative analysis of these motifs revealed similarities with known cis -elements found in PLACE database and led to the discovery of cis -regulatory motifs that were not yet implicated in drought response. Compiled with the proposed synthetic promoter design rationale, three synthetic assemblies were constructed by concatenating multiple copies of drought-inducible cis motifs in a specific order with inter-motif spacing using random bases and placed upstream of 35s minimal core promoter. Each synthetic module substituted 35S promoter in pBI121 and pCAMBIA3301 to drive glucuronidase expression in soybean hairy roots and Arabidopsis thaliana L. Chimeric soybean seedlings and 3-week-old transgenic Arabidopsis plants were treated with simulated with different levels of osmotic stress. Histochemical staining of transgenic soybean hairy roots and Arabidopsis displayed drought-inducible GUS activity of synthetic promoters. Fluorometric assay and expression analysis revealed that SP2 is the better manual combination of cis -elements for stress-inducible expression. qRT-PCR results further demonstrated that designed synthetic promoters are not tissue-specific and thus active in different parts upon treatment with osmotic stress in Arabidopsis plants. This study provides tools for transcriptional upgradation of valuable crops against drought stress and adds to the current knowledge of synthetic biology.

Molecular and expression analysis indicate the role of CBL interacting protein kinases (CIPKs) in abiotic stress signaling and development in chickpea

Calcineurin B-like proteins (CBL)-interacting protein kinases (CIPKs) regulate the developmental processes, hormone signal transduction and stress responses in plants. Although the genome sequence of chickpea is available, information related to the CIPK gene family is missing in this important crop plant. Here, a total of 22 CIPK genes were identified and characterized in chickpea. We found a high degree of structural and evolutionary conservation in the chickpea CIPK family. Our analysis showed that chickpea CIPKs have evolved with dicots such as Arabidopsis and soybean, and extensive gene duplication events have played an important role in the evolution and expansion of the CIPK gene family in chickpea. The three-dimensional structure of chickpea CIPKs was described by protein homology modelling. Most CIPK proteins are localized in the cytoplasm and nucleus, as predicted by subcellular localization analysis. Promoter analysis revealed various cis-regulatory elements related to plant development, hormone signaling, and abiotic stresses. RNA-seq expression analysis indicated that CIPKs are significantly expressed through a spectrum of developmental stages, tissue/organs that hinted at their important role in plant development. The qRT-PCR analysis revealed that several CaCIPK genes had specific and overlapping expressions in different abiotic stresses like drought, salt, and ABA, suggesting the important role of this gene family in abiotic stress signaling in chickpea. Thus, this study provides an avenue for detailed functional characterization of the CIPK gene family in chickpea and other legume crops.

Comprehensive In Silico Analysis of RNA Silencing-Related Genes and Their Regulatory Elements in Wheat (Triticum aestivum L.)

Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are known as the three major gene families that act as the critical components of RNA interference or silencing mechanisms through the noncoding small RNA molecules (miRNA and siRNA) to regulate the expressions of protein-coding genes in eukaryotic organisms. However, most of their characteristics including structures, chromosomal location, subcellular locations, regulatory elements, and gene networking were not rigorously studied. Our analysis identified 7 TaDCL, 39 TaAGO, and 16 TaRDR genes as RNA interference (RNAi) genes from the wheat genome. Phylogenetic analysis of predicted RNAi proteins with the RNAi proteins of Arabidopsis and rice showed that the predicted proteins of TaDCL, TaAGO, and TaRDR groups are clustered into four, eight, and four subgroups, respectively. Domain, 3D protein structure, motif, and exon-intron structure analyses showed that these proteins conserve identical characteristics within groups and maintain differences between groups. The nonsynonymous/synonymous mutation ratio (Ka/Ks) < 1 suggested that these protein sequences conserve some purifying functions. RNAi genes networking with TFs revealed that ERF, MIKC-MADS, C2H2, BBR-BPC, MYB, and Dof are the key transcriptional regulators of the predicted RNAi-related genes. The cis-regulatory element (CREs) analysis detected some important CREs of RNAi genes that are significantly associated with light, stress, and hormone responses. Expression analysis based on an online database exhibited that almost all of the predicted RNAi genes are expressed in different tissues and organs. A case-control study from the gene expression level showed that some RNAi genes significantly responded to the drought and heat stresses. Overall results would therefore provide an excellent basis for in-depth molecular investigation of these genes and their regulatory elements for wheat crop improvement against different stressors.

Functional characterization of the pUceS8.3 promoter and its potential use for ectopic gene overexpression

The pUceS8.3 is a constitutive gene promoter with potential for ectopic and strong genes overexpression or active biomolecules in plant tissues attacked by pests, including nematode-induced giant cells or galls. Soybean (Glycine max) is one of the most important agricultural commodities worldwide and a major protein and oil source. Herein, we identified the soybean ubiquitin-conjugating (E2) enzyme gene (GmUBC4; Glyma.18G216000), which is significantly upregulated in response to Anticarsia gemmatalis attack and Meloidogyne incognita-induced galls during plant parasitism by plant nematode. The GmUBC4 promoter sequence and its different modules were functionally characterized in silico and in planta using transgenic Arabidopsis thaliana and G. max lines. Its full-length transcriptional regulatory region (promoter and 5´-UTR sequences, named pUceS8.3 promoter) was able to drive higher levels of uidA (β-glucuronidase) gene expression in different tissues of transgenic A. thaliana lines compared to its three shortened modules and the p35SdAMV promoter. Notably, higher β-glucuronidase (GUS) enzymatic activity was shown in M. incognita-induced giant cells when the full pUceS8.3 promoter drove the expression of this reporter gene. Furthermore, nematode-specific dsRNA molecules were successfully overexpressed under the control of the pUceS8.3 promoter in transgenic soybean lines. The RNAi gene construct used here was designed to post-transcriptionally downregulate the previously characterized pre-mRNA splicing factor genes from Heterodera glycines and M. incognita. A total of six transgenic soybean lines containing RNAi gene construct were selected for molecular characterization after infection with M. incognita pre-parasitic second-stage (ppJ2) nematodes. A strong reduction in the egg number produced by M. incognita after parasitism was observed in those transgenic soybean lines, ranging from 71 to 92% compared to wild-type control plants. The present data demonstrated that pUceS8.3 is a gene promoter capable of effectively driving dsRNA overexpression in nematode-induced giant cells of transgenic soybean lines and can be successfully applied as an important biotechnological asset to generate transgenic crops with improved resistance to root-knot nematodes as well as other pests.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression analysis under Fusarium kyushuense and drought stress conditions in Passiflora edulis

The NAC gene family is one of the largest plant transcription factors (TFs) families and plays important roles in plant growth, development, metabolism, and biotic and abiotic stresses. However, NAC gene family has not been reported in passion fruit (Passiflora edulis). In this study, a total of 105 NAC genes were identified in the passion fruit genome and were unevenly distributed across all nine-passion fruit chromomere, with a maximum of 48 PeNAC genes on chromosome one. The physicochemical features of all 105 PeNAC genes varied including 120 to 3,052 amino acids, 3 to 8 conserved motifs, and 1 to 3 introns. The PeNAC genes were named (PeNAC001-PeNAC105) according to their chromosomal locations and phylogenetically grouped into 15 clades (NAC-a to NAC-o). Most PeNAC proteins were predicted to be localized in the nucleus. The cis-element analysis indicated the possible roles of PeNAC genes in plant growth, development, light, hormones, and stress responsiveness. Moreover, the PeNAC gene duplications including tandem (11 gene pairs) and segmental (12 gene pairs) were identified and subjected to purifying selection. All PeNAC proteins exhibited similar 3D structures, and a protein-protein interaction network analysis with known Arabidopsis proteins was predicted. Furthermore, 17 putative ped-miRNAs were identified to target 25 PeNAC genes. Potential TFs including ERF, BBR-BPC, Dof, and bZIP were identified in promoter region of all 105 PeNAC genes and visualized in a TF regulatory network. GO and KEGG annotation analysis exposed that PeNAC genes were related to different biological, molecular, and cellular terms. The qRT-PCR expression analysis discovered that most of the PeNAC genes including PeNAC001, PeNAC003, PeNAC008, PeNAC028, PeNAC033, PeNAC058, PeNAC063, and PeNAC077 were significantly upregulated under Fusarium kyushuense and drought stress conditions compared to controls. In conclusion, these findings lay the foundation for further functional studies of PeNAC genes to facilitate the genetic improvement of plants to stress resistance.

Isolation, cloning and functional analysis of a putative constitutive promoter of E3 ubiquitin- protein ligase RF4 from Coleus amboinicus Lour

Promoter is a region in the genome sequence located upstream of the transcription start site comprising cis acting elements, which initiates and regulates the transcription of an associated gene. As the need for genetically engineered plants has widened, the requirement to develop methods to optimize the control of transgene expression has also increased. Therefore, analyzing the functionality of the promoter is very important in understanding the target gene expression. The widespread use of viral constitutive promoters (Cauliflower mosaic virus - CaMV35) has raised concerns about the safety and containment of the transgene in the environment. Hence isolation and characterization of novel promoters using fast and efficient genetic engineering tools is the need of the hour. The present study, for the first time, describes the isolation and characterization of a novel constitutive promoter driving Ubiquitin E3 ligase from the plant Coleus amboinicus, a perennial herb, of Lamiaceae family. The functionality of the isolated promoter was demonstrated using the β Glucuronidase as a reporter in tobacco var Petit havana. Development of blue color in the tobacco leaves indicated the presence of a functional promoter. We describe for the first time the isolation and characterization of E3 ubiquitin- protein ligase RF4 promoter from Coleus amboinicus Lour. In silico analysis revealed the presence of core promoter elements and other responsive elements in the promoter. The functionality of the promoter was demonstrated in tobacco leaf discs via GUS staining. This article is protected by copyright. All rights reserved.

In-Silico Study of Brassinosteroid Signaling Genes in Rice Provides Insight Into Mechanisms Which Regulate Their Expression

Brassinosteroids (BRs) regulate a diverse spectrum of processes during plant growth and development and modulate plant physiology in response to environmental fluctuations and stress factors. Thus, the BR signaling regulators have the potential to be targeted for gene editing to optimize the architecture of plants and make them more resilient to environmental stress. Our understanding of the BR signaling mechanism in monocot crop species is limited compared to our knowledge of this process accumulated in the model dicot species - Arabidopsis thaliana. A deeper understanding of the BR signaling and response during plant growth and adaptation to continually changing environmental conditions will provide insight into mechanisms that govern the coordinated expression of the BR signaling genes in rice (Oryza sativa) which is a model for cereal crops. Therefore, in this study a comprehensive and detailed in silico analysis of promoter sequences of rice BR signaling genes was performed. Moreover, expression profiles of these genes during various developmental stages and reactions to several stress conditions were analyzed. Additionally, a model of interactions between the encoded proteins was also established. The obtained results revealed that promoters of the 39 BR signaling genes are involved in various regulatory mechanisms and interdependent processes that influence growth, development, and stress response in rice. Different transcription factor-binding sites and cis-regulatory elements in the gene promoters were identified which are involved in regulation of the genes’ expression during plant development and reactions to stress conditions. The in-silico analysis of BR signaling genes in O. sativa provides information about mechanisms which regulate the coordinated expression of these genes during rice development and in response to other phytohormones and environmental factors. Since rice is both an important crop and the model species for other cereals, this information may be important for understanding the regulatory mechanisms that modulate the BR signaling in monocot species. It can also provide new ways for the plant genetic engineering technology by providing novel potential targets, either cis-elements or transcriptional factors, to create elite genotypes with desirable traits.

Genome-wide analysis of the CAD gene family reveals two bona fide CAD genes in oil palm

Cinnamyl alcohol dehydrogenase (CAD) is the key enzyme for lignin biosynthesis in plants. In this study, genome-wide analysis was performed to identify CAD genes in oil palm (Elaeis guineensis). Phylogenetic analysis was then conducted to select the bona fide EgCADs. The bona fide EgCAD genes and their respective 5' flanking regions were cloned and analysed. Their expression profiles were evaluated in various organs using RT-PCR. Seven EgCAD genes (EgCAD1-7) were identified and divided into four phylogenetic groups. EgCAD1 and EgCAD2 display high sequence similarities with other bona fide CADs and possess all the signature motifs of the bona fide CAD. They also display similar 3D protein structures. Gene expression analysis showed that EgCAD1 was expressed most abundantly in the root tissues, while EgCAD2 was expressed constitutively in all the tissues studied. EgCAD1 possesses only one transcription start site, while EgCAD2 has five. Interestingly, a TC microsatellite was found in the 5' flanking region of EgCAD2. The 5' flanking regions of EgCAD1 and EgCAD2 contain lignin-associated regulatory elements i.e. AC-elements, and other defence-related motifs, including W-box, GT-1 motif and CGTCA-motif. Altogether, these results imply that EgCAD1 and EgCAD2 are bona fide CAD involved in lignin biosynthesis during the normal development of oil palm and in response to stresses. Our findings shed some light on the roles of the bona fide CAD genes in oil palm and pave the way for manipulating lignin content in oil palm through a genetic approach.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03208-0.

Comprehensive Genome-Wide Identification and Expression Profiling of Eceriferum (CER) Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Plant surfaces are covered with cuticle wax and are the first barrier between a plant and environmental stresses. Eceriferum (CER) is an important gene family involved in wax biosynthesis and stress resistance. In this study, for the first time, 34 CER genes were identified in the passion fruit (Passiflora edulis) genome, and PeCER proteins varied in physicochemical properties. A phylogenetic tree was constructed and divided into seven clades to identify the evolutionary relationship with other plant species. Gene structure analyses revealed that conserved motifs ranged from 1 to 24, and that exons ranged from 1 to 29. The cis-element analysis provides insight into possible roles of PeCER genes in plant growth, development and stress responses. The syntenic analysis revealed that segmental (six gene pairs) and tandem (six gene pairs) gene duplication played an important role in the expansion of PeCER genes and underwent a strong purifying selection. In addition, 12 putative ped-miRNAs were identified to be targeting 16 PeCER genes, and PeCER6 was the most targeted by four miRNAs including ped-miR157a-5p, ped-miR164b-5p, ped-miR319b, and ped-miR319l. Potential transcription factors (TFs) such as ERF, AP2, MYB, and bZIP were predicted and visualized in a TF regulatory network interacting with PeCER genes. GO and KEGG annotation analysis revealed that PeCER genes were highly related to fatty acid, cutin, and wax biosynthesis, plant-pathogen interactions, and stress response pathways. The hypothesis that most PeCER proteins were predicted to localize to the plasma membrane was validated by transient expression assays of PeCER32 protein in onion epidermal cells. qRT-PCR expression results showed that most of the PeCER genes including PeCER1, PeCER11, PeCER15, PeCER17, and PeCER32 were upregulated under drought and Fusarium kyushuense stress conditions compared to controls. These findings provide a foundation for further studies on functions of PeCER genes to further facilitate the genetic modification of passion fruit wax biosynthesis and stress resistance.

Uncovering the role of DNA methyltransferases in grapevine - Plasmopara viticola interaction: From genome-wide characterization to global methylation patterns

Epigenetic regulation has recently gained prominence in the field of plant-pathogen interactions, providing a deeper understanding of the molecular mechanisms associated with plant infection. In grapevine interaction with pathogens, epigenetic regulation still remains a black box. In this work, we characterized grapevine DNA methyltransferase gene family and identified nine DNA methyltransferases genes across eight grapevine chromosomes coding for 17 proteins. We also assessed the modulation of global cytosine methylation and gene expression levels of these genes with the aim of establishing a connection between DNA methylation and grapevine resistance towards downy mildew. Our results revealed that, in the incompatible interaction, an early hypomethylation, coupled with downregulation of DNMT and CMT genes occurs very early after pathogen inoculation. Additionally, the compatible interaction is characterized by a hypermethylation at 6hpi. A temporal delay is evident between the shifts in DNA methyltransferases gene expression in both compatible and incompatible interactions which in turn may be reflected in the global methylation percentage. Overall, we present the first evidence of an epigenetic regulation role in grapevine defense against P. viticola.