In Silico Genome-Wide Analysis of B3 Transcription Factors in Cannabis sativa L

Introduction: The B3 transcription factor has been identified in Arabidopsis thaliana, Oryza sativa, and Solanum lycopersicum, among other species. This family of transcription factors regulates seed growth, development, and stress. Cannabis is a valuable crop with numerous applications; however, no B3 transcription factors have been identified in this plant. Materials and Methods: The cannabis B3 gene family was identified and analyzed using bioinformatics analysis tools, such as the NCBI database, plantTFDB website, TBtools, and MEGA software. Quantitative real-time polymerase chain reaction (qRT-PCR) experiments were used to confirm its function. Results: The cannabis B3 family contains 65 members spread across 10 chromosomes. The isoelectric point ranged from 10.03 to 4.65, and the molecular weight ranged from 99,542.88 to 14,310.9 Da. Most of the members were found in the nucleus. The upstream promoter region of the gene contains a variety of cis-acting elements related to the stress response. RNA-seq data and qRT-PCR results showed that CsB3 genes were expressed differently in five organs of female Diku plants and in glandular hairs of nine distinct types of female cannabis inflorescences. Collinearity analysis revealed that there were more homologous genes between cannabis and dicotyledons than monocotyledonous plants, which was consistent with the evolutionary relationship. Conclusions: Hormones and external environmental factors might influence CsB3 expression. Furthermore, some genes such as CsB3-02, CsB3-07, CsB3-50, CsB3-62, and CsB3-65 may participate in cannabis growth and development and play a role in secondary metabolite synthesis. This study provides a solid foundation for further research into the gene function of the cannabis B3 family.

Genome-Wide Identification and Expression Analysis of Sucrose Nonfermenting 1-Related Protein Kinase (SnRK) Genes in Salvia miltiorrhiza in Response to Hormone

The SnRK gene family is the chief component of plant stress resistance and metabolism through activating the phosphorylation of downstream proteins. S. miltiorrhiza is widely used for the treatment of cardiovascular diseases in Asian countries. However, information about the SnRK gene family of S. miltiorrhiza is not clear. The aim of this study is to comprehensively analyze the SnRK gene family of S. miltiorrhiza and its response to phytohormone. Here, 33 SmSnRK genes were identified and divided into three subfamilies (SmSnRK1, SmSnRK2 and SmSnRK3) according to phylogenetic analysis and domain. SmSnRK genes within same subgroup shared similar protein motif composition and were unevenly distributed on eight chromosomes of S. miltiorrhiza. Cis-acting element analysis showed that the promoter of SmSnRK genes was enriched with ABRE motifs. Expression pattern analysis revealed that SmSnRK genes were preferentially expressed in leaves and roots. Most SmSnRK genes were induced by ABA and MeJA treatment. Correlation analysis showed that SmSnRK3.15 and SmSnRK3.18 might positively regulate tanshinone biosynthesis; SmSnRK3.10 and SmSnRK3.12 might positively regulate salvianolic acid biosynthesis. RNAi-based silencing of SmSnRK2.6 down-regulated the biosynthesis of tanshinones and biosynthetic genes expression. An in vitro phosphorylation assay verified that SmSnRK2.2 interacted with and phosphorylated SmAREB1. These findings will provide a valuable basis for the functional characterization of SmSnRK genes and quality improvement of S. miltiorrhiza.

Genome-Wide Identification and Expression Pattern Analysis of BAHD Acyltransferase Family in Taxus mairei

BAHD acyltransferases are involved in catalyzing and regulating the secondary metabolism in plants. Despite this, the members of BAHD family and their functions have not been reported in the Taxus species. In this study, a total of 123 TwBAHD acyltransferases from Taxus wallichiana var. mairei genome were identified and divided into six clades based on phylogenetic analysis, of which Clade VI contained a Taxus-specific branch of 52 members potentially involved in taxol biosynthesis. Most TwBAHDs from the same clade shared similar conserved motifs and gene structures. Besides the typical conserved motifs within the BAHD family, the YPLAGR motif was also conserved in multiple clades of T. mairei. Moreover, only one pair of tandem duplicate genes was found on chromosome 1, with a Ka/Ks ratio < 1, indicating that the function of duplicate genes did not differentiate significantly. RNA-seq analysis revealed different expression patterns of TwBAHDs in MeJA induction and tissue-specific expression experiments. Several TwBAHD genes in the Taxus-specific branch were highly expressed in different tissues of T. mairei, suggesting an important role in the taxol pathway. This study provides comprehensive information for the TwBAHD gene family and sets up a basis for its potential functions.

Genome-Wide Investigation of Class III Peroxidase Genes in Brassica napus Reveals Their Responsiveness to Abiotic Stresses

Brassica napus (B. napus) is susceptible to multiple abiotic stresses that can affect plant growth and development, ultimately reducing crop yields. In the past, many genes that provide tolerance to abiotic stresses have been identified and characterized. Peroxidase (POD) proteins, members of the oxidoreductase enzyme family, play a critical role in protecting plants against abiotic stresses. This study demonstrated a comprehensive investigation of the POD gene family in B. napus. As a result, a total of 109 POD genes were identified across the 19 chromosomes and classified into five distinct subgroups. Further, 44 duplicate events were identified; of these, two gene pairs were tandem and 42 were segmental. Synteny analysis revealed that segmental duplication was more prominent than tandem duplication among POD genes. Expression pattern analysis based on the RNA-seq data of B. napus indicated that BnPOD genes were expressed differently in various tissues; most of them were expressed in roots rather than in other tissues. To validate these findings, we performed RT-qPCR analysis on ten genes; these genes showed various expression levels under abiotic stresses. Our findings not only furnish valuable insights into the evolutionary dynamics of the BnPOD gene family but also serve as a foundation for subsequent investigations into the functional roles of POD genes in B. napus.

Combined genome and transcriptome analysis of elite fiber quality in Gossypium barbadense

Gossypium barbadense, which is one of several species of cotton, is well-known for its superior fiber quality. However, the genetic basis of its high-quality fiber remains largely unexplored. Here, we re-sequenced 269 G. barbadense accessions. Phylogenetic structure analysis showed that the set of accessions were clustered into three groups, G1 and G2 mainly included modern cultivars from Xinjiang, China, and G3 was related to widely introduced accessions in different regions worldwide. A genome-wide association study of five fiber quality traits across multiple field environments identified a total of 512 qtls (main effect QTLs) and 94 qtlEs (QTL-by-environment interactions) related to fiber quality, of which 292 qtls and 57 qtlEs co-located with previous studies. We extracted the genes located in these loci and performed expression comparison, local association analysis, and introgression segment identification. The results showed that high expression of hormone-related genes during fiber development, introgressions from G. hirsutum, and the recombination of domesticated elite allelic variation, were three major contributors to improve fiber quality of G. barbadense. In total, 839 candidate genes with encoding region variation associated with elite fiber quality were mined. We confirmed that haplotype GB_D03G0092H traced from G. hirsutum introgression, with a 1-bp deletion leading to a frameshift mutation compared with GB_D03G0092B, significantly improved fiber quality. GB_D03G0092H localized in the plasma membrane, while GB_D03G0092B in both the nucleus and plasma membrane. Overexpression of GB_D03G0092H in Arabidopsis (Arabidopsis thaliana) significantly improved the elongation of longitudinal cells. Our study systematically reveals the genetic basis of the superior fiber quality of G. barbadense and provides elite segments and gene resources for breeding high-quality cotton cultivars.

Identification and characterization of two O-methyltransferases involved in methylated 2-(2-phenethyl) chromones biosynthesis in agarwood

The 2-(2-phenethyl)chromones (PECs) are the signature constituents responsible for the fragrance and pharmacological properties of agarwood. O-methyltransferases (OMTs) are necessary for the methylated PECs biosynthesis, however, there is little information known about OMTs in Aquilaria sinensis. In this study, we identified 29 OMT genes from the A. sinensis genome. Expression analysis showed they were differentially expressed in different tissues and responded to drill wounding. Comprehensive analysis of the gene expression and methylated PECs content revealed that the AsOMT2, AsOMT8, AsOMT11, AsOMT16, and AsOMT28 could potentially be involved in methylated PECs biosynthesis. The vitro enzyme assays and functional analysis in Nicotiana benthamiana demonstrated that AsOMT11 and AsOMT16 could methylate 6-hydroxy-2-(2-phenylethyl)chromone to form 6-methoxy-2-(2-phenylethyl) chromone. The transient overexpression experiment in Qi-Nan revealed that AsOMT11 and AsOMT16 could significantly promote the accumulation of three major methylated PECs. Our results provide candidate genes for the mass production of methylated PECs using synthetic biology.

Research Note: SOCS2 contributes to reduction of the third digit during development of the chicken forelimb

