QTL Analysis of Five Silique-Related Traits in Brassica napus L. Across Multiple Environments

As an important physiological and reproductive organ, the silique is a determining factor of seed yield and a breeding target trait in rapeseed (Brassica napus L.). Genetic studies of silique-related traits are helpful for rapeseed marker-assisted high-yield breeding. In this study, a recombinant inbred population containing 189 lines was used to perform a quantitative trait loci (QTLs) analysis for five silique-related traits in seven different environments. As a result, 120 consensus QTLs related to five silique-related traits were identified, including 23 for silique length, 25 for silique breadth, 29 for silique thickness, 22 for seed number per silique and 21 for silique volume, which covered all the chromosomes, except C5. Among them, 13 consensus QTLs, one, five, two, four and one for silique length, silique breadth, silique thickness, seed number per silique and silique volume, respectively, were repeatedly detected in multiple environments and explained 4.38–13.0% of the phenotypic variation. On the basis of the functional annotations of Arabidopsis homologous genes and previously reported silique-related genes, 12 potential candidate genes underlying these 13 QTLs were screened and found to be stable in multiple environments by analyzing the re-sequencing results of the two parental lines. These findings provide new insights into the gene networks affecting silique-related traits at the QTL level in rapeseed.


INTRODUCTION
Rapeseed (Brassica napus, AACC, 2n = 38), a major oil crops worldwide, evolved from the double diploidization of Brassica rapa (AA,2n = 20) and Brassica oleracea (CC, 2n = 18) through interspecific hybridization (Murphy, 1999). B. napus is an oil crop with the highest oil production efficiency and is currently an important source of edible vegetable oil in China (Wang, 2010). Because of the high protein content of rapeseed meal, rapeseed is also a high-quality source of animal feed. Additionally, rape straw undergoes an efficient biomass digestion and can be used as a raw material for bio-energy production .
The global population, and thus, food demand, continues to grow; however, the amount of cultivated land is decreasing owing to a variety of factors, such as human activities and climate change. Therefore, increasing the production of rapeseed is currently an urgent requirement and major goal of global rapeseed producers. Yield is an extremely complex trait, and the product of a series of developmental and physiological processes (Van Camp, 2005). Rapeseed yield is determined by three constituent factors: seed weight, silique number and seed number per silique (SPS) (Clarke and Simpson, 1978). The silique, as an important photosynthetic organ of rapeseed, is closely related to the final grain yield (Allen et al., 1971). Silique-related traits, such as silique length (SL), silique breadth (SB), silique thickness (ST) and silique volume (SV), affect the appearance and morphology of silique, and ultimately affect the production and yield of photosynthetic substances (Ferrándiz et al., 1999). Thus, siliques and their related traits are considered the major contributing factors for increasing rapeseed yield. Understanding their genetic bases is of great significance for breeding high-yield rapeseed.
The silique-related traits of rapeseed are all complex traits controlled by quantitative trait loci (QTLs), which are also easily affected by the environment. QTL mapping has been successfully applied to the genetic analyses of the quantitative traits of a variety of crops. In rapeseed, a large number of QTLs, such as those affecting plant height Dong et al., 2021), branch angle Wang et al., 2019) and flowering time (Yu et al., 2019;Xu et al., 2021), have been genetically mapped. In siliques, SL and SPS are the most important traits. At present, more than 100 QTLs for SL, along with more than 200 QTLs for SPS, have been detected in different mapping populations of B. napus, and QTLs controlling SL and SPS are distributed on all 19 chromosomes, explaining 1.9-65.6% and 0.78-57.8% of the phenotypic variance (PV), respectively (Chen et al., 2007(Chen et al., , 2011Radoev et al., 2008;Shi et al., 2009;Zhang et al., 2011;Yang et al., 2012Yang et al., , 2016Yang et al., , 2017Cai et al., 2014Cai et al., , 2016Qi et al., 2014;Fu et al., 2015;Wang X. et al., 2016;Luo et al., 2017;Zhu et al., 2020). Recently,  successfully isolated qSS.C9, the main locus controlling SPS in B. napus that had been discovered by Zhang et al. (2012). qSS.C9 encodes the predicted protein BnaC9.SMG7b. In Arabidopsis, SMG7 controls female fertility and then SPS by regulating the progression from anaphase to telophase in the second meiotic division (Riehs et al., 2008). As in Arabidopsis, BnaC9.SMG7b determines the formation of functional female gametophytes (FGs) by influencing the progression through meiotic anaphase II. Plants with BnaC9.SMG7b deletions exhibit reduced SPS and SL values, which is caused by the reduced ovule number and shorter siliques because of defects in the development of functional FGs. In addition, two major QTL for SL have been successfully cloned on the A9 chromosome of B. napus Shi et al., 2019). The target genes of these two QTLs were identified as BnaA9.ARF18 and BnaA9.CYP78A9 through a finemapping analysis. BnaA9.ARF18 is a homolog of the Arabidopsis auxin response factor 18 (ARF18) that is expressed differentially in various tissues, including root, leaf, stem, bud and ovule, but mainly in the silique wall. A transcription analysis has shown that ARF18 regulates cell elongation in the silique wall and then SL by acting through the auxin-response pathway . BnaA9.CYP78A9 is an ortholog of Arabidopsis CYP78A9, which regulates reproductive development and floral organ size (Sotelo-Silveira et al., 2013). BnaA9.CYP78A9 affects SL in rapeseed by promoting cell elongation in silique valves during silique growth and development (Shi et al., 2019). Although a few silique genes, such as BnaC9.SMG7b, BnaA9.ARF18 and BnaA9.CYP78A9, have been analyzed in B. napus, the genetic mechanisms behind silique-related traits, especially those other than SPS and SL, are far from understood.
The formation of the SPS is closely related to other silique characteristics (SL, SB, ST, and SV), and SPS is significantly positively correlated with SL Wang X. et al., 2016;Yang et al., 2017). As an important silique-related morphological trait, SL also plays a vital role in seed yield, and it has long been used as an indirect selective indicator in breeding for improved seed yields in B. napus. To date, there are extremely limited QTL mapping studies of SB, ST, and SV in B. napus , and more research is needed to determine the relationships among them and their effects on the yield, which is very important for increasing the efficiency levels of breeding programs.
In the present study, a recombinant inbred line (RIL) population containing 189 lines was used to investigate the QTLs for five silique-related traits in seven environments. The objectives were to identify (i) QTLs associated with SL, SPS, SB, ST, and SV across multiple environments and (ii) candidate genes underlying the QTLs that could be stably detected based on the re-sequencing information of two parental lines. This study will lay a good foundation for studying the molecular mechanisms of silique-related traits and increasing yield through rapeseed breeding.

