Next Article in Journal
The Content of Anthocyanins in Cowpea (Vigna unguiculata (L.) Walp.) Seeds and Contribution of the MYB Gene Cluster to Their Coloration Pattern
Next Article in Special Issue
Novel Pathogen–Plant Host Interaction: Colletotrichum jiangxiense and Fraxinus americana L. (White Ash) in a Sentinel Garden in China
Previous Article in Journal
Genome-Wide Identification and Characterization of the NAC Gene Family and Its Involvement in Cold Response in Dendrobium officinale
Previous Article in Special Issue
Phenotypic Evaluation of Soybean Genotypes for Their Reaction to a Mississippi Isolate of Phakopsora pachyrhizi Causing Soybean Rust
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance- and Resistance-Related Early Responses of Colombian Elite Rubber Tree Genotypes under Low Pressure of South American Leaf Blight: Implications for Disease Management in the Amazon

by
Lyda Constanza Galindo-Rodríguez
1,2,3,*,
Armando Sterling
2,*,
Herminton Muñoz-Ramirez
3 and
Jesica Andrea Fonseca-Restrepo
2,3
1
Doctoral Program in Natural Sciences and Sustainable Development, Faculty of Agricultural Sciences, Universidad de la Amazonia, Florencia 180001, Colombia
2
Laboratory of Phytopathology, Amazonian Scientific Research Institute Sinchi, Faculty of Basic Sciences, Universidad de la Amazonia, Florencia 180001, Colombia
3
Mycology and Phytoprotection Laboratory, Faculty of Basic Sciences, Universidad de la Amazonia, Florencia 180001, Colombia
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(20), 3627; https://doi.org/10.3390/plants12203627
Submission received: 14 September 2023 / Revised: 18 October 2023 / Accepted: 18 October 2023 / Published: 20 October 2023
(This article belongs to the Special Issue Leaf Diseases and Management)

Abstract

:
The cultivation of Hevea brasiliensis, the primary commercial source of natural rubber, is strongly impacted by South American leaf blight (SALB) disease, caused by the fungus Pseudocercospora ulei. Various management strategies have been implemented, including the selection of resistant genotypes and the identification of escape zones. This study evaluated the growth, early yield, and resistance to SALB of nine Colombian elite genotypes from the ECC-100 series and IAN 873 clone (control) in a large-scale clone trial in an area with low SALB pressure in the Colombian Amazon during 2017–2020. Favorable early performance was evident, although there was a significant increase in the severity and sporulation of P. ulei over time, especially in the ECC 35, ECC 60, and IAN 873 genotypes. However, these scores indicate low susceptibility. Genotypes with higher resistance to SALB demonstrated greater growth and early yield compared to more highly susceptible genotypes. The ECC 64, ECC 73, ECC 90, ECC 25, and ECC 29 genotypes were more desirable in low SALB pressure zones due to their higher resistance and early performance. It is important to highlight that this research contributes to the selection of new SALB-resistant Colombian genotypes of H. brasiliensis. However, it is also necessary to evaluate the productivity of these selections in the mature stage and long-term resistance to SALB before recommending and promoting their commercial adoption in the Colombian Amazon.

Graphical Abstract

1. Introduction

Hevea brasilienses (Willd. ex A.Juss.) Müll.Arg. is the primary commercial source of natural rubber (NR) and supports the livelihoods of approximately 20 million people worldwide [1]. This resource is essential as a raw material in the manufacturing of over 50,000 products across various industrial sectors, including healthcare, construction, and automotive. In the latter sector, its usage is particularly prominent in the production of high-strength tires and rubber products [2,3].
Between 1961 and 2020, there was a remarkable increase in the global average yield of rubber cultivation, rising from 0.5 to 1.2 tons per hectare, and annual production surged from 2.1 to 14.8 million tons [4]. This increase as achieved through the implementation of advanced agronomic practices, the careful selection of highly productive and disease-resistant genotypes, the identification of suitable cultivation areas, and the expansion of agricultural land dedicated to rubber tree cultivation [5,6]. The scarcity of other commercial sources of NR offering comparable quality and production capacity has justified research into various aspects of H. brasiliensis cultivation [7,8]. In contrast, rubber tree cultivation faces various phytosanitary limitations that, if not addressed in a timely manner, can lead to significant economic losses. In South America, the primary disease affecting H. brasiliensis is South American leaf blight (SALB), caused by the fungus Pseudocercospora ulei (Henn.) [9].
P. ulei is an exclusive hemibiotrophic fungus of the Hevea genus [10]. This fungus has three developmental stages in its life cycle. The first stage corresponds to the conidial phase (i.e., asexual phase), which is responsible for the infectious stage of the disease. The lesions caused during the fungus’s development primarily appear on the leaves, stems, petioles, and young fruits. On young leaves, there are spots on the underside that spread across the entire leaf surface, varying in coffee-gray color, shape, and size, with a diameter ranging from 8 to 10 mm [11,12,13,14]. The second stage, corresponding to the onset of the sexual phase (i.e., spermogonial phase) in leaves at the beginning of physiological maturity, is not relevant for the pathogen’s dissemination [15,16]. Finally, the third stage (i.e., sexual phase) occurs in leaves at physiological maturity, with structures of resistance known as stromata [17,18]. These stromata are brown, globular, usually abaxial, carbonaceous, and superficially clustered, and sometimes form rings around the perforations of the leaves, with a size of between 200 and 400 mm [18].
This disease results in repeated defoliations and reduced latex production, and, in severe cases, it can lead to the death of the plant [19], thereby limiting production in the Amazon region [20].
Chemical management of SALB is possible; however, due to the height of the trees (20–30 m) and the required treatment frequency, it is economically unfeasible and raises significant ecological and public health concerns [21]. For the management of SALB, the primary approaches implemented have been the selection of fungus-resistant clones [22] and the identification of escape zones (areas where rubber grows with lower fungal pressure) [1]. Escape zones are areas where the pressure of the P. ulei fungus is lower, allowing for the normal development and production of H. brasiliensis [23]. These areas have a well-defined dry season that can last from one to five months, a period during which the fungus’s life cycle is interrupted, as the pathogen cannot find favorable conditions for germination, infect receptive tissue, promote the development of the infection process, and subsequently cause damage. This enables the trees to escape the disease, especially during the natural defoliation–refoliation period [24,25]. In addition to the dry season, other climatic variables are considered, such as annual average temperature, evapotranspiration, water deficit, relative humidity, altitude, and soil suitability in the area [9,23,25,26]. This strategy has been implemented in Brazil [27,28], Ecuador [1,25], and Colombia [9,23,26], where the production of natural rubber has been severely impacted by SALB incidence [29,30,31].
According to Food and Agriculture Organization (FAO) data [4], the yield of rubber tree plantations in Colombia has reached 1.6 tons per hectare, surpassing the global average of 1.2 tons per hectare but falling below the averages of plantations in Mexico, Guatemala, and India, where average yields have reached 4.0, 2.8, and 2.1 tons per hectare, respectively. In the Amazon, despite the introduction of genotypes with good resistance and yield characteristics, the cultivation continues to be heavily impacted by the disease, which diminishes productivity [32,33]. Therefore, it is necessary to select new high-yield clones that are more resistant to SALB.
With the establishment of the rubber tree clone trials network in Latin America [1], clones/genotypes that exhibit partial or complete resistance to SALB have been identified [24]. Previous research conducted in the Colombian Amazon has identified superior genotypes in terms of performance, quality, and adaptability to climatic variations [31,34]. For the Elite Caquetá Colombia (ECC) series, Sterling et al. [30] observed favorable early growth in some genotypes in non-SALB escape zones. However, the genotypes with high resistance showed lower growth rates compared to the more susceptible genotypes. Therefore, further research is necessary for the final selection and commercial adoption of resistant genotypes with desirable productivity traits under high and low SALB pressure.
In Colombia, there has been an effort to expand the genetic base of the rubber tree by assessing the growth, physiology, latex production, and response to SALB of introduced clones and genotypes from the ECC-100 series in small-scale and large-scale clone trials under non-SALB escape conditions [30]. However, in escape zones or areas with low disease pressure, it is expected that the environmental response of these genotypes in agronomic and phytosanitary terms will be more favorable. Therefore, the objective of this study was to evaluate the growth, early yield, and SALB resistance of nine Colombian elite genotypes from the ECC-100 series and IAN 873 (control) under low disease pressure in the Colombian Amazon region.

2. Results

2.1. Climate

