Abstract

Kenaf (Hibiscus cannabinus L.) is a natural fibre crop that can be used for a variety of purposes and has various applications in industry. Despite this, its potential has not been fully exploited because of low yields and a narrow genetic base, limiting hybrids’ development. Based on this background, eight kenaf mutants and one commercial cultivar were selected and crossed in a half-diallel for general and specific combining abilities (GCA and SCA) to get the desired results done in this investigation. The 36 hybrid offspring and their parental lines were tested in the field over two environments. Diallel results based on Griffing B method 2 indicated significant differences for all characters studied except for GCA in top diameter and plant height and top diameter SCA, indicating the existence of both additive and nonadditive gene actions for the inheritance of the traits. The amplitude of GCA variation was much higher than that of SCA variation for all parameters except top diameter and node number, showing the additive gene’s prevalence and the likelihood of genetic advancement through selection. In both conditions, Hayman and Jinks graphical studies demonstrated that partial dominance controlled various fibre yield component parameters such as plant height, middle diameter, stick weight, and fibre weight. On the other hand, fibre yield and the majority of physical features indicated either dominance or overdominance gene action. Plant height, base diameter, core diameter, middle diameter, fresh stem weight, and stick weight all strongly positively correlated with fibre yield. These traits also had a higher proportion of additive effects, a moderate narrow-sense heritability, and a higher baker ratio, indicating successful indirect selection for fibre yield. The parents P1, P3, and P4 had the most dominant alleles for most of the features, while the parents P2, P7, and P9 had the most recessive alleles. The hybrids , , , , , , , , and outperformed the parents in terms of heterotic responses and showed that they have a lot of genetic potential for kenaf enhancement in tropical climates.

1. Introduction

In recent years, kenaf (Hibiscus cannabinus L.) has garnered significant interest as an intriguing multipurpose crop for producing energy, pulp, thermal insulation boards, and fibre-reinforced thermoplastic composites in Mediterranean nations. It has been employed as an alternate raw material for wood, pulp, paper, and textile manufacturing industries. The economic product of kenaf is the stem material, which consists of the outer bark (bast). The bark of the fibre-producing plant’s stem is removed as kenaf through a washing and drying procedure [1]. It creates a high-quality pulp that is also suited for making ropes, cordages, sacs, canvases, carpets, etc. In addition, the central section of the plant, known as the core, is abundant in cellulose and hemicellulose. This makes it a potential source for the synthesis of bioethanol; additionally, it can be used as an absorbent material for animal bedding [2]. According to Dempsey, the productivity of a kenaf F1 generation is 14 to 43 percent more than the parents [3]. As part of the ASEAN Free Trade Area (AFTA), the Malaysian government encouraged kenaf planting to replace tobacco to reduce tobacco import duties by 2010 [4]. Despite its wide distribution in tropical countries, kenaf productivity in Malaysia is low due to a lack of high-yielding varieties with fibre production of around 5 to 10 tonnes per hectare, and research into developing new kenaf varieties is still lacking [5]. Developing new kenaf varieties that produce high biomass is vital for effective kenaf production [6]. Thus, breeding kenaf in Malaysia is essential for producing high fibre content and adaptable growth in the local climates [5]. In some countries, including China, Russia, and Thailand, kenaf hybrids have been employed for commercial farming contributing to increasing production [7]. Due to its improved fibre quality and resistance to force, hybrid kenaf has gained a lot of interest [8]. In Malaysia, hybrid kenaf technology is considered novel, and more research is required to assess genetics, agronomy, and crop management. A hybrid kenaf breeding programme is required to develop a high-yielding, stable performance hybrid kenaf for the Malaysian environment. The availability of high-quality jute and kenaf seed would help to ensure fibre quality to a large extent [9]. As a result, developing locally adapted hybrid kenaf seeds is a viable alternative for enhancing national kenaf yield and increasing kenaf producer income.

The diallel analysis is a good way to screen parents for hybrid creation when compared to other mating designs. Diallel cross is an appropriate genetic design for obtaining genetic information on desired features in a short period. The ability to forecast the efficiency of selection in a population for improving yield in commercially relevant crops requires understanding the nature of gene action for various traits [10]. To better understand the nature of gene action involved in quantitative trait expression, the diallel mating design has been widely used in both self-pollinated and cross-pollinated species [11]. Plant breeders typically want to know how valuable a group of parents is for breeding. Combining ability analysis is a powerful tool for identifying superior hybrid parents with strong general combining ability (GCA) and progenies with stronger specific combining ability (SCA) [12]. In terms of gene activity in character quantitative inheritance, exploitation, and breeding, it is also useful for quantifying the genetic value of parents and crossings [13].

GCA refers to a parent’s average performance over several crosses [14]. On the other hand, SCA refers to a hybrid combination that outperforms or underperforms the parent inbred lines based on their average performance [15]. Cultivar parents with a strong GCA impact have additive gene action but do not always have a good SCA [16]. Meanwhile, genetic studies can benefit from SCA estimate to determine the type of gene action that impacts the phenotypes of interest. Nonadditive gene action is identified by a high SCA [17]. SCA effects were greater than GCA effects when measured in terms of average effects (components), demonstrating the relevance of nonadditive gene action in regulating kenaf yield component expression [18]. SCA and GCA statistics help breeders choose hybrids and parents to produce good offspring [19]. Parents with a high level of GCA generate strong hybrids [20]. The crossing lines with strong GCA are all crucial approaches for improving kenaf production.

