Genome-wide identification and expression profiling of auxin response factor (ARF) gene family in maize

Genome-wide analysis of ARF gene family and miR160 in avocado (Persea americana Mill.) and their roles in somatic embryogenesis from zygotic embryos

MAIN CONCLUSION

The genome-wide analysis revealed that miRNA160 and PaARFs are involved in avocado somatic embryogenesis and play a role in the low efficiency of embryo induction and the ineffective conversion of embryos into plants. The auxin response transcription factors (ARFs) play a role in signaling the auxin phytohormone Indole-3-acetic Acid (IAA) and are involved in plant growth, development, abiotic stress responses and somatic embryogenesis (SE). In the Lauraceae family, particularly in avocado (Persea americana Mill.), propagation via SE remains challenging due to the low efficiency of embryo induction and ineffective conversion of embryos into plants. This study investigates the phylogenetic relationships and evolutionary history of avocado ARFs (PaARFs). This multigenic family consists of at least 20 members that evolved from a now-extinct common ancestor shared by bryophytes and angiosperms. The expression profile of these genes was analyzed in immature zygotic embryos and three SE stages: early globular, late globular and white-opaque. Additionally, we identified six genes that contributed to the formation of a 100% identical single mature miRNA, the miR160. Almost all PaARF genes were upregulated during the embryo-induction stage, while genes such as PaARF1a, PaARF1c, PaARF2a, PaARF2b, and PaARF17 were downregulated at the white-opaque stage. We observed that the expression of miRNA160 differed significantly between the zygotic embryos and the three subsequent development stages. Additionally, free IAA distributions were highly concentrated in immature zygotic embryos. Our results suggest that miR160 and PaARF-mediated auxin signaling play a role in avocado SE, potentially contributing to the low efficiency of SE. This study is the first report on the ARF gene family in avocado. Our findings provide a valuable reference for comparative and functional analyses of ARFs in the context of avocado somatic embryogenesis.

Comprehensive genome-wide analysis of ARF transcription factors in orchardgrass (Dactylis glomerata): the positive regulatory role of DgARF7 in drought resistance

Auxin response factor (ARF), a transcription factor, is crucial in controlling growth, development, and response to environmental stress. Orchardgrass (Dactylis glomerata) is an economically significant, widely cultivated forage grass. However, information on the genome-wide information and functional characterization of ARFs in orchardgrass is limited. This study identified 27 ARF genes based on the orchardgrass genome database. These DgARFs were unevenly distributed across the seven orchardgrass chromosomes and clustered into four classes. Phylogenetic analysis with multispecies of ARF proteins indicated that the ARFs exhibit a relatively conserved evolutionary path. Focusing on hormone signaling responses, DgARF7 demonstrated a potential positive regulatory role in response to 3-indole acetic acid, methyl jasmonate, gibberellin, salicylic acid, and abscisic acid signals. Additionally, exposure to drought stress induced noticeable oscillatory changes in DgARF7 gene. Notably, DgARF7 enhanced drought tolerance through heterologous expression in yeast and overexpression in Arabidopsis. Overexpressed Arabidopsis lines of DgARF7 exhibited a markedly higher relative water content and superoxide dismutase activity, while the malondialdehyde content was significantly decreased compared to wild type under drought stress. DgARF7 also accelerated flowering time by inducing the flowering-related gene expression levels in Arabidopsis. This research provides important insights into the role of DgARF7 in orchardgrass and provides further understanding in molecular breeding.

The Identification of Auxin Response Factors and Expression Analyses of Different Floral Development Stages in Roses

Background/Objectives:Auxin response factors (ARFs) are important in plant growth and development, especially flower development. However, there is limited research on the comprehensive identification and characterization of ARF genes in roses. Methods: We employed bioinformatics tools to identify the ARF genes of roses. These genes were characterized for their phylogenetic relationships, chromosomal positions, conserved motifs, gene structures, and expression patterns. Results: In this study, a total of 17 ARF genes were identified in the genomes of Rosa chinensis 'OB', R. chinensis 'CH', R. rugosa, and R. wichurana. Based on RNA-seq analyses, we found that the ARF genes had diverse transcript patterns in various tissues and cultivars. In 'CH', the expression levels of RcCH_ARFs during different flower-development stages were classified into four clusters. In cluster 3 and cluster 4, RcCH_ARFs were specifically high and low in different stages of floral evocation. Gene expression and phylogenetic analyses showed that RcCH_ARF3, RcCH_ARF4, and RcCH_ARF18 were likely to be the key genes for rose flower development. Conclusions: The identification and characterization of ARF genes in Rosa were investigated. The results presented here provide a theoretical basis for the molecular mechanisms of ARF genes in plant development and flowering for roses, with a broader application for other species in the rose family and for the development of novel cultivars.

Systematic analysis of the ARF gene family in Fagopyrum dibotrys and its potential roles in stress tolerance

The auxin response factor (ARF) is a plant-specific transcription factor that regulates the expression of auxin response genes by binding directly to their promoters. They play an important role in the regulation of plant growth and development, as well as in the response to biotic and abiotic stresses. However, the identification and functional analysis of ARFs in Fagopyrum dibotrys are still unclear. In this study, a total of 26 FdARF genes were identified using bioinformatic methods. Their chromosomal location, gene structure, physical and chemical properties of their encoded protein, subcellular location, phylogenetic tree, conserved motifs and cis-acting elements in FdARF promoters were analyzed. The results showed that 26 FdARF genes were unevenly distributed on 8 chromosomes, with the largest distribution on chromosome 4 and the least distribution on chromosome 3. Most FdARF proteins are located in the nucleus, except for the proteins FdARF7 and FdARF21 located to the cytoplasm and nucleus, while FdARF14, FdARF16, and FdARF25 proteins are located outside the chloroplast and nucleus. According to phylogenetic analysis, 26 FdARF genes were divided into 6 subgroups. Duplication analysis indicates that the expansion of the FdARF gene family was derived from segmental duplication rather than tandem duplication. The prediction based on cis-elements of the promoter showed that 26 FdARF genes were rich in multiple stress response elements, suggesting that FdARFs may be involved in the response to abiotic stress. Expression profiling analysis showed that most of the FdARF genes were expressed in the roots, stems, leaves, and tubers of F. dibotrys, but their expression exhibits a certain degree of tissue specificity. qRT-PCR analysis revealed that most members of the FdARF gene were up- or down-regulated in response to abiotic stress. The results of this study expand our understanding of the functional role of FdARFs in response to abiotic stress and lay a theoretical foundation for further exploration of other functions of FdARF genes.

Functional Divergence of the Closely Related Genes PhARF5 and PhARF19a in Petunia hybrida Flower Formation and Hormone Signaling

The ARF gene family plays a vital role in regulating multiple aspects of plant growth and development. However, detailed research on the role of the ARF family in regulating flower development in petunia and other plants remains limited. This study investigates the distinct roles of PhARF5 and PhARF19a in Petunia hybrida flower development. Phylogenetic analysis identified 29 PhARFs, which were grouped into four clades. VIGS-mediated silencing of PhARF5 and PhARF19a led to notable phenotypic changes, highlighting their non-redundant functions. PhARF5 silencing resulted in reduced petal number and limb abnormalities, while PhARF19a silencing disrupted corolla tube formation and orientation. Both genes showed high expression in the roots, leaves, and corollas, with nuclear localization. The transcriptomic analysis revealed significant overlaps in DEGs between PhARF5 and PhARF19a silencing, indicating shared pathways in hormone metabolism, signal transduction, and stress responses. Phytohormone analysis confirmed their broad impact on phytohormone biosynthesis, suggesting involvement in complex feedback mechanisms. Silencing PhARF5 and PhARF19a led to differential transcription of numerous genes related to hormone signaling pathways beyond auxin signaling, indicating their direct or indirect crosstalk with other phytohormones. However, significant differences in the regulation of these signaling pathways were observed between PhARF5 and PhARF19a. These findings reveal the roles of ARF genes in regulating petunia flower development, as well as the phylogenetic distribution of the PhARFs involved in this process. This study provides a valuable reference for molecular breeding aimed at improving floral traits in the petunia genus and related species.

Genomic identification and expression profiling of DMP genes in oat (Avena sativa) elucidate their responsiveness to seed aging

BACKGROUND

The Domain of unknown function 679 membrane protein (DMP) family, which is unique to plants, plays a crucial role in reproductive development, stress response and aging. A comprehensive study was conducted to identify the DMP gene members of oat (Avena sativa) and to investigate their structural features and tissue-specific expression profiles. Utilizing whole genome and transcriptome data, we analyzed the physicochemical properties, gene structure, cis-acting elements, phylogenetic relationships, conserved structural (CS) domains, CS motifs and expression patterns of the AsDMP family in A. sativa.

RESULTS

The DMP family genes of A. sativa were distributed across 17 chromosomal scaffolds, encompassing a total of 33 members. Based on phylogenetic relationships, the AsDMP genes were classified into five distinct subfamilies. The gene structure also suggests that A. sativa may have undergone an intron loss event during its evolution. Covariance analysis indicates that genome-wide duplication and segmental duplication may be the major contributor to the expansion of the AsDMP gene family. Ka/Ks selective pressure analysis of the AsDMP gene family suggests that DMP gene pairs are generally conserved over evolutionary time. The upstream promoters of these genes contain several cis-acting elements, suggesting a potential role in abiotic stress responses and hormone induction. Transcriptome data revealed that the expression patterns of the DMP genes are involved in tissue and organ development. In this study, the AsDMP genes (AsDMP1, AsDMP19, and AsDMP22) were identified as potential regulators of seed senescence in A. sativa. These genes could serve as candidates for breeding studies focused on seed longevity and anti-aging germplasm in A. sativa. The study provides valuable insights into the regulatory mechanisms of the AsDMP gene family in the aging process of A. sativa germplasm and offers theoretical support for further function investigation into the functions of AsDMP genes and the molecular mechanisms underlying seed anti-aging.

