The grapevine kinome: annotation, classification and expression patterns in developmental processes and stress responses

Comparative genomics analysis provide insights into evolution and stress responses of Lhcb genes in Rosaceae fruit crops

BACKGROUND

Light-harvesting chlorophyll a/b b evelopment of higher plants and in response to abiotic stress. Previous works has demonstrated that that Lhcb genes were involved in the phytochrome regulation and responded to the different light and temperature conditions in Poaceae (such as maize). However, the evolution and functions of Lhcb genes remains poorly characterized in important Rosaceae species.

RESULTS

In this investigation, we conducted a genome-wide analysis and identified a total of 212 Lhcb genes across nine Rosaceae species. Specifically, we found 23 Lhcb genes in Fragaria vesca, 20 in Prunus armeniaca, 33 in Malus domestica 'Gala', 21 in Prunus persica, 33 in Rosa chinensis, 29 in Pyrus bretschneideri, 18 in Rubus occidentalis, 20 in Prunus mume, and 15 in Prunus salicina. Phylogenetic analysis revealed that the Lhcb gene family could be classified into seven major subfamilies, with members of each subfamily sharing similar conserved motifs. And, the functions of each subfamily was predicted based on the previous reports from other species. The Lhcb proteins were highly conserved within their respective subfamilies, suggesting similar functions. Interestingly, we observed similar peaks in Ks values (0.1-0.2) for Lhcb genes in apple and pear, indicating a recent whole genome duplication event (about 30 to 45 million years ago). Additionally, a few Lhcb genes underwent tandem duplication and were located across all chromosomes of nine species of Rosaceae. Furthermore, the analysis of the cis-acting elements in the 2000 bp promoter region upstream of the pear Lhcb gene revealed four main categories: light response correlation, stress response correlation, hormone response correlation, and plant growth. Quantitative expression analysis demonstrated that Lhcb genes exhibited tissue-specific expression patterns and responded differently to low-temperature stress in Rosaceae species.

CONCLUSIONS

These findings shed light on the evolution and phylogeny of Lhcb genes in Rosaceae and highlight the critical role of Lhcb in pear's response to low temperatures. The results obtained provide valuable insights for further investigations into the functions of Lhcb genes in Rosaceae, and these functional genes will be used for further fruit tree breeding and improvement to cope with the current climate changes.

Genome-wide characterization of the common bean kinome: Catalog and insights into expression patterns and genetic organization

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and is the core regulator of cellular signaling. Even considering this substantial importance, the kinome of common bean (Phaseolus vulgaris) has not been profiled yet. Here, we identified and characterised the complete set of kinases of common bean, performing an in-depth investigation with phylogenetic analyses and measurements of gene distribution, structural organization, protein properties, and expression patterns over a large set of RNA-Sequencing data. Being composed of 1,203 PKs distributed across all P. vulgaris chromosomes, this set represents 3.25% of all predicted proteins for the species. These PKs could be classified into 20 groups and 119 subfamilies, with a more pronounced abundance of subfamilies belonging to the receptor-like kinase (RLK)-Pelle group. In addition to provide a vast and rich reservoir of data, our study supplied insights into the compositional similarities between PK subfamilies, their evolutionary divergences, highly variable functional profile, structural diversity, and expression patterns, modeled with coexpression networks for investigating putative interactions associated with stress response.

Trends in the evolution of intronless genes in Poaceae

Intronless genes (IGs), which are a feature of prokaryotes, are a fascinating group of genes that are also present in eukaryotes. In the current study, a comparison of Poaceae genomes revealed that the origin of IGs may have involved ancient intronic splicing, reverse transcription, and retrotranspositions. Additionally, IGs exhibit the typical features of rapid evolution, including recent duplications, variable copy numbers, low divergence between paralogs, and high non-synonymous to synonymous substitution ratios. By tracing IG families along the phylogenetic tree, we determined that the evolutionary dynamics of IGs differed among Poaceae subfamilies. IG families developed rapidly before the divergence of Pooideae and Oryzoideae and expanded slowly after the divergence. In contrast, they emerged gradually and consistently in the Chloridoideae and Panicoideae clades during evolution. Furthermore, IGs are expressed at low levels. Under relaxed selection pressure, retrotranspositions, intron loss, and gene duplications and conversions may promote the evolution of IGs. The comprehensive characterization of IGs is critical for in-depth studies on intron functions and evolution as well as for assessing the importance of introns in eukaryotes.

