Evolution of intron/exon structure of DEAD helicase family genes in Arabidopsis, Caenorhabditis, and Drosophila

The RNA helicase LOS4 regulates pre-mRNA splicing of key genes (EIN2, ERS2, CTR1) in the ethylene signaling pathway

KEY MESSAGE

The Arabidopsis RNA helicase LOS4 plays a key role in regulating pre-mRNA splicing of the genes EIN2, CTR1, and ERS2 in ethylene signaling pathway. The plant hormone ethylene plays diverse roles in plant growth, development, and responses to stress. Ethylene is perceived by the membrane-bound ethylene receptors complex, and then triggers downstream components, such as EIN2, to initiate signal transduction into the nucleus, leading to the activation of ethylene-responsive genes. Over the past decades, substantial information has been accumulated regarding gene cloning, protein-protein interactions, and downstream gene expressions in the ethylene pathway. However, our understanding of mRNA post-transcriptional processing and modification of key genes in the ethylene signaling pathway remains limited. This study aims to provide evidence demonstrating the involvement of the Arabidopsis RNA helicase LOS4 in pre-mRNA splicing of the genes EIN2, CTR1, and ERS2 in ethylene signaling pathway. Various genetic approaches including RNAi gene silencing, CRISPR-Cas9 gene editing, and amino acid mutations were employed in this study. When LOS4 was silenced or knocked down, the ethylene sensitivity of etiolated seedlings was significantly enhanced. Further investigation revealed errors in the EIN2 pre-mRNA splicing when LOS4 was knocked down. In addition, aberrant pre-mRNA splicing was observed in the ERS2 and CTR1 genes in the pathway. Biochemical assays indicated that the los4-2 (E94K) mutant protein exhibited increased ATP binding and enhanced ATP hydrolytic activity. Conversely, the los4-1 (G364R) mutant had reduced substrate RNA binding and lower ATP binding activities. These findings significantly advanced our comprehension of the regulatory functions and molecular mechanisms of RNA helicase in ethylene signaling.

DDX RNA helicases: key players in cellular homeostasis and innate antiviral immunity

RNA helicases are integral in RNA metabolism, performing important roles in cellular homeostasis and stress responses. In particular, the DExD/H-box (DDX) helicase family possesses a conserved catalytic core that binds structural features rather than specific sequences in RNA targets. DDXs have critical roles in all aspects of RNA metabolism including ribosome biogenesis, translation, RNA export, and RNA stability. Importantly, functional specialization within this family arises from divergent N and C termini and is driven at least in part by gene duplications with 18 of the 42 human helicases having paralogs. In addition to their key roles in the homeostatic control of cellular RNA, these factors have critical roles in RNA virus infection. The canonical RIG-I-like receptors (RLRs) play pivotal roles in cytoplasmic sensing of viral RNA structures, inducing antiviral gene expression. Additional RNA helicases function as viral sensors or regulators, further diversifying the innate immune defense arsenal. Moreover, some of these helicases have been coopted by viruses to facilitate their replication. Altogether, DDX helicases exhibit functional specificity, playing intricate roles in RNA metabolism and host defense. This review will discuss the mechanisms by which these RNA helicases recognize diverse RNA structures in cellular and viral RNAs, and how this impacts RNA processing and innate immune responses.

Mining key genes associated with phosphorus deficiency through genome-wide identification and characterization of cucumber SPX family genes

BACKGROUND

Proteins harboring the SPX domain are crucial for the regulation of phosphate (Pi) homeostasis in plants. This study aimed to identify and analyze the entire SPX gene family within the cucumber genome.

RESULTS

The cucumber genome encompassed 16 SPX domain-containing genes, which were distributed across six chromosomes and categorized into four distinct subfamilies: SPX, SPX-MFS, SPX-EXS and SPX-RING, based on their structure characteristics. Additionally, gene duplications and synteny analysis were conducted for CsSPXs, revealing that their promoter regions were enriched with a variety of hormone-responsive, biotic/abiotic stress and typical P1BS-related elements. Tissue expression profiling of CsSPX genes revealed that certain members were specifically expressed in particular organs, suggesting essential roles in cucumber growth and development. Under low Pi stress, CsSPX1 and CsSPX2 exhibited a particularly strong response to Pi starvation. It was observed that the cucumber cultivar Xintaimici displayed greater tolerance to low Pi compared to black-spined cucumber under low Pi stress conditions. Protein interaction networks for the 16 CsSPX proteins were predicted, and yeast two-hybrid assay revealed that CsPHR1 interacted with CsSPX2, CsSPX3, CsSPX4 and CsSPX5, implying their involvement in the Pi signaling pathway in conjunction with CsPHR1.

CONCLUSION

This research lays the foundation for further exploration of the function of the CsSPX genes in response to low Pi stress and for elucidating the underlying mechanism.

Crape myrtle LiGAoxs displaying activities of gibberellin oxidases respond to branching architecture

In the realm of ornamental horticulture, crape myrtle (Lagerstroemia indica) stands out for its aesthetic appeal, attributed largely to its vibrant flowers and distinctive branching architecture. This study embarked on a comprehensive exploration of the gibberellin oxidase (GAox) gene family in crape myrtle, illuminating its pivotal role in regulating GA levels, a key determinant of plant developmental processes. We identified and characterized 36 LiGAox genes, subdivided into GA2ox, GA3ox, GA20ox, and GAox-like subgroups, through genomic analyses. These genes' evolutionary trajectories were delineated, revealing significant gene expansions attributed to segmental duplication events. Functional analyses highlighted the divergent expression patterns of LiGAox genes across different crape myrtle varieties, associating them with variations in flower color and branching architecture. Enzymatic activity assays on selected LiGA2ox enzymes exhibited pronounced GA2 oxidase activity, suggesting a potential regulatory role in GA biosynthesis. Our findings offered a novel insight into the molecular underpinnings of GA-mediated growth and development in L. indica, providing a foundational framework for future genetic enhancements aimed at optimizing ornamental traits.

Evolution of the SPX gene family and its role in the response mechanism to low phosphorus stress in self-rooted apple stock

BACKGROUND

Phosphorus plays a key role in plant adaptation to adversity and plays a positive role in the yield and quality formation of apples. Genes of the SPX domain-containing family are widely involved in the regulation of phosphorus signalling networks. However, the mechanisms controlling phosphorus deficiency are not completely understood in self-rooted apple stock.

RESULTS

In this study, 26 members of the apple SPX gene family were identified by genome-wide analysis, and further divided into four subfamilies (SPX, SPX-MFS, SPX-EXS, and SPX-RING) based on their structural features. The chromosome distribution and gene duplications of MdSPXs were also examined. The promoter regions of MdSPXs were enriched for multiple biotic/abiotic stresses, hormone responses and typical P1BS-related elements. Analysis of the expression levels of 26 MdSPXs showed that some members were remarkably induced when subjected to low phosphate (Pi) stress, and in particular MdSPX2, MdSPX3, and MdPHO1.5 exhibited an intense response to low Pi stress. MdSPX2 and MdSPX3 showed significantly divergent expression levels in low Pi sensitive and insensitive apple species. Protein interaction networks were predicted for 26 MdSPX proteins. The interaction of MdPHR1 with MdSPX2, MdSPX3, MdSPX4, and MdSPX6 was demonstrated by yeast two-hybrid assay, suggesting that these proteins might be involved in the Pi-signaling pathway by interacting with MdPHR1.

CONCLUSION

This research improved the understanding of the apple SPX gene family and contribute to future biological studies of MdSPX genes in self-rooted apple stock.

Characterization, expression pattern, and function analysis of gibberellin oxidases in Salix matsudana

Gibberellin oxidases (GAoxs) identified from many species play indispensable roles in GA biosynthesis and GA signal transduction. However, there has been limited research conducted on the GAox family of Salix matsudana, a tetraploid ornamental tree species. Here, 54 GAox genes were identified from S. matsudana and renamed as SmGA20ox1-22, SmGA2ox1-24, SmGA3ox1-6, and SmGAox-like1/2. Gene structure and conserved motif analysis showed that SmGA3ox members possess the 1 intron and other SmGAoxs contain 2-3 introns, and motif 1/2/7 universally present in all SmGAoxs. A total of 69 gene pairs were identified from SmGAox family members, and the Ka/Ks values indicated the SmGAoxs experience the purifying selection. The intra species collinearity analysis implied S. matsudana, S. purpurea, and Populus trichocarpa have the close genetic relationship. The GO analysis suggested SmGAoxs are dominantly involved in GA metabolic process, ion binding, and oxidoreductase activity. RNA-sequencing demonstrated that some SmGAoxs may play an essential role in salt and submergence stresses. In addition, the SmGA20ox13/21 displayed the dominant vitality of GA20 oxidase, but the SmGA20ox13/21 still possessed low activities of GA2 and GA3 oxidases. This study can contribute to reveal the regulatory mechanism of salt and submergence tolerance in willow.

