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

Comprehensive analysis of pepper (Capsicum annuum) RAV genes family and functional identification of CaRAV1 under chilling stress

BACKGROUND

Despite its known significance in plant abiotic stress responses, the role of the RAV gene family in the response of Capsicum annuum to chilling stress remains largely unexplored.

RESULTS

In this study, we identified and characterized six members of the CaRAV gene subfamily in pepper plants through genome-wide analysis. Subsequently, the CaRAV subfamily was classified into four branches based on homology with Arabidopsis thaliana, each exhibiting relatively conserved domains within the branch. We discovered that light response elements accounted for the majority of CaRAVs, whereas low-temperature response elements were specific to the NGA gene subfamily. After pepper plants were subjected to chilling stress, qRT‒PCR analysis revealed that CaRAV1, CaRAV2 and CaNGA1 were significantly induced in response to chilling stress, indicating that CaRAVs play a role in the response to chilling stress. Using virus-induced gene silencing (VIGS) vectors, we targeted key members of the CaRAV gene family. Under normal growth conditions, the MDA content and SOD enzyme activity of the silenced plants were slightly greater than those of the control plants, and the REC activity was significantly greater than that of the control plants. The levels of MDA and electrolyte leakage were greater in the silenced plants after they were exposed to chilling stress, and the POD and CAT enzyme activities were significantly lower than those in the control, which was particularly evident under repeated chilling stress. In addition, the relative expression of CaPOD and CaCAT was greater in V2 plants upon repeated chilling stress, especially CaCAT was significantly greater in V2 plants than in the other two silenced plants, with 3.29 and 1.10 increases within 12 and 24 h. These findings suggest that CaRAV1 and CaNGA1 positively regulate the response to chilling stress.

CONCLUSIONS

Silencing of key members of the CaRAV gene family results in increased susceptibility to chilling damage and reduced antioxidant enzyme activity in plants, particularly under repeated chilling stress. This study provides valuable information for understanding the classification and putative functions of RAV transcription factors in pepper plants.

Identification and expression analysis of XIP gene family members in rice

Xylanase inhibitor proteins (XIP) are widely distributed in the plant kingdom, and also exist in rice. However, a systematic bioinformatics analysis of this gene family in rice (OsXIP) has not been conducted to date. In this study, we identified 32 members of the OsXIP gene family and analyzed their physicochemical properties, chromosomal localization, gene structure, protein structure, expression profiles, and interaction networks. Our results indicated that OsXIP genes exhibit an uneven distribution across eight rice chromosomes. These genes generally feature a low number of introns or are intronless, all family members, except for OsXIP20, contain two highly conserved motifs, namely Motif 8 and Motif 9. In addition, it is worth noting that the promoter regions of OsXIP gene family members feature a widespread presence of abscisic acid response elements (ABRE) and gibberellin response elements (GARE-motif and TATC-box). Quantitative Real-time PCR (qRT-PCR) analysis unveiled that the expression of OsXIP genes exhibited higher levels in leaves and roots, with considerable variation in the expression of each gene in these tissues both prior to and following treatments with abscisic acid (ABA) and gibberellin (GA3). Protein interaction studies and microRNA (miRNA) target prediction showed that OsXIP engages with key elements within the hormone-responsive and drought signaling pathways. The qRT-PCR suggested osa-miR2927 as a potential key regulator in the rice responding to drought stress, functioning as tissue-specific and temporally regulation. This study provides a theoretical foundation for further analysis of the functions within the OsXIP gene family.

Genome-wide characterization of LEA gene family reveals a positive role of BnaA.LEA6.a in freezing tolerance in rapeseed (Brassica napus L.)

BACKGROUND

Freezing stress is one of the major abiotic stresses that causes extensive damage to plants. LEA (Late embryogenesis abundant) proteins play a crucial role in plant growth, development, and abiotic stress. However, there is limited research on the function of LEA genes in low-temperature stress in Brassica napus (rapeseed).

RESULTS

Total 306 potential LEA genes were identified in B. rapa (79), B. oleracea (79) and B. napus (148) and divided into eight subgroups. LEA genes of the same subgroup had similar gene structures and predicted subcellular locations. Cis-regulatory elements analysis showed that the promoters of BnaLEA genes rich in cis-regulatory elements related to various abiotic stresses. Additionally, RNA-seq and real-time PCR results indicated that the majority of BnaLEA family members were highly expressed in senescent tissues of rapeseed, especially during late stages of seed maturation, and most BnaLEA genes can be induced by salt and osmotic stress. Interestingly, the BnaA.LEA6.a and BnaC.LEA6.a genes were highly expressed across different vegetative and reproductive organs during different development stages, and showed strong responses to salt, osmotic, and cold stress, particularly freezing stress. Further analysis showed that overexpression of BnaA.LEA6.a increased the freezing tolerance in rapeseed, as evidenced by lower relative electrical leakage and higher survival rates compared to the wild-type (WT) under freezing treatment.

CONCLUSION

This study is of great significance for understanding the functions of BnaLEA genes in freezing tolerance in rapeseed and offers an ideal candidate gene (BnaA.LEA6.a) for molecular breeding of freezing-tolerant rapeseed cultivars.

Genome-Wide Identification and Expression Analysis of Bx Involved in Benzoxazinoids Biosynthesis Revealed the Roles of DIMBOA during Early Somatic Embryogenesis in Dimocarpus longan Lour

Benzoxazinoids (BXs) are tryptophan-derived indole metabolites and play a role in various physiological processes, such as auxin metabolism. Auxin is essential in the process of somatic embryogenesis (SE) in plants. In this study, we used bioinformatics, transcriptome data, exogenous treatment experiments, and qPCR analysis to study the evolutionary pattern of Bx genes in green plants, the regulatory mechanism of DlBx genes during early SE, and the effect of 2,4-dihydroxy-7-methoxy-1,4-benzoxazine-3-one (DIMBOA) on the early SE in Dimocarpus longan Lour. The results showed that 27 putative DlBxs were identified in the longan genome; the Bx genes evolved independently in monocots and dicots, and the main way of gene duplication for the DlBx was tandem duplication (TD) and the DlBx were strongly constrained by purification selection during evolution. The transcriptome data indicated varying expression levels of DlBx during longan early SE, and most DlBxs responded to light, temperature, drought stress, and 2,4-dichlorophenoxyacetic acid (2,4-D) treatment; qRT-PCR results showed DlBx1, DlBx6g and DlBx6h were responsive to auxin, and treatment with 0.1mg/L DIMBOA for 9 days significantly upregulated the expression levels of DlBx1, DlBx3g, DlBx6c, DlBx6f, DlB6h, DlBx7d, DlBx8, and DlBx9b. The correlation analysis showed a significantly negative correlation between the expression level of DlBx1 and the endogenous IAA contents; DIMBOA significantly promoted the early SE and significantly changed the endogenous IAA content, and the IAA content increased significantly at the 9th day and decreased significantly at the 13th day. Therefore, the results suggested that DIMBOA indirectly promote the early SE by changing the endogenous IAA content via affecting the expression level of DlBx1 and hydrogen peroxide (H2O2) content in longan.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C. japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, Ks, and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

Comparative transcriptional analysis of Persea americana MYB, WRKY and AP2/ERF transcription factors following Phytophthora cinnamomi infection

Plant cells undergo extensive transcriptional reprogramming following pathogen infection, with these reprogramming patterns becoming more complex when pathogens, such as hemibiotrophs, exhibit different lifestyles. These transcriptional changes are often orchestrated by MYB, WRKY and AP2/ERF transcription factors (TFs), which modulate both growth and defence-related gene expression. Transcriptional analysis of defence-related genes in avocado (Persea americana) infected with Phytophthora cinnamomi indicated differential immune response activation when comparing a partially resistant and susceptible rootstock. This study identified 226 MYB, 82 WRKY, and 174 AP2/ERF TF-encoding genes in avocado, using a genome-wide approach. Phylogenetic analysis revealed substantial sequence conservation within TF groups underscoring their functional significance. RNA-sequencing analysis in a partially resistant and susceptible avocado rootstock infected with P. cinnamomi was indicative of an immune response switch occurring in either rootstock after 24 and 6 h post-inoculation, respectively. Different clusters of co-expressed TF genes were observed at these times, suggesting the activation of necrotroph-related immune responses at varying intervals between the two rootstocks. This study aids our understanding of avocado immune response activation following P. cinnamomi infection, and the role of the TFs therein, elucidating the transcriptional reprogramming disparities between partially resistant and susceptible rootstocks.