Dehydration response in Stylosanthes scabra: Transcriptional, biochemical, and physiological modulations

Stylosanthes scabra, popularly known as stylo, is native to the Brazilian Caatinga semiarid region and stands out as a drought-tolerant shrub forage crop. This work provides information about the plant response during the first 48 h of water deficit, followed by a rehydration treatment. Besides root transcriptomics data, 13 physiological or biochemical parameters were scrutinized. Additionally, RNA-Seq annotated transcripts not associated with the "Viridiplantae" clade were taxonomically categorized. It was found that S. scabra quickly perceives and recovers from the oscillations of the imposed water regime. Physiologically, mechanisms that minimize evapotranspiration or protect the photosynthetic apparatus stood out. Biochemically, it was found that the root tissue invests in synthesizing compounds that can act as osmolytes (proline and sugars), emphasizing the importance of osmoregulation to water deficit acclimation. Consistently, transcriptome and qPCR analyses showed that a set of enriched biological processes with upregulated (UR) transcripts were involved in protective functions against reactive oxygen species or encoding enzymes of important metabolic pathways, which might contribute to S. scabra response to water deficit. Additionally, several UR kinases and transcription factors were identified. Finally, in an innovative approach, some naturally occurring microbial groups (such as Schizosaccharomyces, Bradyrhizobium, etc.) were identified in the S. scabra roots. This study reveals insights into the physiological, biochemical, and molecular mechanisms underlying the S. scabra response to water deficit and provides candidate genes that may be useful in developing drought-tolerant crop varieties through biotechnological applications.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

The emergence and evolution of intron-poor and intronless genes in intron-rich plant gene families

Eukaryotic genes can be classified into intronless (no introns), intron-poor (three or fewer introns per gene) or intron-rich. Early eukaryotic genes were mostly intron-rich, and their alternative splicing into multiple transcripts, giving rise to different proteins, might have played pivotal roles in adaptation and evolution. Interestingly, extant plant genomes contain many gene families with one or sometimes few sub-families with genes that are intron-poor or intronless, and it remains unknown when and how these intron-poor or intronless genes have originated and evolved, and what their possible functions are. In this study, we identified 33 such gene families that contained intronless and intron-poor sub-families. Intronless genes seemed to have first emerged in early land plant evolution, while intron-poor sub-families seemed first to have appeared in green algae. In contrast to intron-rich genes, intronless genes in intron-poor sub-families occurred later, and were subject to stronger functional constraints. Based on RNA-seq analyses in Arabidopsis and rice, intronless or intron-poor genes in AP2, EF-hand_7, bZIP, FAD_binding_4, STE_STE11, CAMK_CAMKL-CHK1 and C2 gene families were more likely to play a role in response to drought and salt stress, compared with intron-rich genes in the same gene families, whereas intronless genes in the B_lectin and S_locus_glycop gene family were more likely to participate in epigenetic processes and plant development. Understanding the origin and evolutionary trajectory, as well as the potential functions, of intronless and intron-poor sub-families provides further insight into plant genome evolution and the functional divergence of genes.