Between July 2017 and June 2020, the study area recorded an average monthly precipitation of 209.5 mm and an annual average precipitation of 2514.6 mm. The average annual temperature was 25.1 °C. The maximum temperature ranged from 28.5 °C in June and July 2019 to 34.1 °C in January 2020, whereas the minimum temperature ranged from 21.84 °C in December 2019 to 23.6 °C in February 2020. The average thermal amplitude was 8.7 °C, with an average maximum temperature of 31.6 °C and an average minimum temperature of 22.7 °C. The lowest relative humidity was recorded in January 2020 with a value of 69.3%, whereas the highest occurred in June 2019 with a value of 88.8%, and the annual average was of 81.5%. The dew point ranged from 19.4 °C in January 2020 to 22.9 °C in May 2018 (Figure 1).

2.2. Effects of the Fixed Factors

The Akaike (AIC) criterion was minimized in the unstructured with heterogeneous variances (UNH) model for all measured variables except for the type of reaction (TR), for which the first-order autoregressive with heterogeneous variances (ARH1) model performed better in terms of the AIC criterion. For the Bayesian criterion (BIC), the best model was UNH for the circumference of the trunk in the pre-tapping phase (CTpt), the severity of the attack on young leaves (AT1), and TR variables. The ARH1 model was superior for the BIC criterion in terms of the severity of the attack on mature leaves (AT2), and compound symmetry with heterogeneous variances (CSH) model was the best for the BIC in terms of stroma density (ST) (Table 1). Based on the above, the unstructured model with heterogeneous variances (UNH) was the best model. Consequently, significant differences were observed between the genotypes, between the years, and in the genotype × year interaction (p < 0.05) (Table 2).

2.3. Growth and Early Yield

The annual tree growth data are shown in Table 3. Significant differences were found between the 10 genotypes (p < 0.05). Higher mean values of CTpt were observed in the ECC 64, ECC 25, and ECC 83 genotypes, where the ECC 64 genotype had the highest average CTpt over the three years of evaluation (p < 0.05). The ECC 25 genotype showed higher means than the other genotypes in years 2 and 3 (p < 0.05), where the ECC 64 and ECC 25 genotypes had the greatest vigor (Table 3).
The CTo ranged from 25.00 cm (ECC 35) to 27.92 cm (ECC 64), where the ECC 64 genotype was significantly higher than the other genotypes (p < 0.05). The ECC 64, ECC 90, and ECC 25 genotypes had higher values of CTo, with 27.92 cm, 26.94 cm, and 26.86 cm, respectively (Table 4).
In the early tapping phase, the ECC 25, ECC 64, and ECC 83 genotypes had higher means of dry rubber yield (DRY), with 5.57, 3.76, and 3.30 g.tree−1.tap−1, respectively. These values were higher than that of IAN 873 (p < 0.05). On the other hand, the ECC 29, ECC 35, and ECC 66 genotypes showed lower DRY means, with values of 0.74, 0.42, and 1.09 g.tree−1.tap−1, respectively, although these values did not differ significantly from that of IAN 873, which had a mean of 1.88 g.tree−1.tap−1 (Table 4).
The dry rubber content (DRC) was significantly higher in the IAN 873, ECC 25, ECC 73, ECC 29, ECC 64, and ECC 66 genotypes, with means of 32.83, 32.45, 32.42, 31.44, 31.38, and 20.75%, respectively. The ECC 35 and ECC 90 genotypes had lower DRC values (Table 4).

2.4. Resistance to SALB

The severity of a P. ulei attack on H. brasiliesis is evident in Figure 2. Figure 2a,b represent minor and moderate attacks, respectively, which are observed in genotypes with low susceptibility to SALB. Figure 2c,d represent severe and very severe attacks, respectively. Severe attacks occurred in clones with high susceptibility to SALB, as was the case with IAN 873.
According to Figure 3, there was a low SALB incidence in the first year, but in years two and three, there was a significant increase in susceptibility to P. ulei. In young leaves, mean scores of AT1 and TR were lower than 1. Some genotypes had mean scores > 1 and >3 for AT2 and ST, respectively. In terms of susceptibility, it was evident that the ECC 35 genotype had the highest SALB mean scores, whereas ECC 73 had the lowest scores.
For AT1, the ECC 64, ECC 29, and ECC 90 genotypes had lower attack intensity, with AT1 < 1. Despite the significant increase in the attacks during the third year, the scores were lower than 1 (Figure 3a). The TR means were statistically different between genotypes, with a significant increase in the second and third years. The ECC 64, ECC 90, and ECC 29 genotypes were the least susceptible, with TR = 0. On the other hand, the ECC 60, ECC 83, and ECC 35 genotypes showed higher sporulation intensities, especially in the third year. However, TR means < 1 were observed in all genotypes over the years (Figure 3b).
In mature leaves, a significant increase in the AT2 means were observed between the genotypes and between the years. The ECC 35 genotype was the most susceptible, with an average AT2 score > 1, followed by IAN 873, ECC 83, ECC 60, ECC 25, and ECC 66. However, the highest score was observed in the third year, where ECC 35 had the highest severity scores (AT2 > 2) and was significantly superior to other genotypes in all years. Furthermore, for the second year, ECC 35 exhibited the highest susceptibility to SALB. On the other hand, the ECC 64, ECC 73, ECC 29, and ECC 90 genotypes were the least susceptible, with an average AT2 < 1 (Figure 3c).
With respect to the sexual phase of P. ulei, it was found that the ECC 35 and IAN 873 genotypes had higher formation of stroma (ST > 2). In the third year, ECC 35 was the most susceptible, with mean ST scores > 3, followed by IAN 873, ECC 60, ECC 25, ECC 66, ECC 83, and ECC 60 (all, ST > 2). In the second year, ECC 25, ECC 83, and ECC 60 had mean ST scores > 2. Overall, ECC 64 was the least susceptible, with an average ST < 1. The ECC 73, ECC 64, and ECC 90 genotypes had mean ST scores <1 and were the least susceptible in the third year (Figure 3d).
The co-inertia analysis between the data matrices did not show significant covariation (p > 0.05). This indicates that the co-variation structure between early performance and SALB resistance was not entirely similar to the structures obtained in the individual analyses of each variable matrix. However, the first two axes of the co-inertia analysis explained 88.36% and 11.1% of the total variability, respectively (Figure 4a,b, respectively), indicating that the simultaneous analysis of both matrices was suitable for explaining the observed ordering pattern in the 10 genotypes of H. brasiliensis.
The most susceptible genotypes positioned themselves in the lower-left quadrant of Figure 4c, characterized by lower vigor and lower yield (DRY and DRC, respectively) (Figure 4d), as well as higher scores of AT2 and ST. Additionally, a positive association was observed between the ECC 60 genotype and the AT1 and TR variables (i.e., a higher susceptibility in young leaves) (Figure 4e). On the other hand, the ECC 64 genotype showed higher values of CTo and CTpty3 (Figure 4d). This genotype had greater vigor and lower SALB susceptibility, positioning itself in the lower-right quadrant of Figure 4c. The ECC 25, ECC 29, ECC 73, ECC 64, and ECC 90 genotypes were located towards the right side of the plot and had the best results in terms of growth and early yield, along with low susceptibility to SALB.

3. Discussion