The development of the avian wing pattern has been the subject of heated debate due to its special shape. The Suppressor of cytokine signaling 2 (SOCS2) gene encodes a negative regulator of growth hormone (GH) signaling and bone growth and is known to be strongly expressed in the third digit of chicken forelimbs. These observations suggest that SOCS2 might regulate the morphology of the avian wing, however, the function of SOCS2 in avian limb development remains unknown. Here, we reexamined SOCS2 expression in successive developmental stages of chicken limb development by in situ hybridization (ISH) and describe extended expression from the posterior of the stypolod to the third digit of the forelimbs. We used the RCAS avian retrovirus to overexpress SOCS2 in the developing chicken limb buds, which resulted in reduced or malformed chicken wings while hindlimbs developed normally. Transcriptome sequencing (mRNA-Seq) revealed changes in expression of genes known to be associated with growth and development in forelimbs with overexpressed SOCS2. This study highlights a pivotal role for SOCS2 during the development of the wing in the chicken and provides new insight into molecular mechanisms regulating avian limb development.

Genome-Wide Identification and Co-Expression Networks of WOX Gene Family in Nelumbo nucifera

WUSCHEL-related homeobox (WOX) genes are a class of plant-specific transcription factors, regulating the development of multiple tissues. However, the genomic characterizations and expression patterns of WOX genes have not been analyzed in lotus. In this study, 15 NnWOX genes were identified based on the well-annotated reference genome of lotus. According to the phylogenetic analysis, the NnWOX genes were clustered into three clades, i.e., ancient clade, intermediate clade, and WUS clade. Except for the conserved homeobox motif, we further found specific motifs of NnWOX genes in different clades and divergence gene structures, suggesting their distinct functions. In addition, two NnWOX genes in the ancient clade have conserved expression patterns and other NnWOX genes exhibit different expression patterns in lotus tissues, suggesting a low level of functional redundancy in lotus WOX genes. Furthermore, we constructed the gene co-expression networks for each NnWOX gene. Based on weighted gene co-expression network analysis (WGCNA), ten NnWOX genes and their co-expressed genes were assigned to the modules that were significantly related to the cotyledon and seed coat. We further performed RT-qPCR experiments, validating the expression levels of ten NnWOX genes in the co-expression networks. Our study reveals comprehensive genomic features of NnWOX genes in lotus, providing a solid basis for further function studies.

Identification and Analysis of the MIR399 Gene Family in Grapevine Reveal Their Potential Functions in Abiotic Stress

MiR399 plays an important role in plant growth and development. The objective of the present study was to elucidate the evolutionary characteristics of the MIR399 gene family in grapevine and investigate its role in stress response. To comprehensively investigate the functions of miR399 in grapevine, nine members of the Vvi-MIR399 family were identified based on the genome, using a miRBase database search, located on four chromosomes (Chr 2, Chr 10, Chr 15, and Chr 16). The lengths of the Vvi-miR399 precursor sequences ranged from 82 to 122 nt and they formed stable stem-loop structures, indicating that they could produce microRNAs (miRNAs). Furthermore, our results suggested that the 2 to 20 nt region of miR399 mature sequences were relatively conserved among family members. Phylogenetic analysis revealed that the Vvi-MIR399 members of dicots (Arabidopsis, tomato, and sweet orange) and monocots (rice and grapevine) could be divided into three clades, and most of the Vvi-MIR399s were closely related to sweet orange in dicots. Promoter analysis of Vvi-MIR399s showed that the majority of the predicted cis-elements were related to stress response. A total of 66.7% (6/9) of the Vvi-MIR399 promoters harbored drought, GA, and SA response elements, and 44.4% (4/9) of the Vvi-MIRR399 promoters also presented elements involved in ABA and MeJA response. The expression trend of Vvi-MIR399s was consistent in different tissues, with the lowest expression level in mature and young fruits and the highest expression level in stems and young leaves. However, nine Vvi-MIR399s and four target genes showed different expression patterns when exposed to low light, high light, heat, cold, drought, and salt stress. Interestingly, a putative target of Vvi-MIR399 targeted multiple genes; for example, seven Vvi-MIR399s simultaneously targeted VIT_213s0067g03280.1. Furthermore, overexpression of Vvi_MIR399e and Vvi_MIR399f in Arabidopsis enhanced tolerance to drought compared with wild-type (WT). In contrast, the survival rate of Vvi_MIR399d-overexpressed plants were zero after drought stress. In conclusion, Vvi-MIR399e and Vvi-MIR399f, which are related to drought tolerance in grapevine, provide candidate genes for future drought resistance breeding.

Transcriptome-Based Identification of the SaR2R3-MYB Gene Family in Sophora alopecuroides and Function Analysis of SaR2R3-MYB15 in Salt Stress Tolerance

As one of the most prominent gene families, R2R3-MYB transcription factors significantly regulate biochemical and physiological processes under salt stress. However, in Sophora alopecuroides, a perennial herb known for its exceptional saline alkali resistance, the comprehensive identification and characterization of SaR2R3-MYB genes and their potential functions in response to salt stress have yet to be determined. We investigated the expression profiles and biological functions of SaR2R3-MYB transcription factors in response to salt stress, utilizing a transcriptome-wide mining method. Our analysis identified 28 SaR2R3-MYB transcription factors, all sharing a highly conserved R2R3 domain, which were further divided into 28 subgroups through phylogenetic analysis. Some SaR2R3-MYB transcription factors showed induction under salt stress, with SaR2R3-MYB15 emerging as a potential regulator based on analysis of the protein-protein interaction network. Validation revealed the transcriptional activity and nuclear localization of SaR2R3-MYB15. Remarkably, overexpression of SaR2R3-MYB15 in transgenic plants could increase the activity of antioxidant enzymes and the accumulation of proline but decrease the content of malondialdehyde (MDA), compared with wild-type plants. Moreover, several salt stress-related genes showed higher expression levels in transgenic plants, implying their potential to enhance salt tolerance. Our findings shed light on the role of SaR2R3-MYB genes in salt tolerance in S. alopecuroides.

Comprehensive analysis of the aldehyde dehydrogenase gene family in Phaseolus vulgaris L. and their response to saline-alkali stress

Background

Aldehyde dehydrogenase (ALDH) scavenges toxic aldehyde molecules by catalyzing the oxidation of aldehydes to carboxylic acids. Although ALDH gene family members in various plants have been extensively studied and were found to regulate plant response to abiotic stress, reports on ALDH genes in the common bean (Phaseolus vulgaris L.) are limited. In this study, we aimed to investigate the effects of neutral (NS) and basic alkaline (AS) stresses on growth, physiological and biochemical indices, and ALDH activity, ALDH gene expression of common bean. In addition, We used bioinformatics techniques to analyze the physical and chemical properties, phylogenetic relationships, gene replication, collinearity, cis-acting elements, gene structure, motifs, and protein structural characteristics of PvALDH family members.

Results

We found that both NS and AS stresses weakened the photosynthetic performance of the leaves, induced oxidative stress, inhibited common bean growth, and enhanced the antioxidative system to scavenge reactive oxygen species. Furthermore, we our findings revealed that ALDH in the common bean actively responds to NS or AS stress by inducing the expression of PvALDH genes. In addition, using the established classification criteria and phylogenetic analysis, 27 PvALDHs were identified in the common bean genome, belonging to 10 ALDH families. The primary expansion mode of PvALDH genes was segmental duplication. Cis-acting elemental analysis showed that PvALDHs were associated with abiotic stress and phytohormonal responses. Gene expression analysis revealed that the PvALDH gene expression was tissue-specific. For instance, PvALDH3F1 and PvALDH3H1 were highly expressed in flower buds and flowers, respectively, whereas PvALDH3H2 and PvALDH2B4 were highly expressed in green mature pods and young pods, respectively. PvALDH22A1 and PvALDH11A2 were highly expressed in leaves and young trifoliates, respectively; PvALDH18B2 and PvALDH18B3 were highly expressed in stems and nodules, respectively; and PvALDH2C2 and PvALDH2C3 were highly expressed in the roots. PvALDHs expression in the roots responded positively to NS-AS stress, and PvALDH2C3, PvALDH5F1, and PvALDH10A1 were significantly (P < 0.05) upregulated in the roots.

Conclusion

These results indicate that AS stress causes higher levels of oxidative damage than NS stress, resulting in weaker photosynthetic performance and more significant inhibition of common bean growth. The influence of PvALDHs potentially modulates abiotic stress response, particularly in the context of saline-alkali stress. These findings establish a basis for future research into the potential roles of ALDHs in the common bean.

Genome-Wide Analysis of Aux/IAA Gene Family in Artemisia argyi: Identification, Phylogenetic Analysis, and Determination of Response to Various Phytohormones

Artemisia argyi is a traditional herbal medicine plant, and its folium artemisia argyi is widely in demand due to moxibustion applications globally. The Auxin/indole-3-acetic acid (Aux/IAA, or IAA) gene family has critical roles in the primary auxin-response process, with extensive involvement in plant development and stresses, controlling various essential traits of plants. However, the systematic investigation of the Aux/IAA gene family in A. argyi remains limited. In this study, a total of 61 Aux/IAA genes were comprehensively identified and characterized. Gene structural analysis indicated that 46 Aux/IAA proteins contain the four typical domains, and 15 Aux/IAA proteins belong to non-canonical IAA proteins. Collinear prediction and phylogenetic relationship analyses suggested that Aux/IAA proteins were grouped into 13 distinct categories, and most Aux/IAA genes might experience gene loss during the tandem duplication process. Promoter cis-element investigation indicated that Aux/IAA promoters contain a variety of plant hormone response and stress response cis-elements. Protein interaction prediction analysis demonstrated that AaIAA26/29/7/34 proteins are possibly core members of the Aux/IAA family interaction. Expression analysis in roots and leaves via RNA-seq data indicated that the expression of some AaIAAs exhibited tissue-specific expression patterns, and some AaIAAs were involved in the regulation of salt and saline-alkali stresses. In addition, RT-qPCR results indicated that AaIAA genes have differential responses to auxin, with complex response patterns in response to other hormones, indicating that Aux/IAA may play a role in connecting auxin and other hormone signaling pathways. Overall, these findings shed more light on AaIAA genes and offer critical foundational knowledge toward the elucidation of their function during plant growth, stress response, and hormone networking of Aux/IAA family genes in A. argyi.