The Wild Sugarcane and Sorghum Kinomes: Insights Into Expansion, Diversification, and Expression Patterns

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and the core regulator of cellular signaling. Despite this substantial importance, the kinomes of sugarcane and sorghum have not been profiled. Here, we identified and profiled the complete kinomes of the polyploid Saccharum spontaneum (Ssp) and Sorghum bicolor (Sbi), a close diploid relative. The Sbi kinome was composed of 1,210 PKs; for Ssp, we identified 2,919 PKs when disregarding duplications and allelic copies, and these were related to 1,345 representative gene models. The Ssp and Sbi PKs were grouped into 20 groups and 120 subfamilies and exhibited high compositional similarities and evolutionary divergences. By utilizing the collinearity between the species, this study offers insights into Sbi and Ssp speciation, PK differentiation and selection. We assessed the PK subfamily expression profiles via RNA-Seq and identified significant similarities between Sbi and Ssp. Moreover, coexpression networks allowed inference of a core structure of kinase interactions with specific key elements. This study provides the first categorization of the allelic specificity of a kinome and offers a wide reservoir of molecular and genetic information, thereby enhancing the understanding of Sbi and Ssp PK evolutionary history.

Confirmation and Fine Mapping of the Resistance Locus Ren9 from the Grapevine Cultivar 'Regent'

Grapevine (Vitis vinifera ssp. vinifera) is a major fruit crop with high economic importance. Due to its susceptibility towards fungal and oomycete pathogens such as Erysiphe necator and Plasmopara viticola, the causal agents of powdery and downy mildew (PM and DM, respectively), grapevine growers annually face a major challenge in coping with shortfalls of yield caused by these diseases. Here we report the confirmation of a genetic resource for grapevine resistance breeding against PM. During the delimitation process of Ren3 on chromosome 15 from the cultivar ‘Regent’, a second resistance-encoding region on chromosome 15 termed Ren9 was characterized. It mediates a trailing necrosis associated with the appressoria of E. necator and restricts pathogen growth. In this study, we confirm this QTL in a related mapping population of ‘Regent’ × ‘Cabernet Sauvignon’. The data show that this locus is located at the upper arm of chromosome 15 between markers GF15-58 (0.15 Mb) and GF15-53 (4 Mb). The efficiency of the resistance against one of the prominent European PM isolates (EU-B) is demonstrated. Based on fine-mapping and literature knowledge we propose two possible regions of interest and supply molecular markers to follow both regions in marker-assisted selection.

Grapevine Responses to Heat Stress and Global Warming

The potential effects of the forthcoming climate change include the rising of the average annual temperature and the accumulation of extreme weather events, like frequent and severe heatwaves, a phenomenon known as global warming. Temperature is an important environmental factor affecting almost all aspects of growth and development in plants. The grapevine (Vitis spp.) is quite sensitive to extreme temperatures. Over the current century, temperatures are projected to continue rising with negative impacts on viticulture. These consequences range from short-term effects on wine quality to long-term issues such as the suitability of certain varieties and the sustainability of viticulture in traditional wine regions. Many viticultural zones, particularly in Mediterranean climate regions, may not be suitable for growing winegrapes in the near future unless we develop heat-stress-adapted genotypes or identify and exploit stress-tolerant germplasm. Grapevines, like other plants, have developed strategies to maintain homeostasis and cope with high-temperature stress. These mechanisms include physiological adaptations and activation of signaling pathways and gene regulatory networks governing heat stress response and acquisition of thermotolerance. Here, we review the major impacts of global warming on grape phenology and viticulture and focus on the physiological and molecular responses of the grapevine to heat stress.

The wild strawberry kinome: identification, classification and transcript profiling of protein kinases during development and in response to gray mold infection

Background

Protein kinases (PKs) play an important role in signaling cascades and are one of the largest and most conserved protein super families in plants. Despite their importance, the woodland strawberry (Fragaria vesca) kinome and expression patterns of PK genes remain to be characterized.

Results

Here, we report on the identification and classification of 954 Fragaria vesca PK genes, which were classified into nine groups and 124 gene families. These genes were distributed unevenly among the seven chromosomes, and the number of introns per gene varied from 0 to 47. Almost half of the putative PKs were predicted to localize to the nucleus and 24.6% were predicted to localize to the cell membrane. The expansion of the woodland strawberry PK gene family occurred via different duplication mechanisms and tandem duplicates occurred relatively late as compared to other duplication types. Moreover, we found that tandem and transposed duplicated PK gene pairs had undergone stronger diversifying selection and evolved relatively faster than WGD genes. The GO enrichment and transcriptome analysis implicates the involvement of strawberry PK genes in multiple biological processes and molecular functions in differential tissues, especially in pollens. Finally, 109 PKs, mostly the receptor-like kinases (RLKs), were found transcriptionally responsive to Botrytis cinerea infection.