CONCLUSIONS

This study identified the AsDMP genes as being involved in the aging process of A. sativa seeds, marking the first report on the potential role of DMP genes in seed aging for A. sativa.

Characterizing seed dormancy in Epimedium brevicornu Maxim.: Development of novel chill models and determination of dormancy release mechanisms by transcriptomics

PURPOSE

Epimedium brevicornu Maxim. is a perennial persistent C3 plant of the genus Epimedium Linn. in the family Berberaceae that exhibits severe physiological and morphological seed dormancy.We placed mature E. brevicornu seeds under nine stratification treatment conditions and explored the mechanisms of influence by combining seed embryo growth status assessment with related metabolic pathways and gene co-expression analysis.

RESULTS

We identified 3.9 °C as the optimum cold-stratification temperature of E. brevicornu seeds via a chilling unit (CU) model. The best treatment was variable-temperature stratification (10/20 °C, 12/12 h) for 4 months followed by low-temperature stratification (4 °C) for 3 months (4-3). A total of 63801 differentially expressed genes were annotated to 2587 transcription factors (TFs) in 17 clusters in nine treatments (0-0, 0-3, 1-3, 2-3, 3-3, 4-3, 4-2, 4-1, 4-0). Genes specifically highly expressed in the dormancy release treatment group were significantly enriched in embryo development ending in seed dormancy and fatty acid degradation, indicating the importance of these two processes. Coexpression analysis implied that the TF GRF had the most reciprocal relationships with genes, and multiple interactions centred on zf-HD and YABBY as well as on MYB, GRF, and TCP were observed.

CONCLUSION

In this study, analyses of plant hormone signal pathways and fatty acid degradation pathways revealed changes in key genes during the dormancy release of E. brevicornu seeds, providing evidence for the filtering of E. brevicornu seed dormancy-related genes.

Ultrasonic treatment can improve maize seed germination and abiotic stress resistance

Constant-frequency ultrasonic treatment helped to improve seed germination. However, variable-frequency ultrasonic treatment on maize seed germination were rarely reported. In this study, maize seeds were exposed to 20-40 kHz ultrasonic for 40 s. The germination percentage and radicle length of maize seeds increased by 10.4% and 230.5%. Ultrasonic treatment also significantly increased the acid protease, α-amylase, and β-amylase contents by 96.4%, 73.8%, and 49.1%, respectively. Transcriptome analysis showed that 11,475 differentially expressed genes (DEGs) were found in the ultrasonic treatment and control groups, including 5,695 upregulated and 5,780 downregulated. Metabolic pathways and transcription factors (TFs) were significantly enriched among DEGs after ultrasonic treatment. This included metabolism and genetic information processing, that is, ribosome, proteasome, and pyruvate metabolism, sesquiterpenoid, triterpenoid, and phenylpropanoid biosynthesis, and oxidative phosphorylation, as well as transcription factors in the NAC, MYB, bHLH, WRKY, AP2, bZIP, and ARF families. Variable-frequency ultrasonic treatment increased auxin, gibberellin, and salicylic acid by 5.5%, 37.3%, and 28.9%, respectively. Abscisic acid significantly decreased by 33.2%. The related DEGs were upregulated and downregulated to varying degrees. Seed germination under the abiotic stress conditions of salt stress (NaCl solution), drought (PEG solution), and waterlogging (water-saturated sand bed) under ultrasonic treatment were promoted, radicle length was significantly increased by 30.2%, 30.5%, and 27.3%, respectively; and germination percentage by 14.8%, 20.1%, and 21.6%, respectively. These findings provide new insight into the mechanisms through ultrasonic to promote maize seed germination.

Divergent Roles of the Auxin Response Factors in Lemongrass (Cymbopogon flexuosus (Nees ex Steud.) W. Watson) during Plant Growth

Auxin Response Factors (ARFs) make up a plant-specific transcription factor family that mainly couples perception of the phytohormone, auxin, and gene expression programs and plays an important and multi-faceted role during plant growth and development. Lemongrass (Cymbopogon flexuosus) is a representative Cymbopogon species widely used in gardening, beverages, fragrances, traditional medicine, and heavy metal phytoremediation. Biomass yield is an important trait for several agro-economic purposes of lemongrass, such as landscaping, essential oil production, and phytoremediation. Therefore, we performed gene mining of CfARFs and identified 26 and 27 CfARF-encoding genes in each of the haplotype genomes of lemongrass, respectively. Phylogenetic and domain architecture analyses showed that CfARFs can be divided into four groups, among which groups 1, 2, and 3 correspond to activator, repressor, and ETTN-like ARFs, respectively. To identify the CfARFs that may play major roles during the growth of lemongrass plants, RNA-seq was performed on three tissues (leaf, stem, and root) and four developmental stages (3-leaf, 4-leaf, 5-leaf. and mature stages). The expression profiling of CfARFs identified several highly expressed activator and repressor CfARFs and three CfARFs (CfARF3, 18, and 35) with gradually increased levels during leaf growth. Haplotype-resolved transcriptome analysis revealed that biallelic expression dominance is frequent among CfARFs and contributes to their gene expression patterns. In addition, co-expression network analysis identified the modules enriched with CfARFs. By establishing orthologous relationships among CfARFs, sorghum ARFs, and maize ARFs, we showed that CfARFs were mainly expanded by whole-genome duplications, and that the duplicated CfARFs might have been divergent due to differential expression and variations in domains and motifs. Our work provides a detailed catalog of CfARFs in lemongrass, representing a first step toward characterizing CfARF functions, and may be useful in molecular breeding to enhance lemongrass plant growth.

Genome-Wide Characterization of Fructose 1,6-Bisphosphate Aldolase Genes and Expression Profile Reveals Their Regulatory Role in Abiotic Stress in Cucumber

The fructose-1,6-bisphosphate aldolase (FBA) gene family exists in higher plants, with the genes of this family playing significant roles in plant growth and development, as well as response to abiotic stresses. However, systematic reports on the FBA gene family and its functions in cucumber are lacking. In this study, we identified five cucumber FBA genes, named CsFBA1-5, that are distributed randomly across chromosomes. Phylogenetic analyses involving these cucumber FBAs, alongside eight Arabidopsis FBA proteins and eight tomato FBA proteins, were conducted to assess their homology. The CsFBAs were grouped into two clades. We also analyzed the physicochemical properties, motif composition, and gene structure of the cucumber FBAs. This analysis highlighted differences in the physicochemical properties and revealed highly conserved domains within the CsFBA family. Additionally, to explore the evolutionary relationships of the CsFBA family further, we constructed comparative syntenic maps with Arabidopsis and tomato, which showed high homology but only one segmental duplication event within the cucumber genome. Expression profiles indicated that the CsFBA gene family is responsive to various abiotic stresses, including low temperature, heat, and salt. Taken together, the results of this study provide a theoretical foundation for understanding the evolution of and future research into the functional characterization of cucumber FBA genes during plant growth and development.

PbARF19-mediated auxin signaling regulates lignification in pear fruit stone cells

The stone cells in pear fruits cause rough flesh and low juice, seriously affecting the taste. Lignin has been demonstrated as the main component of stone cells. Auxin, one of the most important plant hormone, regulates most physiological processes in plants including lignification. However, the concentration effect and regulators of auxin on pear fruits stone cell formation remains unclear. Here, endogenous indole-3-acetic acid (IAA) and stone cells were found to be co-localized in lignified cells by immunofluorescence localization analysis. The exogenous treatment of different concentrations of IAA demonstrated that the application of 200 µM IAA significantly reduced stone cell content, while concentrations greater than 500 µM significantly increased stone cell content. Besides, 31 auxin response factors (ARFs) were identified in pear genome. Putative ARFs were predicted as critical regulators involved in the lignification of pear flesh cells by phylogenetic relationship and expression analysis. Furthermore, the negative regulation of PbARF19 on stone cell formation in pear fruit was demonstrated by overexpression in pear fruitlets and Arabidopsis. These results illustrated that the PbARF19-mediated auxin signal plays a critical role in the lignification of pear stone cell by regulating lignin biosynthetic genes. This study provides theoretical and practical guidance for improving fruit quality in pear production.

Comprehensive analyses of the ARF gene family in cannabis reveals their potential roles in regulating cannabidiol biosynthesis and male flower development

Background

Cannabidiol (CBD), as an important therapeutic property of the cannabis plants, is mainly produced in the flower organs. Auxin response factors (ARFs) are play a crucial role in flower development and secondary metabolite production. However, the specific roles of ARF gene family in cannabis remain unknown.

Methods

In this study, various bioinformatics analysis of CsARF genes were conducted using online website and bioinformatics, quantitative real time PCR technology was used to investigate the expression patterns of the CsARF gene family in different tissues of different cannabis varieties, and subcellular localization analysis was performed in tobacco leaf.

Results

