Genome-Wide Identification and Expression Analysis of Two-Component System Genes in Tomato

Computational analysis and expression profiling of two-component system (TCS) gene family members in mango (Mangifera indica) indicated their roles in stress response

The two-component system (TCS) gene family is among the most important signal transduction families in plants and is involved in the regulation of various abiotic stresses, cell growth and division. To understand the role of TCS genes in mango (Mangifera indica ), a comprehensive analysis of TCS gene family was carried out in mango leading to identification of 65 MiTCS genes. Phylogenetic analysis divided MiTCSs into three groups (histidine kinases, histidine-containing phosphotransfer proteins, and response regulators) and 11 subgroups. One tandem duplication and 23 pairs of segmental duplicates were found within the MiTCSs . Promoter analysis revealed that MiTCSs contain a large number of cis -elements associated with environmental stresses, hormone response, light signalling, and plant development. Gene ontology analysis showed their involvement in various biological processes and molecular functions, particularly signal transduction. Protein-protein interaction analysis showed that MiTCS proteins interacted with each other. The expression pattern in various tissues and under many stresses (drought, cold, and disease) showed that expression levels varied among various genes in different conditions. MiTCSs 3D structure predictions showed structural conservation among members of the same groups. This information can be further used to develop improved cultivars and will serve as a foundation for gaining more functional insights into the TCS gene family.

Surviving a Double-Edged Sword: Response of Horticultural Crops to Multiple Abiotic Stressors

Climate change-induced weather events, such as extreme temperatures, prolonged drought spells, or flooding, pose an enormous risk to crop productivity. Studies on the implications of multiple stresses may vary from those on a single stress. Usually, these stresses coincide, amplifying the extent of collateral damage and contributing to significant financial losses. The breadth of investigations focusing on the response of horticultural crops to a single abiotic stress is immense. However, the tolerance mechanisms of horticultural crops to multiple abiotic stresses remain poorly understood. In this review, we described the most prevalent types of abiotic stresses that occur simultaneously and discussed them in in-depth detail regarding the physiological and molecular responses of horticultural crops. In particular, we discussed the transcriptional, posttranscriptional, and metabolic responses of horticultural crops to multiple abiotic stresses. Strategies to breed multi-stress-resilient lines have been presented. Our manuscript presents an interesting amount of proposed knowledge that could be valuable in generating resilient genotypes for multiple stressors.

Genome-wide identification of the melon (Cucumis melo L.) response regulator gene family and functional analysis of CmRR6 and CmPRR3 in response to cold stress

The response regulator (RR) gene family play crucial roles in cytokinin signal transduction, plant development, and resistance to abiotic stress. However, there are no reports on the identification and functional characterization of RR genes in melon. In this study, a total of 18 CmRRs were identified and classified into type A, type B, and clock PRRs, based on phylogenetic analysis. Most of the CmRRs displayed tissue-specific expression patterns, and some were induced by cold stress according to two RNA-seq datasets. The expression patterns of CmRR2/6/11/15 and CmPRR2/3 under cold treatment were confirmed by qRT-PCR. Subcellular localization assays indicated that CmRR6 and CmPRR3 were primarily localized in the nucleus and chloroplast. Furthermore, when either CmRR6 or CmPRR3 were silenced using tobacco ringspot virus (TRSV), the cold tolerance of the virus-induced gene silencing (VIGS) melon plants were significantly enhanced, as evidenced by measurements of chlorophyll fluorescence, ion leakage, reactive oxygen, proline, and malondialdehyde levels. Additionally, the expression levels of CmCBF1, CmCBF2, and CmCBF3 were significantly increased in CmRR6-silenced and CmPRR3-silenced plants under cold treatment. Our findings suggest that CmRRs contribute to cold stress responses and provide new insights for further pursuing the molecular mechanisms underlying CmRRs-mediated cold tolerance in melon.

Phenomic and transcriptomic analyses reveal the sequential synthesis of Fe3O4 nanoparticles in Acidithiobacillus ferrooxidans BYM

IMPORTANCE

As the most important non-magnetotactic magnetosome-producing bacteria, Acidithiobacillus ferrooxidans only requires very mild conditions to produce Fe3O4 nanoparticles, thus conferring greater flexibility and potential application in biomagnetic nanoparticle production. However, the available information cannot explain the mechanism of Fe3O4 nanoparticle formation in A. ferrooxidans. In this study, we applied phenomic and transcriptomic analyses to reveal this mechanism. We found that different treatment condition factors notably affect the phenomic data of Fe3O4 nanoparticle in A. ferrooxidans. Using transcriptomic analyses, the gene network controlling/regulating Fe3O4 nanoparticle biogenesis in A. ferrooxidans was proposed, excavating the candidate hub genes for Fe3O4 nanoparticle formation in A. ferrooxidans. Based on this information, a sequential model for Fe3O4 nanoparticle synthesis in A. ferrooxidans was hypothesized. It lays the groundwork for further clarifying the feature of Fe3O4 nanoparticle synthesis.

Two-Component System Genes in Brassica napus: Identification, Analysis, and Expression Patterns in Response to Abiotic and Biotic Stresses

The two-component system (TCS), consisting of histidine kinases (HKs), histidine phosphotransfer proteins (HPs) and response regulators (RRs) in eukaryotes, plays pivotal roles in regulating plant growth, development, and responses to environment stimuli. However, the TCS genes were poorly characterized in rapeseed, which is an important tetraploid crop in Brassicaceae. In this work, a total of 182 BnaTCS genes were identified, including 43 HKs, 16 HPs, and 123 RRs, which was more than that in other crops due to segmental duplications during the process of polyploidization. It was significantly different in genetic diversity between the three subfamilies, and some members showed substantial genetic differentiation among the three rapeseed ecotypes. Several hormone- and stress-responsive cis-elements were identified in the putative promoter regions of BnaTCS genes. Furthermore, the expression of BnaTCS genes under abiotic stresses, exogenous phytohormone, and biotic stresses was analyzed, and numerous candidate stress-responsive genes were screened out. Meanwhile, using a natural population with 505 B. napus accessions, we explored the genetic effects of BnaTCS genes on salt tolerance by association mapping analysis and detected some significant association SNPs/genes. The result will help to further understand the functions of TCS genes in the developmental and stress tolerance improvement in B. napus.

Tomato salt-responsive pseudo-response regulator 1, SlSRP1, negatively regulates the high-salt and dehydration stress responses

