Genome-Wide Identification and Expression Profiling of KCS Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in passion fruit (Passiflora edulis) and their involvement in flower and fruit development

BACKGROUND

The calmodulin (CaM) and calmodulin-like (CML) proteins play regulatory roles in plant growth and development, responses to biotic and abiotic stresses, and other biological processes. As a popular fruit and ornamental crop, it is important to explore the regulatory mechanism of flower and fruit development of passion fruit.

RESULTS

In this study, 32 PeCaM/PeCML genes were identified from passion fruit genome and were divided into 9 groups based on phylogenetic analysis. The structural analysis, including conserved motifs, gene structure and homologous modeling, illustrates that the PeCaM/PeCML in the same subgroup have relative conserved structural features. Collinearity analysis suggested that the expansion of the CaM/CML gene family likely took place mainly by segmental duplication, and the whole genome replication events were closely related with the rapid expansion of the gene group. PeCaM/PeCMLs were potentially required for different floral tissues development. Significantly, PeCML26 had extremely high expression levels during ovule and fruit development compared with other PeCML genes, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. The co-presence of various cis-elements associated with growth and development, hormone responsiveness, and stress responsiveness in the promoter regions of these PeCaM/PeCMLs might contribute to their diverse regulatory roles. Furthermore, PeCaM/PeCMLs were also induced by various abiotic stresses. This work provides a comprehensive understanding of the CaM/CML gene family and valuable clues for future studies on the function and evolution of CaM/CML genes in passion fruit.

CONCLUSION

A total of 32 PeCaM/PeCML genes were divided into 9 groups. The PeCaM/PeCML genes showed differential expression patterns in floral tissues at different development stages. It is worth noting that PeCML26, which is highly homologous to AtCaM2, not only interacts with multiple BBR-BPC TFs, but also has high expression levels during ovule and fruit development, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. This research lays the foundation for future investigations and validation of the potential function of PeCaM/PeCML genes in the growth and development of passion fruit.

Class I TCP in fruit development: much more than growth

Fruit development can be viewed as the succession of three main steps consisting of the fruit initiation, growth and ripening. These processes are orchestrated by different factors, notably the successful fertilization of flowers, the environmental conditions and the hormones whose action is coordinated by a large variety of transcription factors. Among the different transcription factor families, TEOSINTE BRANCHED 1, CYCLOIDEA, PROLIFERATING CELL FACTOR (TCP) family has received little attention in the frame of fruit biology despite its large effects on several developmental processes and its action as modulator of different hormonal pathways. In this respect, the comprehension of TCP functions in fruit development remains an incomplete puzzle that needs to be assembled. Building on the abundance of genomic and transcriptomic data, this review aims at collecting available TCP expression data to allow their integration in the light of the different functional genetic studies reported so far. This reveals that several Class I TCP genes, already known for their involvement in the cell proliferation and growth, display significant expression levels in developing fruit, although clear evidence supporting their functional significance in this process remains scarce. The extensive expression data compiled in our study provide convincing elements that shed light on the specific involvement of Class I TCP genes in fruit ripening, once these reproductive organs acquire their mature size. They also emphasize their putative role in the control of specific biological processes such as fruit metabolism and hormonal dialogue.

Genome-Wide Comprehensive Identification and In Silico Characterization of Lectin Receptor-Like Kinase Gene Family in Barley (Hordeum vulgare L.)