In this study, 22 CsARF genes were identified and found to be unevenly distributed across 9 chromosomes of the cannabis genome. Phylogenetic analysis revealed that the ARF proteins were divided into 4 subgroups. Duplication analysis identified one pair of segmental/whole-genome duplicated CsARF, and three pairs of tandemly duplicated CsARF. Collinearity analysis revealed that two CsARF genes, CsARF4 and CsARF19, were orthologous in both rice and soybean. Furthermore, subcellular localization analysis showed that CsARF2 was localized in the nucleus. Tissue-specific expression analysis revealed that six genes were highly expressed in cannabis male flowers, and among these genes, 3 genes were further found to be highly expressed at different developmental stages of male flowers. Meanwhile, correlation analysis between the expression level of CsARF genes and CBD content in two cultivars 'H8' and 'Y7' showed that the expression level of CsARF13 was negatively correlated with CBD content, while the expression levels of six genes were positively correlated with CBD content. In addition, most of CsARF genes were responsive to IAA treatment.

Conclusion

Our study laid a foundation for the further studies of CsARFs function in cannabis, and provides candidate genes for breeding varieties with high CBD yield in cannabis production.

The auxin response factor (ARF) gene family in Cyclocarya paliurus: genome-wide identification and their expression profiling under heat and drought stresses

Auxin response factors (ARFs), as the main components of auxin signaling, play a crucial role in various processes of plant growth and development, as well as in stress response. So far, there have been no reports on the genome-wide identification of the ARF transcription factor family in Cyclocarya paliurus, a deciduous tree plant in the family Juglaceae. In this study, a total of 34 CpARF genes were identified based on whole genome sequence, and they were unevenly distributed on 16 chromosomes, with the highest distribution on chromosome 6. Domain analysis of CpARF proteins displayed that 31 out of 34 CpARF proteins contain a typical B3 domain (DBD domain), except CpARF12/ CpARF14/CpARF31, which all belong to Class VI. And 20 CpARFs (58.8%) contain an auxin_IAA binding domain, and are mainly distributed in classes I, and VI. Phylogenetic analysis showed that CpARF was divided into six classes (I-VI), each containing 4, 4, 1, 8, 4, and 13 members, respectively. Gene duplication analysis showed that there are 14 segmental duplications and zero tandem repeats were identified in the CpARF gene family of the C. paliurus genome. The Ka/Ks ratio of duplicate gene pairs indicates that CpARF genes are subjected to strong purification selection pressure. Synteny analysis showed that C. paliurus shared the highest homology in 74 ARF gene pairs with Juglans regia, followed by 73, 51, 25, and 11 homologous gene pairs with Populus trichocarpa, Juglans cathayensis, Arabidopsis, and rice, respectively. Promoter analysis revealed that 34 CpARF genes had cis-elements related to hormones, stress, light, and growth and development except for CpARF12. The expression profile analysis showed that almost all CpARF genes were differentially expressed in at least one tissue, and several CpARF genes displayed tissue-specific expression. Furthermore, 24 out of the 34 CpARF genes have significantly response to drought stress (P < 0.05), and most of them (16) being significantly down-regulated under moderate drought treatment. Meanwhile, the majority of CpARF genes (28) have significantly response to drought stress (P < 0.05), and most of them (26) are significantly down-regulated under severe drought treatment. Furthermore, 32 out of the 34 CpARF genes have significantly response to high, middle, and low salt stress under salt treatment (P < 0.05). Additionally, subcellular localization analysis confirmed that CpARF16 and CpARF32 were all localized to nucleus. Thus, our findings expand the understanding of the function of CpARF genes and provide a basis for further functional studies on CpARF genes in C. paliurus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01474-1.

Comparative genomic profiling of transport inhibitor Response1/Auxin signaling F-box (TIR1/AFB) genes in eight Pyrus genomes revealed the intraspecies diversity and stress responsiveness patterns

In the genomics of plants and the phytoecosystem, Pyrus (pear) is among the most nutritious fruits and contains fiber that has great health benefits to humans. It is mostly cultivated in temperate regions and is one of the most cultivated pome fruits globally. Pears are highly subjected to biotic and abiotic stresses that affect their yield. TIR1/AFB proteins act as auxin co-receptors during the signaling of nuclear auxins and play a primary role in development-related regulatory processes and responses to biotic and abiotic stresses. However, this gene family and its members have not been explored in Pyrus genomes, and understanding these genes will help obtain useful insights into stress tolerance and ultimately help maintain a high yield of pears. This study reports a pangenome-wide investigation of TIR1/AFB genes from eight Pyrus genomes: Cuiguan (Pyrus pyrifolia), Shanxi Duli (P. betulifolia), Zhongai 1 [(P. ussuriensis × communis) × spp.], Nijisseiki (P. pyrifolia), Yunhong No.1 (P. pyrifolia), d'Anjou (P. communis), Bartlett v2.0 (P. communis), and Dangshansuli v.1.1 (P. bretschneideri). These genes were randomly distributed on 17 chromosomes in each genome. Based on phylogenetics, the identified TIR1/AFB genes were divided into six groups. Their gene structure and motif pattern showed the intraspecific structural conservation as well as evolutionary patterns of Pyrus TIR1/AFBs. The expansion of this gene family in Pyrus is mainly caused by segmental duplication; however, a few genes showed tandem duplication. Moreover, positive and negative selection pressure equally directed the gene's duplication process. The GO and PPI analysis showed that Pyrus TIR1/AFB genes are associated with abiotic stress- and development-related signaling pathways. The promoter regions of Pyrus TIR1/AFB genes were enriched in hormone-, light-, development-, and stress-related cis elements. Furthermore, publicly available RNA-seq data analysis showed that DaTIR1/AFBs have varied levels of expression in various tissues and developmental stages, fruit hardening disease conditions, and drought stress conditions. This indicated that DaTIR1/AFB genes might play critical roles in response to biotic and abiotic stresses. The DaTIR1/AFBs have similar protein structures, which show that they are involved in the same function. Hence, this study will broaden our knowledge of the TIR1/AFB gene family in Pyrus, elucidating their contribution to conferring resistance against various environmental stresses, and will also provide valuable insights for future researchers.

Unraveling the role of PlARF2 in regulating deed formancy in Paeonia lactiflora

MAIN CONCLUSION

PlARF2 can positively regulate the seed dormancy in Paeonia lactiflora Pall. and bind the RY cis-element. Auxin, a significant phytohormone influencing seed dormancy, has been demonstrated to be regulated by auxin response factors (ARFs), key transcriptional modulators in the auxin signaling pathway. However, the role of this class of transcription factors (TFs) in perennials with complex seed dormancy mechanisms remains largely unexplored. Here, we cloned and characterized an ARF gene from Paeonia lactiflora, named PlARF2, which exhibited differential expression levels in the seeds during the process of seed dormancy release. The deduced amino acid sequence of PlARF2 had high homology with those of other plants and contained typical conserved Auxin_resp domain of the ARF family. Phylogenetic analysis revealed that PlARF2 was closely related to VvARF3 in Vitis vinifera. The subcellular localization and transcriptional activation assay showed that PlARF2 is a nuclear protein possessing transcriptional activation activity. The expression levels of dormancy-related genes in transgenic callus indicated that PlARF2 was positively correlated with the contents of PlABI3 and PlDOG1. The germination assay showed that PlARF2 promoted seed dormancy. Moreover, TF Centered Yeast one-hybrid assay (TF-Centered Y1H), electrophoretic mobility shift assay (EMSA) and dual-luciferase reporter assay analysis (Dual-Luciferase) provided evidence that PlARF2 can bind to the 'CATGCATG' motif. Collectively, our findings suggest that PlARF2, as TF, could be involved in the regulation of seed dormancy and may act as a repressor of germination.

Genome-wide identification and expression analysis of ARF gene family in embryonic development of Korean pine (Pinus koraiensis)

BACKGROUND

The Auxin Responsive Factor (ARF) family plays a crucial role in mediating auxin signal transduction and is vital for plant growth and development. However, the function of ARF genes in Korean pine (Pinus koraiensis), a conifer species of significant economic value, remains unclear.

RESULTS

This study utilized the whole genome of Korean pine to conduct bioinformatics analysis, resulting in the identification of 13 ARF genes. A phylogenetic analysis revealed that these 13 PkorARF genes can be classified into 4 subfamilies, indicating the presence of conserved structural characteristics within each subfamily. Protein interaction prediction indicated that Pkor01G00962.1 and Pkor07G00704.1 may have a significant role in regulating plant growth and development as core components of the PkorARFs family. Additionally, the analysis of RNA-seq and RT-qPCR expression patterns suggested that PkorARF genes play a crucial role in the development process of Korean pine.

CONCLUSION

Pkor01G00962.1 and Pkor07G00704.1, which are core genes of the PkorARFs family, play a potentially crucial role in regulating the fertilization and developmental process of Korean pine. This study provides a valuable reference for investigating the molecular mechanism of embryonic development in Korean pine and establishes a foundation for cultivating high-quality Korean pine.

An ARF gene mutation creates flint kernel architecture in dent maize

Dent and flint kernel architectures are important characteristics that affect the physical properties of maize kernels and their grain end uses. The genes controlling these traits are unknown, so it is difficult to combine the advantageous kernel traits of both. We found mutation of ARFTF17 in a dent genetic background reduces IAA content in the seed pericarp, creating a flint-like kernel phenotype. ARFTF17 is highly expressed in the pericarp and encodes a protein that interacts with and inhibits MYB40, a transcription factor with the dual functions of repressing PIN1 expression and transactivating genes for flavonoid biosynthesis. Enhanced flavonoid biosynthesis could reduce the metabolic flux responsible for auxin biosynthesis. The decreased IAA content of the dent pericarp appears to reduce cell division and expansion, creating a shorter, denser kernel. Introgression of the ARFTF17 mutation into dent inbreds and hybrids improved their kernel texture, integrity, and desiccation, without affecting yield.