Genome-Wide Analysis of C-Repeat Binding Factor Gene Family in Capsicum baccatum and Functional Exploration in Low-Temperature Response

Capsicum baccatum is a close relative of edible chili peppers (Capsicum annuum) with high economic value. The CBF gene family plays an important role in plant stress resistance physiology. We detected a total of five CBF genes in the C. baccatum genome-wide sequencing data. These genes were scattered irregularly across four chromosomes. The genes were categorized into three groupings according to their evolutionary relationships, with genes in the same category showing comparable principles for motif composition. The 2000 bp upstream of CbCBF contains many resistance-responsive elements, hormone-responsive elements, and transcription factor binding sites. These findings emphasize the crucial functions of these genes in responding to challenging conditions and physiological regulation. Analysis of tissue-specific expression revealed that CbCBF3 exhibited the greatest level of expression among all tissues. Under conditions of low-temperature stress, all CbCBF genes exhibited different levels of responsiveness, with CbCBF3 showing a considerable up-regulation after 0.25 h of cold stress, indicating a high sensitivity to low-temperature response. The importance of the CbCBF3 gene in the cold response of C. baccatum was confirmed by the use of virus-induced gene silencing (VIGS) technology, as well as the prediction of its protein interaction network. To summarize, this study conducts a thorough bioinformatics investigation of the CbCBF gene family, showcases the practicality of employing VIGS technology in C. baccatum, and confirms the significance of the CbCBF3 gene in response to low temperatures. These findings provide significant references for future research on the adaptation of C. baccatum to low temperatures.

Genome-Wide Mining of CULLIN E3 Ubiquitin Ligase Genes from Uncaria rhynchophylla

CULLIN (CUL) protein is a subtype of E3 ubiquitin ligase that is involved in a variety of biological processes and responses to stress in plants. In Uncaria rhynchophylla, the CUL gene family has not been identified and its role in plant development, stress response and secondary metabolite synthesis has not been studied. In this study, 12 UrCUL gene members all contained the typical N-terminal domain and C-terminal domain identified from the U. rhynchophylla genome and were classified into four subfamilies based on the phylogenetic relationship with CULs in Arabidopsis thaliana. They were unevenly distributed on eight chromosomes but had a similar structural composition in the same subfamily, indicating that they were relatively conserved and potentially had similar gene functions. An interspecific and intraspecific collinearity analysis showed that fragment duplication played an important role in the evolution of the CUL gene family. The analysis of the cis-acting elements suggests that the UrCULs may play an important role in various biological processes, including the abscisic acid (ABA) response. To investigate this hypothesis, we treated the roots of U. rhynchophylla tissue-cultured seedlings with ABA. The expression pattern analysis showed that all the UrCUL genes were widely expressed in roots with various expression patterns. The co-expression association analysis of the UrCULs and key enzyme genes in the terpenoid indole alkaloid (TIA) synthesis pathway revealed the complex expression patterns of 12 UrCUL genes and some key TIA enzyme genes, especially UrCUL1, UrCUL1-likeA, UrCUL2-likeA and UrCUL2-likeB, which might be involved in the biosynthesis of TIAs. The results showed that the UrCULs were involved in the response to ABA hormones, providing important information for elucidating the function of UrCULs in U. rhynchophylla. The mining of UrCULs in the whole genome of U. rhynchophylla provided new information for understanding the CUL gene and its function in plant secondary metabolites, growth and development.

Genome-Wide Analysis of the Oat (Avena sativa) HSP90 Gene Family Reveals Its Identification, Evolution, and Response to Abiotic Stress

Oats (Avena sativa) are an important cereal crop and cool-season forage worldwide. Heat shock protein 90 (HSP90) is a protein ubiquitously expressed in response to heat stress in almost all plants. To date, the HSP90 gene family has not been comprehensively reported in oats. Herein, we have identified twenty HSP90 genes in oats and elucidated their evolutionary pathways and responses to five abiotic stresses. The gene structure and motif analyses demonstrated consistency across the phylogenetic tree branches, and the groups exhibited relative structural conservation. Additionally, we identified ten pairs of segmentally duplicated genes in oats. Interspecies synteny analysis and orthologous gene identification indicated that oats share a significant number of orthologous genes with their ancestral species; this implies that the expansion of the oat HSP90 gene family may have occurred through oat polyploidization and large fragment duplication. The analysis of cis-acting elements revealed their influential role in the expression pattern of HSP90 genes under abiotic stresses. Analysis of oat gene expression under high-temperature, salt, cadmium (Cd), polyethylene glycol (PEG), and abscisic acid (ABA) stresses demonstrated that most AsHSP90 genes were significantly up-regulated by heat stress, particularly AsHSP90-7, AsHSP90-8, and AsHSP90-9. This study offers new insights into the amplification and evolutionary processes of the AsHSP90 protein, as well as its potential role in response to abiotic stresses. Furthermore, it lays the groundwork for understanding oat adaptation to abiotic stress, contributing to research and applications in plant breeding.

Genome-Wide Identification and Expression Profile Analysis of Sugars Will Eventually Be Exported Transporter (SWEET) Genes in Zantedeschia elliottiana and Their Responsiveness to Pectobacterium carotovora subspecies Carotovora (Pcc) Infection

SWEET, sugars will eventually be exported transporter, is a novel class of sugar transporter proteins that can transport sugars across membranes down a concentration gradient. It plays a key role in plant photosynthetic assimilates, phloem loading, nectar secretion from nectar glands, seed grouting, pollen development, pathogen interactions, and adversity regulation, and has received widespread attention in recent years. To date, systematic analysis of the SWEET family in Zantedeschia has not been documented, although the genome has been reported in Zantedeschia elliottiana. In this study, 19 ZeSWEET genes were genome-wide identified in Z. elliottiana, and unevenly located in 10 chromosomes. They were further clustered into four clades by a phylogenetic tree, and almost every clade has its own unique motifs. Synthetic analysis confirmed two pairs of segmental duplication events of ZeSWEET genes. Heatmaps of tissue-specific and Pectobacterium carotovora subsp. Carotovora (Pcc) infection showed that ZeSWEET genes had different expression patterns, so SWEETs may play widely varying roles in development and stress tolerance in Zantedeschia. Moreover, quantitative reverse transcription-PCR (qRT-PCR) analysis revealed that some of the ZeSWEETs responded to Pcc infection, among which eight genes were significantly upregulated and six genes were significantly downregulated, revealing their potential functions in response to Pcc infection. The promoter sequences of ZeSWEETs contained 51 different types of the 1380 cis-regulatory elements, and each ZeSWEET gene contained at least two phytohormone responsive elements and one stress response element. In addition, a subcellular localization study indicated that ZeSWEET07 and ZeSWEET18 were found to be localized to the plasma membrane. These findings provide insights into the characteristics of SWEET genes and contribute to future studies on the functional characteristics of ZeSWEET genes, and then improve Pcc infection tolerance in Zantedeschia through molecular breeding.

Genome-Wide Identification and Expression Analysis of DWARF53 Gene in Response to GA and SL Related to Plant Height in Banana

Dwarfing is one of the common phenotypic variations in asexually reproduced progeny of banana, and dwarfed banana is not only windproof and anti-fallout but also effective in increasing acreage yield. As a key gene in the strigolactone signalling pathway, DWARF53 (D53) plays an important role in the regulation of the height of plants. In order to gain insight into the function of the banana D53 gene, this study conducted genome-wide identification of banana D53 gene based on the M. acuminata, M. balbisiana and M. itinerans genome database. Analysis of MaD53 gene expression under high temperature, low temperature and osmotic stress based on transcriptome data and RT-qPCR was used to analyse MaD53 gene expression in different tissues as well as in different concentrations of GA and SL treatments. In this study, we identified three MaD53, three MbD53 and two MiD53 genes in banana. Phylogenetic tree analysis showed that D53 Musa are equally related to D53 Asparagales and Poales. Both high and low-temperature stresses substantially reduced the expression of the MaD53 gene, but osmotic stress treatments had less effect on the expression of the MaD53 gene. GR24 treatment did not significantly promote the height of the banana, but the expression of the MaD53 gene was significantly reduced in roots and leaves. GA treatment at 100 mg/L significantly promoted the expression of the MaD53 gene in roots, but the expression of this gene was significantly reduced in leaves. In this study, we concluded that MaD53 responds to GA and SL treatments, but "Yinniaijiao" dwarf banana may not be sensitive to GA and SL.

