Genome-wide analysis of Dof transcription factors reveals functional characteristics during development and response to biotic stresses in pepper

Exploring CDF gene family in wild potato under salinity stress unveils promising candidates for developing climate-resilient crops

The DNA-binding with one finger (Dof) gene family is a class of plant-specific transcription factors involved in diverse biological processes, including response to biotic and abiotic stresses. Members of this family have been reported in the cultivated potato Solanum tuberosum, but clues to the roles of several Dof genes are still lacking. Potato wild relatives represent a genetic reservoir for breeding as they could provide useful alleles for adaptation to the environment and tolerance to biotic and abiotic stresses. We performed an in silico analysis to identify genes belonging to the Dof family in the wild potato S. commersonii, confirming that the identified Dof genes can be grouped in four classes (A, B, C, D), as reported for cultivated potato. A special focus was dedicated to Cycling Dof Factors (CDFs), which play a crucial role in plant responses to abiotic stresses. Analysis of available RNA-seq data confirmed CDF genes as regulated by stresses and often in a tissue specific manner. To ascertain their involvement in the stress response, S. tuberosum and S. commersonii plantlets growing in vitro were subjected to salt stress (80mM NaCl) for short (2 days) and prolonged (7 days) times. Analysis of phenotypic traits and qRT-PCR expression profiles of target CDF genes in aerial and root tissues showed differences between the two species. In addition, after saline treatment, changes in total phenols, proline, and malondialdehyde suggested a diverse perception of saline stress in S. commersonii vs. S. tuberosum. Overall, this study provided useful clues to the involvement of CDF genes in salt response and promoted the identification of potential candidate genes for further functional studies.

Genome-wide identification and molecular evolution of Dof gene family in Camellia oleifera

DNA binding with one finger(Dof) gene family is a class of transcription factors which play an important role on plant growth and development. Genome-wide identification results indicated that there were 45 Dof genes(ColDof) in C.oleifera genome. All 45 ColDof proteins were non-transmembrane and non-secretory proteins. Phosphorylation site analysis showed that biological function of ColDof proteins were mainly realized by phosphorylation at serine (Ser) site. The secondary structure of 44 ColDof proteins was dominated by random coil, and only one ColDof protein was dominated by α-helix. ColDof genes' promoter region contained a variety of cis-acting elements, including light responsive regulators, gibberellin responsive regulators, abscisic acid responsive regulators, auxin responsive regulators and drought induction responsive regulators. The SSR sites analysis showed that the proportion of single nucleotide repeats and the frequency of A/T in ColDof genes were the largest. Non-coding RNA analysis showed that 45 ColDof genes contained 232 miRNAs. Transcription factor binding sites of ColDof genes showed that ColDof genes had 5793 ERF binding sites, 4381 Dof binding sites, 2206 MYB binding sites, 3702 BCR-BPC binding sites. ColDof9, ColDof39 and ColDof44 were expected to have the most TFBSs. The collinearity analysis showed that there were 40 colinear locis between ColDof proteins and AtDof proteins. Phylogenetic analysis showed that ColDof gene family was most closely related to that of Camellia sinensis var. sinensis cv.Biyun and Camellia lanceoleosa. Protein-protein interaction analysis showed that ColDof34, ColDof20, ColDof28, ColDof35, ColDof42 and ColDof26 had the most protein interactions. The transcriptome analysis of C. oleifera seeds showed that 21 ColDof genes were involved in the growth and development process of C. oleifera seeds, and were expressed in 221 C. oleifera varieties. The results of qRT-PCR experiments treated with different concentrations NaCl and PEG6000 solutions indicated that ColDof1, ColDof2, ColDof14 and ColDof36 not only had significant molecular mechanisms for salt stress tolerance, but also significant molecular functions for drought stress tolerance in C. oleifera. The results of this study provide a reference for further understanding of the function of ColDof genes in C.oleifera.

Identification and functional analysis of the DOF gene family in Populus simonii: implications for development and stress response

Background

Populus simonii, a notable native tree species in northern China, demonstrates impressive resistance to stress, broad adaptability, and exceptional hybridization potential. DOF family is a class of specific transcription factors that only exist in plants, widely participating in plant growth and development, and also playing an important role in abiotic stress response. To date, there have been no reported studies on the DOF gene family in P. simonii, and the expression levels of this gene family in different tissues of poplar, as well as its expression patterns under cold, heat, and other stress conditions, remain unclear.

Methods

In this study, DOF gene family were identified from the P. simonii genome, and various bioinformatics data on the DOF gene family, gene structure, gene distribution, promoters and regulatory networks were analyzed. Quantitative real time PCR technology was used to verify the expression patterns of the DOF gene family in different poplar tissues.

Results

This research initially pinpointed 41 PSDOF genes in P. simonii genome. The chromosomal localization results revealed that the PSDOF genes is unevenly distributed among 19 chromosomes, with the highest number of genes located on chromosomes 4, 5, and 11. A phylogenetic tree was constructed based on the homology between Arabidopsis thaliana and P. simonii, dividing the 41 PSDOF genes into seven subgroups. The expression patterns of PSDOF genes indicated that specific genes are consistently upregulated in various tissues and under different stress conditions, suggesting their pivotal involvement in both plant development and response to stress. Notably, PSDOF35 and PSDOF28 serve as pivotal hubs in the interaction network, playing a unique role in coordinating with other genes within the family.

Conclusion

The analysis enhances our comprehension of the functions of the DOF gene family in tissue development and stress responses within P. simonii. These findings provide a foundation for future exploration into the biological roles of DOF gene family.

Identification of Dof transcription factors in Dendrobium huoshanense and expression pattern under abiotic stresses

Introduction: DNA-binding with one finger (Dof) transcription factors (TFs) are a unique family of TFs found in higher plants that regulate plant responses to light, hormones, and abiotic stresses. The specific involvement of Dof genes in the response to environmental stresses remains unknown in D. huoshanense. Methods: A total of 22 Dof family genes were identified from the D. huoshanense genome. Results: Chromosome location analysis showed that DhDof genes were distributed on 12 chromosomes, with the largest number of Dof genes located on chromosome 8. The phylogenetic tree revealed that DhDofs could be categorized into 11 distinct subgroups. In addition to the common groups, DhDof4, DhDof5, DhDof17, and the AtDof1.4 ortholog were clustered into the B3 subgroup. Group E was a newly identified branch, among which DhDof6, DhDof7, DhDof8, and DhDof9 were in an independent branch. The conserved motifs and gene structure revealed the differences in motif number and composition of DhDofs. The dof domain near the N-terminus was highly conserved and contained a C2-C2-type zinc finger structure linked with four cysteines. Microsynteny and interspecies collinearity revealed gene duplication events and phylogenetic tree among DhDofs. Large-scale gene duplication had not occurred among the DhDofs genes and only in one pair of genes on chromosome 13. Synteny blocks were found more often between D. huoshanense and its relatives and less often between Oryza sativa and Arabidopsis thaliana. Selection pressure analysis indicated that DhDof genes were subject to purifying selection. Expression profiles and correlation analyses revealed that the Dof gene under hormone treatments showed several different expression patterns. DhDof20 and DhDof21 had the highest expression levels and were co-expressed under MeJA induction. The cis-acting element analysis revealed that each DhDof had several regulatory elements involved in plant growth as well as abiotic stresses. qRT-PCR analysis demonstrated that DhDof2 was the main ABA-responsive gene and DhDof7 was the main cold stress-related gene. IAA suppressed the expression of some Dof candidates, and SA inhibited most of the candidate genes. Discussion: Our results may provide new insights for the further investigation of the Dof genes and the screening of the core stress-resistance genes.

Identification and expression characteristics of NLP (NIN-like protein) gene family in pepper (Capsicum annuum L.)

BACKGROUND

Pepper (Capsicum annum L.) is the main crop in the vegetable industry. The growth and development of peppers are regulated by nitrate, but there is limited research on the molecular mechanisms of nitrate absorption and assimilation in peppers. A plant specific transcription factor NLP plays an important role in nitrate signal transduction.

METHODS AND RESULTS

