Genome-wide investigation of WRKY gene family in pineapple: evolution and expression profiles during development and stress

Genome-Wide Systematic Characterization of the NPF Family Genes and Their Transcriptional Responses to Multiple Nutrient Stresses in Allotetraploid Rapeseed

NITRATE TRANSPORTER 1 (NRT1)/PEPTIDE TRANSPORTER (PTR) family (NPF) proteins can transport various substrates, and play crucial roles in governing plant nitrogen (N) uptake and distribution. However, little is known about the NPF genes in Brassica napus. Here, a comprehensive genome-wide systematic characterization of the NPF family led to the identification of 193 NPF genes in the whole genome of B. napus. The BnaNPF family exhibited high levels of genetic diversity among sub-families but this was conserved within each subfamily. Whole-genome duplication and segmental duplication played a major role in BnaNPF evolution. The expression analysis indicated that a broad range of expression patterns for individual gene occurred in response to multiple nutrient stresses, including N, phosphorus (P) and potassium (K) deficiencies, as well as ammonium toxicity. Furthermore, 10 core BnaNPF genes in response to N stress were identified. These genes contained 6–13 transmembrane domains, located in plasma membrane, that respond discrepantly to N deficiency in different tissues. Robust cis-regulatory elements were identified within the promoter regions of the core genes. Taken together, our results suggest that BnaNPFs are versatile transporters that might evolve new functions in B. napus. Our findings benefit future research on this gene family.

Genome-wide identification and characterization of DCL, AGO and RDR gene families in Saccharum spontaneum

RNA silencing is a conserved mechanism in eukaryotic organisms to regulate gene expression. Argonaute (AGO), Dicer-like (DCL) and RNA-dependent RNA polymerase (RDR) proteins are critical components of RNA silencing, but how these gene families’ functions in sugarcane were largely unknown. Most stress-resistance genes in modern sugarcane cultivars (Saccharum spp.) were originated from wild species of Saccharum, for example S. spontaneum. Here, we used genome-wide analysis and a phylogenetic approach to identify four DCL, 21 AGO and 11 RDR genes in the S. spontaneum genome (termed SsDCL, SsAGO and SsRDR, respectively). Several genes, particularly some of the SsAGOs, appeared to have undergone tandem or segmental duplications events. RNA-sequencing data revealed that four SsAGO genes (SsAGO18c, SsAGO18b, SsAGO10e and SsAGO6b) and three SsRDR genes (SsRDR2b, SsRDR2d and SsRDR3) tended to have preferential expression in stem tissue, while SsRDR5 was preferentially expressed in leaves. qRT-PCR analysis showed that SsAGO10c, SsDCL2 and SsRDR6b expressions were strongly upregulated, whereas that of SsAGO18b, SsRDR1a, SsRDR2b/2d and SsRDR5 was significantly depressed in S. spontaneum plants exposed to PEG-induced dehydration stress or infected with Xanthomonas albilineans, causal agent of leaf scald disease of sugarcane, suggesting that these genes play important roles in responses of S. spontaneum to biotic and abiotic stresses.

Genome-Wide Analysis, Characterization, and Expression Profile of the Basic Leucine Zipper Transcription Factor Family in Pineapple

This study identified 57 basic leucine zipper (bZIP) genes from the pineapple genome, and the analysis of these bZIP genes was focused on the evolution and divergence after multiple duplication events in relation to the pineapple genome fusion. According to bioinformatics analysis of a phylogenetic tree, the bZIP gene family was divided into 11 subgroups in pineapple, Arabidopsis, and rice; gene structure and conserved motif analyses showed that bZIP genes within the same subgroup shared similar intron-exon organizations and motif composition. Further synteny analysis showed 17 segmental duplication events with 27 bZIP genes. The study also analyzed the pineapple gene expression of bZIP genes in different tissues, organs, and developmental stages, as well as in abiotic stress responses. The RNA-sequencing data showed that AcobZIP57 was upregulated in all tissues, including vegetative and reproductive tissues. AcobZIP28 and AcobZIP43 together with the other 25 bZIP genes did not show high expression levels in any tissue. Six bZIP genes were exposed to abiotic stress, and the relative expression levels were detected by quantitative real-time PCR. A significant response was observed for AcobZIP24 against all kinds of abiotic stresses at 24 and 48 h in pineapple root tissues. Our study provides a perspective for the evolutionary history and general biological involvement of the bZIP gene family of pineapple, which laid the foundation for future functional characterization of the bZIP genes in pineapple.

Phylogenetic and Transcription Analysis of Hibiscus hamabo Sieb. et Zucc. WRKY Transcription Factors

WRKY transcription factors are known to play important roles in the regulation of various aspects of plant growth and development, including germination, stress resistance, and senescence. Nevertheless, there is little information about the WRKY genes in Hibiscus hamabo Sieb. et Zucc., an important semimangrove plant. In this study, HhWRKY genes in H. hamabo were identificated based on Illumina RNA-sequencing and isoform sequencing from salt-treated roots. Then phylogenetic analysis and conserved motif analysis of the WRKY family in H. hamabo and Arabidopsis thaliana were used to classify WRKY genes. Sixteen HhWRKY genes were selected from different groups to detect their expression patterns using real-time quantitative PCR in different organ (root, old leaf, tender leaf, receptacle, petal, or stamen) from 10-year-old H. hamabo plants grown in their natural environment and in seedlings with 8 to 10 true leaves challenged by phytohormone (salicylic acid, methyl jasmonate, or abscisic acid) and abiotic stress (drought, salt, or high temperature). As a result, the identified 78 HhWRKY genes were divided into two major groups and several subgroups based on their structural and phylogenetic features. Most transcripts of the selected 16 HhWRKY genes were more abundant in old than in tender leaves of H. hamabo. HhWRKY genes were regulated in reaction to abiotic stresses and phytohormone treatments and may participate in signaling networks to improve plant stress resistance. Some of HhWRKY genes behaved as would be predicted based on their homology with A. thaliana WRKY genes, but others showed divergent behavior. This systematic analysis lays the foundation for further identification of WRKY gene functions, with the aim of improving woody plants.

Heavy Metal Stress-Associated Proteins in Rice and Arabidopsis: Genome-Wide Identification, Phylogenetics, Duplication, and Expression Profiles Analysis

