Comprehensive genomics and expression analysis of eceriferum (CER) genes in sunflower (Helianthus annuus)

Decoding drought resilience: a comprehensive exploration of the cotton Eceriferum (CER) gene family and its role in stress adaptation

BACKGROUND

The cuticular wax serves as a primary barrier that protects plants from environmental stresses. The Eceriferum (CER) gene family is associated with wax production and stress resistance.

RESULTS

In a genome-wide identification study, a total of 52 members of the CER family were discovered in four Gossypium species: G. arboreum, G. barbadense, G. raimondii, and G. hirsutum. There were variations in the physicochemical characteristics of the Gossypium CER (GCER) proteins. Evolutionary analysis classified the identified GCERs into five groups, with purifying selection emerging as the primary evolutionary force. Gene structure analysis revealed that the number of conserved motifs ranged from 1 to 15, and the number of exons varied from 3 to 13. Closely related GCERs exhibited similar conserved motifs and gene structures. Analyses of chromosomal positions, selection pressure, and collinearity revealed numerous fragment duplications in the GCER genes. Additionally, nine putative ghr-miRNAs targeting seven G. hirsutum CER (GhCER) genes were identified. Among them, three miRNAs, including ghr-miR394, ghr-miR414d, and ghr-miR414f, targeted GhCER09A, representing the most targeted gene. The prediction of transcription factors (TFs) and the visualization of the regulatory TF network revealed interactions with GhCER genes involving ERF, MYB, Dof, bHLH, and bZIP. Analysis of cis-regulatory elements suggests potential associations between the CER gene family of cotton and responses to abiotic stress, light, and other biological processes. Enrichment analysis demonstrated a robust correlation between GhCER genes and pathways associated with cutin biosynthesis, fatty acid biosynthesis, wax production, and stress response. Localization analysis showed that most GCER proteins are localized in the plasma membrane. Transcriptome and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) expression assessments demonstrated that several GhCER genes, including GhCER15D, GhCER04A, GhCER06A, and GhCER12D, exhibited elevated expression levels in response to water deficiency stress compared to control conditions. The functional identification through virus-induced gene silencing (VIGS) highlighted the pivotal role of the GhCER04A gene in enhancing drought resistance by promoting increased tissue water retention.

CONCLUSIONS

This investigation not only provides valuable evidence but also offers novel insights that contribute to a deeper understanding of the roles of GhCER genes in cotton, their role in adaptation to drought and other abiotic stress and their potential applications for cotton improvement.

Comparative genomics of light harvesting chlorophyll (LHC) gene family and impact of chlorophyll-A contents under drought stress in Helianthus annuus

Drought is one of the main environmental stressors that can alter the water status of plants; negatively affect growth, assimilation, and photosynthesis; and eventually reduce crop yield. We explored the dependence of drought tolerance traits on chlorophyll-A content. Local sunflower cultivars (FH-01, FH-628, FH-633, FH-572, and FH-653) were grown in pots and subjected to drought by withholding water for 10, 15, or 20 d. One month after germination, the leaves of the treated and non-treated plants were collected and subjected to biochemical analyses. Under different water stress levels, the levels of peroxidase (POD), superoxide dismutase (SOD), catalase (CAT), and proline increased, whereas those of chlorophyll-A decreased. Regression analysis clearly found that proline (-0.442), POD (-0.528), SOD (-0.532), and CAT (-0.814) have negative beta coefficient values. Phylogenetic analysis revealed that the LHC gene family is divided into six clades. Subcellular locations indicated that most LHC genes were located in the chloroplast; however, only few genes were present in the peroxisomes and endoplasmic reticulum. Our research found that Arabidopsis thaliana LHC genes were highly homologous to the LHC genes of Helianthus annuus. Furthermore, the LHC genes of both species are located in the chloroplasts; therefore, they play a role in photosynthesis and renewable energy production. This study opens a new horizon for discussing the role of chlorophyll-A in the drought-related traits of sunflowers.

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 of the CER Gene Family and Significant Features in Climate Adaptation of Castanea mollissima