In this study, a total of 7 NLP members were identified based on pepper genome data. Two nitrogen transport elements (GCN4) were found in the CaNLP5 promoter. In the phylogenetic tree, CaNLP members are divided into three branches, with pepper NLP and tomato NLP having the closest genetic relationship. The expression levels of CaNLP1, CaNLP3, and CaNLP4 are relatively high in the roots, stems, and leaves. The expression level of CaNLP7 gene is relatively high during the 5-7 days of pepper fruit color transformation. After various non-Biotic stress and hormone treatments, the expression of CaNLP1 was at a high level. The expression of CaNLP3 and CaNLP4 was down regulated in leaves, but up regulated in roots. Under conditions of nitrogen deficiency and sufficient nitrate, the expression patterns of NLP genes in pepper leaves and roots were determined.

CONCLUSION

These results provide important insights into the multiple functions of CaNLPs in regulating nitrate absorption and transport.

Genome-Wide Analysis of Barley bHLH Transcription Factors and the Functional Characterization of HvbHLH56 in Low Nitrogen Tolerance in Arabidopsis

Improvement of low nitrogen (LN) tolerance or nitrogen use efficiency (NUE) in crops is imperative for environment-friendly agriculture development. The basic helix-loop-helix (bHLH) transcription factors are involved in multiple abiotic stresses and are suitable as candidate genes for improving LN tolerance. Few studies were performed on the characterization of the HvbHLH gene family and their function in response to LN stress in barley. In this study, 103 HvbHLH genes were identified through genome-wide analysis. HvbHLH proteins were classified into 20 subfamilies based on phylogenetic analysis in barley, which was supported by conserved motifs and gene structure analysis. The stress-related cis-element analysis in the promoters showed that HvbHLHs are probably involved in multiple stress responses. By phylogenetic analysis of HvbHLHs and bHLHs in other plants, some HvbHLHs were predicted to play roles in response to nutrition deficiency stress. Furthermore, at least 16 HvbHLHs were differentially expressed in two barley genotypes differing in LN tolerance under LN stress. Finally, overexpression of HvbHLH56 enhanced LN stress tolerance in transgenic Arabidopsis, suggesting it is an important regulator in LN stress response. The differentially expressed HvbHLHs identified herein may be valuable for the breeding of barley cultivars with LN tolerance.

Expression of RsPORB Is Associated with Radish Root Color

Radish (Raphanus sativus) plants exhibit varied root colors due to the accumulation of chlorophylls and anthocyanins compounds that are beneficial for both human health and visual quality. The mechanisms of chlorophyll biosynthesis have been extensively studied in foliar tissues but remain largely unknown in other tissues. In this study, we examined the role of NADPH:protochlorophyllide oxidoreductases (PORs), which are key enzymes in chlorophyll biosynthesis, in radish roots. The transcript level of RsPORB was abundantly expressed in green roots and positively correlated with chlorophyll content in radish roots. Sequences of the RsPORB coding region were identical between white (948) and green (847) radish breeding lines. Additionally, virus-induced gene silencing assay with RsPORB exhibited reduced chlorophyll contents, verifying that RsPORB is a functional enzyme for chlorophyll biosynthesis. Sequence comparison of RsPORB promoters from white and green radishes showed several insertions and deletions (InDels) and single-nucleotide polymorphisms. Promoter activation assays using radish root protoplasts verified that InDels of the RsPORB promoter contribute to its expression level. These results suggested that RsPORB is one of the key genes underlying chlorophyll biosynthesis and green coloration in non-foliar tissues, such as roots.

Identification and characterization of Dof genes in Cerasus humilis

Introduction

Dof genes encode plant-specific transcription factors, which regulate various biological processes such as growth, development, and secondary metabolite accumulation.

Methods

We conducted whole-genome analysis of Chinese dwarf cherry (Cerasus humilis) to identify ChDof genes and characterize the structure, motif composition, cis-acting elements, chromosomal distribution, and collinearity of these genes as well as the physical and chemical properties, amino acid sequences, and phylogenetic evolution of the encoded proteins.

Results

The results revealed the presence of 25 ChDof genes in C. humilis genome. All 25 ChDof genes could be divided into eight groups, and the members of the same group had similar motif arrangement and intron-exon structure. Promoter analysis showed that cis-acting elements responsive to abscisic acid, low temperature stress, and light were dominant. Transcriptome data revealed that most ChDof genes exhibited tissue-specific expression. Then, we performed by qRT-PCR to analyze the expression patterns of all 25 ChDof genes in fruit during storage. The results indicated that these genes exhibited different expression patterns, suggesting that they played an important role in fruit storage.

Discussion

The results of this study provide a basis for further investigation of the biological function of Dof genes in C. humilis fruit.

Identification of Photoperiod- and Phytohormone-Responsive DNA-Binding One Zinc Finger (Dof) Transcription Factors in Akebia trifoliata via Genome-Wide Expression Analysis

As a kind of plant-specific transcription factor (TF), DNA-Binding One Zinc Finger (Dof) is widely involved in the response to environmental change, and as an evolutionarily important perennial plant species, Akebia trifoliata is ideal for studying environmental adaptation. In this study, a total of 41 AktDofs were identified in the A. trifoliata genome. First, the characteristics, including the length, exon number, and chromosomal distribution of the AktDofs and the isoelectric point (PI), amino acid number, molecular weight (MW), and conserved motifs of their putative proteins, were reported. Second, we found that all AktDofs evolutionarily underwent strong purifying selection, and many (33, 80.5%) of them were generated by whole-genome duplication (WGD). Third, we outlined their expression profiles by the use of available transcriptomic data and RT-qPCR analysis. Finally, we identified four candidate genes (AktDof21, AktDof20, AktDof36, and AktDof17) and three other candidate genes (AktDof26, AktDof16, and AktDof12) that respond to long day (LD) and darkness, respectively, and that are closely associated with phytohormone-regulating pathways. Overall, this research is the first to identify and characterize the AktDofs family and is very helpful for further research on A. trifoliata adaptation to environmental factors, especially photoperiod changes.

Genome-wide survey and expression analysis of Dof transcription factor family in sweetpotato shed light on their promising functions in stress tolerance

DNA-binding with one finger (Dof) transcription factors play a crucial role in plant abiotic stress regulatory networks, although massive Dofs have been systematically characterized in plants, they have not been identified in the hexaploid crop sweetpotato. Herein, 43 IbDof genes were detected to be disproportionally dispersed across 14 of the 15 chromosomes of sweetpotato, and segmental duplications were discovered to be the major driving force for the expansion of IbDofs. The collinearity analysis of IbDofs with their related orthologs from eight plants revealed the potential evolutionary history of Dof gene family. Phylogenetic analysis displayed that IbDof proteins were assigned into nine subfamilies, and the regularity of gene structures and conserved motifs was consistent with the subgroup classification. Additionally, five chosen IbDof genes were shown to be substantially and variably induced under various abiotic conditions (salt, drought, heat, and cold), as well as hormone treatments (ABA and SA), according to their transcriptome data and qRT-PCR experiments. Consistently, the promoters of IbDofs contained a number of cis-acting elements associated with hormone and stress responses. Besides, it was noted that IbDof2 had transactivation activity in yeasts, while IbDof-11/-16/-36 did not, and protein interaction network analysis and yeast two-hybrid experiments revealed a complicated interaction connection amongst IbDofs. Collectively, these data lay a foundation for further functional explorations of IbDof genes, especially with regards to the possible application of multiple IbDof members in breeding the tolerant plants.

Genome-Wide Association Study (GWAS) of White Mold Resistance in Snap Bean

White mold can result in snap bean yield losses of 90 to 100% when field conditions favor the pathogen. A genome-wide association study (GWAS) was conducted to detect loci significantly associated with white mold resistance in a panel of snap bean (Phaseolus vulgaris L.) cultivars. Two populations of snap bean were used in this study. The first population was the BeanCAP (Coordinated Agriculture Project) Snap Bean Diversity Panel (SBDP) (n = 136), and the second population was the Snap Bean Association Panel (SnAP) (n = 378). SBDP was evaluated for white mold reaction in the field in 2012 and 2013, and SnAP was screened in a greenhouse only using the seedling straw test in 2016. Two reference genomes representing the Andean and Middle American centers of domestication were utilized to align the genotyping-by-sequencing (GBS) data. A GWAS was performed using FarmCPU with one principal component after comparing five models. Thirty-four single-nucleotide polymorphisms (SNPs) significantly associated with white mold resistance were detected. Eleven significant SNPs were identified by the seedling straw test, and 23 significant SNPs were identified by field data. Fifteen SNPs were identified within a 100 kb window containing pentatricopeptide repeat (PPR)-encoding genes, and eleven were close to leucine-rich repeat (LRR)-encoding genes, suggesting that these two classes are of outsized importance for snap bean resistance to white mold.