Molecular characterization and evolutionary relationships of DOFs in four cherry species and functional analysis in sweet cherry

The DOF (DNA binding with one finger) has multiple functions in plants. However, it has received little attention in the research field of cherries. In this study, the evolutionary relationship and molecular characterization of DOF in four cherry species were analyzed, revealing its expression pattern in sweet cherry. There are 23 members in Prunus avium cv. 'Tieton', 88 in Prunus cerasus, 53 in Cerasus × yedoensis, and 27 in Cerasus serrulata. Most of these genes are intron-less or non-intron, with a conserved C2-C2 domain. Due to heterozygosity and chromosomal ploidy, whole-genome duplication (WGD) events occur to varying degrees, and DOF genes are contracted during evolution. Furthermore, these genes are affected by purifying selection pressure. Under low-temperature treatment, the expression of PavDOF2 and PavDOF18 were significantly up-regulated, while that of PavDOF16 is significantly down-regulated. The expression of PavDOF9, PavDOF12, PavDOF14, PavDOF16, PavDOF17, PavDOF18, and PavDOF19 exhibits an increasing trend during flower development and varies during sweet cherry fruit development. PavDOF1, PavDOF8, PavDOF9, and PavDOF15 are localized in the nucleus but is not transcriptionally active. The findings systemically demonstrate the molecular characteristics of DOF in different cherry varieties, providing a basis for further research on the functions of these genes.

CiAP2/ERF65 and CiAP2/ERF106, a pair of homologous genes in pecan (Carya illinoensis), regulate plant responses during submergence in transgenic Arabidopsis thaliana

When plants are entirely submerged, photosynthesis and respiration are severely restricted, affecting plant growth and potentially even causing plant death. The AP2/ERF superfamily has been widely reported to play a vital role in plant growth, development and resistance to biotic and abiotic stresses. However, no relevant studies exist on flooding stress in pecan. In this investigation, we observed that CiAP2/ERF65 positively modulated the hypoxia response during submergence, whereas CiAP2/ERF106 was sensitive to submergence. The levels of physiological and biochemical indicators, such as POD, CAT and among others, in CiAP2/ERF65-OE lines were significantly higher than those in wild-type Arabidopsis thaliana, indicating that the antioxidant capacity of CiAP2/ERF65-OE lines was enhanced under submergence. The RNA-seq results revealed that the maintenance of the expression levels of the antenna protein gene, different signaling pathways for regulation, as well as the storage and consumption of ATP, might account for the opposite phenotypes of CiAP2/ERF65 and CiAP2/ERF106. Furthermore, the expression of some stress-related genes was altered during submergence and reoxygenation. Overall, these findings enhance our understanding of submergence stress in pecan, providing important candidate genes for the molecular design and breeding of hypoxia resistant in plants.

Genome-wide identification of the PFK gene family and their expression analysis in Quercus rubra

The glycolytic pathway involves phosphofructokinase (PFK), a rate-limiting enzyme that catalyzes the phosphorylation of fructose-6-phosphate. In plants, the two PFK members are ATP-dependent phosphofructokinase (PFK) and pyrophosphate-fructose-6-phosphate phosphotransferase (PFP). However, the functions of the PFK family members in Quercus rubra are not well understood. The purpose of this study was to investigate the genome-wide distribution of the PFK family members and their roles in Q. rubra by performing a systematic study of the phylogenetic relationships, molecular characteristics, motifs, chromosomal and subcellular locations, and cis-elements of QrPFKs. We identified 14 QrPFK genes in the genome of Q. rubra, followed by examining their expression in different tissues, including the roots, stems, and leaves. The phylogenetic tree divided the 14 QrPFK genes into two groups: 11 belonging to PFK and three belonging to PFP. The expression profiles of all 14 proteins were relatively the same in leaves but differed between stems and roots. Four genes (Qurub.02G189400.1, Qurub.02G189400.2, Qurub.09G134300.1, and Qurub.09G134300.2) were expressed at very low levels in both stems and roots, while two (Qurub.05G235500.1 and Qurub.05G235500.1) were expressed at low levels and the others showed relatively high expression in all tissues.

Introns and Their Therapeutic Applications in Biomedical Researches

Context

Although for a long time, it was thought that intervening sequences (introns) were junk DNA without any function, their critical roles and the underlying molecular mechanisms in genome regulation have only recently come to light. Introns not only carry information for splicing, but they also play many supportive roles in gene regulation at different levels. They are supposed to function as useful tools in various biological processes, particularly in the diagnosis and treatment of diseases. Introns can contribute to numerous biological processes, including gene silencing, gene imprinting, transcription, mRNA metabolism, mRNA nuclear export, mRNA localization, mRNA surveillance, RNA editing, NMD, translation, protein stability, ribosome biogenesis, cell growth, embryonic development, apoptosis, molecular evolution, genome expansion, and proteome diversity through various mechanisms.

Evidence Acquisition

In order to fulfill the objectives of this study, the following databases were searched: Medline, Scopus, Web of Science, EBSCO, Open Access Journals, and Google Scholar. Only articles published in English were included.

Results & Conclusions

The intervening sequences of eukaryotic genes have critical functions in genome regulation, as well as in molecular evolution. Here, we summarize recent advances in our understanding of how introns influence genome regulation, as well as their effects on molecular evolution. Moreover, therapeutic strategies based on intron sequences are discussed. According to the obtained results, a thorough understanding of intron functional mechanisms could lead to new opportunities in disease diagnosis and therapies, as well as in biotechnology applications.

Hemibiotrophic Phytophthora infestans Modulates the Expression of SWEET Genes in Potato (Solanum tuberosum L.)

Sugar Efflux transporters (SWEET) are involved in diverse biological processes of plants. Pathogens have exploited them for nutritional gain and subsequently promote disease progression. Recent studies have implied the involvement of potato SWEET genes in the most devastating late blight disease caused by Phytophthora infestans. Here, we identified and designated 37 putative SWEET genes as StSWEET in potato. We performed detailed in silico analysis, including gene structure, conserved domains, and phylogenetic relationship. Publicly available RNA-seq data was harnessed to retrieve the expression profiles of SWEET genes. The late blight-responsive SWEET genes were identified from the RNA-seq data and then validated using quantitative real-time PCR. The SWEET gene expression was studied along with the biotrophic (SNE1) and necrotrophic (PiNPP1) marker genes of P. infestans. Furthermore, we explored the co-localization of P. infestans resistance loci and SWEET genes. The results indicated that nine transporter genes were responsive to the P. infestans in potato. Among these, six transporters, namely StSWEET10, 12, 18, 27, 29, and 31, showed increased expression after P. infestans inoculation. Interestingly, the observed expression levels aligned with the life cycle of P. infestans, wherein expression of these genes remained upregulated during the biotrophic phase and decreased later on. In contrast, StSWEET13, 14, and 32 didn't show upregulation in inoculated samples suggesting non-targeting by pathogens. This study underscores these transporters as prime P. infestans targets in potato late blight, pivotal in disease progression, and potential candidates for engineering blight-resistant potato genotypes.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Genome-Wide Comparison and Functional Characterization of HMGR Gene Family Associated with Shikonin Biosynthesis in Lithospermum erythrorhizon

3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), as the rate-limiting enzyme in the mevalonate pathway, is essential for the biosynthesis of shikonin in Lithospermum erythrorhizon. However, in the absence of sufficient data, the principles of a genome-wide in-depth evolutionary exploration of HMGR family members in plants, as well as key members related to shikonin biosynthesis, remain unidentified. In this study, 124 HMGRs were identified and characterized from 36 representative plants, including L. erythrorhizon. Vascular plants were found to have more HMGR family genes than nonvascular plants. The phylogenetic tree revealed that during lineage and species diversification, the HMGRs evolved independently and intronless LerHMGRs emerged from multi-intron HMGR in land plants. Among them, Pinus tabuliformis and L. erythrorhizon had the most HMGR gene duplications, with 11 LerHMGRs most likely expanded through WGD/segmental and tandem duplications. In seedling roots and M9 cultured cells/hairy roots, where shikonin biosynthesis occurs, LerHMGR1 and LerHMGR2 were expressed significantly more than other genes. The enzymatic activities of LerHMGR1 and LerHMGR2 further supported their roles in catalyzing the conversion of HMG-CoA to mevalonate. Our findings provide insight into the molecular evolutionary properties and function of the HMGR family in plants and a basis for the genetic improvement of efficiently produced secondary metabolites in L. erythrorhizon.