Heavy metal exposure is a serious environmental stress in plants. However, plants have evolved several strategies to improve their heavy metal tolerance. Heavy metal-associated proteins (HMPs) participate in heavy metal detoxification. Here, we identified 46 and 55 HMPs in rice and Arabidopsis, respectively, and named them OsHMP 1–46 and AtHMP 1–55 according to their chromosomal locations. The HMPs from both plants were divided into six clades based on the characteristics of their heavy metal-associated domains (HMA). The HMP gene structures and motifs varied greatly among the different classifications. The HMPs had high collinearity and were segmentally duplicated. A cis-element analysis revealed that the HMPs may be regulated by different transcription factors. An expression profile analysis disclosed that only eight OsHMPs were constitutive in rice tissues. Of these, the expression of OsHMP37 was far higher than that of the other seven genes while OsHMP28 was expressed exclusively in the roots. For Arabidopsis, nine AtHMPs presented with very high transcript levels in all organs. Most of the selected OsHMPs were differentially expressed in various tissues under different heavy metal stresses. Only OsHMP09, OsHMP18, and OsHMP22 showed higher expression levels in all tissues under different heavy metal stresses. In contrast, most of the selected AtHMPs had nearly constant expression levels in different tissues under various heavy metal stresses. The AtHMP20, AtHMP23, AtHMP25, AtHMP31, AtHMP35, AtHMP46 expression levels under different heavy metal stresses were higher in the leaves and roots. The foregoing discoveries elucidated HMP evolution in monocotyledonous and dicotyledonous plants and may helpful functionally characterize HMPs in the future.

Identification and expression analysis of WRKY gene family under drought stress in peanut (Arachis hypogaea L.)

WRKY transcription factors play crucial roles in regulation mechanism leading to the adaption of plants to the complex environment. In this study, AhWRKY family was comprehensively analyzed using bioinformatic approaches in combination with transcriptome sequencing data of the drought-tolerant peanut variety ‘L422’. A total of 158 AhWRKY genes were identified and named according to their distribution on the chromosomes. Based on the structural features and phylogenetic analysis of AhWRKY proteins, the AhWRKY family members were classified into three (3) groups, of which group II included five (5) subgroups. Results of structure and conserved motifs analysis for the AhWRKY genes confirmed the accuracy of the clustering analysis. In addition, 12 tandem and 136 segmental duplication genes were identified. The results indicated that segmental duplication events were the main driving force in the evolution of AhWRKY family. Collinearity analysis found that 32 gene pairs existed between Arachis hypogaea and two diploid wild ancestors (Arachis duranensis and Arachis ipaensis), which provided valuable clues for phylogenetic characteristics of AhWRKY family. Furthermore, 19 stress-related cis-acting elements were found in the promoter regions. During the study of gene expression level of AhWRKY family members in response to drought stress, 73 differentially expressed AhWRKY genes were obtained to have been influenced by drought stress. These results provide fundamental insights for further study of WRKY genes in peanut drought resistance.

Genome-wide investigation of WRKY transcription factors in sweet osmanthus and their potential regulation of aroma synthesis

WRKY transcription factors, one of the largest transcription factor families, play important roles in regulating the synthesis of secondary metabolites. In sweet osmanthus (Osmanthus fragrans), the monoterpenes have been demonstrated as the most important volatile compounds, and the W-box, which is the cognate binding site of WRKY transcription factors, could be identified in most of the terpene-synthesis-related genes' promoters. However, the role of the WRKY family in terpene synthesis in sweet osmanthus has rarely been examined. In this study, 154 WRKY genes with conserved WRKY domain were identified and classified into three groups. The group II was further divided into five subgroups, and almost all members of IId contained a plant zinc cluster domain. Eight OfWRKYs (OfWRKY7/19/36/38/42/84/95/139) were screened from 20 OfWRKYs for their flower-specific expression patterns in different tissues. Simultaneously, the expression patterns of OfWRKYs and emission patterns of volatile compounds during the flowering process were determined and gas chromatography-mass spectrometry results showed that monoterpenes, such as linalool and ocimene, accounted for the highest proportion, contributing to the floral scent of sweet osmanthus in two cultivars. In addition, correlation analysis revealed the expression patterns of OfWRKYs (OfWRKY7/19/36/139) were each correlated with distinct monoterpenes (linalool, linalool derivatives, ocimene and ocimene derivatives). Subcellular localization analysis showed that p35S::GFP-OfWRKY7/38/95/139 were localized in the nucleus and OfWRKY139 had very strong transactivation activity. Collectively, the results indicated potential roles of OfWRKY139 and OfWRKYs with plant zinc cluster domain in regulating synthesis of aromatic compounds in sweet osmanthus, laying the foundation for use of OfWRKYs to improve the aroma of ornamental plants.

Genome-wide identification and expression analysis of WRKY transcription factors in pearl millet (Pennisetum glaucum) under dehydration and salinity stress

Background

Plants have developed various sophisticated mechanisms to cope up with climate extremes and different stress conditions, especially by involving specific transcription factors (TFs). The members of the WRKY TF family are well known for their role in plant development, phytohormone signaling and developing resistance against biotic or abiotic stresses. In this study, we performed a genome-wide screening to identify and analyze the WRKY TFs in pearl millet (Pennisetum glaucum; PgWRKY), which is one of the most widely grown cereal crops in the semi-arid regions.

Results

A total number of 97 putative PgWRKY proteins were identified and classified into three major Groups (I-III) based on the presence of WRKY DNA binding domain and zinc-finger motif structures. Members of Group II have been further subdivided into five subgroups (IIa-IIe) based on the phylogenetic analysis. In-silico analysis of PgWRKYs revealed the presence of various cis-regulatory elements in their promoter region like ABRE, DRE, ERE, EIRE, Dof, AUXRR, G-box, etc., suggesting their probable involvement in growth, development and stress responses of pearl millet. Chromosomal mapping evidenced uneven distribution of identified 97 PgWRKY genes across all the seven chromosomes of pearl millet. Synteny analysis of PgWRKYs established their orthologous and paralogous relationship among the WRKY gene family of Arabidopsis thaliana, Oryza sativa and Setaria italica. Gene ontology (GO) annotation functionally categorized these PgWRKYs under cellular components, molecular functions and biological processes. Further, the differential expression pattern of PgWRKYs was noticed in different tissues (leaf, stem, root) and under both drought and salt stress conditions. The expression pattern of PgWRKY33, PgWRKY62 and PgWRKY65 indicates their probable involvement in both dehydration and salinity stress responses in pearl millet.

Conclusion

Functional characterization of identified PgWRKYs can be useful in delineating their role behind the natural stress tolerance of pearl millet against harsh environmental conditions. Further, these PgWRKYs can be employed in genome editing for millet crop improvement.

Genome-Wide Identification and Expression Profile Analysis of WRKY Family Genes in the Autopolyploid Saccharum spontaneum