In-Depth Characterization of bZIP Genes in the Context of Endoplasmic Reticulum (ER) Stress in Brassica campestris ssp. chinensis

Numerous studies have been conducted to investigate the genomic characterization of bZIP genes and their involvement in the cellular response to endoplasmic reticulum (ER) stress. These studies have provided valuable insights into the coordinated cellular response to ER stress, which is mediated by bZIP transcription factors (TFs). However, a comprehensive and systematic investigations regarding the role of bZIP genes and their involvement in ER stress response in pak choi is currently lacking in the existing literature. To address this knowledge gap, the current study was initiated to elucidate the genomic characteristics of bZIP genes, gain insight into their expression patterns during ER stress in pak choi, and investigate the protein-to-protein interaction of bZIP genes with the ER chaperone BiP. In total, 112 members of the BcbZIP genes were identified through a comprehensive genome-wide analysis. Based on an analysis of sequence similarity, gene structure, conserved domains, and responsive motifs, the identified BcbZIP genes were categorized into 10 distinct subfamilies through phylogenetic analysis. Chromosomal location and duplication events provided insight into their genomic context and evolutionary history. Divergence analysis estimated their evolutionary history with a predicted divergence time ranging from 0.73 to 80.71 million years ago (MYA). Promoter regions of the BcbZIP genes were discovered to exhibit a wide variety of cis-elements, including light, hormone, and stress-responsive elements. GO enrichment analysis further confirmed their roles in the ER unfolded protein response (UPR), while co-expression network analysis showed a strong relationship of BcbZIP genes with ER-stress-responsive genes. Moreover, gene expression profiles and protein-protein interaction with ER chaperone BiP further confirmed their roles and capacity to respond to ER stress in pak choi.

Genome-wide systematic survey and analysis of the RNA helicase gene family and their response to abiotic stress in sweetpotato

Sweetpotato (Ipomoea batatas (L.) Lam.) holds a crucial position as one of the staple foods globally, however, its yields are frequently impacted by environmental stresses. In the realm of plant evolution and the response to abiotic stress, the RNA helicase family assumes a significant role. Despite this importance, a comprehensive understanding of the RNA helicase gene family in sweetpotato has been lacking. Therefore, we conducted a comprehensive genome-wide analysis of the sweetpotato RNA helicase family, encompassing aspects such as chromosome distribution, promoter elements, and motif compositions. This study aims to shed light on the intricate mechanisms underlying the stress responses and evolutionary adaptations in sweetpotato, thereby facilitating the development of strategies for enhancing its resilience and productivity. 300 RNA helicase genes were identified in sweetpotato and categorized into three subfamilies, namely IbDEAD, IbDEAH and IbDExDH. The collinearity relationship between the sweetpotato RNA helicase gene and 8 related homologous genes from other species was explored, providing a reliable foundation for further study of the sweetpotato RNA helicase gene family's evolution. Furthermore, through RNA-Seq analysis and qRT-PCR verification, it was observed that the expression of eight RNA helicase genes exhibited significant responsiveness to four abiotic stresses (cold, drought, heat, and salt) across various tissues of ten different sweetpotato varieties. Sweetpotato transgenic lines overexpressing the RNA helicase gene IbDExDH96 were generated using A.rhizogenes-mediated technology. This approach allowed for the preliminary investigation of the role of sweetpotato RNA helicase genes in the response to cold stress. Notably, the promoters of RNA helicase genes contained numerous cis-acting elements associated with temperature, hormone, and light response, highlighting their crucial role in sweetpotato abiotic stress response.

Analysis of the HD-Zip I transcription factor family in Salvia miltiorrhiza and functional research of SmHD-Zip12 in tanshinone synthesis

Background

The homeodomain-leucine zipper I (HD-Zip I) transcription factor is a plant-specific protein that plays an essential role in the abiotic stress response of plants. Research on the HD-Zip I family in Salvia miltiorrhiza is still lacking.

Methods and Results

In this study, a total of 25 SmHD-Zip I proteins were identified. Their characterizations, phylogenetic relationships, conserved motifs, gene structures, and cis-elements were analyzed comprehensively using bioinformatics methods. Expression profiling revealed that SmHD-Zip I genes exhibited distinctive tissue-specific patterns and divergent responses to ABA, PEG, and NaCl stresses. SmHD-Zip12 responded the most strongly to ABA, PEG, and NaCl, so it was used for transgenic experiments. The overexpression of SmHD-Zip12 significantly increased the content of cryptotanshinone, dihydrotanshinone I, tanshinone I, and tanshinone IIA by 2.89-fold, 1.85-fold, 2.14-fold, and 8.91-fold compared to the wild type, respectively. Moreover, in the tanshinone biosynthetic pathways, the overexpression of SmHD-Zip12 up-regulated the expression levels of SmAACT, SmDXS, SmIDS, SmGGPPS, SmCPS1, SmCPS2, SmCYP76AH1, SmCYP76AH3, and SmCYP76AK1 compared with the wild type.

Conclusions

This study provides information the possible functions of the HD-Zip I family and lays a theoretical foundation for clarifying the functional mechanism of the SmHD-Zip12 gene in regulating the synthesis of tanshinone in S. miltiorrhiza.

Genome-wide identification glutathione-S-transferase gene superfamily in Daphnia pulex and its transcriptional response to nanoplastics

Glutathione S-transferases (GSTs) are key multifunctional phase II detoxification enzymes involved in the regulation of growth, development, and stress responses. However, the knowledge of GSTs in the model invertebrate organism Daphnia pulex at the genomic level remains limited. In the present study, 35 GST genes were identified in D. pulex (Dp-GST), belonging to eight subfamilies, with the sigma, mu, and delta/epsilon subfamilies constituting approximately 29 %, 20 %, and 20 % of the GST superfamily, respectively. Chromosome tandem duplication of genes within the same subfamily was observed, which may be the main force driving GST expansion in D. pulex. DpGST genes showed different expression patterns in response to nanoplastic exposure for 96 h and 21 days. Some homologous GST genes in D. pulex showed similar expression patterns in response to nanoplastic exposure, likely owing to their unique motifs. For example, motif 9 is found in all delta/epsilon GST genes, whereas motifs 1, 2, 3, 5, and 7 are highly conserved in sigma GST genes. The characterization of D. pulex GSTs extending the knowledge of GST-mediated environmental contaminants, especially nanoplastics.

Comprehensive analysis of cucumber RAV family genes and functional characterization of CsRAV1 in salt and ABA tolerance in cucumber

The RAV (related to ABI3 and VP1) transcription factors are specific and exist in plants, which contain a B3 DNA binding domain and/or an APETALA2 (AP2) DNA binding domain. RAVs have been extensively studied in plants, and more and more evidences show that RAVs are involved in various aspects of plant growth and development, stress resistance and hormone signal transduction. However, the systematic analysis of RAV family in cucumber is rarely reported. In this study, eight CsRAV genes were identified in cucumber genome and we further comprehensively analyzed their protein physicochemical properties, conserved domains, gene structure and phylogenetic relationships. The synteny analysis and gene duplications of CsRAV genes were also analysed. Cis-element analysis revealed that the CsRAVs promoter contained several elements related to plant hormones and abiotic stress. Expression analysis showed that NaCl and ABA could significantly induce CsRAV genes expression. Subcellular localization revealed that all CsRAVs were localized in the nucleus. In addition, 35S:CsRAV1 transgenic Arabidopsis and cucumber seedlings enhanced NaCl and ABA tolerance, revealing CsRAV1 may be an important regulator of abiotic stress response. In conclusion, comprehensive analysis of CsRAVs would provide certain reference for understanding the evolution and function of the CsRAV genes.

Genome-wide characterization of C2H2 zinc-finger gene family provides insight into the mechanisms and evolution of the dehydration-rehydration responses in Physcomitrium and Arabidopsis