Conclusions

The findings of this research expand the understanding of the evolutionary dynamics of PK genes in plant species and provide a potential link between cell signaling pathways and pathogen attack.

Genome-wide analysis of lectin receptor-like kinases family from potato (Solanum tuberosum L.)

Lectin receptor-like kinases (LecRLKs) are involved in responses to diverse environmental stresses and pathogenic microbes. A comprehensive acknowledgment of the family members in potato (Solanum tuberosum) genome is largely limited until now. In total, 113 potato LecRLKs (StLecRLKs) were first identified, including 85 G-type, 26 L-type and 2 C-type members. Based on phylogenetic analysis, StLecRLKs were sub-grouped into seven clades, including C-type, L-type, G-I, G-II, G-III G-IV and G-V. Chromosomal distribution and gene duplication analysis revealed the expansion of StLecRLKs occurred majorly through tandem duplication although the whole-genome duplication (WGD)/segmental duplication events were found. Cis-elements in the StLecRLKs promoter region responded mainly to signals of defense and stress, phytohormone, biotic or abiotic stress. Moreover, expressional investigations indicated that the family members of the clades L-type, G-I, G-IV and G-V were responsive to both bacterial and fungal infection. Based on qRT-PCR analysis, the expressions of PGSC0003DMP400055136 and PGSC0003DMP400067047 were strongly induced in all treatments by both Fusarium sulphureum (Fs) and Phytophthora infestans (Pi) inoculation. The present study provides valuable information for LecRLKs gene family in potato genome, and establishes a foundation for further research into the functional analysis.

The grapevine NaE sodium exclusion locus encodes sodium transporters with diverse transport properties and localisation

Grapevine (Vitis vinifera L.) is a valuable crop for human consumption and wine production, and is prone to suffering from salinity stress in arid regions or when exposed to low quality irrigation water. A previous study identified a quantitative trait locus (QTL) NaE, containing six High-affinity Potassium Transporter 1 genes, that was associated with shoot Na⁺ exclusion in grapevine. While HKT1;1 was predicted to be the most likely gene within this QTL to encode for this important salinity tolerance sub-trait, four other HKTs within the QTL remained uncharacterised; VviHKT1;2 encodes a truncated transcript unlikely to form a functional transporter. In this study, two allelic variants for each of VviHKT1;6, VviHKT1;7 and VviHKT1;8 from the heterozygous grapevine variety Cabernet Sauvignon were functionally characterised. Using the Xenopus laevis oocyte heterologous expression system, as well as transient expression in tobacco leaves, we found that the VviHKT1;6 and VviHKT1;7 alleles encoded plasma membrane localised proteins that facilitated significant non-rectifying Na⁺ transport. Conversely, proteins encoded by the VviHKT1;8 alleles were inwardly-rectifying, weak Na⁺ transporters that localised to intracellular organelles. Mining of previous RNA-seq gene expression data suggested that VviHKT1;6-8 are weakly expressed in grapevine roots, flower buds, and seeds under normal conditions and different nutrient regimes. We propose that VviHKT1;6 and VviHKT1;7 are likely to have a less significant role in grapevine leaf Na⁺ exclusion than VviHKT1;1, and that VviHKT1;8 is involved in endomembrane Na⁺ transport.

Selection of reliable reference genes for quantitative RT-PCR in garlic under salt stress