The frequency of dominant and recessive genes and the number or cluster of genes affecting the quality of interest can all be used to improve breeding programme planning. Plant breeders use GCA and SCA data to choose appropriate genotypes and generate new high-yielding cultivars [21, 22]. Heritability is a fundamental genetic trait calculated in broad and narrow sense [23]. Selection can increase a trait with a high heritability estimate [2]. It will be extremely beneficial to improve the kenaf crop as a variety in Malaysia’s tropical environment by selecting genotypes with unique genetic histories [24]. The main purpose of this study was to find genotypes (parents and offspring) with good combining ability, which could help future kenaf improvement with high fibre production. The goal of this research was to use a diallel analysis to quantify the genetic variation of fibre yield components and examine if different kenaf mutants in tropical climates have different patterns of morpho-physiological parameters combining abilities.

2. Materials and Methods

2.1. Experimental Location

The research was conducted at Universiti Putra Malaysia, Selangor, Malaysia, at a latitude of 2°59 N, a longitude of 101°42 E, and 48 metres above sea level [2]. The field trials took place in two environments: one from June to September 2020 and the other from March to June 2021, both in a humid tropical climate. The maximum and minimum mean monthly temperatures (°C), total rainfall (mm), and relative humidity for the experimental periods are shown in Table 1.

2.2. Breeding Materials

Based on morpho molecular characterization [25, 26], eight kenaf mutants and one commercial cultivar were chosen as parents in half diallel crosses (Table 2). From February to May 2020, the parents were mated in all feasible combinations, except reciprocals, at University Putra Malaysia’s Field 10 to produce 36 F1 hybrids (Supplementary figure 1).

2.3. Development of F1 Using Half Diallel Mating Design

Although the genotype clustering differed greatly across the two methodologies, the genotypes were placed in similar groups. Consequently, when choosing parents for hybridized activity, integrating phenotypic and molecular data will provide a more accurate summary. Pure and healthy seeds were transplanted in three batches at 10-day intervals in the trial field to achieve flowering synchronisation. After the fruits had ripened, normal emasculation and pollination methods were used to obtain F1 seeds.

2.4. Experimental Design and Layout

Mechanical ploughing and laddering were used to prepare the soil for the plant’s cultivation. Kenaf seeds from 45 entries comprising nine parents and 36 F1s were sown in peat moss soil in germination trays. Twenty-one seedlings were transplanted into a plot with a distance of  cm inter- and intrarow spacing after two weeks, to a depth of two to 2.5 cm. There were 135 plots, and the entire plot measured . The experiment was carried out in a Randomised Complete Block Design (RCBD) with three replications using a table of random numbers [27].

2.5. Crop Husbandry

Standard kenaf production provided a good crop with the best management. NPK Green (15 : 15 : 15) and NPK Blue (12 : 12 : 17) were applied at 450 kg per hectare shortly after seeding and 40 days after planting. Intercultural activities like weeding, thinning, additional watering, and plant protection measure were carried out as needed throughout the cropping environment. The recommended cultural guidelines for growing a healthy crop were followed as described by Wong et al. [28].

2.6. Data Collection

A total of 10 fibre yield components were discovered. For each feature, observations were obtained from 10 randomly selected plants for each genotype for each replication. Plant height, base diameter, core diameter, middle diameter, top diameter, node number, days to 50 percent flowering, fresh stem weight, stick weight, and fibre weight were among the quantitative data collected in the same direction as the procedure of Al-Mamun et al. [25]. Days to 50 percent flowering were confirmed in the field, and remaining attributes were measured 90 days after transplantation in the laboratory.

2.7. Statistical Analysis
2.7.1. Performance of Hybrids

An analysis of variance (ANOVA) was used with the PROC GLM function of the Statistical Analysis System (SAS) computer software version 9.4 to determine the significance of variation among genotypes and blocks. A combined ANOVA from diallel (F1 and parents) was performed for each growth environment, eliminating reciprocals as indicated by Steel et al. [29]. The least significant difference (LSD) method was used to compare the mean performance among hybrids and paternal inbred lines.

2.7.2. Combining Ability Analysis

Griffing’s and Hayman’s techniques [30] were used to analyze the genetics of the diallel population. According to Baker [31], the combining ability ratio was determined. A modified Hayman [32] ANOVA was calculated after Morley Jones [33] adjustment. The genetic system influencing important kenaf features and general and specific combining abilities has been studied.

(1) Diallel Analysis by Griffing’s Method. Griffing’s method 2 (one set of F1 progenies and parental lines) was used to analyze the diallel [34]. The GCA of parents and SCA of hybrids were determined following Griffing’s method 2 model 1 (fixed effects). According to Zhang and Kang [35], the linear model for data from the combined analysis of performance and environments was as follows: where is the observed value of each experimental unit, is the mean of the population, is the environmental effects of th environment, is the replication effect in each environment, is the F1 hybrid effect, is the interaction effect between F1 hybrids and environments, and is the residual effects. where il denotes the interaction effect between GCA and environments for the th parent, denotes the interaction effect between GCA and environments for the th parent, and denotes the interaction effect between SCA and environments for the th F1 hybrid.