iFLAS: positive-unlabeled learning facilitates full-length transcriptome-based identification and functional exploration of alternatively spliced isoforms in maize

The advent of full-length transcriptome sequencing technologies has accelerated the discovery of novel splicing isoforms. However, existing alternative splicing (AS) tools are either tailored for short-read RNA-Seq data or designed for human and animal studies. The disparities in AS patterns between plants and animals still pose a challenge to the reliable identification and functional exploration of novel isoforms in plants. Here, we developed integrated full-length alternative splicing analysis (iFLAS), a plant-optimized AS toolkit that introduced a semi-supervised machine learning method known as positive-unlabeled (PU) learning to accurately identify novel isoforms. iFLAS also enables the investigation of AS functions from various perspectives, such as differential AS, poly(A) tail length, and allele-specific AS (ASAS) analyses. By applying iFLAS to three full-length transcriptome sequencing datasets, we systematically identified and functionally characterized maize (Zea mays) AS patterns. We found intron retention not only introduces premature termination codons, resulting in lower expression levels of isoforms, but may also regulate the length of 3'UTR and poly(A) tail, thereby affecting the functional differentiation of isoforms. Moreover, we observed distinct ASAS patterns in two genes within heterosis offspring, highlighting their potential value in breeding. These results underscore the broad applicability of iFLAS in plant full-length transcriptome-based AS research.

The ZmbHLH32-ZmIAA9-ZmARF1 module regulates salt tolerance in maize

The growth and productivity of maize (Zea mays), along with other crop plants, can be significantly hindered by salt stress. Nevertheless, the precise molecular mechanism underlying salt tolerance in maize has yet to be fully elucidated. Hence, it was attempted to identify ZmIAA9, a member of the maize Aux/IAA gene family, as a positive regulator of salt tolerance in maize, which was accompanied by the increased ROS detoxification and elevated transcript abundances of ROS scavenging genes. Molecular and biochemical assays have provided compelling evidence that ZmbHLH32, a transcription factor belonging to the bHLH family, was capable of binding directly to the promoter region of ZmIAA9, thereby activating its expression. This interaction between ZmbHLH32 and ZmIAA9 could be critical for the regulation of salt tolerance in maize. As expected, overexpression of ZmbHLH32 led to the enhanced salt tolerance. In contrast, decreased salt tolerance was attained after application of knockout mutants of ZmbHLH32. Furthermore, ZmARF1, which could act as a downstream of ZmIAA9, was found to physically interact with ZmIAA9 and repress the expression levels of ROS scavenging genes. Thus, our work uncovers a novel mechanism of ZmbHLH32-ZmIAA9-ZmARF1 module-mediated salt tolerance in maize, which can be exploited for breeding salt-tolerant maize varieties.

C4 leaf development and evolution

C4 photosynthesis is more efficient than C3 photosynthesis for two reasons. First, C4 plants have evolved efficient C4 enzymes to suppress wasteful photorespiration and enhance CO2 fixation. Second, C4 leaves have Kranz anatomy in which the veins are surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. The BS and M cells are functionally well differentiated and also well coordinated for rapid assimilation of atmospheric CO2 and transport of photo-assimilates between the two types of cells. Recent comparative transcriptomics of developing M and BS cells in young maize embryonic leaves revealed not only potential regulators of BS and M cell differentiation but also rapid early BS cell differentiation whereas slower, more prolonged M cell differentiation, contrary to the traditional view of a far simpler process of M cell development. Moreover, new upstream regulators of Kranz anatomy development have been identified and a number of gene co-expression modules for early vascular development have been inferred. Also, a candidate gene regulatory network associated with Kranz anatomy and vascular development has been constructed. Additionally, how whole genome duplication (WGD) may facilitate C4 evolution has been studied and the reasons for why the same WGD event led to successful C4 evolution in Gynandropsis gynandra but not in the sister species Tarenaya hassleriana have been proposed. Finally, new future research directions are suggested.

Genome-Wide Analysis of the DC1 Domain Protein Gene Family in Tomatoes under Abiotic Stress

DC1 (Divergent C1) domain proteins are a new class of proteins that have been discovered in recent years, which play an important role in plant growth, development, and stress response. In order to better study the distribution and function of DC1 domain proteins in tomatoes, a genome-wide identification was conducted. It was found that there are twenty-one DC1 domain protein genes distributed on nine chromosomes of tomatoes, named SlCHP1-21. Phylogenetic analysis shows that twenty-one SlCHP genes are divided into six subfamilies. Most of the SlCHP genes in tomatoes have no or very short introns. All SlCHP proteins, with the exception of SlCHP8 and SlCHP17, contain variable amounts of C1 domain. Analysis of the SlCHP gene promoter sequence revealed multiple cis-elements responsive to plant stress. qRT-CR analysis showed that most members of SlCHP gene expressed in the roots. The SlCHP11, 13, 16, 17, and SlCHP20 genes showed specific responses to high temperature, low temperature, salt, and drought stress. In addition, the subcellular localization and interaction proteins of SlCHP were analyzed and predicted. Together, these results provides a theoretical basis for further exploration of the function and mechanism of the SlCHP gene in tomatoes.

Identification and Evolutionary Analysis of the Auxin Response Factor (ARF) Family Based on Transcriptome Data from Caucasian Clover and Analysis of Expression Responses to Hormones

Caucasian clover (Trifolium ambiguum M. Bieb.) is an excellent perennial plant in the legume family Fabaceae, with a well-developed rhizome and strong clonal growth. Auxin is one of the most important phytohormones in plants and plays an important role in plant growth and development. Auxin response factor (ARF) can regulate the expression of auxin-responsive genes, thus participating in multiple pathways of auxin transduction signaling in a synergistic manner. No genomic database has been established for Caucasian clover. In this study, 71 TaARF genes were identified through a transcriptomic database of Caucasian clover rhizome development. Phylogenetic analysis grouped the TaARFs into six (1-6) clades. Thirty TaARFs contained a complete ARF structure, including three relatively conserved regions. Physical and chemical property analysis revealed that TaARFs are unstable and hydrophilic proteins. We also analyzed the expression pattern of TaARFs in different tissues (taproot, horizontal rhizome, swelling of taproot, rhizome bud and rhizome bud tip). Quantitative real-time RT-PCR revealed that all TaARFs were responsive to phytohormones (indole-3-acetic acid, gibberellic acid, abscisic acid and methyl jasmonate) in roots, stems and leaves. These results helped elucidate the role of ARFs in responses to different hormone treatments in Caucasian clover.

Identification of ARF transcription factor gene family and its defense responses to bacterial infection and salicylic acid treatment in sugarcane

Auxin response factor (ARF) is a critical regulator in the auxin signaling pathway, involved in a variety of plant biological processes. Here, gene members of 24 SpapARFs and 39 SpnpARFs were identified in two genomes of Saccharum spontaneum clones AP85-441 and Np-X, respectively. Phylogenetic analysis showed that all ARF genes were clustered into four clades, which is identical to those ARF genes in maize (Zea mays) and sorghum (Sorghum bicolor). The gene structure and domain composition of this ARF family are conserved to a large degree across plant species. The SpapARF and SpnpARF genes were unevenly distributed on chromosomes 1-8 and 1-10 in the two genomes of AP85-441 and Np-X, respectively. Segmental duplication events may also contribute to this gene family expansion in S. spontaneum. The post-transcriptional regulation of ARF genes likely involves sugarcane against various stressors through a miRNA-medicated pathway. Expression levels of six representative ShARF genes were analyzed by qRT-PCR assays on two sugarcane cultivars [LCP85-384 (resistant to leaf scald) and ROC20 (susceptible to leaf scald)] triggered by Acidovorax avenae subsp. avenae (Aaa) and Xanthomonas albilineans (Xa) infections and salicylic acid (SA) treatment. ShARF04 functioned as a positive regulator under Xa and Aaa stress, whereas it was a negative regulator under SA treatment. ShARF07/17 genes played positive roles against both pathogenic bacteria and SA stresses. Additionally, ShARF22 was negatively modulated by Xa and Aaa stimuli in both cultivars, particularly LCP85-384. These findings imply that sugarcane ARFs exhibit functional redundancy and divergence against stressful conditions. This work lays the foundation for further research on ARF gene functions in sugarcane against diverse environmental stressors.

Genome-wide investigation of ARF transcription factor gene family and its responses to abiotic stress in Coix (Coix lacryma-jobi L.)