WRKY is one of the largest transcription factor families in plants and plays important roles in the regulation of developmental and physiological processes. To date, the WRKY gene family has not been identified in Saccharum species because of its complex polyploid genome. In this study, a total of 294 sequences for 154 SsWRKY genes were identified in the polyploid Saccharum spontaneum genome and then named on the basis of their chromosome locations, including 13 (8.4%) genes with four alleles, 29 (18.8%) genes with three alleles and 41 (26.6%) genes with two alleles. Among them, 73.8% and 16.0% of the SsWRKY genes originated from segmental duplications and tandem duplications, respectively. The WRKY members exhibited conserved gene structures and amino acid sequences among the allelic haplotypes, which were accompanied by variations in intron sizes. Phylogenetic and collinearity analyses revealed that 27 SsWRKYs originated after the split of sorghum and Saccharum, resulting in a significantly higher number of WRKYs in sugarcane than in the proximal diploid species sorghum. The analysis of RNA-seq data revealed that SsWRKYs' expression profiles in 46 different samples including different developmental stages revealed distinct temporal and spatial patterns with 52 genes expressed in all tissues, four genes not expressed in any tissues and 21 SsWRKY genes likely to be involved in photosynthesis. The comprehensive analysis of SsWRKYs' expression will provide an important and valuable foundation for further investigation of the regulatory mechanisms of WRKYs in physiological roles in sugarcane S. spontaneum.

Genome-Wide Identification and Characterization of the NAC Transcription Factor Family in Musa Acuminata and Expression Analysis during Fruit Ripening

Banana (Musa acuminata, AAA group) is a representative climacteric fruit with essential nutrients and pleasant flavors. Control of its ripening determines both the fruit quality and the shelf life. NAC (NAM, ATAF, CUC2) proteins, as one of the largest superfamilies of transcription factors, play crucial roles in various functions, especially developmental processes. Thus, it is important to conduct a comprehensive identification and characterization of the NAC transcription factor family at the genomic level in M. acuminata. In this article, a total of 181 banana NAC genes were identified. Phylogenetic analysis indicated that NAC genes in M. acuminata, Arabidopsis, and rice were clustered into 18 groups (S1–S18), and MCScanX analysis disclosed that the evolution of MaNAC genes was promoted by segmental duplication events. Expression patterns of NAC genes during banana fruit ripening induced by ethylene were investigated using RNA-Seq data, and 10 MaNAC genes were identified as related to fruit ripening. A subcellular localization assay of selected MaNACs revealed that they were all localized to the nucleus. These results lay a good foundation for the investigation of NAC genes in banana toward the biological functions and evolution.

Decrypting tubby-like protein gene family of multiple functions in starch root crop cassava

Tubby-like proteins (TLPs) are ubiquitous in eukaryotes and function in abiotic stress tolerance of some plants. Cassava (Manihot esculenta Crantz) is a high-yield starch root crop and has a high tolerance to poor soil conditions and abiotic stress. However, little is known about TLP gene characteristics and their expression in cassava. We identified cassava TLP genes, MeTLPs, and further analysed structure, duplication, chromosome localization and collinearity, cis-acting elements in the promoter regions and expression patterns of MeTLPs, and three-dimensional structure of the encoded proteins MeTLPs. In conclusion, there is a MeTLP family containing 13 members, which are grouped into A and C subfamilies. There are 11 pairs of MeTLPs that show the duplication which took place between 10.11 and 126.69 million years ago. Two MeTLPs 6 and 9 likely originate from one gene in an ancestral species, may be common ancestors for other MeTLPs and would most likely not be eligible for ubiquitin-related protein degradation because their corresponding proteins (MeTLPs 6 and 9) have no the F-box domain in the N-terminus. MeTLPs feature differences in the number from TLPs in wheat, apple, Arabidopsis, poplar and maize, and are highlighted by segmental duplication but more importantly by the chromosomal collinearity with potato StTLPs. MeTLPs are at least related to abiotic stress tolerance in cassava. However, the subtle differences in function among MeTLPs are predictable partly because of their differential expression profiles, which are coupled with various cis‑acting elements existing in the promoter regions depending on genes.

Genome-Wide Analysis of the YABBY Transcription Factor Family in Pineapple and Functional Identification of AcYABBY4 Involvement in Salt Stress

The plant-specific transcription factor gene family, YABBY, belongs to the subfamily of zinc finger protein superfamily and plays an essential regulatory role in lateral organ development. In this study, nine YABBY genes were identified in the pineapple genome. Seven of them were located on seven different chromosomes and the remaining two were located on scaffold 1235. Through protein structure prediction and protein multiple sequence alignment, we found that AcYABBY3, AcYABBY5 and AcYABBY7 lack a C2 structure in their N-terminal C2C2 zinc finger protein structure. Analysis of the cis-acting element indicated that all the seven pineapple YABBY genes contain multiple MYB and MYC elements. Further, the expression patterns analysis using the RNA-seq data of different pineapple tissues indicated that different AcYABBYs are preferentially expressed in various tissues. RT-qPCR showed that the expression of AcYABBY2, AcYABBY3, AcYABBY6 and AcYABBY7 were highly sensitive to abiotic stresses. Subcellular localization in pineapple protoplasts, tobacco leaves and Arabidopsis roots showed that all the seven pineapple YABBY proteins were nucleus localized. Overexpression of AcYABBY4 in Arabidopsis resulted in short root under NaCl treatment, indicating a negative regulatory role of AcYABBY4 in plant resistance to salt stress. This study provides valuable information for the classification of pineapple AcYABBY genes and established a basis for further research on the functions of AcYABBY proteins in plant development and environmental stress response.

Comprehensive analysis of multiprotein bridging factor 1 family genes and SlMBF1c negatively regulate the resistance to Botrytis cinerea in tomato

BACKGROUND

Multiprotein bridging factor 1 s (MBF1s) are members of the transcriptional co-activator family that have involved in plant growth, development and stress responses. However, little is known about the Solanum lycopersicum MBF1 (SlMBF1) gene family.

RESULTS

In total, five SlMBF1 genes were identified based on the tomato reference genome, and these genes were mapped to five chromosomes. All of the SlMBF1 proteins were highly conserved, with a typical MBF1 domain and helix-turn-helix_3 domain. In addition, the promoter regions of the SlMBF1 genes have various stress and hormone responsive cis-regulatory elements. Encouragingly, the SlMBF1 genes were expressed with different expression profiles in different tissues and responded to various stress and hormone treatments. The biological function of SlMBF1c was further identified through its overexpression in tomato, and the transgenic tomato lines showed increased susceptibility to Botrytis cinerea (B. cinerea). Additionally, the expression patterns of salicylic acid (SA)-, jasmonic acid (JA)- and ethylene (ET)- mediated defense related genes were altered in the transgenic plants.

CONCLUSIONS