Lectin receptor-like kinases (LecRLKs) are a significant subgroup of the receptor-like kinases (RLKs) protein family. They play crucial roles in plant growth, development, immune responses, signal transduction, and stress tolerance. However, the genome-wide identification and characterization of LecRLK genes and their regulatory elements have not been explored in a major cereal crop, barley (Hordeum vulgare L.). Therefore, in this study, integrated bioinformatics tools were used to identify and characterize the LecRLK gene family in barley. Based on the phylogenetic tree and domain organization, a total of 113 LecRLK genes were identified in the barley genome (referred to as HvlecRLK) corresponding to the LecRLK genes of Arabidopsis thaliana. These putative HvlecRLK genes were classified into three groups: 62 G-type LecRLKs, 1 C-type LecRLK, and 50 L-type LecRLKs. They were unevenly distributed across eight chromosomes, including one unknown chromosome, and were predominantly located in the plasma membrane (G-type HvlecRLK (96.8%), C-type HvlecRLK (100%), and L-type HvlecRLK (98%)). An analysis of motif composition and exon-intron configuration revealed remarkable homogeneity with the members of AtlecRLK. Notably, most of the HvlecRLKs (27 G-type, 43 L-type) have no intron, suggesting their rapid functionality. The Ka/Ks and syntenic analysis demonstrated that HvlecRLK gene pairs evolved through purifying selection and gene duplication was the major factor for the expansion of the HvlecRLK gene family. Exploration of gene ontology (GO) enrichment indicated that the identified HvlecRLK genes are associated with various cellular processes, metabolic pathways, defense mechanisms, kinase activity, catalytic activity, ion binding, and other essential pathways. The regulatory network analysis identified 29 transcription factor families (TFFs), with seven major TFFs including bZIP, C2H2, ERF, MIKC_MADS, MYB, NAC, and WRKY participating in the regulation of HvlecRLK gene functions. Most notably, eight TFFs were found to be linked to the promoter region of both L-type HvleckRLK64 and HvleckRLK86. The promoter cis-acting regulatory element (CARE) analysis of barley identified a total of 75 CARE motifs responsive to light responsiveness (LR), tissue-specific (TS), hormone responsiveness (HR), and stress responsiveness (SR). The maximum number of CAREs was identified in HvleckRLK11 (25 for LR), HvleckRLK69 (17 for TS), and HvleckRLK80 (12 for HR). Additionally, HvleckRLK14, HvleckRLK16, HvleckRLK33, HvleckRLK50, HvleckRLK52, HvleckRLK56, and HvleckRLK110 were predicted to exhibit higher responses in stress conditions. In addition, 46 putative miRNAs were predicted to target 81 HvlecRLK genes and HvlecRLK13 was the most targeted gene by 8 different miRNAs. Protein-protein interaction analysis demonstrated higher functional similarities of 63 HvlecRLKs with 7 Arabidopsis STRING proteins. Our overall findings provide valuable information on the LecRLK gene family which might pave the way to advanced research on the functional mechanism of the candidate genes as well as to develop new barley cultivars in breeding programs.

Genome-wide identification and expression analysis of the KCS gene family in soybean (Glycine max) reveal their potential roles in response to abiotic stress

Very long chain fatty acids (VLCFAs) are fatty acids with chain lengths of 20 or more carbon atoms, which are the building blocks of various lipids that regulate developmental processes and plant stress responses. 3-ketoacyl-CoA synthase encoded by the KCS gene is the key rate-limiting enzyme in VLCFA biosynthesis, but the KCS gene family in soybean (Glycine max) has not been adequately studied thus far. In this study, 31 KCS genes (namely GmKCS1 - GmKCS31) were identified in the soybean genome, which are unevenly distributed on 14 chromosomes. These GmKCS genes could be phylogenetically classified into seven groups. A total of 27 paralogous GmKCS gene pairs were identified with their Ka/Ks ratios indicating that they had undergone purifying selection during soybean genome expansion. Cis-acting element analysis revealed that GmKCS promoters contained multiple hormone- and stress-responsive elements, indicating that GmKCS gene expression levels may be regulated by various developmental and environmental stimuli. Expression profiles derived from RNA-seq data and qRT-PCR experiments indicated that GmKCS genes were diversely expressed in different organs/tissues, and many GmKCS genes were found to be differentially expressed in the leaves under cold, heat, salt, and drought stresses, suggesting their critical role in soybean resistance to abiotic stress. These results provide fundamental information about the soybean KCS genes and will aid in their further functional elucidation and exploitation.

Cotyledon peeling method for passion fruit protoplasts: a versatile cell system for transient gene expression in passion fruit (Passiflora edulis)