Under severe environmental stress conditions, plants inhibit their growth and development and initiate various defense mechanisms to survive. The pseudo-response regulator (PRRs) genes have been known to be involved in fruit ripening and plant immunity in various plant species, but their role in responses to environmental stresses, especially high salinity and dehydration, remains unclear. Here, we focused on PRRs in tomato plants and identified two PRR2-like genes, SlSRP1 and SlSRP1H, from the leaves of salt-treated tomato plants. After exposure to dehydration and high-salt stresses, expression of SISRP1, but not SlSRP1H, was significantly induced in tomato leaves. Subcellular localization analysis showed that SlSRP1 was predominantly located in the nucleus, while SlSRP1H was equally distributed in the nucleus and cytoplasm. To further investigate the potential role of SlSRP1 in the osmotic stress response, we generated SISRP1-silenced tomato plants. Compared to control plants, SISRP1-silenced tomato plants exhibited enhanced tolerance to high salinity, as evidenced by a high accumulation of proline and reduced chlorosis, ion leakage, and lipid peroxidation. Moreover, SISRP1-silenced tomato plants showed dehydration-tolerant phenotypes with enhanced abscisic acid sensitivity and increased expression of stress-related genes, including SlRD29, SlAREB, and SlDREB2. Overall, our findings suggest that SlSRP1 negatively regulates the osmotic stress response.

Characterization of Two-Component System gene (TCS) in melatonin-treated common bean under salt and drought stress

The two-component system (TCS) generally consists of three elements, namely the histidine kinase (HK), response regulator (RR), and histidine phosphotransfer (HP) gene families. This study aimed to assess the expression of TCS genes in P. vulgaris leaf tissue under salt and drought stress and perform a genome-wide analysis of TCS gene family members using bioinformatics methods. This study identified 67 PvTCS genes, including 10 PvHP, 38 PvRR, and 19 PvHK, in the bean genome. PvHK2 had the maximum number of amino acids with 1261, whilst PvHP8 had the lowest number with 87. In addition, their theoretical isoelectric points were between 4.56 (PvHP8) and 9.15 (PvPRR10). The majority of PvTCS genes are unstable. Phylogenetic analysis of TCS genes in A. thaliana, G. max, and bean found that PvTCS genes had close phylogenetic relationships with the genes of other plants. Segmental and tandem duplicate gene pairs were detected among the TCS genes and TCS genes have been subjected to purifying selection pressure in the evolutionary process. Furthermore, the TCS gene family, which has an important role in abiotic stress and hormonal responses in plants, was characterized for the first time in beans, and its expression of TCS genes in bean leaves under salt and drought stress was established using RNAseq and qRT-PCR analyses. The findings of this study will aid future functional and genomic studies by providing essential information about the members of the TCS gene family in beans.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01406-5.

Cytokinin-inducible response regulator SlRR6 controls plant height through gibberellin and auxin pathways in tomato

Plant height is a key agronomic trait regulated by several phytohormones such as gibberellins (GAs) and auxin. However, little is known about how cytokinin (CK) participates in this process. Here, we report that SlRR6, a type-A response regulator in the CK signaling pathway, positively regulates plant height in tomato. SlRR6 was induced by exogenous kinetin and GA3, but inhibited by indole-3-acetic acid (IAA). Knock out of SlRR6 reduced tomato plant height through shortening internode length, while overexpression of SlRR6 caused taller plants due to increased internode number. Cytological observation of longitudinal stems showed that both knock out and overexpression of SlRR6 generated larger cells, but significantly reduced cell numbers in each internode. Further studies demonstrated that overexpression of SlRR6 enhanced GA accumulation and lowered IAA content, along with expression changes in GA- and IAA-related genes. Exogenous paclobutrazol and IAA treatments restored the increased plant height phenotype in SlRR6-overexpressing lines. Yeast two-hybrid, bimolecular fluorescence complementation, and co-immunoprecipitation assays showed that SlRR6 interacts with a small auxin up RNA protein, SlSAUR58. Moreover, SlSAUR58-overexpressing plants were dwarf with decreased internode length. Overall, our findings establish SlRR6 as a vital component in the CK signaling, GA, and IAA regulatory network that controls plant height.

Genome-wide identification and expression profiling of two-component system (TCS) genes in Brassica oleracea in response to shade stress

The Two-component system (TCS) consists of Histidine kinases (HKs), Phosphotransfers (HPs), and response regulator (RR) proteins. It has an important role in signal transduction to respond to a wide variety of abiotic stresses and hence in plant development. Brassica oleracea (cabbage) is a leafy vegetable, which is used for food and medicinal purposes. Although this system was identified in several plants, it had not been identified in Brassica oleracea yet. This genome-wide study identified 80 BoTCS genes consisting of 21 HKs, 8 HPs, 39 RRs, and 12 PRRs. This classification was done based on conserved domains and motif structure. Phylogenetic relationships of BoTCS genes with Arabidopsis thaliana, Oryza sativa, Glycine max, and Cicer arietinum showed conservation in TCS genes. Gene structure analysis revealed that each subfamily had conserved introns and exons. Both tandem and segmental duplication led to the expansion of this gene family. Almost all of the HPs and RRs were expanded through segmental duplication. Chromosomal analysis showed that BoTCS genes were dispersed across all nine chromosomes. The promoter regions of these genes were found to contain a variety of cis-regulatory elements. The 3D structure prediction of proteins also confirmed the conservation of structure within subfamilies. MicroRNAs (miRNAs) involved in the regulation of BoTCSs were also predicted and their regulatory roles were also evaluated. Moreover, BoTCSs were docked with abscisic acid to evaluate their binding. RNA-seq-based expression analysis and validation by qRT-PCR showed significant variation of expression for BoPHYs, BoERS1.1, BoERS2.1, BoERS2.2, BoRR10.2, and BoRR7.1 suggesting their importance in stress response. These genes showing unique expression can be further used in manipulating the plant's genome to make the plant more resistant the environmental stresses which will ultimately help in the increase of plant's yield. More specifically, these genes have altered expression in shade stress which clearly indicates their importance in biological functions. These findings are important for future functional characterization of TCS genes in generating stress-responsive cultivars.

Omics-assisted characterization of two-component system genes from Gossypium Raimondii in response to salinity and molecular interaction with abscisic acid