Our comprehensive analysis provides valuable information for clarifying the evolutionary relationship of the SlMBF1 members and their expression patterns in different tissues and under different stresses. The overexpression of SlMBF1c decreased the resistance of tomato to B. cinerea through enhancing the gene expression of the SA-mediated signaling pathway and depressing JA/ET-mediated signaling pathways. These results will facilitate future functional studies of the transcriptional co-activator family.

An Efficient Agrobacterium Mediated Transformation of Pineapple with GFP-Tagged Protein Allows Easy, Non-Destructive Screening of Transgenic Pineapple Plants

Quite a few studies have been conducted to improve the Agrobacterium-mediated transformation of pineapple, which is the second most important commercial tropical fruit crop worldwide. However, pineapple transformation remains challenging, due to technical difficulties, the lengthy regeneration process, and a high labor requirement. There have not been any studies specifically addressing the introduction of GFP-tagged genes into pineapples through Agrobacterium-mediated transformation, which would enable easy, non-destructive expression detection. It would also allow expression localization at the organelle level, which is not possible with GUS a reporter gene that encodes β-glucuronidase or a herbicide resistance reporter gene. Here, we report a method for the introduction of GFP-tagged genes into pineapples through Agrobacterium-mediated transformation. We used embryonic calli for transformation, and plants were regenerated through somatic embryogenesis. In this study, we optimized the incubation time for Agrobacterium infection, the co-cultivation time, the hygromycin concentration for selection, and the callus growth conditions after selection. Our strategy reduced the time required to obtain transgenic plants from 7.6 months to 6.1 months. The expression of GFP-tagged AcWRKY28 was observed in the nuclei of transgenic pineapple root cells. This method allows easy, non-destructive expression detection of transgenic constructs at the organelle level. These findings on pineapple transformation will help accelerate pineapple molecular breeding efforts to introduce new desirable traits.

Genome-Wide Identification and Characterization of Cucumber BPC Transcription Factors and Their Responses to Abiotic Stresses and Exogenous Phytohormones

BASIC PENTACYSTEINE (BPC) is a small transcription factor family that functions in diverse growth and development processes in plants. However, the roles of BPCs in plants, especially cucumber (Cucumis sativus L.), in response to abiotic stress and exogenous phytohormones are still unclear. Here, we identified four BPC genes in the cucumber genome, and classified them into two groups according to phylogenetic analysis. We also investigated the gene structures and detected five conserved motifs in these CsBPCs. Tissue expression pattern analysis revealed that the four CsBPCs were expressed ubiquitously in both vegetative and reproductive organs. Additionally, the transcriptional levels of the four CsBPCs were induced by various abiotic stress and hormone treatments. Overexpression of CsBPC2 in tobacco (Nicotiana tabacum) inhibited seed germination under saline, polyethylene glycol, and abscisic acid (ABA) conditions. The results suggest that the CsBPC genes may play crucial roles in cucumber growth and development, as well as responses to abiotic stresses and plant hormones. CsBPC2 overexpression in tobacco negatively affected seed germination under hyperosmotic conditions. Additionally, CsBPC2 functioned in ABA-inhibited seed germination and hypersensitivity to ABA-mediated responses. Our results provide fundamental information for further research on the biological functions of BPCs in development and abiotic stress responses in cucumber and other plant species.

Comprehensive Analysis and Functional Studies of WRKY Transcription Factors in Nelumbo nucifera

The WRKY family is one of the largest transcription factor (TF) families in plants and plays central roles in modulating plant stress responses and developmental processes, as well as secondary metabolic regulations. Lotus (Nelumbo nucifera) is an aquatic crop that has significant food, ornamental and pharmacological values. Here, we performed an overview analysis of WRKY TF family members in lotus, and studied their functions in environmental adaptation and regulation of lotus benzylisoquinoline alkaloid (BIA) biosynthesis. A total of 65 WRKY genes were identified in the lotus genome and they were well clustered in a similar pattern with their Arabidopsis homologs in seven groups (designated I, IIa-IIe, and III), although no lotus WRKY was clustered in the group IIIa. Most lotus WRKYs were functionally paired, which was attributed to the recently occurred whole genome duplication in lotus. In addition, lotus WRKYs were regulated dramatically by salicilic acid (SA), jasmonic acid (JA), and submergence treatments, and two lotus WRKYs, NnWRKY40a and NnWRKY40b, were significantly induced by JA and promoted lotus BIA biosynthesis through activating BIA biosynthetic genes. The investigation of WRKY TFs for this basal eudicot reveals new insights into the evolution of the WRKY family, and provides fundamental information for their functional studies and lotus breeding.

Genome-Wide Identification of WRKY Transcription Factors in the Asteranae

The WRKY transcription factors family, which participates in many physiological processes in plants, constitutes one of the largest transcription factor families. The Asterales and the Apiales are two orders of flowering plants in the superorder Asteranae. Among the members of the Asterales, globe artichoke (Cynara cardunculus var. scolymus L.), sunflower (Helianthus annuus L.), and lettuce (Lactuca sativa L.) are important economic crops worldwide. Within the Apiales, ginseng (Panax ginseng C. A. Meyer) and Panax notoginseng (Burk.) F.H. Chen are important medicinal plants, while carrot (Daucus carota subsp. carota L.) has significant economic value. Research involving genome-wide identification of WRKY transcription factors in the Asterales and the Apiales has been limited. In this study, 490 WRKY genes, 244 from three species of the Apiales and 246 from three species of the Asterales, were identified and categorized into three groups. Within each group, WRKY motif characteristics and gene structures were similar. WRKY gene promoter sequences contained light responsive elements, core regulatory elements, and 12 abiotic stress cis-acting elements. WRKY genes were evenly distributed on each chromosome. Evidence of segmental and tandem duplication events was found in all six species in the Asterales and the Apiales, with segmental duplication inferred to play a major role in WRKY gene evolution. Among the six species, we uncovered 54 syntenic gene pairs between globe artichoke and lettuce. The six species are thus relatively closely related, consistent with their traditional taxonomic placement in the Asterales. This study, based on traditional species classifications, was the first to identify WRKY transcription factors in six species from the Asteranae. Our results lay a foundation for further understanding of the role of WRKY transcription factors in species evolution and functional differentiation.

MYB Gene Family in Potato (Solanum tuberosum L.): Genome-Wide Identification of Hormone-Responsive Reveals Their Potential Functions in Growth and Development