The combining ability ratio was calculated according to Baker [31] as below: where and are the mean squares for GCA and SCA, respectively.

(2) Diallel Analysis by Hayman’s Approach. Most plant breeding programs that attempt to improve yield and other related parameters start with parallel analysis. Breeders utilize it the most to determine the utility of varieties as parents and to analyze gene action in diverse qualities [34, 36]. The following are the main features of Hayman’s approach: (i) Hayman’s ANOVA; (ii) Vr-Wr analysis with graphical representation; and (iii) components of variation and genetic parameters.

(3) Hayman’s ANOVA and Morley Jones Modification. Hayman [32] analyzed variance for the entire diallel table, expanding on Yates’ work [37] in one direction. Reciprocal differences are typically absent; therefore, just one reciprocal cross from each pair is elevated. Morley Jones [33] modified Hayman’s technique in response to this situation. The sum of squares correlating to the additive impact (a) in Hayman’s model indicates dominance (b1), and additional dominance effects that genes can explain with one allele present in only one line (b2) (the residual n-1 lines were supposed to bring the same substitute allele), and remaining dominance effects (b3) are essentially a straightforward application of least squares suitable constants.

(4) Vr-Wr Regression Analysis. The regression coefficient was calculated by using the following formula: where is the variance of each array, is the covariance between parents and their offspring’s, is the variance of , is the variance of , and is the covariance between and .

Significant difference of “” from zero and unity was tested as follows:

These results were compared to the “” values for degrees of freedom.

The Vr-Wr analysis shows the significant variation of the regression coefficient () from unity. The absence of nonallelic interaction is demonstrated if the regression coefficient does not significantly depart from unity. The “” test for the uniformity of Vr-Wr values was used to verify the validity of diallel theory assumptions.

(5) Components of Variation and Genetic Parameters. Hayman [38] and Jinks [39] were used to determine the genetic and environmental components of variation and allied or related genetic factors in F1. The various components estimated were as follows: where is the covariance of additive and dominant effects in a single array. where and represent the percentage of positive and negative genes in the parents, respectively. as algebraic sum over the loci in heterozygous phase in the crosses.

is the variance of parents, is the mean variance of the arrays, is the mean covariance between parents and arrays, is the variance of the means of arrays, is the difference between the parents’ mean, and the progenies’ mean is known as the dominance relationship.

The estimation of genetic ratios was as follows: . The ratio of dominat to recessive allele were categorized as 0 = no dominance, >0<1 = partial dominance, 1 = complete dominance, and >1 = over dominance is indicated.

In the analysis of combining ability, the variance components in the ANOVA table were used to estimate broad-sense and narrow-sense heritability. The following equation was used to determine the heritability estimates for each kenaf attribute in each environment suggested by Rojas and Sprague [40]:

Broad-sense heritability:

Narrow-sense heritability:

(5)1. Vr-Wr Graph. The Vr-Wr graph can be used to interpret a diallel cross [32, 39]. By calculating the array variance (), parent-offspring covariance (), and regression of on , it is possible to test the adequacy of the simple additive dominance genetic model, discern the relative proportion of dominant and recessive alleles present in the common parents of each array, and determine the average level of dominance.

The (covariance) values for each array were calculated by using the following formula: where is the variance of th array and is the variance of parents.

By plotting values against values, the external limits of the parabola were found. The limiting parabola is drawn using the values. The expected values of are needed to create a regression line. The following formula is used to acquire these values for each array separately: where is the array means of covariances, is the array means of variances, is the individual array variance, and is the regression coefficient.

(5)2. Inference from the Vr-Wr Graph. The average degree of dominance is shown by the regression line’s position on the Vr-Wr graph. When the regression line crosses across the origin, it indicates complete dominance (). Partial dominance () exists when the regression line crosses the axis above the origin. It represents the absence of dominance when it crosses the origin, cuts the axis, and contacts the limiting parabola. Overdominance is shown when the regression line crosses the axis below the origin.

The order of dominance for each parent is determined by the position of parental points on the regression line. Parents with more dominant genes are closer to their ancestors than parents with more recessive genes. Parents with the equal amount of dominant and recessive genes fall in the middle of the range.

3. Results and Discussion

3.1. Variation in Pooled Environments across All Genotypes

All traits tested had highly significant () variations in environment and genotypes (parents and offspring) except for the middle and top diameters (Table 3). In the pooled quantitative data of the two environments for nine parents and their crosses, genotypes by environment interaction () were significant ( or 0.05) except for plant height. The CV percent for fibre yield components varies from 12.02 to 69.05, indicating a wide range of variability in the features studied.