Dehydration tolerance is a vital factor for land plant evolution and world agricultural production. Numerous studies enlightened that the plant-specific C2H2-type zinc-finger proteins (C2H2-ZFPs) as master regulators played pivotal roles in the abiotic stress responses of plants. However, a comprehensive understanding of the evolution of C2H2-ZFPs in terrestrial plants and its regulatory mechanism in dehydration and rehydration response remains a mystery. In this study, the genome-wide identification of C2H2-ZFP genes revealed 549 homologs in the representatives of terrestrial plant lineages from liverwort to angiosperms. Based on the characteristics of the conserved C2H2-ZF domains, four major C2H2-ZF types (M-, Z-, Q-, and D-type) were identified in the C2H2-ZFPs, with the dominants of M-type in all selected species and followed by Z-type in non-seed plants and Q-type in seed plants, respectively. Phylogenetic analyses of the identified C2H2-ZFPs supported four major groups in the land plant representatives, among which the members from the desiccation-tolerant Physcomitrium patens and the dehydration-sensitive Arabidopsis thaliana displayed different topological relationships in the phylogenies reconstructed for a single species. C2H2-ZFPs clustered in the same subclades shared similar features in their conserved domains and gene structures. Approximately, 81% of the C2H2-ZFP promoters of all 549 identified C2H2-ZFPs harbored the conserved ABA-responsive elements (ABREs) and/or dehydration-responsive elements (DREs). Comparative transcriptomic analyses showed that 50 PpZFPs and 56 AtZFPs significantly changed their transcripts abundance. Interestingly, most of the dehydration- and rehydration-responsive PpZPFs and AtZFPs had been predicted to contain the ABRE and DRE elements in their promoter regions and with over half of which phylogenetically belonging to group III. The differences in the expression patterns of C2H2-ZFPs in responses to dehydration and rehydration between P. patens and A. thaliana reflected their different strategies to adapt to dehydration. The identified candidate PpZFPs were specifically induced by moderate dehydration and reached the peak transcript abundance in severe dehydration. Our study lays the foundations for further functional investigation of C2H2-ZFPs in dehydration responses from an evolutionary perspective in land plants. The findings will provide us with genetic resources and potential targets for drought tolerance breeding in crops and beyond.

Genome-wide comprehensive characterization and expression analysis of TLP gene family revealed its responses to hormonal and abiotic stresses in watermelon (Citrullus lanatus)

Thaumatin-like protein (TLP) is the well-known sweetest protein which plays a crucial role in diverse developmental processes and different stress conditions in plants, fungi and animals. The TLP gene family is extensively studied in different plant species including crop plants. Watermelon (Citrullus lanatus) is an important cucurbit crop cultivated worldwide; however, the comprehensive information about the TLP gene family is not available in watermelon. In the present study, we identified the 29 TLP genes as gene family members in watermelon using various computational methods to understand its role in different developmental processes and stress conditions. ClaTLP gene family members were not uniformly distributed on 22 chromosomes. Phylogenetic analysis revealed that the ClaTLP gene family members were grouped into 10 sub-groups. Further, gene duplication analysis showed thirteen gene duplication events which included one tandem and twelve segmental duplications. Amino acid sequence alignment has shown that ClaTLP proteins shared 16 conserved cysteine residues in their THN domain. Furthermore, cis-acting regulatory elements analysis also displayed that ClaTLP gene family members contain diverse phytohormone, various defense, and stress-responsive elements in their promoter region. The expression profile of the ClaTLP gene family revealed the differential expression of gene family members in different tissues and abiotic stresses conditions. Moreover, the expression profile of ClaTLP genes was further validated by semi-quantitative reverse transcriptase PCR. Taken together, these results indicate that ClaTLP genes might play an important role in developmental processes and diverse stress conditions. Therefore, the outcome of this study brings forth the valuable information for further interpret the precise role of ClaTLP gene family members in watermelon.

Genome-wide screening of lectin putative genes from Sorghum bicolor L., distribution in QTLs and a probable implications of lectins in abiotic stress tolerance

BACKGROUND

Sorghum bicolor is one of the most important crops worldwide with the potential to provide resilience when other economic staples might fail against the continuous environmental changes. Many physiological, developmental and tolerance traits in plants are either controlled or influenced by lectins; carbohydrate binding proteins. Hence, we aimed at providing a comprehensive in silico account on sorghum's lectins and study their possible implication on various desired agronomical traits.

RESULTS

We have searched sorghum's genome from grain and sweet types for lectins putative genes that encode proteins with domains capable of differentially binding carbohydrate moieties and trigger various physiological responses. Of the 12 known plant lectin families, 8 were identified regarding their domain architectures, evolutionary relationships, physiochemical characteristics, and gene expansion mechanisms, and they were thoroughly addressed. Variations between grain and sweet sorghum lectin homologs in term of the presence/absence of certain other joint domains like dirigent and nucleotide-binding adaptor shared by APAF-1, R-proteins, and CED-4 (NB-ARC) indicate a possible neofunctionalization. Lectin sequences were found to be preferentially overrepresented in certain quantitative trait loci (QTLs) related to various traits under several subcategories such as cold, drought, salinity, panicle/grain composition, and leaf morphology. The co-localization and distribution of lectins among multiple QTLs provide insights into the pleiotropic effects that could be played by one lectin gene in numerous traits.

CONCLUSION

Our study offers a first-time inclusive details on sorghum lectins and their possible role in conferring tolerance against abiotic stresses and other economically important traits that can be informative for future functional analysis and breeding studies.

Genome-wide identification, in silico characterization and expression analysis of the RNA helicase gene family in chickpea (C. arietinum L.)

The RNA helicases are an important class of enzymes which are known to influence almost every aspect of RNA metabolism. The majority of RNA helicases belong to the SF2 (superfamily 2) superfamily, members of which are further categorized into three separate subfamilies i.e., the DEAD, DEAH and DExD/H-box subfamilies. In chickpea, these RNA helicases have not been characterized until now. A genome-wide analysis across the chickpea genome led to the identification of a total of 150 RNA helicase genes which included 50 DEAD, 33 DEAH and 67 DExD/H-box genes. These were distributed across all the eight chromosomes, with highest number on chromosome 4 (26) and least on chromosome 8 (8). Gene duplication analysis resulted in identification of 15 paralogous gene pairs with Ka/Ks values < 1, indicating towards the genes being under purifying selection during the course of evolution. The promoter regions of the RNA helicase genes were enriched in cis-acting elements like the light and ABA-responsive elements. The drought responsiveness of the genes was analysed by studying the expression profiles of few of these genes, in two different genotypes, the cultivated variety ICC 8261 (kabuli, C. arietinum) and the wild accession ILWC 292 (C. reticulatum), through qRT-PCR. These genotypes were selected based on their drought responsiveness in a field experiment, where it was observed that the percentage (%) reduction in relative water content (RWC) and membrane stability index (MSI) for the drought stressed plants after withholding water for 24 days, over the control or well-watered plants, was least for both the genotypes. The genes CaDEAD50 and CaDExD/H66 were identified as drought-responsive RNA helicase genes in chickpea. The protein encoded by the CaDExD/H66 gene shares a high degree of homology with one of the CLSY (CLASSY) proteins of A. thaliana. We hypothesize that this gene could possibly be involved in regulation of DNA methylation levels in chickpea by regulating siRNA production, in conjunction with other proteins like the Argonaute, RNA dependent RNA polymerases and Dicer-like proteins.

Identification and Expression Analysis of Calcium-Dependent Protein Kinases Gene Family in Potato Under Drought Stress

Calcium-dependent protein kinases (CDPKs) are a class of serine/threonine protein kinases encoded by several gene families that play key roles in stress response and plant growth and development. In this study, the BLAST method was used to search for protein sequences of the potato Calcium-dependent protein kinase gene family. The chromosome location, phylogeny, gene structures, gene duplication, cis-acting elements, protein-protein interaction, and expression profiles were analyzed. Twenty-five CDPK genes in the potato genome were identified based on RNA-seq data and were clustered into four groups (I-IV) based on their structural features and phylogenetic analysis. The result showed the composition of the promoter region of the StCDPKs gene, including light-responsive elements such as Box4, hormone-responsive elements such as ABRE, and stress-responsive elements such as MBS. Four pairs of segmental duplications were found in StCDPKs genes and the Ka/Ks ratios were below 1, indicating a purifying selection of the genes. The protein-protein interaction network revealed defense-related proteins such as; respiratory burst oxidase homologs (RBOHs) interacting with potato CDPKs. Transcript abundance was measured via RT-PCR between the two cultivars and their relative expression of CDPK genes was analyzed after 15, 20, and 25 days of drought. There were varied expression patterns of StCDPK3/13/21 and 23, between the two potato cultivars under mannitol induced-drought conditions. Correlation analysis showed that StCDPK21/22 and StCDPK3 may be the major differentially expressed genes involved in the regulation of malondialdehyde (MDA) and proline content in response to drought stress, opening a new research direction for genetic improvement of drought resistance in potato.

Genome-wide identification and characterization of the Prunus persica ferredoxin gene family and its role in improving heat tolerance

Ferredoxin is involved in many biological processes, such as carbon fixation, nitrogen assimilation, chlorophyll metabolism, and fatty acid synthesis, and it plays a role in plant resistance to stress. However, the functions of Fds in peach during stress are unclear. In this study, 11 members of the peach Fd gene family were identified and divided into six groups (I- VI). We carried out bioinformatics analysis on these sequences, analyzed the physical and chemical properties of PpFd protein and the cis-elements in its promoter region, and predicted and compared the differences in gene structure and conserved protein motifs among groups. The results showed that the PpFd protein was highly conserved in plant species. In addition, overexpression of PpFd08 significantly increased the tolerance of transgenic tomato to high-temperature stress. The transcriptome analysis and qRT-PCR results of PpFd08 transgenic apple calli showed that PpFd08 might enhance heat resistance by modulating the expression of heat tolerance related genes. The results of this study provide a new understanding for the further study of the function of PpFd protein in peach and a candidate gene for improving the heat resistance of peach.