Genome-Wide Identification of the DOF Gene Family Involved in Fruitlet Abscission in Areca catechu L

Fruitlet abscission frequently occurs in Areca catechu L. and causes considerable production loss. However, the inducement mechanism of fruitlet abscission remains mysterious. In this study, we observed that the cell architecture in the abscission zone (AZ) was distinct with surrounding tissues, and varied obviously before and after abscission. Transcriptome analysis of the “about-to-abscise” and “non-abscised” AZs were performed in A. catechu, and the genes encoding the plant-specific DOF (DNA-binding with one finger) transcription factors showed a uniform up-regulation in AZ, suggesting a role of the DOF transcription in A. catechu fruitlet abscission. In total, 36 members of the DOF gene family distributed in 13 chromosomes were identified from the A. catechu genome. The 36 AcDOF genes were classified into nine subgroups based on phylogenic analysis. Six of them showed an AZ-specific expression pattern, and their expression levels varied according to the abscission process. In total, nine types of phytohormone response cis-elements and five types of abiotic stress related cis-elements were identified in the promoter regions of the AcDOF genes. In addition, histochemical staining showed that lignin accumulation of vascular bundles in AZ was significantly lower than that in pedicel and mesocarp, indicating the specific characteristics of the cell architecture in AZ. Our data suggests that the DOF transcription factors might play a role in fruitlet abscission regulation in A. catechu.

Root Cap to Soil Interface: A Driving Force Toward Plant Adaptation and Development

Land plants have developed robust roots to grow in diverse soil ecosystems. The distal end of the root tip has a specialized organ called the 'root cap'. The root cap assists the roots in penetrating the ground, absorbing water and minerals, avoiding heavy metals and regulating the rhizosphere microbiota. Furthermore, root-cap-derived auxin governs the lateral root patterning and directs root growth under varying soil conditions. The root cap formation is hypothesized as one of the key innovations during root evolution. Morphologically diversified root caps in early land plant lineage and later in angiosperms aid in improving the adaptation of roots and, thereby, plants in diverse soil environments. This review article presents a retrospective view of the root cap's important morphological and physiological characteristics for the root-soil interaction and their response toward various abiotic and biotic stimuli. Recent single-cell RNAseq data shed light on root cap cell-type-enriched genes. We compiled root cap cell-type-enriched genes from Arabidopsis, rice, maize and tomato and analyzed their transcription factor (TF) binding site enrichment. Further, the putative gene regulatory networks derived from root-cap-enriched genes and their TF regulators highlight the species-specific biological functions of root cap genes across the four plant species.

Identification of the 12-oxo-phytoeienoic acid reductase (OPR) gene family in pepper (Capsicum annuum L.) and functional characterization of CaOPR6 in pepper fruit development and stress response

The 12-oxo-phytoeienoic acid reductase (OPR) is a kind of enzyme in octadecanoid biosynthesis pathway, which determines the biosynthesis of jasmonic acid. Although the roles of OPRs have been extensively studied in several crop plants, little is known about the biological functions of OPR encoding genes in Capsicum annuum plants. In this study, seven OPR family genes (CaOPR1-7) were identified from the C. annuum genome. The physical and chemical properties of CaOPR1-7 were further analyzed, including gene expression patterns, promoter elements and chromosomal locations. The results showed that the seven CaOPR homologous could be divided into two subgroups, and CaOPR6 was highly similar to AtOPR3 in Arabidopsis. The expression of CaOPR6 was significantly induced by various stresses such as cold, salt and pathogen infection, indicating that CaOPR6 plays important roles in response to abiotic and biotic stresses. Overall, these findings improve the understanding of the biological functions of CaOPR6 in the development of pepper fruit and stress response of pepper plants, and facilitate further studies on the molecular biology of OPR proteins in Solanaceae vegetables.

Genome-Wide Identification and Expression Analysis of Dof Transcription Factors in Lotus (Nelumbo nucifera Gaertn.)

Lotus (Nelumbo nucifera Gaertn.) is a traditional Chinese aquatic flower with high ornamental and economic value, but water salinity seriously affects lotus cultivation and distribution. The Dof transcription factors (TFs) play a crucial function in the regulatory network of growth and defense in plants. However, no systematic investigations of the Dof TFs in lotus have been performed. In this study, comprehensive searches of the lotus genome yielded 29 potential NnDofs. We carried out a series of standardized analyses, which include physical properties, multiple sequence alignment, phylogenetic analysis, gene structure, motif composition, cis-acting element prediction, chromosome distribution, and synteny analysis. The results showed that segment duplication probably caused the NnDofs gene family expansion. The potential functions of NnDofs in lotus development and stress conditions are speculated by promoter analysis. Furthermore, a complete expression investigation of NnDofs utilizing an RNA-seq atlas and quantitative real-time polymerase chain reaction (qRT-PCR) was performed. The majority of the NnDofs exhibit tissue-specific expression patterns, and many genes have been identified as being extremely sensitive to salt stressors. Overall, this study is the first to report a genome-wide assessment of the Dof family in lotus, and the findings offer vital insights for prospective functional studies on lotus salinity stress.

GBS-Based SNP Map Pinpoints the QTL Associated With Sorghum Downy Mildew Resistance in Maize (Zea mays L.)

Sorghum downy mildew (SDM), caused by the biotrophic fungi Peronosclerospora sorghi, threatens maize production worldwide, including India. To identify quantitative trait loci (QTL) associated with resistance to SDM, we used a recombinant inbred line (RIL) population derived from a cross between resistant inbred line UMI936 (w) and susceptible inbred line UMI79. The RIL population was phenotyped for SDM resistance in three environments [E1-field (Coimbatore), E2-greenhouse (Coimbatore), and E3-field (Mandya)] and also utilized to construct the genetic linkage map by genotyping by sequencing (GBS) approach. The map comprises 1516 SNP markers in 10 linkage groups (LGs) with a total length of 6924.7 cM and an average marker distance of 4.57 cM. The QTL analysis with the phenotype and marker data detected nine QTL on chromosome 1, 2, 3, 5, 6, and 7 across three environments. Of these, QTL namely qDMR1.2, qDMR3.1, qDMR5.1, and qDMR6.1 were notable due to their high phenotypic variance. qDMR3.1 from chromosome 3 was detected in more than one environment (E1 and E2), explaining the 10.3% and 13.1% phenotypic variance. Three QTL, qDMR1.2, qDMR5.1, and qDMR6.1 from chromosomes 1, 5, and 6 were identified in either E1 or E3, explaining 15.2%–18% phenotypic variance. Moreover, genome mining on three QTL (qDMR3.1, qDMR5.1, and qDMR6.1) reveals the putative candidate genes related to SDM resistance. The information generated in this study will be helpful for map-based cloning and marker-assisted selection in maize breeding programs.

Transcriptomics Advancement in the Complex Response of Plants to Viroid Infection

Viroids are the smallest plant pathogens, consisting of a single-stranded circular RNA of less than 500 ribonucleotides in length. Despite their noncoding nature, viroids elicit disease symptoms in many economically important plant hosts, and are, thus, a class of pathogens of great interest. How these viroids establish disease within host plants, however, is not yet fully understood. Recent transcriptomic studies have revealed that viroid infection influences the expression of genes in several pathways and processes in plants, including defence responses, phytohormone signalling, cell wall modification, photosynthesis, secondary metabolism, transport, gene expression and protein modification. There is much debate about whether affected pathways signify a plant response to viroid infection, or are associated with the appearance of disease symptoms in these interactions. In this review, we consolidate the findings of viroid–host transcriptome studies to provide an overview of trends observed in the data. When considered together, changes in the gene expression of different hosts upon viroid infection reveal commonalities and differences in diverse interactions. Here, we discuss whether trends in host gene expression can be correlated to plant defence or disease development during viroid infection, and highlight avenues for future research in this field.