The genotypes performed differently across the environments studied, based on the results of the analysis on the combined data. These results indicate that the effects of genes controlling these traits were expressed differently in different environments. Most of the traits showed significant differences across genotypes and replications in the combined analyses of variance, indicating that the materials had enough genetic variation to improve these traits. The crop’s extensive variation in botanical and agromorphological properties suggests that the genotypes are genetically diverse [41]. For all studied traits, the most significant environments, mean squares, and their interactions with genotypes were detected, indicating that the environment had enough seasonal variability to cause fluctuations in all population component rankings, i.e., different genotypes ranked differently from environment to environment. Alza and Fernandez-Martinez [42] and Abd El-Satar et al. [30] reported similar findings.

3.2. Genotype Performance Averaged over Two Environments

Table 4 shows the average comparison between all genotypes. The heights of the plants ranged from 230.44 to 280.59 cm, with parent P4 having the highest mean value and P2 having the lowest. The base diameters ranged from 20.88 mm (P6) to 25.74 mm (P1), with P4 having the highest core diameter (21.88 mm) and P6 having the lowest (17.15 mm). The middle diameter ranged from 9.69 mm (P2) to 12.49 mm (P3), with P2 and P9 having the largest and lowest top diameters, respectively. The most nodes were found in P5 (7.38), whereas the fewest were found in P2 (4.15). The period it takes to achieve 50 percent flowering varied between 55.67 and 69.17 days. The hybrid P7 was the first to reach maturity, while the parent P2 was the last (Table 4). In addition to other fibre yield components, Golam et al. [43] discovered that 50 percent flowering and days to maturity may be the two most relevant traits in identifying kenaf accessions.

The highest and lowest fresh stem weight was recorded for P4 (351.42 g) and P6 (222.41 g), respectively. The stick weights varied from 71.35 g (P6) to 123.24 g (P3). Among the parental lines, P1 had the highest mean fibre weight per plant (26.71 g), followed by P4 (24.78 g) and P8 (23.59 g), and P2 had the lowest (17.55 g) (Table 4). According to the findings, environmental factors influenced the expression of kenaf genetic characteristics. Parent P3 had the highest mean stick weight per plant (123.24 g), followed by P7 (96.26 g) and P8 (95.06 g). However, as the node number increased, the fibre production per plant decreased.

3.3. Combining Ability Analysis of Kenaf Hybrids for Fibre Yield Components
3.3.1. Combining Ability Effects by Griffing’s Method

For all analyzed traits in the combined data (Table 5) for both environments, analysis of variance for combining ability was performed using Griffing’s technique [34]. GCA is well understood to be a function of additive gene effects and additive sections of epistatic variance. In contrast, SCA is a function of nonadditive gene effects and the remaining epistatic variance [44]. GCA effects were highly significant for all parameters measured except top diameter, according to the agronomic performance of the parental inbred lines. The SCA effects of the parental inbred lines were highly significant () for all traits except plant height and top diameter and significant () for days to 50 percent flowering, showing that additive gene effects predominantly influenced variability between those traits. In contrast, both additive and nonadditive gene effects controlled plant height and top diameter variations. The interaction of GCA with environment effects was highly significant () for all traits except plant height and fresh stem weight and significant () for base diameter and days to 50 percent flowering, indicating that the inbred lines’ GCA was influenced by the environmental conditions surrounding the hybrids.

On the other hand, the SCA by environment interaction effects were highly significant for all traits except plant height and top diameter, while significant for days to 50 percent flowering, indicating that environments altered the effects of specific hybrid combinations on the expression of these traits. As a result, the impacts of nonadditive genes on phenotypes interacted with the environment more. For all traits except top diameter and node number, the mean squares of GCA were bigger than the mean squares of SCA, indicating that additive gene action predominated in our study across both environments (Table 5). Furthermore, the Baker ratio for the examined traits was altered from 0.52 to 0.93, indicating that the additive effect plays a larger role in determining the traits.

The combining ability analysis revealed significant GCA for all traits except top diameter and significant SCA for all traits except plant height and top diameter from the combined data of the two environments. The results indicate that variations for those traits were controlled mainly by additive gene effects, while both additive and nonadditive gene effects controlled variations for plant height and top diameter. For quantitative traits, the predominance of additive over nonadditive gene effects was rather prevalent. Similar findings were found in kenaf by Jianmin et al. [45] and Heliyanto et al. [46]. The analyzed traits (except top diameter and node number) had large mean square GCA values, implying that the parental materials examined had a high level of genetic variability. The investigated characteristics (except for plant height and fresh stem weight) significantly affected GCA by environment interaction, demonstrating that environmental variation influenced additive gene action. Likewise, the characters’ environments had a significant impact on how they evolved [47]. Mostofa et al. [48] and Cai et al. [49] reported comparable outcomes in Hibiscus cannabinus, Sobhan [50] in Hibiscus sabdariffa, and Khatun [51] in Corchorus capsularis by combining ability studies.

For all traits except top diameter and node number, the mean squares of GCA were larger than SCA’s, showing that additive gene action effects played a significant role in their inheritance. Although additive gene effects contributed the most to the variability in most traits, dominance and overdominance significantly influenced the genetic system that controls yield components. Because of the strong additive gene effects, this conclusion usually favours the breeding selection technique [52]. Mostofa et al. [53] discovered that fibre weight correlated with the presence of additive gene influences in the development of the characteristics. One dominant gene pair was responsible for days to first flowering and plant height, whereas three genes were responsible for raw fibre yield [54]. Pace et al. [18] discovered additive gene action to be more important for yield components such as plant height, fresh and dry bark weight, and usable stick in kenaf. Baker ratios [31] greater than 0.80 indicated that additive gene effects played a larger role in the genetic control of plant height, middle diameter, days to 50 percent flowering, fresh stem weight, stick weight, and fibre weight (Table 5). Additive influences contributed more to genetic variation, according to the Baker ratio. Additive variances are linked to heredity and react well to selection for improving desired traits. Anwar et al. [55] and Hassan et al. [56] both found similar findings in wheat research.