Insight into the CBL and CIPK gene families in pecan (Carya illinoinensis): identification, evolution and expression patterns in drought response

BACKGROUND

Calcium (Ca²⁺) serves as a ubiquitous second messenger and plays a pivotal role in signal transduction. Calcineurin B-like proteins (CBLs) are plant-specific Ca²⁺ sensors that interact with CBL-interacting protein kinases (CIPKs) to transmit Ca²⁺ signals. CBL-CIPK complexes have been reported to play pivotal roles in plant development and response to drought stress; however, limited information is available about the CBL and CIPK genes in pecan, an important nut crop.

RESULTS

In the present study, a total of 9 CBL and 30 CIPK genes were identified from the pecan genome and divided into four and five clades based on phylogeny, respectively. Gene structure and distribution of conserved sequence motif analysis suggested that family members in the same clade commonly exhibited similar exon-intron structures and motif compositions. The segmental duplication events contributed largely to the expansion of pecan CBL and CIPK gene families, and Ka/Ks values revealed that all of them experienced strong negative selection. Phylogenetic analysis of CIPK proteins from 14 plant species revealed that CIPKs in the intron-poor clade originated in seed plants. Tissue-specific expression profiles of CiCBLs and CiCIPKs were analysed, presenting functional diversity. Expression profiles derived from RNA-Seq revealed distinct expression patterns of CiCBLs and CiCIPKs under drought treatment in pecan. Moreover, coexpression network analysis helped to elucidate the relationships between these genes and identify potential candidates for the regulation of drought response, which were verified by qRT-PCR analysis.

CONCLUSIONS

The characterization and analysis of CBL and CIPK genes in pecan genome could provide a basis for further functional analysis of CiCBLs and CiCIPKs in the drought stress response of pecan.

Genome-Wide Identification of bZIP Transcription Factor Genes and Functional Analyses of Two Members in Cytospora chrysosperma

The basic leucine zipper (bZIP) transcription factor (TF) family, one of the largest and the most diverse TF families, is widely distributed across the eukaryotes. It has been described that the bZIP TFs play diverse roles in development, nutrient utilization, and various stress responses in fungi. However, little is known of the bZIP members in Cytospora chrysosperma, a notorious plant pathogenic fungus, which causes canker disease on over 80 woody plant species. In this study, 26 bZIP genes were systematically identified in the genome of C. chrysosperma, and two of them (named CcbZIP05 and CcbZIP23) significantly down-regulated in CcPmk1 deletion mutant (a pathogenicity-related mitogen-activated protein kinase) were selected for further analysis. Deletion of CcbZIP05 or CcbZIP23 displayed a dramatic reduction in fungal growth but showed increased hypha branching and resistance to cell wall inhibitors and abiotic stresses. The CcbZIP05 deletion mutants but not CcbZIP23 deletion mutants were more sensitive to the hydrogen peroxide compared to the wild-type and complemented strains. Additionally, the CcbZIP23 deletion mutants produced few pycnidia but more pigment. Remarkably, both CcbZIP05 and CcbZIP23 deletion mutants were significantly reduced in fungal virulence. Further analysis showed that CcbZIP05 and CcbZIP23 could regulate the expression of putative effector genes and chitin synthesis-related genes. Taken together, our results suggest that CcbZIP05 and CcbZIP23 play important roles in fungal growth, abiotic stresses response, and pathogenicity, which will provide comprehensive information on the CcbZIP genes and lay the foundation for further research on the bZIP members in C. chrysosperma.

Genome-Wide Identification of CBL-CIPK Gene Family in Honeysuckle (Lonicera japonica Thunb.) and Their Regulated Expression Under Salt Stress

In plants, calcineurin B-like proteins (CBLs) are a unique group of Ca²⁺ sensors that decode Ca²⁺ signals by activating a family of plant-specific protein kinases known as CBL-interacting protein kinases (CIPKs). CBL-CIPK gene families and their interacting complexes are involved in regulating plant responses to various environmental stimuli. To gain insight into the functional divergence of CBL-CIPK genes in honeysuckle, a total of six LjCBL and 17 LjCIPK genes were identified. The phylogenetic analysis along with the gene structure analysis divided both CBL and CBL-interacting protein kinase genes into four subgroups and validated by the distribution of conserved protein motifs. The 3-D structure prediction of proteins shown that most LjCBLs shared the same Protein Data Bank hit 1uhnA and most LjCIPKs shared the 6c9Da. Analysis of cis-acting elements and gene ontology implied that both LjCBL and LjCIPK genes could be involved in hormone signal responsiveness and stress adaptation. Protein-protein interaction prediction suggested that LjCBL4 is hypothesized to interact with LjCIPK7/9/15/16 and SOS1/NHX1. Gene expression analysis in response to salinity stress revealed that LjCBL2/4, LjCIPK1/15/17 under all treatments gradually increased over time until peak expression at 72 h. These results demonstrated the conservation of salt overly sensitive pathway genes in honeysuckle and a model of Ca²⁺-LjCBL4/LjSOS3-LjCIPK16/LjSOS2 module-mediated salt stress signaling in honeysuckle is proposed. This study provides insight into the characteristics of the CBL-CIPK gene families involved in honeysuckle salt stress responses, which could serve as a foundation for gene transformation technology, to obtain highly salt-tolerant medicinal plants in the context of the global reduction of cultivated land.

The temperature-regulated DEAD-box RNA helicase CrhR interactome: Autoregulation and photosynthesis-related transcripts

RNA helicases play crucial functions in RNA biology. In plants, RNA helicases are encoded by large gene families, performing roles in abiotic stress responses, development, the post-transcriptional regulation of gene expression as well as house-keeping functions. Several of these RNA helicases are targeted to the organelles, mitochondria and chloroplasts. Cyanobacteria are the direct evolutionary ancestors of plant chloroplasts. The cyanobacterium Synechocystis 6803 encodes a single DEAD-box RNA helicase, CrhR, that is induced by a range of abiotic stresses, including low temperature. Though the ΔcrhR mutant exhibits a severe cold-sensitive phenotype, the physiological function(s) performed by CrhR have not been described. To identify transcripts interacting with CrhR, we performed RNA co-immunoprecipitation with extracts from a Synechocystis crhR deletion mutant expressing the FLAG-tagged native CrhR or a K57A mutated version with an anticipated enhanced RNA binding. The composition of the interactome was strikingly biased towards photosynthesis-associated and redox-controlled transcripts. A transcript highly enriched in all experiments was the crhR mRNA, suggesting an auto-regulatory molecular mechanism. The identified interactome explains the described physiological role of CrhR in response to the redox poise of the photosynthetic electron transport chain and characterizes CrhR as an enzyme with a diverse range of transcripts as molecular targets.

Identification and Characterization of Abiotic Stress Responsive CBL-CIPK Family Genes in Medicago

The calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) play important roles in plant signal transduction and response to abiotic stress. Plants of Medicago genus contain many important forages, and their growth is often affected by a variety of abiotic stresses. However, studies on the CBL and CIPK family member and their function are rare in Medicago. In this study, a total of 23 CBL and 58 CIPK genes were identified from the genome of Medicago sativa as an important forage crop, and Medicaog truncatula as the model plant. Phylogenetic analysis suggested that these CBL and CIPK genes could be classified into five and seven groups, respectively. Moreover, these genes/proteins showed diverse exon-intron organizations, architectures of conserved protein motifs. Many stress-related cis-acting elements were found in their promoter region. In addition, transcriptional analyses showed that these CBL and CIPK genes exhibited distinct expression patterns in various tissues, and in response to drought, salt, and abscisic acid treatments. In particular, the expression levels of MtCIPK2 (MsCIPK3), MtCIPK17 (MsCIPK11), and MtCIPK18 (MsCIPK12) were significantly increased under PEG, NaCl, and ABA treatments. Collectively, our study suggested that CBL and CIPK genes play crucial roles in response to various abiotic stresses in Medicago.

Genome-wide identification and characterization of DEAD-box helicase family associated with early somatic embryogenesis in Dimocarpus longan Lour