Emerging Roles of Plant DNA-Binding With One Finger Transcription Factors in Various Hormone and Stress Signaling Pathways

Coordinated transcriptional regulation of stress-responsive genes orchestrated by a complex network of transcription factors (TFs) and the reprogramming of metabolism ensure a plant's continued growth and survival under adverse environmental conditions (e.g., abiotic stress). DNA-binding with one finger (Dof) proteins, a group of plant-specific TF, were identified as one of several key components of the transcriptional regulatory network involved in abiotic stress responses. In many plant species, Dofs are often activated in response to a wide range of adverse environmental conditions. Dofs play central roles in stress tolerance by regulating the expression of stress-responsive genes via the DOFCORE element or by interacting with other regulatory proteins. Moreover, Dofs act as a key regulatory hub of several phytohormone pathways, integrating abscisic acid, jasmonate, SA and redox signaling in response to many abiotic stresses. Taken together, we highlight a unique role of Dofs in hormone and stress signaling that integrates plant response to adverse environmental conditions with different aspects of plant growth and development.

Magnesium- a Forgotten Element: Phenotypic Variation and Genome Wide Association Study in Turkish Common Bean Germplasm

Magnesium (Mg) is the fourth most abundant element in the human body and plays the role of cofactor for more than 300 enzymatic reactions. In plants, Mg is involved in various key physiological and biochemical processes like growth, development, photophosphorylation, chlorophyll formation, protein synthesis, and resistance to biotic and abiotic stresses. Keeping in view the importance of this element, the present investigation aimed to explore the Mg contents diversity in the seeds of Turkish common bean germplasm and to identify the genomic regions associated with this element. A total of 183 common bean accessions collected from 19 provinces of Turkey were used as plant material. Field experiments were conducted according to an augmented block design during 2018 in two provinces of Turkey, and six commercial cultivars were used as a control group. Analysis of variance depicted that Mg concentration among common bean accessions was statistically significant (p < 0.05) within each environment, however genotype × environment interaction was non-significant. A moderate level (0.60) of heritability was found in this study. Overall mean Mg contents for both environments varied from 0.33 for Nigde-Dermasyon to 1.52 mg kg⁻¹ for Nigde-Derinkuyu landraces, while gross mean Mg contents were 0.92 mg kg⁻¹. At the province level, landraces from Bolu were rich while the landraces from Bitlis were poor in seed Mg contents respectively. The cluster constellation plot divided the studied germplasm into two populations on the basis of their Mg contents. Marker-trait association was performed using a mixed linear model (Q + K) with a total of 7,900 DArTseq markers. A total of six markers present on various chromosomes (two at Pv01, and one marker at each chromosome i.e., Pv03, Pv07, Pv08, Pv11) showed statistically significant association for seed Mg contents. Among these identified markers, the DArT-3367607 marker present on chromosome Pv03 contributed to maximum phenotypic variation (7.5%). Additionally, this marker was found within a narrow region of previously reported markers. We are confident that the results of this study will contribute significantly to start common bean breeding activities using marker assisted selection regarding improved Mg contents.

Universal gene co-expression network reveals receptor-like protein genes involved in broad-spectrum resistance in pepper (Capsicum annuum L.)

Receptor-like proteins (RLPs) on plant cells have been implicated in immune responses and developmental processes. Although hundreds of RLP genes have been identified in plants, only a few RLPs have been functionally characterized in a limited number of plant species. Here, we identified RLPs in the pepper (Capsicum annuum) genome and performed comparative transcriptomics coupled with the analysis of conserved gene co-expression networks (GCNs) to reveal the role of core RLP regulators in pepper-pathogen interactions. A total of 102 RNA-seq datasets of pepper plants infected with four pathogens were used to construct CaRLP-targeted GCNs (CaRLP-GCNs). Resistance-responsive CaRLP-GCNs were merged to construct a universal GCN. Fourteen hub CaRLPs, tightly connected with defense-related gene clusters, were identified in eight modules. Based on the CaRLP-GCNs, we evaluated whether hub CaRLPs in the universal GCN are involved in the biotic stress response. Of the nine hub CaRLPs tested by virus-induced gene silencing, three genes (CaRLP264, CaRLP277, and CaRLP351) showed defense suppression with less hypersensitive response-like cell death in race-specific and non-host resistance response to viruses and bacteria, respectively, and consistently enhanced susceptibility to Ralstonia solanacearum and/or Phytophthora capsici. These data suggest that key CaRLPs are involved in the defense response to multiple biotic stresses and can be used to engineer a plant with broad-spectrum resistance. Together, our data show that generating a universal GCN using comprehensive transcriptome datasets can provide important clues to uncover genes involved in various biological processes.

Metabolic and Transcriptomic Profiling of Lilium Leaves Infected With Botrytis elliptica Reveals Different Stages of Plant Defense Mechanisms

Botrytis elliptica, the causal agent of gray mold disease, poses a major threat to commercial Lilium production, limiting its ornamental value and yield. The molecular and metabolic regulation mechanisms of Lilium's defense response to B. elliptica infection have not been completely elucidated. Here, we performed transcriptomic and metabolomic analyses of B. elliptica resistant Lilium oriental hybrid "Sorbonne" to understand the molecular basis of gray mold disease resistance in gray mold disease. A total of 115 differentially accumulated metabolites (DAMs) were detected by comparing the different temporal stages of pathogen infection. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed the differentially expressed genes (DEGs) and DAMs were enriched in the phenylpropanoid and flavonoid pathways at all stages of infection, demonstrating the prominence of these pathways in the defense response of "Sorbonne" to B. elliptica. Network analysis revealed high interconnectivity of the induced defense response. Furthermore, time-course analysis of the transcriptome and a weighted gene coexpression network analysis (WGCNA) led to the identification of a number of hub genes at different stages, revealing that jasmonic acid (JA), salicylic acid (SA), brassinolide (BR), and calcium ions (Ca2+) play a crucial role in the response of "Sorbonne" to fungal infection. Our work provides a comprehensive perspective on the defense response of Lilium to B. elliptica infection, along with a potential transcriptional regulatory network underlying the defense response, thereby offering gene candidates for resistance breeding and metabolic engineering of Lilium.

The Cowpea Kinome: Genomic and Transcriptomic Analysis Under Biotic and Abiotic Stresses

The present work represents a pioneering effort, being the first to analyze genomic and transcriptomic data from Vigna unguiculata (cowpea) kinases. We evaluated the cowpea kinome considering its genome-wide distribution and structural characteristics (at the gene and protein levels), sequence evolution, conservation among Viridiplantae species, and gene expression in three cowpea genotypes under different stress situations, including biotic (injury followed by virus inoculation-CABMV or CPSMV) and abiotic (root dehydration). The structural features of cowpea kinases (VuPKs) indicated that 1,293 bona fide VuPKs covered 20 groups and 118 different families. The RLK-Pelle was the largest group, with 908 members. Insights on the mechanisms of VuPK genomic expansion and conservation among Viridiplantae species indicated dispersed and tandem duplications as major forces for VuPKs' distribution pattern and high orthology indexes and synteny with other legume species, respectively. K _a /K _s ratios showed that almost all (91%) of the tandem duplication events were under purifying selection. Candidate cis-regulatory elements were associated with different transcription factors (TFs) in the promoter regions of the RLK-Pelle group. C2H2 TFs were closely associated with the promoter regions of almost all scrutinized families for the mentioned group. At the transcriptional level, it was suggested that VuPK up-regulation was stress, genotype, or tissue dependent (or a combination of them). The most prominent families in responding (up-regulation) to all the analyzed stresses were RLK-Pelle_DLSV and CAMK_CAMKL-CHK1. Concerning root dehydration, it was suggested that the up-regulated VuPKs are associated with ABA hormone signaling, auxin hormone transport, and potassium ion metabolism. Additionally, up-regulated VuPKs under root dehydration potentially assist in a critical physiological strategy of the studied cowpea genotype in this assay, with activation of defense mechanisms against biotic stress while responding to root dehydration. This study provides the foundation for further studies on the evolution and molecular function of VuPKs.

Importance of inositols and their derivatives in cowpea under root dehydration: An omics perspective