3.3.2. General and Specific Combining Ability of Kenaf Genotypes for Fibre Yield Components

Plant height in a taller stature combination should have a positive GCA effect, while the node number should have a negative GCA effect [57]. In the pooled data (Table 6), the parent P7 had the maximum plant height (12.89) and the lowest negative GCA values for node number (-0.42), suggesting that they were strong general combiners for quality fibre yield and might be exploited in future breeding attempts. Parent P4 showed the most favourably significant GCA effect () of plant height (9.11), base diameter (2.16), core diameter (1.93), middle diameter (0.70), fresh stem weight (54.54), stick weight (14.14), and fibre weight (3.02). With the exception of the P4, parent P1 (1.49) and P3 (1.38) were also positive and significant general combiners () for fibre weight, with parent P3 having a positive and very significant GCA effect () for stick weight (12.66). The parent P2 had the most positively significant GCA effect () for days to 50 percent flowering (2.99). Parent P7 had the smallest GCA effect for days to 50 percent flowering (-3.89). The parent P2 showed the lowest negative and highly significant GCA effect () for plant height (-16.83). For base diameter (-1.08), core diameter (-1.01), and fresh stem weight (-21.63), the parent P9 exhibited the lowest negative GCA values.

The effects of SCA on hybrids in various contexts in pooled environments are shown in Table 6. SCA impacts were seen in all 36 hybrids investigated, with 20 showing positive (desired direction) SCA influences on plant height. The hybrids with significant and beneficial SCA effects () were produced by the cross (23.63), which was evaluated as a good specific combiner for tallness, followed by (19.91). The SCA effects on the base diameter ranged from -4.24 to 5.61. Of the 23 positive SCA effects, two ( and ) were found to be highly significant () for base diameter, while six (, , , , , and ) were reported to be significant (). Six of the 36 cross combinations had significantly positive SCA values for kenaf core diameter, showing heterotic performance compared to their parents’ mean. The cross had the highest positive SCA effect (5.97), followed by , , and , indicating that they were found to be highly significant () for the trait. The SCA effects for middle diameter stem varied from -1.52 to 2.97. The hybrid had the best SCA effects, followed by and , considered the best specific combiners for this trait. Another 18 crosses had positive but insignificant SCA values, implying that they could be used as average specific combiners. The SCA effects on top diameter ranged from -0.71 to 1.06. Three crossings exhibited significant positive SCA effects, , , and were identified as good specific combiners for top diameter.

SCA impacts ranged from -1.94 to 2.52 for a given node number. Negative SCA effects were found in 14 crosses, with one () being highly significant (-1.94) and another () being significant (-1.31), showing that these hybrids had good SCA for lower branch stems and improved fibre yield. SCA impacts ranged from -5.55 to 4.91 days to 50 percent flowering. The significance of four of the 22 positive SCA effects was determined to be quite high. The cross of the 15 positive SCA effects, which displayed a good specific combining capacity for this characteristic, had the most significant SCA effects. Five more crossings of the 23 positive SCA effects produced positive and significant SCA effects, making them good specific combiners for fresh stem weight. Three crosses, (103.70), (91.15), and (78.19) had highly significant positive SCA effects, whereas two crosses, (66.50) and (58.02), had significant positive SCA effects, indicating that these hybrids had good SCA for fresh stem weight.

In kenaf, 23 cross pairings had positive SCA effects, while the other 13 had negative for stick weight per plant, SCA impacts ranging from -19.01 to 31.08. had the largest significant positive SCA impact (31.08), followed by (30.16), (29.98), and (29.54); it was an excellent specific combiner for the stick weight per plant attribute, meaning that it was a good specific combiner. Three crosses, (11.70), (5.31), and (4.43), showed highly significant positive SCA effects for fibre weight, while four crosses, (7.93), (4.06), (3.53), and (3.45), showed significant positive SCA effects. This suggests that hybrids with higher fibre weights than their parents’ means are the best specific combiners for increased kenaf fibre weight.

The research demonstrated that parents’ GCA effects were linked to their crossings’ SCA impacts, which had the highest significant positive intensity. Additive genetic variance is a major contributor to the GCA component. Because additive variance may be fixed, selecting qualities regulated by additive variance is a very effective strategy [58]. As a result, each parent’s GCA variation has a major impact on the parents’ decisions. A parent with a significant positive GCA effect is a good general combiner [59]. The high value of GCA for the traits of interest was dispersed across genotypes in this study, indicating that none of the genotypes employed had the best combination of GCA values for the various traits of interest. Parent P4 was among the best parental lines for plant height, base diameter, core diameter, middle diameter, fresh stem weight, stick weight, and fibre weight content, showing the accumulation of favourable additive genes for these traits in the hybrids. The parent P3 had favourable impacts on base diameter, fresh stem weight, stick weight, and fibre weight; on the other hand, the parent P1 was the greatest general combiner for fibre weight and contributed positively to the hybrid for these traits. The main purpose of this breeding effort is to create high-yielding hybrids with high potential fibre production comparable to or equal to existing cultivars.