Quantitative real-time reverse-transcriptase PCR (qRT-PCR) has been frequently used for detecting gene expression. To obtain reliable results, selection of suitable reference genes is a fundamental and necessary step. Garlic (Allium sativum), a member from Alliaceae family, has been used both as a food flavoring and as a traditional medicine. In the present study, garlic plants were exposed to salt stress (200 mM NaCl) for 0, 1, 4 and 12 h, and garlic roots, bulbs, and leaves were harvested for subsequent analysis. The expression stability of eight candidate reference genes, eukaryotic translation initiation factor 4α (eIF-4α), actin (ACTIN), tubulin β-7 (TUB7), TAP42-interacting protein of 41 kDa (TIP41), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), SAND family protein (SAND), elongation factor 1 alpha (EF-1α), and protein phosphatase 2A (PP2A) were evaluated by geNorm, NormFinder, and BestKeeper. All genes tested displayed variable expression profiles under salt stress. In the leaf and root group, ACTIN was the best reference gene for normalizing gene expression. In garlic clove, ACTIN and SAND were the least variable, and were suitable for gene expression studies under salt stress; these two genes also performed well in all samples tested. Based on our results, we recommend that it is essential to use specific reference genes in different situations to obtain accurate results. Using a combination of multiple stable reference genes, such as ACTIN and SAND, to normalize gene expression is encouraged. The results from the study will be beneficial for accurate determination of gene expression in garlic and other plants.

Genome-Wide Identification, Characterization, and Expression Analysis of the Grapevine Superoxide Dismutase (SOD) Family

Superoxide dismutase (SOD) is an essential enzyme of the plant antioxidant system that responds to oxidative damage caused by adverse conditions. However, little is known about the SOD gene family in Vitis vinifera (Vv). In the present study, ten SOD genes, including 6 copper/zinc SODs, 2 iron SODs, and 2 manganese SODs, were identified in the grapevine genome where they were unevenly distributed on 12 chromosomes. Ten VvSOD genes were divided into three main groups based on phylogenetic analysis, subcellular localization, and the distribution of conserved protein motifs. Additionally, many cis-elements related to different stresses were found in the promoters of the 10 VvSOD genes. Syntenic analysis revealed that VvMSD1 and VvMSD2 were derived from segmental duplication, and VvCSD4 and VvCSD5 belong to a pair of tandemly duplicated genes. Gene expression analysis based on microarray data showed that the 10 VvSOD genes were expressed in all the tested tissues. Interestingly, the segmentally duplicated gene pair (VvMSD1 and VvMSD2) exhibited differential expression patterns in various organs. In contrast, the tandemly duplicated gene pair (VvCSD4 and VvCSD5) displayed similar expression patterns in the tested organs. Our results provide a basis for further functional research on the SOD gene family in grapevine.

The kinome of pineapple: catalog and insights into functions in crassulacean acid metabolism plants

BACKGROUND

Crassulacean acid metabolism (CAM) plants use water 20-80% more efficiently by shifting stomata opening and primary CO₂ uptake and fixation to the nighttime. Protein kinases (PKs) play pivotal roles in this biological process. However, few PKs have been functionally analyzed precisely due to their abundance and potential functional redundancy (caused by numerous gene duplications).

RESULTS

In this study, we systematically identified a total of 758 predicted PK genes in the genome of a CAM plant, pineapple (Ananas comosus). The pineapple kinome was classified into 20 groups and 116 families based on the kinase domain sequences. The RLK was the largest group, containing 480 members, and over half of them were predicted to locate at the plasma membrane. Both segmental and tandem duplications make important contributions to the expansion of pineapple kinome based on the synteny analysis. Ka/Ks ratios showed all of the duplication events were under purifying selection. The global expression analysis revealed that pineapple PKs exhibit different tissue-specific and diurnal expression patterns. Forty PK genes in a cluster performed higher expression levels in green leaf tip than in white leaf base, and fourteen of them had strong differential expression patterns between the photosynthetic green leaf tip and the non-photosynthetic white leaf base tissues.

CONCLUSIONS

Our findings provide insights into the evolution and biological function of pineapple PKs and a foundation for further functional analysis of PKs in CAM plants. The gene duplication, expression, and coexpression analysis helped us to rapidly identify the key candidates in pineapple kinome, which may play roles in the carbon fixation process in pineapple and help engineering CAM pathway into C3 crops for improved drought tolerance.