DEAD-box (DDX) proteins belong to the largest subfamily of RNA helicase SF2, which contributes to all biological processes of RNA metabolism in the plant kingdom. Till now, no significant data are available regarding studies on DDX in Somatic Embryogenesis (SE) of woody plants. It is important to investigate the biological function of the DlDDX family in longan SE. Thus, a comprehensive analysis of 58 longan DEAD-box (DlDDX) genes characterization was performed by genome-wide identification and transcript abundance validation analysis. Homologous evolution has revealed that some DlDDXs in longan had high sequence similarity with Mus musculus, Citrus and Saccharomyces cerevisiae, indicating that DlDDXs were highly conservative in the animal, plant, and microorganism. Remarkably, gene duplication, purifying selection, and alternative splicing events, and new auxiliary domains have likely contributed to the functional evolution of DlDDX, indicating that DlDDX appeared neofunctionalization in longan. Besides, DlDDX3, 15, 28, 36 might interact with protein complex (MAC3A, MAC3B, CDC5, CBP20) of miRNA biosynthesis. Notably, DlDDX28 contained a novel auxiliary domain (CAF-1 p150), which might contribute to DNA demethylation in longan early SE. 4 DlDDX genes significantly expressed not only in early SE and zygotic embryogenesis (ZE) but also up-regulated at high levels in 'Honghezi' and 'Quanlongbaihe' with abortive seeds, which are of great significance. Moreover, some DlDDXs presented abiotic stress-response dynamic expression patterns by ABA, SA, JA, and NaCl treatments during early SE. Hence, DEAD-box is essential to SE development and seed abortive in longan.

Systematic Analysis of Gibberellin Pathway Components in Medicago truncatula Reveals the Potential Application of Gibberellin in Biomass Improvement

Gibberellins (GAs), a class of phytohormones, act as an essential natural regulator of plant growth and development. Many studies have shown that GA is related to rhizobial infection and nodule organogenesis in legume species. However, thus far, GA metabolism and signaling components are largely unknown in the model legume Medicago truncatula. In this study, a genome-wide analysis of GA metabolism and signaling genes was carried out. In total 29 components, including 8 MtGA20ox genes, 2 MtGA3ox genes, 13 MtGA2ox genes, 3 MtGID1 genes, and 3 MtDELLA genes were identified in M. truncatula genome. Expression profiles revealed that most members of MtGAox, MtGID1, and MtDELLA showed tissue-specific expression patterns. In addition, the GA biosynthesis and deactivation genes displayed a feedback regulation on GA treatment, respectively. Yeast two-hybrid assays showed that all the three MtGID1s interacted with MtDELLA1 and MtDELLA2, suggesting that the MtGID1s are functional GA receptors. More importantly, M. truncatula exhibited increased plant height and biomass by ectopic expression of the MtGA20ox1, suggesting that enhanced GA response has the potential for forage improvement.

Genome-Wide Identification and Expression Analysis of the bZIP Transcription Factors in the Mycoparasite Coniothyrium minitans

The basic leucine zipper (bZIP) proteins family is one of the largest and most diverse transcription factors, widely distributed in eukaryotes. However, no information is available regarding the bZIP gene family in Coniothyrium minitans, an important biocontrol agent of the plant pathogen Sclerotinia sclerotiorum. In this study, we identified 34 bZIP genes from the C. minitans genome, which were classified into 8 groups based on their phylogenetic relationships. Intron analysis showed that 28 CmbZIP genes harbored a variable number of introns, and 15 of them shared a feature that intron inserted into the bZIP domain. The intron position in bZIP domain was highly conserved, which was related to recognize the arginine (R) and could be treated as a genomic imprinting. Expression analysis of the CmbZIP genes in response to abiotic stresses indicated that they might play distinct roles in abiotic stress responses. Results showed that 22 CmbZIP genes were upregulated during the later stage of conidial development. Furthermore, transcriptome analysis indicated that CmbZIP genes are involved in different stages of mycoparasitism. Among deletion mutants of four CmbZIPs (CmbZIP07, -09, -13, and -16), only ΔCmbZIP16 mutants significantly reduced its tolerance to the oxidative stress. The other mutants exhibited no significant effects on colony morphology, mycelial growth, conidiation, and mycoparasitism. Taken together, our results suggested that CmbZIP genes play important roles in the abiotic stress responses, conidial development, and mycoparasitism. These results provide comprehensive information of the CmbZIP gene family and lay the foundation for further research on the bZIP gene family regarding their biological functions and evolutionary history.

Gametophyte-specific DEAD-box RNA helicase 29 is required for functional maturation of male and female gametophytes in Arabidopsis

The maturation of male and female gametophytes together with its impact on plant sexual reproduction has not received much attention, and the molecular mechanisms underlying the process are largely unknown. Here, we show that Arabidopsis DEAD-box RNA helicase 29 (RH29) is critical for the functional maturation of both male and female gametophytes. Homozygous rh29 mutants could not be obtained, and heterozygous mutant plants were semi-sterile. Progression of the cell cycle in rh29 female gametophytes was delayed. Delayed pollination experiments showed that rh29 female gametophytes underwent cell-fate specification but were unable to develop into functional gametophytes. Functional specification but not morphogenesis was also disrupted in rh29 male gametophytes, causing defective pollen tube growth in the pistil. RH29 was highly and specifically expressed in gametophytic cells. RH29 shares high amino acid sequence identity with yeast Dbp10p, which partially rescues the aborted-ovules phenotype of rh29/RH29 plants. RH29 is essential for the synthesis of REGULATORY PARTICLE TRIPLE A ATPase 5a (RPT5a), a subunit of the regulatory particle of the 26S proteasome. Our results suggest that gametophyte functional maturation is a necessary process for successful fertilization and that RH29 is essential for the functional maturation of both male and female gametophytes.

Genome-wide identification, and characterization of the CDPK gene family reveal their involvement in abiotic stress response in Fragaria x ananassa

Calcium-dependent protein kinases (CDPKs) are encoded by a large gene family and play important roles against biotic and abiotic stresses and in plant growth and development. To date, little is known about the CDPK genes in strawberry (Fragaria x ananassa). In this study, analysis of Fragaria x ananassa CDPK gene family was performed, including gene structures, phylogeny, interactome and expression profiles. Nine new CDPK genes in Fragaria x ananassa were identified based on RNA-seq data. These identified strawberry FaCDPK genes were classified into four main groups, based on the phylogenetic analysis and structural features. FaCDPK genes were differentially expressed during fruit development and ripening, as well as in response to abiotic stress (salt and drought), and hormone (abscisic acid) treatment. In addition, the interaction network analysis pointed out proteins involved in the ABA-dependent response to plant stress via Ca2+ signaling, especially RBOHs. To our knowledge, this is the first report on CDPK families in Fragaria x ananassa, and it will provide valuable information for development of biofortified fruits and stress tolerant plants.

Evolutionary Analysis of Calcium-Dependent Protein Kinase in Five Asteraceae Species

Calcium-dependent protein kinase (CPK) is crucial in Ca²⁺ signal transduction, and is a large gene family in plants. In our previous work, we reported Hevea brasiliensis CPKs were important for natural rubber biosynthesis. However, this CPK gene family in other rubber producing plants has not been investigated. Here, we report the CPKs in five representative Asteraceae species, including three rubber-producing and two non-rubber species. A total of 34, 34, 40, 34 and 30 CPKs were identified from Taraxacum koksaghyz, Lactuca sativa, Helianthus annuus, Chrysanthemum nankingense and Cynara cardunculus, respectively. All CPKs were classified into four individual groups (group I to IV). In addition, 10 TkCPKs, 11 LsCPKs, 20 HaCPKs, 13 CnCPKs and 7 CcCPKs duplicated paralogs were identified. Further evolutionary analysis showed that, compared to other subfamilies, the group III had been expanded in the Asteraceae species, especially in the rubber-producing species. Meanwhile, the CPKs in group III from Asteraceae species tend to expand with low calcium binding capacity. This study provides a systematical evolutionary investigation of the CPKs in five representative Asteraceae species, suggesting that the sub-family specific expansion of CPKs might be related to natural rubber producing.

Genome-Wide Identification and Expression Profiling of CBL-CIPK Gene Family in Pineapple (Ananas comosus) and the Role of AcCBL1 in Abiotic and Biotic Stress Response

Ca²⁺ serves as a ubiquitous second messenger regulating several aspects of plant growth and development. A group of unique calcium sensor proteins, calcineurin B-like (CBL), interact with CBL-interacting protein kinases (CIPKs) to decode the Ca²⁺ signature inside the cell. Although CBL-CIPK signaling toolkit has been shown to play significant roles in the responses to numerous stresses in different plants, the information about pineapple CBL-CIPK remains obscure. In the present study, a total of eight AcCBL and 21 AcCIPK genes were identified genome-wide in pineapple. The identified genes were renamed on the basis of gene ID in ascending order and phylogenetic analysis divided into five groups. Transcriptomic data analysis showed that AcCBL and AcCIPK genes were expressed differentially in different tissues. Further, the expression analysis of AcCBL1 in different tissues showed significant changes under various abiotic stimuli. Additionally, the ectopic expression of AcCBL1 in Arabidopsis resulted in enhanced tolerance to salinity, osmotic, and fungal stress. The present study revealed the crucial contribution of the CBL-CIPK gene in various biological and physiological processes in pineapple.