Genome-Wide Analyses of SlFWL Family Genes and Their Expression Profiles under Cold, Heat, Salt and Drought Stress in Tomato

PLAC8 is a cysteine-rich protein that serves as a central mediator of tumor evolution in mammals. PLAC8 motif-containing proteins widely distribute in fungi, algae, higher plants and animals that have been described to be implicated in fruit size, cell number and the transport of heavy metals such as cadmium or zinc. In tomatoes, FW2.2 is a PLAC8 motif-containing gene that negatively controls fruit size by regulating cell division and expansion in the carpel ovary during fruit development. However, despite FW2.2, other FWL (FW2.2-Like) genes in tomatoes have not been investigated. In this study, we identified the 21 SlFWL genes, including FW2.2, examined their expression profiles under various abiotic adversity-related conditions. The SlFWL gene structures and motif compositions are conserved, indicating that tomato SlFWL genes may have similar roles. Cis-acting element analysis revealed that the SlFWL genes may participate in light and abiotic stress responses, and they also interacted with a variety of phytohormone-responsive proteins and plant development elements. Phylogenetic analyses were performed on five additional plant species, including Arabidopsis, pepper, soybean, rice and maize, these genes were classified into five subfamilies. Based on the results of collinearity analyses, the SlFWL genes have a tighter homologous evolutionary relationship with soybean, and these orthologous FWL gene pairs might have the common ancestor. Expression profiling of SlFWL genes show that they were all responsive to abiotic stresses, each subgroup of genes exhibited a different expression trend. Our findings provide a strong foundation for investigating the function and abiotic stress responses of the SlFWL family genes.

DREB1 and DREB2 Genes in Garlic (Allium sativum L.): Genome-Wide Identification, Characterization, and Stress Response

Dehydration-responsive element-binding (DREB) transcription factors (TFs) of the A1 and A2 subfamilies involved in plant stress responses have not yet been reported in Allium species. In this study, we used bioinformatics and comparative transcriptomics to identify and characterize DREB A1 and A2 genes redundant in garlic (Allium sativum L.) and analyze their expression in A. sativum cultivars differing in the sensitivity to cold and Fusarium infection. Eight A1 (AsaDREB1.1-1.8) and eight A2 (AsaDREB2.1-2.8) genes were identified. AsaDREB1.1-1.8 genes located in tandem on chromosome 1 had similar expression patterns, suggesting functional redundancy. AsaDREB2.1-2.8 were scattered on different chromosomes and had organ- and genotype-specific expressions. AsaDREB1 and AsaDREB2 promoters contained 7 and 9 hormone- and stress-responsive cis-regulatory elements, respectively, and 13 sites associated with TF binding and plant development. In both Fusarium-resistant and -sensitive cultivars, fungal infection upregulated the AsaDREB1.1-1.5, 1.8, 2.2, 2.6, and 2.8 genes and downregulated AsaDREB2.5, but the magnitude of response depended on the infection susceptibility of the cultivar. Cold exposure strongly upregulated the AsaDREB1 genes, but downregulated most AsaDREB2 genes. Our results provide the foundation for further functional analysis of the DREB TFs in Allium crops and could contribute to the breeding of stress-tolerant varieties.

Intra-Varietal Diversity and Its Contribution to Wheat Evolution, Domestication, and Improvement in Wheat

Crop genetic diversity is essential for adaptation and productivity in agriculture. A previous study revealed that poor allele diversity in wheat commercial cultivars is a major barrier to its further improvement. Homologs within a variety, including paralogs and orthologs in polyploid, account for a large part of the total genes of a species. Homolog diversity, intra-varietal diversity (IVD), and their functions have not been elucidated. Common wheat, an important food crop, is a hexaploid species with three subgenomes. This study analyzed the sequence, expression, and functional diversity of homologous genes in common wheat based on high-quality reference genomes of two representative varieties, a modern commercial variety Aikang 58 (AK58) and a landrace Chinese Spring (CS). A total of 85,908 homologous genes, accounting for 71.9% of all wheat genes, including inparalogs (IPs), outparalogs (OPs), and single-copy orthologs (SORs), were identified, suggesting that homologs are an important part of the wheat genome. The levels of sequence, expression, and functional variation in OPs and SORs were higher than that of IPs, which indicates that polyploids have more homologous diversity than diploids. Expansion genes, a specific type of OPs, made a great contribution to crop evolution and adaptation and endowed crop with special characteristics. Almost all agronomically important genes were from OPs and SORs, demonstrating their essential functions for polyploid evolution, domestication, and improvement. Our results suggest that IVD analysis is a novel approach for evaluating intra-genomic variations, and exploitation of IVD might be a new road for plant breeding, especially for polyploid crops, such as wheat.

Unraveling the genomic reorganization of polygalacturonase-inhibiting proteins in chickpea

Polygalacturonase-inhibiting proteins (PGIPs) are cell wall proteins that inhibit pathogen polygalacturonases (PGs). PGIPs, like other defense-related proteins, contain extracellular leucine-rich repeats (eLRRs), which are required for pathogen PG recognition. The importance of these PGIPs in plant defense has been well documented. This study focuses on chickpea (Cicer arietinum) PGIPs (CaPGIPs) owing to the limited information available on this important crop. This study identified two novel CaPGIPs (CaPGIP3 and CaPGIP4) and computationally characterized all four CaPGIPs in the gene family, including the previously reported CaPGIP1 and CaPGIP2. The findings suggest that CaPGIP1, CaPGIP3, and CaPGIP4 proteins possess N-terminal signal peptides, ten LRRs, theoretical molecular mass, and isoelectric points comparable to other legume PGIPs. Phylogenetic analysis and multiple sequence alignment revealed that the CaPGIP1, CaPGIP3, and CaPGIP4 amino acid sequences are similar to the other PGIPs reported in legumes. In addition, several cis-acting elements that are typical of pathogen response, tissue-specific activity, hormone response, and abiotic stress-related are present in the promoters of CaPGIP1, CaPGIP3, and CaPGIP4 genes. Localization experiments showed that CaPGIP1, CaPGIP3, and CaPGIP4 are located in the cell wall or membrane. Transcript levels of CaPGIP1, CaPGIP3, and CaPGIP4 genes analyzed at untreated conditions show varied expression patterns analogous to other defense-related gene families. Interestingly, CaPGIP2 lacked a signal peptide, more than half of the LRRs, and other characteristics of a typical PGIP and subcellular localization indicated it is not located in the cell wall or membrane. The study's findings demonstrate CaPGIP1, CaPGIP3, and CaPGIP4's similarity to other legume PGIPs and suggest they might possess the potential to combat chickpea pathogens.

Whole-genome and dispersed duplication, including transposed duplication, jointly advance the evolution of TLP genes in seven representative Poaceae lineages

BACKGROUND

In the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages.

RESULTS

Our research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations.

CONCLUSIONS

This study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.

Genome-Wide Analysis of Proline-Rich Extensin-Like Receptor Kinases (PERKs) Gene Family Reveals Their Roles in Plant Development and Stress Conditions in Oryza sativa L