Passion fruit (Passiflora edulis) is a perennial evergreen vine that grows mainly in tropical and subtropical regions due to its nutritional, medicinal and ornamental values. However, the molecular biology study of passion fruit is extremely hindered by the lack of an easy and efficient method for transformation. The protoplast transformation system plays a vital role in plant regeneration, gene function analysis and genome editing. Here, we present a new method ('Cotyledon Peeling Method') for simple and efficient passion fruit protoplast isolation using cotyledon as the source tissue. A high yield (2.3 × 107 protoplasts per gram of fresh tissues) and viability (76%) of protoplasts were obtained upon incubation in the enzyme solution [1% (w/v) cellulase R10, 0.25% (w/v) macerozyme R10, 0.4 M mannitol, 10 mM CaCl2, 20 mM KCl, 20 mM MES and 0.1% (w/v) BSA, pH 5.7] for 2 hours. In addition, we achieved high transfection efficiency of 83% via the polyethylene glycol (PEG)-mediated transformation with a green fluorescent protein (GFP)-tagged plasmid upon optimization. The crucial factors affecting transformation efficiency were optimized as follows: 3 μg of plasmid DNA, 5 min transfection time, PEG concentration at 40% and protoplast density of 100 × 104 cells/ml. Furthermore, the established protoplast system was successfully applied for subcellular localization analysis of multiple fluorescent organelle markers and protein-protein interaction study. Taken together, we report a simple and efficient passion fruit protoplast isolation and transformation system, and demonstrate its usage in transient gene expression for the first time in passion fruit. The protoplast system would provide essential support for various passion fruit biology studies, including genome editing, gene function analysis and whole plant regeneration.

Comprehensive analysis of KCS gene family in pear reveals the involvement of PbrKCSs in cuticular wax and suberin synthesis and pear fruit skin formation

Cuticular wax, cutin and suberin polyesters covering the surface of some fleshy fruit are tightly associated with skin color and appearance. β-Ketoacyl-CoA synthase (KCS) is a rate-limiting enzyme participating in the synthesis of very-long-chain fatty acids (VLCFAs), the essential precursors of cuticular waxes and aliphatic monomers of suberin. However, information on the KCS gene family in pear genome and the specific members involved in pear fruit skin formation remain unclear. In the present study, we performed an investigation of the composition and amount of cuticular waxes, cutin and aliphatic suberin in skins of four sand pear varieties with distinct colors (russet, semi-russet, and green) and demonstrated that the metabolic shifts of cuticular waxes and suberin leading to the significant differences of sand pear skin color. A genome-wide identification of KCS genes from the pear genome was conducted and 35 KCS coding genes were characterized and analyzed. Expression profile analysis revealed that the KCS genes had diverse expression patterns among different pear skins and the transcript abundance of PbrKCS15, PbrKCS19, PbrKCS24, and PbrKCS28 were consistent with the accumulation of cuticular waxes and suberin in fruit skin respectively. Subcellular localization analysis demonstrated that PbrKCS15, PbrKCS19, PbrKCS24 and PbrKCS28 located on the endoplasmic reticulum (ER). Further, transient over-expression of PbrKCS15, PbrKCS19, and PbrKCS24 in pear fruit skins significantly increased cuticular wax accumulation, whereas PbrKCS28 notably induced suberin deposition. In conclusion, pear fruit skin color and appearance are controlled in a coordinated way by the deposition of the cuticular waxes and suberin. PbrKCS15, PbrKCS19, and PbrKCS24 are involved in cuticular wax biosynthesis, and PbrKCS28 is involved in suberin biosynthesis, which play essential roles in pear fruit skin formation. Moreover, this work provides a foundation for further understanding the functions of KCS genes in pear.

Genome-wide identification and expression analysis of 3-ketoacyl-CoA synthase gene family in rice (Oryza sativa L.) under cadmium stress