Genome-wide identification and expression analysis of sucrose synthase genes in allotetraploid Brassica juncea

Sucrose plays pivotal role in energy metabolism and regulating gene expression of several physiological processes in higher plants. Here, fourteen sucrose synthase (SUS) genes have been identified in the allotetraploid genome of Indian mustard, Brassica juncea. The identified SUS genes in B. juncea (BjSUS) were derived from the two-progenitor species, B. rapa and B. nigra. Intron-exon analysis indicated loss or gain of 1-3 introns in diversification of SUS gene family. Phylogenetic analysis revealed discrete evolutionary paths for the BjSUS genes, originating from three ancestor groups, SUS I, SUS II and SUS III. Gene expression study revealed significant variability in expression of the BjSUS paralogs across the different tissues. BjSUS genes showed transcriptional activation in response to defense hormones and a late response to wounding. Tissue and temporal specificity of expression revealed importance of specific SUS paralogs at different developmental stages and under different stress conditions. The study highlighted differential involvement of SUS paralogs in sucrose metabolism across the tissues and stress-responses, in a major oilseed crop B. juncea.

Transcriptional and Post-transcriptional Regulation of Organellar Gene Expression (OGE) and Its Roles in Plant Salt Tolerance

Given their endosymbiotic origin, chloroplasts and mitochondria genomes harbor only between 100 and 200 genes that encode the proteins involved in organellar gene expression (OGE), photosynthesis, and the electron transport chain. However, as the activity of these organelles also needs a few thousand proteins encoded by the nuclear genome, a close coordination of the gene expression between the nucleus and organelles must exist. In line with this, OGE regulation is crucial for plant growth and development, and is achieved mainly through post-transcriptional mechanisms performed by nuclear genes. In this way, the nucleus controls the activity of organelles and these, in turn, transmit information about their functional state to the nucleus by modulating nuclear expression according to the organelles' physiological requirements. This adjusts organelle function to plant physiological, developmental, or growth demands. Therefore, OGE must appropriately respond to both the endogenous signals and exogenous environmental cues that can jeopardize plant survival. As sessile organisms, plants have to respond to adverse conditions to acclimate and adapt to them. Salinity is a major abiotic stress that negatively affects plant development and growth, disrupts chloroplast and mitochondria function, and leads to reduced yields. Information on the effects that the disturbance of the OGE function has on plant tolerance to salinity is still quite fragmented. Nonetheless, many plant mutants which display altered responses to salinity have been characterized in recent years, and interestingly, several are affected in nuclear genes encoding organelle-localized proteins that regulate the expression of organelle genes. These results strongly support a link between OGE and plant salt tolerance, likely through retrograde signaling. Our review analyzes recent findings on the OGE functions required by plants to respond and tolerate salinity, and highlights the fundamental role that chloroplast and mitochondrion homeostasis plays in plant adaptation to salt stress.

Genome-wide identification and functional analysis of the splicing component SYF2/NTC31/p29 across different plant species

This study systematically identifies plant SYF2/NTC31/p29 genes from 62 plant species by a combinatory bioinformatics approach, revealing the importance of this gene family in phylogenetics, duplication, transcriptional, and post-transcriptional regulation. Alternative splicing is a post-transcriptional regulatory mechanism, which is critical for plant development and stress responses. The entire process is strictly attenuated by a complex of splicing-related proteins, designated splicing factors. Human p29, also referred to as synthetic lethal with cdc forty 2 (SYF2) or the NineTeen complex 31 (NTC31), is a core protein found in the NTC complex of humans and yeast. This splicing factor participates in a variety of biological processes, including DNA damage repair, control of the cell cycle, splicing, and tumorigenesis. However, its function in plants has been seldom reported. Thus, we have systematically identified 89 putative plant SYF2s from 62 plant species among the deposited entries in the Phytozome database. The phylogenetic relationships and evolutionary history among these plant SYF2s were carefully examined. The results revealed that plant SYF2s exhibited distinct patterns regarding their gene structure, promoter sequences, and expression levels, suggesting their functional diversity in response to developmental cues or stress treatments. Although local duplication events, such as tandem duplication and retrotransposition, were found among several plant species, most of the plant species contained only one copy of SYF2, suggesting the existence of additional mechanisms to confer duplication resistance. Further investigation using the model dicot and monocot representatives Arabidopsis and rice SYF2s indicated that the splicing pattern and resulting protein isoforms might play an alternative role in the functional diversity.

Comprehensive Analysis of Cucumber Gibberellin Oxidase Family Genes and Functional Characterization of CsGA20ox1 in Root Development in Arabidopsis

Cucumber (Cucumis sativus L.) is an important vegetable crop worldwide and gibberellins (GAs) play important roles in the regulation of cucumber developmental and growth processes. GA oxidases (GAoxs), which are encoded by different gene subfamilies, are particularly important in regulating bioactive GA levels by catalyzing the later steps in the biosynthetic pathway. Although GAoxs are critical enzymes in GA synthesis pathway, little is known about GAox genes in cucumber, in particular about their evolutionary relationships, expression profiles and biological function. In this study, we identified 17 GAox genes in cucumber genome and classified them into five subfamilies based on a phylogenetic tree, gene structures, and conserved motifs. Synteny analysis indicated that the tandem duplication or segmental duplication events played a minor role in the expansion of cucumber GA2ox, GA3ox and GA7ox gene families. Comparative syntenic analysis combined with phylogenetic analysis provided deep insight into the phylogenetic relationships of CsGAox genes and suggested that protein homology CsGAox are closer to AtGAox than OsGAox. In addition, candidate transcription factors BBR/BPC (BARLEY B RECOMBINANT/BASIC PENTACYSTEINE) and GRAS (GIBBERELLIC ACID-INSENSITIVE, REPRESSOR of GAI, and SCARECROW) which may directly bind promoters of CsGAox genes were predicted. Expression profiles derived from transcriptome data indicated that some CsGAox genes, especially CsGA20ox1, are highly expressed in seedling roots and were down-regulated under GA₃ treatment. Ectopic over-expression of CsGA20ox1 in Arabidopsis significantly increased primary root length and lateral root number. Taken together, comprehensive analysis of CsGAoxs would provide a basis for understanding the evolution and function of the CsGAox family.

Genome-wide identification and analysis of the basic leucine zipper (bZIP) transcription factor gene family in Ustilaginoidea virens

The basic leucine zipper (bZIP) transcription factor (TF) family is one of the largest and most diverse TF families widely distributed across the eukaryotes. The bZIP TF family plays an important role in growth, development, and response to abiotic or biotic stresses, which have been well characterized in plants, but not in plant pathogenic fungi. In this study, we performed genome-wide and systematic bioinformatics analysis of bZIP genes in the fungus Ustilaginoidea virens, the causal agent of rice false smut disease. We identified 28 bZIP family members in the U. virens genome by searching for the bZIP domain in predicted genes. The gene structures, motifs, and phylogenetic relationships were analyzed for bZIP genes in U. virens (UvbZIP). Together with bZIP proteins from two other fungi, the bZIP genes can be divided into eight groups according to their phylogenetic relationships. Based on RNA-Seq data, the expression profiles of UvbZIP genes at different infection stages were evaluated. Results showed that 17 UvbZIP genes were up-regulated during the infection period. Furthermore, 11 infection-related UvbZIP genes were investigated under H2O2 stress and the expression level of eight genes were changed, which confirmed their role in stress tolerance and pathogenicity. In summary, our genome-wide systematic characterization and expression analysis of UvbZIP genes provided insight into the molecular function of these genes in U. virens and provides a reference for other pathogens.

Genome-wide survey indicates diverse physiological roles of the turnip (Brassica rapa var. rapa) calcium-dependent protein kinase genes

Calcium-dependent protein kinases (CDPKs) as crucial sensors of calcium concentration changes play important roles in responding to abiotic and biotic stresses. In this study, 55 BrrCDPK genes, which were phylogenetically clustered into four subfamilies, were identified. Chromosome locations indicated that the CDPK family in turnip expanded by segmental duplication and genome rearrangement. Moreover, gene expression profiles showed that different BrrCDPKs were expressed in specific tissues or stages. Transcript levels of BrrCDPKs indicated that they were involved in abiotic and biotic stresses and that paralogs exhibited functional divergence. Additionally, we identified 15 Rboh genes in turnip; the results of yeast two-hybrid analysis suggested that BrrRbohD1 interacted only with BrrCDPK10 and that BrrRbohD2 interacted with BrrCDPK4/7/9/10/17/22/23. Most of the genes play an important role in pst DC3000 defense by regulating the accumulation of H2O2 and stomatal closure. Our study may provide an important foundation for future functional analysis of BrrCDPKs and reveal further biological roles.

Genome and transcriptome-wide analyses of cellulose synthase gene superfamily in soybean