Plant Materials
An RIL population containing 189 lines was constructed by a cross between "APL01" (female parent) and "Holly" (male parent), and it was named the AH population (Wang et al., 2015). The two parental lines showed diversity in silique-related traits. The AH population was previously used for high-density SNP genetic map construction (Wang et al., 2015). The map contains 2,755 SNP-bins, including 11,458 SNP markers and 57 simple sequence repeats, spanning a genetic distance of 2,027.53 cM, with an average distance of 0.72 cM between markers. This population was previously used for QTL mapping and the analysis of apetalous characteristics (Wang et al., 2015), seed fatty acid composition  and seed-related traits . In this study, it was used for the QTL mapping of five silique-related traits.

Field Experiment and Trait Measurement
The AH population, and the parental lines, were grown in four locations in China in 2014-2016. Yangling, Shaanxi Province was planted in September 2015, harvested in May 2016 and recorded as 15YL; Sunan, Gansu Province, was planted in April and harvested in September 2016 and recorded as 16GS; Dali, Shaanxi Province was planted in September 2014 and 2016, harvested in May 2015 and 2017, and recorded as 14DL and 16DL, respectively; and Nanjing, Jiangsu Province was planted in September 2014, 2015 and 2016, harvested in the May of the following years and recorded as 14NJ, 15NJ, and 16NJ, respectively. The field planting of the parents and AH population followed a randomized complete block design with two replications. Each repetition contained two rows with an average spacing of 40 cm between rows and 20 cm between individual plants. The planting, management and harvesting of field materials followed local field-breeding practices.
At the mature stage, three plants growing uniformly from each line were selected, and three well-developed siliques of the first branch adjacent to the main inflorescence from each plant were collected to analyze the silique-related traits . Here, five silique-related traits were investigated: (1) SL, the average length of nine sampled siliques; (2) SB, the widest horizontal length of the space occupied by the silique; (3) ST, the thickness at the same position as SB; (4) SV, the volume measured using the drainage method, which takes the difference between the volume reading of the water after silique immersion and the volume reading of the water prior to the immersion; and (5) SPS, the average seed number of nine sampled siliques. A Vernier caliper was used to measure SL, SB and ST, and a volumetric cylinder was used to measure SV.