This study evaluated the resistance to SALB, growth, and early yield of nine Colombian elite genotypes of H. brasiliensis from the ECC-100 series (Elite, Caquetá Colombia) and IAN 873 (control) under conditions of low SALB pressure.
The climatic characteristics of the study area (i.e., average temperature of 25.1 °C, average relative humidity in the driest months ranging from 69 to 75%, and precipitation of 2514.6 mm) are typical of regions with low SALB pressure [9,23,26,31,38]. These climatic parameters are unfavorable for the production of P. ulei spores [39] but they are optimal for rubber tree cultivation [9,23,24,25].
For the three years of the study, between November and March, there was a reduction in precipitation and relative humidity, along with an increase in temperature compared to the period from April to October. According to Sterling et al. [30], these climatic variations are optimal for the growth of H. brasiliensis in areas with no escape from SALB, resulting in more vigorous trees in the first few years after planting. On the other hand, Rivano et al. [24] found that climatic variations allowed for the foliage recovery of some clones during the dry season, and, under suitable phytosanitary conditions, they were not affected by SALB in a suboptimal rubber cultivation area in Ecuador. Furthermore, Sterling et al. [31] found that the climatic conditions in San Vicente del Caguán contribute to SALB resistance for introduced rubber tree genotypes from Brazil, Peru, and Guatemala. Our results support the idea that these climatic conditions are favorable, promoting better vigor and greater resistance to SALB for Colombian rubber tree genotypes.
In this study, the circumference of the trunk (CT) was higher than that reported in previous studies in the three first years of growth [30,31]. In escape SALB conditions, Sterling et al. [31] found CTpt values that ranged from 6.50 to 8.61 cm in the first year of planting. In contrast, in non-escape SALB conditions, Sterling et al. [30] reported CT values that ranged from 3.64 to 5.98 cm in the first year for ECC-100 series genotypes. Our results showed CTpt values of between 8.73 and 10.69 cm in the first year, higher than those reported in the Colombian Amazon region.
In this study, in the second year, the average CTpt values ranged from 13.33 to 20.44 cm, similar to those reported by Sterling et al. [40]. For the third year, the average CTpt values ranged from 20.23 cm to 27.48 cm, which were higher than those reported by Gireesh et al. [41] in India, where 20 promising genotypes had a girth of between 15.3 and 23.0 cm. The ECC 64 and ECC 25 genotypes were the most vigorous, with an average annual trunk circumference increase of 7 to 8 cm, similar to the results of Rivano et al. [24], who observed an 8 cm increase in circumference for the two most vigorous clones in a suboptimal rubber cultivation zone in Ecuador. This demonstrates the high potential of ECC-100 series genotypes for zones with low SALB pressure.
The rubber tree immaturity period ranges from four to eight years before conventional tapping is applied. Genotypes with high CT in the early years of planting reach the tapping period faster, allowing for early exploitation and early cost recovery for tree maintenance [41,42,43]. Therefore, one of the goals of rubber tree breeders is to identify genotypes that combine high yield and vigorous growth in early evaluation [44]. In this regard, the Colombian elite genotypes evaluated in this study showed a high CT in the first three years, similar to the highest-performing clones evaluated by Rivano et al. [24] and superior to those evaluated by Gireesh et al. [41].
Climate conditions and clonal variation have a direct impact on the phenotypic expression of H. brasiliensis [45,46]. This indicates that genotype performance depends not only on the genetic nature of the trees but also on the environmental variations of the planting site [47]. In this study, a greater growth of ECC-100 series genotypes was observed compared to previous studies conducted in areas with favorable conditions for SALB [30]. These findings highlight the importance of considering both the selection of suitable genotypes and environmental conditions to maximize the performance of H. brasiliensis.
The dry rubber content (DRC) is an important parameter for assessing latex quality, as it is one of the main components that directly or indirectly contribute to rubber yield, and it determines the amount of rubber in a latex sample [48,49]. In this study, genotypes had a DRC that ranged from 26.24% to 32.83%. These results are lower than those reported by Mendoza-Vargas and Jiménez-Forero [50] (i.e., DRC of 34.3%) in a promising clonal collection from the FX series in Cundinamarca, Colombia. Furthermore, the percentages obtained are lower than those reported by Ajith et al. [51] in two locations in India (i.e., DRC ranged from 35.23% to 48.9%). The DRC values of IAN 873 (control) were also lower than those reported by Quesada-Méndez et al. [48] in genotypes three years after planting (i.e., DRC of 40.91%) and lower than those reported by Mendoza-Vargas and Jiménez-Forero [50] (i.e., DRC of 41.3%).
The differences observed in DRC between the genotypes were mainly due to variations in the colloidal composition of latex, which depends on soil conditions, climate, tapping duration, diseases, and plant material. These factors primarily affect the DRC composition. Overall, the DRC values ranged between 27 and 45% [52,53]. This indicates that the results of our study fall within optimal composition ranges.
Higher DRC values represent greater productivity in terms of dry rubber yield. However, elevated percentages also imply an increase in the viscosity of the colloidal latex composition, which can lead to blockages in the laticifer mantle and obstruction of the latex flow from the tapping panel [48,54], which can reduce dry rubber yield (DRY). In this study, DRY values ranged from 0.42 to 5.57 g.tree−1.tap−1, which are lower than those reported by Ajith et al. [51] (i.e., DRY of between 2.88 and 54.5 g.tree−1.tap−1).
Assessing the temporal dynamics of SALB resistance is crucial to determining optimal areas for crop establishment [30]. This study revealed different susceptibility patterns depending on the year and genotype. The severity of the attacks on both young and mature leaves increased significantly over the years, especially in the third year and in the ECC 35, ECC 60, and IAN 873 genotypes. These results are consistent with those of Sterling et al. [30], who found that the ECC 35 and ECC 60 genotypes were affected the most in non-escape SALB areas. Furthermore, the ECC 64, ECC 73, ECC 90, ECC 25, and ECC 29 genotypes showed lower susceptibility, positively correlating with higher vigor and early yield. These results differ from those of Sterling et al. [30], who found that the genotypes with higher resistance did not necessarily exhibit greater vigor.
With respect to asexual or conidial sporulation, it occurs on young leaves, and during this phase, conidia can disperse for short distances during dry periods [55], but it does not ensure pathogen survival [56]. However, in very rainy seasons, conidia can be particularly abundant [55]. The San Vicente del Caguán site in this study is characterized by low precipitation and well-defined dry periods (i.e., semi-humid warm climate), which is unfavorable for the production of P. ulei conidia (i.e., low TR scores).
In relation to sexual sporulation, which provides quantitative information about the intensity of ascospores on mature leaves [36], the attack severity increased in the second and third years. In the second year, a moderate level of attack was evident for all genotypes, with the exception of the ECC 29 and ECC 64 genotypes, which had lower attack levels. In the third year, there was an increase in susceptibility; however, moderate attack levels were still recorded for all genotypes, with the exception of the ECC 35 genotype, which suffered severe attacks. Overall, we observed an increase in ascospore production over the years. However, the ECC 64, ECC 73, ECC 90, ECC 25, and ECC 29 genotypes were the least susceptible. The SALB scores in the first year were similar to those reported by Rivano et al. [36] in the CDC 312, FDR 4575, FDR 5597, and MDF 180 clones, which were classified as completely resistant, and in the third year, the SALB scores were similar to those reported by Sterling et al. [34].
Considering that clones with high yield and SALB resistance are rare [57], this dual-choice criterion is very important in rubber tree breeding programs. Based on these results, the ECC 64, ECC 73, ECC 90, ECC 25, and ECC 29 genotypes are considered potential candidates for establishment in areas with low SALB pressure in the Colombian Amazon owing to their low susceptibility to SALB and high performance in the first three years of growth.

4. Materials and Methods

4.1. Study Area

4.1.1. Location