This work presents a robust analysis of the inositols (INSs) and raffinose family oligosaccharides (RFOs) pathways, using genomic and transcriptomic tools in cowpea under root dehydration. Nineteen (~70%) of the 26 scrutinized enzymes presented transcriptional up-regulation in at least one treatment time. The transcriptional orchestration allowed categorization of the analyzed enzymes as time-independent (those showing the same regulation throughout the assay) and time-dependent (those showing different transcriptional regulation over time). It is suggested that up-regulated time-independent enzymes (INSs: myo-inositol oxygenase, inositol-tetrakisphosphate 1-kinase 3, phosphatidylinositol 4-phosphate 5-kinase 4-like, 1-phosphatidylinositol-3-phosphate 5-kinase, phosphoinositide phospholipase C, and non-specific phospholipase C; RFOs: α-galactosidase, invertase, and raffinose synthase) actively participate in the reorganization of cowpea molecular physiology under the applied stress. In turn, time-dependent enzymes, especially those up-regulated in some of the treatment times (INSs: inositol-pentakisphosphate 2-kinase, phosphatidylinositol 4-kinase, phosphatidylinositol synthase, multiple inositol polyphosphate phosphatase 1, methylmalonate-semialdehyde dehydrogenase, triosephosphate isomerase, myo-inositol-3-phosphate synthase, phosphatidylinositol 3,4,5-trisphosphate 3-phosphatase and protein-tyrosine-phosphatase, and phosphatidylinositol 3-kinase; RFOs: galactinol synthase) seem to participate in fine-tuning of the molecular physiology, helping the cowpea plants to acclimatize under dehydration stress. Not all loci encoding the studied enzymes were expressed during the assay; most of the expressed ones exhibited a variable transcriptional profile in the different treatment times. Genes of the INSs and RFOs pathways showed high orthology with analyzed Phaseoleae members, suggesting a relevant role within this legume group. Regarding the promoter regions of INSs and RFOs genes, some bona fide cis-regulatory elements were identified in association with seven transcription factor families (AP2-EFR, Dof-type, MADS-box, bZIP, CPP, ZF-HD, and GATA-type). Members of INSs and RFOs pathways potentially participate in other processes regulated by these proteins in cowpea.

Tissue-Specific RNA-Seq Analysis and Identification of Receptor-Like Proteins Related to Plant Growth in Capsicum annuum

Receptor-like proteins (RLPs) are a gene family of cell surface receptors that are involved in plant growth, development, and disease resistance. In a recent study, 438 pepper RLP genes were identified in the Capsicum annuum genome (CaRLPs) and determined to be present in response to multiple biotic stresses. To further understand the role of CaRLPs in plant growth and development, we analyzed expression patterns of all CaRLPs from various pepper tissues and developmental stages using RNA-seq. Ten CaRLP genes were selected for further analysis according to transcript levels with hierarchical clustering. The selected CaRLP genes displayed similarity of motifs within the same groups and structures typical of RLPs. To examine RLP function in growth and development, we performed loss-of-function analysis using a virus-induced gene silencing system. Three of the ten tested CaRLPs (CaRLP238, 253, and 360) in silenced plants exhibited phenotypic alteration with growth retardation compared to controls. All three gene-silenced peppers showed significant differences in root dry weight. Only CaRLP238 had significant differences in both root and shoot dry weight. Our results suggest that CaRLPs may play important roles in regulation of plant growth and development as well as function in defense responses to biotic stresses in the RLP gene family.

Genome-wide analysis and functional characterization of the Dof transcription factor family in rice (Oryza sativa L.)

Exhaustive searches of the rice genome have revealed 30 different potential OsDof (Oryza sativa DNA binding with One Finger) genes. Their subcellular localization, phylogenetic relationship, conserved motifs identification, chromosomal allocation, expression patterns, and interaction networks were analyzed. The Dof (DNA binding with One Finger) family of transcription factors represents a particular class of plant-specific transcriptional regulators, contain a highly conserved region of 50-52 amino acids (Dof domain) and involved in various plant developmental processes and response to various environmental stresses. Few (Oryza sativa) OsDof genes have been demonstrated previously for their biological functions but there is no comprehensive study on most of the Dof genes of rice. In the current study, exhaustive searches of the rice genome revealed 30 different potential OsDof genes, and then their subcellular localization, phylogenetic relationship, conserved motifs identification, chromosomal allocation, expression patterns, and interaction networks were analyzed. Phylogenetic analysis of Dof proteins in rice showed that they are distributed in 4 groups. By genome-wide observation of gene expression profiles, we found that OsDof genes showed significant variances in expression levels in different tissues across multiple developmental stages. Protein-protein correlation network analysis, shows a statically significant overlap of some OsDofs, suggesting their similar functions and a high degree of co-expression. The Dof family transcription factors have been reported for their involvement in the regulation of various gene expression processes in rice but still, most of the Dof genes are not characterized for their specific physiological functions. This study revealed useful information and clues about predicting the potential roles of OsDofs in rice by combining their genome-wide characterization, expression profiling, protein-protein interactions, and for further studies to develop high-quality rice varieties.

Genome-wide identification and characterization of basic helix-loop-helix genes in nine molluscs

The basic helix-loop-helix (bHLH) transcription factors form a large superfamily that plays an important role in numerous physiological processes, including development and response to environmental stresses. In this study, the distribution of bHLH genes in nine molluscs was systematically investigated (including five bivalves, three gastropods and one cephalopod). Finally, 53-85 bHLH genes were identified from each genome and classified into corresponding families by using phylogenetic analysis. The results of gene structure and conserved motif analysis illustrated the hereditary conservation of bHLH transcription factors during evolution but showed low similarity in group C. Through transcription profile analysis of C. gigas and T. granosa, we found a important role of bHLH genes in responding to multiple external challenges and development; meanwhile, they also exhibited tissue-specific expression. Interestingly, we were also surprised to find different bHLH genes from the same group generally possess similar patterns expression that tends to simultaneously present high or lower expression of multiple challenges and different tissues in this study. In summary, this study lays the foundation for further investigation of the biological functions and evolution of molluscan bHLH genes.

Genome-Wide Identification and Expression Analysis of the Dof Gene Family in Medicago sativa L. Under Various Abiotic Stresses

The Dof transcription factor is a plant-specific transcriptional regulator that plays important roles in plant development and acts as a mediator in plant external stress responses. However, Dofs have previously been identified in several plants but not in alfalfa (Medicago sativa L.), one of the most widely cultivated forage legumes. In the present study, a total of 40 MsDof genes were identified, and the phylogenetic reconstruction, classification, conserved motifs, and expression patterns under abscisic acid (ABA), cold, heat, drought and salt stresses of these Dof genes were comprehensively analyzed. The Dof genes family in alfalfa could be classified into eight classes. Gene ontology (GO) and tissue-specific analysis indicated that most MsDof genes may be involved in biological functions during plant growth. Moreover, the expression profiles and quantitative real-time PCR analysis indicated that eight candidate abiotic tolerance genes were induced in response to four abiotic stresses. This study identified the possibility of abiotic tolerance candidate genes playing various roles in stress resistance at the whole genome level, which would provide new information on the Dof family in alfalfa.

Capsidiol-related genes are highly expressed in response to Colletotrichum scovillei during Capsicum annuum fruit development stages

Capsicum annuum is one of the most important horticultural crops worldwide. Anthracnose disease (Colletotrichum spp.) is a major constraint for chili production, causing substantial losses. Capsidiol is a sesquiterpene phytoalexin present in pepper fruits that can enhance plant resistance. The genetic mechanisms involved in capisidiol biosynthesis are still poorly understood. In this study, a 3′ RNA sequencing approach was used to develop the transcriptional profile dataset of C. annuum genes in unripe (UF) and ripe fruits (RF) in response to C. scovillei infection. Results showed 4,845 upregulated and 4,720 downregulated genes in UF, and 2,560 upregulated and 1,762 downregulated genes in RF under fungus inoculation. Four capsidiol-related genes were selected for RT-qPCR analysis, two 5-epi-aristolochene synthase (CA12g05030, CA02g09520) and two 5-epi-aristolochene-1,3-dihydroxylase genes (CA12g05070, CA01g05990). CA12g05030 and CA01g05990 genes showed an early response to fungus infection in RF (24 h post-inoculation—HPI), being 68-fold and 53-fold more expressed at 96 HPI, respectively. In UF, all genes showed a late response, especially CA12g05030, which was 700-fold more expressed at 96 HPI compared to control plants. We are proving here the first high-throughput expression dataset of pepper fruits in response to anthracnose disease in order to contribute for future pepper breeding programs.