Statistical Analysis
SPSS 22.0 software (SPSS Inc., Chicago, IL, United SA) was used for descriptive statistical and correlation analyses of each trait in each environment, and the mean value of the seven environments for each trait was used to carry out correlation analyses. The coefficient of variation was calculated as σ/µ, where σ represents the standard deviation and µ represents the average. The broadsense heritability (h 2 ) was obtained using the R package lme4 (Merk et al., 2012). The computational formula for broad-sense heritability was as follows: h 2 = σ 2 G /(σ 2 G + σ 2 GE /e + σ 2 e /re), where σ 2 G , σ 2 GE , and σ 2 e , respectively, represent genotypic variance, the interaction variance of genotype-by-environment and error variance, n represents the number of environments and r represents the number of replications (Sun et al., 2016).

QTL Mapping and Integration
The composite interval mapping procedure of the Windows QTL Cartographer 2.5 was used to detect the QTLs associated with the five silique-related traits (Wang et al., 2012). The window size was set to 10 cM, a walking speed of 2 cM was selected, and five markers were set as background cofactors. The permutation analysis with 1,000 repetitions (P = 0.05) was performed, and a LOD threshold (2.5-2.8) was used to determine the existence of a QTL and named as "identified QTL." The QTL intervals were established by 2-LOD as approximately 95% QTL confidence intervals (CIs), which is automatically generated by Windows QTL Cartographer 2.5. Identified QTLs for each trait have undergone two rounds of meta-analysis using BioMercator 2.1 software, which can be used to decide the best fitting QTL in accordance with the Akaike criterion (Arcade et al., 2004). The QTL meta-analysis was performed following the description of Wang X. et al. (2016). In the first round, identified QTLs for specific traits with overlapping CIs detected in two or more environments were integrated and named as "consensus QTLs." In the second round, consensus QTLs for different traits located in the same chromosomal region were integrated and named as "unique QTLs" . The designation of the QTL refers to the method of McCouch et al. (1997) with modifications. The name of identified QTL starts with its abbreviation "iq," followed by the environment/trait abbreviation and linkage group (A1-10 and C1-9). When there are multiple QTLs within the same linkage group, the identified QTLs are numbered in accordance with their physical locations. Similarly, the names of consensus QTLs and unique QTLs start with their abbreviations "cq" and "uq," respectively.