Genome-Wide Identification and Characterization of the TIFY Gene Family and Their Expression Patterns in Response to MeJA and Aluminum Stress in Centipedegrass (Eremochloa ophiuroides)

The TIFY family is a group of novel plant-specific transcription factors involved in plant development, signal transduction, and responses to stress and hormones. TIFY genes have been found and functionally characterized in a number of plant species. However, there is no information about this family in warm-season grass plants. The current study identified 24 TIFY genes in Eremochloa ophiuroides, a well-known perennial warm-season grass species with a high tolerance to aluminum toxicity and good adaptability to the barren acidic soils. All of the 24 EoTIFYs were unevenly located on six out of nine chromosomes and could be classified into two subfamilies (ZIM/ZML and JAZ), consisting of 3 and 21 genes, respectively, with the JAZ subfamily being further divided into five subgroups (JAZ I to JAZ V). The amino acids of 24 EoTIFYs showed apparent differences between the two subfamilies based on the analysis of gene structures and conserved motifs. MCScanX analysis revealed the tandem duplication and segmental duplication of several EoTIFY genes occurred during E. ophiuroides genome evolution. Syntenic analyses of TIFY genes between E. ophiuroides and other five plant species (including A. thaliana, O. sativa, B. distachyon, S. biocolor, and S. italica) provided valuable clues for understanding the potential evolution of the EoTIFY family. qRT-PCR analysis revealed that EoTIFY genes exhibited different spatial expression patterns in different tissues. In addition, the expressions of EoTIFY genes were highly induced by MeJA and all of the EoTIFY family members except for EoJAZ2 displayed upregulated expression by MeJA. Ten EoTIFY genes (EoZML1, EoZML1, EoJAZ1, EoJAZ3, EoJAZ5, EoJAZ6, EoJAZ8, EoJAZ9, EoJAZ10, and EoJAZ21) were observed to be highly expressed under both exogenous MeJA treatment and aluminum stress, respectively. These results suggest that EoTIFY genes play a role in the JA-regulated pathway of plant growth and aluminum resistance as well. The results of this study laid a foundation for further understanding the function of TIFY genes in E. ophiuroides, and provided useful information for future aluminum tolerance related breeding and gene function research in warm-season grass plants.

Genome-wide identification and characterization of TCP gene family in Dendrobium nobile and their role in perianth development

TCP is a widely distributed, essential plant transcription factor that regulates plant growth and development. An in-depth study of TCP genes in Dendrobium nobile, a crucial parent in genetic breeding and an excellent model material to explore perianth development in Dendrobium, has not been conducted. We identified 23 DnTCP genes unevenly distributed across 19 chromosomes and classified them as Class I PCF (12 members), Class II: CIN (10 members), and CYC/TB1 (1 member) based on the conserved domain and phylogenetic analysis. Most DnTCPs in the same subclade had similar gene and motif structures. Segmental duplication was the predominant duplication event for TCP genes, and no tandem duplication was observed. Seven genes in the CIN subclade had potential miR319 and -159 target sites. Cis-acting element analysis showed that most DnTCP genes contained many developmental stress-, light-, and phytohormone-responsive elements in their promoter regions. Distinct expression patterns were observed among the 23 DnTCP genes, suggesting that these genes have diverse regulatory roles at different stages of perianth development or in different organs. For instance, DnTCP4 and DnTCP18 play a role in early perianth development, and DnTCP5 and DnTCP10 are significantly expressed during late perianth development. DnTCP17, 20, 21, and 22 are the most likely to be involved in perianth and leaf development. DnTCP11 was significantly expressed in the gynandrium. Specially, MADS-specific binding sites were present in most DnTCP genes putative promoters, and two Class I DnTCPs were in the nucleus and interacted with each other or with the MADS-box. The interactions between TCP and the MADS-box have been described for the first time in orchids, which broadens our understanding of the regulatory network of TCP involved in perianth development in orchids.

The Molecular and Functional Characterization of Sensory Neuron Membrane Protein 1b (SNMP1b) from Cyrtotrachelus buqueti (Coleoptera: Curculionidae)

Sensory neuron membrane proteins (SNMPs) play important roles in insect chemoreception and SNMP1s have been reported to be essential in detecting sex pheromones in Drosophila and some lepidopteran species. However, SNMPs for Cyrtotrachelus buqueti (Coleoptera: Curculionidae), a major insect pest of bamboo plantations, remain uncharacterized. In this study, a novel SNMP gene, CbuqSNMP1b, from C. buqueti was functionally characterized. The expression of CbuqSNMP1b was significantly higher in antennae than in other tissues of both sexes and the expression level was significantly male-biased. Additionally, CbuqSNMP1b showed significantly higher transcription levels in the adult stage and very low transcription levels in other stages, suggesting that CbuqSNMP1b is involved in the process of olfaction. Fluorescence binding assays indicated that CbuqSNMP1b displayed the strongest binding affinity to dibutyl phthalate (Ki = 9.03 μM) followed by benzothiazole (Ki = 11.59 μM) and phenol (Ki = 20.95 μM) among fourteen C. buqueti volatiles. Furthermore, molecular docking revealed key residues in CbuqSNMP1b that interact with dibutyl phthalate, benzothiazole, and phenol. In conclusion, these findings will lay a foundation to further understand the olfactory mechanisms of C. buqueti and promote the development of novel methods for controlling this pest.

Role of TRAF6 from obscure puffer (Takifugu obscurus) in immune response against Edwardsiella tarda infection

Tumor necrosis factor receptor-associated factor 6 (TRAF6) is a vital molecule of inflammatory signaling pathways in innate immune response against pathogens. To elucidate its role in defense against Edwardsiella tarda infection in teleost fish, TRAF6 homologue was identified from obscure puffer (Takifugu obscurus) and functionally analyzed in this study. The obscure puffer TRAF6 (ToTRAF6) is a protein of 565 amino acids containing conserved RING domain, zinc finger-TRAF and MATH_TRAF6 domain. ToTRAF6 mRNA distributed in various healthy tissues of obscure puffer and was upregulated in the immune related tissues after E. tarda infection. ToTRAF6 protein was localized in the cytoplasm and aggregate as dots around the nuclei in FHM cells. The overexpression of ToTRAF6 in FHM cells decreased the quantity of E. tarda and induced the significant upregulation of downstream MAPK signaling pathway genes. These data suggest that ToTRAF6 is a key molecule of MAPK signaling pathway in defense against E. tarda infection.

Transcriptome-Wide Identification and Expression Analysis of bHLH Family Genes in Iris domestica under Drought and Cu Stress

The role of bHLH transcription factors in plant response to abiotic stress and regulation of flavonoid metabolism is well documented. However, to date, the bHLH transcription factor family in Iris domestica remains unreported, impeding further research on flavonoid metabolism in this plant. To address this knowledge gap, we employed bioinformatics to identify 39 IdbHLH genes and characterised their phylogenetic relationships and gene expression patterns under both drought and copper stress conditions. Our evolutionary tree analysis classified the 39 IdbHLHs into 17 subfamilies. Expression pattern analysis revealed that different IdbHLH transcription factors had distinct expression trends in various organs, suggesting that they might be involved in diverse biological processes. We found that IdbHLH36 was highly expressed in all organs (Transcripts Per Million (TPM) > 10), while only 12 IdbHLH genes in the rhizome and four in the root were significantly upregulated under drought stress. Of these, four genes (IdbHLH05, -37, -38, -39) were co-upregulated in both the rhizome and root, indicating their potential role in drought resistance. With regards to copper stress, we found that only 12 genes were upregulated. Further co-expression analysis revealed that most bHLH genes were significantly correlated with key enzyme genes involved in isoflavone biosynthesis. Thereinto, IdbHLH06 showed a significant positive correlation with IdC4H1 and Id4CL1 (p < 0.05). Furthermore, a transient expression assay confirmed that the IdbHLH06 protein was localised in the nucleus. Our findings provide new insights into the molecular basis and regulatory mechanisms of bHLH transcription factors in isoflavone biosynthesis in I. domestica.

Rice tetraspanins express in specific domains of diverse tissues and regulate plant architecture and root growth

Tetraspanins (TETs) are small transmembrane scaffold proteins that distribute proteins into highly organized microdomains, consisting of adaptors and signaling proteins, which play important roles in various biological events. In plants, understanding of tetraspanin is limited to the Arabidopsis TET genes' expression pattern and their function in leaf and root growth. Here, we comprehensively analyzed all rice tetraspanin (OsTET) family members, including their gene expression pattern, protein topology, and subcellular localization. We found that the core domain of OsTETs is conserved and shares a similar topology of four membrane-spanning domains with animal and plant TETs. OsTET genes are partially overlapping expressed in diverse tissue domains in vegetative and reproductive organs. OsTET proteins preferentially targeted the endoplasmic reticulum. Mutation analysis showed that OsTET5, OsTET6, OsTET9, and OsTET10 regulated plant height and tillering, and that OsTET13 controlled root growth in association with the jasmonic acid pathway. In summary, our work provides systematic new insights into the function of OsTETs in rice growth and development, and the data provides valuable resources for future research.

A Call for Gene Expression Analysis in Whole Blood of Patients With Rheumatoid Arthritis (RA) as a Biomarker for RA-Associated Interstitial Lung Disease