The plant cuticle is the outermost layer of the aerial organs and an important barrier against biotic and abiotic stresses. The climate varies greatly between the north and south of China, with large differences in temperature and humidity, but Chinese chestnut is found in both regions. This study investigated the relationship between the wax layer of chestnut leaves and environmental adaptation. Firstly, semi-thin sections were used to verify that there is a significant difference in the thickness of the epicuticular wax layer between wild chestnut leaves in northwest and southeast China. Secondly, a whole-genome selective sweep was used to resequence wild chestnut samples from two typical regional populations, and significant genetic divergence was identified between the two populations in the CmCER1-1, CmCER1-5 and CmCER3 genes. Thirty-four CER genes were identified in the whole chestnut genome, and a series of predictive analyses were performed on the identified CmCER genes. The expression patterns of CmCER genes were classified into three trends-upregulation, upregulation followed by downregulation and continuous downregulation-when chestnut seedlings were treated with drought stress. Analysis of cultivars from two resource beds in Beijing and Liyang showed that the wax layer of the northern variety was thicker than that of the southern variety. For the Y-2 (Castanea mollissima genome sequencing material) cultivar, there were significant differences in the expression of CmCER1-1, CmCER1-5 and CmCER3 between the southern variety and the northern one-year-grafted variety. Therefore, this study suggests that the CER family genes play a role in environmental adaptations in chestnut, laying the foundation for further exploration of CmCER genes. It also demonstrates the importance of studying the adaptation of Chinese chestnut wax biosynthesis to the southern and northern environments.

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.

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.

Genome-Wide Identification and Expression Profiling of DUF221 Gene Family Provides New Insights Into Abiotic Stress Responses in Potato

The domain of the unknown function 221 proteins regulate several processes in plants, including development, growth, hormone transduction mechanism, and abiotic stress response. Therefore, a comprehensive analysis of the potato genome was conducted to identify the deafness-dystonia peptide (DDP) proteins' role in potatoes. In the present study, we performed a genome-wide analysis of the potato domain of the unknown function 221 (DUF221) genes, including phylogenetic inferences, chromosomal locations, gene duplications, gene structures, and expression analysis. In our results, we identified 10 DDP genes in the potato genome. The phylogenetic analysis results indicated that StDDPs genes were distributed in all four clades, and clade IV was the largest clade. The gene duplication under selection pressure analysis indicated various positive and purifying selections in StDDP genes. The putative stu-miRNAs from different families targeting StDDPs were also predicted in the present study. Promoter regions of StDDP genes contain different cis-acting components involved in multiple stress responses, such as phytohormones and abiotic stress-responsive factors. The analysis of the tissue-specific expression profiling indicated the StDDPs gene expression in stem, root, and leaf tissues. We subsequently observed that StDDP4, StDDP5, and StDDP8 showed higher expressions in roots, stems, and leaves. StDDP5 exhibited high expression against heat stress response, and StDDP7 showed high transcript abundance against salt stress in potatoes. Under abscisic acid (ABA) and indole acetic acid (IAA) treatments, seven StDDP genes' expressions indicated that ABA and IAA performed important roles in immunity response. The expression profiling and real-time qPCR of stems, roots, and leaves revealed StDDPs' significant role in growth and development. These expression results of DDPs are primary functional analysis and present basic information for other economically important crops.

Genetic Dissection of Mature Root Characteristics by Genome-Wide Association Studies in Rapeseed (Brassica napus L.)

Roots are complicated quantitative characteristics that play an essential role in absorbing water and nutrients. To uncover the genetic variations for root-related traits in rapeseed, twelve mature root traits of a Brassica napus association panel were investigated in the field within three environments. All traits showed significant phenotypic variation among genotypes, with heritabilities ranging from 55.18% to 79.68%. Genome-wide association studies (GWAS) using 20,131 SNPs discovered 172 marker-trait associations, including 103 significant SNPs (-log10 (p) > 4.30) that explained 5.24-20.31% of the phenotypic variance. With the linkage disequilibrium r² > 0.2, these significant associations were binned into 40 quantitative trait loci (QTL) clusters. Among them, 14 important QTL clusters were discovered in two environments and/or with phenotypic contributions greater than 10%. By analyzing the genomic regions within 100 kb upstream and downstream of the peak SNPs within the 14 loci, 334 annotated genes were found. Among these, 32 genes were potentially associated with root development according to their expression analysis. Furthermore, the protein interaction network using the 334 annotated genes gave nine genes involved in a substantial number of interactions, including a key gene associated with root development, BnaC09g36350D. This research provides the groundwork for deciphering B. napus' genetic variations and improving its root system architecture.