Prediction of Underlying Candidate Genes for SL, SB, ST, SPS, and SV
The Illumina HiSeq 2500 platform (Illumina, Inc., San Diego, CA, United States) was used to re-sequence the parental APL01 and Holly genomes with a sequencing depth of 30 × coverage in our previous study (Yu et al., 2021). A total of 74.89 G bp clean data were obtained, with Q30 reaching an average of 92.0%. For candidate gene prediction the following steps were taken: First, the SNP/Indel loci with differences and homozygous between the parents were detected. The clean reads of APL01 and Holly were aligned against B. napus reference genome "Darmor-bzh" (Chalhoub et al., 2014) using BWA software (Li and Durbin, 2009). The SAMtools command "mpileup" was used to genotypes calling ). The locus were screened and retained with criteria: (1) Depth ≥ 6; (2) Quality ≥ 20; (3) Allele mutation ratio of homozygous loci <0.1 or >0.9. After filtration, the genotypes of the two parents were combined by the Bcftools command "merge" (Danecek and McCarthy, 2017), and the SNP/Indel loci with differences and homozygous between the parents were retained; Second, the probe sequences of SNPs on the genetic map were used to map markers on both sides of the QTL interval to the physical position of the reference genome "Darmor-bzh" using BLAST (E value ≤ 1e-10) (Chalhoub et al., 2014), and the variation loci within the QTL interval were found using the results of the first step; and Third, on the basis of the functional annotations of Arabidopsis homologous genes and previously reported silique-related genes, candidate genes for silique-related traits of rapeseed were screened from the mutant loci of these QTLs.

Phenotypic Analysis of Five Silique-Related Traits
The phenotypic performance and h 2 estimates for the five siliquerelated traits in the two parental lines and AH population are presented in Table 1. The five silique-related traits of the two parents, "APL01" and "Holly, " differed significantly in most of the investigated environments, and only two environments, 16DL and 16NJ, showed no significant differences in ST. Compared with "Holly, " "APL01" had significantly greater SB, ST, SV, A wide range of variation, as well as transgressive segregation, were observed for the five silique-related traits, suggesting that alleles with positive effects were distributed in both parents (Table 1 and Figure 1). In addition, the distributions of the five traits were continuous, and most skewness and kurtosis values for the distributions of these traits were <1.0 (Table 1 and Figure 1), which are characteristic of the normal distribution model, indicating that the AH population is suitable for QTL mapping.

Correlation Analysis Among Five Silique-Related Traits
The correlation coefficients among the five silique-related traits were calculated ( Table 2), and most were extremely significantly correlated with each other. For instance, highly