The plant cellulose synthase gene superfamily belongs to the category of type-2 glycosyltransferases, and is involved in cellulose and hemicellulose biosynthesis. These enzymes are vital for maintaining cell-wall structural integrity throughout plant life. Here, we identified 78 putative cellulose synthases (CS) in the soybean genome. Phylogenetic analysis against 40 reference Arabidopsis CS genes clustered soybean CSs into seven major groups (CESA, CSL A, B, C, D, E and G), located on 19 chromosomes (except chromosome 18). Soybean CS expansion occurred in 66 duplication events. Additionally, we identified 95 simple sequence repeat makers related to 44 CSs. We next performed digital expression analysis using publically available datasets to understand potential CS functions in soybean. We found that CSs were highly expressed during soybean seed development, a pattern confirmed with an Affymatrix soybean IVT array and validated with RNA-seq profiles. Within CS groups, CESAs had higher relative expression than CSLs. Soybean CS models were designed based on maximum average RPKM values. Gene co-expression networks were developed to explore which CSs could work together in soybean. Finally, RT-PCR analysis confirmed the expression of 15 selected CSs during all four seed developmental stages.

Chlorosis caused by two recessively interacting genes reveals a role of RNA helicase in hybrid breakdown in Arabidopsis thaliana

Hybrids often differ in fitness from their parents. They may be superior, translating into hybrid vigour or heterosis, but they may also be markedly inferior, because of hybrid weakness or incompatibility. The underlying genetic causes for the latter can often be traced back to genes that evolve rapidly because of sexual or host-pathogen conflicts. Hybrid weakness may manifest itself only in later generations, in a phenomenon called hybrid breakdown. We have characterized a case of hybrid breakdown among two Arabidopsis thaliana accessions, Shahdara (Sha, Tajikistan) and Lövvik-5 (Lov-5, Northern Sweden). In addition to chlorosis, a fraction of the F₂ plants have defects in leaf and embryo development, and reduced photosynthetic efficiency. Hybrid chlorosis is due to two major-effect loci, of which one, originating from Lov-5, appears to encode an RNA helicase (AtRH18). To examine the role of the chlorosis allele in the Lövvik area, in addition to eight accessions collected in 2009, we collected another 240 accessions from 15 collections sites, including Lövvik, from Northern Sweden in 2015. Genotyping revealed that Lövvik collection site is separated from the rest. Crosses between 109 accessions from this area and Sha revealed 85 cases of hybrid chlorosis, indicating that the chlorosis-causing allele is common in this area. These results suggest that hybrid breakdown alleles not only occur at rapidly evolving loci, but also at genes that code for conserved processes.

Chloroplast- or Mitochondria-Targeted DEAD-Box RNA Helicases Play Essential Roles in Organellar RNA Metabolism and Abiotic Stress Responses

The yields and productivity of crops are greatly diminished by various abiotic stresses, including drought, cold, heat, and high salinity. Chloroplasts and mitochondria are cellular organelles that can sense diverse environmental stimuli and alter gene expression to cope with adverse environmental stresses. Organellar gene expression is mainly regulated at posttranscriptional levels, including RNA processing, intron splicing, RNA editing, RNA turnover, and translational control, during which a variety of nucleus-encoded RNA-binding proteins (RBPs) are targeted to chloroplasts or mitochondria where they play essential roles in organellar RNA metabolism. DEAD-box RNA helicases (RHs) are enzymes that can alter RNA structures and affect RNA metabolism in all living organisms. Although a number of DEAD-box RHs have been found to play important roles in RNA metabolism in the nucleus and cytoplasm, our understanding on the roles of DEAD-box RHs in the regulation of RNA metabolism in chloroplasts and mitochondria is only at the beginning. Considering that organellar RNA metabolism and gene expression are tightly regulated by anterograde signaling from the nucleus, it is imperative to determine the functions of nucleus-encoded organellar RBPs. In this review, we summarize the emerging roles of nucleus-encoded chloroplast- or mitochondria-targeted DEAD-box RHs in organellar RNA metabolism and plant response to diverse abiotic stresses.

The CBL and CIPK Gene Family in Grapevine (Vitis vinifera): Genome-Wide Analysis and Expression Profiles in Response to Various Abiotic Stresses

Calcium plays a central role in regulating signal transduction pathways. Calcineurin B-like proteins (CBLs), which harbor a crucial region consisting of EF hands that capture Ca²⁺, interact in a specific manner with CBL-interacting protein kinases (CIPKs). This two gene families or their interacting-complex widely respond to various environment stimuli and development processes. The genome-wide annotation and specific expression patterns of CBLs and CIPKs, however, in grapevine remain unclear. In the present study, eight CBL and 20 CIPK genes were identified in grapevine genome, and divided into four and five subfamilies, respectively, based on phylogenetic analysis, and validated by gene structure and the distribution of conserved protein motifs. Four (50%) out of eight VvCBLs and eight (40%) out of 20 VvCIPKs were found to be derived from tandem duplication, and five (25%) out of 20 VvCIPKs were derived from segmental duplication, indicating that the expansion of grapevine CBL and CIPK gene families were mainly contributed by gene duplication, and all duplication events between VvCIPK genes only detected in intron poor clade. Estimating of synonymous and non-synonymous substitution rates of both gene families suggested that VvCBL genes seems more conserved than VvCIPK genes, and were derived by positive selection pressure, whereas VvCIPK genes were mainly derived by purifying selection pressure. Expressional analyses of VvCBL and VvCIPK genes based on microarray and qRT-PCR data performed diverse expression patterns of VvCBLs and VvCIPKs in response to both various abiotic stimuli and at different development stages. Furthermore, the co-expression analysis of grapevine CBLs and CIPKs suggested that CBL-CIPK complex seems to be more responsive to abiotic stimuli than during different development stages. VvCBLs may play an important and special role in regulating low temperature stress. The protein interaction analysis suggested divergent mechanisms might exist between Arabidopsis and grapevine. Our results will facilitate the future functional characterization of individual VvCBLs and VvCIPKs.

eIF4A RNA Helicase Associates with Cyclin-Dependent Protein Kinase A in Proliferating Cells and Is Modulated by Phosphorylation

Eukaryotic initiation factor 4A (eIF4A) is a highly conserved RNA-stimulated ATPase and helicase involved in the initiation of messenger RNA translation. Previously, we found that eIF4A interacts with cyclin-dependent kinase A (CDKA), the plant ortholog of mammalian CDK1. Here, we show that this interaction occurs only in proliferating cells where the two proteins coassociate with 5'-cap-binding protein complexes, eIF4F or the plant-specific eIFiso4F. CDKA phosphorylates eIF4A on a conserved threonine residue (threonine-164) within the RNA-binding motif 1b TPGR. In vivo, a phospho-null (APGR) variant of the Arabidopsis (Arabidopsis thaliana) eIF4A1 protein retains the ability to functionally complement a mutant (eif4a1) plant line lacking eIF4A1, whereas a phosphomimetic (EPGR) variant fails to complement. The phospho-null variant (APGR) rescues the slow growth rate of roots and rosettes, together with the ovule-abortion and late-flowering phenotypes. In vitro, wild-type recombinant eIF4A1 and its phospho-null variant both support translation in cell-free wheat germ extracts dependent upon eIF4A, but the phosphomimetic variant does not support translation and also was deficient in ATP hydrolysis and helicase activity. These observations suggest a mechanism whereby CDK phosphorylation has the potential to down-regulate eIF4A activity and thereby affect translation.

Recruitment of Arabidopsis RNA Helicase AtRH9 to the Viral Replication Complex by Viral Replicase to Promote Turnip Mosaic Virus Replication

Positive-sense RNA viruses have a small genome with very limited coding capacity and are highly dependent on host components to fulfill their life cycle. Recent studies have suggested that DEAD-box RNA helicases play vital roles in many aspects of RNA metabolism. To explore the possible role of the RNA helicases in viral infection, we used the Turnip mosaic virus (TuMV)-Arabidopsis pathosystem. The Arabidopsis genome encodes more than 100 putative RNA helicases (AtRH). Over 41 Arabidopsis T-DNA insertion mutants carrying genetic lesions in the corresponding 26 AtRH genes were screened for their requirement in TuMV infection. TuMV infection assays revealed that virus accumulation significantly decreased in the Arabidopsis mutants of three genes, AtRH9, AtRH26, and PRH75. In the present work, AtRH9 was further characterized. Yeast two-hybrid and bimolecular fluorescence complementation (BiFC) assays showed that AtRH9 interacted with the TuMV NIb protein, the viral RNA-dependent RNA polymerase. Moreover, the subcellular distribution of AtRH9 was altered in the virus-infected cells, and AtRH9 was recruited to the viral replication complex. These results suggest that Arabidopsis AtRH9 is an important component of the TuMV replication complex, possibly recruited via its interaction with NIb.

Analysis of basic leucine zipper genes and their expression during bud dormancy in peach (Prunus persica)