Proline-rich extensin-like receptor kinases (PERKs) play a crucial role in a wide range of biological processes in plants. In model plants like Arabidopsis, the PERK gene family has been well investigated. Conversely, no information available on the PERK gene family and their biological functions largely remained unknown in rice. This study analyzed the basic physicochemical properties, phylogeny, gene structure, cis-acting elements, Gene ontology (GO) annotation and protein-protein interaction of OsPERK gene family members using various bioinformatics tools based on the whole-genome data of O. sativa. Thus, in this work, 8 PERK genes in rice were identified, and their roles in plant development, growth, and response to various stresses were studied. A phylogenetic study revealed that OsPERKs are grouped into seven classes. Chromosomal mapping also displayed that 8 PERK genes were unevenly distributed on 12 chromosomes. Further, the prediction of subcellular localization indicated that OsPERKs were mainly located at the endomembrane system. Gene structure analysis of OsPERKs has shown a distinctive evolutionary path. In addition, synteny analysis exhibited the 40 orthologous gene pairs in Arabidopsis thaliana, Triticum aestivum, Hordeum vulgare and Medicago truncatula. Furthermore, Ka to Ks proportion shows that most OsPERK genes experienced resilient purifying selection during evolutionary processes. The OsPERK promoters contained several cis-acting regulatory, which are crucial for plant development processes, phytohormone signaling, stress, and defense response. Moreover, the expression pattern of OsPERK family members showed differential expression patterns in different tissues and various stress conditions. Taken together, these results provide clear messages for a better understanding the roles of OsPERK genes in various development stages, tissues, and multifactorial stress as well as enriched the related research of OsPERK family members in rice.

Genome-Wide Identification and Expression Analysis of the Dof Transcription Factor in Annual Alfalfa Medicago polymorpha

The Dof transcription factor is a plant-specific transcription gene family that plays various biological functions in plant development and stress response. However, no relevant research has been conducted on Medicago polymorpha. Here, 36 MpDof genes were identified in the M. polymorpha genome and further divided into 10 groups based on the comparative phylogenetic analysis. The essential information of MpDof genes, such as chromosomal localization, gene structure, conserved motifs, and selective pressures were systematically analyzed. All 36 MpDof genes were predicted to contain more cis-acting elements related to hormone response. MpDof24 and MpDof25 were predicted to interact with MpDof11 and MpDof26 to involve in the photoperiod blooms process. The MpDof genes showed a diverse expression pattern in different tissues. Notably, MpDof29 and MpDof31 were specifically expressed in the large pod and root, respectively, suggesting their crucial role in the pod and root development. qRT-PCR analysis indicated that the expression levels of MpDof10, MpDof25, MpDof26, and MpDof29 were obviously up-regulated under drought, salt, and cold stress. Collectively, genome-wide identification, evolutionary, and expression analysis of the Dof transcription gene family in M. polymorpha will provide new information to further understand and utilize the function of these Dof genes in Medicago plants.

Genome-wide characterization and comparative analysis of the OSCA gene family and identification of its potential stress-responsive members in legumes

Cicer arietinum, Cajanus cajan, Vigna radiata, and Phaseolus vulgaris are economically important legume crops with high nutritional value. They are negatively impacted globally by different biotic and abiotic stresses. Hyperosmolality-gated calcium-permeable channels (OSCA) have been characterized as osmosensors in Arabidopsis thaliana but have not previously reported in legumes. This study provides a genome-wide identification, characterization, and comparative analysis of OSCA genes in legumes. Our study identified and characterized 13 OSCA genes in C. cajan, V. radiata, P. vulgaris, and 12 in C. arietinum, classified into four distinct clades. We found evidence to suggest that the OSCAs might be involved in the interaction between hormone signalling pathways and stress signalling pathways. Furthermore, they play a major role in plant growth and development. The expression levels of the OSCAs vary under different stress conditions in a tissue-specific manner. Our study can be used to develop a detailed understanding of stress regulatory mechanisms of the OSCA gene family in legumes.

Mutational Bias and Natural Selection Driving the Synonymous Codon Usage of Single-Exon Genes in Rice (Oryza sativa L.)

The relative abundance of single-exon genes (SEGs) in higher plants is perplexing. Uncovering the synonymous codon usage pattern of SEGs will benefit for further understanding their underlying evolutionary mechanism in plants. Using internal correspondence analysis (ICA), we reveal a significant difference in synonymous codon usage between SEGs and multiple-exon genes (MEGs) in rice. But the effect is weak, accounting for only 2.61% of the total codon usage variability. SEGs and MEGs contain remarkably different base compositions, and are under clearly differential selective constraints, with the former having higher GC content, and evolving relatively faster during evolution. In the group of SEGs, the variability in synonymous codon usage among genes is partially due to the variations in GC content, gene function, and gene expression level, which accounts for 22.03%, 5.99%, and 3.32% of the total codon usage variability, respectively. Therefore, mutational bias and natural selection should work on affecting the synonymous codon usage of SEGs in rice. These findings may deepen our knowledge for the mechanisms of origination, differentiation and regulation of SEGs in plants.

Trends in the evolution of intronless genes in Poaceae

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

Decipher the molecular evolution and expression patterns of Cupin family genes in oilseed rape

Cupin proteins are involved in plant growth and development as well as in response to various stresses. Here, a total of 173 Cupin genes were identified in Brassica napus, and their molecular evolution and expression patterns were analyzed. These genes were classified into ten groups. Motif and exon-intron structure indicated a high degree of conservation within each group during evolution. BnaCupins were distributed on 19 chromosomes and their expansion is mainly contributed by whole-genome duplication (WGD) and segmental duplication events. BnaCupins have undergone severe purifying selection during a long evolutionary process. Meanwhile, some positive selection sites were identified. Expression patterns and cis-element analysis indicated that BnaCupins play significant roles in plant growth and stress responses. In addition, the expression levels of some BnCupins were significantly altered when treated with different conditions (cold, salt, drought, IAA, ABA, and 6-BA). Some BnaCupin interacting proteins, such as glycosyl hydrolase5 (GHs5), carbohydrate kinase (CHKs), ATP-dependent 6-phosphofructokinase (ATP-PFK), S-adenosylmethionine synthase (S-MAT), and aldolase class II (ALD II), were identified by the protein-protein interaction network. It will contribute to enriching our knowledge of the Cupin gene family in B. napus and provide a basis for further studies of their functions.

A deeply conserved protease, acylamino acid-releasing enzyme (AARE), acts in ageing in Physcomitrella and Arabidopsis

Reactive oxygen species (ROS) are constant by-products of aerobic life. In excess, ROS lead to cytotoxic protein aggregates, which are a hallmark of ageing in animals and linked to age-related pathologies in humans. Acylamino acid-releasing enzymes (AARE) are bifunctional serine proteases, acting on oxidized proteins. AARE are found in all domains of life, albeit under different names, such as acylpeptide hydrolase (APEH/ACPH), acylaminoacyl peptidase (AAP), or oxidized protein hydrolase (OPH). In humans, AARE malfunction is associated with age-related pathologies, while their function in plants is less clear. Here, we provide a detailed analysis of AARE genes in the plant lineage and an in-depth analysis of AARE localization and function in the moss Physcomitrella and the angiosperm Arabidopsis. AARE loss-of-function mutants have not been described for any organism so far. We generated and analysed such mutants and describe a connection between AARE function, aggregation of oxidized proteins and plant ageing, including accelerated developmental progression and reduced life span. Our findings complement similar findings in animals and humans, and suggest a unified concept of ageing may exist in different life forms.

Genome-Wide Identification of DUF668 Gene Family and Expression Analysis under Drought and Salt Stresses in Sweet Potato [Ipomoea batatas (L.) Lam]

The domain of unknown function 668 (DUF668) is a gene family that plays a vital role in responses to adversity coercion stresses in plant. However, the function of the DUF668 gene family is not fully understood in sweet potato. In this study, bioinformatics methods were used to analyze the number, physicochemical properties, evolution, structure, and promoter cis-acting elements of the IbDUF668 family genes, and RNA-seq and qRT-PCR were performed to detect gene expression and their regulation under hormonal and abiotic stress. A total of 14 IbDUF668 proteins were identified in sweet potato, distributed on nine chromosomes. By phylogenetic analysis, IbDUF668 proteins can be divided into two subfamilies. Transcriptome expression profiling revealed that many genes from DUF668 in sweet potato showed specificity and differential expression under cold, heat, drought, salt and hormones (ABA, GA3 and IAA). Four genes (IbDUF668-6, 7, 11 and 13) of sweet potato were significantly upregulated by qRT-PCR under ABA, drought and NaCl stress. Results suggest that the DUF668 gene family is involved in drought and salt tolerance in sweet potato, and it will further provide the basic information of DUF668 gene mechanisms in plants.