QTL Mapping of Five Silique-Related Traits
Genome-wide QTL analyses were performed separately using the phenotypic data of SL, SB, ST, SPS, and SV grown in seven environments, and a total of 137 identified QTLs distributed across 18 chromosomes were obtained, which explained 4.23-15.31% of the PV (Supplementary Additional  File 1). Among the QTLs, 24 were identified for SL, 33 for SB, 31 for ST, 27 for SPS, and 22 for SV ( Supplementary  Additional File 1). For each trait, QTLs with overlapping CIs detected in different environments were integrated into one QTL through a meta-analysis. As a result, 137 identified QTLs were integrated into 120 consensus QTLs (Figure 2, Supplementary Additional Files 1, 2). Among them, 107, 10, 2 and 1 consensus QTL were detected in 1, 2, 3, and 4 environments, respectively. The consensus QTLs repeatedly detected in multiple environments are listed in Table 3.
For the 33 identified QTLs of SB, each QTL accounted for 4.23-15.3% of the PV, and these QTLs were mapped onto 12 chromosomes. Finally, the 33 identified QTLs were integrated into 25 consensus QTLs, including 3, 1 and 1 QTLs stably expressed in 2, 3 and 4 environments, respectively ( Table 3). The stable QTL cqSB.C6-1 showed a larger effect in the DL than in the NJ environment, with PVs of 12.96 and 9.53% in 14DL and 16DL, respectively but PVs of 4.38 and 5.18% in 14NJ and 15NJ, respectively. The QTL cqSB.C6-1 was repeatedly detected in four environments and may be used in markerassisted selection. Four of the five consensus QTLs (cqSB.C1, cqSB.C7-2, cqSB.C7-4, and cqSB.C6-1) repeatedly detected in different environments had positive additive effects, indicating that the alleles responsible for increasing SB were inherited from female parent "APL01." For the 31 identified QTLs of ST, each QTL explained 4.5-9.5% of the PV and were located on ten chromosomes. The 31 QTLs were integrated into 29 consensus QTLs, including two QTLs that were repeatedly detected in two environments. QTL cqST.A5-2 was integrated from iq16DL.A5 and iq16GS.A5, whereas cqST.A5-3 was integrated from iq14DL.A5-1 and iq14NJ.A5-2 ( Table 3). These two consensus QTLs had negative additive effects, indicating that the favorable alleles were derived from the male parent "Holly." For SPS, 27 identified QTLs were obtained, explaining 4.5-13.8% of the PV. These QTLs were located on nine chromosomes. As a result of meta-QTL analyses, 22 consensus QTLs were obtained, including 18, 3 (cqSPS.A7-3, cqSPS.C3-2, and cqSPS.C6-2) and 1 (cqSPS.C3-3) QTL that were expressed in 1, 2, and 3 environments, respectively ( Table 3). Among the 22 consensus QTLs, only four had negative additive effects and exhibited minor effects, whereas the other 18 QTLs had positive additive effects, indicating that the female parent "APL01" contributed favorable alleles.
For SV, 22 identified QTLs were detected, with each QTL explaining 4.62-11.9% of the PV. These QTLs were distributed across 13 chromosomes. There are two identified QTLs, iq14DL.A10-1 and iq15NJ.A10, with overlapping CIs that were integrated into one consensus QTL, cqSV.A10-1 (Table 3), and the other QTLs were only specifically expressed in a single environment (Supplementary Additional File 1). Of the 21 consensus QTLs, 15 had positive additive effects, including the two QTLs, cqSV.C7-2 and cqSV.C7-3, which explained 10.3 and 11.9%, respectively, of the PV in 16GS, indicating that the female parent "APL01" contributed favorable alleles.

Pleiotropic Unique QTLs for the Five Silique-Related Traits
As shown in Table 2, most of the five silique-related traits showed highly significant correlations. The genetic correlations may be caused by pleiotropism, in which, a single QTL affects the phenotypic variation of multiple traits. Consensus QTLs for the five silique-related traits with overlapping CIs were integrated into unique QTLs using a meta-analysis. In total, 120 consensus QTLs were integrated into 89 unique QTLs (Supplementary  Additional File 3), including 25 QTLs that had pleiotropic effects on two to three traits ( Table 4). These pleiotropic QTLs can partly explain the significant correlations between silique-related traits. For instance, SB had a significant positive correlation with SV, with the highest positive correlation coefficient (0.718) among the five silique-related traits ( Table 2). Five unique QTLs simultaneously controlled SB and SV, in which uqC6-4, uqC7-7, uqC7-8, and uqC8 had positive additive effects, while uqA10-8 had a negative effect on both SB and SV (Table 4).  Additional File 4). We analyzed these 2-271 genes one by one based on the functional annotations of Arabidopsis homologous genes and previously reported siliquerelated genes. Twelve genes underlying 12 consensus except for cqSV.A10-1 were found to have known functions related to the silique-related traits. The CI of cqSL.A10-3 contained the putative gene SMG7 (BnaA10g15730D), which have been successively cloned by . The CI of cqSB.A10-2 contained the putative genes LNG1 (BnaA10g18650D) and CNGC18 (BnaA10g19030D), lng1 dominant mutant displayed longitudinally elongated and transversely narrowed in cells of the siliques (Lee et al., 2006), the CNGC18 point mutations resulted in shorter siliques, reduced male fertility and fewer seeds per silique (Gao Q. et al., 2016). The CI of cqSB.C1 contained the putative gene STY2 (BnaC01g02360D), the sty2 mutation prevents siliques from elongating and showed short silique in Arabidopsis (Kuusk et al., 2002). The CI of cqSB.C7-2 contained the putative gene KINβγ (BnaC07g33610D), its mutant displayed defects in organogenesis and growth, including shorter stature and silique . The CI of cqSB.C7-4 contained the putative gene ABCC13 (BnaC07g35090D), ABCC13 loss function in Arabidopsis leads to decreased silique length and seed yield (El Guizani et al., 2014). The CI of cqST.A5-2 contained the putative gene ARID1 (BnaA05g05310D), mutation in ARID1 showed reduced seed set and short siliques (Zheng et al., 2014). The CI of cqST.A5-3 contained the putative gene LAC4 (BnaA05g06610D), the overexpression of miR397b in Arabidopsis showed increased silique length and seed yield via modulating gene LAC4 . The CI of cqSPS.A7-3 contained the putative gene HTH (BnaA07g22900D), hth mutant showed reduced seed set and short siliques (Krolikowski et al., 2003). The CI of cqSPS.C3-2 contained the putative gene AOG1 (BnaC03g12770D), the siliques of aog1 mutant were much smaller than those of wild type and had a dramatically reduced seed set . The CI of cqSPS.C3-3 contained the putative gene PYL8 (BnaC03g15210D), the pyl8 T-DNA mutant had an extremely low seed yield due to fewer siliques with lower length (Gonzalez-Guzman et al., 2012). In this study, no candidate genes were discovered in the CI of cqSV.A10-1. The CI of cqSB.C6-1 and cqSPS.C6-2 contained the putative gene HMG1 (BnaC06g26470D), T-DNA insertion hmg1 mutant showed shorter siliques and fewer seeds per silique in Arabidopsis (Suzuki et al., 2004).