Auxin response factor (ARF) is an important transcription factor that regulates the expression of auxin-responsive genes by direct binding to their promoters, which play a central role in plant growth, development, and response to abiotic stresses. The availability of the entire Coix (Coix lacryma-jobi L.) genome sequence provides an opportunity to investigate the characteristics and evolutionary history of the ARF gene family in this medicine and food homology plant for the first time. In this study, a total of 27 ClARF genes were identified based on the genome-wide sequence of Coix. Twenty-four of the 27 ClARF genes were unevenly distributed on 8 chromosomes except Chr 4 and 10, and the remaining three genes (ClARF25-27) were not assigned to any chromosome. Most of the ClARF proteins were predicted to be localized to the nucleus, except ClARF24, which was localized to both the plasma membrane and nucleus. Twenty-seven ClARFs were clustered into six subgroups based on the phylogenetic analysis. Duplication analysis showed that segmental duplication, rather than tandem duplications promoting the expansion of the ClARF gene family. Synteny analysis showed that purifying selection might have been a primary driving force in the development of the ARF gene family in Coix and other investigated cereal plants. The prediction of the cis element of the promoter showed that 27 ClARF genes contain several stress response elements, suggesting that ClARFs might be involved in the abiotic stress response. Expression profile analysis shows that 27 ClARF genes were all expressed in the root, shoot, leaf, kernel, glume, and male flower of Coix with varying expression levels. Furthermore, qRT-PCR analyses revealed that the majority of ClARFs members were upregulated or downregulated in response to hormone treatment and abiotic stress. The current study expands our understanding of the functional roles of ClARFs in stress responses and provides basic information for the ClARF genes.

Genome-wide identification and expression analysis of the Auxin-Response factor (ARF) gene family in Medicago sativa under abiotic stress

BACKGROUND

Alfalfa (Medicago sativa) is the most widely planted legume forage and one of the most economically valuable crops in the world. The periodic changes in its growth and development and abiotic stress determine its yield and economic benefits. Auxin controls many aspects of alfalfa growth by regulating gene expression, including organ differentiation and stress response. Auxin response factors (ARF) are transcription factors that play an essential role in auxin signal transduction and regulate the expression of auxin-responsive genes. However, the function of ARF transcription factors is unclear in autotetraploid-cultivated alfalfa.

RESULT

A total of 81 ARF were identified in the alfalfa genome in this study. Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed, identifying that ARF genes are mainly involved in transcriptional regulation and plant hormone signal transduction pathways. Phylogenetic analysis showed that MsARF was divided into four clades: I, II, III, and IV, each containing 52, 13, 7, and 9 genes, respectively. The promoter region of the MsARF gene contained stress-related elements, such as ABRE, TC-rich repeats, MBS, LTR. Proteins encoded by 50 ARF genes were localized in the nucleus without guide peptides, signal peptides, or transmembrane structures, indicating that most MsARF genes are not secreted or transported but only function in the nucleus. Protein structure analysis revealed that the secondary and tertiary structures of the 81 MsARF genes varied. Chromosomal localization analysis showed 81 MsARF genes were unevenly distributed on 25 chromosomes, with the highest distribution on chromosome 5. Furthermore, 14 segmental duplications and two sets of tandem repeats were identified. Expression analysis indicated that the MsARF was differentially expressed in different tissues and under various abiotic stressors. The quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that the expression profiles of 23 MsARF genes were specific to abiotic stresses such as drought, salt, high temperature, and low temperature, as well as tissue-specific and closely related to the duration of stress.

CONCLUSION

This study identified MsARF in the cultivated alfalfa genome based on the autotetraploid level, which GO, KEGG analysis, phylogenetic analysis, sequence characteristics, and expression pattern analysis further confirmed. Together, these findings provide clues for further investigation of MsARF functional verification and molecular breeding of alfalfa. This study provides a novel approach to systematically identify and characterize ARF transcription factors in autotetraploid cultivated alfalfa, revealing 23 MsARF genes significantly involved in response to various stresses.

The Knockdown of AUXIN RESPONSE FACTOR 2 Confers Enhanced Tolerance to Salt and Drought Stresses in Tomato (Solanum lycopersicum L.)

Auxin response factors (ARFs) act as key elements of the auxin-signaling pathway and play important roles in the process of a plant's growth, development, and response to environmental conditions. We studied the implication of the SlARF2 gene in the tomato response to salt (150 mM of NaCl) and drought (15% PEG 20000) stresses. The functional characterization of SlARF2 knockdown tomato mutants revealed that the downregulation of this gene enhanced primary root length and root branching and reduced plant wilting. At the physiological level, the arf2 mutant line displayed higher chlorophyll, soluble sugars, proline, and relative water contents as well as lower stomatal conductance and a decreased malondialdehyde content. Moreover, SlARF2 knockdown tomato mutants demonstrated higher activities of the antioxidant enzymes superoxide dismutase (SOD) and catalase (CAT) under salt and drought stresses than the wild type. Indeed, the stress tolerance of the arf2 mutant was also reflected by the upregulation of stress-related genes involved in ROS scavenging and plant defense, including SOD, CAT, dehydration-responsive element-binding protein, and early responsive to dehydration, which can ultimately result in a better resistance to salt and drought stresses. Furthermore, the transcriptional levels of the Δ1-pyrroline-5-carboxylate synthase (P5CS) gene were upregulated in the arf2 mutant after stress, in correlation with the higher levels of proline. Taken together, our findings reveal that SlARF2 is implicated in salt and drought tolerance in tomato and provides some considerable elements for improving the abiotic stress tolerance and increasing the crop yields of tomato.

Genome-Wide Identification of MADS-Box Genes in Taraxacum kok-saghyz and Taraxacum mongolicum: Evolutionary Mechanisms, Conserved Functions and New Functions Related to Natural Rubber Yield Formation

MADS-box transcription regulators play important roles in plant growth and development. However, very few MADS-box genes have been isolated in the genus Taraxacum, which consists of more than 3000 species. To explore their functions in the promising natural rubber (NR)-producing plant Taraxacum kok-saghyz (TKS), MADS-box genes were identified in the genome of TKS and the related species Taraxacum mongolicum (TM; non-NR-producing) via genome-wide screening. In total, 66 TkMADSs and 59 TmMADSs were identified in the TKS and TM genomes, respectively. From diploid TKS to triploid TM, the total number of MADS-box genes did not increase, but expansion occurred in specific subfamilies. Between the two genomes, a total of 11 duplications, which promoted the expansion of MADS-box genes, were identified in the two species. TkMADS and TmMADS were highly conserved, and showed good collinearity. Furthermore, most TkMADS genes exhibiting tissue-specific expression patterns, especially genes associated with the ABCDE model, were preferentially expressed in the flowers, suggesting their conserved and dominant functions in flower development in TKS. Moreover, by comparing the transcriptomes of different TKS lines, we identified 25 TkMADSs related to biomass formation and 4 TkMADSs related to NR content, which represented new targets for improving the NR yield of TKS.

Characterization of myo-inositol oxygenase from rice (OsMIOX): influence of salinity stress in different indica rice cultivars

Myo-inositol oxygenase (MIOX), the only catabolic enzyme of the inositol pathway, catalyzes conversion of myo-inositol to D-GlcA (glucuronic acid). The present study encompasses bioinformatic analysis of MIOX gene across phylogenetically related plant lineages and representative animal groups. Comparative motif analysis of the MIOX gene(s) across various plant groups suggested existence of abiotic- stress related cis-acting elements such as, DRE, MYB, MYC, STRE, MeJa among others. A detailed analysis revealed a single isoform of MIOX gene, located in chromosome 6 of indica rice (Oryza sativa) with an open reading frame of 938 bp coding for 308 amino acids producing a protein of ~ 35 kD. Secondary structure prediction of the protein gave the predicted number of 144 alpha helices and 154 random coils. The three-dimensional structure suggested it to be a monomeric protein with a single domain. Bacterial overexpression of the protein, purification and enzyme assay showed optimal catalytic activity at pH 7.5-8 at an optimal temperature of 37 °C with Michaelis constant of 40.92 mM. The range of Km was determined as 22.74-28.7 mM and the range of Vmax was calculated as 3.51-3.6 µM/min, respectively. Four salt-tolerant and salt-sensitive rice cultivars displayed differential gene expression of OsMIOX at different time points in different tissues under salinity and drought stress as observed from qRT-PCR data, microarray results and protein expression profile in immunoblot analysis. Gel volumetric analysis confirmed a very high expression of MIOX in roots and leaves on 7th day following germination. Microarray data showed high expression of MIOX at all developmental stages including seedling growth and reproduction. These data suggest that OsMIOX might have a role to play in rice abiotic stress responses mediated through the myo-inositol oxidation pathway.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01340-6.

Identification and transcriptome data analysis of ARF family genes in five Orchidaceae species

The Orchidaceae is a large family of perennial herbs especially noted for the exceptional diversity of specialized flowers. Elucidating the genetic regulation of flowering and seed development of orchids is an important research goal with potential utility in orchid breeding programs. Auxin Response Factor (ARF) genes encode auxin-responsive transcription factors, which are involved in the regulation of diverse morphogenetic processes, including flowering and seed development. However, limited information on the ARF gene family in the Orchidaceae is available. In this study, 112 ARF genes were identified in the genomes of 5 orchid species (Apostasia shenzhenica, Dendrobium catenatum, Phalaenopsis aphrodite, Phalaenopsis equestris and Vanilla planifolia,). These genes were grouped into 7 subfamilies based on their phylogenetic relationships. Compared with the ARF family in model plants, such as Arabidopsis thaliana and Oryza sativa, one group of ARF genes involved in pollen wall synthesis has been lost during evolution of the Orchidaceae. This loss corresponds with absence of the exine in the pollinia. Through mining of the published genomic and transcriptomic data for the 5 orchid species: the ARF genes of subfamily 4 may play an important role in flower formation and plant growth, whereas those of subfamily 3 are potentially involved in pollen wall development. the study results provide novel insights into the genetic regulation of unique morphogenetic phenomena of orchids, which lay a foundation for further analysis of the regulatory mechanisms and functions of sexual reproduction-related genes in orchids.