Identification of Callose Synthases in Stinging Nettle and Analysis of Their Expression in Different Tissues

Callose is an important biopolymer of β-1,3-linked glucose units involved in different phases of plant development, reproduction and response to external stimuli. It is synthesized by glycosyltransferases (GTs) known as callose synthases (CalS) belonging to family 48 in the Carbohydrate-Active enZymes (CAZymes) database. These GTs are anchored to the plasma membrane via transmembrane domains. Several genes encoding CalS have been characterized in higher plants with 12 reported in the model organism Arabidopsis thaliana. Recently, the de novo transcriptome of a fibre-producing clone of stinging nettle (Urtica dioica L.) was published and here it is mined for CalS genes with the aim of identifying members differentially expressed in the core and cortical tissues of the stem. The goal is to understand whether specific CalS genes are associated with distinct developmental stages of the stem internodes (elongation, thickening). Nine genes, eight of which encoding full-length CalS, are identified in stinging nettle. The phylogenetic analysis with CalS proteins from other fibre crops, namely textile hemp and flax, reveals grouping into 6 clades. The expression profiles in nettle tissues (roots, leaves, stem internodes sampled at different heights) reveal differences that are most noteworthy in roots vs. leaves. Two CalS are differentially expressed in the internodes sampled at the top and middle of the stem. Implications of their role in nettle stem tissue development are discussed.

Genome-Wide Identification, Structure Characterization, and Expression Profiling of Dof Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

DNA binding with one finger (Dof) proteins are plant-specific transcription factors with crucial roles in plant growth and stress response. Even so, little is known about them in wheat. In this study, 108 wheat Dof (TaDof) genes across 21 chromosomes were detected. Although variable in sequence length, molecular weight, and isoelectric point, all TaDof proteins contained conserved zinc-finger structures and were phylogenetically divided into 7 sub-groups. Exon/intron and motif analyses suggested that TaDof structures and conserved motifs were similar within sub-groups but diverse among sub-groups. Many segmental duplications were identified and Ka/Ks and inter-species synthetic analyses indicated that polyploidization was main reason for increased number of TaDofs. Prediction and experimental confirmation revealed that TaDofs functioned as transcription factors in the nucleus. Expression pattern profiling showed that TaDofs specifically affected growth and development, and biotic and abiotic stress responses. Wheat miRNAs and cis-regulator were predicted as essential players in molding TaDofs expression patterns. qRT-PCR analysis revealed that TaDofs were induced by salt and drought stresses. Customized annotation revealed that TaDofs were widely involved in phytohormone response, defense, growth and development, and metabolism. Our study provided a comprehensive understanding to wheat TaDofs.

Genome-wide characterization and expression analysis of the Dof gene family related to abiotic stress in watermelon

The plant DNA-binding with one finger (Dof) gene family is a class of plant-specific transcription factors that play vital roles in many biological processes and stress responses. In the present study, a total of 36 ClDof genes were identified in the watermelon genome, which were unevenly distributed on 10 chromosomes. Phylogenetic analysis showed that the ClDof proteins could be divided into nine groups, and the members in a particular group had similar motif arrangement and exon-intron structure. Synteny analysis indicated the presence of a large number of syntenic relationship events between watermelon and cucumber. In promoter analysis, five kinds of stress-related and nine kinds of hormone-related cis-elements were identified in the promoter regions of ClDof genes. We then analyzed the expression patterns of nine selected ClDof genes in eight specific tissues by qRT-PCR, and the results showed that they have tissue-specific expression patterns. We also evaluated the expression levels of 12 selected ClDof genes under salt stress and ABA treatments using qRT-PCR. As a result, they showed differential expression under these treatments, suggesting their important roles in stress response. Taken together, our results provide a basis for future research on the biological functions of Dof genes in watermelon.

Transcriptome profiling of abiotic responses to heat, cold, salt, and osmotic stress of Capsicum annuum L

Peppers (Capsicum annuum L.), belonging to the Solanaceae family, are one of the most economically important crops globally. Like other crops, peppers are threatened by diverse environmental conditions due to different pathogens and abiotic stresses. High-quality reference genomes with massive datasets of transcriptomes from various conditions can provide clues to preferred agronomic traits for breeding. However, few global gene expression profiling datasets have been published to examine the environmental stress-resistant mechanisms in peppers. In this study, we report the RNA-seq analyses of peppers treated with heat, cold, salinity, and osmotic stress at six different time points. RNA-seq libraries from 78 RNA samples containing three biological replicates per time point for each of the abiotic stresses and a mock control were constructed. A total of 204.68 Gb of transcriptome data were verified by differentially expressed genes and gene ontology enrichment analysis. Analyses of the transcriptome data in this study will provide useful information for basic studies of various stimuli to facilitate the development of stress-resistant pepper cultivars.

The Role of Major Transcription Factors in Solanaceous Food Crops under Different Stress Conditions: Current and Future Perspectives

Plant growth, development, and productivity are adversely affected by environmental stresses such as drought (osmotic stress), soil salinity, cold, oxidative stress, irradiation, and diverse diseases. These impacts are of increasing concern in light of climate change. Noticeably, plants have developed their adaptive mechanism to respond to environmental stresses by transcriptional activation of stress-responsive genes. Among the known transcription factors, DoF, WRKY, MYB, NAC, bZIP, ERF, ARF and HSF are those widely associated with abiotic and biotic stress response in plants. Genome-wide identification and characterization analyses of these transcription factors have been almost completed in major solanaceous food crops, emphasizing these transcription factor families which have much potential for the improvement of yield, stress tolerance, reducing marginal land and increase the water use efficiency of solanaceous crops in arid and semi-arid areas where plant demand more water. Most importantly, transcription factors are proteins that play a key role in improving crop yield under water-deficient areas and a place where the severity of pathogen is very high to withstand the ongoing climate change. Therefore, this review highlights the role of major transcription factors in solanaceous crops, current and future perspectives in improving the crop traits towards abiotic and biotic stress tolerance and beyond. We have tried to accentuate the importance of using genome editing molecular technologies like CRISPR/Cas9, Virus-induced gene silencing and some other methods to improve the plant potential in giving yield under unfavorable environmental conditions.

Knockdown of CaHSP60-6 confers enhanced sensitivity to heat stress in pepper (Capsicum annuum L.)

HSP60 gene family in pepper was analyzed through bioinformatics along with transcriptional regulation against multiple abiotic and hormonal stresses. Furthermore, the knockdown of CaHSP60-6 increased sensitivity to heat stress. The 60 kDa heat shock protein (HSP60) also known as chaperonin (cpn60) is encoded by multi-gene family that plays an important role in plant growth, development and in stress response as a molecular chaperone. However, little is known about the HSP60 gene family in pepper (Capsicum annuum L.). In this study, 16 putative pepper HSP60 genes were identified through bioinformatic tools. The phylogenetic tree revealed that eight of the pepper HSP60 genes (50%) clustered into group I, three (19%) into group II, and five (31%) into group III. Twelve (75%) CaHSP60 genes have more than 10 introns, while only a single gene contained no introns. Chromosomal mapping revealed that the tandem and segmental duplication events occurred in the process of evolution. Gene ontology enrichment analysis predicted that CaHSP60 genes were responsible for protein folding and refolding in an ATP-dependent manner in response to various stresses in the biological processes category. Multiple stress-related cis-regulatory elements were found in the promoter region of these CaHSP60 genes, which indicated that these genes were regulated in response to multiple stresses. Tissue-specific expression was studied under normal conditions and induced under 2 h of heat stress measured by RNA-Seq data and qRT-PCR in different tissues (roots, stems, leaves, and flowers). The data implied that HSP60 genes play a crucial role in pepper growth, development, and stress responses. Fifteen (93%) CaHSP60 genes were induced in both, thermo-sensitive B6 and thermo-tolerant R9 lines under heat treatment. The relative expression of nine representative CaHSP60 genes in response to other abiotic stresses (cold, NaCl, and mannitol) and hormonal applications [ABA, methyl jasmonate (MeJA), and salicylic acid (SA)] was also evaluated. Knockdown of CaHSP60-6 increased the sensitivity to heat shock treatment as documented by a higher relative electrolyte leakage, lipid peroxidation, and reactive oxygen species accumulation in silenced pepper plants along with a substantial lower chlorophyll content and antioxidant enzyme activity. These results suggested that HSP60 might act as a positive regulator in pepper defense against heat and other abiotic stresses. Our results provide a basis for further functional analysis of HSP60 genes in pepper.