Dormancy is a biological characteristic developed to resist the cold conditions in winter. The bZIP transcription factors are present exclusively in eukaryotes and have been identified and classified in many species. bZIP proteins are known to regulate numerous biological processes, however, the role of bZIP in bud dodormancy has not been studied extensively. In total, 50 PpbZIP transcription factor-encoding genes were identified and categorized them into 10 groups (A-I and S). Similar intron/exon structures, additional conserved motifs, and DNA-binding site specificity supported our classification scheme. Additionally, chromosomal distribution and collinearity analyses suggested that expansion of the PpbZIP transcription factor family was due to segment/chromosomal duplications. We also predicted the dimerization properties based on characteristic features of the leucine zipper and classified PpbZIP proteins into 23 subfamilies. Furthermore, qRT-PCR results indicated that PpbZIPs genes may be involved in regulating dormancy. The same gene of different species might participate in different regulating networks through interactions with specific partners. Our expression profiling results complemented the microarray data, suggesting that co-expression patterns of bZIP transcription factors during dormancy differed among deciduous fruit trees. Our findings further clarify the molecular characteristics of the PpbZIP transcription factor family, including potential gene functions during dormancy. This information may facilitate further research on the evolutionary history and biological functions of bZIP proteins in peach and other rosaceae plants.

A Cold-Inducible DEAD-Box RNA Helicase from Arabidopsis thaliana Regulates Plant Growth and Development under Low Temperature

DEAD-box RNA helicases comprise a large family and are involved in a range of RNA processing events. Here, we identified one of the Arabidopsis thaliana DEAD-box RNA helicases, AtRH7, as an interactor of Arabidopsis COLD SHOCK DOMAIN PROTEIN 3 (AtCSP3), which is an RNA chaperone involved in cold adaptation. Promoter:GUS transgenic plants revealed that AtRH7 is expressed ubiquitously and that its levels of the expression are higher in rapidly growing tissues. Knockout mutant lines displayed several morphological alterations such as disturbed vein pattern, pointed first true leaves, and short roots, which resemble ribosome-related mutants of Arabidopsis. In addition, aberrant floral development was also observed in rh7 mutants. When the mutants were germinated at low temperature (12°C), both radicle and first leaf emergence were severely delayed; after exposure of seedlings to a long period of cold, the mutants developed aberrant, fewer, and smaller leaves. RNA blots and circular RT-PCR revealed that 35S and 18S rRNA precursors accumulated to higher levels in the mutants than in WT under both normal and cold conditions, suggesting the mutants are partially impaired in pre-rRNA processing. Taken together, the results suggest that AtRH7 affects rRNA biogenesis and plays an important role in plant growth under cold.

The DEAD-Box RNA Helicase AtRH7/PRH75 Participates in Pre-rRNA Processing, Plant Development and Cold Tolerance in Arabidopsis

DEAD-box RNA helicases belong to an RNA helicase family that plays specific roles in various RNA metabolism processes, including ribosome biogenesis, mRNA splicing, RNA export, mRNA translation and RNA decay. This study investigated a DEAD-box RNA helicase, AtRH7/PRH75, in Arabidopsis. Expression of AtRH7/PRH75 was ubiquitous; however, the levels of mRNA accumulation were increased in cell division regions and were induced by cold stress. The phenotypes of two allelic AtRH7/PRH75-knockout mutants, atrh7-2 and atrh7-3, resembled auxin-related developmental defects that were exhibited in several ribosomal protein mutants, and were more severe under cold stress. Northern blot and circular reverse transcription-PCR (RT-PCR) analyses indicated that unprocessed 18S pre-rRNAs accumulated in the atrh7 mutants. The atrh7 mutants were hyposensitive to the antibiotic streptomycin, which targets ribosomal small subunits, suggesting that AtRH7 was also involved in ribosome assembly. In addition, the atrh7-2 and atrh7-3 mutants displayed cold hypersensitivity and decreased expression of CBF1, CBF2 and CBF3, which might be responsible for the cold intolerance. The present study indicated that AtRH7 participates in rRNA biogenesis and is also involved in plant development and cold tolerance in Arabidopsis.

Genome-wide analysis of the basic leucine zipper (bZIP) transcription factor gene family in six legume genomes

BACKGROUND

Plant bZIP proteins characteristically harbor a highly conserved bZIP domain with two structural features: a DNA-binding basic region and a leucine (Leu) zipper dimerization region. They have been shown to be diverse transcriptional regulators, playing crucial roles in plant development, physiological processes, and biotic/abiotic stress responses. Despite the availability of six completely sequenced legume genomes, a comprehensive investigation of bZIP family members in legumes has yet to be presented.

RESULTS

In this study, we identified 428 bZIP genes encoding 585 distinct proteins in six legumes, Glycine max, Medicago truncatula, Phaseolus vulgaris, Cicer arietinum, Cajanus cajan, and Lotus japonicus. The legume bZIP genes were categorized into 11 groups according to their phylogenetic relationships with genes from Arabidopsis. Four kinds of intron patterns (a-d) within the basic and hinge regions were defined and additional conserved motifs were identified, both presenting high group specificity and supporting the group classification. We predicted the DNA-binding patterns and the dimerization properties, based on the characteristic features in the basic and hinge regions and the Leu zipper, respectively, which indicated that some highly conserved amino acid residues existed across each major group. The chromosome distribution and analysis for WGD-derived duplicated blocks revealed that the legume bZIP genes have expanded mainly by segmental duplication rather than tandem duplication. Expression data further revealed that the legume bZIP genes were expressed constitutively or in an organ-specific, development-dependent manner playing roles in multiple seed developmental stages and tissues. We also detected several key legume bZIP genes involved in drought- and salt-responses by comparing fold changes of expression values in drought-stressed or salt-stressed roots and leaves.

CONCLUSIONS

In summary, this genome-wide identification, characterization and expression analysis of legume bZIP genes provides valuable information for understanding the molecular functions and evolution of the legume bZIP transcription factor family, and highlights potential legume bZIP genes involved in regulating tissue development and abiotic stress responses.

Divergence and adaptive evolution of the gibberellin oxidase genes in plants

BACKGROUND

The important phytohormone gibberellins (GAs) play key roles in various developmental processes. GA oxidases (GAoxs) are critical enzymes in GA synthesis pathway, but their classification, evolutionary history and the forces driving the evolution of plant GAox genes remain poorly understood.

RESULTS

This study provides the first large-scale evolutionary analysis of GAox genes in plants by using an extensive whole-genome dataset of 41 species, representing green algae, bryophytes, pteridophyte, and seed plants. We defined eight subfamilies under the GAox family, namely C19-GA2ox, C20-GA2ox, GA20ox,GA3ox, GAox-A, GAox-B, GAox-C and GAox-D. Of these, subfamilies GAox-A, GAox-B, GAox-C and GAox-D are described for the first time. On the basis of phylogenetic analyses and characteristic motifs of GAox genes, we demonstrated a rapid expansion and functional divergence of the GAox genes during the diversification of land plants. We also detected the subfamily-specific motifs and potential sites of some GAox genes, which might have evolved under positive selection.

CONCLUSIONS

GAox genes originated very early-before the divergence of bryophytes and the vascular plants and the diversification of GAox genes is associated with the functional divergence and could be driven by positive selection. Our study not only provides information on the classification of GAox genes, but also facilitates the further functional characterization and analysis of GA oxidases.

SlDEAD31, a Putative DEAD-Box RNA Helicase Gene, Regulates Salt and Drought Tolerance and Stress-Related Genes in Tomato

The DEAD-box RNA helicases are involved in almost every aspect of RNA metabolism, associated with diverse cellular functions including plant growth and development, and their importance in response to biotic and abiotic stresses is only beginning to emerge. However, none of DEAD-box genes was well characterized in tomato so far. In this study, we reported on the identification and characterization of two putative DEAD-box RNA helicase genes, SlDEAD30 and SlDEAD31 from tomato, which were classified into stress-related DEAD-box proteins by phylogenetic analysis. Expression analysis indicated that SlDEAD30 was highly expressed in roots and mature leaves, while SlDEAD31 was constantly expressed in various tissues. Furthermore, the expression of both genes was induced mainly in roots under NaCl stress, and SlDEAD31 mRNA was also increased by heat, cold, and dehydration. In stress assays, transgenic tomato plants overexpressing SlDEAD31 exhibited dramatically enhanced salt tolerance and slightly improved drought resistance, which were simultaneously demonstrated by significantly enhanced expression of multiple biotic and abiotic stress-related genes, higher survival rate, relative water content (RWC) and chlorophyll content, and lower water loss rate and malondialdehyde (MDA) production compared to wild-type plants. Collectively, these results provide a preliminary characterization of SlDEAD30 and SlDEAD31 genes in tomato, and suggest that stress-responsive SlDEAD31 is essential for salt and drought tolerance and stress-related gene regulation in plants.