LkARF7 and LkARF19 overexpression promote adventitious root formation in a heterologous poplar model by positively regulating LkBBM1

Cuttage propagation involves adventitious root formation induced by auxin. In our previous study, Larix kaempferi BABY BOOM 1 (LkBBM1), which is known to regulate adventitious root formation, was affected by auxin. However, the relationship between LkBBM1 and auxin remains unclear. Auxin response factors (ARFs) are a class of important transcription factors in the auxin signaling pathway and modulate the expression of early auxin-responsive genes by binding to auxin response elements. In the present study, we identified 14 L. kaempferi ARFs (LkARFs), and found LkARF7 and LkARF19 bound to LkBBM1 promoter and enhanced its transcription using yeast one-hybrid, ChIP-qPCR, and dual-luciferase assays. In addition, the treatment with naphthalene acetic acid promoted the expression of LkARF7 and LkARF19. We also found that overexpression of these two genes in poplar promoted adventitious root formation. Furthermore, LkARF19 interacted with the DEAD-box ATP-dependent RNA helicase 53-like protein to form a heterodimer to regulate adventitious root formation. Altogether, our results reveal an additional regulatory mechanism underlying the control of adventitious root formation by auxin.

Advances in structure and function of auxin response factor in plants

Auxin is a crucial phytohormone that has various effects on the regulators of plant growth and development. Auxin signal transduction is mainly controlled by two gene families: auxin response factor (ARF) and auxin/indole-3-acetic acid (Aux/IAA). ARFs are plant-specific transcription factors that bind directly to auxin response elements in the promoters of auxin-responsive genes. ARF proteins contain three conserved regions: a conserved N-terminal B3 DNA-binding domain, a variable intermediate middle region domain that functions in activation or repression, and a C-terminal domain including the Phox and Bem1p region for dimerization, similar to the III and IV elements of Aux/IAA, which facilitate protein-protein interaction through homodimerization of ARF proteins or heterodimerization of ARF and Aux/IAA proteins. In the two decades following the identification of the first ARF, 23 ARF members have been identified and characterized in Arabidopsis. Using whole-genome sequencing, 22, 25, 23, 25, and 36 ARF genes have been identified in tomato, rice, wheat, sorghum, and maize, respectively, in addition to which the related biofunctions of some ARFs have been reported. ARFs play crucial roles in regulating the growth and development of roots, leaves, flowers, fruits, seeds, responses to biotic and abiotic stresses, and phytohormone signal crosstalk. In this review, we summarize the research progress on the structures and functions of ARFs in Arabidopsis, tomato, and cereal crops, to provide clues for future basic research on phytohormone signaling and the molecular design breeding of crops.

Phylogeny, gene structures, and expression patterns of the auxin response factor (GhARF2) in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Auxin response factors (ARFs) are a class of transcription factors that regulate the expression of auxin-responsive genes and play important functions in plant growth and development. To understand the biological functions of the auxin response factor GhARF2 gene in upland cotton, the coding sequence (CDS) of GhARF2 gene was cloned, and its protein sequence, evolutionary relationship, subcellular localization and expression pattern were analysed.

METHODS

The CDS sequence of GhARF2 gene was cloned from upland cotton variety Baimian No.1, and its protein sequence was analyzed by bioinformatics method. The subcellular localization of GhARF2 protein was detected by tobacco epidermal transient transformation system, and the tissue expression and stress expression pattern of GhARF2 were analyzed by quantitative Real‑Time PCR (qRT-PCR).

RESULTS

The full-length CDS of GhARF2 gene was 2583 bp, encoded 860 amino acids, and had a molecular weight and an isoelectric point of 95.46 KDa and 6.02, respectively. The GhARF2 protein had multiple phosphorylation sites, no transmembrane domain, and secondary structures dominated by random coils and alpha helix. The GhARF2 protein had 3 conserved typical domains of ARF gene family members, including the B3 DNA binding domain, the Auxin_resp domain, and the Aux/IAA domain. Phylogenetic analysis revealed that ARF2 proteins in different species were clustered in the Group A subgroup, in which GhARF2 was closely related to TcARF2 of Theobroma cacao L. (Malvaceae). The subcellular localization results showed that the GhARF2 protein was localized in the nucleus. Analysis of tissue expression pattern showed that the GhARF2 gene was expressed in all tested tissues, with the highest expression levels in sepal, followed by leaf, and the lowest expression levels in fiber. Further stress expression analysis showed that the GhARF2 gene was induced by drought, high-temperature, low-temperature and salt stress, and had different expression patterns under different stress conditions.

CONCLUSION

These results established a foundation for understanding the functions of GhARF2 and breeding varieties with high-stress tolerance in cotton.

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.

PyuARF16/33 Are Involved in the Regulation of Lignin Synthesis and Rapid Growth in Populus yunnanensis

(1) Background: Lignin is a unique component of the secondary cell wall, which provides structural support for perennial woody plants. ARFs are the core factors of the auxin-signaling pathway, which plays an important role in promoting plant growth, but the specific relationship between auxin response factors (ARFs) and lignin has not been fully elucidated with regard to rapid plant growth in forest trees. (2) Objectives: This study aimed to investigate the relationship between ARFs and lignin with regard to rapid plant growth in forest trees. (3) Methods: We used bioinformatics analysis to investigate the PyuARF family, find genes homologous to ARF6 and ARF8 in Populus yunnanensis, and explore the changes in gene expression and lignin content under light treatment. (4) Results: We identified and characterized 35 PyuARFs based on chromosome-level genome data from P. yunnanensis. In total, we identified 92 ARF genes in P. yunnanensis, Arabidopsis thaliana, and Populus trichocarpa, which were subsequently divided into three subgroups based on phylogenetic analysis and classified the conserved exon-intron structures and motif compositions of the ARFs into the same subgroups. Collinearity analysis suggested that segmental duplication and whole-genome duplication events were majorly responsible for the expansion of the PyuARF family, and the analysis of Ka/Ks indicated that the majority of the duplicated PyuARFs underwent purifying selection. The analysis of cis-acting elements showed that PyuARFs were sensitive to light, plant hormones, and stress. We analyzed the tissue-specific transcription profiles of PyuARFs with transcriptional activation function and the transcription profiles of PyuARFs with high expression under light in the stem. We also measured the lignin content under light treatment. The data showed that the lignin content was lower, and the gene transcription profiles were more limited under red light than under white light on days 1, 7, and 14 of the light treatments. The results suggest that PyuARF16/33 may be involved in the regulation of lignin synthesis, thereby promoting the rapid growth of P. yunnanensis. (5) Conclusions: Collectively, this study firstly reports that PyuARF16/33 may be involved in the regulation of lignin synthesis and in promoting the rapid growth in P. yunnanensis.

Genome-wide analysis of the SBT gene family involved in drought tolerance in cotton

The subtilisin-like proteases (SBTs) are a large family of serine peptidases that are unique to plants. Previous studies have shown that SBTs are associated with developmental processes and environmental responses. However, comprehensive identification and systematic analysis of the SBT family have not been conducted in cotton. We used bioinformatics methods to analyze the structural characteristics, phylogenetic relationships, gene structures, expression modes, evolutionary relationships, selection pressures and stress responses of SBT gene family members in upland cotton. In this study, we identified 120 and 112 SBTs in the tetraploid cotton species G. hirsutum and G. barbadense, while 67 and 69 SBTs were identified in the diploid species G. arboreum and G. raimondii, respectively; these SBTs were divided into five distinct subfamilies. We identified the SBT gene GhSBT27A, and explore its function through virus-induced gene silencing and transmission electron microscopy. These results suggested that the GhSBT27A gene was involved in the response to drought stress. These results lay a foundation for further study on the drought stress mechanism of cotton.

Genome-Wide Analysis of Auxin Response Factors in Lettuce (Lactuca sativa L.) Reveals the Positive Roles of LsARF8a in Thermally Induced Bolting

Warm temperatures induce plant bolting accompanied by flower initiation, where endogenous auxin is dynamically associated with accelerated growth. Auxin signaling is primarily regulated by a family of plant-specific transcription factors, AUXIN RESPONSE FACTORS (ARFs), which either activate or repress the expression of downstream genes in response to developmental and environmental cues. However, the relationship between ARFs and bolting has not been completely understood in lettuce yet. Here, we identified 24 LsARFs (Lactuca sativa ARFs) in the lettuce genome. The phylogenetic tree indicated that LsARFs could be classified into three clusters, which was well supported by the analysis of exon-intron structure, consensus motifs, and domain compositions. RNA-Seq analysis revealed that more than half of the LsARFs were ubiquitously expressed in all tissues examined, whereas a small number of LsARFs responded to UV or cadmium stresses. qRT-PCR analysis indicated that the expression of most LsARFs could be activated by more than one phytohormone, underling their key roles as integrative hubs of different phytohormone signaling pathways. Importantly, the majority of LsARFs displayed altered expression profiles under warm temperatures, implying that their functions were tightly associated with thermally accelerated bolting in lettuce. Importantly, we demonstrated that silencing of LsARF8a, expression of which was significantly increased by elevated temperatures, resulted in delayed bolting under warm temperatures, suggesting that LsARF8a might conduce to the thermally induced bolting. Together, our results provide molecular insights into the LsARF gene family in lettuce, which will facilitate the genetic improvement of the lettuce in an era of global warming.