The two-component system (TCS) genes are involved in a wide range of physiological processes in prokaryotes and eukaryotes. In plants, the TCS elements help in a variety of functions, including cell proliferation, response to abiotic and biotic stresses, leaf senescence, nutritional signaling, and division of chloroplasts. Three different kinds of proteins make up the TCS system in plants. These are known as HKs (histidine kinases), HPs (histidine phosphotransfer), and RRs (response regulators). We investigated the genome of Gossypium raimondii and discovered a total of 59 GrTCS candidates, which include 23 members of the HK family, 8 members of the HP family, and 28 members of the RR family. RR candidates are further classified as type-A (6 members), type-B (11 members), type-C (2 members), and pseudo-RRs (9 members). The GrTCS genes were analyzed in comparison with the TCS components of other plant species such as Arabidopsis thaliana, Cicer arietinum, Sorghum bicolor, Glycine max, and Oryza sativa. This analysis revealed both conservation and changes in their structures. We identified 5 pairs of GrTCS syntenic homologs in the G. raimondii genome. All 59 TCS genes in G. raimondii are located on all thirteen chromosomes. The GrTCS promoter regions have several cis-regulatory elements, which function as switches and respond to a wide variety of abiotic stresses. RNA-seq and real-time qPCR analysis showed that the majority of GrTCS genes are differentially regulated in response to salt and cold stress. 3D structures of GrTCS proteins were predicted to reveal the specific function. GrTCSs were docked with abscisic acid to assess their binding interactions. This research establishes the groundwork for future functional studies of TCS elements in G. raimondii, which will further focus on stress resistance and overall development.

The mutation of CaCKI1 causes seedless fruits in chili pepper (Capsicum annuum)

The seedless mutant tn-1 in chili pepper is caused by a mutation in CaCKI1 (CA12g21620), which encodes histidine kinase involving female gametophyte development. An amino acid insertion in the receiver domain of CaCKI1 may be the mutation responsible for tn-1. Seedlessness is a desirable trait in fruit crops because the removal of seeds is a troublesome step for consumers and processing industries. However, little knowledge is available to develop seedless chili peppers. In a previous study, a chili pepper mutant tn-1, which stably produces seedless fruits, was isolated. In this study, we report characterization of tn-1 and identification of the causative gene. Although pollen germination was normal, confocal laser microscopy observations revealed deficiency in embryo sac development in tn-1. By marker analysis, the tn-1 locus was narrowed down to a 313 kb region on chromosome 12. Further analysis combined with mapping-by-sequencing identified CA12g21620, which encodes histidine kinase as a candidate gene. Phylogenetic analysis revealed CA12g21620 was the homolog of Arabidopsis CKI1 (Cytokinin Independent 1), which plays an important role in female gametophyte development, and CA12g21620 was designated as CaCKI1. Sequence analysis revealed that tn-1 has a 3-bp insertion in the 6th exon resulting in one lysine (K) residue insertion in receiver domain of CaCKI1, and the sequence nearby the insertion is widely conserved among CKI1 orthologs in various plants. This suggested that one K residue insertion may reduce the phosphorylation relay downstream of CaCKI1 and impair normal development of female gametophyte, resulting in seedless fruits production in tn-1. Furthermore, we demonstrated that virus-induced gene silencing of CaCKI1 reduced normally developed female gametophyte in chili pepper. This study describes the significant role of CaCKI1 in seed development in chili pepper and the possibility of developing seedless cultivars using its mutation.

Genome-wide identification and expression analysis of the response regulator gene family in alfalfa (Medicago sativa L.) reveals their multifarious roles in stress response

As important components of the two-component regulatory system, response regulatory proteins (RRPs) play a crucial role in histidine phosphorylation-mediated signal transduction in response to environmental fluctuations. Accumulating evidence has revealed that RRPs play important roles in plant growth and stress response. However, the specific functions of RR genes (RRs) in cultivated alfalfa remain ambiguous. Therefore, in this study, we identified and characterized the RR family genes in the alfalfa genome using bioinformatics methods. Our analysis revealed 37 RRs in the alfalfa genome of Zhongmu No.1 that were unevenly distributed on the chromosomes. Cis-elements analysis revealed the involvement of RRs in responses to light, stress, and various plant hormones. Expression analysis of RRs in different tissues revealed their distinct tissue expression patterns. These findings provide preliminary insights into the roles of RRs in plant responses to abiotic stress, which can be used to improve the stress tolerance of autotetraploid-cultivated alfalfa plants via genetic engineering.

Prunus Knotted-like Genes: Genome-Wide Analysis, Transcriptional Response to Cytokinin in Micropropagation, and Rootstock Transformation

Knotted1-like homeobox (KNOX) transcription factors are involved in plant development, playing complex roles in aerial organs. As Prunus species include important fruit tree crops of Italy, an exhaustive investigation of KNOX genes was performed using genomic and RNA-seq meta-analyses. Micropropagation is an essential technology for rootstock multiplication; hence, we investigated KNOX transcriptional behavior upon increasing 6-benzylaminopurine (BA) doses and the effects on GF677 propagules. Moreover, gene function in Prunus spp. was assessed by Gisela 6 rootstock transformation using fluorescence and peach KNOX transgenes. Based on ten Prunus spp., KNOX proteins fit into I-II-M classes named after Arabidopsis. Gene number, class member distribution, and chromosome positions were maintained, and exceptions supported the diversification of Prunus from Cerasus subgenera, and that of Armeniaca from the other sections within Prunus. Cytokinin (CK) cis-elements occurred in peach and almond KNOX promoters, suggesting a BA regulatory role in GF677 shoot multiplication as confirmed by KNOX expression variation dependent on dose, time, and interaction. The tripled BA concentration exacerbated stress, altered CK perception genes, and modified KNOX transcriptions, which are proposed to concur in in vitro anomalies. Finally, Gisela 6 transformation efficiency varied (2.6-0.6%) with the genetic construct, with 35S:GFP being more stable than 35S:KNOPE1 lines, which showed leaf modification typical of KNOX overexpression.

Jasmonic acid enhances osmotic stress responses by MYC2-mediated inhibition of protein phosphatase 2C1 and response regulators 26 transcription factor in tomato

The jasmonic acid (JA) signaling pathway is involved in the plant response to drought stress. JA and other hormones synergistically regulate the drought response in plants. However, the molecular mechanism underlying this synergism remains poorly defined. In the present study, transcriptome analyses of guard cells and quantitative PCR experiments revealed that MYC2 negatively regulated the negative regulator of ABA signaling, SlPP2C1, and the type-B response regulator in the cytokinin pathway, SlRR26, and this negative regulation was direct. SlRR26 overexpression reduced drought tolerance in transgenic tomatoes, whereas slrr26^cr lines were more tolerant to drought. SlRR26 negatively modulated reactive oxygen species levels in stomata and stomatal closure through RobhB. Moreover, SlRR26 overexpression counteracted JA-mediated stomatal closure, suggesting that SlRR26 played a negative role in the JA-mediated drought response. These findings suggest that MYC2 plays a key role in JA-regulated stomatal closure under drought stress by inhibiting SlPP2C1 and SlRR26.