Crosses with high vs. low general combiners outperform others in terms of yield components in general. According to an investigation of combining ability impacts [60], high-specific combiners involved high vs. high, high vs. average, high vs. low, average vs. average, average vs. low, and low vs. low combining parents. According to Jinks, overdominance and epistasis induced severe SCA effects in crosses with high vs. low and low general combiners [61]. In crosses comprising high vs. low general combiners for yield components, mutual cancellation of heterosis components, especially dominance and its association, resulted in unfavourable SCA effects [62]. Crossing two parents with low general combiners produces high performance, which is attributable to complementary gene action [63].

SCA effects were found to be significant for most yield traits, as mentioned in Table 6. For fibre weight, the hybrids , , , and were shown to be the most effective specific combiners. The cross is good for plant height, base diameter, middle diameter, top diameter, stick weight, and fibre weight. As the finest specific combiners for stick weight, the hybrids , , and were selected. Days to 50 percent flowering for hybrids were also picked as unique combiners.

Complementing gene effects could explain the strong SCA effects of these crosses. Crossings of , , , , , , , , and showed promising heterotic responses and could be beneficial in future breeding program. Hybrid vigour can be induced by dominant, overdominant, or epistatic gene action in any combination of parents, according to Moll and Stuber [64]. In this study, both additive and nonadditive genetic components influenced fibre yield components, with nonadditive gene action dominating most of the characters.

3.3.3. Hayman’s Method of Analysis of Variance and Genetic Component Estimate

(1) Hayman’s Analysis of Variance. After Morley Jones’ modification [33], Hayman’s ANOVA results for all tested traits in mixed environments were identical for additive effect (a), dominance effect (b), and error variance (Table 7). In combined environments, additive genetic effects (a) were significant for most of the traits, implying that both additive and dominance gene actions are involved in the inheritance of these traits. Most traits had a much larger a than b, indicating that additive effects were more important. For the top diameter and node number, a highly significant mean square was detected due to the interaction of bl with the environment, and a significant mean square was detected for fibre weight and fresh stem weight, revealing that changes in soil types and climatic factors at each environment influenced the mean departure of the F1s from their mid-parental values for the two traits [30]. The remaining traits, on the other hand, exhibited negligible environment interaction mean squares with bi, indicating that they were stable in both environments. In addition, all traits had insignificant mean squares of b2 with environment interaction, indicating that b2 and variety heterosis components were consistent across the two environments [30]. All traits had insignificant mean squares of b3-environment interaction, demonstrating that b3 and specific heterosis components were stable across the two environments [30].

In combined environments, additive genetic effects (a) were significantly larger than dominance genetic effects (b), suggesting the presence of both additive and dominant gene actions in the inheritance of these traits [65]. The findings agree with those of Akter [66]. Xu [67] discovered a larger magnitude of the additive genetic variance than dominance in the inheritance of Hibiscus cannabinus flowering days. Significant mean square due to the interaction of bl with the environment was detected for top diameter, node number, fibre weight, and fresh stem weight, indicating that the mean deviation of the F1s from their midparental values for the two traits was most likely influenced by differences in soil types and climate conditions at each environment [30].

(2) Components of Variation and Genetic Parameters. The genetic component analysis of Hayman [33] differed little from the analyses for additive (D), dominance (H1), and environmental error in most parameters (Table 8). The data were submitted to Hayman’s diallel analysis [37] to isolate the genetic variance components and their ratios for all traits studied. For the eight examined traits out of 10, including fibre weight, positive dominance (H1 and H2) values were found, demonstrating the relevance of nonadditive components in the inheritance of these traits [68]. For the 10 traits, the magnitude of dominance (H1 and H2) was greater than the additive component (D), indicating overdominance [30]. H1 was greater than H2 for six of 10 traits, showing that the frequency of gene distribution in the parents was asymmetrical, implying dominance for all of these traits [30]. These findings matched with Akter [66], who observed that positive and negative alleles were not equal in the parents for all features in tossa jute. The value was negative for all traits suggesting the presence of dominant recessive genes in the parents influencing these traits [30].

Because the balance of positive and negative alleles (UV parameter) is less than 0.25 for 6 out of the 10 traits and the b2 component is significant for a few, the frequency of dominant and recessive alleles cannot be assumed equal. Furthermore, the importance of nonadditive gene action in influencing the traits was demonstrated because h2B and h2N for the examined traits varied from 0 to 37.86 percent and 0 to 23.05 percent, respectively. Fibre weight shows a moderate heritability of 37.86 percent. Plant height had the highest narrow-sense heritability of 23.05 percent in the combined environment, followed by fresh stem weight (20.06 percent). On the other hand, poor narrow-sense heritability indicated that nonadditive gene effects predominated in trait inheritance [30]. The average degree of dominance at overall loci, as estimated by the (H1/D)0.5 ratio, was found to be greater than one for all traits except base diameter, core diameter, node number, and fibre weight, indicated by the Vr-Wr graph; overdominance gene effects have a role in the inheritance of most examined traits. On the other hand, the remaining attributes had a zero (H1/D)0.5 ratio, indicating no dominance in the control of the traits. Patil and Thombre [54] found a higher proportion of additive genetic components in kenaf in terms of days to flowering, plant height, and fibre length. The dominant impact is important for all traits except plant height, middle diameter, top diameter, and days to 50 percent flowering, according to the results of the ratio dominant to recessive alleles (h2/H1) (Table 8).