As an important nongrain crop, the growth and yield of potato (Solanum tuberosum L.) is often affected by an unfavorable external environment in the process of cultivation. The MYB family is one of the largest and most important gene families, participating in the regulation of plant growth and development and response to abiotic stresses. Several MYB genes in potato that regulate anthocyanin synthesis and participate in abiotic stress responses have been identified. To identify all Solanum tuberosum L. MYB (StMYB) genes involved in hormone or stress responses to potentially regulate potato growth and development, we identified the MYB gene family at the genome-wide level. In this work, 158 StMYB genes were found in the potato genome. According to the amino acid sequence of the MYB domain and gene structure, the StMYB genes were divided into R2R3-MYB and R1R2R3-MYB families, and the R2R3-MYB family was divided into 20 subgroups (SGs). The expression of 21 StMYB genes from different SGs in roots, stems, leaves, flowers, shoots, stolons, young tubers, and mature tubers was determined by quantitative real-time polymerase chain reaction (qRT-PCR). The expression patterns of StMYB genes in potatoes treated with abscisic acid (ABA), indole-3-acetic acid (IAA), gibberellin acid 3 (GA3), NaCl, mannitol, and heat were also measured. We have identified several potential candidate genes that regulate the synthesis of potato flavonoids or participate in hormone or stress responses. This work provides a comprehensive understanding of the MYB family in potato and will lay a foundation for the future investigation of the potential functions of StMYB genes in the growth and development of potato.

Identification of SWI2/SNF2-Related 1 Chromatin Remodeling Complex (SWR1-C) Subunits in Pineapple and the Role of Pineapple SWR1 COMPLEX 6 (AcSWC6) in Biotic and Abiotic Stress Response

Chromatin remodeling complex orchestrates numerous aspects of growth and development in eukaryotes. SWI2/SNF2-Related 1 chromatin remodeling complex (SWR1-C) is a member of the SWI/SNF ATPase-containing chromatin remodeling complex responsible for the exchange of H2A for H2A.Z. In plants, SWR1-C plays a crucial role by transcriptionally regulating numerous biological and developmental processes. However, SWR1-C activity remains obscure in pineapple. Here, we aim to identify the SWR1-C subunits in pineapple. By genome-wide identification, we found a total of 11 SWR1-C subunits in the pineapple. The identified SWR1-C subunits were named and classified based on the sequence similarity and phylogenetic analysis. RNA-Seq analysis showed that pineapple SWR1-C subunits are expressed differentially in different organs and at different stages. Additionally, the qRT-PCR of pineapple SWR1-C subunits during abiotic stress exposure showed significant changes in their expression. We further investigated the functions of pineapple SWR1 COMPLEX 6 (AcSWC6) by ectopically expressing it in Arabidopsis. Interestingly, transgenic plants ectopically expressing AcSWC6 showed susceptibility to fungal infection and enhanced resistance to salt and osmotic stress, revealing its involvement in biotic and abiotic stress. Moreover, the complementation of mutant Arabidopsisswc6 by pineapple SWC6 suggested the conserved function of SWC6 in plants.

Genome-wide exploration and characterization of miR172/euAP2 genes in Brassica napus L. for likely role in flower organ development

BACKGROUND

APETALA2-like genes encode plant-specific transcription factors, some of which possess one microRNA172 (miR172) binding site. The miR172 and its target euAP2 genes are involved in the process of phase transformation and flower organ development in many plants. However, the roles of miR172 and its target AP2 genes remain largely unknown in Brassica napus (B. napus).

RESULTS

In this study, 19 euAP2 and four miR172 genes were identified in the B. napus genome. A sequence analysis suggested that 17 euAP2 genes were targeted by Bna-miR172 in the 3' coding region. EuAP2s were classified into five major groups in B.napus. This classification was consistent with the exon-intron structure and motif organization. An analysis of the nonsynonymous and synonymous substitution rates revealed that the euAP2 genes had gone through purifying selection. Whole genome duplication (WGD) or segmental duplication events played a major role in the expansion of the euAP2 gene family. A cis-regulatory element (CRE) analysis suggested that the euAP2s were involved in the response to light, hormones, stress, and developmental processes including circadian control, endosperm and meristem expression. Expression analysis of the miR172-targeted euAP2s in nine different tissues showed diverse spatiotemporal expression patterns. Most euAP2 genes were highly expressed in the floral organs, suggesting their specific functions in flower development. BnaAP2-1, BnaAP2-5 and BnaTOE1-2 had higher expression levels in late-flowering material than early-flowering material based on RNA-seq and qRT-PCR, indicating that they may act as floral suppressors.

CONCLUSIONS

Overall, analyses of the evolution, structure, tissue specificity and expression of the euAP2 genes were peformed in B.napus. Based on the RNA-seq and experimental data, euAP2 may be involved in flower development. Three euAP2 genes (BnaAP2-1, BnaAP2-5 and BnaTOE1-2) might be regarded as floral suppressors. The results of this study provide insights for further functional characterization of the miR172 /euAP2 module in B.napus.

Systematic Analysis of MYB Family Genes in Potato and Their Multiple Roles in Development and Stress Responses

The MYB proteins represent a large family of transcription factors and play important roles in development, senescence, and stress responses in plants. In the current study, 233 MYB transcription factor-encoding genes were identified and analyzed in the potato genome, including 119 R1-MYB, 112 R2R3-MYB, and two R1R2R3-MYB members. R2R3-MYB is the most abundant MYB subclass and potato R2R3-MYB members together with their Arabidopsis homologs were divided into 35 well-supported subgroups as the result of phylogenetic analyses. Analyses on gene structure and protein motif revealed that members from the same subgroup shared similar exon/intron and motif organization, further supporting the results of phylogenetic analyses. Evolution of the potato MYB family was studied via syntenic analysis. Forty-one pairs of StMYB genes were predicted to have arisen from tandem or segmental duplication events, which played important roles in the expansion of the StMYB family. Expression profiling revealed that the StMYB genes were expressed in various tissues and several StMYB genes were identified to be induced by different stress conditions. Notably, StMYB030 was found to act as the homolog of AtMYB44 and was significantly up-regulated by salt and drought stress treatments. Furthermore, overexpression of StMYB030 in Arabidopsis enhanced salt stress tolerance of transgenic plants. The results from this study provided information for further functional analysis and for crop improvements through genetic manipulation of these StMYB genes.

Molecular identification of the phosphate transporter family 1 (PHT1) genes and their expression profiles in response to phosphorus deprivation and other abiotic stresses in Brassica napus

