1. Introduction
The ubiquitin-protease system plays a crucial role in protein degradation in eukaryotes. Within the cell, ubiquitin E3 ligase participates in varous of physiological processes by regulating the ubiquitination of regulatory proteins [
1,
2]. The Skp-Cullin-F box (SCF) protein, a subset of ligase E3, is characterized by the involvement of SKP-1 and SKP-1 like proteins in substrate recognition. SKP1 further interacts with specific N-terminal domain proteins of the F-box, ultimately assisting in protein degradation [
3,
4]. The F-box protein family is one of the largest protein families in plants and its encoded genes play criitical roles in plant growth and development. They are also involved in the transduction of plant hormone and light signals [
5]. F-box proteins possess a specific structural pattern consisting of approximately 40-50 amino acids. Notably, this family is distinguished by the presence of at least one F-box structural domain at the N-terminal end and an interacting secondary structure at the C-terminal end of the protein [
6]. These features allow for the categorization of the protein into different subfamilies, each with divergent functions. To date, a large number of F-box protein family members have been identified in various plant species, such as
Arabidopsis thaliana (694), rice (687), maize (359), alfalfa (972), soybean (725), ,tomato (139), eggplant (389), and poplar (337), among others [
7,
8,
9,
10,
11,
12,
13,
14].
The expression pattern and gene function of some F–box protein genes in plant growth and development had been clarified. For instance, in
Arabidopsis thaliana,
MAX2 has been identified as a positive regulator of photomorphogenesis. By constructing a
pps mutant, the expression of
MAX2 was reduced, leading to repression of shoot lateral branching [
15]. Song et al. demonstrated that the development of plant leaf morphology is influenced by the expression of
At1g27340 (
LCR), using a constructed
Arabidopsis miR394 mutant [
16]. Interestingly, abnormal leaf development was observed when this gene was expressed either excessively low or high. Similarly, Marrocco et al. confirmed the crucial role of
EID1 in shaping plant photomorphogenesis by manipulating photopigments in
Arabidopsis thaliana [
17]. Additionally, Qin et al. substantiated that the F-box protein family gene
QDtbn1 acts as a negative regulatory factor, determining the number of tassel branches in maize [
18]. Currently, the availability of the peanut genomes allows for the exploration of various gene families in the entire genome [
19,
20,
21,
22], including the SNARE [
23], SAUR [
24], PIF [
25], and aquaporin [
26] families. However, a comprehensive study on the of F-box protein family in peanut has yet to be conducted.
Peanut is a vital oilseed crop that provides high-quality vegetable oil and proteins, making a significant contribution to economic value. The variations in plant varieties depend on the nutritional and reproductive growth patterns, as well as the differences in the angle of lateral branches and main stem [
27]. Therefore, peanut plants can be classified into four branching habits: erect, bunch, spreading and prostrate [
28,
29]. The plant ideotype plays a crucial role in optimizing the utilization of light energy, which greatly affects yield, cultivation methods, and suitability for machine harvesting [
30]. Therefore, studying plant branching is essential for cultivating optimal crop varieties. Previous research on staple cereal crops has successfully developed high-yield versions by studying genes related to branch development [
31,
32]. Similarly, the growth pattern of lateral branches, particularly the first lateral branch, is an important agronomic trait that determines the plant architecture and production in peanuts. However, the study of branch habit genes in peanuts is still limited. Fortunately, our previous research has identified three candidate genes of the F-box family associated with lateral branch angle (LAB) using genome-wide association study (GWAS) and bulk segregant analysis (BSA) [
33]. Among these genes,
Araip.E64SW could be negatively regulated the lateral branch development, potentially leading to a spreading or prostrate appearance. Therefore, the aim of this study is to investigate the expression patterns of F-box protein genes in various tissues of cultivated peanuts, with a focus on their patterns during lateral branch development.
To comprehensively study F-box protein genes and their potential involvement in lateral branch development, we have identified F-box protein genes using bioinformatics approaches, and integrated the expression data from various tissues of "Tifrunner" cultivar and previous transcriptomic data from different stages of branching development. This process highllights the potential AhFBX genes related to lateral branch development in peanut. The main objectives of this research are: (1) assessing the quantity, physicochemical properties, conserved domain features, gene structure, and evolutionary origins of AhFBXs; (2) explicating the computationally predicted expression patterns of AhFBXs aross different peanut tissues; (3) pinpointing the potential AhFBX candidates related to lateral branch development; and (4) elucidating the expression patterns of these candidate genes in branching development across diverse plant types. With this investigation, our goal is to provide gene resources that enhance our understanding of the molecular mechanisms in which F-box associated genes contribute to peanut branching development and impact on plant architecture. This study holds significant value in the formulation of plant architecture.
4. Discussion
The F-box protein superfamily, primarily found in plants, plays a crucial role in a wide range of physiological and biochemical processes. The heterogeneity of F-box protein structures is predominantly due to variations in their C-terminal domains which are primarily responsible for substrate recognition and binding via interaction with SCF subunits. These proteins partake in different signal transduction pathways, thereby helping regulate vital cellular processes, rendering the F-box proteins one of the largest families of regulatory proteins. The identification of F-box proteins through whole-genome sequencing, however, can vary substantially among species, largely due to differences in chromosome number and evolutionary lineage. Moreover, the subfamily composition within the F-box protein family also exhibit diversity across species. Previous researches reported the discovery of 694 F-box proteins in
Arabidopsis, 687 in rice, 359 in maize, 927 in alfalfa, and 725 in soybean [
5,
6,
7,
8,
9,
10,
11,
46]. However, our study has identified a markedly lower count of F-box proteins, totaling merely 95, a figure that falls drastically short compared to crops like rice, maize, soybean, and alfalfa. A previous investigation on maize revealed the highest number of F-box protein subfamilies with a total of 12 subfamilies, which include FBX, FBT (Tub), FBD, FBW (WD40), FBL (LRR), FBK (Kelch-type), FBP (PP2), FBDUF (DUF295), TLH, FBA, JmjC, and ALH [
9]. Alfalfa, on the other hand, boasts the highest number of F-box protein subfamilies among legumes, with 15 subfamilies. These comprised of F-box, FBA, LRR, FBD, Kelch, DUF, PP2, TUB, WD40, PAS, Actin, GSH-synth ATP, ARM, JmjC, and LysM [
10]. In our study, we observed the presence of merely eight subfamilies within the F-box protein family in peanut, a substantially smaller quantity than in maize and alfalfa. We speculate that the restricted count of F-box protein family members and subfamilies in peanut may be accredited to the significant homology between the A and B subgenomes of cultivated peanut, the abundance of repetitive sequences, and the diminished genome heterozygosity [
19].
Distinct subfamilies demonstrating individual domains illustrate diverse functions integral to plant growth and development. Prior research predominantly delved into the role of the F-box protein family members in stress response processes, while investigations related to growth and development have not been extensively conducted. Consequently, it is essential to examine the equilibrium between target genes involved in stress response and those related to growth and development. For instance, a study on the FBK subfamily revealed that OsFBK12, an F-box protein encompassing the Kelch repeats domain in rice, collaborates with OSK1 to form the SCF complex. This results in the degradation of SAMs impacting the ETH content and influencing the leaf senescence process in rice [
47]. Similarly, in
Arabidopsis, the Kelch repeats F-box (KFB) proteins KFB01, KFB20, and KFB50 interact with PAL isoenzymes. The regulation, either through an increase or decrease of KFB expression, impacted the activity of PAL, which subsequently affected the rate-limiting process in the phenylpropanoid pathway. This eventually determined the lignin content in cells [
48]. The Tub domain, earmarked as a distinct feature of the FBT subfamily, is recognized in numerous plant proteins. Overexpressed
AtTLP9 mutants in
Arabidopsis emonstrated ABA sensitivity, hinting at a potential involvement of the Tub domain in the ABA signaling pathway [
49]. The PP2 (PHLOEM PROTEIN 2) domain is a characteristic recognition domain of the FBP subfamily. PP2 proteins are involved in vascular formation, transport of plant nutrients and macromolecules, and signal transduction processes, and are closely related to nutrient transport. PP2 proteins are involved in wound healing, resistance to biotic stress, and nutrient transport in plants. It is speculated that the FBP subfamily could be involved in regulating the formation of phloem in plants and participating in nutrient transport pathways [
50,
51]. The LRR (leucine-rich repeats) domain has been confirmed to be involved in plant root development and immunity in
Arabidopsis thaliana [
52]. In our study, FBL subfamily contained the LRR domain, and the molecular mechanisms underlying their regulation of lateral branch development is not yet understood. The FBO, FBD, FBA, FBU subfamilies identified in this study also have unclear functions in other crops, and other subfamilies, apart from the already known functions, may also participate in regulating lateral branch development.
Whole-genome duplication (WGD) events can provide an explanation for the large number of members within the F-box protein family and their diverse range of functions. Additionally, collinearity and the estimation of
Ka/Ks ratios can assist researchers in indirectly confirming the allotetraploid origin of the peanut genome.By analyzing the Ka/Ks ratios of various F-box genes, it can be inferred that the cultivated peanut has undergone two significant WGD events throughout its evolutionary history. The first event occurred around 60 MYA when legume crops diversified into different species [
53]. The second event took place approximately 2.16 MYA, resulting in the formation of the current allotetraploid peanut [
19].In this study, 19 gene pairs were identified through intraspecific collinearity analysis. Through
Ka/Ks analysis of these gene pairs, it was determined that the duplication timeframe of family members ranged from 13.4 MYA to 0.34 MYA. Among the selected gene pairs, 11 pairs demonstrated differentiation occurring after 2.16 MYA, while 8 pairs indicated differentiation occurring before 2.16 MYA. These findings suggest that the cultivated peanut differentiation also experienced the WGD event. Thus, the F-box protein family members have undergone both the differentiation of the legume crop genome and the duplication event within the peanut gene during the process of evolution.
The F-box protein family displays a wide range of structural and functional diversity, which is evident in the diversity of their expression patterns. Our investigation of
AhFBX the sequences located 2,000 bp upstream of the CDS of cultivated peanuts revealed significant variations, indicating their diverse functionalities. We analyzed the promoter of the
AhFBXs genes and found numerous
cis-elements, including ones associated with light, plant hormone, and stress responses. This suggests that
AhFBX genes play roles in various physiological and biochemical processes. Some commonly found
cis-elements include ARE, G-box, AuxRE, MYB, BOX, and TCA-element. Of particular interest is the BOX element, a light-responsive element composed of TAATTA. The promoter region of the rbcS-3A gene in pea contains several BOX elements, which contribute to the regulation of the light-controlled molecular switch [
54]. In our study, we identified a total of 39 BOX elements in the promoter region of
AhFBX genes associated with lateral branch development. This suggested that peanut may regulate lateral branch development through a light-controlled molecular switch. Another important
cis-regulatory element is the G-box element, which is widely involved in light responses and has a core sequence of CCACGTGG. The G-box binding protein GBF in the rbcS-1A gene of
Arabidopsis thaliana regulates photomorphogenesis [
55]. Additionally, the G-box element in
Arabidopsis thaliana participates in the regulation of circadian rhythm by binding to the pseudo response regulator PRR [
56].
AhFBXs were a multifunctional gene family, which was confirmed via the RNA-seq data and qRT-PCR results in this study. In this study, we identified 9 DEGs for RNA-Seq. Through phylogenetic analysis, we classified these DEGs into four groups based on their protein sequences in
Arabidopsis and peanuts. Group I consisted of
AhFBX46,
AhFBX91,
AtFBX12 (At2g24250), and
AtFBX56 (At2g16220). Group III consisted of
AhFBX40,
AhFBX45,
AhFBX55,
AhFBX90,
AtFBX58 (At3g50080), and
AtFBX5 (At1g23390). Group IV consisted of
AhFBX8,
AhFBX31,
AhFBX57, and
AtFBX79(At4g35930).
AhFBX8,
AhFBX57,
AhFBX31 and
AtFBX79 were divided into Group IV. In previous studies have linked the FBS mutant
fbs4-1 of
Arabidopsis (At4g35930) to stomata development [
57]. Considering the transcriptome data, the specific expression of AhFBX46 and AhFBX91 in the seeds of peanuts, suggested that their involve in seed growth and development. A similar pattern has been observed for
AtFBX12 (At2g24250) in
Arabidopsis, which exhibits specific expression during seed development and was associated with the mitochondrial protein pathway [
58]. Furthermore,
AhFBX40,
AhFBX45,
AhFBX55,
AhFBX90,
AtFBX58 (At3g50080), and
AtFBX5 (At1g23390) were divided into Group III, and potential involved in root development and secondary metabolism. Previous studies showed that
AtFBX58 (At3g50080) in
Arabidopsis regulates lateral root formation [
59], and
AtFBX5 (At1g23390) was a negative regulator involved in the synthesis of brassinosteroids and flavonoids, as well as responsive to leaf, stem, pod development, and environmental stresses [
60]. These findings implied that these genes in cultivated peanuts may be had similar functions. In our study, 9 DEGs displayed different expression pattern in developmental changes of lateral branches, 15DAP-25DAP had significant change trend between the two cultivars. In this study, we speculated that the period of 15-25 DAP may be a critical stage for plant type differentiation in peanut. These DEGs could be provided a basis for future research on the gene function of lateral branch development and differentiation.