Genome-wide identification of the auxin response factor (ARF) gene family in Magnolia sieboldii and functional analysis of MsARF5

Auxin plays an essential role in flowering, embryonic development, seed dormancy, and germination. Auxin response factors (ARFs) are plant-specific key transcriptional factors in mediating the gene expression network of auxin signaling. Although ARFs in model plants such as Arabidopsis had been well characterized, their identities and potential roles in non-model plants are less studied. Here, we performed genome-wide identification of ARFs in Magnolia sieboldii K. Koch, a primitive species with high taxonomic importance and medicinal values. We found 25 ARF genes in M. sieboldii, which were widely distributed across multiple chromosomes. Based on sequence similarity, the encoded proteins could be either transcriptional repressors or activators. Gene expression analysis showed a dynamic pattern for many ARFs including MsARF5 during seed germination. In addition, overexpressing of MsARF5 showed that it restores many developmental defects in the Arabidopsis mutant. Moreover, two phenotypically distinct transgenic Arabidopsis lines were obtained, indicating a link between gene expression levels and developmental phenotypes. Taken together, we provided a systematic investigation of the ARF gene family in M. sieboldii and revealed an important role of MsARF5 in mediating auxin signaling.

Temporal and spatial auxin responsive networks in maize primary roots

Auxin is a key regulator of root morphogenesis across angiosperms. To better understand auxin-regulated networks underlying maize root development, we have characterized auxin-responsive transcription across two time points (30 and 120 min) and four regions of the primary root: the meristematic zone, elongation zone, cortex and stele. Hundreds of auxin-regulated genes involved in diverse biological processes were quantified in these different root regions. In general, most auxin-regulated genes are region unique and are predominantly observed in differentiated tissues compared with the root meristem. Auxin gene regulatory networks were reconstructed with these data to identify key transcription factors that may underlie auxin responses in maize roots. Additionally, Auxin-Response Factor subnetworks were generated to identify target genes that exhibit tissue or temporal specificity in response to auxin. These networks describe novel molecular connections underlying maize root development and provide a foundation for functional genomic studies in a key crop.

Regulators of early maize leaf development inferred from transcriptomes of laser capture microdissection (LCM)-isolated embryonic leaf cells

The superior photosynthetic efficiency of C₄ leaves over C₃ leaves is owing to their unique Kranz anatomy, in which the vein is surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. Kranz anatomy development starts from three contiguous ground meristem (GM) cells, but its regulators and underlying molecular mechanism are largely unknown. To identify the regulators, we obtained the transcriptomes of 11 maize embryonic leaf cell types from five stages of pre-Kranz cells starting from median GM cells and six stages of pre-M cells starting from undifferentiated cells. Principal component and clustering analyses of transcriptomic data revealed rapid pre-Kranz cell differentiation in the first two stages but slow differentiation in the last three stages, suggesting early Kranz cell fate determination. In contrast, pre-M cells exhibit a more prolonged transcriptional differentiation process. Differential gene expression and coexpression analyses identified gene coexpression modules, one of which included 3 auxin transporter and 18 transcription factor (TF) genes, including known regulators of Kranz anatomy and/or vascular development. In situ hybridization of 11 TF genes validated their expression in early Kranz development. We determined the binding motifs of 15 TFs, predicted TF target gene relationships among the 18 TF and 3 auxin transporter genes, and validated 67 predictions by electrophoresis mobility shift assay. From these data, we constructed a gene regulatory network for Kranz development. Our study sheds light on the regulation of early maize leaf development and provides candidate leaf development regulators for future study.

Identification and characterization of auxin response factor (ARF) family members involved in fig (Ficus carica L.) fruit development

The auxin response factor (ARF) combines with AuxREs cis-acting elements in response to auxin to regulate plant development. To date, no comprehensive analysis of ARF genes expressed during fruit development has been conducted for common fig (Ficus carica L.). In this study, members of the FcARF gene family were screened, identified in the fig genome database and their features characterized using bioinformatics. Twenty FcARF genes were clustered into three classes, with almost similar highly conserved DBD (B3-like DNA binding domain), AUX/IAA (auxin/indole-3-acetic acid gene family) and MR domain structure among class members. Analysis of amino acid species in MR domain revealed 10 potential transcription activators and 10 transcription inhibitors, and 17 FcARF members were predicted to be located in the nucleus. DNA sequence analysis showed that the ARF gene family consisted of 4-25 exons, and the promoter region contained 16 cis-acting elements involved in stress response, hormone response and flavonoid biosynthesis. ARF genes were expressed in most tissues of fig, especially flower and peel. Transcriptomics analysis results showed that FcARF2, FcARF11 and FcARF12, belonging to class-Ia, were stably and highly expressed in the early development stage of flower and peel of 'Purple peel' fig. However, their expression levels decreased after maturity. Expression of class-Ic member FcARF3 conformed to the regularity of fig fruit development. These four potential transcription inhibitors may regulate fruit growth and development of 'Purple Peel' fig. This study provides comprehensive information on the fig ARF gene family, including gene structure, chromosome position, phylogenetic relationship and expression pattern. Our work provides a foundation for further research on auxin-mediated fig fruit development.

Systematic Identification and Expression Analysis of the Auxin Response Factor (ARF) Gene Family in Ginkgo biloba L

Auxin participates in various physiological and molecular response-related developmental processes and is a pivotal hormone that regulates phenotypic formation in plants. Auxin response factors (ARFs) are vital transcription factors that mediate downstream auxin signaling by explicitly binding to auxin-responsive genes' promoters. Here, to investigate the possible developmental regulatory functions of ARFs in Ginkgo biloba, through employing comprehensive bioinformatics, we recognized 15 putative GbARF members. Conserved domains and motifs, gene and protein structure, gene duplication, GO enrichment, transcriptome expression profiles, and qRT-PCR all showed that Group I and III members were highly conserved. Among them, GbARF10b and GbARF10a were revealed as transcriptional activators in the auxin response for the development of Ginkgo male flowers through sequences alignment, cis-elements analysis and GO annotation; the results were corroborated for the treatment of exogenous SA. Moreover, the GbARFs expansion occurred predominantly by segmental duplication, and most GbARFs have undergone purifying selection. The Ka/Ks ratio test identified the functional consistence of GbARF2a and GbARF2c, GbARF10b, and GbARF10a in tissue expression profiles and male flower development. In summary, our study established a new research basis for exploring Ginkgo GbARF members' roles in floral organ development and hormone response.

Genome-wide investigation of histone acetyltransferase gene family and its responses to biotic and abiotic stress in foxtail millet (Setaria italica [L.] P. Beauv)

BACKGROUND

Modification of histone acetylation is a ubiquitous and reversible process in eukaryotes and prokaryotes and plays crucial roles in the regulation of gene expression during plant development and stress responses. Histone acetylation is co-regulated by histone acetyltransferase (HAT) and histone deacetylase (HDAC). HAT plays an essential regulatory role in various growth and development processes by modifying the chromatin structure through interactions with other histone modifications and transcription factors in eukaryotic cells, affecting the transcription of genes. Comprehensive analyses of HAT genes have been performed in Arabidopsis thaliana and Oryza sativa. However, little information is available on the HAT genes in foxtail millet (Setaria italica [L.] P. Beauv).

RESULTS

In this study, 24 HAT genes (SiHATs) were identified and divided into four groups with conserved gene structures via motif composition analysis. Phylogenetic analysis of the genes was performed to predict functional similarities between Arabidopsis thaliana, Oryza sativa, and foxtail millet; 19 and 2 orthologous gene pairs were individually identified. Moreover, all identified HAT gene pairs likely underwent purified selection based on their non-synonymous/synonymous nucleotide substitutions. Using published transcriptome data, we found that SiHAT genes were preferentially expressed in some tissues and organs. Stress responses were also examined, and data showed that SiHAT gene transcription was influenced by drought, salt, low nitrogen, and low phosphorus stress, and that the expression of four SiHATs was altered as a result of infection by Sclerospora graminicola.

CONCLUSIONS

Results indicated that histone acetylation may play an important role in plant growth and development and stress adaptations. These findings suggest that SiHATs play specific roles in the response to abiotic stress and viral infection. This study lays a foundation for further analysis of the biological functions of SiHATs in foxtail millet.

Maize miR167-ARF3/30-polyamine oxidase 1 module-regulated H2O2 production confers resistance to maize chlorotic mottle virus

Maize chlorotic mottle virus (MCMV) is the key pathogen causing maize lethal necrosis (MLN). Due to the sharply increased incidence of MLN in many countries, there is an urgent need to identify resistant lines and uncover the underlying resistance mechanism. Here, we showed that the abundance of maize (Zea mays) microR167 (Zma-miR167) positively modulates the degree of resistance to MCMV. Zma-miR167 directly targets Auxin Response Factor3 (ZmARF3) and ZmARF30, both of which negatively regulate resistance to MCMV. RNA-sequencing coupled with gene expression assays revealed that both ZmARF3 and ZmARF30 directly bind the promoter of Polyamine Oxidase 1 (ZmPAO1) and activate its expression. Knockdown or inhibition of enzymatic activity of ZmPAO1 suppressed MCMV infection. Nevertheless, MCMV-encoded p31 protein directly targets ZmPAO1 and enhances the enzyme activity to counteract Zma-miR167-mediated defense to some degree. We uncovered a role of the Zma-miR167-ZmARF3/30 module for restricting MCMV infection by regulating ZmPAO1 expression, while MCMV employs p31 to counteract this defense.