The revealing of a novel lipid transfer protein lineage in green algae

BACKGROUND

Non-specific lipid transfer proteins (nsLTPs) are a group of small and basic proteins that can bind and transfer various lipid molecules to the apoplastic space. A typical nsLTP carries a conserved architecture termed eight-cysteine motif (8CM), a scaffold of loop-linked helices folding into a hydrophobic cavity for lipids binding. Encoded by a multigene family, nsLTPs are widely distributed in terrestrial plants from bryophytes to angiosperms with dozens of gene members in a single species. Although the nsLTPs in the most primitive plants such as Marchantia already reach 14 members and are divergent enough to form separate groups, so far none have been identified in any species of green algae.

RESULTS

By using a refined searching strategy, we identified putative nsLTP genes in more than ten species of green algae as one or two genes per haploid genome but not in red and brown algae. The analyses show that the algal nsLTPs carry unique characteristics, including the extended 8CM spacing, larger molecular mass, lower pI value and multiple introns in a gene, which suggests that they could be a novel nsLTP lineage. Moreover, the results of further investigation on the two Chlamydomonas nsLTPs using transcript and protein assays demonstrated their late zygotic stage expression patterns and the canonical nsLTP properties were also verified, such as the fatty acids binding and proteinase resistance activities.

CONCLUSIONS

In conclusion, a novel nsLTP lineage is identified in green algae, which carries some unique sequences and molecular features that are distinguishable from those in land plants. Combined with the results of further examinations of the Chlamydomonas nsLTPs in vitro, possible roles of the algal nsLTPs are also suggested. This study not only reveals the existence of the nsLTPs in green algae but also contributes to facilitating future studies on this enigmatic protein family.

Functional Characterization of TkSRPP Promoter in Response to Hormones and Wounding Stress in Transgenic Tobacco

Taraxacum kok-saghyz is a model species for studying natural rubber biosynthesis because its root can produce high-quality rubber. Small rubber particle protein (SRPP), a stress-related gene to multiple stress responses, involves in natural rubber biosynthesis. To investigate the transcriptional regulation of the TkSRPP promoter, the full-length promoter PR0 (2188 bp) and its four deletion derivatives, PR1 (1592 bp), PR2 (1274 bp), PR3 (934 bp), and PR4 (450 bp), were fused to β-glucuronidase (GUS) reporter gene and transformed into tobacco. The GUS tissue staining showed that the five promoters distinctly regulated GUS expression utilizing transient transformation of tobacco. The GUS activity driven by a PR0 promoter was detected in transgenic tobacco leaves, stem and roots, suggesting that the TkSRPP promoter was not tissue-specific. Deletion analyses in transgenic tobacco have demonstrated that the PR3 from -934 bp to -450 bp core region responded strongly to the hormones, methyl jasmonate (MeJA), abscisic acid (ABA), and salicylic acid (SA), and also to injury induction. The TkSRPP gene was highly expressed under hormones and wound-induced conditions. This study reveals the regulation pattern of the SRPP promoter, and provides valuable information for studying natural rubber biosynthesis under hormones and wounding stress.

Genome-Wide Identification and Expression Profiling of Aconitase Gene Family Members Reveals Their Roles in Plant Development and Adaptation to Diverse Stress in Triticum aestivum L

Global warming is a serious threat to food security and severely affects plant growth, developmental processes, and, eventually, crop productivity. Respiratory metabolism plays a critical role in the adaptation of diverse stress in plants. Aconitase (ACO) is the main enzyme, which catalyzes the revocable isomerization of citrate to isocitrate in the Krebs cycle. The function of ACO gene family members has been extensively studied in model plants, for instance Arabidopsis. However, their role in plant developmental processes and various stress conditions largely remained unknown in other plant species. Thus, we identified 15 ACO genes in wheat to elucidate their function in plant developmental processes and different stress environments. The phylogenetic tree revealed that TaACO genes were classified into six groups. Further, gene structure analysis of TaACOs has shown a distinctive evolutionary path. Synteny analysis showed the 84 orthologous gene pairs in Brachypodium distachyon, Aegilops tauschii, Triticum dicoccoides, Oryza sativa, and Arabidopsis thaliana. Furthermore, Ka/Ks ratio revealed that most TaACO genes experienced strong purifying selection during evolution. Numerous cis-acting regulatory elements were detected in the TaACO promoters, which play a crucial role in plant development processes, phytohormone signaling, and are related to defense and stress. To understand the function of TaACO genes, the expression profiling of TaACO genes were investigated in different tissues, developmental stages, and stress conditions. The transcript per million values of TaACOs genes were retrieved from the Wheat Expression Browser Database. We noticed the differential expression of the TaACO genes in different tissues and various stress conditions. Moreover, gene ontology analysis has shown enrichment in the tricarboxylic acid metabolic process (GO:0072350), citrate metabolic process (GO:0006101), isocitrate metabolic process GO:0006102, carbohydrate metabolic (GO:0005975), and glyoxylate metabolic process (GO:0046487). Therefore, this study provided valuable insight into the ACO gene family in wheat and contributed to the further functional characterization of TaACO during different plant development processes and various stress conditions.

Dissection of leucine-rich repeat receptor-like protein kinases: insight into resistance to Fusarium wilt in tung tree

The tung tree is a woody oil plant native to China and widely distributed in the subtropics. The three main species commonly known as Vernicia are V. fordii, V. montana, and V. cordata. The growth and development of V. fordii are affected by a large number of plant pathogens, such as Fusarium wilt caused by Fusarium sp. In contrast, V. montana shows significant resistance to Fusarium wilt. The leucine-rich repeat receptor-like protein kinase (LRR-RLK) is the largest class of receptor-like kinases associated with plant resistance to Fusarium wilt. Here, we identified 239 VmLRR-RLKs in V. montana, and found that there were characteristic domains of resistance to Fusarium wilt in them. Phylogenetic analysis suggested that the VmLRR-RLKs are divided into 14 subfamilies, indicating that homologous genes in the same group may have similar functions. Chromosomal localization analysis showed that VmLRR-RLKs were unevenly distributed on chromosomes, and segment duplications were the main reason for the expansion of VmLRR-RLK family members. The transcriptome data showed that six orthologous pairs were up-regulated in V. montana in response to Fusarium wilt, while the corresponding orthologous genes showed low or no expression in V. fordii in resistance Fusarium wilt, further indicating the important role of LRR-RLKs in V. montana's resistance to infection by Fusarium spp. Our study provides important reference genes for the future use of molecular breeding to improve oil yield and control of Fusarium wilt in tung tree.

Comprehensive Identification and Functional Analysis of Stress-Associated Protein (SAP) Genes in Osmotic Stress in Maize

Stress-associated proteins (SAPs) are a kind of zinc finger protein with an A20/AN1 domain and contribute to plants' adaption to various abiotic and biological stimuli. However, little is known about the SAP genes in maize (Zea mays L.). In the present study, the SAP genes were identified from the maize genome. Subsequently, the protein properties, gene structure and duplication, chromosomal location, and cis-acting elements were analyzed by bioinformatic methods. Finally, their expression profiles under osmotic stresses, including drought and salinity, as well as ABA, and overexpression in Saccharomyces cerevisiae W303a cells, were performed to uncover the potential function. The results showed that a total of 10 SAP genes were identified and named ZmSAP1 to ZmSAP10 in maize, which was unevenly distributed on six of the ten maize chromosomes. The ZmSAP1, ZmSAP4, ZmSAP5, ZmSAP6, ZmSAP7, ZmSAP8 and ZmSAP10 had an A20 domain at N terminus and AN1 domain at C terminus, respectively. Only ZmSAP2 possessed a single AN1 domain at the N terminus. ZmSAP3 and ZmSAP9 both contained two AN1 domains without an A20 domain. Most ZmSAP genes lost introns and had abundant stress- and hormone-responsive cis-elements in their promoter region. The results of quantitative real-time PCR showed that all ZmSAP genes were regulated by drought and saline stresses, as well as ABA induction. Moreover, heterologous expression of ZmSAP2 and ZmSAP7 significantly improved the saline tolerance of yeast cells. The study provides insights into further underlying the function of ZmSAPs in regulating stress response in maize.