OBJECTIVE

Rheumatoid arthritis (RA)-associated interstitial lung disease (ILD) is one of the most common and prognostic organ manifestations of RA. Therefore, to allow effective treatment, it is of crucial importance to diagnose RA-ILD at the earliest possible stage. So far, the gold standard of early detection has been high-resolution computed tomography (HRCT) of the lungs. This procedure involves considerable radiation exposure for the patient and is therefore unsuitable as a routine screening measure for ethical reasons. Here, we propose the analysis of characteristic gene expression patterns as a biomarker to aid in the early detection and initiation of appropriate, possibly antifibrotic, therapy.

METHODS

To investigate unique molecular patterns of RA-ILD, whole blood samples were taken from 12 female patients with RA-ILD (n = 7) or RA (n = 5). The RNA was extracted, sequenced by RNA-Seq, and analyzed for characteristic differences in the gene expression patterns between patients with RA-ILD and those with RA without ILD.

RESULTS

The differential gene expression analysis revealed 9 significantly upregulated genes in RA-ILD compared to RA without ILD: arginase 1 (ARG1), thymidylate synthetase (TYMS), sortilin 1 (SORT1), marker of proliferation Ki-67 (MKI67), olfactomedin 4 (OLFM4), baculoviral inhibitor of apoptosis repeat containing 5 (BIRC5), membrane spanning 4-domains A4A (MS4A4A), C-type lectin domain family 12 member A (CLEC12A), and the long intergenic nonprotein coding RNA (LINC02967).

CONCLUSION

All gene products of these genes (except for LINC02967) are known from the literature to be involved in the pathogenesis of fibrosis. Further, for some, a contribution to the development of pulmonary fibrosis has even been demonstrated in experimental studies. Therefore, the results presented here provide an encouraging perspective for using specific gene expression patterns as biomarkers for the early detection and differential diagnosis of RA-ILD as a routine screening test.

Involvement of Sep38β in the Insecticidal Activity of Bacillus thuringiensis against Beet Armyworm, Spodoptera exigua (Lepidoptera)

The p38 mitogen-activated protein kinases (MAPKs) are associated with insect immunity, tissue repair, and the insecticidal activity of Bacillus thuringiensis (Bt). Here, a p38 MAPK family gene (Sep38β) was identified from Spodoptera exigua. Among the developmental stages, the transcription level of Sep38β was the highest in egg, followed by that in prepupa and pupa. Sep38β expression peaked in Malpighian tubules and the hemolymph of fifth instar larvae. Knockdown of Sep38β or injection of SB203580 (a p38 MAPK inhibitor) significantly downregulated the SeDUOX expression and reactive oxygen species (ROS) level in the midgut, accounting for deterioration of the midgut to scavenge pathogens and enhancement of Bt insecticidal activity. In conclusion, all the results demonstrate that Sep38β regulates the immune-related ROS level in the insect midgut, which suppresses the insecticidal activity of Bt against S. exigua by 17-22%. Our study highlights that Sep38β is essential for insect immunity and the insecticidal activity of Bt to S. exigua and is a potential target for pest control.

[Identification and expression analysis of the YABBY gene family in strawberry]

YABBY proteins are important transcription factors that regulate morphogenesis and organ development in plants. In order to study the YABBY of strawberry, bioinformatic technique were used to identify the YABBY gene families in Fragaria vesca (diploid) and Fragaria×ananassa (octoploid), and then analyze the sequence characters, phylogeny and collinearity of the family members. The RNA-seq data and the quantitative reverse transcription-polymerase chain reaction (qRT-PCR) technique were used to assay the expression patterns of the family members. A green fluorescent protein (GFP) was fused with FvYABBYs and transiently expressed in tobacco leaf cells for the subcellular localization. As the results, six FvYABBY genes and 26 FxaYABBY genes were identified from F. vesca and F.×ananassa, respectively. The FvYABBY genes were grouped into five clades, and five family members were orthologous with AtYABBY genes of Arabidopsis. In F. vesca, all of the FvYABBYs were basically not expressed not expressed in root and receptacle, while FvYABBY1, FvYABBY2, FvYABBY5 and FvYABBY6 were highly expressed in leaf, shoot, flower and achene. In F.×ananassa, FxaYABBY1, FxaYABBY2, FxaYABBY5 and FxaYABBY6 were expressed in achene, and all FxaYABBY were poorly or not expressed in receptacle. Additionally, under the abiotic stresses of low temperature, high salt and drought, the expression of FvYABBY1, FvYABBY3, FvYABBY4 and FvYABBY6 were down-regulated, FvYABBY5 was up-regulated, and FvYABBY2 was up-regulated and then down-regulated. In tobacco leaf cells, the subcellular localization of FvYABBY proteins were in the nucleus. These results provides a foundation for the functional researches of YABBY gene in strawberry.

hirsutum

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

The Identification, Characterization, and Functional Analysis of the Sugar Transporter Gene Family of the Rice False Smut Pathogen, Villosiclava virens

False smut, caused by Villosiclava virens, is becoming increasingly serious in modern rice production systems, leading to yield losses and quality declines. Successful infection requires efficient acquisition of sucrose, abundant in rice panicles, as well as other sugars. Sugar transporters (STPs) may play an important role in this process. STPs belong to a major facilitator superfamily, which consists of large multigenic families necessary to partition sugars between fungal pathogens and their hosts. This study identified and characterized the STP family of V. viren, and further analyzed their gene functions to uncover their roles in interactions with rice. Through genome-wide and systematic bioinformatics analyses, 35 STPs were identified from V.virens and named from VvSTP1 to VvSTP35. Transmembrane domains, gene structures, and conserved motifs of VvSTPs have been identified and characterized through the bioinformatic analysis. In addition, a phylogenetic analysis revealed relationship between VvSTPs and STPs from the other three reference fungi. According to a qRT-PCR and RNA-sequencing analysis, VvSTP expression responded differently to different sole carbon sources and H2O2 treatments, and changed during the pathogenic process, suggesting that these proteins are involved in interactions with rice and potentially functional in pathogenesis. In total, 12 representative VvSTPs were knocked out through genetic recombination in order to analyze their roles in pathogenicity of V. virens. The knock-out mutants of VvSTPs showed little difference in mycelia growth and conidiation, indicating a single gene in this family cannot influence vegetative growth of V. virens. It is clear, however, that these mutants result in a change in infection efficiency in a different way, indicating that VvSTPs play an important role in the pathogenicity of virens. This study is expected to contribute to a better understanding of how host-derived sugars contribute to V. virens pathogenicity.

Genome-wide identification of the glutamate receptor-like gene family in Vanilla planifolia and their response to Fusarium oxysporum infection

Glutamate receptor-like genes (GLRs) are essential for plant growth and development and for coping with environmental (biological and non-biological) stresses. In this study, 13 GLR members were identified in the Vanilla planifolia genome and attributed to two subgroups (Clade I and Clade III) based on their physical relationships. Cis-acting element analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations indicated the GLR gene regulation's complexity and their functional diversity. Expression analysis revealed a relatively higher and more general expression pattern of Clade III members compared to the Clade I subgroup in tissues. Most GLRs showed significant differences in expression during Fusarium oxysporum infection. This suggested that GLRs play a critical role in the response of V. planifolia to pathogenic infection. These results provide helpful information for further functional research and crop improvement of VpGLRs.

The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)

Micro RNAs (miRNAs) are short non-coding RNAs that act as post-transcriptional gene expression regulators. Genes regulated in vertebrates include those affecting growth and development or stress and immune response. Pikeperch (Sander lucioperca) is a species that is increasingly being considered for farming in recirculation aquaculture systems. We characterized the pikeperch miRNA repertoire to increase the knowledge of the genomic mechanisms affecting performance and health traits by applying small RNA sequencing to different developmental stages and organs. There were 234 conserved and 8 novel miRNA genes belonging to 104 families. A total of 375 unique mature miRNAs were processed from these genes. Many mature miRNAs showed high relative abundances or were significantly more expressed at early developmental stages, like the miR-10 and miR-430 family, let-7, the miRNA clusters 106-25-93, and 17-19-92. Several miRNAs associated with immune responses (e.g., slu-mir-731-5p, slu-mir-2188-5p, and slu-mir-8159-5p) were enriched in the spleen. The mature miRNAs slu-mir-203a-3p and slu-mir-205-5p were enriched in gills. These miRNAs are similarly abundant in many vertebrates, indicating that they have shared regulatory functions. There was also a significantly increased expression of the disease-associated miR-462/miR-731 cluster in response to hypoxia stress. This first pikeperch miRNAome reference resource paves the way for future functional studies to identify miRNA-associated variations that can be utilized in marker-assisted breeding programs.

Genome-Wide Identification and Expression Analysis of WNK Kinase Gene Family in Acorus