Genome–Wide Identification and Functional Characterization of Auxin Response Factor (ARF) Genes in Eggplant

Auxin response factors (ARFs) are important plant transcription factors that are differentially expressed in response to auxin and various abiotic stresses. ARFs play important roles in mediating plant growth and stress responses; however, these factors have not been studied in eggplants. In this study, genome–wide identification and the functional analysis of the ARF gene family in eggplants (Solanum melongena L.) were performed. A total of 20 ARF (SmARF) genes were identified and phylogenetically classified into three groups. Our analysis revealed four functional domains and 10 motifs in these proteins. Subcellular localization showed that the SmARFs localized in the nucleus. To investigate the biological functions of the SmARFs under 2,4–D and salt stress treatments, quantitative real–time RT–PCR (qRT–PCR) was conducted. Most SmARF genes exhibited changes in expression in response to 2,4–D treatments in the flowers, especially SmARF4 and 7B. All SmARF genes quickly responded to salt stress, except SmARF17 and 19 in leaves, SmARF1A and 7B in roots, and SmARF2A, SmARF7B, and SmARF16B in stems. These results helped to elucidate the role of ARFs in auxin signaling under 2,4–D and salt stress in eggplants.

Genome-Wide Identification of Auxin Response Factors in Peanut (Arachis hypogaea L.) and Functional Analysis in Root Morphology

Auxin response factors (ARFs) play important roles in plant growth and development; however, research in peanut (Arachis hypogaea L.) is still lacking. Here, 63, 30, and 30 AhARF genes were identified from an allotetraploid peanut cultivar and two diploid ancestors (A. duranensis and A. ipaensis). Phylogenetic tree and gene structure analysis showed that most AhARFs were highly similar to those in the ancestors. By scanning the whole-genome for ARF-recognized cis-elements, we obtained a potential target gene pool of AhARFs, and the further cluster analysis and comparative analysis showed that numerous members were closely related to root development. Furthermore, we comprehensively analyzed the relationship between the root morphology and the expression levels of AhARFs in 11 peanut varieties. The results showed that the expression levels of AhARF14/26/45 were positively correlated with root length, root surface area, and root tip number, suggesting an important regulatory role of these genes in root architecture and potential application values in peanut breeding.

TOPLESS promotes plant immunity by repressing auxin signaling and is targeted by the fungal effector Naked1

In plants, the antagonism between growth and defense is hardwired by hormonal signaling. The perception of pathogen-associated molecular patterns (PAMPs) from invading microorganisms inhibits auxin signaling and plant growth. Conversely, pathogens manipulate auxin signaling to promote disease, but how this hormone inhibits immunity is not fully understood. Ustilago maydis is a maize pathogen that induces auxin signaling in its host. We characterized a U. maydis effector protein, Naked1 (Nkd1), that is translocated into the host nucleus. Through its native ethylene-responsive element binding factor-associated amphiphilic repression (EAR) motif, Nkd1 binds to the transcriptional co-repressors TOPLESS/TOPLESS-related (TPL/TPRs) and prevents the recruitment of a transcriptional repressor involved in hormonal signaling, leading to the de-repression of auxin and jasmonate signaling and thereby promoting susceptibility to (hemi)biotrophic pathogens. A moderate upregulation of auxin signaling inhibits the PAMP-triggered reactive oxygen species (ROS) burst, an early defense response. Thus, our findings establish a clear mechanism for auxin-induced pathogen susceptibility. Engineered Nkd1 variants with increased expression or increased EAR-mediated TPL/TPR binding trigger typical salicylic-acid-mediated defense reactions, leading to pathogen resistance. This implies that moderate binding of Nkd1 to TPL is a result of a balancing evolutionary selection process to enable TPL manipulation while avoiding host recognition.

Maize Transcription Factor ZmARF4 Confers Phosphorus Tolerance by Promoting Root Morphological Development

Plant growth and development are closely related to phosphate (Pi) and auxin. However, data regarding auxin response factors (ARFs) and their response to phosphate in maize are limited. Here, we isolated ZmARF4 in maize and dissected its biological function response to Pi stress. Overexpression of ZmARF4 in Arabidopsis confers tolerance of Pi deficiency with better root morphology than wild-type. Overexpressed ZmARF4 can partially restore the absence of lateral roots in mutant arf7 arf19. The ZmARF4 overexpression promoted Pi remobilization and up-regulated AtRNS1, under Pi limitation while it down-regulated the expression of the anthocyanin biosynthesis genes AtDFR and AtANS. A continuous detection revealed higher activity of promoter in the Pi-tolerant maize P178 line than in the sensitive 9782 line under low-Pi conditions. Meanwhile, GUS activity was specifically detected in new leaves and the stele of roots in transgenic offspring. ZmARF4 was localized to the nucleus and cytoplasm of the mesophyll protoplast and interacted with ZmILL4 and ZmChc5, which mediate lateral root initiation and defense response, respectively. ZmARF4 overexpression also conferred salinity and osmotic stress tolerance in Arabidopsis. Overall, our findings suggest that ZmARF4, a pleiotropic gene, modulates multiple stress signaling pathways, and thus, could be a candidate gene for engineering plants with multiple stress adaptation.

Transcriptomic Variations and Network Hubs Controlling Seed Size and Weight During Maize Seed Development

To elucidate the mechanisms underlying seed development in maize, comprehensive RNA-seq analyses were conducted on Zhengdan1002 (ZD1002), Zhengdan958 (ZD958), and their parental lines during seven seed developmental stages. We found that gene expression levels were largely nonadditive in hybrids and that cis-only or trans × cis pattern played a large role in hybrid gene regulation during seed developmental stage. Weighted gene co-expression network (WGCNA) analysis showed that 36 modules were highly correlated (r = -0.90-0.92, p < 0.05) with kernel weight, length, and width during seed development. Forty-five transcription factors and 38 ribosomal protein genes were identified as major hub genes determining seed size/weight. We also described a network hub, Auxin Response Factor 12 of maize (ZmARF12), a member of a family of transcription factor that mediate gene expression in response to auxin, potentially links auxin signal pathways, cell division, and the size of the seeds. The ZmARF12 mutant exhibited larger seed size and higher grain weight. ZmARF12 transcription was negatively associated with cell division during seed development, which was confirmed by evaluating the yield of protoplasts that isolated from the kernels of the mutant and other inbred lines. Transient knock-down of ZmARF12 in maize plants facilitated cell expansion and division, whereas transient silencing of its potential interactor ZmIAA8 impaired cell division. ZmIAA8 expression was repressed in the ZmARF12 over-expressed protoplasts. The mutant phenotype and the genetics studies presented here illustrated evidence that ZmARF12 is a cell division repressor, and potentially determines the final seed size.

Chinese Cherry (Cerasus pseudocerasus Lindl.) ARF7 Participates in Root Development and Responds to Drought and Low Phosphorus

In this paper, an auxin-responsive transcription factor, CpARF7, was isolated from the roots of Chinese cherry (Cerasus pseudocerasus Lindl. Cv. “Manao Hong”). CpARF7 is highly homologous to AtARF7 or AtARF19 in Arabidopsis, and PavARF1 or PavARF14 in sweet cherry. However, in the phenotype of transgenic tomatoes, the root morphology changed, the main root elongated, and the lateral root increased. Both drought treatment and low-phosphorus conditions can elongate the roots of transgenic tomatoes. In addition, the drought resistance and low-phosphorus tolerance of the transgenic lines are improved, and the POD, SOD, and CAT activities under drought and low-phosphorus environments are increased. There is an effect on the tomato somatotropin suppressor gene, SlIAAs, in which SlIAA1/14/19/29 are up-regulated and SlIAA2/11/12/16 are down-regulated. These results indicate that CpARF7 plays an essential regulatory role in root formation and abiotic stress response, and deepens the understanding of auxin-responsive genes in root growth and abiotic stress.

Genome-Wide Identification and Co-Expression Analysis of ARF and IAA Family Genes in Euscaphis konishii: Potential Regulators of Triterpenoids and Anthocyanin Biosynthesis

Euscaphis konishii is an evergreen plant that is widely planted as an industrial crop in Southern China. It produces red fruits with abundant secondary metabolites, giving E. konishii high medicinal and ornamental value. Auxin signaling mediated by members of the AUXIN RESPONSE FACTOR (ARF) and auxin/indole-3-acetic acid (Aux/IAA) protein families plays important roles during plant growth and development. Aux/IAA and ARF genes have been described in many plants but have not yet been described in E. konishii. In this study, we identified 34 EkIAA and 29 EkARF proteins encoded by the E. konishii genome through database searching using HMMER. We also performed a bioinformatic characterization of EkIAA and EkARF genes, including their phylogenetic relationships, gene structures, chromosomal distribution, and cis-element analysis, as well as conserved motifs in the proteins. Our results suggest that EkIAA and EkARF genes have been relatively conserved over evolutionary history. Furthermore, we conducted expression and co-expression analyses of EkIAA and EkARF genes in leaves, branches, and fruits, which identified a subset of seven EkARF genes as potential regulators of triterpenoids and anthocyanin biosynthesis. RT-qPCR, yeast one-hybrid, and transient expression analyses showed that EkARF5.1 can directly interact with auxin response elements and regulate downstream gene expression. Our results may pave the way to elucidating the function of EkIAA and EkARF gene families in E. konishii, laying a foundation for further research on high-yielding industrial products and E. konishii breeding.