DISCUSSION
Most agricultural production-related traits are quantitative traits controlled by multiple genes that are also sensitive to the environmental conditions (Abe et al., 2012). Therefore, QTL detection across different environments is necessary for the genetic dissection of complex traits in crops. In this study, phenotypic observations of five silique-related traits replicated in seven environments were used for the QTL analysis, and they might improve the accuracy of QTL detection.
In the present study, most of the five silique-related traits have significant positive correlations with each other ( Table 2). The significant correlations can be explained by unique pleiotropic QTLs. For example, SB had significant positive correlations with SV (coefficient 0.718), whereas, five unique QTLs simultaneously effect SB and SV with the same direction of additive effects (Supplementary Additional File 3). ST showed a negative but low correlations with SL (coefficient -0.026) and SPS (coefficient −0.027) ( Table 4). Accordingly, no unique QTL had effects on both ST and SL, whereas only two unique QTLs controlled ST and SPS. These were uqA4-2, with the same direction of additive effects, and uqA5-4, with the opposite direction of additive effects. In addition, SPS had significant positive correlations with SL, SB and SV, but not ST, which is consistent with previous results Wang X. et al., 2016;Yang et al., 2017). In this study, 25, 29, and 21 consensus QTLs were obtained for SB, ST and SV, respectively. Among these QTLs, five for SB, two for ST and one for SV were repeatedly identified in two to four different environments. To our knowledge, QTLs for SB, ST, and SV in B. napus have rarely been reported . In this study, the QTLs for SB were first mapped on chromosomes A1, A10, C1, C3, C4, C8, and C9; QTLs for ST were first mapped on chromosomes A2, A4, C3, and C7; and QTLs for SV were first mapped on chromosomes A2, A4, A7, A8, A10, C2, C3, C4, C7, C8, and C9. Meanwhile, in previously published results, two QTLs related to SB located at 60.47-62.39 Mb and 68.69-72.5 Mb of C6 had overlapping CIs with the results of our study (cqSB.C6-2, 60.6-61.9 Mb and cqSB.C6-3, 70.1-83.3 Mb), two QTLs related to ST located at 61.3-69.2 Mb of A9 and 34.1-60.4 Mb of A10 had overlapping CIs with the results of our study (cqST.A9-2, 68.1-75 Mb and cqST.A10-1, 38.9-52 Mb) and one QTL related to SV located at 54.8-60.7 Mb of C6 had an overlapping CI with the results of our study (cqSV.C6-2, 50.9-55.9 Mb). This indicated that these QTLs were stably expressed in different genetic backgrounds and that they may be useful in molecular marker-assisted selection breeding.
Siliques represent a fruit type specific to members of the Brassicaceae family that form from the gynoecium after flowering (Seymour et al., 2008). The transformation of the gynoecium into a silique depends on whether the ovule has been successfully fertilized, and this signal may be produced by the pollen grains (Seymour et al., 2013). Successfully fertilized ovules result in developmental alterations of pistils from senescence to growing fruit. Thus, the failure of any part of these successive processes, including gynoecium formation, fertilization and silique growth, that involving cell proliferation, differentiation and expansion, will eventually affect the development of siliques and seeds, which are the determining factors of rapeseed yield. For example, plants with a BnaC9.SMG7b deletion exhibit reduced SPS and SL values, which are caused by developmental defects in the formation of functional FGs owing to the interruption of meiotic anaphase II . BnaC9.SMG7b, the successfully cloned rapeseed silique gene, also underlies the CI of cqSL.A10-3 in this study. LNG1 regulates longitudinal cell elongation in Arabidopsis, and lng1 dominant mutant plants are characterized by elongated siliques owing to longitudinally elongated and transversely narrowed cells in the siliques (Lee et al., 2006). CNGC18 is an essential Ca 2+ channel for pollen tube guidance in Arabidopsis. CNGC18 point mutations result in shorter siliques, reduced male fertility and fewer seeds per silique (Gao Q. et al., 2016). The Arabidopsis gene STY2 promotes the formation of the apical tissues of the gynoecium. The sty2 mutant lines exhibit many developmental defects in reproductive tissues, including shortened siliques and aborted ovules (Kuusk et al., 2002). KINβγ, a component of the regulatory subunit of the SNF1related protein kinase, is required for pollen germination on the stigma surface. KINβγ mutants display defects in organogenesis and growth, including shorter statures and siliques . ABCC13 is expressed in the seed coat and embryo, and its loss of function in Arabidopsis leads to decreased silique lengths and seed yields (El Guizani et al., 2014). ARID1 is required for sperm cell formation in Arabidopsis, and the ARID1 mutant shows reduced seed set and short siliques, which are caused by defects in gametophyte formation owing to an arrested mitotic cell cycle (Zheng et al., 2014). LAC4, a laccase gene regulated by miR397b, controls both lignin biosynthesis and seed yield in Arabidopsis, and overexpressing miR397b may increase silique lengths and seed sizes by modulating LAC4 . HTH is involved in regulating floral organ fusion in Arabidopsis, and the hth mutant shows reduced seed set and short siliques, which are caused by reduced pollen fertility and aborted ovules (Krolikowski et al., 2003). The homozygous aog1 in Arabidopsis shows reduced seed set and short siliques compared with the wild type, and this caused by a significant reduction in fertility owing to reduced pollen formation and severe defects in embryo sacs . PYL8 is a regulatory component of the ABA receptor, which is vital for regulating seed germination, root and shoot development and abiotic stress responses (Garcia-Maquilon et al., 2021), and a pyl8 mutant shows an ABA-insensitive phenotype, including a lower seed yield owing to fewer siliques having shorter lengths (Gonzalez-Guzman et al., 2012). Arabidopsis HMG1 encodes a 3-hydroxy-3-methylglutaryl coenzyme A reductase, which is a key element of the sterol biosynthetic process that is required for cell viability and growth. The T-DNA insertion hmg1 mutant in Arabidopsis shows shorter siliques and fewer seeds per silique compared with the wild type, and this was caused by reduced cell elongation and fertility owing to a reduced sterol level (Suzuki et al., 2004).

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.