WNK (With No Lysine) kinases are members of serine/threonine protein kinase family, which lack conserved a catalytic lysine (K) residue in protein kinase subdomain II and this residue is replaced by either asparagine, serine, or glycine residues. They are involved in various physiological regulations of flowering time, circadian rhythms, and abiotic stresses in plants. In this study, we identified the WNK gene family in two species of Acorus, and analyzed their phylogenetic relationship, physiochemical properties, subcellular localization, collinearity, and cis-elements. The results showed twenty-two WNKs in two Acorus (seven in Ac. gramineus and fifteen in Ac. calamus) have been identified and clustered into five main clades phylogenetically. Gene structure analysis showed all WNKs possessed essential STKc_WNK or PKc_like superfamily domains, and the gene structures and conserved motifs of the same clade were similar. All the WNKs harbored a large number of light response elements, plant hormone signaling elements, and stress resistance elements. Through a collinearity analysis, two and fourteen segmental duplicated gene pairs were identified in the Ac. gramineus and Ac. calamus, respectively. Moreover, we observed tissue-specificity of WNKs in Acorus using transcriptomic data, and their expressions in response to salt stress and cold stress were analyzed by qRT-PCR. The results showed WNKs are involved in the regulation of abiotic stresses. There were significant differences in the expression levels of most of the WNKs in the leaves and roots of Acorus under salt stress and cold stress, among which two members in Ac. gramineus (AgWNK3 and AgWNK4) and two members in Ac. calamus (AcWNK8 and AcWNK12) were most sensitive to stress. In summary, this paper will significantly contribute to the understanding of WNKs in monocots and thus provide a set up for functional genomics studies of WNK protein kinases.

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

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

Genome-Wide Identification and Expression Pattern Profiling of the Aquaporin Gene Family in Papaya (Carica papaya L.)

Aquaporins (AQPs) are mainly responsible for the transportation of water and other small molecules such as CO2 and H2O2, and they perform diverse functions in plant growth, in development, and under stress conditions. They are also active participants in cell signal transduction in plants. However, little is known about AQP diversity, biological functions, and protein characteristics in papaya. To better understand the structure and function of CpAQPs in papaya, a total of 29 CpAQPs were identified and classified into five subfamilies. Analysis of gene structure and conserved motifs revealed that CpAQPs exhibited a degree of conservation, with some differentiation among subfamilies. The predicted interaction network showed that the PIP subfamily had the strongest protein interactions within the subfamily, while the SIP subfamily showed extensive interaction with members of the PIP, TIP, NIP, and XIP subfamilies. Furthermore, the analysis of CpAQPs' promoters revealed a large number of cis-elements participating in light, hormone, and stress responses. CpAQPs exhibited different expression patterns in various tissues and under different stress conditions. Collectively, these results provided a foundation for further functional investigations of CpAQPs in ripening, as well as leaf, flower, fruit, and seed development. They also shed light on the potential roles of CpAQP genes in response to environmental factors, offering valuable insights into their biological functions in papaya.

Comparative Transcriptomic Assessment of Chemosensory Genes in Adult and Larval Olfactory Organs of Cnaphalocrocis medinalis

The rice leaf folder, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae), is a notorious pest of rice in Asia. The larvae and adults of C. medinalis utilize specialized chemosensory systems to adapt to different environmental odors and physiological behaviors. However, the differences in chemosensory genes between the olfactory organs of these two different developmental stages remain unclear. Here, we conducted a transcriptome analysis of larvae heads, male antennae, and female antennae in C. medinalis and identified 131 putative chemosensory genes, including 32 OBPs (8 novel OBPs), 23 CSPs (2 novel CSPs), 55 ORs (17 novel ORs), 19 IRs (5 novel IRs) and 2 SNMPs. Comparisons between larvae and adults of C. medinalis by transcriptome and RT-qPCR analysis revealed that the number and expression of chemosensory genes in larval heads were less than that of adult antennae. Only 17 chemosensory genes (7 OBPs and 10 CSPs) were specifically or preferentially expressed in the larval heads, while a total of 101 chemosensory genes (21 OBPs, 9 CSPs, 51 ORs, 18 IRs, and 2 SNMPs) were specifically or preferentially expressed in adult antennae. Our study found differences in chemosensory gene expression between larvae and adults, suggesting their specialized functions at different developmental stages of C. medinalis. These results provide a theoretical basis for screening chemosensory genes as potential molecular targets and developing novel management strategies to control C. medinalis.

Genome-Wide Analysis of Family-1 UDP-Glycosyltransferases in Potato (Solanum tuberosum L.): Identification, Phylogenetic Analysis and Determination of Response to Osmotic Stress

Family-1 UDP-glycosyltransferases (UGTs) are the most common and functional glycosyltransferases in the plant world. UGT is closely related to plant growth and the response to abiotic stress. However, despite systematic research, our understanding of potato UGT genes is still unclear. In this study, we identified 174 potato UGT proteins based on their conserved plant secondary product glycosyltransferase (PSPG) motifs. Phylogenetic analyses were used to compare these proteins with Arabidopsis UGTs and other plant UGTs, and it was found that they could be clustered into 18 distinct groups. Patterns of intron gain/loss and intron phases within potato UGTs revealed highly conserved intron insertion events. The promoter cis-elements of these 174 UGT genes were systematically investigated. The promoter regions of these UGT genes are known to contain various classes of cis-acting compounds. These include elements that are light-responsive, phytohormone-responsive, and stress-responsive. Transcriptome data analysis established that 25, 10, 6, and 4 of these 174 UGT genes were specifically expressed in leaves, roots, stolons, and young tubers, respectively. The mannitol-treated transcriptomic data showed thirty-eight UGT genes were significantly upregulated. The quantitative real-time PCR results showed that the four genes were all responsive to osmotic stress under a 10% PEG6000 treatment. The results of our study provide a basis for clarifying the molecular mechanism of potato osmotic stress resistance and better understanding its function in the future.

Evaluation of Reference Genes for Normalizing RT-qPCR and Analysis of the Expression Patterns of WRKY1 Transcription Factor and Rhynchophylline Biosynthesis-Related Genes in Uncaria rhynchophylla

Uncaria rhynchophylla (Miq.) Miq. ex Havil, a traditional medicinal herb, is enriched with several pharmacologically active terpenoid indole alkaloids (TIAs). At present, no method has been reported that can comprehensively select and evaluate the appropriate reference genes for gene expression analysis, especially the transcription factors and key enzyme genes involved in the biosynthesis pathway of TIAs in U. rhynchophylla. Reverse transcription quantitative PCR (RT-qPCR) is currently the most common method for detecting gene expression levels due to its high sensitivity, specificity, reproducibility, and ease of use. However, this methodology is dependent on selecting an optimal reference gene to accurately normalize the RT-qPCR results. Ten candidate reference genes, which are homologues of genes used in other plant species and are common reference genes, were used to evaluate the expression stability under three stress-related experimental treatments (methyl jasmonate, ethylene, and low temperature) using multiple stability analysis methodologies. The results showed that, among the candidate reference genes, S-adenosylmethionine decarboxylase (SAM) exhibited a higher expression stability under the experimental conditions tested. Using SAM as a reference gene, the expression profiles of 14 genes for key TIA enzymes and a WRKY1 transcription factor were examined under three experimental stress treatments that affect the accumulation of TIAs in U. rhynchophylla. The expression pattern of WRKY1 was similar to that of tryptophan decarboxylase (TDC) under ETH treatment. This research is the first to report the stability of reference genes in U. rhynchophylla and provides an important foundation for future gene expression analyses in U. rhynchophylla. The RT-qPCR results indicate that the expression of WRKY1 is similar to that of TDC under ETH treatment. It may coordinate the expression of TDC, providing a possible method to enhance alkaloid production in the future through synthetic biology.

The Basic/Helix-Loop-Helix Transcription Factor Family Gene RcbHLH112 Is a Susceptibility Gene in Gray Mould Resistance of Rose (Rosa Chinensis)

The basic/helix-loop-helix (bHLH) family is a major family of transcription factors in plants. Although it has been reported that bHLH plays a defensive role against pathogen infection in plants, there is no comprehensive study on the bHLH-related defence response in rose (Rosa sp.). In this study, a genome-wide analysis of bHLH family genes (RcbHLHs) in rose was carried out, including their phylogenetic relationships, gene structure, chromosome localization and collinearity analysis. Via phylogenetic analysis, a total of 121 RcbHLH genes in the rose genome were divided into 21 sub-groups. These RcbHLHs are unevenly distributed in all 7 chromosomes of rose. The occurrence of gene duplication events indicates that whole-genome duplication and segmental duplication may play a key role in gene duplication. Ratios of non-synonymous to synonymous mutation frequency (Ka/Ks) analysis showed that the replicated RcbHLH genes mainly underwent purification selection, and their functional differentiation was limited. Gene expression analysis showed that 46 RcbHLHs were differentially expressed in rose petals upon B. cinerea infection. It is speculated that these RcbHLHs are candidate genes that regulate the response of rose plants to B. cinerea infection. Virus-induced gene silencing (VIGS) confirmed that RcbHLH112 in rose is a susceptibility factor for infection with B. cinerea. This study provides useful information for further study of the functions of the rose bHLH gene family.

Genome-Wide Identification and Characterization of CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Foxtail Millet (Setaria italica L.)