Genome-wide identification and expression analysis of two-component system genes in sweet potato (Ipomoea batatas L.)

Two-component system (TCS), which comprises histidine kinases (HKs), histidine phosphotransfer proteins (HPs), and response regulators (RRs), plays essential roles in regulating plant growth, development, and response to various environmental stimuli. TCS genes have been comprehensively identified in various plants, while studies on the genome-wide identification and analysis of TCS in sweet potato were still not reported. Therefore, in this study, a total of 90 TCS members consisting of 20 HK(L)s, 11 HPs, and 59 RRs were identified in the genome of Ipomoea batatas. Furthermore, their gene structures, conserved domains, and phylogenetic relationships were analyzed in detail. Additionally, the gene expression profiles in various organs were analyzed, and response patterns to adverse environmental stresses were investigated. The results showed that these 90 TCS genes were mapped on 15 chromosomes with a notably uneven distribution, and the expansion of TCS genes in sweet potato was attributed to both segmental and tandem duplications. The majority of the TCS genes showed distinct organ-specific expression profiles, especially in three types of roots (stem roots, fibrous roots, tuberous roots). Moreover, most of the TCS genes were either induced or suppressed upon treatment with abiotic stresses (drought, salinity, cold, heat) and exogenous phytohormone abscisic acid (ABA). In addition, the yeast-two hybrid system was used to reveal the HK-HP-RR protein-protein interactions. IbHP1, IbHP2, IbHP4, and IbHP5 could interact with three HKs (IbHK1a, IbHK1b, and IbHK5), and also interact with majority of the type-B RRs (IbRR20-IbRR28), while no interaction affinity was detected for IbHP3. Our systematic analyses could provide insights into the characterization of the TCS genes, and further the development of functional studies in sweet potato.

Pangenome-wide analysis of cyclic nucleotide-gated channel (CNGC) gene family in citrus Spp. Revealed their intraspecies diversity and potential roles in abiotic stress tolerance

Cyclic nucleotide-gated channels (CNGC) gene family has been found to be involved in physiological processes including signaling pathways, environmental stresses, plant growth, and development. This gene family of non-selective cation channels is known to regulate the uptake of calcium and is reported in several plant species. The pangenome-wide studies enable researchers to understand the genetic diversity comprehensively; as a comparative analysis of multiple plant species or member of a species at once helps to better understand the evolutionary relationships and diversity present among them. In the current study, pangenome-wide analysis of the CNGC gene family has been performed on five Citrus species. As a result, a total of 32 genes in Citrus sinensis, 27 genes in Citrus recticulata, 30 genes in Citrus grandis, 31 genes in Atalantia buxfolia, and 30 genes in Poncirus trifoliata were identified. In addition, two unique genes CNGC13 and CNGC14 were identified, which may have potential roles. All the identified CNGC genes were unevenly distributed on 9 chromosomes except P. trifoliata had genes distributed on 7 chromosomes and were classified into four major groups and two sub-groups namely I, II, III, IV-A, and IV-B. Cyclic nucleotide binding (CNB) motif, calmodulin-binding motif (CaMB), and motif for IQ-domain were conserved in Citrus Spp. Intron exon structures of citrus species were not exactly as same as the gene structures of Arabidopsis. The majority of cis-regulatory elements (CREs) were light responsive and others include growth, development, and stress-related indicating potential roles of the CNGC gene family in these functions. Both segmental and tandem duplication were involved in the expansion of the CNGC gene family in Citrus Spp. The miRNAs are involved in the response of CsCNGC genes towards drought stress along with having regulatory association in the expression of these genes. Protein- Protein interaction (PPI) analysis also showed the interaction of CNGC proteins with other CNGCs which suggested their potential role in pathways regulating different biological processes. GO enrichment revealed that CNGC genes were involved in the transport of ions across membranes. Furthermore, tissue-specific expression patterns of leaves sample of C. sinensis were studied under drought stress. Out of 32 genes of C. sinensis 3 genes i.e., CsCNGC1.4, CsCNGC2.1, and CsCNGC4.2 were highly up-regulated, and only CsCNGC4.6 was highly down-regulated. The qRT-PCR analysis also showed that CNGC genes were highly expressed after treatment with drought stress, while gene expression was lower under controlled conditions. This work includes findings based on multiple genomes instead of one, therefore, this will provide more genomic information rather than single genome-based studies. These findings will serve as a basis for further functional insights into the CNGC gene family.

Genome wide identification and characterization of nodulation related genes in Arachis hypogaea

Nitrogen is an important plant nutrient that has a significant role in crop yield. Hence, to fulfill the needs of sustainable agriculture, it is necessary to improve biological nitrogen fixation in leguminous crops. Nod inducing gene families plays a crucial role in the interaction between rhizobia and legumes, leading to biological nitrogen fixation. However, nod inducing genes identification and characterization has not yet been performed in Arachis hypogaea. In this study, identification and genome-wide analysis of nod inducing genes are performed so that to explore their potential functions in the Arachis hypogaea for the first time. Nod genes were comprehensively analyzed by phylogenetic clustering analysis, gene structure determination, detection of conserved motifs, subcellular localization, conserved motifs, cis-acting elements and promoter region analysis. This study identified 42 Nod inducing genes in Arachis hypogaea, their sequences were submitted to NCBI and accession numbers were obtained. Potential involvement of these genes in biological nitrogen fixation has been unraveled, such as, phylogenetic analysis revealed that nod inducing genes evolved independently in Arachis hypogaea, the amino acid structures exhibited 20 highly conserved motifs, the proteins are present at different locations in cells and the gene structures revealed that all the genes are full-length genes with upstream intronic regions. Further, the promoter analysis determined a large number of cis-regulatory elements involved in nodulation. Moreover, this study not only provides identification and characterization of genes underlying developmental and functional stages of nodulation and biological nitrogen fixation but also lays the foundation for further revelation of nod inducing gene family. Besides, identification and structural analysis of these genes in Arachis hypogaea may provide a theoretical basis for the study of evolutionary relationships in future analysis.

Two-Component System Genes in Sorghum bicolor: Genome-Wide Identification and Expression Profiling in Response to Environmental Stresses

The two-component signal transduction system (TCS) acts in a variety of physiological processes in lower organisms and has emerged as a key signaling system in both prokaryotes and eukaryotes, including plants. TCS genes assist plants in processes such as stress resistance, cell division, nutrition signaling, leaf senescence, and chloroplast division. In plants, this system is composed of three types of proteins: response regulators (RRs), histidine kinases (HKs), and histidine phosphotransfer proteins (HPs). We aimed to study the Sorghum bicolor genome and identified 37 SbTCS genes consisting of 13 HKs, 5 HPs, and 19 RRs (3 type-A RRs, 7 type-B RRs, 2 type-C RRs, and 7 pseudo-RRs). The structural and phylogenetic comparison of the SbTCS members with their counterparts in Arabidopsis thaliana, Oryza sativa, Cicer arietinum, and Glycine max showed group-specific conservations and variations. Expansion of the gene family members is mostly a result of gene duplication, of both the tandem and segmental types. HKs and RRs were observed to be originated from segmental duplication, while some HPs originated from tandem duplication. The nuclear genome of S. bicolor contain 10 chromosomes and these SbTCS genes are randomly distributed on all the chromosomes. The promoter sequences of the SbTCS genes contain several abiotic stress-related cis-elements. RNA-seq and qRT-PCR-based expression analysis demonstrated most of the TCS genes were responsive to drought and salt stresses in leaves, which suggest their role in leaf development. This study lays a foundation for further functional study of TCS genes for stress tolerance and developmental improvement in S. bicolor.

Genome-wide analysis of ARF transcription factors reveals HcARF5 expression profile associated with the biosynthesis of β-ocimene synthase in Hedychium coronarium

Herein, 37 ARF genes were identified and analyzed in Hedychium coronarium and HcARF5 showed a potential role in the regulation of HcTPS3. Auxin is an important plant hormone, implicated in various aspects of plant growth and development processes especially in the biosynthesis of various secondary metabolites. Auxin response factors (ARF) belong to the transcription factors (TFs) gene family and play a crucial role in transcriptional activation/repression of auxin-responsive genes by directly binding to their promoter region. Nevertheless, whether ARF genes are involved in the regulatory mechanism of volatile compounds in flowering plants is largely unknown. β-ocimene is a key floral volatile compound synthesized by terpene synthase 3 (HcTPS3) in Hedychium coronarium. A comprehensive analysis of H. coronarium genome reveals 37 candidate ARF genes in the whole genome. Tissue-specific expression patterns of HcARFs family members were assessed using available transcriptome data. Among them, HcARF5 showed a higher expression level in flowers, and significantly correlated with the key structural β-ocimene synthesis gene (HcTPS3). Furthermore, transcript levels of both genes were associated with the flower development. Under hormone treatments, the response of HcARF5 and HcTPS3, and the emission level of β-ocimene contents were evaluated. Subcellular and transcriptional activity assay showed that HcARF5 localizes to the nucleus and possesses transcriptional activity. Yeast one-hybrid (Y1H) and dual-luciferase assays revealed that HcARF5 directly regulates the transcriptional activity of HcTPS3. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that HcARF5 interacts with scent-related HcIAA4, HcIAA6, and HcMYB1 in vivo. Overall, these results indicate that HcARF5 is potentially involved in the regulation of β-ocimene synthesis in H. coronarium.

Genome-wide analysis, identification, evolution and genomic organization of dehydration responsive element-binding (DREB) gene family in Solanum tuberosum

Background

The dehydration responsive element-binding (DREB) gene family plays a crucial role as transcription regulators and enhances plant tolerance to abiotic stresses. Although the DREB gene family has been identified and characterized in many plants, knowledge about it in Solanum tuberosum (Potato) is limited.

Results

In the present study, StDREB gene family was comprehensively analyzed using bioinformatics approaches. We identified 66 StDREB genes through genome wide screening of the Potato genome based on the AP2 domain architecture and amino acid conservation analysis (Valine at position 14th). Phylogenetic analysis divided them into six distinct subgroups (A1–A6). The categorization of StDREB genes into six subgroups was further supported by gene structure and conserved motif analysis. Potato DREB genes were found to be distributed unevenly across 12 chromosomes. Gene duplication proved that StDREB genes experienced tandem and segmental duplication events which led to the expansion of the gene family. The Ka/Ks ratios of the orthologous pairs also demonstrated the StDREB genes were under strong purification selection in the course of evolution. Interspecies synteny analysis revealed 45 and 36 StDREB genes were orthologous to Arabidopsis and Solanum lycopersicum, respectively. Moreover, subcellular localization indicated that StDREB genes were predominantly located within the nucleus and the StDREB family’s major function was DNA binding according to gene ontology (GO) annotation.

Conclusions

This study provides a comprehensive and systematic understanding of precise molecular mechanism and functional characterization of StDREB genes in abiotic stress responses and will lead to improvement in Solanum tuberosum.

Histidine Kinases: Diverse Functions in Plant Development and Responses to Environmental Conditions

The two-component system (TCS), which is one of the most evolutionarily conserved signaling pathway systems, has been known to regulate multiple biological activities and environmental responses in plants. Significant progress has been made in characterizing the biological functions of the TCS components, including signal receptor histidine kinase (HK) proteins, signal transducer histidine-containing phosphotransfer proteins, and effector response regulator proteins. In this review, our scope is focused on the diverse structure, subcellular localization, and interactions of the HK proteins, as well as their signaling functions during development and environmental responses across different plant species. Based on data collected from scientific studies, knowledge about acting mechanisms and regulatory roles of HK proteins is presented. This comprehensive summary ofthe HK-related network provides a panorama of sophisticated modulating activities of HK members and gaps in understanding these activities, as well as the basis for developing biotechnological strategies to enhance the quality of crop plants.

Ectopic expression of GmHP08 enhances resistance of transgenic Arabidopsis toward drought stress

Ectopic expression of Glycine max two-component system member GmHP08 in Arabidopsis enhanced drought tolerance of transgenic plants, possibly via ABA-dependent pathways. Phosphorelay by two-component system (TCS) is a signal transduction mechanism which has been evolutionarily conserved in both prokaryotic and eukaryotic organisms. Previous studies have provided lines of evidence on the involvement of TCS genes in plant perception and responses to environmental stimuli. In this research, drought-associated functions of GmHP08, a TCS member from soybean (Glycine max L.), were investigated via its ectopic expression in Arabidopsis system. Results from the drought survival assay showed that GmHP08-transgenic plants exhibited higher survival rates compared with their wild-type (WT) counterparts, indicating better drought resistance of the former group. Analyses revealed that the transgenic plants outperformed the WT in various regards, i.e. capability of water retention, prevention of hydrogen peroxide accumulation and enhancement of antioxidant enzymatic activities under water-deficit conditions. Additionally, the expression of stress-marker genes, especially antioxidant enzyme-encoding genes, in the transgenic plants were found greater than that of the WT plants. In contrary, the expression of SAG13 gene, one of the senescence-associated genes, and of several abscisic acid (ABA)-related genes was repressed. Data from this study also revealed that the ectopic expression lines at germination and early seedling development stages were hypersensitive to exogenous ABA treatment. Taken together, our results demonstrated that GmHP08 could play an important role in mediating plant response to drought, possibly via an ABA-dependent manner.

The prohibitins (PHB) gene family in tomato: Bioinformatic identification and expression analysis under abiotic and phytohormone stresses

The prohibitins (PHB) are SPFH domain-containing proteins found in the prokaryotes to eukaryotes. The plant PHBs are associated with a wide range of biological processes, including senescence, development, and responses to biotic and abiotic stresses. The PHB proteins are identified and characterized in the number of plant species, such as Arabidopsis, rice, maize, and soybean. However, no systematic identification of PHB proteins was performed in Solanum lycopersicum. In this study, we identified 16 PHB proteins in the tomato genome. The analysis of conserved motifs and gene structure validated the phylogenetic classification of tomato PHB proteins. It was observed that various members of tomato PHB proteins undergo purifying selection based on the Ka/Ks ratio and are targeted by four families of miRNAs. Moreover, SlPHB proteins displayed a very unique expression pattern in different plant parts including fruits at various development stages. It was found that SlPHBs processed various development-related and phytohormone responsive cis-regulatory elements in their promoter regions. Furthermore, the exogenous phytohormones treatments (Abscisic acid, indole-3-acetic acid, gibberellic acid, methyl jasmonate) salt and drought stresses induce the expression of SlPHB. Moreover, the subcellular localization assay revealed that SlPHB5 and SlPHB10 were located in the mitochondria. This study systematically summarized the general characterization of SlPHBs in the tomato genome and provides a foundation for the functional characterization of PHB genes in tomato and other plant species.

The DUF221 domain-containing (DDP) genes identification and expression analysis in tomato under abiotic and phytohormone stress

The domain of unknown function (DUF221 domain-containing) proteins regulates various aspects of plant growth, development, responses to abiotic stresses, and hormone transduction pathways. To understand the role of DDP proteins in tomato, a comprehensive genome-wide analysis was performed in the tomato genome. A total of 12 DDP genes were identified and distributed in 8 chromosomes in the tomato genome. Phylogenetically all SlDDPs were clustered into four clades, subsequently supported by their gene structure and conserved motifs distribution. The SlDDPs contained various cis-acting elements involved in plant responses to abiotic and various phytohormones stresses. The tissue-specific expression profile analysis revealed the constitutive expression of SlDDPs in roots, leaves, and developmental phases of fruit. It was found that SlDDP1, SlDDP3, SlDDP4, SlDDP9, SlDDP10, and SlDDP12 exhibited high expression levels in fruits at different development stages. Of these genes, SlDDP12 contained ethylene (ERE) responsive elements in their promoter regions, suggesting its role in ethylene-dependent fruit ripening. It was found that a single SlDDP induced by two or more abiotic and phytohormone stresses. These include, SlDDP1, SlDDP2, SlDDP3, SlDDP4, SlDDP7, SlDDP8, and SlDDP10 was induced under salt, drought, ABA, and IAA stresses. Moreover, tomato SlDDPs were targeted by multiple miRNA gene families as well. In conclusion, this study predicted that the putative DDP genes might help improve abiotic and phytohormone tolerance in plants, particularly tomato, rice, and other economically important crop plant species.

Genome-Wide Characterization and Expression of Two-Component System Genes in Cytokinin-Regulated Gall Formation in Zizania latifolia

The thickening of Zizania latifolia shoots, referred to as gall formation, depends on infection with the fungal endophyte Ustilago esculenta. The swollen and juicy shoots are a popular vegetable in Asia. A key role for cytokinin action in this process was postulated. Here, trans-zeatin stimulated swelling in fungi-infected Z. latifolia. A two-component system (TCS) linked cytokinin binding to receptors with transcriptional regulation in the nucleus and played important roles in diverse biological processes. We characterized 69 TCS genes encoding for 25 histidine kinase/histidine-kinase-like (HK(L)) (21 HKs and 4 HKLs), 8 histidine phosphotransfer proteins (HP) (5 authentic and 3 pseudo), and 36 response regulators (RR; 14 type A, 14 type B, 2 type C, and 6 pseudo) in the genome of Z. latifolia. These TCS genes have a close phylogenetic relationship with their rice counterparts. Nineteen duplicated TCS gene pairs were found and the ratio of nonsynonymous to synonymous mutations indicated that a strong purifying selection acted on these duplicated genes, leading to few mutations during evolution. Finally, ZlCHK1, ZlRRA5, ZIRRA9, ZlRRA10, ZlPRR1, and ZlPHYA expression was associated with gall formation. Among them, ARR5, ARR9, and ZlPHYA are quickly induced by trans-zeatin, suggesting a role for cytokinin signaling in shoot swelling of Z. latifolia.

Genome-Wide Characterization, Evolution, and Expression Analysis of the Leucine-Rich Repeat Receptor-Like Protein Kinase (LRR-RLK) Gene Family in Medicago truncatula

Leucine-rich repeat receptor-like kinases (LRR-RLKs) constitute the largest subfamily of receptor-like kinases (RLKs) in plants. They play roles in plant growth and developmental and physiological processes, but less is known about the functions of LRR-RLKs in Medicago truncatula. Our genome-wide analysis revealed 329 LRR-RLK genes in the M.truncatula genome. Phylogenetic and classification analysis suggested that these genes could be classified into 15 groups and 24 subgroups. A total of 321 genes were mapped onto all chromosomes, and 23 tandem duplications (TDs) involving 56 genes were distributed on each chromosome except 4. Twenty-seven M.truncatula LRR-RLK segmental duplication gene pairs were colinearly related. The exon/intron organization, motif composition and arrangements were relatively conserved among members of the same groups or subgroups. Using publicly available RNAseq data and quantitative real-time polymerase chain reaction (qRT-PCR), expression profiling suggested that LRR-RLKs were differentially expressed among different tissues, while some were expressed specifically in the roots and nodules. The expression of LRR-RLKs in A17 and 4 nodule mutants under rhizobial infection showed that 36 LRR-RKLs were highly upregulated in the sickle (skl) mutant [an ethylene (ET)-insensitive, Nod factor-hypersensitive mutant] after 12 h of rhizobium inoculation. Among these LRR-RLKs, six genes were also expressed specifically in the roots and nodules, which might be specific to the Nod factor and involved in autoregulation of the nodulation signal. Our results provide information on the LRR-RLK gene family in M. truncatula and serve as a guide for functional research of the LRR-RLKs.

The Interaction Network and Signaling Specificity of Two-Component System in Arabidopsis

Two-component systems (TCS) in plants have evolved into a more complicated multi-step phosphorelay (MSP) pathway, which employs histidine kinases (HKs), histidine-containing phosphotransfer proteins (HPts), and response regulators (RRs) to regulate various aspects of plant growth and development. How plants perceive the external signals, then integrate and transduce the secondary signals specifically to the desired destination, is a fundamental characteristic of the MSP signaling network. The TCS elements involved in the MSP pathway and molecular mechanisms of signal transduction have been best understood in the model plant Arabidopsis thaliana. In this review, we focus on updated knowledge on TCS signal transduction in Arabidopsis. We first present a brief description of the TCS elements; then, the protein–protein interaction network is established. Finally, we discuss the possible molecular mechanisms involved in the specificity of the MSP signaling at the mRNA and protein levels.

Genome-wide characterization of two-component system (TCS) genes in melon (Cucumis melo L.)

To better understand cytokinin signaling in melon (Cucumis melo L.), one of the most important fruit crops in the Cucurbitaceae family, we identified and characterized melon two-component system (TCS) genes in this study. The results showed that there were 51 genes encoding putative TCS proteins in melon, and these TCS genes were classified into 3 subgroups, with 17 HK(L)s (histidine kinase/histidine-kinase like; 9 HKs and 8 HKLs), 9 HPs (histidine phosphotransfer proteins; 6 authentic and 3 pseudo), and 25 RRs (response regulators; 8 Type-A, 11 Type-B and 6 pseudo). The identity values of these cytokinin signaling proteins were revealed by analyzing their conserved motifs, domains and amino acid sequences. By analyzing TCS genes in different plant species, we found that melon HK(L)s, HPs and RRs had closer phylogenetic relationships with cucumber genes than with the genes of other plants, and the expansion of melon cytokinin signaling genes might be attributed to segmental duplication events. Analysis of the putative promoter regions (2-kb upstream regions of the start codon) revealed the enrichment of stress- and hormone-response cis-elements. The involvement of these putative TCS genes in melon cytokinin signaling was further supported by qRT-PCR data.

Genome-wide identification and characterization of TCP family genes in Brassica juncea var. tumida

Background

Teosinte branched1/Cycloidea/proliferating cell factors (TCPs) are plant-specific transcription factors widely involved in leaf development, flowering, shoot branching, the circadian rhythm, hormone signaling, and stress responses. However, the TCP function in Brassica juncea var. tumida, the tumorous stem mustard, has not yet been reported. This study identified and characterized the entire TCP family members in B. juncea var. tumida.

Methods

We identified 62 BjTCP genes from the B. juncea var. tumida genome and analyzed their phylogenetic relationship, gene structure, protein motifs, chromosome location, and expression profile in different tissues.

Results

Of the 62 BjTCP genes we identified in B. juncea var. tumida, containing 34 class I and 28 class II subfamily members, 61 were distributed on 18 chromosomes. Gene structure and conserved motif analysis showed that the same clade genes displayed a similar exon/intron gene structure and conserved motifs. Cis-acting element results showed that the same clade genes also had a similar cis-acting element; however, subtle differences implied a different regulatory pathway. The BjTCP18s members were low-expressed in Dayejie strains and the unswelling stage of Yonganxiaoye strains. Treatment with gibberellin (GA) and salicylic acid (SA) showed that GA and SA affect the expression levels of multiple TCP genes.

Conclusion

We performed the first genome-wide analysis of the TCP gene family of B. juncea var. tumida. Our results have provided valuable information for understanding the classification and functions of TCP genes in B. juncea var. tumida.

Heterologous Expression of a Soybean Gene RR34 Conferred Improved Drought Resistance of Transgenic Arabidopsis

Two-component systems (TCSs) have been identified as participants in mediating plant response to water deficit. Nevertheless, insights of their contribution to plant drought responses and associated regulatory mechanisms remain limited. Herein, a soybean response regulator (RR) gene RR34, which is the potential drought-responsive downstream member of a TCS, was ectopically expressed in the model plant Arabidopsis for the analysis of its biological roles in drought stress response. Results from the survival test revealed outstanding recovery ratios of 52%-53% in the examined transgenic lines compared with 28% of the wild-type plants. Additionally, remarkedly lower water loss rates in detached leaves as well as enhanced antioxidant enzyme activities of catalase and superoxide dismutase were observed in the transgenic group. Further transcriptional analysis of a subset of drought-responsive genes demonstrated higher expression in GmRR34-transgenic plants upon exposure to drought, including abscisic acid (ABA)-related genes NCED3, OST1, ABI5, and RAB18. These ectopic expression lines also displayed hypersensitivity to ABA treatment at germination and post-germination stages. Collectively, these findings indicated the ABA-associated mode of action of GmRR34 in conferring better plant performance under the adverse drought conditions.

Comparative analysis of basic helix-loop-helix gene family among Brassica oleracea, Brassica rapa, and Brassica napus

Background

The basic helix–loop–helix (bHLH) is the second largest gene family in the plant, some members play important roles in pistil development and response to drought, waterlogging, cold stress and salt stress. The bHLH gene family has been identified in many species, except for Brassica oleracea and B. napus thus far. This study aims to identify the bHLH family members in B. oleracea, B. rapa and B. napus, and elucidate the expression, duplication, phylogeny and evolution characters of them.

Result

A total of 268 bHLH genes in B. oleracea, 440 genes in B. napus, and 251 genes in B. rapa, including 21 new bHLH members, have been identified. Subsequently, the analyses of the phylogenetic trees, conserved motifs and gene structures showed that the members in the same subfamily were highly conserved. Most Ka/Ks values of homologous gene were < 1, which indicated that these genes suffered from strong purifying selection for retention. The retention rates of BrabHLH and BolbHLH genes were 51.6 and 55.1%, respectively. The comparative expression patterns between B. rapa and B. napus showed that they had similar expression patterns in the root and contrasting patterns in the stems, leaves, and reproductive tissues. In addition, there were 41 and 30 differential expression bHLH genes under the treatments of ABA and JA, respectively, and the number of down regulation genes was significantly more than up regulation genes.

Conclusion

In the present study, we identified and performed the comparative genomics analysis of bHLH gene family among B. oleracea, B. rapa and B. napus, and also investigated their diversity. The expression patterns between B. rapa and B. napus shows that they have the similar expression pattern in the root and opposite patterns in the stems, leaves, and reproduction tissues. Further analysis demonstrated that some bHLH gene members may play crucial roles under the abiotic and biotic stress conditions. This is the first to report on the bHLH gene family analysis in B. oleracea and B. napus, which can offer useful information on the functional analysis of the bHLH gene in plants.

Genome-wide identification and interactome analysis of members of two-component system in Banana

BACKGROUND

Ethylene signal transduction in plants is conducted by the two-component system (TCS) which consists of histidine kinase (HK), histidine phosphotransferase (HPT) and response regulators (RRs). This system plays an important role in signal transduction during various cellular processes, including fruit ripening and response to multiple environmental cues. Though members of TCS have been identified in a few plants, no detailed analysis has been carried out in banana.

RESULTS

Through genome-wide analysis, we identified a total of 80 (25 HK, 10 HPT and 45 RR) and 72 (25 HK, 5 HPT and 42 RR) TCS genes in Musa acuminata and Musa balbisiana respectively. The analysis of identified genes revealed that most of the genes are highly conserved however; there are subtle divergences among various members. Comparative expression analysis revealed an involvement of a set of TCS members during banana fruit ripening. Co-expression network analysis identified a working TCS module with direct interactions of HK-HPT and RR members. The molecular dynamics analysis of TCS module showed a significant change in structural trajectories of TCS proteins in the presence of ethylene. Analysis suggests possible interactions between the HK-HPTs and RRs as well as other members leading to banana fruit ripening.

CONCLUSIONS

In this study, we identified and compared the members of TCS gene family in two banana species and showed their diversity, within groups on the basis of whole-genome duplication events. Our analysis showed that during banana fruit ripening TCS module plays a crucial role. We also demonstrated a possible interaction mechanism of TCS proteins in the presence and absence of ethylene by molecular dynamics simulations. These findings will help in understanding the functional mechanism of TCS proteins in plants in different conditions.

Comprehensive genomic survey, structural classification and expression analysis of C2H2 zinc finger protein gene family in Brassica rapa L

C2H2 zinc finger protein (ZFP) genes have been extensively studied in many organisms and can function as transcription factors and be involved in many biological processes including plant growth and development and stress responses. In the current study, a comprehensive genomics analysis of the C2H2-ZFP genes in B. rapa was performed. A total of 301 B. rapa putative C2H2-ZFP (BrC2H2-ZFP) genes were identified from the available Brassica genome databases, and further characterized through analysis of conserved amino acid residues in C2H2-ZF domains and their organization, subcellular localization, phylogeny, additional domain, chromosomal location, synteny relationship, Ka/Ks ratio, and expression pattern. We also analyzed the expression patterns of eight B. rapa C2H2-ZFP genes under salt and drought stress conditions by using qRT-PCR technique. Our results showed that about one-third of these B. rapa C2H2-ZFP genes were originated from segmental duplication caused by the WGT around 13 to 17 MYA, one-third of them were highly and consecutively expressed in all tested tissues, and 92% of them were located in nucleus by prediction supporting then their functional roles as transcription factors, of which some may play important roles in plant growth and development. The Ka/Ks ratios of 264 orthologous C2H2-ZFP gene pairs between A. thaliana and B. rapa were all, except two, inferior to 1 (varied from 0.0116 to 1.4919, with an average value of 0.3082), implying that these genes had mainly experienced purifying selection during species evolution. The estimated divergence times of the same set of gene pairs ranged from 6.23 to 38.60 MY, with an average value of 18.29 MY, indicating that these gene members have undergone different selective pressures resulting in different evolutionary rates during species evolution. In addition, a few of these B. rapa C2H2-ZFPs were shown to be involved in stress responses in a similar way as their orthologs in A. thaliana. Comparison between A. thaliana and B. rapa orthologous C2H2-ZFP genes showed that the majority of these C2H2-ZFP gene members encodes proteins with conserved subcellular localization and functional domains between the two species but differed in their expression patterns in five tissues or organs. Thus, our study provides valuable information for further functional determination of each C2H2-ZFP gene across the Brassica species, and may help to select the appropriate gene targets for further in-depth studies, and genetic engineering and improvement of Brassica crops.

Diversity and Evolution of Sensor Histidine Kinases in Eukaryotes

Histidine kinases (HKs) are primary sensor proteins that act in cell signaling pathways generically referred to as "two-component systems" (TCSs). TCSs are among the most widely distributed transduction systems used by both prokaryotic and eukaryotic organisms to detect and respond to a broad range of environmental cues. The structure and distribution of HK proteins are now well documented in prokaryotes, but information is still fragmentary for eukaryotes. Here, we have taken advantage of recent genomic resources to explore the structural diversity and the phylogenetic distribution of HKs in the prominent eukaryotic supergroups. Searches of the genomes of 67 eukaryotic species spread evenly throughout the phylogenetic tree of life identified 748 predicted HK proteins. Independent phylogenetic analyses of predicted HK proteins were carried out for each of the major eukaryotic supergroups. This allowed most of the compiled sequences to be categorized into previously described HK groups. Beyond the phylogenetic analysis of eukaryotic HKs, this study revealed some interesting findings: 1) characterization of some previously undescribed eukaryotic HK groups with predicted functions putatively related to physiological traits; 2) discovery of HK groups that were previously believed to be restricted to a single kingdom in additional supergroups, and 3) indications that some evolutionary paths have led to the appearance, transfer, duplication, and loss of HK genes in some phylogenetic lineages. This study provides an unprecedented overview of the structure and distribution of HKs in the Eukaryota and represents a first step toward deciphering the evolution of TCS signaling in living organisms.

Genome-Wide Identification, Molecular Evolution, and Expression Profiling Analysis of Pectin Methylesterase Inhibitor Genes in Brassica campestris ssp. chinensis

Pectin methylesterase inhibitor genes (PMEIs) are a large multigene family and play crucial roles in cell wall modifications in plant growth and development. Here, a comprehensive analysis of the PMEI gene family in Brassicacampestris, an important leaf vegetable, was performed. We identified 100 BrassicacampestrisPMEI genes (BcPMEIs), among which 96 BcPMEIs were unevenly distributed on 10 chromosomes and nine tandem arrays containing 20 BcPMEIs were found. We also detected 80 pairs of syntenic PMEI orthologs. These findings indicated that whole-genome triplication (WGT) and tandem duplication (TD) were the main mechanisms accounting for the current number of BcPMEIs. In evolution, BcPMEIs were retained preferentially and biasedly, consistent with the gene balance hypothesis and two-step theory, respectively. The molecular evolution analysis of BcPMEIs manifested that they evolved through purifying selection and the divergence time is in accordance with the WGT data of B. campestris. To obtain the functional information of BcPMEIs, the expression patterns in five tissues and the cis-elements distributed in promoter regions were investigated. This work can provide a better understanding of the molecular evolution and biological function of PMEIs in B. campestris.