3-Ketoacyl-CoA synthase (KCS) is the key rate-limiting enzyme for the synthesis of very long-chain fatty acids (VLCFAs) in plants, which determines the carbon chain length of VLCFAs. However, a comprehensive study of KCSs in Oryza sativa has not been reported yet. In this study, we identified 22 OsKCS genes in rice, which are unevenly distributed on nine chromosomes. The OsKCS gene family is divided into six subclasses. Many cis-acting elements related to plant growth, light, hormone, and stress response were enriched in the promoters of OsKCS genes. Gene duplication played a crucial role in the expansion of the OsKCS gene family and underwent a strong purifying selection. Quantitative Real-time polymerase chain reaction (qRT-PCR) results revealed that most KCS genes are constitutively expressed. We also revealed that KCS genes responded differently to exogenous cadmium stress in japonica and indica background, and the KCS genes with higher expression in leaves and seeds may have functions under cadmium stress. This study provides a basis for further understanding the functions of KCS genes and the biosynthesis of VLCFA in rice.

Integrated Analysis of lncRNA–mRNA Regulatory Networks Related to Lipid Metabolism in High-Oleic-Acid Rapeseed

A high oleic acid content is considered an essential characteristic in the breeding of high-quality rapeseed in China. Long-chain non-coding RNA (lncRNA) molecules play an important role in the plant’s growth and its response to stress. To better understand the role of lncRNAs in regulating plant reproductive development, we analyzed whole-transcriptome and physiological data to characterize the dynamic changes in lncRNA expression during the four representative times of seed development of high- and low-oleic-acid rapeseed in three regions. We identified 21 and 14 lncRNA and mRNA modules, respectively. These modules were divided into three types related to region, development stages, and material. Next, we analyzed the key modules related to the oil content and the oleic acid, linoleic acid, and linolenic acid contents with physiological data and constructed the key functional network analysis on this basis. Genes related to lipid metabolism, such as 3-ketoacyl-CoA synthase 16 (KCS16) and acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1), were present in the co-expression network, suggesting that the effect of these genes on lipid metabolism might be embodied by the expression of these lncRNAs. Our results provide a fresh insight into region-, development-stage-, and material-biased changes in lncRNA expression in the seeds of Brassica napus. Some of these lncRNAs may participate in the regulatory network of lipid accumulation and metabolism, together with regulated genes. These results may help elucidate the regulatory system of lncRNAs in the lipid metabolism of high-oleic-acid rapeseed seeds.

Genome-wide characterization of phospholipase D family genes in allotetraploid peanut and its diploid progenitors revealed their crucial roles in growth and abiotic stress responses

Abiotic stresses such as cold, drought and salinity are the key environmental factors that limit the yield and quality of oil crop peanut. Phospholipase Ds (PLDs) are crucial hydrolyzing enzymes involved in lipid mediated signaling and have valuable functions in plant growth, development and stress tolerance. Here, 22, 22 and 46 PLD genes were identified in Arachis duranensis, Arachis ipaensis and Arachis hypogaea, respectively, and divided into α, β, γ, δ, ε, ζ and φ isoforms. Phylogenetic relationships, structural domains and molecular evolution proved the conservation of PLDs between allotetraploid peanut and its diploid progenitors. Almost each A. hypogaea PLD except for AhPLDα6B had a corresponding homolog in A. duranensis and A. ipaensis genomes. The expansion of Arachis PLD gene families were mainly attributed to segmental and tandem duplications under strong purifying selection. Functionally, the most proteins interacting with AhPLDs were crucial components of lipid metabolic pathways, in which ahy-miR3510, ahy-miR3513-3p and ahy-miR3516 might be hub regulators. Furthermore, plenty of cis-regulatory elements involved in plant growth and development, hormones and stress responses were identified. The tissue-specific transcription profiling revealed the broad and unique expression patterns of AhPLDs in various developmental stages. The qRT-PCR analysis indicated that most AhPLDs could be induced by specific or multiple abiotic stresses. Especially, AhPLDα3A, AhPLDα5A, AhPLDβ1A, AhPLDβ2A and AhPLDδ4A were highly up-regulated under all three abiotic stresses, whereas AhPLDα9A was neither expressed in 22 peanut tissues nor induced by any abiotic stresses. This genome-wide study provides a systematic analysis of the Arachis PLD gene families and valuable information for further functional study of candidate AhPLDs in peanut growth and abiotic stress responses.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression analysis under Fusarium kyushuense and drought stress conditions in Passiflora edulis