The CLAVATA3/EMBRYO-SURROUNDING REGION (CLE) genes encode signaling peptides that play important roles in various developmental and physiological processes. However, the systematic identification and characterization of CLE genes in foxtail millet (Setaria italica L.) remain limited. In this study, we identified and characterized 41 SiCLE genes in the foxtail millet genome. These genes were distributed across nine chromosomes and classified into four groups, with five pairs resulting from gene duplication events. SiCLE genes within the same phylogenetic group shared similar gene structure and motif patterns, while 34 genes were found to be single-exon genes. All SiCLE peptides harbored the conserved C-terminal CLE domain, with highly conserved positions in the CLE core sequences shared among foxtail millet, Arabidopsis, rice, and maize. The SiCLE genes contained various cis-elements, including five plant hormone-responsive elements. Notably, 34 SiCLE genes possessed more than three types of phytohormone-responsive elements on their promoters. Comparative analysis revealed higher collinearity between CLE genes in maize and foxtail millet, which may be because they are both C4 plants. Tissue-specific expression patterns were observed, with genes within the same group exhibiting similar and specific expression profiles. SiCLE32 and SiCLE41, classified in Group D, displayed relatively high expression levels in all tissues except panicles. Most SiCLE genes exhibited low expression levels in young panicles, while SiCLE6, SiCLE24, SiCLE25, and SiCLE34 showed higher expression in young panicles, with SiCLE24 down-regulated during later panicle development. Greater numbers of SiCLE genes exhibited higher expression in roots, with SiCLE7, SiCLE22, and SiCLE36 showing the highest levels and SiCLE36 significantly down-regulated after abscisic acid (ABA) treatment. Following treatments with ABA, 6-benzylaminopurine (6-BA), and gibberellic acid 3 (GA3), most SiCLE genes displayed down-regulation followed by subsequent recovery, while jasmonic acid (JA) and indole-3-acetic acid (IAA) treatments led to upregulation at 30 min in leaves. Moreover, identical hormone treatments elicited different expression patterns of the same genes in leaves and stems. This comprehensive study enhances our understanding of the SiCLE gene family and provides a foundation for further investigations into the functions and evolution of SiCLE genes in foxtail millet.

Transcriptome-Wide Identification of TCP Transcription Factor Family Members in Pinus massoniana and Their Expression in Regulation of Development and in Response to Stress

Pinus massoniana is an important coniferous tree species for barren mountain afforestation with enormous ecological and economic significance. It has strong adaptability to the environment. TEOSINTE BRANCHED 1/CYCLOIDEA/PCF (TCP) transcription factors (TFs) play crucial roles in plant stress response, hormone signal transduction, and development processes. At present, TCP TFs have been widely studied in multiple plant species, but research in P. massoniana has not been carried out. In this study, 13 PmTCP TFs were identified from the transcriptomes of P. massoniana. The phylogenetic results revealed that these PmTCP members were divided into two categories: Class I and Class II. Each PmTCP TF contained a conserved TCP domain, and the conserved motif types and numbers were similar in the same subgroup. According to the transcriptional profiling analysis under drought stress conditions, it was found that seven PmTCP genes responded to drought treatment to varying degrees. The qRT-PCR results showed that the majority of PmTCP genes were significantly expressed in the needles and may play a role in the developmental stage. Meanwhile, the PmTCPs could respond to several stresses and hormone treatments at different levels, which may be important for stress resistance. In addition, PmTCP7 and PmTCP12 were nuclear localization proteins, and PmTCP7 was a transcriptional suppressor. These results will help to explore the regulatory factors related to the growth and development of P. massoniana, enhance its stress resistance, and lay the foundation for further exploration of the physiological effects on PmTCPs.

Genome-wide investigation and expression analysis of OSCA gene family in response to abiotic stress in alfalfa

Alfalfa is an excellent leguminous forage crop that is widely cultivated worldwide, but its yield and quality are often affected by drought and soil salinization. Hyperosmolality-gated calcium-permeable channel (OSCA) proteins are hyperosmotic calcium ion (Ca2+) receptors that play an essential role in regulating plant growth, development, and abiotic stress responses. However, no systematic analysis of the OSCA gene family has been conducted in alfalfa. In this study, a total of 14 OSCA genes were identified from the alfalfa genome and classified into three groups based on their sequence composition and phylogenetic relationships. Gene structure, conserved motifs and functional domain prediction showed that all MsOSCA genes had the same functional domain DUF221. Cis-acting element analysis showed that MsOSCA genes had many cis-regulatory elements in response to abiotic or biotic stresses and hormones. Tissue expression pattern analysis demonstrated that the MsOSCA genes had tissue-specific expression; for example, MsOSCA12 was only expressed in roots and leaves but not in stem and petiole tissues. Furthermore, RT-qPCR results indicated that the expression of MsOSCA genes was induced by abiotic stress (drought and salt) and hormones (JA, SA, and ABA). In particular, the expression levels of MsOSCA3, MsOSCA5, MsOSCA12 and MsOSCA13 were significantly increased under drought and salt stress, and MsOSCA7, MsOSCA10, MsOSCA12 and MsOSCA13 genes exhibited significant upregulation under plant hormone treatments, indicating that these genes play a positive role in drought, salt and hormone responses. Subcellular localization results showed that the MsOSCA3 protein was localized on the plasma membrane. This study provides a basis for understanding the biological information and further functional analysis of the MsOSCA gene family and provides candidate genes for stress resistance breeding in alfalfa.

Identification of the Maize LEA Gene Family and Its Relationship with Kernel Dehydration

Maize, the most widely planted and highest yielding of the three major crops in the world, requires the development and breeding of new varieties to accommodate the shift towards mechanized harvesting. However, the moisture content of kernels during harvest poses a significant challenge to mechanized harvesting, leading to seed breakage and increased storage costs. Previous studies highlighted the importance of LEA (Late Embryogenesis Abundant) members in regulating kernel dehydration. In this study, we aimed to gain a better understanding of the relationship between the LEA family and grain dehydration in maize. Through expression pattern analysis of maize, we identified 52 LEA genes (ZmLEAs) distributed across 10 chromosomes, organized into seven subgroups based on phylogenetic analysis, gene structure, and conserved motifs. Evolutionary and selective pressure analysis revealed that the amplification of ZmLEA genes primarily resulted from whole-genome or fragment replication events, with strong purifying selection effects during evolution. Furthermore, the transcriptome data of kernels of two maize inbred lines with varying dehydration rates at different developmental stages showed that 14 ZmLEA genes were expressed differentially in the two inbreds. This suggested that the ZmLEA genes might participate in regulating the kernel dehydration rate (KDR) in maize. Overall, this study enhances our understanding of the ZmLEA family and provides a foundation for further research into its role in regulating genes associated with grain dehydration in maize.

Characterization and Expression Analysis of Genes from Megalobrama amblycephala Encoding Hemoglobins with Extracellular Microbicidal Activity

Hemoglobin (Hb) usually comprises two α and two β subunits, forming a tetramer responsible for oxygen transportation and storage. Few studies have elucidated fish hemoglobin immune functions. Megalobrama amblycephala is a freshwater-cultured fish prevalent in China. We identified two M. amblycephala hemoglobin subunits and analyzed their expression patterns and antibacterial activities. The respective full-length cDNA sequences of the M. amblycephala Hb α (MaHbα) and β (MaHbβ) subunits were 588 and 603 bp, encoding 143 and 148 amino acids. MaHbα and MaHbβ were highly homologous to hemoglobins from other fish, displaying typical globin-like domains, most heme-binding sites, and tetramer interface regions highly conserved in teleosts. In phylogenetic analyses, the hemoglobin genes from M. amblycephala and other cypriniformes clustered into one branch, and those from other fishes and mammals clustered into other branches, revealing fish hemoglobin conservation. These M. amblycephala Hb subunits exhibit different expression patterns in various tissues and during development. MaHbα is mainly expressed in the blood and brain, while MaHbβ gene expression is highest in the muscle. MaHbα expression was detectable and abundant post-fertilization, with levels fluctuating during the developmental stages. MaHbβ expression began at 3 dph and gradually increased. Expression of both M. amblycephala Hb subunits was down-regulated in most examined tissues and time points post-Aeromonas hydrophila infection, which might be due to red blood cell (RBC) and hematopoietic organ damage. Synthetic MaHbα and MaHbβ peptides showed excellent antimicrobial activities, which could inhibit survival and growth in five aquatic pathogens. Two M. amblycephala hemoglobin subunits were identified, and their expression patterns and antibacterial activities were analyzed, thereby providing a basis for the understanding of evolution and functions of fish hemoglobins.

Comprehensive Functional Analysis of the bZIP Family in Bletilla striata Reveals That BsbZIP13 Could Respond to Multiple Abiotic Stresses

Basic leucine zipper (bZIP) transcription factors (TFs) are one of the largest families involved in plant physiological processes such as biotic and abiotic responses, growth, and development, etc. In this study, 66 members of the bZIP family were identified in Bletilla striata, which were divided into 10 groups based on their phylogenetic relationships with AtbZIPs. A structural analysis of BsbZIPs revealed significant intron-exon differences among BsbZIPs. A total of 63 bZIP genes were distributed across 16 chromosomes in B. striata. The tissue-specific and germination stage expression patterns of BsbZIPs were based on RNA-seq. Stress-responsive expression analysis revealed that partial BsbZIPs were highly expressed under low temperatures, wounding, oxidative stress, and GA treatments. Furthermore, subcellular localization studies indicated that BsbZIP13 was localized in the nucleus. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays suggested that BsbZIP13 could interact with multiple BsSnRK2s. The results of this study provide insightful data regarding bZIP TF as one of the stress response regulators in B. striata, while providing a theoretical basis for transgenic and functional studies of the bZIP gene family in B. striata.

Phylogenomic Analysis of Cytochrome P450 Gene Superfamily and Their Association with Flavonoids Biosynthesis in Peanut (Arachis hypogaea L.)

Cytochrome P450s (CYPs) constitute extensive enzyme superfamilies in the plants, playing pivotal roles in a multitude of biosynthetic and detoxification pathways essential for growth and development, such as the flavonoid biosynthesis pathway. However, CYPs have not yet been systematically studied in the cultivated peanuts (Arachis hypogaea L.), a globally significant cash crop. This study addresses this knowledge deficit through a comprehensive genome-wide analysis, leading to the identification of 589 AhCYP genes in peanuts. Through phylogenetic analysis, all AhCYPs were systematically classified into 9 clans, 43 gene families. The variability in the number of gene family members suggests specialization in biological functions. Intriguingly, both tandem duplication and fragment duplication events have emerged as pivotal drivers in the evolutionary expansion of the AhCYP superfamily. Ka/Ks analysis underscored the substantial influence of strong purifying selection on the evolution of AhCYPs. Furthermore, we selected 21 genes encoding 8 enzymes associated with the flavonoid pathway. The results of quantitative real-time PCR (qRT-PCR) experiments unveiled stage-specific expression patterns during the development of peanut testa, with discernible variations between pink and red testa. Importantly, we identified a direct correlation between gene expression levels and the accumulation of metabolites. These findings offer valuable insights into elucidating the comprehensive functions of AhCYPs and the underlying mechanisms governing the divergent accumulation of flavonoids in testa of different colors.

Comprehensive Analysis of GH3 Gene Family in Potato and Functional Characterization of StGH3.3 under Drought Stress

As an important hormone response gene, Gretchen Hagen 3 (GH3) maintains hormonal homeostasis by conjugating excess auxin with amino acids during plant stress-related signaling pathways. GH3 genes have been characterized in many plant species, but they are rarely reported in potato. Here, 19 StGH3 genes were isolated and characterized. Phylogenetic analysis indicated that StGH3s were divided into two categories (group I and group III). Analyses of gene structure and motif composition showed that the members of a specific StGH3 subfamily are relatively conserved. Collinearity analysis of StGH3 genes in potato and other plants laid a foundation for further exploring the evolutionary characteristics of the StGH3 genes. Promoter analysis showed that most StGH3 promoters contained hormone and abiotic stress response elements. Multiple transcriptome studies indicated that some StGH3 genes were responsive to ABA, water deficits, and salt treatments. Moreover, qRT-PCR analysis indicated that StGH3 genes could be induced by phytohormones (ABA, SA, and MeJA) and abiotic stresses (water deficit, high salt, and low temperature), although with different patterns. Furthermore, transgenic tobacco with transient overexpression of the StGH3.3 gene showed positive regulation in response to water deficits by increasing proline accumulation and reducing the leaf water loss rate. These results suggested that StGH3 genes may be involved in the response to abiotic stress through hormonal signal pathways. Overall, this study provides useful insights into the evolution and function of StGH3s and lays a foundation for further study on the molecular mechanisms of StGH3s in the regulation of potato drought resistance.

Genome-Wide Identification and Characterization of WRKY Transcription Factors in Betula platyphylla Suk. and Their Responses to Abiotic Stresses

The WRKY transcription factor (TF) family is one the largest plant-specific transcription factor families. It has been proven to play significant roles in multiple plant biological processes, especially stress response. Although many WRKY TFs have been identified in various plant species, WRKYs in white birch (Betula platyphylla Suk.) remain to be studied. Here, we identified a total of 68 BpWRKYs, which could be classified into four main groups. The basic physiochemical properties of these TFs were analyzed using bioinformatics tools, including molecular weight, isoelectric point, chromosome location, and predicted subcellular localization. Most BpWRKYs were predicted to be located in the nucleus. Synteny analysis found 17 syntenic gene pairs among BpWRKYs and 52 syntenic gene pairs between BpWRKYs and AtWRKYs. The cis-acting elements in the promoters of BpWRKYs could be enriched in multiple plant biological processes, including stress response, hormone response, growth and development, and binding sites. Tissue-specific expression analysis using qRT-PCR showed that most BpWRKYs exhibited highest expression levels in the root. After ABA, salt (NaCl), or cold treatment, different BpWRKYs showed different expression patterns at different treatment times. Furthermore, the results of the Y2H assay proved the interaction between BpWRKY17 and a cold-responsive TF, BpCBF7. By transient expression assay, BpWRKY17 and BpWRKY67 were localized in the nucleus, consistent with the previous prediction. Our study hopes to shed light for research on WRKY TFs and plant stress response.

Corrigendum: Genome-wide characterization of NAC transcription factors in Camellia sinensis and the involvement of CsNAC28 in drought tolerance

[This corrects the article DOI: 10.3389/fpls.2022.1065261.].

Genome-Wide Identification and Expression Analysis of the TCP Gene Family Related to Developmental and Abiotic Stress in Ginger

Ginger (Zingiber officinale Roscoe), a widely consumed edible and medicinal plant, possesses significant nutritional and economic value. Abiotic stresses such as drought and low temperatures can impact the growth and development of ginger. The plant-specific transcription factor Teosinte branched1/cycloidea/proliferating cell factor (TCP) has progressively been identified in various plants for its role in regulating plant growth and development as well as conferring resistance to abiotic stresses. However, limited information on the TCP family is available in ginger. In this study, we identified 20 TCP members in the ginger genome, which were randomly distributed across 9 chromosomes. Based on phylogenetic analysis, these ginger TCP were classified into two subfamilies: Class I (PCF) and Class II (CIN, CYC/TB). The classification of the identified ginger TCPs was supported by a multi-species phylogenetic tree and motif structure analysis, suggesting that the amplification of the ginger TCP gene family occurred prior to the differentiation of angiosperms. The promoter region of ginger TCP genes was found to contain numerous cis-acting elements associated with plant growth, development, and abiotic stress response. Among these elements, the stress response element, anaerobic induction, and MYB binding site play a dominant role in drought responsiveness. Additionally, expression pattern analysis revealed variations in the expression of ginger TCP gene among different tissues and in response to diverse abiotic stresses (drought, low temperature, heat, and salt). Our research offers a thorough examination of TCP members within the ginger plant. This analysis greatly contributes to the understanding of how TCP genes regulate tissue development and response to stress, opening up new avenues for further exploration in this field.

Genome-wide identification, expression analysis, and potential roles under low-temperature stress of bHLH gene family in Prunus sibirica

The basic helix-loop-helix (bHLH) family is one of the most well-known transcription factor families in plants, and it regulates growth, development, and abiotic stress responses. However, systematic analyses of the bHLH gene family in Prunus sibirica have not been reported to date. In this study, 104 PsbHLHs were identified and classified into 23 subfamilies that were unevenly distributed on eight chromosomes. Nineteen pairs of segmental replication genes and ten pairs of tandem replication genes were identified, and all duplicated gene pairs were under purifying selection. PsbHLHs of the same subfamily usually share similar motif compositions and exon-intron structures. PsbHLHs contain multiple stress-responsive elements. PsbHLHs exhibit functional diversity by interacting and coordinating with other members. Twenty PsbHLHs showed varying degrees of expression. Eleven genes up-regulated and nine genes down-regulated in -4°C. The majority of PsbHLHs were highly expressed in the roots and pistils. Transient transfection experiments demonstrated that transgenic plants with overexpressed PsbHLH42 have better cold tolerance. In conclusion, the results of this study have significant implications for future research on the involvement of bHLH genes in the development and stress responses of Prunus sibirica.

Genome-Wide Identification of TCP Gene Family in Dendrobium and Their Expression Patterns in Dendrobium chrysotoxum

The TCP gene family are plant-specific transcription factors that play important roles in plant growth and development. Dendrobium chrysotoxum, D. nobile, and D. huoshanense are orchids with a high ornamental value, but few studies have investigated the specific functions of TCPs in Dendrobium flower development. In this study, we used these three Dendrobium species to analyze TCPs, examining their physicochemical properties, phylogenetic relationships, gene structures, and expression profiles. A total of 50 TCPs were identified across three Dendrobium species; they were divided into two clades-Class-I (PCF subfamily) and Class-II (CIN and CYC/TB1 subfamilies)-based on their phylogenetic relationships. Our sequence logo analysis showed that almost all Dendrobium TCPs contain a conserved TCP domain, as well as the existence of fewer exons, and the cis-regulatory elements of the TCPs were mostly related to light response. In addition, our transcriptomic data and qRT-PCR results showed that DchTCP2 and DchTCP13 had a significant impact on lateral organs. Moreover, changes in the expression level of DchTCP4 suggested its important role in the phenotypic variation of floral organs. Therefore, this study provides a significant reference for the further exploration of TCP gene functions in the regulation of different floral organs in Dendrobium orchids.