Phosphate (Pi) transporters play critical roles in Pi acquisition and homeostasis. However, little is known about these transporters in oilseed rape. Therefore, the aim of the present study was to characterize the members of the PHT1 gene family in allotetraploid Brassica napus and to analyze their expression profiles in response to environmental stresses. In total, 49 PHT1 family members were identified in B. napus, including 27 genes in the A subgenome and 22 in the C subgenome. Most of the PHT1 proteins were predicted to localize to the plasma membrane. Phylogenetic analysis suggested that the members of the PHT1 gene family can be divided into seven clades, with the introns/exons and protein motifs conserved in each clade. Collinearity analysis revealed that most of the BnaPHT1 genes shared syntenic relationships with PHT1 members in Arabidopsis thaliana, B. rapa, and B. oleracea, and that whole-genome duplication (polyploidy) played a major driving force for BnaPHT1 evolution in addition to segmental duplication. Transcript abundance analysis showed that a broad range of expression patterns of individual BnaPHT1 genes occurred in response to phosphorus (P) deficiency. In addition, the expression levels of BnaPHT1 genes can be regulated by different nutrient stresses, including nitrogen (N), potassium (K), sulfur (S) and iron (Fe) stresses. Moveover, salt and drought stresses can regulate the transcript abundances of BnaPHT1s, as well as phytohormones including auxin and cytokinin. Gene coexpression analysis based on the RNA-seq data implied that BnaPHT1s might cooperate with each other as well as with other genes to regulate nutrient homeostasis in B. napus. Further analysis of the promoters revealed that GT-1, DRE and P1BS elements are widely distributed within the promoter regions of BnaPHT1 genes. Our results indicate that BnaPHT1s might be involved in cross-talk for sensing the external status of P, N, K, S and Fe, as well as salt and drought stresses. Moreover, these processes might be mediated by phytohormones. Our findings provide the first step in the complex genetic dissection of the Pi transport system in plants and implicate multiple transcriptional regulation, which probably refers to new roles of PHT1 genes in B. napus.

Genome-Wide Identification and Expression Profiling of CBL-CIPK Gene Family in Pineapple (Ananas comosus) and the Role of AcCBL1 in Abiotic and Biotic Stress Response

Ca²⁺ serves as a ubiquitous second messenger regulating several aspects of plant growth and development. A group of unique calcium sensor proteins, calcineurin B-like (CBL), interact with CBL-interacting protein kinases (CIPKs) to decode the Ca²⁺ signature inside the cell. Although CBL-CIPK signaling toolkit has been shown to play significant roles in the responses to numerous stresses in different plants, the information about pineapple CBL-CIPK remains obscure. In the present study, a total of eight AcCBL and 21 AcCIPK genes were identified genome-wide in pineapple. The identified genes were renamed on the basis of gene ID in ascending order and phylogenetic analysis divided into five groups. Transcriptomic data analysis showed that AcCBL and AcCIPK genes were expressed differentially in different tissues. Further, the expression analysis of AcCBL1 in different tissues showed significant changes under various abiotic stimuli. Additionally, the ectopic expression of AcCBL1 in Arabidopsis resulted in enhanced tolerance to salinity, osmotic, and fungal stress. The present study revealed the crucial contribution of the CBL-CIPK gene in various biological and physiological processes in pineapple.

Genome-wide identification of GRF transcription factors in soybean and expression analysis of GmGRF family under shade stress

BACKGROUND

The Growth-regulating factor (GRF) family encodes plant-specific transcription factors which contain two conserved domains, QLQ and WRC. Members of this family play vital roles in plant development and stress response processes. Although GRFs have been identified in various plant species, we still know little about the GRF family in soybean (Glycine max).

RESULTS

In the present study, 22 GmGRFs distributed on 14 chromosomes and one scaffold were identified by searching soybean genome database and were clustered into five subgroups according to their phylogenetic relationships. GmGRFs belonging to the same subgroup shared a similar motif composition and gene structure. Synteny analysis revealed that large-scale duplications played key roles in the expansion of the GmGRF family. Tissue-specific expression data showed that GmGRFs were strongly expressed in growing tissues, including the shoot apical meristems, developing seeds and flowers, indicating that GmGRFs play critical roles in plant growth and development. On the basis of expression analysis of GmGRFs under shade conditions, we found that all GmGRFs responded to shade stress. Most GmGRFs were down-regulated in soybean leaves after shade treatment.

CONCLUSIONS

Taken together, this research systematically analyzed the characterization of the GmGRF family and its primary roles in soybean development and shade stress response. Further studies of the function of the GmGRFs in the growth, development and stress tolerance of soybean, especially under shade stress, will be valuable.

Genome-wide characterization of the WRKY gene family in cultivated strawberry (Fragaria × ananassa Duch.) and the importance of several group III members in continuous cropping

WRKY transcription factors play important roles in many plant developmental processes and adaptation to the environment. However, little knowledge is available about the WRKY gene family in cultivated strawberry (Fragaria × ananassa Duch.), an important soft fruit worldwide. In this study, a total of 47 WRKY gene members were identified and renamed on the basis of their order on the chromosomes. According to their evolutionary events and conserved structure, the 47 FaWRKYs were divided into three major groups with several subgroups. A cis-element analysis showed that all FaWRKYs possessed at least one stress response-related cis-element. Comprehensive analysis, including phylogenetic analysis and expression profiling, based on real-time qPCR analysis in root, stem, leaf and fruit was performed on group III FaWRKY genes. The phylogenetic tree of the WRKY III genes in cultivated strawberry, wild Strawberry, Arabidopsis, tomato, and rice was divided into five clades. Additionally, the expression profiles of the FaWRKY genes in response to continuous cropping were further investigated based on RNA-seq data. FaWRKY25, FaWRKY32, and FaWRKY45, which are group III FaWRKY genes, were upregulated after continuous cropping. The level of reactive oxygen species (ROS) and the expression levels of PR1 and peroxidase were higher in continuous cropping (CC) than in non-continuous cropping (NCC). The results indicated that group III FaWRKYs might play an important role in continuous cropping. These results provide a foundation for genetic improvements for continuous cropping tolerance in cultivated strawberry.

Ectopic expression of the Aechmea fasciata APETALA2 gene AfAP2-2 reduces seed size and delays flowering in Arabidopsis

The Bromeliaceae family, which is distributed pantropically, is one of the most morphologically diverse families. Except for the edible pineapple (Ananas comosus), the vast majority of bromeliads cultivated worldwide are appreciated mainly for their ornamental value. As subtropical and tropical flowering plants, these bromeliads, among with Aechmea fasciata, have significant economic importance. However, the molecular mechanism of flowering in bromeliads remains unrevealed. In this study, an APETALA2 (AP2) homologue, AfAP2-2, which belongs to the AP2/ethylene response element binding protein (AP2/EREBP) transcription factor superfamily, was identified in A. fasciata. AfAP2-2 contains two conserved AP2 domains and is a nuclear-localized transactivator. The expression level of AfAP2-2 was predominantly higher in vegetative organs of the reproductive phase than in those of the vegetative phase. Ectopic expression of AfAP2-2 in Arabidopsis specifically delayed flowering in short-day (SD) conditions. Furthermore, the size and weight of seeds of AfAP2-2-overexpressing Arabidopsis plants were significantly reduced compared to those of the wild type (WT). Our findings suggest that AfAP2-2 might be a negative regulator of flowering and seed size and weight. These results may help facilitate the molecular breeding of bromeliads.