Genome-Wide Analysis of Basic Helix-Loop-Helix Superfamily Members Reveals Organization and Chilling-Responsive Patterns in Cabbage (Brassica oleracea var. capitata L.)

Basic helix–loop–helix (bHLH) transcription factor (TF) family is commonly found in eukaryotes, which is one of the largest families of regulator proteins. It plays an important role in plant growth and development, as well as various biotic and abiotic stresses. However, a comprehensive analysis of the bHLH family has not been reported in Brassica oleracea. In this study, we systematically describe the BobHLHs in the phylogenetic relationships, expression patterns in different organs/tissues, and in response to chilling stress, and gene and protein characteristics. A total of 234 BobHLH genes were identified in the B. oleracea genome and were further clustered into twenty-three subfamilies based on the phylogenetic analyses. A large number of BobHLH genes were unevenly located on nine chromosomes of B. oleracea. Analysis of RNA-Seq expression profiles revealed that 21 BobHLH genes exhibited organ/tissue-specific expression. Additionally, the expression of six BobHLHs (BobHLH003, -048, -059, -093, -109, and -148) were significantly down-regulated in chilling-sensitive cabbage (CS-D9) and chilling-tolerant cabbage (CT-923). At 24 h chilling stress, BobHLH054 was significantly down-regulated and up-regulated in chilling-treated CS-D9 and CT-923. Conserved motif characterization and exon/intron structural patterns showed that BobHLH genes had similar structures in the same subfamily. This study provides a comprehensive analysis of BobHLH genes and reveals several candidate genes involved in chilling tolerance of B. oleracea, which may be helpful to clarify the roles of bHLH family members and understand the regulatory mechanisms of BobHLH genes in response to the chilling stress of cabbage.

Identification, expression, alternative splicing and functional analysis of pepper WRKY gene family in response to biotic and abiotic stresses

WRKY proteins are a large group of plant transcription factors that are involved in various biological processes, including biotic and abiotic stress responses, hormone response, plant development, and metabolism. WRKY proteins have been identified in several plants, but only a few have been identified in Capsicum annuum. Here, we identified a total of 62 WRKY genes in the latest pepper genome. These genes were classified into three groups (Groups 1–3) based on the structural features of their proteins. The structures of the encoded proteins, evolution, and expression under normal growth conditions were analyzed and 35 putative miRNA target sites were predicted in 20 CaWRKY genes. Moreover, the response to cold or CMV treatments of selected WRKY genes were examined to validate the roles under stresses. And alternative splicing (AS) events of some CaWRKYs were also identified under CMV infection. Promoter analysis confirmed that CaWRKY genes are involved in growth, development, and biotic or abiotic stress responses in hot pepper. The comprehensive analysis provides fundamental information for better understanding of the signaling pathways involved in the WRKY-mediated regulation of developmental processes, as well as biotic and abiotic stress responses.

Bioinformatical analysis and prediction of Nicotiana benthamiana bHLH transcription factors in Phytophthora parasitica resistance

The basic helix-loop-helix (bHLH) family, one of the largest transcription factor groups in plants, regulates many critical developmental processes. However, their functions in plant defense have not been extensively studied in Nicotiana benthamiana, an important model plant species for phytopathology. Here, we identified N. benthamiana bHLH genes (NbbHLHs) using a whole-genome searching approach, and found that the NbbHLHs are highly enriched and some subfamilies are selectively expanded in N. benthamiana. The results showed that gene duplication may be responsible for bHLH family expansion in this plant. Furthermore, we analyzed their expression profiles upon infection with Phytophthora parasitica. Finally, 28 candidate NbbHLHs may play important roles in Phytophthora pathogen resistance using cis-element analysis and protein-interaction network prediction. Taken together, our results established a platform for future studies of the gene family and provide molecular insights into plant immune responses against P. parasitica.

Transcriptomic analysis of cultivated cotton Gossypium hirsutum provides insights into host responses upon whitefly-mediated transmission of cotton leaf curl disease

Cotton is a commercial and economically important crop that generates billions of dollars in annual revenue worldwide. However, cotton yield is affected by a sap-sucking insect Bemisia tabaci (whitefly), and whitefly-borne cotton leaf curl disease (CLCuD). The causative agent of devastating CLCuD is led by the viruses belonging to the genus Begomovirus (family Geminiviridae), collectively called cotton leaf curl viruses. Unfortunately, the extensively cultivated cotton (Gossypium hirsutum) species are highly susceptible and vulnerable to CLCuD. Yet, the concomitant influence of whitefly and CLCuD on the susceptible G. hirsutum transcriptome has not been interpreted. In the present study we have employed an RNA Sequencing (RNA-Seq) transcriptomics approach to explore the differential gene expression in susceptible G. hirsutum variety upon infection with viruliferous whiteflies. Comparative RNA-Seq of control and CLCuD infected plants was done using Illumina HiSeq 2500. This study yielded 468 differentially expressed genes (DEGs). Among them, we identified 220 up and 248 downregulated DEGs involved in disease responses and pathogen defense. We selected ten genes for downstream RT-qPCR analyses on two cultivars, Karishma and MNH 786 that are susceptible to CLCuD. We observed a similar expression pattern of these genes in both susceptible cultivars that was also consistent with our transcriptome data further implying a wider application of our global transcription study on host susceptibility to CLCuD. We next performed weighted gene co-expression network analysis that revealed six modules. This analysis also identified highly co-expressed genes as well as 55 hub genes that co-express with ≥ 50 genes. Intriguingly, most of these hub genes are shown to be downregulated and enriched in cellular processes. Under-expression of such highly co-expressed genes suggests their roles in favoring the virus and enhancing plant susceptibility to CLCuD. We also discuss the potential mechanisms governing the establishment of disease susceptibility. Overall, our study provides a comprehensive differential gene expression analysis of G. hirsutum under whitefly-mediated CLCuD infection. This vital study will advance the understanding of simultaneous effect of whitefly and virus on their host and aid in identifying important G. hirsutum genes which intricate in its susceptibility to CLCuD.

Genome-wide Identification, Classification, and Expression Analysis of the Receptor-Like Protein Family in Tomato

Receptor-like proteins (RLPs) are involved in plant development and disease resistance. Only some of the RLPs in tomato (Solanum lycopersicum L.) have been functionally characterized though 176 genes encoding RLPs, which have been identified in the tomato genome. To further understand the role of RLPs in tomato, we performed genome-guided classification and transcriptome analysis of these genes. Phylogenic comparisons revealed that the tomato RLP members could be divided into eight subgroups and that the genes evolved independently compared to similar genes in Arabidopsis. Based on location and physical clustering analyses, we conclude that tomato RLPs likely expanded primarily through tandem duplication events. According to tissue specific RNA-seq data, 71 RLPs were expressed in at least one of the following tissues: root, leaf, bud, flower, or fruit. Several genes had expression patterns that were tissue specific. In addition, tomato RLP expression profiles after infection with different pathogens showed distinguish gene regulations according to disease induction and resistance response as well as infection by bacteria and virus. Notably, Some RLPs were highly and/or unique expressed in susceptible tomato to pathogen, suggesting that the RLP could be involved in disease response, possibly as a host-susceptibility factor. Our study could provide an important clues for further investigations into the function of tomato RLPs involved in developmental and response to pathogens.

Genome-wide identification of the AP2/ERF transcription factor family in pepper (Capsicum annuum L.)