The NAC gene family is one of the largest plant transcription factors (TFs) families and plays important roles in plant growth, development, metabolism, and biotic and abiotic stresses. However, NAC gene family has not been reported in passion fruit (Passiflora edulis). In this study, a total of 105 NAC genes were identified in the passion fruit genome and were unevenly distributed across all nine-passion fruit chromomere, with a maximum of 48 PeNAC genes on chromosome one. The physicochemical features of all 105 PeNAC genes varied including 120 to 3,052 amino acids, 3 to 8 conserved motifs, and 1 to 3 introns. The PeNAC genes were named (PeNAC001-PeNAC105) according to their chromosomal locations and phylogenetically grouped into 15 clades (NAC-a to NAC-o). Most PeNAC proteins were predicted to be localized in the nucleus. The cis-element analysis indicated the possible roles of PeNAC genes in plant growth, development, light, hormones, and stress responsiveness. Moreover, the PeNAC gene duplications including tandem (11 gene pairs) and segmental (12 gene pairs) were identified and subjected to purifying selection. All PeNAC proteins exhibited similar 3D structures, and a protein-protein interaction network analysis with known Arabidopsis proteins was predicted. Furthermore, 17 putative ped-miRNAs were identified to target 25 PeNAC genes. Potential TFs including ERF, BBR-BPC, Dof, and bZIP were identified in promoter region of all 105 PeNAC genes and visualized in a TF regulatory network. GO and KEGG annotation analysis exposed that PeNAC genes were related to different biological, molecular, and cellular terms. The qRT-PCR expression analysis discovered that most of the PeNAC genes including PeNAC001, PeNAC003, PeNAC008, PeNAC028, PeNAC033, PeNAC058, PeNAC063, and PeNAC077 were significantly upregulated under Fusarium kyushuense and drought stress conditions compared to controls. In conclusion, these findings lay the foundation for further functional studies of PeNAC genes to facilitate the genetic improvement of plants to stress resistance.

The Homeodomain-Leucine Zipper Genes Family Regulates the Jinggangmycin Mediated Immune Response of Oryza sativa to Nilaparvata lugens, and Laodelphax striatellus

The homeodomain-leucine zipper (HDZIP) is an important transcription factor family, instrumental not only in growth but in finetuning plant responses to environmental adversaries. Despite the plethora of literature available, the role of HDZIP genes under chewing and sucking insects remains elusive. Herein, we identified 40 OsHDZIP genes from the rice genome database. The evolutionary relationship, gene structure, conserved motifs, and chemical properties highlight the key aspects of OsHDZIP genes in rice. The OsHDZIP family is divided into a further four subfamilies (i.e., HDZIP I, HDZIP II, HDZIP III, and HDZIP IV). Moreover, the protein–protein interaction and Gene Ontology (GO) analysis showed that OsHDZIP genes regulate plant growth and response to various environmental stimuli. Various microRNA (miRNA) families targeted HDZIP III subfamily genes. The microarray data analysis showed that OsHDZIP was expressed in almost all tested tissues. Additionally, the differential expression patterns of the OsHDZIP genes were found under salinity stress and hormonal treatments, whereas under brown planthopper (BPH), striped stem borer (SSB), and rice leaf folder (RLF), only OsHDZIP3, OsHDZIP4, OsHDZIP40, OsHDZIP10, and OsHDZIP20 displayed expression. The qRT-PCR analysis further validated the expression of OsHDZIP20, OsHDZIP40, and OsHDZIP10 under BPH, small brown planthopper (SBPH) infestations, and jinggangmycin (JGM) spraying applications. Our results provide detailed knowledge of the OsHDZIP gene family resistance in rice plants and will facilitate the development of stress-resilient cultivars, particularly against chewing and sucking insect pests.