According to Hayman’s analysis of variance and components of variation in gene action, all traits were influenced by both additive and dominant (nonadditive) gene actions. It also revealed that positive and negative alleles were not equally common in the parents and asymmetrical gene distribution at the loci, indicating that all these traits were dominant. However, the value of additive (D) and dominance (H1) of sqr (H1/D) were zero, which indicated the presence of no dominance effects for the base diameter, core diameter, node number, and fibre weight, while there were dominance effects for the other traits showing 1.34 to 2.06. Heritability in the narrow-sense assessed from genetic components appeared to be moderate to low for most traits. The moderate narrow-sense heritability values of plant height and fresh stem weight were indicative of their early generation selection success for breeding improvement [7]. All parameters of kenaf (fresh plant yield, defoliated stalk yield, plant height, basal diameter, middle diameter, and dry stalk yield) had a low narrow-sense heritability, according to Liu [7]. According to Falconer and Mackay [68], low additive effects and high dominant gene action caused reduction in narrow-sense heritability. Pedigree selection based on plant height and fresh stem weight can improve fibre yield.

3.4. Vr-Wr Regression Analysis

For the traits and their related statistics, Table 9 and Figure 1 show a graphical assessment of parent-offspring covariance (Wr) and array variance (Vr). The regression analysis showed that days to 50 percent flowering fitted with a simple additive dominance genetic model (Ho: and Ho: ) involving the nine parents studied. This is considerably different from zero but not from unity, showing no nonallelic interaction for these traits’ inheritance. However, plant height and fresh stem weight showed highly significant “” values in F1. It therefore did not follow the models clearly, which were significantly different both from zero and unity, indicating the presence of an interallelic interaction in the inheritance of these traits. The remaining of the nine traits in combined environments were considered to follow epistatic or nonallelic gene interaction as the regression coefficient significantly differed from unity but not from zero involving all nine parents studied.

The regression analysis revealed that days to 50 percent flowering differ considerably from zero but not from unity, indicating that nonallelic interactions do not inherit these traits. The findings were comparable to those of Sobhan [50], who found that in Hibiscus sabdariffa, a close relative of H. cannabinus, the days to flower suited the model well and the regression coefficient was substantially different from zero but not from unity. Plant height and fresh stem weight differed significantly both from zero and unity, indicating the presence of an inter-allelic interaction in the inheritance of these traits. The remaining seven traits in combined environments were considered to follow epistatic or nonallelic gene interaction as the regression coefficient significantly differed from unity but not from zero involving all nine parents studied. The regression coefficient was not statistically significant when compared to zero. The deviation from unity for all 10 characters except for the top diameter indicated the validity of the hypothesis proposed by Hayman [62]. The nonsignificant values satisfied the uniformity of covariance and variance (Wr-Vr) and thus supported the validation of the assumptions of Hayman [32] for these traits.

3.4.1. Average Degree of Dominance

The deviation from the origin of the point where the regression line cuts the Wr axis, one of the information points acquired from the graph, measures the average level of dominance. In view of this, the intercept of the regression line on the covariance axis (Figures 1(b), 1(c), 1(f), 1(g), and 1(h)) is below the origin of base diameter, core diameter, node number, days to 50 percent flowering, and fresh stem weight, implying that these traits show overdominance. However, for plant height, middle diameter, stick weight, and fibre weight (Figures 1(a), 1(d), 1(i), and 1(j)), partial dominance was indicated by the regression line cutting the Wr axis above the place of origin. These findings match those of Vr-Wr regression analysis by Akter [66] for Corchorus olitorius and Khatun [51] for Corchorus capsularis.

3.4.2. Parental Inheritance of Dominant and Recessive Genes

The array points along the regression line reflect the distribution of dominant and recessive genes between the parents. The parents with the most dominant genes are closest to the origin, whereas the parents with the most recessive genes are the furthest away. It is noticed that in Figures 1(a)1(d) and 1(j), the parent P4 falls near the point of origin, in plant height, base diameter, core diameter, middle diameter, and fibre weight, as well as in Figures 1(d), 1(e), and 1(i), the parent P3 also falls near the point of origin, in middle diameter, top diameter, and stick weight. In Figures 1(b), 1(c), 1(g), and 1(f), the parent P1 also falls near the point of origin, in base diameter, core diameter, days to 50 percent flowering, and fibre weight. Results demonstrated the presence of dominant genes for the traits as mentioned earlier in the parents P1, P3, and P4. Those parents that are located closest to the regression line’s start had more dominant genes for that trait, whereas those located farthest away had more recessive genes [69]. In contrast, the parents P2, P8, and P9 are the furthest from the origin in most features, implying that they have the most recessive genes for these traits.