The AP2/ERF family is one of the largest transcription factor families in the plant kingdom. AP2/ERF genes contributing to various processes including plant growth, development, and response to various stresses have been identified. In this study, 175 putative AP2/ERF genes were identified in the latest pepper genome database and classified into AP2, RAV, ERF, and Soloist subfamilies. Their chromosomal localization, gene structure, conserved motif, cis-acting elements within the promoter region, and subcellular locations were analyzed. Transient expression of CaAP2/ERF proteins in tobacco revealed that CaAP2/ERF064, CaAP2/ERF109, and CaAP2/ERF127 were located in the nucleus, while CaAP2/ERF171 was located in the nucleus and cytoplasm. Most of the CaAP2/ERF genes contained cis-elements within their promoter regions that responded to various stresses (HSE, LTR, MBS, Box-W1/W-box, and TC-rich repeats) and phytohormones (ABRE, CGTCA-motif, and TCA-element). Furthermore, RNA-seq analysis revealed that CaAP2/ERF genes showed differential expression profiles in various tissues as well as under biotic stresses. Moreover, qRT-PCR analysis of eight selected CaAP2/ERF genes also showed differential expression patterns in response to infection with Phytophthora capsici (HX-9) and in response to phytohormones (SA, MeJA, and ETH). This study will provide basic insights for further studies of the CaAP2/ERF genes involved in the interaction between pepper and P. capsici.

Genome-wide identification of the auxin/indole-3-acetic acid (Aux/IAA) gene family in pepper, its characterisation, and comprehensive expression profiling under environmental and phytohormones stress

Auxin is an essential phytohormone that plays a crucial role in the growth and development of plants in stressful environments. Here, we analysed the auxin/indole-3-acetic acid (Aux/IAA) gene family, which produces auxin in pepper, and succeeded in identifying 27 putative members containing four conserved domains (I. II. III and IV) in their protein sequences. Sequence analysis, chromosomal mapping and motif prediction of all identified CaAux/IAA genes were performed. It was observed that these genes contained four conserved motifs divided into nine different groups and distributed across nine chromosomes in pepper plants. RNA-seq analysis revealed the organ specific expression of many CaAux/IAA genes. However, the majority of genes were expressed with high expression levels in the early stages of fruit development. However, the maximum expression level of the CA03g34540 gene was observed in the breaker stage. Moreover, thirteen CaAux/IAA genes were labelled as early responsive genes to various phytohormone and abiotic stresses. Furthermore, RNA-seq analysis in response to pathogen inoculation (PepMoV, TMV strains P0/P1, and Phytophthora capsici) showed distinct expression profiles of all identified genes, suggesting the diverse expression nature of genes under these stress conditions. Overall, this study provides insight into the dynamic response of CaAux/IAA genes under environmental and phytohormones stress conditions, providing bases to further explore the importance of these genes through mutant/transgenic analysis in pepper.

Classification and Genome-Wide Analysis of Chitin-Binding Proteins Gene Family in Pepper (Capsicum annuum L.) and Transcriptional Regulation to Phytophthora capsici, Abiotic Stresses and Hormonal Applications

Chitin-binding proteins are pathogenesis-related gene family, which play a key role in the defense response of plants. However, thus far, little is known about the chitin-binding family genes in pepper (Capsicum annuum L.). In current study, 16 putative chitin genes (CaChi) were retrieved from the latest pepper genome database, and were classified into four distinct classes (I, III, IV and VI) based on their sequence structure and domain architectures. Furthermore, the structure of gene, genome location, gene duplication and phylogenetic relationship were examined to clarify a comprehensive background of the CaChi genes in pepper. The tissue-specific expression analysis of the CaChi showed the highest transcript levels in seed followed by stem, flower, leaf and root, whereas the lowest transcript levels were noted in red-fruit. Phytophthora capsici post inoculation, most of the CaChi (CaChiI3, CaChiIII1, CaChiIII2, CaChiIII4, CaChiIII6, CaChiIII7, CaChiIV1, CaChiVI1 and CaChiVI2) were induced by both strains (PC and HX-9). Under abiotic and exogenous hormonal treatments, the CaChiIII2, CaChiIII7, CaChiVI1 and CaChiVI2 were upregulated by abiotic stress, while CaChiI1, CaChiIII7, CaChiIV1 and CaChiIV2 responded to hormonal treatments. Furthermore, CaChiIV1-silenced plants display weakened defense by reducing (60%) root activity and increase susceptibility to NaCl stress. Gene ontology (GO) enrichment analysis revealed that CaChi genes primarily contribute in response to biotic, abiotic stresses and metabolic/catabolic process within the biological process category. These results exposed that CaChi genes are involved in defense response and signal transduction, suggesting their vital roles in growth regulation as well as response to stresses in pepper plant. In conclusion, these finding provide basic insights for functional validation of the CaChi genes in different biotic and abiotic stresses.

Global gene expression profiling for fruit organs and pathogen infections in the pepper, Capsicum annuum L

Hot pepper (Capsicum annuum) is one of the most consumed vegetable crops in the world and useful to human as it has many nutritional and medicinal values. Genomic resources of pepper are publically available since the pepper genomes have been completed and massive data such as transcriptomes have been deposited. Nevertheless, global transcriptome profiling is needed to identify molecular mechanisms related to agronomic traits in pepper, but limited analyses are published. Here, we report the comprehensive analysis of pepper transcriptomes during fruit ripening and pathogen infection. For the ripening, transcriptome data were obtained from placenta and pericarp at seven developmental stages. To reveal global transcriptomic landscapes during infection, leaves at six time points post-infection by one of three pathogens (Phytophthora infestans, Pepper mottle virus, and Tobacco mosaic virus P0 strain) were profiled. The massive parallel transcriptome profiling in this study will serve as a valuable resource for detection of molecular networks of fruit development and disease resistance in Capsicum annuum.

Genome-wide analysis of dirigent gene family in pepper (Capsicum annuum L.) and characterization of CaDIR7 in biotic and abiotic stresses

The dirigent (DIR and DIR-like) proteins involved in lignification, play a pivotal role against biotic and abiotic stresses in plants. However, no information is available about DIR gene family in pepper (Capsicum annuum L.). In this study, 24 putative dirigent genes (CaDIRs) were identified, their gene structure, genome location, gene duplication and phylogenetic relationship were elucidated. Tissue-specific expression analysis displayed the highest transcription levels in flower, stem and leaf. Some CaDIRs were up-regulated by virulent (CaDIR2, 3, 6, 7, 11, 14, 16, 22 and 23) and avirulent (CaDIR3, 5, 7, 16, 20, 22, 23 and 24) Phytophthora capsici strains, as well as by Methyl jasmonate, salicylic acid, NaCl and mannitol stresses. Acid-soluble lignin content increased (103.21%) after P. capsici inoculation (48-hour). Silencing of CaDIR7 weakened plant defense by reducing (~50%) root activity and made plants more susceptible (35.7%) to P. capsici and NaCl (300 mM). Leaf discs of the CaDIR7:silenced plants exposed to NaCl and mannitol (300 mM each), exhibited a significant decrease (56.25% and 48% respectively) in the chlorophyll content. These results suggested that CaDIR7 is involved in pepper defense response against pathogen and abiotic stresses and the study will provide basic insights for future research regarding CaDIRs.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

BACKGROUND

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

RESULTS

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

CONCLUSIONS

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Evolution analysis of Dof transcription factor family and their expression in response to multiple abiotic stresses in Malus domestica

As a family of transcription factors, DNA binding with one figure (Dof) proteins play important roles in various biological processes in plants. Here, a total of 60 putative apple (Malus domestica) Dof genes (MdDof) were identified and mapped to different chromosomes. Chromosomal distribution and synteny analysis indicated that the expansion of the MdDof genes came primarily from segmental and duplication events, and from whole genome duplication, which lead to more Dof members in apples than in other plants. All 60 MdDof genes were classified into thirteen groups, according to multiple sequence alignment and the phylogenetic tree constructed of Dof genes from apple, peach (Prunus persica), Arabidopsis and rice. Within each group, the members shared a similar exon/intron and motif compositions, although the sizes of the MdDof genes and encoding proteins were quite different. Several Dof genes from the apple and peach were identified to be homologues based on their close synteny relationship, which suggested that these genes bear similar functions. Half of the MdDof genes were randomly selected to determine their responses to different stresses. The majority of MdDof genes were quite sensitive to PEG, NaCl, cold and exogenous ABA treatment. Our results suggested that MdDof family members may play important roles in plant tolerance to abiotic stress.