Revisiting the Role of Transcription Factors in Coordinating the Defense Response Against Citrus Bark Cracking Viroid Infection in Commercial Hop (Humulus Lupulus L.)

Transcription factors (TFs) play a major role in controlling gene expression by intricately regulating diverse biological processes such as growth and development, the response to external stimuli and the activation of defense responses. The systematic identification and classification of TF genes are essential to gain insight into their evolutionary history, biological roles, and regulatory networks. In this study, we performed a global mining and characterization of hop TFs and their involvement in Citrus bark cracking viroid CBCVd infection by employing a digital gene expression analysis. Our systematic analysis resulted in the identification of a total of 3,818 putative hop TFs that were classified into 99 families based on their conserved domains. A phylogenetic analysis classified the hop TFs into several subgroups based on a phylogenetic comparison with reference TF proteins from Arabidopsis thaliana providing glimpses of their evolutionary history. Members of the same subfamily and subgroup shared conserved motif compositions. The putative functions of the CBCVd-responsive hop TFs were predicted using their orthologous counterparts in A. thaliana. The analysis of the expression profiling of the CBCVd-responsive hop TFs revealed a massive differential modulation, and the expression of the selected TFs was validated using qRT-PCR. Together, the comprehensive integrated analysis in this study provides better insights into the TF regulatory networks associated with CBCVd infections in the hop, and also offers candidate TF genes for improving the resistance in hop against viroids.

Genome-Wide Characterization, Evolution, and Expression Profiling of VQ Gene Family in Response to Phytohormone Treatments and Abiotic Stress in Eucalyptus grandis

VQ genes play important roles in plant development, growth, and stress responses. However, little information regarding the functions of VQ genes is available for Eucalyptus grandis. In our study, genome-wide characterization and identification of VQ genes were performed in E. grandis. Results showed that 27 VQ genes, which divided into seven sub-families (I-VII), were found, and all but two VQ genes showed no intron by gene structure and conserved motif analysis. To further identify the function of EgrVQ proteins, gene expression analyses were also developed under hormone treatments (brassinosteroids, methyl jasmonate, salicylic acid, and abscisic acid) and abiotic conditions (salt stress, cold 4 °C, and heat 42 °C). The results of a quantitative real-time PCR analysis indicated that the EgrVQs were variously expressed under different hormone treatments and abiotic stressors. Our study provides a comprehensive overview of VQ genes in E. grandis, which will be beneficial in the molecular breeding of E. grandis to promote its resistance to abiotic stressors; the results also provide a basis from which to conduct further investigation into the functions of VQ genes in E. grandis.

Transcriptome Analysis Reveals Key Cold-Stress-Responsive Genes in Winter Rapeseed (Brassica rapa L.)

Low ambient air temperature limits the growth and selection of crops in cold regions, and cold tolerance is a survival strategy for overwintering plants in cold winters. Studies of differences in transcriptional levels of winter rapeseed (Brassica rapa L.) under cold stress can improve our understanding of transcript-mediated cold stress responses. In this study, two winter rapeseed varieties, Longyou-7 (cold-tolerant) and Lenox (cold-sensitive), were used to reveal morphological, physiological, and transcriptome levels after 24 h of cold stress, and 24 h at room temperature, to identify the mechanism of tolerance to cold stress. Compared to Lenox, Longyou-7 has a shorter growth period and greater belowground mass, and exhibits stronger physiological activity after cold stress. Subsequently, more complete genomic annotation was obtained by sequencing. A total of 10,251 and 10,972 differentially expressed genes (DEG) were identified in Longyou-7 and Lenox, respectively. Six terms closely related to cold stress were found by the Gene Ontology (GO) function annotation. Some of these terms had greater upregulated expression in Longyou-7, and the expression of these genes was verified by qRT-PCR. Most of these DEGs are involved in phenylpropanoid biosynthesis, plant hormone signal transduction, ribosome biogenesis, MAPK signaling pathway, basal transcription factors, and photosynthesis. Analysis of the genes involved in the peroxisome pathway revealed that Longyou-7 and Lenox may have different metabolic patterns. Some transcription factors may play an important role in winter rapeseed tolerance to cold stress, and Longyou-7 is slightly slower than Lenox. Our results provide a transcriptome database and candidate genes for further study of winter rapeseed cold stress.

Evolution and Identification of the WRKY Gene Family in Quinoa (Chenopodium quinoa)

The WRKY gene family plays a unique role in plant stress tolerance. Quinoa is a cultivated crop worldwide that is known for its high stress tolerance. The WRKY gene family in quinoa has not yet been studied. Using a genome-wide search method, we identified 1226 WRKY genes in 15 plant species, seven animal species, and seven fungi species. WRKY proteins were not found in animal species and five fungi species, but were, however, widespread in land plants. A total of 92 CqWRKY genes were identified in quinoa. Based on the phylogenetic analysis, these CqWRKY genes were classified into three groups. The CqWRKY proteins have a highly conserved heptapeptide WRKYGQK with 15 conserved elements. Furthermore, a total of 25 CqWRKY genes were involved in the co-expression pathway of organ development and osmotic stress. The expression level of more than half of these CqWRKY genes showed significant variation under salt or drought stress. This study reports, for the first time, the findings of the CqWRKY gene family in quinoa at the genome-wide level. This information will be beneficial for our understanding of the molecular mechanisms of stress tolerance in crops, such as quinoa.

Combining bioinformatics and conventional PCR optimization strategy for one-time design of high-specificity primers for WRKY gene family using unigene database

Gene families, like the conserved transcription factor families, evolve through gene duplications and share moderate similarity between member genes. Lack of genomic data makes it difficult to design high-specificity primers to the target genes. Furthermore, many primers under-perform in highly sensitive assays like quantitative PCR due to issues of thermodynamic nature, thereby increasing the cost and time for analysis. A methodology involving intra-species and inter-generic bioinformatic sequence comparison combined with thermodynamic estimation of primer performance was used for one-time design of gene specific primers for different WRKYs, Mitogen Activated Protein-kinases and N-methyltransferases of Coffea canephora without the aid of genome sequence resources. Out of a total 37 primer sets including 31 pairs of primers for WRKY from 34 mined WRKY Unigenes/ESTs and six pairs for genes coding for MAP kinases and NBS-LRR proteins, 32 sets exhibited high specificity of amplification upon genome analysis as well as in the high-resolution melt analysis. Furthermore, PCR optimization strategies-both in silico and experimental-indicated a superior performance of the primer sets for different applications like quantitative PCR and rapid amplification of cDNA ends. Only one set of primer resulted in mis-priming upon confirmation by DNA sequencing of the cloned amplicons. The intra-species differences and inter-generic similarities ensure high specificity of primers in all cases studied. The procedure allowed design of primers for the use in different downstream applications with high performance, specificity, yield and ease-of-use.