The study was conducted from July 2017 to June 2020 on the Parcela No. 12 farm, located in the Buenos Aires rural settlement area in the Municipality of San Vicente del Caguán, Caquetá, Colombia. This site is located in the southern part of the Colombian Amazon (02°01′42.62″ N and 74°54′38.95″ W, at an elevation of 346 m above sea level. According to the classification by IGAC [58], the landscape in the study zona is hilly, characterized by hills with heights of less than 300 m and slopes ranging from 7 to 12 percent. The region is a zone with low SALB pressure, making it a promising region for rubber cultivation [26,33].

4.1.2. Soils

The soils in Caquetá have unfavorable characteristics such as high compaction, low fertility and infiltration, poor drainage, shallow depth, and a pH range of between 4.5 and 5.8. Furthermore, these soils have high levels of aluminum saturation and low contents of sodium (Na), phosphorus (P), potassium (K), and magnesium (Mg) [58].
The study area had a very acidic pH of 5.05, an organic matter (OM) content of 1.21%, and an organic carbon (OC) content of 0.715%. Saturation levels were 17.51% (Ca), 1.18% (K), 4.99% (Mg), and 2.23% (Na). Elemental contents were as follows: iron (Fe) (161.59 mg kg−1), manganese (Mn) (4.27 mg kg−1), copper (Cu) (1.71 mg kg−1), zinc (Zn) (0.54 mg kg−1), boron (B) (0.19 mg kg−1), phosphorus (P) (0.84 mg kg−1), sulfur (S) (6.91 mg kg−1), and total nitrogen (N) (0.06%).

4.1.3. Climate

The climatic conditions in the study area correspond to a semi-humid warm climate, typical of the tropical wet region. It records an average temperature of 25.4 °C, a relative humidity of 82.3%, and an annual precipitation of 2503 mm [58]. The climate follows a monomodal regime with a transition to a dry period from July to October, followed by a dry period with less rain between November and February and, finally, a rainy period from March to June [59].
Microclimatic data for the study area were collected with a Portable Microclimate Station (Decagon Devices Inc., Pullman, WA, USA), from which monthly averages were calculated between July 2017 and June 2020 (Figure 1).

4.2. Experimental Design and Plot Management

A large-scale clonal trial (LSCT) [60] was established using a randomized complete block design (RCBD) with four replications (Fisher blocks). Each treatment consisted of planting 60 trees in single rows, with a planting distance of 7.0 × 3.0 m, equivalent to a planting density of 476 trees per hectare, and the total area of the LSCT was 5.04 hectares. The measured plot unit was 1260 square meters, containing 60 trees organized in 3 rows of 20 trees [30,31,34].
With respect to the management of the experimental plots, the fertilization plan employed by Sterling et al. [61] was applied every six months. Thus, a compound fertilizer (N (15%), P2O5 (15%), K2O (15%), CaO (2.2%), and SSO4 (1.7%)) was used at a rate of 150 g per plant. Minor elements included N (8%), P2O5 (5%), CaO (18%), MgO (6%), S (1.6%), B (1%), Cu (0.14%), Mo (0.005%), and Zn (2.5%), applied at a rate of 75 g per plant. Additionally, organic compost was applied at a dosage of 1000 g per plant. Weed control was carried out every three months, and no phytosanitary controls were performed [30,31,34].

4.3. Plant Material

Nine Colombian elite H. brasiliensis genotypes from the ECC-100 series (Elite Caquetá Colombia) were evaluated: ECC 25, ECC 29, ECC 35, ECC 60, ECC 64, ECC 66, ECC 73, ECC 83, and ECC 90. These genotypes were obtained from asexual propagation (cloning) of plus trees of H. brasiliensis from rubber-producing farms in Caquetá, Colombia, through natural cross-pollination and were characterized using molecular markers [30].
Evaluations carried out in small-scale [62] and large-scale clonal trials [30] in non-SALB escape zones demonstrated favorable performance in terms of growth and partial resistance to SALB compared to traditional cultivars introduced in the Colombian Amazon.
In this study, the cultivar IAN 873 was used as a control, which is a clone widely cultivated in Colombia due to its good agronomic performance [63], which contrasts with its high susceptibility to P. ulei, as evidenced in the last 15 years mainly in non-SALB escape zones [33].

4.4. Evaluations and Data Analysis

4.4.1. Growth and Early Yield

The annual trunk circumference at 1.2 m above the ground was recorded during the pre-tapping phase (CTpt) of all trees from July 2017 to June 2020 and in the early tapping phase in the third year (CTpty3). Then, the CT mean was calculated only for the trees that reached the threshold for extraction (CT ≥ 25 cm) at 1.2 m above the ground at the opening of the tapping panel (circumference of the trunk at opening, CTo) of the trees suitable for early tapping in the third year after planting (modified from Silva et al. [64]).
The dry rubber yield (DRY) (g.tree−1.tap−1) was recorded during the first six months of tapping, following the methods described by Meenakumari et al. [37] and Silva et al. [64]. The dry rubber content (DRC) as a percentage was measured according to the method described by Lacote et al. [65] and Mendoza-Vargas and Jiménez-Forero [50].

4.4.2. Resistance to SALB

The resistance of H. brasiliensis to SALB was assessed monthly in 50% of the trees in each plot (30 trees) between July 2017 and June 2020 using a method adapted from Rivano et al. [36]. The data were analyzed by year. To evaluate resistance to P. ulei, both young and mature leaves (stages C and D, respectively) were examined [31].
The severity of the attacks on young leaves (AT1) and mature leaves (AT2) was assessed using a scale adapted by Rivano et al. [36] based on the scale proposed by Chee and Holliday [18]. This scale defines five categories based on the percentage of leaf surface with symptoms or signs of the disease: 0 represents no attack (<1%), 1 represents a minor attack (1–5%), 2 represents a moderate attack (6–15%), 3 represents a severe attack (16–30%), and 4 represents a very severe attack (>30%).
In young leaves, the intensity of conidial sporulation (asexual phase of the fungus) on the lesions formed on the leaf surface at stage C was visually assessed using the variable “type of reaction” (TR). This variable is based on a scale from 1 to 6 proposed by Junqueira et al. [66] and adapted by Mattos et al. [67]. The scores on this scale are as follows: (1) necrotic lesions without spores; (2) non-necrotic lesions without spores; (3) very scarce and heterogeneous sporulation on the underside of the leaf; (4) high, heterogeneous, or partial sporulation on the underside of the leaf; (5) very high and uniform sporulation covering the entire lesion on the underside of the leaf; and (6) very high sporulation covering the entire lesion on both the underside and upper side of the leaf.
The intensity of sexual sporulation of the fungus on the upper side of the leaves at stage D was assessed using the variable “stroma density” (ST). A scale proposed by Rivano et al. [36] was used, defining four scores based on the number of stromata per leaflet: (0) absence of stromata, (1) fewer than 5 stromata, (2) between 5 and 10 stromata, (3) between 11 and 30 stromata, and (4) more than 30 stromata. For all of these variables, the most severely affected leaves on each plant were evaluated and the average of each variable was calculated for each replicate of each genotype [34].

4.4.3. Statistical Analysis

A linear mixed-effects model (LME) was used to analyze the fixed factors genotype, year, and their interaction for the variables CTpt, AT1, AT2, TR, and ST. QQ-plots of the model residuals were used to assess the normality assumption and fitted plot (residual vs. predicted) for the homogeneity of variances. A transformation was applied to the AT1, AT2, TR, and ST variables to normalize the data, using the formula (value + 0.5)1/2. The effect of genotype on the growth variables was evaluated at the end of the evaluation cycle. Block and plot effects associated with genotypes within blocks were considered random effects. Residual variance was modeled to adjust the heteroscedasticity by genotype, and residual correlation for successive observations was considered using models for longitudinal data: compound symmetry with heterogeneous variances (CSH), first-order autoregressive with heterogeneous variances (ARH1), continuous first-order autoregressive with heterogeneous variances (CARH1), and unstructured with heterogeneous variances (UNH) [35]. The Akaike (AIC), Bayesian (BIC), and log-likelihood (Log Lik) criteria were used to select the structure of variances and residual correlations [34,35].
To analyze the effect of genotype on variables related to early yield (CTo, DRY, and DRC), a linear mixed-effects model (LME) with independent errors was used. The analyses were performed with the “lme” function (linear mixed-effects models) in the “nlme” package (linear and nonlinear mixed-effects models) [68] in the R language software, version 3.5.1 [69], and the interface in InfoStat v. 2018 [70]. Differences between means of the variables were analyzed using Fisher’s LSD test (α = 0.05).
To explore the relationships between performance (growth and early yield) and SALB resistance in the 10 genotypes, two matrices were constructed with each group of variables, and a Principal Component Analysis (PCA) was performed as an ordination method. The results of the ordination of the two groups of variables were carried out with a co-inertia analysis with the Monte Carlo test [71]. The software used for this analysis was the ADE-4 package [72] in R 3.5.1 [69].

5. Conclusions

The results of this study support the idea that, under favorable climatic conditions (i.e., zones with low SALB pressure), the genotypes from the ECC-100 series had higher yields and lower SALB susceptibility, and, therefore, could be successfully established in regions such as the Colombian Amazon. This research contributes to the development of more effective management strategies based on genetic resistance to the disease.
Most tested genotypes showed low susceptibility to P. ulei. Over time, the impact of the disease increased significantly; however, when compared to other studies, the susceptibility level was very low (partial resistance to SALB, with mean scores of TR < 2 and ST < 3). Furthermore, the genotypes with lower susceptibility to SALB expressed higher early performance indices (CTpt > 23 cm, DRY > 3 g.tree−1.tap− 1, and DRC > 31% in the third year after planting). The optimal early responses in yield, growth, and SALB resistance reported in this study, although preliminary, are of interest for the selection of new desirable H. brasiliensis genotypes for the Amazon region. However, it is necessary to analyze the productive performance and SALB resistance for a longer period comprising the crop maturity stage before generating a commercial recommendation under the Colombian Amazon region conditions.

Author Contributions

Conceptualization, L.C.G.-R. and A.S.; methodology, L.C.G.-R., A.S. and J.A.F.-R.; software, A.S. and H.M.-R.; validation, L.C.G.-R. and A.S.; formal analysis, L.C.G.-R. and A.S.; investigation, L.C.G.-R., A.S., H.M.-R. and J.A.F.-R.; resources, L.C.G.-R. and A.S.; data curation, L.C.G.-R. and A.S.; writing—original draft preparation, L.C.G.-R., A.S. and H.M.-R.; writing—review and editing, L.C.G.-R., A.S. and H.M.-R.; visualization, L.C.G.-R., A.S. and H.M.-R.; supervision, L.C.G.-R. and A.S.; project administration, L.C.G.-R. and A.S.; funding acquisition, L.C.G.-R. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research is part of the “Ampliación de la base genética de caucho natural, Caquetá, Amazonía” project, funded by the Fondo de Ciencia, Tecnología e Innovación (FCTeI)—SGR, Contract 59/2013; the Instituto Amazónico de Investigaciones Científicas SINCHI—Gobernación del Caquetá; the Universidad de la Amazonia; and the Asociación de Reforestadores y Cultivadores de Caucho del Caquetá Asoheca.

Data Availability Statement

Data are available from the authors upon request.

Acknowledgments

The authors thank ASOHECA, the owners of the experimental farms, and the project’s technical personnel, who made the field work possible.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rivano, F.; Mattos, C.R.R.; Cardoso, S.E.A.; Martinez, M.; Cevallos, V.; Le Guen, V.; Garcia, D. Breeding Hevea brasiliensis for Yield, Growth and SALB resistance for high disease environments. Ind. Crops Prod. 2013, 44, 659–670. [Google Scholar] [CrossRef]
  2. Mazlan, S.; Jaafar, N.; Wahab, A.; Sulaiman, Z.; Rajandas, H.; Zulperi, D. Major diseases of rubber (Hevea brasiliensis) in Malaysia. Pertanika J. Sch. Res. Rev. 2019, 5, 10–21. [Google Scholar]
  3. Prabhakaran, K.P. The Agronomy and Economy of Important Tree Crops of the Developing World. Nair, K.P.P., Ed.; Elsevier: New Delhi, India, 2010; ISBN 978-0-12-384677-8. [Google Scholar]
  4. Food and Agriculture Organization—FAO FAOSTAT. 2023. Available online: http://www.fao.org/faostat/ (accessed on 13 September 2023).
  5. de Castro Sant’Anna, I.; Gouvêa, L.R.L.; Martins, M.A.; Scaloppi, E.J., Jr.; de Freitas, R.S.; Gonçalves, P. Genetic diversity associated with natural rubber quality in elite genotypes of the rubber tree. Sci. Rep. 2021, 11, 1081. [Google Scholar] [CrossRef]
  6. Van Truong, V.; Van Lam, L.; Thuy, L.M.; Dinh, H.D.; Thao, N.T.; Rivallan, R.; Van Biet, H.; Guen, V. Le population genetic structure of a thousand rubber tree accessions from wild Rondônia populations conserved in Vietnam. Genet. Resour. Crop Evol. 2020, 67, 475–487. [Google Scholar]
  7. Maiden, N.A.; Syd Ali, N.; Ahmad, K.; Atan, S.; Wong, M.Y. Growth and physiological responses of Hevea brasiliensis to Rigidoporus microporus infection. J. Rubber Res. 2022, 25, 213–221. [Google Scholar] [CrossRef]
  8. Malmonge, J.A.; Camillo, E.C.; Moreno, R.M.B.; Mattoso, L.H.C.; McMahan, C.M. Comparative study on the technological properties of latex and natural rubber from Hancornia speciosa Gomes and Hevea brasiliensis. J. Appl. Polym. Sci. 2008, 111, 2986–2991. [Google Scholar] [CrossRef]
  9. Jaimes, Y.; Rojas, J.; Cilas, C.; Furtado, E.L. Suitable climate for rubber trees affected by the South American Leaf Blight (SALB): Example for identification of escape zones in the Colombian Middle Magdalena. Crop Prot. 2016, 81, 99–114. [Google Scholar] [CrossRef]
  10. Tapiero, A.L. Enfermedades y su control. In Modelo Productivo Para el Cultivo de Caucho Natural (Hevea brasiliensis) en el Bajo Cauca Antioqueño y en el sur de Córdoba; Gutiérrez Vanegas, A.J., Correa Pinilla, D.E., Gómez Valderrama, J.A., Amado Saavedra, G.M., Martínez Atencia, J., Álvarez-Láinez, M.L., Córdoba Gaona, O.d.J., Tapiero Ortiz, A.L., Monsalve García, D.A., Gil Restrepo, J.P., Eds.; AGROSAVIA: Mosquera, Colombia, 2022; pp. 145–162. [Google Scholar]
  11. de Araújo, A.E.; Filho, A.Ν.; Nóbrega, M.Β.M.; Sousa, N.R.; Santos, J.W. dos Reação de Dez Clones de Seringueira (Hevea benthamiana) a três isolados de Microcyclus ulei. Acta Amaz. 2001, 31, 349. [Google Scholar] [CrossRef]
  12. Junqueira, N.T.V.; Gasparotto, L.; Lieberei, R.; Normando, M.C.S.; Lima, M.I.P.M. Especialização fisiológica de Microcyclus ulei em diferentes espécies de Seringueira: Identificação de grupos de patótipos. Fitopatol. Bras. 1989, 14, 205. [Google Scholar]
  13. Gasparotto, L.; Lieberei, R.; Trindade, D. In vitro conidia germination of Microcyclus ulei and its sensitivity to fungicides. Fitopatol. Bras. 1984, 9, 505–511. [Google Scholar]
  14. Gasparotto, L.; Pereira, J.C.R.; Furtado, E.L. Capítulo 3: Doenças das folhas. In Doenças da Seringueira no Brasil; Embrapa: Brasilia, Brazil, 2012; pp. 39–176. [Google Scholar]
  15. Sterling, A.; Rodríguez, O.L.; Correa, J. Evaluación fitosanitaria y resistencia en campo a Microcyclus ulei de Hevea Brasiliensis en la Amazonia Colombiana. Rev. Colomb. Amaz. 2009, 5, 80–90. [Google Scholar]
  16. Garzón, C. Principales Enfermedades y Plagas en el Cultivo de Caucho (Hevea brasiliensis) con Énfasis en la Amazonia Colombiana; Instituto Amazónico de Investigaciones Científicas SINCHI: Bogotá, Colombia, 2000. [Google Scholar]
  17. Jaimes, Y.Y.; Rojas, J.; Furtado, E.L. Clones de Hevea brasiliensis de alta productividad caracterizados por resistencia a Microcyclus ulei en jardin clonal en el Magdalena Medio Colombiano. Summa Phytopathol. 2015, 41, 115–120. [Google Scholar]
  18. Chee, K.H.; Holliday, P. South American Leaf Blight of Hevea rubber. Malaysian Rubber Res. Dev. Board Monogr. 1986, 13, 50. [Google Scholar]
  19. Lieberei, R. South American Leaf Blight of the Rubber Tree (Hevea spp.): New steps in plant domestication using physiological features and molecular markers. Ann. Bot. 2007, 100, 1125–1142. [Google Scholar] [CrossRef]
  20. da Hora Júnior, B.T.; de Macedo, D.M.; Barreto, R.W.; Evans, H.C.; Mattos, C.R.R.; Maffia, L.A.; Mizubuti, E.S.G. Erasing the Past: A New Identity for the samoclean pathogen causing South American Leaf Blight of rubber. PLoS ONE 2014, 9, e104750. [Google Scholar] [CrossRef]
  21. Le Guen, V.; Garcia, D.; Mattos, C.R.R.; Doaré, F.; Lespinasse, D.; Seguin, M. Bypassing of a polygenic Microcyclus ulei resistance in rubber tree, analyzed by QTL detection. New Phytol. 2007, 173, 335–345. [Google Scholar] [CrossRef] [PubMed]
  22. Sobha, S.; Rekha, K.; Uthup, T.K. Abstract Biotechnological Advances in rubber tree (Hevea brasiliensis Muell. Arg.) breeding. In Advances in Plant Breeding Strategies: Industrial and Food Crops; Al-Khayri, J.M., Jain, S.M., Johnson, D.V., Eds.; Springer Nature: Cham, Switzerland, 2019; Volume 6, pp. 179–236. [Google Scholar]
  23. Correa-Pinilla, D.E.; Gutiérrez-Vanegas, A.J.; Gil-Restrepo, J.P.; Martínez-Atencia, J.; Córdoba-Gaona, O.d.J. Agroecological and South American Leaf Blight Escape Zones for rubber cultivation in Colombia. Agron. J. 2022, 114, 2830–2844. [Google Scholar] [CrossRef]
  24. Rivano, F.; Vera, J.; Cevallos, V.; Almeida, D.; Maldonado, L.; Flori, A. Performance of 10 Hevea brasiliensis Clones in Ecuador, under South American Leaf Blight escape conditions. Ind. Crops Prod. 2016, 94, 762–773. [Google Scholar] [CrossRef]
  25. Rivano, F.; Maldonado, L.; Simbaña, B.; Lucero, R.; Gohet, E.; Cevallos, V.; Yugcha, T. Suitable Rubber Growing in Ecuador: An approach to South American Leaf Blight. Ind. Crops Prod. 2015, 66, 262–270. [Google Scholar] [CrossRef]
  26. Castañeda, A. Zonas Aptas Para el Cultivo de Caucho en Colombia; Serie Técnica No. 39; CONIF: Bogotá, Colombia, 1997. [Google Scholar]
  27. Pilau, F.G.; Marin, F.R.; de Camargo, M.B.P.; Assad, E.D.; Pinto, H.S.; Barbarisi, B.F. Zoneamento agroclimático da heveicultura para as regiões Sudeste e Centro-Oeste Do Brasil. Rev. Bras. Agrometeorol. Piracicaba 2010, 15, 161–168. [Google Scholar]
  28. Carmo, C.A.F.d.S.d.; Lumbreras, J.F.; Naime, U.J.; Gonçalves, A.O.; Fidalgo, E.C.C.; Ágilo, M.L.D.; Lima, J.A.d.S. Aspectos culrurais e zoneamento da Seringueirano Estado Do Rio de Janeiro. Bol. Pesqui. Desenvolv. 2004, 60, 49. [Google Scholar]
  29. Sterling, A.; Pimentel-Parra, G.A.; Virguez-Díaz, Y.R.; Suárez-Córdoba, Y.D.; Hoyos-Duarte, J.D.; Fonseca-Restrepo, J.A. Long-Term resistance in promising rubber tree genotypes as a breeding source for improving South American Leaf Blight management under high disease incidence in the Colombian Amazon. Crop Prot. 2021, 150, 105817. [Google Scholar] [CrossRef]
  30. Sterling, A.; Galindo-Rodríguez, L.C.; Suárez-Córdoba, Y.D.; Velasco-Anacona, G.; Andrade-Ramírez, T.; Gómez-Torres, A.K. Early assessing performance and tesistance of Colombian Rubber Tree Genotypes under High South American Leaf Blight Pressure in Amazon. Ind. Crops Prod. 2019, 141, 111775. [Google Scholar] [CrossRef]
  31. Sterling, A.; Martínez-Viuche, E.J.; Suárez-Córdoba, Y.D.; Agudelo-Sánchez, A.A.; Fonseca-Restrepo, J.A.; Andrade-Ramírez, T.K.; Virguez-Díaz, Y.R. Assessing Growth, early yielding and resistance in rubber tree clones under low South American Leaf Blight pressure in the Amazon Region, Colombia. Ind. Crops Prod. 2020, 158, 112958. [Google Scholar] [CrossRef]
  32. Sterling, A.; Rodríguez, C.H. Valoración de Nuevos Clones de Hevea brasiliensis con Proyección Para la Amazonia Colombiana: Fases de pre y Post-Sangría Temprana en el Caquetá; Sterling, A., Rodríguez, C.H., Eds.; Instituto Amazónico de Investigaciones Científicas—SINCHI: Bogotá, Colombia, 2020. [Google Scholar]
  33. Sterling, A.; Rodríguez, C.H. Estrategias de Manejo Para las Principales Enfermedades y Plagas del Cultivo del Caucho con Énfasis en La Amazonia Colombiana; Instituto Amazónico de Investigaciones Científicas—SINCHI: Florencia Caquetá, Colombia, 2018. [Google Scholar]
  34. Sterling, A.; Martínez-Viuche, E.J.; Pimentel-Parra, G.A.; Suárez-Córdoba, Y.D.; Fonseca-Restrepo, J.A.; Virguez-Díaz, Y.R. Dynamics of adaptive responses in growth and resistance of rubber tree clones under South American Leaf Blight non-escape conditions in the Colombian Amazon. Ind. Crops Prod. 2019, 141, 111811. [Google Scholar] [CrossRef]
  35. Verbeke, G.; Molenberghs, G. Linear Mixed Models for Longitudinal Data; Springer: New York, NY, USA, 2000. [Google Scholar]
  36. Rivano, F.; Martinez, M.; Cevallos, V.; Cilas, C. Assessing resistance of rubber tree clones to Microcyclus ulei in large-scale clone trials in Ecuador: A less time-consuming field method. Eur. J. Plant Pathol. 2010, 126, 541–552. [Google Scholar] [CrossRef]
  37. Meenakumari, T.; Meenattoor, J.R.; Thirunavoukkarasu, M.; Vinod, K.K.; Krishan, B.; Gireesh, T.; Thomas, V.; Mydin, K.K.; Jacob, J. Dynamics of long-term adaptive responses in growth and rubber yield among Hevea brasiliensis genotypes introduced to a dry sub-humid climate of eastern India. Ind. Crops Prod. 2018, 119, 294–303. [Google Scholar] [CrossRef]
  38. Da Silva, K.R.; Cecílio, R.A.; Xavier, A.C.; Pezzopane, J.R.M.; de Oliveira Garcia, E. Zoneamento edafoclimático para a cultura da Seringueira No Espírito Santo. IRRIGA 2013, 18, 1–12. [Google Scholar] [CrossRef]
  39. Guyot, J.; Condina, V.; Doaré, F.; Cilas, C.; Sache, I. Segmentation Applied to Weather-Disease Relationships in South American Leaf Blight of the Rubber Tree. Eur. J. Plant Pathol. 2010, 126, 349–362. [Google Scholar] [CrossRef]
  40. Sterling, A.; Suárez, J.C.; Rodríguez, C.H.; Caicedo, D.F.; Salas, Y.M.; Virguez, Y.R. Crecimiento inicial de clones promisorios de Hevea brasiliensis (Willd. Ex A. Juss.) Müll. Arg. en sistema agroforestal en Caquetá, Colombia. Colomb. For. 2015, 18, 175. [Google Scholar] [CrossRef]
  41. Gireesh, T.; Meenakumari, T.; Mydin, K.K. Fast Track Evaluation and selection of Hevea brasiliensis clones from a clonal nursery. Ind. Crops Prod. 2017, 103, 195–201. [Google Scholar] [CrossRef]
  42. Gonçalves, P.D.S.; Bortoletto, N.; Furtado, E.L. Desempenho de Clones de Seringueira da Série IAC 300 Selecionados Para a Região Noroeste Do Estado de São Paulo; Embrapa Secretaria de Pesquisa e Desenvolvimento, Pesquisa Agropecuária Brasileira: Brasília, Brazil, 2001; pp. 589–599. [Google Scholar]
  43. Vrignon-Brenas, S.; Gay, F.; Ricard, S.; Snoeck, D.; Perron, T.; Mareschal, L.; Laclau, J.P.; Gohet, É.; Malagoli, P. Nutrient management of immature rubber plantations. A review. Agron. Sustain. Dev. 2019, 39, 11. [Google Scholar] [CrossRef]
  44. Gouvêa, L.R.L.; Silva, G.A.P.; Verardi, C.K.; Oliveira, A.L.B.; Gonçalves, P. Simultaneous selection of rubber yield and girth growth in young rubber trees. Ind. Crops Prod. 2013, 50, 39–43. [Google Scholar] [CrossRef]
  45. Gonçalves, P.d.S.; Silva, M.d.A.; Gouvêa, L.R.L.; Scaloppi, E.J., Jr. Genetic variability for girth growth and rubber yield in Hevea brasiliensis. Sci. Agric. 2006, 63, 246–254. [Google Scholar] [CrossRef]
  46. Vélez-Sánchez, D.F.; Nieto-Rodríguez, V.M. Desempeño inicial de clones de caucho en el municipio de Puerto Carreño, Vichada, Colombia. Colomb. For. 2013, 16, 186–199. [Google Scholar]
  47. Vinod, K.K.; Rajeswari, M.J.; Pothen, J.; Krishnakumar, A.K.; Sethuraj, M.R. Performance Analysis for Wintering Pattern in Hevea brasiliensis Clones. Ind. J. Nat. Rubber Res. 1996, 9, 44–47. [Google Scholar]
  48. Quesada-Méndez, I.; Aristizábal-Gutiérrez, F. Characterization of two latex parameters from Hevea Brasiliensis (Willd. Ex A. Juss.) Müll. Arg. Clones in the Colombian Altillanura. Colomb. For. 2012, 15, 139–146. [Google Scholar] [CrossRef]
  49. Sethuraj, M.R. Yield Components in Hevea brasiliensis. In Natural Rubber Biology, Cultivation and Technology; Sethuraj, M.R., Mathew, N.M., Eds.; Developments in Crop Science; Elsevier: Kerala, India, 1992; Volume 23, pp. 137–163. [Google Scholar]
  50. Mendoza-Vargas, L.P.; Jiménez-Forero, J.A. Effect of inclusion of anticoagulants during the benefit of Hevea brasiliensis (Willd. Ex a.Juss.) Müll.Arg. Colomb. For. 2018, 21, 91–101. [Google Scholar] [CrossRef]
  51. Ajith, S.V.; Nair, N.N.; Sathik, M.B.M.; Meenakumari, T. Breeding for low-temperature stress tolerance in Hevea brasiliensis: Screening of newly developed clones using latex biochemical parameters. J. Rubber Res. 2021, 24, 819–828. [Google Scholar] [CrossRef]
  52. Cáceres, S.A.P.; Gauthier-Maradei, P. Thermogravimetric analysis as a new method to determine of total solid content (TSC) and dry rubber content (DRC) of natural latex. Rev. ION 2012, 25, 57–65. [Google Scholar]
  53. Támara, A.M.; Reyes, H. Physicochemical analysis and characterization of the Hevea brasiliensis species. Entre Cienc. Ingienería 2015, 18, 9–15. [Google Scholar]
  54. Jacob, J.; Prévöt, J.C.; Roussel, D.; Lacrotte, R.; Serres, E.; d’Auzac, J.; Eschbach, J.M.; Omont, H. Yield limiting factors for latex yield of Hevea. In Physiology of Rubber Tree Latex. The Laticiferous Cell and Latex: A Model of Cytoplasm; D’Auzac, J., Jacob, J.-L., Chrestin, H., Eds.; CRC Press: Boca Raton, FL, USA, 1989; Volume 1, pp. 345–397. [Google Scholar]
  55. Guyot, J.; Condina, V.; Doaré, F.; Cilas, C.; Sache, I. Role of Ascospores and Conidia in the initiation and spread of South American Leaf Blight in a rubber tree plantation. Plant Pathol. 2014, 63, 510–518. [Google Scholar] [CrossRef]
  56. Guyot, J.; Le Guen, V. A Review of a century of studies on South American Leaf Blight of the rubber tree. Plant Dis. 2018, 102, 1–14. [Google Scholar] [CrossRef] [PubMed]
  57. Le Guen, V.; Garcia, D.; Mattos, C.; Fouet, O.; Doaré, F.; Condina, V.; Seguin, M. A newly identified locus controls complete resistance to Microcyclus ulei in the Fx2784 rubber clone. Tree Genet. Genomes 2013, 9, 805–812. [Google Scholar] [CrossRef]
  58. Instituto Geográfico Agustín Codazzi (IGAC) Estudio General de Suelos y Zonificación de Tierras Departamento de Caquetá Escala 1.100.000; Imprenta Nacional de Colombia: Bogotá, Colombia, 2014.
  59. Instituto Geográfico Agustín Codazzi-IGAC Caquetá. Características Geográficas; IGAC: Bogotá, Colombia, 2010. [Google Scholar]
  60. Clément-Demange, A.; Nicolas, D.; Legnaté, H.; Rivano, F.; Le Guen, V.; Gnagne, M.; Chapuset, T. Hévéa: Stratégies de sélection. Plant. Rech. Dev. 1995, 2, 5–14. [Google Scholar]
  61. Sterling, A.; Andrade-Ramírez, T.K.; Clavijo-Arias, E.A.; Guaca-Cruz, L.; Celis-Hernández, M.D.; Noriega, J.E.; Torres-Barrero, R.; Suárez-Salazar, J.C. evaluación de los rasgos de crecimiento, nutricionales y fisiológicos de genotipos caquetenses selección ECC 1 de caucho (Hevea brasiliensis) en campos clonales a gran escala en el Caquetá. In Valoración de Nuevos Clones de Hevea brasiliensis con Proyección Para la Amazonia Colombiana: Fases de pre y Post-Sangría Temprana en el Caquetá; Sterling, A., Rodríguez, C.H., Eds.; Instituto Amazónico de Investigaciones Científicas SINCHI: Bogotá, Colombia, 2020; pp. 127–152. [Google Scholar]
  62. Sterling, A.; Rodríguez, C.H. Nuevos Clones de Cuacho Natural Para la Amazonía Colombiana: Énfasis en la Resistencia al mal Suramericano de las Hojas (Mycrocyclus ilei), 1st ed.; Sterling, A., Rodríguez, C.H., Eds.; Instituto Amazónico de Investigaciones Científicas SINCHI: Bogotá, Colombia, 2011. [Google Scholar]
  63. Confederación Cauchera Colombiana (CCC). Estado Actual Del Gremio Cauchero Colombiano; CAELCA: Bogotá, Colombia, 2015. [Google Scholar]
  64. Silva, G.A.P.; Gouvêa, L.R.L.; Verardi, C.K.; de Oliveira, A.L.B.; Gonçalves, P.d.S. Annual growth increment and stability of rubber yield in the tapping phase in rubber tree clones: Implications for early selection. Ind. Crops Prod. 2014, 52, 801–808. [Google Scholar] [CrossRef]
  65. Lacote, R.; Gabla, O.; Obouayeba, S.; Eschbach, J.M.; Rivano, F.; Dian, K.; Gohet, E. Long-Term Effect of ethylene stimulation on the yield of rubber trees is linked to latex cell biochemistry. Field Crop Res. 2010, 115, 94–98. [Google Scholar] [CrossRef]
  66. Junqueira, N.T.V.; Chaves, G.M.; Zambolim, L.; Gasparotto, L.; Alfenas, A.C. Variabilidade fisiológica de Microcyclus ulei. Fitopatol. Bras. 1986, 11, 823–833. [Google Scholar]
  67. Mattos, C.; Garcia, D.; Pinard, F.; Le Guen, V. Variabilidade de isolados de Microcyclus ulei No Sudeste de Bahia. Fitopatol. Bras. 2003, 5, 502–507. [Google Scholar] [CrossRef]
  68. Pinheiro, J.; Bates, D.; Saikat, D.; Sarkar, D.; R Development Core Team. Nlme: Linear and Nonlinear Mixed Effects Models; R Core Team: Vienna, Austria, 2013; Volume 3, pp. 1–109. [Google Scholar]
  69. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2018. [Google Scholar]
  70. Di Rienzo, J.A.; Casanoves, F.; Balzarine, M.G.; Gonzales, L.; Tablada, M.; Robledo, C.W. InfoStat, versión 2018; National University of Cordoba: Córdoba, Argentina, 2018.
  71. Doledec, S.; Chessel, D. Co-inertia Analysis: An alternative method for studying species–environment relationships. Freshw. Biol. 1994, 31, 277–294. [Google Scholar] [CrossRef]
  72. Thioulouse, J.; Chessel, D.; Dolédec, S.; Olivier, J.-M. ADE-4 a Multivariate analysis and graphical display software. Stat. Comput. 1997, 7, 75–83. [Google Scholar] [CrossRef]
Figure 1. Monthly precipitation, T-max and T-min (°C): maximum and minimum temperature; RH (%): relative humidity (San Vicente del Caguán, Caquetá, Colombia) between July 2017 and June 2020.
Figure 1. Monthly precipitation, T-max and T-min (°C): maximum and minimum temperature; RH (%): relative humidity (San Vicente del Caguán, Caquetá, Colombia) between July 2017 and June 2020.
Plants 12 03627 g001
Figure 2. Representation of the attack of P. ulei on H. brasiliensis. (a) AT1 represents a minor attack (1–5%). (b) AT2 represents a moderate attack (6–15%). (c) AT3 represents a severe attack (16–30%). (d) AT4 represents a very severe attack (>30%).
Figure 2. Representation of the attack of P. ulei on H. brasiliensis. (a) AT1 represents a minor attack (1–5%). (b) AT2 represents a moderate attack (6–15%). (c) AT3 represents a severe attack (16–30%). (d) AT4 represents a very severe attack (>30%).
Plants 12 03627 g002
Figure 3. Results of susceptibility to SALB over three years after planting. The following parameters were evaluated: (a) attacks of P. ulei on young leaves at stage C (AT1); (b) type of reaction (TR) to P. ulei on young leaves at stage C; (c) attacks of P. ulei on mature leaves at stage D (AT2); (d) stromal density on mature leaves at stage D (ST). Means between the genotypes followed by the same letter in the same year do not differ statistically (Fisher’s LSD test, p < 0.05).
Figure 3. Results of susceptibility to SALB over three years after planting. The following parameters were evaluated: (a) attacks of P. ulei on young leaves at stage C (AT1); (b) type of reaction (TR) to P. ulei on young leaves at stage C; (c) attacks of P. ulei on mature leaves at stage D (AT2); (d) stromal density on mature leaves at stage D (ST). Means between the genotypes followed by the same letter in the same year do not differ statistically (Fisher’s LSD test, p < 0.05).
Plants 12 03627 g003
Figure 4. Co-inertia analysis between SALB resistance of 10 genotypes of Hevea brasiliensis: (a) PCA (Principal Component Analysis); (b) captured inertia (%); (c) genotype ordination plot; (d) growth and early yield projection; (e) SALB resistance projection.
Figure 4. Co-inertia analysis between SALB resistance of 10 genotypes of Hevea brasiliensis: (a) PCA (Principal Component Analysis); (b) captured inertia (%); (c) genotype ordination plot; (d) growth and early yield projection; (e) SALB resistance projection.
Plants 12 03627 g004
Table 1. Selection criteria for the longitudinal data analysis model for assessing the growth and resistance to SALB.
Table 1. Selection criteria for the longitudinal data analysis model for assessing the growth and resistance to SALB.
ModelCompound Symmetric Model with Heterogeneous Variances (CSH)First-Order Autoregressive Model with Heterogeneous Variances (ARH1)First-Order Continuous Autoregressive Model with Heterogeneous Variances (CARH1)Unstructured Model with Heterogeneous Variances (UNH)
VariableLog LikAICBICLog LikAICBICLog LikAICBICLog LikAICBIC
CTpt−175.38420.77508.26−170.39410.77498.27−170.39410.77498.27−165.18404.35496.84
AT1159.78−235.6−130.6159.97−236−131159.97−235.95−131161.4−234.8−124.81
AT263.18−56.3631.1363.29−56.5930.9163.29−56.5930.9164.22−54.4438.06
ST48.8−27.5959.948.68−27.3760.1348.68−27.3660.1349.65−25.367.2
TR97.9−125.8−38.3197.1−124.2−36.797.1−124.19−36.7101.28−128.57−36.07
Log Lik, log-likelihood; AIC, Akaike criterion; BIC, Bayesian criterion [35]; CTpt, circumference of the trunk in the pre-tapping phase; AT1, attack severity on young leaves; AT2, attack severity on mature leaves; TR, type of reaction on young leaves (conidial sporulation intensity, TR); ST, stroma density on mature leaves [31,34,36].
Table 2. Analysis of variance (ANOVA) in growth, SALB resistance, and early yield for effects of genotype, year, and genotype × year interaction.
Table 2. Analysis of variance (ANOVA) in growth, SALB resistance, and early yield for effects of genotype, year, and genotype × year interaction.
VariableDF1DF2DF3GenotypeYearGenotype × Year
CTpt92184.42 (0.0001)714.72 (<0.0001)2.23 (0.0074)
AT192185.35 (<0.0001)89.39 (<0.0001)2.98 (0.0004)
AT2921816.19 (<0.0001)256.90 (<0.0001)2.68 (0.0012)
ST921819.30 (<0.0001)215.16 (<0.0001)3.08 (0.0002)
TR92183.50 (0.0010)116.21 (<0.0001)2.09 (0.0125)
DRC9--10.14 (<0.0001)--
DRY9--18.85 (<0.0001)--
CTo9--5.69 (<0.0001)--
CTpt, circumference of the trunk in the pre-tapping phase; AT1, attack severity on young leaves; AT2, attack severity on mature leaves; TR, type of reaction on young leaves (conidial sporulation intensity, TR); ST, stroma density on mature leaves [31,34,36]. DRC, dry rubber content; DRY, dry rubber yield; CTo, circumference of the trunk at opening [31,37]; DF1, degrees of freedom for the factor genotype; DF2, degrees of freedom for the factor year; DF3, degrees of freedom for interaction genotype × year; -, not applicable.
Table 3. Annual circumference of the trunk (cm) in the pre-tapping phase (CTpt) for 10 genotypes of Hevea brasiliensis.
Table 3. Annual circumference of the trunk (cm) in the pre-tapping phase (CTpt) for 10 genotypes of Hevea brasiliensis.
GenotypeYear 1Year 2Year 3
ECC 259.37 ± 0.49 a,b k18.33 ± 0.88 def25.75 ± 1.24 ab
ECC 299.43 ± 0.49 k16.98 ± 0.88 fgh22.80 ± 1.24 bc
ECC 357.78 ± 0.49 l13.33 ± 0.88 i20.23 ± 1.24 cde
ECC 609.03 ± 0.49 k14.89 ± 0.88 hi21.44 ± 1.24 c
ECC 6410.69 ± 0.49 j20.44 ± 0.88 cd27.48 ± 1. 24 a
ECC 668.73 ± 0.49 kl15.31 ± 0.88 ghi22.54 ± 1.24 bc
ECC 739.44 ± 0.49 k16.66 ± 0.88 fgh22.84 ± 1.24 bc
ECC 839.62 ± 0.49 jk17.31 ± 0.88 efg23.61 ± 1.24 bc
ECC 909.02 ± 0.49 k15.83 ± 0.88 gh23.63 ± 1.24 bc
IAN 8739.82 ± 0.49 jk16.28 ± 0.88 fgh22.41 ± 1.24 bc
a Standard error; b Means in each column with the same letter do not differ statistically (Fisher’s LSD test, p < 0.05).
Table 4. Circumference of the trunk (cm) at opening (time of commencement of tapping or 3rd year) (CTo); DRY, dry rubber yield; DRC, dry rubber content in the early tapping phase for 10 genotypes of Hevea brasiliensis.
Table 4. Circumference of the trunk (cm) at opening (time of commencement of tapping or 3rd year) (CTo); DRY, dry rubber yield; DRC, dry rubber content in the early tapping phase for 10 genotypes of Hevea brasiliensis.
GenotypeCTo (cm)DRY (g.tree−1.tap−1)DRC (%)
ECC 2526.86 ± 0.35 a,b b5.57 ± 0.43 a32.45 ± 0.44 a
ECC 2926.49 ± 0.37 b0.74 ± 0.18 c31.44 ± 2.29 a
ECC 3525.00 ± 2.45 b0.42 ± 1.14 c26.24 ± 0.76 b
ECC 6025.03 ± 0.93 b1.21 ± 0.43 c30.43 ± 0.70 ab
ECC 6427.92 ± 0.26 a3.76 ± 0.66 b31.38 ± 1.01 a
ECC 6625.54 ± 0.78 b1.09 ± 0.43 c30.75 ± 1.10 a
ECC 7325.62 ± 0.56 b2.23 ± 0.43 bc32.42 ± 1.18 a
ECC 8326.09 ± 0.45 b3.30 ± 0.24 b30.53 ± 1.98 ab
ECC 9026.94 ± 0.57 b2.96 ± 0.81 bc28.83 ± 3.05 ab
IAN 87325.47 ± 0.60 b1.88 ± 0.28 c32.83 ± 0.10 a
a Standard error. b Means in each column with the same letter do not differ statistically (Fisher’s LSD test, p < 0.05).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Galindo-Rodríguez, L.C.; Sterling, A.; Muñoz-Ramirez, H.; Fonseca-Restrepo, J.A. Performance- and Resistance-Related Early Responses of Colombian Elite Rubber Tree Genotypes under Low Pressure of South American Leaf Blight: Implications for Disease Management in the Amazon. Plants 2023, 12, 3627. https://doi.org/10.3390/plants12203627

AMA Style

Galindo-Rodríguez LC, Sterling A, Muñoz-Ramirez H, Fonseca-Restrepo JA. Performance- and Resistance-Related Early Responses of Colombian Elite Rubber Tree Genotypes under Low Pressure of South American Leaf Blight: Implications for Disease Management in the Amazon. Plants. 2023; 12(20):3627. https://doi.org/10.3390/plants12203627

Chicago/Turabian Style

Galindo-Rodríguez, Lyda Constanza, Armando Sterling, Herminton Muñoz-Ramirez, and Jesica Andrea Fonseca-Restrepo. 2023. "Performance- and Resistance-Related Early Responses of Colombian Elite Rubber Tree Genotypes under Low Pressure of South American Leaf Blight: Implications for Disease Management in the Amazon" Plants 12, no. 20: 3627. https://doi.org/10.3390/plants12203627

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop