Comparative Genomic Analysis Reveals Extensive Genetic Variations of WRKYs in Solanaceae and Functional Variations of CaWRKYs in Pepper

As a conserved protein family, WRKY has been shown to be involved in multiple biological processes in plants. However, the mechanism of functional diversity for WRKYs in pepper has not been well elucidated. Here, a total of 223 WRKY members from solanaceae crops including pepper, tomato and potato, were analyzed using comparative genomics. A tremendous genetic variation among WRKY members of different solanaceous plants or groups was demonstrated by the comparison of some WRKY features, including number/size, group constitution, gene structure, and domain composition. The phylogenetic analysis showed that except for the known WRKY groups (I, IIa/b/c/d/e and III), two extra WRKY subgroups specifically existed in solanaceous plants, which were named group IIf and group IIg in this study, and their genetic variations were also revealed by the characteristics of some group IIf and IIg WRKYs. Except for the extensive genetic variations, certain degrees of conservatism for solanaceae WRKYs were also revealed. Moreover, the variant zinc-finger structure (CX4,7CX22-24HXC) in group III of solanaceae WRKYs was identified. Expression profiles of CaWRKY genes suggested their potential roles in pepper development and stress responses, and demonstrated a functional division pattern for pepper CaWRKYs. Furthermore, functional analysis using virus induced gene silencing (VIGS) revealed critical roles of two CaWRKYs (CaWRKY45 and CaWRKY58) in plant responses to disease and drought, respectively. This study provides a solid foundation for further dissection of the evolutionary and functional diversity of solanaceae WRKYs in crop plants

The PAP Gene Family in Tomato: Comprehensive Comparative Analysis, Phylogenetic Relationships and Expression Profiles

Purple acid phosphatase (PAP) plays a vital role in plant phosphate acquisition and utilization, as well as cell wall synthesis and redox reactions. In this study, comprehensive comparative analyses of PAP genes were carried out using the integration of phylogeny, chromosomal localization, intron/exon structural characteristics, and expression profiling. It was shown that the number of introns of the PAP genes, which were distributed unevenly on 12 chromosomes, ranged from 1 to 12. These findings pointed to the existence of complex structures. Phylogenetic analyses revealed that PAPs from tomato, rice, and Arabidopsis could be divided into three groups (Groups I, II, and III). It was assumed that the diversity of these PAP genes occurred before the monocot–dicot split. RNA-seq analysis revealed that most of the genes were expressed in all of the tissues analyzed, with the exception of SlPAP02, SlPAP11, and SlPAP14, which were not detected. It was also found that expression levels of most of the SlPAP gene family of members were changed under phosphorus stress conditions, suggesting potential functional diversification. The findings of this work will help us to achieve a better insight into the function of SlPAP genes in the future, as well as enhance our understanding of their evolutionary relationships in plants

Molecular Evolution and Functional Divergence of Stress-Responsive Cu/Zn Superoxide Dismutases in Plants

Superoxide dismutases (SODs), a family of antioxidant enzymes, are the first line of defense against oxidative damage and are ubiquitous in every cell of all plant types. The Cu/Zn SOD, one of three types of SODs present in plant species, is involved in many of the biological functions of plants in response to abiotic and biotic stresses. Here, we carried out a comprehensive analysis of the Cu/Zn SOD gene family in different plant species, ranging from lower plants to higher plants, and further investigated their organization, sequence features, and expression patterns in response to biotic and abiotic stresses. Our results show that plant Cu/Zn SODs can be divided into two subfamilies (group I and group II). Group II appeared to be conserved only as single- or low-copy genes in all lineages, whereas group I genes underwent at least two duplication events, resulting in multiple gene copies and forming three different subgroups (group Ia, group Ib, and group Ic). We also found that, among these genes, two important events—the loss of introns and the loss of and variation in signal peptides—occurred over the long course of their evolution, indicating that they were involved in shifts in subcellular localization from the chloroplast to cytosol or peroxisome and underwent functional divergence. In addition, expression patterns of Cu/Zn SOD genes from Arabidopsis thaliana and Solanum lycopersicum were tested in different tissues/organs and developmental stages and under different abiotic stresses. The results indicate that the Cu/Zn SOD gene family possesses potential functional divergence and may play vital roles in ROS scavenging in response to various stresses in plants. This study will help establish a foundation for further understanding these genes’ function during stress responses

Divergent Retention of Sucrose Metabolism Genes after Whole Genome Triplication in the Tomato (<i>Solanum lycopersicum</i>)

Sucrose, the primary carbon transport mode and vital carbohydrate for higher plants, significantly impacts plant growth, development, yield, and quality formation. Its metabolism involves three key steps: synthesis, transport, and degradation. Two genome triplication events have occurred in Solanaceae, which have resulted in massive gene loss. In this study, a total of 48 and 65 genes from seven sucrose metabolism gene families in Vitis vinifera and Solanum lycopersicum were identified, respectively. The number of members comprising the different gene families varied widely. And there were significant variations in the pattern of gene duplication and loss in the tomato following two WGD events. Tandem duplication is a major factor in the expansion of the SWEET and Acid INV gene families. All the genes are irregularly distributed on the chromosomes, with the majority of the genes showing collinearity with the grape, particularly the CIN family. And the seven gene families were subjected to a purifying selection. The expression patterns of the different gene families exhibited notable variations. This study presents basic information about the sucrose metabolism genes in the tomato and grape, and paves the way for further investigations into the impact of SCT events on the phylogeny, gene retention duplication, and function of sucrose metabolism gene families in the tomato or Solanaceae, and the adaptive evolution of the tomato

Genome-wide identification and expression profiling of tomato Hsp20 gene family in response to biotic and abiotic stresses

The Hsp20 genes are involved in the response of plants to environment stresses including heat shock and also play a vital role in plant growth and development. They represent the most abundant small heat shock proteins (sHsps) in plants, but little is known about this family in tomato (Solanum lycopersicum), an important vegetable crop in the world. Here, we characterized heat shock protein 20 (SlHsp20) gene family in tomato through integration of gene structure, chromosome location, phylogenetic relationship and expression profile. Using bioinformatics-based methods, we identified at least 42 putative SlHsp20 genes in tomato. Sequence analysis revealed that most of SlHsp20 genes possessed no intron or a relatively short intron in length. Chromosome mapping indicated that inter-arm and intra-chromosome duplication events contributed remarkably to the expansion of SlHsp20 genes. Phylogentic tree of Hsp20 genes from tomato and other plant species revealed that SlHsp20 genes were grouped into 13 subfamilies, indicating that these genes may have a common ancestor that generated diverse subfamilies prior to the mono-dicot split. In addition, expression analysis using RNA-seq in various tissues and developmental stages of cultivated tomato and the wild relative Solanum pimpinellifolium revealed that most of these genes (83%) were expressed in at least one stage from at least one genotype. Out of 42 genes, 4 genes were expressed constitutively in almost all the tissues analyzed, implying that these genes might have specific housekeeping function in tomato cell under normal growth conditions. Two SlHsp20 genes displayed differential expression levels between cultivated tomato and S. pimpinellifolium in vegetative (leaf and root) and reproductive organs (floral bud and flower), suggesting inter-species diversification for functional specialization during the process of domestication. Based on genome-wide microarray analysis, we showed that the transcript levels of SlHsp20 genes could be induced profusely by abiotic and biotic stresses such as heat, drought, salt, Botrytis cinerea and Tomato Spotted Wilt Virus, indicating their potential roles in mediating the response of tomato plants to environment stresses. In conclusion, these results provide valuable information for elucidating the evolutionary relationship of Hsp20 gene family and functional characterization of the SlHsp20 gene family in the future

Analysis of TIR- and non-TIR-NBS-LRR disease resistance gene analogous in pepper: characterization, genetic variation, functional divergence and expression patterns

<p>Abstract</p> <p>Background</p> <p>Pepper (<it>Capsicum annuum</it> L.) is one of the most important vegetable crops worldwide. However, its yield and fruit quality can be severely threatened by several pathogens. The plant nucleotide-binding site (NBS)-leucine-rich repeat (LRR) gene family is the largest class of known disease resistance genes (R genes) effective against such pathogens. Therefore, the isolation and identification of such R gene homologues from pepper will provide a critical foundation for improving disease resistance breeding programs.</p> <p>Results</p> <p>A total of 78 R gene analogues (CaRGAs) were identified in pepper by degenerate PCR amplification and database mining. Phylogenetic tree analysis of the deduced amino acid sequences for 51 of these CaRGAs with typically conserved motifs ( P-loop, kinase-2 and GLPL) along with some known R genes from <it>Arabidopsis</it> and tomato grouped these CaRGAs into the non-Toll interleukin-1 receptor (TIR)-NBS-LRR (CaRGAs I to IV) and TIR-NBS-LRR (CaRGAs V to VII) subfamilies. The presence of consensus motifs (i.e. P-loop, kinase-2 and hydrophobic domain) is typical of the non-TIR- and TIR-NBS-LRR gene subfamilies. This finding further supports the view that both subfamilies are widely distributed in dicot species. Functional divergence analysis provided strong statistical evidence of altered selective constraints during protein evolution between the two subfamilies. Thirteen critical amino acid sites involved in this divergence were also identified using DIVERGE version 2 software. Analyses of non-synonymous and synonymous substitutions per site showed that purifying selection can play a critical role in the evolutionary processes of non-TIR- and TIR-NBS-LRR RGAs in pepper. In addition, four specificity-determining positions were predicted to be responsible for functional specificity. qRT-PCR analysis showed that both salicylic and abscisic acids induce the expression of CaRGA genes, suggesting that they may primarily be involved in defence responses by activating signaling pathways.</p> <p>Conclusion</p> <p>The identified CaRGAs are a valuable resource for discovering R genes and developing RGA molecular markers for genetic map construction. They will also be useful for improving disease resistance in pepper. The findings of this study provide a better understanding of the evolutionary mechanisms that drive the functional diversification of non-TIR- and TIR-NBS-LRR R genes in pepper.</p

The SOD gene family in tomato: identification, phylogenetic relationships and expression patterns

Superoxide dismutases (SODs) are critical antioxidant enzymes that protect organisms from reactive oxygen species (ROS) caused by adverse conditions, and have been widely found in the cytoplasm, chloroplasts, and mitochondria of eukaryotic and prokaryotic cells. Tomato (Solanum lycopersicum L.) is an important economic crop and is cultivated worldwide. However, abiotic and biotic stresses severely hinder growth and development of the plant, which affects the production and quality of the crop. To reveal the potential roles of SOD genes under various stresses, we performed a systematic analysis of the tomato SOD gene family and analyzed the expression patterns of SlSOD genes in response to abiotic stresses at the whole-genome level. The characteristics of the SlSOD gene family were determined by analyzing gene structure, conserved motifs, chromosomal distribution, phylogenetic relationships, and expression patterns. We determined that there are at least nine SOD genes in tomato, including four Cu/ZnSODs, three FeSODs, and one MnSOD, and they are unevenly distributed on 12 chromosomes. Phylogenetic analyses of SOD genes from tomato and other plant species were separated into two groups with a high bootstrap value, indicating that these SOD genes were present before the monocot-dicot split. Additionally, many cis-elements that respond to different stresses were found in the promoters of nine SlSOD genes. Gene expression analysis based on RNA-seq data showed that most genes were expressed in all tested tissues, with the exception of SlSOD6 and SlSOD8, which were only expressed in young fruits. Microarray data analysis showed that most members of the SlSOD gene family were altered under salt- and drought-stress conditions. This genome-wide analysis of SlSOD genes helps to clarify the function of SlSOD genes under different stress conditions and provides information to aid in further understanding the evolutionary relationships of SOD genes in plants

Identification of Optimal Reference Genes for Normalization of qPCR Analysis during Pepper Fruit Development

Due to its high sensitivity and reproducibility, quantitative real-time PCR (qPCR) is practiced as a useful research tool for targeted gene expression analysis. For qPCR operations, the normalization with suitable reference genes (RGs) is a crucial step that eventually determines the reliability of the obtained results. Although pepper is considered an ideal model plant for the study of non-climacteric fruit development, at present no specific RG have been developed or validated for the qPCR analyses of pepper fruit. Therefore, this study aimed to identify stably expressed genes for their potential use as RGs in pepper fruit studies. Initially, a total of 35 putative RGs were selected by mining the pepper transcriptome data sets derived from the PGP (Pepper Genome Platform) and PGD (Pepper Genome Database). Their expression stabilities were further measured in a set of pepper (Capsicum annuum L. var. 007e) fruit samples, which represented four different fruit developmental stages (IM: Immature; MG: Mature green; B: Break; MR: Mature red) using the qPCR analysis. Then, based on the qPCR results, three different statistical algorithms, namely geNorm, Normfinder, and boxplot, were chosen to evaluate the expression stabilities of these putative RGs. It should be noted that nine genes were proven to be qualified as RGs during pepper fruit development, namely CaREV05 (CA00g79660); CaREV08 (CA06g02180); CaREV09 (CA06g05650); CaREV16 (Capana12g002666); CaREV21 (Capana10g001439); CaREV23 (Capana05g000680); CaREV26 (Capana01g002973); CaREV27 (Capana11g000123); CaREV31 (Capana04g002411); and CaREV33 (Capana08g001826). Further analysis based on geNorm suggested that the application of the two most stably expressed genes (CaREV05 and CaREV08) would provide optimal transcript normalization in the qPCR experiments. Therefore, a new and comprehensive strategy for the identification of optimal RGs was developed. This strategy allowed for the effective normalization of the qPCR analysis of the pepper fruit development at the whole pepper genome level. This study not only explored the optimal RGs specific for studying pepper fruit development, but also introduced a referable strategy of RG mining which could potentially be implicated in other plant species

The Tryptophan Decarboxylase in Solanum lycopersicum

Melatonin plays an important role in plant growth, development, and environmental stress. In this study, a systematic analysis of tomato tryptophan decarboxylase (SlTrpDC), which is the first enzyme of melatonin biosynthesis, was conducted by integrating structural features, phylogenetic relationships, an exon/intron feature, and a divergent expression profile. The results determined that the tomato genome encoded five members (SlTrpDC1-SlTrpDC5). The phylogenetic relationships indicated that gene expansion was proposed as the major mode of evolution of the TrpDC genes from the different plant algae species to the higher plants species. The analyses of the exon/intron configurations revealed that the intron loss events occurred during the structural evolution of the TrpDCs in plants. Additionally, the RNA-seq and qRT-PCR analysis revealed that the expression of the SlTrpDC3 was high in all of the tested tissues, while the SlTrpDC4 and SlTrpDC5 were not expressed. The expression patterns of the remaining two (SlTrpDC1 and SlTrpDC2) were tissue-specific, which indicated that these genes may play important roles within the different tissues. No expression difference was observed in the tomato plants in response to the biotic stresses. This study will expand the current knowledge of the roles of the TrpDC genes in tomato growth and development