GmWRKY16 Enhances Drought and Salt Tolerance Through an ABA-Mediated Pathway in Arabidopsis thaliana

The WRKY transcription factors (TFs) are one of the largest families of TFs in plants and play multiple roles in plant development and stress response. In the present study, GmWRKY16 encoding a WRKY transcription factor in soybean was functionally characterized in Arabidopsis. GmWRKY16 is a nuclear protein that contains a highly conserved WRKY domain and a C2H2 zinc-finger structure, and has the characteristics of transcriptional activation ability, presenting a constitutive expression pattern with relative expression levels of over fourfold in the old leaves, flowers, seeds and roots of soybean. The results of quantitative real time polymerase chain reaction (qRT-PCR) showed that GmWRKY16 could be induced by salt, alkali, ABA, drought and PEG-6000. As compared with the control, overexpression of GmWRKY16 in Arabidopsis increased the seed germination rate and root growth of seedlings in transgenic lines under higher concentrations of mannitol, NaCl and ABA. In the meantime, GmWRKY16 transgenic lines showed over 75% survival rate after rehydration and enhanced Arabidopsis tolerance to salt and drought with higher proline and lower MDA accumulation, less water loss of the detached leaves, and accumulated more endogenous ABA than the control under stress conditions. Further studies showed that AtWRKY8, KIN1, and RD29A were induced in GmWRKY16 transgenic plants under NaCl treatment. The expressions of the ABA biosynthesis gene (NCED3), signaling genes (ABI1, ABI2, ABI4, and ABI5), responsive genes (RD29A, COR15A, COR15B, and RD22) and stress-related marker genes (KIN1, LEA14, LEA76, and CER3) were regulated in transgenic lines under drought stress. In summary, these results suggest that GmWRKY16 as a WRKY TF may promote tolerance to drought and salt stresses through an ABA-mediated pathway.

The WRKY Transcription Factor GmWRKY12 Confers Drought and Salt Tolerance in Soybean

WRKYs are important regulators in plant development and stress responses. However, knowledge of this superfamily in soybean is limited. In this study, we characterized the drought- and salt-induced gene GmWRKY12 based on RNA-Seq and qRT-PCR. GmWRKY12, which is 714 bp in length, encoded 237 amino acids and grouped into WRKY II. The promoter region of GmWRKY12 included ABER4, MYB, MYC, GT-1, W-box and DPBF cis-elements, which possibly participate in abscisic acid (ABA), drought and salt stress responses. GmWRKY12 was minimally expressed in different tissues under normal conditions but highly expressed under drought and salt treatments. As a nucleus protein, GmWRKY12 was responsive to drought, salt, ABA and salicylic acid (SA) stresses. Using a transgenic hairy root assay, we further characterized the roles of GmWRKY12 in abiotic stress tolerance. Compared with control (Williams 82), overexpression of GmWRKY12 enhanced drought and salt tolerance, increased proline (Pro) content and decreased malondialdehyde (MDA) content under drought and salt treatment in transgenic soybean seedlings. These results may provide a basis to understand the functions of GmWRKY12 in abiotic stress responses in soybean.

Comprehensive tissue-specific transcriptome profiling of pineapple (Ananas comosus) and building an eFP-browser for further study

Pineapple is one of the most economically important tropical or subtropical fruit trees. However, few studies focus on the development of its unique collective fruit. In this study, we generated a genome-wide developmental transcriptomic profile of 14 different tissues of the collective fruit of the pineapple covering each of the three major fruit developmental stages. In total, 273 tissue-specific and 1,051 constitutively expressed genes were detected. We also performed gene co-expression analysis and 18 gene modules were classified. Among these, we found three interesting gene modules; one was preferentially expressed in bracts and sepals and was likely involved in plant defense; one was highly expressed at the beginning of fruit expansion and faded afterward and was probably involved in endocytosis; Another gene module increased expression level with pineapple fruit development and was involved in terpenoid and polyketide metabolism. In addition, we built a pineapple electronic fluorescent pictograph (eFP) browser to facilitate exploration of gene expression during pineapple fruit development. With this tool, users can visualize expression data in this study in an intuitive way. Together, the transcriptome profile generated in this work and the corresponding eFP browser will facilitate further study of fruit development in pineapple.

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.

Identification of WRKY Gene Family from Dimocarpus longan and Its Expression Analysis during Flower Induction and Abiotic Stress Responses

Longan is an important fruit tree in the subtropical region of Southeast Asia and Australia. However, its blooming and its yield are susceptible to stresses such as droughts, high salinity, and high and low temperature. To date, the molecular mechanisms of abiotic stress tolerance and flower induction in longan have not been elucidated. WRKY transcription factors (TFs), which have been studied in various plant species, play important regulatory roles in plant growth, development, and responses to stresses. However, there is no report about WRKYs in longan. In this study, we identified 55 WRKY genes with the conserved WRKY domain and zinc finger motif in the longan genome. Based on the structural features of WRKY proteins and topology of the phylogenetic tree, the longan WRKY (DlWRKY) family was classified into three major groups (I–III) and five subgroups (IIa–IIe) in group II. Tissue expression analysis showed that 25 DlWRKYs were highly expressed in almost all organs, suggesting that these genes may be important for plant growth and organ development in longan. Comparative RNA-seq and qRT-PCR-based gene expression analysis revealed that 18 DlWRKY genes showed a specific expression during three stages of flower induction in “Sijimi” (“SJ”), which exhibited the “perpetual flowering” (PF) habit, indicating that these 18 DlWRKY genes may be involved in the flower induction and the genetic control of the perpetual flowering trait in longan. Furthermore, the RT-qPCR analysis illustrated the significant variation of 27, 18, 15, 17, 27, and 23 DlWRKY genes under SA (Salicylic acid), MeJA (Methyl Jasmonate), heat, cold, drought, or high salinity treatment, respectively, implicating that they might be stress- or hormone-responsive genes. In summary, we systematically and comprehensively analyzed the structure, evolution, and expression pattern of the DlWRKY genes. The results presented here increase our understanding of the WRKY family in fruit trees and provide a basis for the further elucidation of the biological function of DlWRKY genes in longan.