3.5. Correlation among Fibre Yield Components

As demonstrated in Table 10, simple Pearson phenotypic correlation coefficients for fibre yield components were computed using proc corr SAS software version 9.4. Most traits do not exist in isolation; rather, they establish a complex connection with one another that ultimately determines the yield. The -value of a correlation coefficient conveys the idea of a relationship between two unique traits by making a connection. Traits’ correlation coefficients ranged from 0.05 to 0.99, indicating that the phenotypic variation is greater. Negative linear relationship, no linear relationship, and perfect positive linear relationship are shown by -values of 1, 0, and +1, respectively. According to Field [70], values ranging from 0.5 to 1, 0.3 to 0.5, and 0 to 0.3 indicate strong, moderate, and low negative linear relationships, respectively, while values ranging from 0.5 to 1, 0.3 to 0.5, and 0 to 0.3 indicate strong, moderate, and low positive linear relationships, respectively.

Plant height, base diameter, core diameter, middle diameter, fresh stem weight, and stick/fibre weight all demonstrated a positive and highly significant correlation with dry fibre and stick weight. This was the most significant correlation yet observed for a single trait. For at least one of the types, in combined environments, a negative phenotypic correlation was seen in days to 50 percent flowering with fibre weight and stick weight. The fact that these attributes had a negative correlation with fibre weight suggested that they would not improve the trait. Plant height, base diameter, core diameter, middle diameter, node number, fresh stem weight, and stick weight all showed significant positive correlations with fibre weight, while top diameter and days to 50 percent flowering showed significant negative correlations. In contrast, there were significant correlations with all traits in the case of combined environment. Between the types, there was no discernible pattern or similarity. This suggested that progenies within each type may have a variety of strategies for expressing a certain phenotype. Significant and strong correlations will aid in selecting traits and qualities that influence commercial yield traits, making future selection easier and more reliable.

The negative relationship between fibre weight, stick weight, plant height, base diameter, core diameter, node number, and fresh stem weight suggests that early maturity reduces fibre yield (bast and core fibre) per plant [71]. Furthermore, the length of the growth environment may impact stem yield and other yield components like plant height. Early flowering plants have a shorter life cycle and produce less fibre than late flowering plants [69]. Several other significant correlations were also discovered between the various qualities investigated in this research. In the future, these correlations could be used to speed up the selection of superior progenies without quantifying all traits. The correlation analysis allows each trait’s performance to be modelled with other major correlated traits, saving time and money in future investigations.

4. Conclusions

The fibre yield components evaluated in this study were genetically controlled by additive and nonadditive variants. From the combined data of the two environments, the combining ability analysis revealed strong GCA for all traits except top diameter and considerable SCA for all traits except plant height and top diameter. Additive gene effects mostly governed variations in those parameters. In contrast, according to the findings, plant height and top diameter were influenced by both additive and nonadditive gene effects. Except for top diameter and node number, GCA effects were higher than SCA effects, as shown by mean squares, indicating that additive gene action predominates for these traits. In conclusion, the parental lines P1 (ML5), P3 (ML36-10), and P4 (ML36-24) were outstanding general combiners for fibre yield components. In contrast, the hybrids P1 (ML5) × P4 (ML36-24), P1 (ML5) × P9 (ML36-21(2)), P2 (ML9) × P3 (ML36-10), P2 (ML9) × P5 (ML36-25), P4 (ML36-24) × P6 (ML36-27), P4 (ML36-24) × P7 (BJRI Kenaf4), P4 (ML36-24) × P9 (ML36-21(2)), P5 (ML36-25) × P8 (MLRing4P2), and P7 (BJRI Kenaf4) × P9 (ML36-21(2)) showed promising heterotic responses and could be beneficial in future breeding program. Plant height, middle diameter, days to first flowering, stick weight, and fibre weight indicated a clear partial dominance in the variance and covariance graphs. For most of the phenotypes evaluated, parents P2 (ML9), P7 (BJRI Kenaf4), and P9 (ML36-21(2)) had the most recessive genes, while parents P1 (ML5), P3 (ML36-10), and P4 (ML36-24) had the most dominant alleles. The Baker ratio showed that selection-based approaches highly responsive to additive effects enhance genetic improvement in the fibre output of hybridization-based breeding programmes. Because the fibre yield components studied have a high to moderate broad-sense heritability, a high Baker ratio, and moderate to low narrow-sense heritability, it was determined that selection should be done in advanced generations after homozygosity and genetically fixed.

Data Availability

All data are provided in the results section of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This study was funded by the Bangladesh Agricultural Research Council (BARC-Project of NATP Phase-II), People’s Republic of Bangladesh; World Bank; IFAD; and Universiti Putra Malaysia (research grant: vote number 6282507).

Supplementary Materials

Supplementary figure 1: the nine genotypes were mated in half diallel at Field 10 in University Putra Malaysia, excluding the reciprocals, to produce 36 F1 hybrids. Comparison of photographs showing stem pigmentation and leaf shape of parents and F1 population. (Supplementary Materials)