ZmDREB2.9 Gene in Maize (Zea mays L.): Genome-Wide Identification, Characterization, Expression, and Stress Response

Dehydration-responsive element-binding (DREB) transcription factors of the A2 subfamily play key roles in plant stress responses. In this study, we identified and characterized a new A2-type DREB gene, ZmDREB2.9, in the Zea mays cv. B73 genome and compared its expression profile with those of the known A2-type maize genes ZmDREB2.1-2.8. ZmDREB2.9 was mapped to chromosome 8, contained 18 predicted hormone- and stress-responsive cis-elements in the promoter, and had two splice isoforms: short ZmDREB2.9-S preferentially expressed in the leaves, embryos, and endosperm and long ZmDREB2.9-L expressed mostly in the male flowers, stamens, and ovaries. Phylogenetically, ZmDREB2.9 was closer to A. thaliana DREB2A than the other ZmDREB2 factors. ZmDREB2.9-S, ZmDREB2.2, and ZmDREB2.1/2A were upregulated in response to cold, drought, and abscisic acid and may play redundant roles in maize stress resistance. ZmDREB2.3, ZmDREB2.4, and ZmDREB2.6 were not expressed in seedlings and could be pseudogenes. ZmDREB2.7 and ZmDREB2.8 showed similar transcript accumulation in response to cold and abscisic acid and could be functionally redundant. Our results provide new data on Z. mays DREB2 factors, which can be used for further functional studies as well as in breeding programs to improve maize stress tolerance.

AP2/ERF genes associated with superfast fig (Ficus carica L.) fruit ripening

Fig fruits have significant health value and are culturally important. Under suitable climatic conditions, fig fruits undergo a superfast ripening process, nearly doubling in size, weight, and sugar content over three days in parallel with a sharp decrease in firmness. In this study, 119 FcAP2/ERF genes were identified in the fig genome, namely 95 ERFs, 20 AP2s, three RAVs, and one soloist. Most of the ERF subfamily members (76) contained no introns, whereas the majority of the AP2 subfamily members had at least two introns each. Three previously published transcriptome datasets were mined to discover expression patterns, encompassing the fruit peel and flesh of the 'Purple Peel' cultivar at six developmental stages; the fruit receptacle and flesh of the 'Brown Turkey' cultivar after ethephon treatment; and the receptacle and flesh of parthenocarpic and pollinated fruits of the 'Brown Turkey' cultivar. Eighty-three FcAP2/ERFs (68 ERFs, 13 AP2s, one RAV, and one soloist) were expressed in the combined transcriptome dataset. Most FcAP2/ERFs were significantly downregulated (|log₂(fold change) | ≥ 1 and p-adjust < 0.05) during both normal fruit development and ethephon-induced accelerated ripening, suggesting a repressive role of these genes in fruit ripening. Five significantly downregulated ERFs also had repression domains in the C-terminal. Seven FcAP2/ERFs were identified as differentially expressed during ripening in all three transcriptome datasets. These genes were strong candidates for future functional genetic studies to elucidate the major FcAP2/ERF regulators of the superfast fig fruit ripening process.

Genome-wide identification and characterization of AP2/ERF gene superfamily during flower development in Actinidia eriantha

Background

As one of the largest transcription factor families in plants, AP2/ERF gene superfamily plays important roles in plant growth, development, fruit ripening and biotic and abiotic stress responses. Despite the great progress has been made in kiwifruit genomic studies, little research has been conducted on the AP2/ERF genes of kiwifruit. The increasing kiwifruit genome resources allowed us to reveal the tissue expression profiles of AP2/ERF genes in kiwifruit on a genome-wide basis.

Results

In present study, a total of 158 AP2/ERF genes in A. eriantha were identified. All genes can be mapped on the 29 chromosomes. Phylogenetic analysis divided them into four main subfamilies based on the complete protein sequences. Additionally, our results revealed that the same subfamilies contained similar gene structures and conserved motifs. Ka/Ks calculation indicated that AP2/ERF gene family was undergoing a strong purifying selection and the evolutionary rates were slow. RNA-seq showed that the AP2/ERF genes were expressed differently in different flower development stages and 56 genes were considered as DEGs among three contrasts. Moreover, qRT-PCR suggested partial genes showed significant expressions as well, suggesting they could be key regulators in flower development in A. eriantha. In addition, two genes (AeAP2/ERF061, AeAP2/ERF067) had abundant transcription level based on transcriptomes, implying that they may play a crucial role in plant flower development regulation and flower tissue forming.

Conclusions

We identified AP2/ERF genes and demonstrated their gene structures, conserved motifs, and phylogeny relationships of AP2/ERF genes in two related species of kiwifruit, A. eriantha and A. chinensis, and their potential roles in flower development in A. eriantha. Such information would lay the foundation for further functional identification of AP2/ERF genes involved in kiwifruit flower development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08871-4.

The ABCISIC ACID INSENSITIVE (ABI) 4 Transcription Factor Is Stabilized by Stress, ABA and Phosphorylation

The Arabidopsis transcription factor ABSCISIC ACID INSENSITIVE 4 (ABI4) is a key player in the plant hormone abscisic acid (ABA) signaling pathway and is involved in plant response to abiotic stress and development. Expression of the ABI4 gene is tightly regulated, with low basal expression. Maximal transcript levels occur during the seed maturation and early seed germination stages. Moreover, ABI4 is an unstable, lowly expressed protein. Here, we studied factors affecting the stability of the ABI4 protein using transgenic Arabidopsis plants expressing 35S::HA-FLAG-ABI4-eGFP. Despite the expression of eGFP-tagged ABI4 being driven by the highly active 35S CaMV promoter, low steady-state levels of ABI4 were detected in the roots of seedlings grown under optimal conditions. These levels were markedly enhanced upon exposure of the seedlings to abiotic stress and ABA. ABI4 is degraded rapidly by the 26S proteasome, and we report on the role of phosphorylation of ABI4-serine 114 in regulating ABI4 stability. Our results indicate that ABI4 is tightly regulated both post-transcriptionally and post-translationally. Moreover, abiotic factors and plant hormones have similar effects on ABI4 transcripts and ABI4 protein levels. This double-check mechanism for controlling ABI4 reflects its central role in plant development and cellular metabolism.

Molecular evolutionary analysis of the SHI/STY gene family in land plants: A focus on the Brassica species

The plant-specific SHORT INTERNODES/STYLISH (SHI/STY) proteins belong to a family of transcription factors that are involved in the formation and development of early lateral roots. However, the molecular evolution of this family is rarely reported. Here, a total of 195 SHI/STY genes were identified in 21 terrestrial plants, and the Brassica species is the focus of our research. Their physicochemical properties, chromosome location and duplication, motif distribution, exon-intron structures, genetic evolution, and expression patterns were systematically analyzed. These genes are divided into four clades (Clade 1/2/3/4) based on phylogenetic analysis. Motif distribution and gene structure are similar in each clade. SHI/STY proteins are localized in the nucleus by the prediction of subcellular localization. Collinearity analysis indicates that the SHI/STYs are relatively conserved in evolution. Whole-genome duplication is the main factor for their expansion. SHI/STYs have undergone intense purifying selection, but several positive selection sites are also identified. Most promoters of SHI/STY genes contain different types of cis-elements, such as light, stress, and hormone-responsive elements, suggesting that they may be involved in many biological processes. Protein-protein interaction predicted some important SHI/STY interacting proteins, such as LPAT4, MBOATs, PPR, and UBQ3. In addition, the RNA-seq and qRT-PCR analysis were studied in detail in rape. As a result, SHI/STYs are highly expressed in root and bud, and can be affected by Sclerotinia sclerotiorum, drought, cold, and heat stresses. Moreover, quantitative real-time PCR (qRT-PCR) analyses indicates that expression levels of BnSHI/STYs are significantly altered in different treatments (cold, salt, drought, IAA, auxin; ABA, abscisic acid; 6-BA, cytokinin). It provides a new understanding of the evolution and expansion of the SHI/STY family in land plants and lays a foundation for further research on their functions.

CFVisual: an interactive desktop platform for drawing gene structure and protein architecture

BACKGROUND

When researchers perform gene family analysis, they often analyze the structural characteristics of the gene, such as the distribution of introns and exons. At the same time, characteristic structural analysis of amino acid sequence is also essential, for example, motif and domain features. Researchers often integrate these analyses into one image to dig out more information, but the tools responsible for this integration are lacking.

RESULTS

Here, we developed a tool (CFVisual) for drawing gene structure and protein architecture. CFVisual can draw the phylogenetic tree, gene structure, and protein architecture in one picture, and has rich interactive capabilities, which can meet the work needs of researchers. Furthermore, it also supports arbitrary stitching of the above analysis images. It has become a useful helper in gene family analysis. The CFVisual package was implemented in Python and is freely available from https://github.com/ChenHuilong1223/CFVisual/ .

CONCLUSION

CFVisual has been used by some researchers and cited by some articles. In the future, CFVisual will continue to serve as a good helper for researchers in the study of gene structure and protein architecture.

Genome-Wide Characterization, Evolutionary Analysis of ARF Gene Family, and the Role of SaARF4 in Cd Accumulation of Sedum alfredii Hance

Auxin response factors (ARFs) play important roles in plant development and environmental adaption. However, the function of ARFs in cadmium (Cd) accumulation are still unknown. Here, 23 SaARFs were detected in the genome of hyperaccumulating ecotype of Sedum alfredii Hance (HE), and they were not evenly distributed on the chromosomes. Their protein domains remained highly conservative. SaARFs in the phylogenetic tree can be divided into three groups. Genes in the group Ⅰ contained three introns at most. However, over ten introns were found in other two groups. Collinearity relationships were exhibited among ten SaARFs. The reasons for generating SaARFs may be segmental duplication and rearrangements. Collinearity analysis among different species revealed that more collinear genes of SaARFs can be found in the species with close relationships of HE. A total of eight elements in SaARFs promoters were related with abiotic stress. The qRT-PCR results indicated that four SaARFs can respond to Cd stress. Moreover, that there may be functional redundancy among six SaARFs. The adaptive selection and functional divergence analysis indicated that SaARF4 may undergo positive selection pressure and an adaptive-evolution process. Overexpressing SaARF4 effectively declined Cd accumulation. Eleven single nucleotide polymorphism (SNP) sites relevant to Cd accumulation can be detected in SaARF4. Among them, only one SNP site can alter the sequence of the SaARF4 protein, but the SaARF4 mutant of this site did not cause a significant difference in cadmium content, compared with wild-type plants. SaARFs may be involved in Cd-stress responses, and SaARF4 may be applied for decreasing Cd accumulation of plants.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Thaumatin-like Protein (TLP) Genes in Garlic (Allium sativum L.): Genome-Wide Identification, Characterization, and Expression in Response to Fusarium proliferatum Infection

Plant antifungal proteins include the pathogenesis-related (PR)-5 family of fungi- and other stress-responsive thaumatin-like proteins (TLPs). However, the information on the TLPs of garlic (Allium sativum L.), which is often infected with soil Fusarium fungi, is very limited. In the present study, we identified 32 TLP homologs in the A. sativum cv. Ershuizao genome, which may function in the defense against Fusarium attack. The promoters of A. sativumTLP (AsTLP) genes contained cis-acting elements associated with hormone signaling and response to various types of stress, including those caused by fungal pathogens and their elicitors. The expression of AsTLP genes in Fusarium-resistant and -susceptible garlic cultivars was differently regulated by F. proliferatum infection. Thus, in the roots the mRNA levels of AsTLP7–9 and 21 genes were increased in resistant and decreased in susceptible A. sativum cultivars, suggesting the involvement of these genes in the garlic response to F. proliferatum attack. Our results provide insights into the role of TLPs in garlic and may be useful for breeding programs to increase the resistance of Allium crops to Fusarium infections.

Genome-wide identification and expression analysis of the VQ gene family in Cucurbita pepo L

VQ protein is a plant specific protein, which plays an important role in plant growth and development and biological and abiotic stress response. This study aimed to systematically analyze for the first time the VQ of Cucurbita pepo and understand their expression patterns in response to different stimuli. Herein, 44 VQ genes were identified, which were divided into eight groups (I-VIII) based on phylogenetic analysis. Two genes (CpVQ1 and CpVQ2) could not be located on the chromosome, whereas the remaining CpVQ genes were randomly distributed on the chromosomes, except for chromosomes 15 and 18. Noteworthy, the main event driving the expansion of the VQ gene family was chromosome fragment duplication. Based on qRT-PCR analysis, VQ genes are expressed in different tissues, and VQ genes are differentially regulated under a variety of abiotic stresses and powdery mildew stress, indicating that they play an important role in plant stress response and other aspects. This report presents the first systematic analysis of VQ genes from C. pepo and provides a solid foundation for further research of the specific functions of VQ proteins.

Genome-Wide Analysis and Expression Profiles of Ethylene Signal Genes and Apetala2/Ethylene-Responsive Factors in Peanut (Arachis hypogaea L.)

Peanut is an important oil and economic crop widely cultivated in the world. It has special characteristics such as blooming on the ground but bearing fruits underground. During the peg penetrating into the ground, it is subjected to mechanical stress from the soil at the same time. It has been proved that mechanical stress affects plant growth and development by regulating the ethylene signaling-related genes. In this study, we identified some genes related to ethylene signal of peanut, including 10 ethylene sensors, two constitutive triple responses (CTRs), four ethylene insensitive 2 (EIN2s), four ethylene insensitive 3 (EIN3s), six EIN3-binding F-box proteins (EBFs), and 188 Apetala2/ethylene-responsive factors (AP2/ERFs). One hundred and eighty-eight AP2/ERFs were further divided into four subfamilies, 123 ERFs, 56 AP2s, 6 Related to ABI3/VP1 (RAVs), and three Soloists, of them one hundred and seventy AP2/ERF gene pairs were clustered into segmental duplication events in genome of Arachis hypogaea. A total of 134, 138, 97, and 150 AhAP2/ERF genes formed 210, 195, 166, and 525 orthologous gene pairs with Arachis duranensis, Arachis ipaensis, Arabidopsis thaliana, and Glycine max, respectively. Our transcriptome results showed that two EIN3s (Arahy.J729H0 and Arahy.S7XF8N) and one EBFs (Arahy.G4JMEM) were highly expressed when mechanical stress increased. Among the 188 AhAP2/ERF genes, there were 31 genes with the fragments per kilobase of exon model per million mapped fragments (FPKM) ≥ 100 at least one of the 15 samples of Tifrunner. Among them, three AhAP2/ERFs (Arahy.15RATX, Arahy.FAI7YU, and Arahy.452FBF) were specifically expressed in seeds and five AhAP2/ERFs (Arahy.HGAZ7D, Arahy.ZW7540, Arahy.4XS3FZ, Arahy.QGFJ76, and Arahy.AS0C7C) were highly expressed in the tissues, which responded mechanical stress, suggesting that they might sense mechanical stress. Mechanical stress simulation experiment showed that three AhAP2/ERFs (Arahy.QGFJ76, Arahy.AS0C7C, and Arahy.HGAZ7D) were sensitive to mechanical stress changes and they all had the conservative repressor motif (DLNXXP) in the C-terminus, indicated that they might transmit mechanical stress signals through transcriptional inhibition. This study reveals the regulatory landscape of ethylene signal-related genes in peanut, providing valuable information for the mining of target genes for further study.

Genome-wide characterization of chalcone synthase genes in sweet cherry and functional characterization of CpCHS1 under drought stress

Cherries are one of the important fruit trees. The growth of cherry is greatly affected by abiotic stresses such as drought, which hinders its development. Chalcone synthase (CHS, EC 2.3.1.74) is a crucial rate-limiting enzyme in the flavonoid biosynthetic pathway that plays an important role in regulating plant growth, development, and abiotic stress tolerance. In the current study, three genes encoding chalcone synthase were identified in the genome of sweet cherry (Prunus avium L.). The three genes contained fewer introns and showed high homology with CHS genes of other Rosaceae members. All members are predicted to localize in the cytoplasm. The conserved catalytic sites may be located at the Cys163, Phe214, His302, and Asn335 residues. These genes were differentially expressed during flower bud dormancy and fruit development. The total flavonoid content of Chinese cherry (Cerasus pseudocerasus Lindl.) was highest in the leaves and slightly higher in the pulp than in the peel. No significant difference in total flavonoid content was detected between aborted kernels and normally developing kernels. Overexpression of Chinese cherry CpCHS1 in tobacco improved the germination frequency of tobacco seeds under drought stress, and the fresh weight of transgenic seedlings under drought stress was higher than that of the wild type, and the contents of SOD, POD, CAT, and Pro in OE lines were significantly increased and higher than WT under drought stress. These results indicate cherry CHS genes are conserved and functionally diverse and will assist in elucidating the functions of flavonoid synthesis pathways in cherry and other Rosaceae species under drought stress.

Differential Functions of Pepper Stress-Associated Proteins in Response to Abiotic Stresses

Stress-associated proteins (SAPs), a group of zinc-finger-type proteins, have been identified as novel regulators of plant abiotic and biotic stresses. However, although they have been discovered in different plant species, their precise functional roles remain unclear. Here, we identified 14 SAP subfamily genes in the pepper genome. An investigation of the promoter regions of these genes for cis-regulatory elements associated with abiotic stress responses revealed the presence of multiple stress-related elements. Domain and phylogenetic analyses using the corresponding protein sequences revealed that the CaSAP genes can be classified into six groups (I-VI) and sorted into two broad types. Expression levels of the CaSAP genes were found to be differentially induced by low temperature, the dehydration stress, or exogenous abscisic acid. Group II and IV genes were highly induced by the low temperature and dehydration treatments, respectively. Moreover, subcellular localization analysis indicated that the proteins in these two groups are distributed in the nucleus, cytoplasm, and plasma membrane. Among the pepper plants silenced with the three identified group II CaSAP genes, the CA02g10410-silenced plants showed tolerance to low temperature, whereas the CA03g17080-silenced plants were found to have temperature-sensitive phenotypes. Interestingly, group IV CaSAP-silenced pepper plants showed drought-tolerant phenotypes. These findings contribute to a preliminary characterization of CaSAP genes and provide directions for future research on the biological role of CaSAPs in response to different abiotic stresses.

Genome-Wide Identification and Evolutionary Analysis of the SRO Gene Family in Tomato

SRO (SIMILAR TO RCD ONE) is a family of plant-specific small molecule proteins that play an important role in plant growth and development and environmental responses. However, SROs still lack systematic characterization in tomato. Based on bioinformatics methods, SRO family genes were identified and characterized from cultivated tomatoes and several wild tomatoes. qRT-PCR was used to study the expression of SRO gene in cultivated tomatoes. Phylogenetic and evolutionary analyses showed that SRO genes in angiosperms share a common ancestor and that the number of SRO family members changed as plants diverged and evolved. Cultivated tomato had six SRO members, five of which still shared some degree of identity with the ancestral SRO genes. Genetic structure and physicochemical properties showed that tomato SRO genes were highly conserved with chromosomal distribution. They could be divided into three groups based on exon-intron structure, and cultivated tomato contained only two of these subclades. A number of hormonal, light and abiotic stress-responsive cis-regulatory elements were identified from the promoter of the tomato SRO gene, and they also interacted with a variety of stress-responsive proteins and microRNAs. RNA-seq analysis showed that SRO genes were widely expressed in different tissues and developmental stages of tomato, with significant tissue-specific features. Expression analysis also showed that SRO genes respond significantly to high temperature and salt stress and mediate the tomato hormone regulatory network. These results provide a theoretical basis for further investigation of the functional expression of tomato SRO genes and provide potential genetic resources for tomato resistance breeding.

Identification and Analysis of GhEXO Gene Family Indicated That GhEXO7_At Promotes Plant Growth and Development Through Brassinosteroid Signaling in Cotton (Gossypium hirsutum L.)

Brassinosteroids (BRs), an efficient plant endogenous hormone, significantly promotes plant nutrient growth adapting to biological and abiotic adversities. BRs mainly promote plant cell elongation by regulating gene expression patterns. EXORDIUM (EXO) genes have been characterized as the indicators of BR response genes. Cotton, an ancient crop, is of great economic value and its fibers can be made into all kinds of fabrics. However, EXO gene family genes have not been full identified in cotton. 175 EXO genes were identified in nine plant species, of which 39 GhEXO genes in Gossypium hirsutum in our study. A phylogenetic analysis grouped all of the proteins encoded by the EXO genes into five major clades. Sequence identification of conserved amino acid residues among monocotyledonous and dicotyledonous species showed a high level of conservation across the N and C terminal regions. Only 25% the GhEXO genes contain introns besides conserved gene structure and protein motifs distribution. The 39 GhEXO genes were unevenly distributed on the 18 At and Dt sub-genome chromosomes. Most of the GhEXO genes were derived from gene duplication events, while only three genes showed evidence of tandem duplication. Homologous locus relationships showed that 15 GhEXO genes are located on collinear blocks and that all orthologous/paralogous gene pairs had Ka > Ks values, indicating purifying selection pressure. The GhEXO genes showed ubiquitous expression in all eight tested cotton tissues and following exposure to three phytohormones, IAA, GA, and BL. Furthermore, GhEXO7_At was mainly expressed in response to BL treatment, and was predominantly expressed in the fibers. GhEXO7_At was found to be a plasma membrane protein, and its ectopic expression in Arabidopsis mediated BR-regulated plant growth and development with altered expression of DWF4, CPD, KCS1, and EXP5. Additionally, the functions of GhEXO7_At were confirmed by virus-induced gene silencing (VIGS) in cotton. This study will provide important genetic resources for future cotton breeding programs.

Membrane localized thaumatin-like protein from tea (CsTLP) enhanced seed yield and the plant survival under drought stress in Arabidopsis thaliana

Thaumatin-like proteins (TLPs) are pathogenesis-related (PR5) proteins, which are induced in response to various biotic and abiotic stresses. The present work was carried out to clone TLP of Camellia sinensis (CsTLP) and to evaluate the response of transgenic lines of Arabidopsis constitutively expressing CsTLP under drought conditions. Data showed that transgenic lines exhibited lower relative electrolyte leakage and higher water retention capacity as compared to the wild-type (WT) plants under drought stress. In addition, results with confocal microscopy showed CsTLP + GFP fusion protein to be localized in the cell membrane which moved to the intercellular spaces under prolonged drought stress. Expression of CsTLP enhanced seed yield and the plant survival in transgenic lines as compared to the WT plants under drought stress. Results suggested the importance of CsTLP in improving drought tolerance in Arabidopsis.

Cytokinin glucosyl transferases, key regulators of cytokinin homeostasis, have potential value for wheat improvement

The cytokinins, which are N⁶ -substituted adenine derivatives, control key aspects of crop productivity. Cytokinin levels are controlled via biosynthesis by isopentenyl transferase (IPT), destruction by cytokinin oxidase/dehydrogenase (CKX), and inactivation via glucosylation by cytokinin glucosyl transferases (CGTs). While both yield components and tolerance to drought and related abiotic stressors have been positively addressed via manipulation of IPT and/or CKX expression, much less attention has been paid to the CGTs. As naming of the CGTs has been unclear, we suggest COGT, CNGT, CONGT and CNOGT to describe the O-, N- and dual function CGTs. As specific CGT mutants of both rice and arabidopsis showed impacts on yield components, we interrogated the wheat genome database, IWGSC RefSeq v1.0 & v2.0, to investigate wheat CGTs. Besides providing unambiguous names for the 53 wheat CGTs, we show their expression patterns in 70 developmental tissues and their response characteristics to various stress conditions by reviewing more than 1000 RNA-seq data sets. These revealed various patterns of responses and showed expression generally being more limited in reproductive tissues than in vegetative tissues. Multiple cis-regulatory elements are present in the 3 kb upstream of the start codons of the 53 CGTs. Elements associated with abscisic acid, light and methyl jasmonate are particularly over-represented, indicative of the responsiveness of CGTs to the environment. These data sets indicate that CGTs have potential value for wheat improvement and that these could be targeted in TILLING or gene editing wheat breeding programmes.