Managing Potato Biodiversity to Cope with Frost Risk in the High Andes: A Modeling Perspective

Austral summer frosts in the Andean highlands are ubiquitous throughout the crop cycle, causing yield losses. In spite of the existing warming trend, climate change models forecast high variability, including freezing temperatures. As the potato center of origin, the region has a rich biodiversity which includes a set of frost resistant genotypes. Four contrasting potato genotypes –representing genetic variability- were considered in the present study: two species of frost resistant native potatoes (the bitter Solanum juzepczukii, var. Luki, and the non-bitter Solanum ajanhuiri, var. Ajanhuiri) and two commercial frost susceptible genotypes (Solanum tuberosum ssp. tuberosum var. Alpha and Solanum tuberosum ssp. andigenum var. Gendarme). The objective of the study was to conduct a comparative growth analysis of four genotypes and modeling their agronomic response under frost events. It included assessing their performance under Andean contrasting agroecological conditions. Independent subsets of data from four field experiments were used to parameterize, calibrate and validate a potato growth model. The validated model was used to ascertain the importance of biodiversity, represented by the four genotypes tested, as constituents of germplasm mixtures in single plots used by local farmers, a coping strategy in the face of climate variability. Also scenarios with a frost routine incorporated in the model were constructed. Luki and Ajanhuiri were the most frost resistant varieties whereas Alpha was the most susceptible. Luki and Ajanhuiri, as monoculture, outperformed the yield obtained with the mixtures under severe frosts. These results highlight the role played by local frost tolerant varieties, and featured the management importance –e.g. clean seed, strategic watering- to attain the yields reported in our experiments. The mixtures of local and introduced potatoes can thus not only provide the products demanded by the markets but also reduce the impact of frosts and thus the vulnerability of the system to abiotic stressors.


Introduction
The Altiplano is a high tropical plateau located at 3600-4300 m above sea level in the Andes of Bolivia and Peru. Most of the cropland is located below 4000 masl; above that elevation land is mainly covered by natural grasslands and is only used for growing bitter potato landraces, which are adapted to cold conditions. Potato is by far the most important crop in the region, accounting for 44% of the gross value of crop production [1] from a cropping area of about 88,000 ha [2]. Potato production is limited by abiotic and biotic factors; Andean farmers manage these constraints mainly by the use of a high diversity of native species and cultivars that are often grown as mixtures in single plots [3,4,5]. As potato originated in the Andes [6], local genetic diversity in cultivated potato is large and includes several species, comprising both bitter -Solanum juzepckzukii (triploid), and l nonbitter frost resistant potatoes: Solanum ajanhuiri (diploid), but also the non-bitter frost susceptible conventional Solanum tuberosum subspecies tuberosum (tetraploid), and Solanum tuberosum ssp. andigenum (tetraploid), which are present in the Altiplano [7,8].
The principal role played by the diversity of potatoes grown in the Altiplano is related to smallholder's food security. Potato fresh yields in the area are low. In Peru and the northernmost part of the Bolivian Altiplano, yield average varies from 4 to 5.2 t/ha whilst in the southern Bolivian section [1] the average yield is 3.6 t/ha. The growing season in the Altiplano extends from October to March, when maximum annual temperature coincides with the rainy season. In the agricultural zones of the Altiplano, average maximum temperature is around 18uC whereas minimum temperature is around 4uC during the growing season. Precipitation is around 800 mm/year in the northeast of the Altiplano whereas in the southwestern Altiplano, it is about 200 mm/year, mostly occurring during the same growing season. Production risk for potato is high due to several recurrent factors, particularly drought, hail, and frost. Frost-free period averages 140 days in the northern Altiplano and 110 days in the Southern areas [9]. The high production risks presented by frost and other factors may also lead to reduced investment in agriculture, resulting in low production which in turn affects food availability.
The varieties of the species Solanum tuberosum ssp. andigenum are the most widely cultivated in the Andes. The Solanum juzepczukii stands out for its high frost and drought tolerance and its capacity to grow at 4000 masl and above [10]. However its tubers are bitter due to a high content of glycoalkaloids, requiring processing for direct human consumption [11]. This processing is an old Andean's strategy for conserving food -chuñ o: dehydrated potatoes -for several years [12]. It has been estimated that at least 25% of total area under potato in the Altiplano is planted with bitter varieties [13]. This assertion is supported by the estimate that bitter potatoes make up 15% of total potato area in Bolivia [11]. In fact, more bitter potatoes are found in the Altiplano than in most other zones. The varieties of Solanum ajanhuiri have characteristics of tolerance to frost and drought similar to those of the Solanum juzepczukii but they do not have high glycoalkaloid contents and are therefore non bitter. At intermediate and lower altitudes the Solanum tuberosum ssp. tuberosum varieties are the most widespread whereas the other species predominate at altitudes higher than 3500 masl. Notwithstanding, recent findings showed that the actual upper limit for all potato varieties has increased to around 4300 masl [14] in response to increased temperatures and disease pressure at lower altitudes brought about by climate change.
The Andes represent the largest and highest mountain range in the tropics, and thus a suitable ecosystem to study changes in climate and how they affect the natural resources and the livelihood they sustain. Recent studies based on local meteorological networks [15,16] have shown a significant warming trend after 1979 (0.32-0.34uC/decade). This warming trend shows some interesting features. Below 1,000 m there is a substantial difference between the eastern (Amazonian) and the western (Pacific watershed) slopes i.e. no warming trend on the eastern slope whereas on the western side there is a warming trend of 0.39uC/ decade. On this flank, there is an almost linear decrease of the warming trend with altitude, reaching down to 0.16uC/decade above 4000 m. Thus there are differences between the Eastern and Western slopes and a definite vertical structure. This behavior differs from other mountain ranges in other latitudes e.g. European Alps and Tibet. A warming trend is also projected for the future [17], but differences as functions of position and altitudes are yet to be studied.
Summer frost events are caused by radiative cooling and are common throughout the Andean highlands and can occur at any time during the growing season [9,18]. Frost can cause partial or complete loss of foliage, leading to a reduction in photosynthate production and hence yield. In turn, crop failure caused by frost damage may lead to a decrease in the total area planted to potato in the subsequent season due to seed shortage [18]. The temperature at which leafage frost damage occurs depends on the species and the cultivar. For Solanum tuberosum subsp. andigenum, frost damage is likely to occur when the temperature drops to 2uC or lower [19]. Higher frost resistance exists in other cultivars and in wild potato species. For example, cultivated potato species such as S. ajanhuiri and S. curtilobum are damaged at 23 to 25uC, whereas S. juzepzuckii generally resists temperatures down to 25uC and perhaps even lower [20].
Crop growth simulation models can be used to analyze constraints and opportunities for crop yields in complex production systems (e.g. [4,5]). One of the limitations for using simulation models in this context is that most available models are calibrated only for varieties of S. tuberosum subspecies. Almost all modeling work has been done for so-called 'modern' varieties from formal breeding programs; with a few exceptions [21]. This can be particularly problematic when working in areas of high genetic diversity, where landraces with rather distinct growth characteristics may coexist with modern cultivars. A more in-depth understanding of the most important crop growth limiting factors as affected by genetic diversity is needed to arrive at more robust recommendations to reduce the vulnerability and improve the productivity of potato in environmentally challenged complex cropping systems. Comparative growth analysis and its translation into potato growth models can constitute a useful tool for scientists, particularly under climate variability and global climate change conditions in areas of rich genetic diversity.
The objective of the present study was to conduct a comparative growth analysis of four contrasting genotypes from the species S. juzepczukii, S. ajanhuiri and S. tuberosum subspecies andigenum and tuberosum and translate it into a modeling of their agronomic response. It included assessing their performance under contrasting agroecological conditions in the high Andes as well as the impact of frost events.

Materials and Methods
This study did not involve neither human nor animal subjects. Field experiments were conducted in research stations and farmer fields and no law regulated permits were required because all materials used are commonly used by farmers with no environmental negative impacts. No protected species were sampled. Permission to work on PROINPA's field research sites (Patacamaya 2, Patacamaya 3, Puchuni and Laurani) was granted by Mr. Enrique Carrasco, then Leader of PROINPA.
Permission to work on farmers' fields was granted by Mr. Vitaliano Mamani, Mrs. Maria Laura and Mr. Juan Carlos Huanca, under a collaboration agreement with PROINPA.
All authorizations were granted as verbal agreements on a bona fide basis. Germplasm Four contrasting germplasm (species, cultivars, varieties, and landraces) -three native landraces and one introduced cultivarrepresenting the genetic diversity cultivated in the high Andes were studied. The four genotypes were Bola luki, referred to as Luki in this paper, a bitter potato variety of Solanum juzepczukii (36 or triploid); Chiar ajanhuiri, referred to as Ajanhuiri, a diploid (26) variety of Solanum ajanhuiri; Gendarme, a tetraploid (46) variety of Solanum tuberosum ssp. andigenum; and Alpha, a European tetraploid variety of Solanum tuberosum ssp. tuberosum [7,22].

Experiments and Data Collection
Four field trials, comparing four potato genotypes, were carried out in farmer's fields or on experimental stations in different locations (Patacamaya 2, Patacamaya 3, Puchuni and Laurani) in high and semi-arid zones in Bolivia (Table 1), where potato is the dominant crop, grown as the first crop after a fallow period. Within each site, complete randomized block designs with 3 or 4 replications were implemented in plots of 25.2 m 2 . The planting density was 4.76 plants/m 2 , a density suitable for the potato plow used in the zone. The tuber seeds used were from certified quality and were stored homogeneously during five months before planting. In addition, another experiment was conducted at Wichukollu locality, where the effect of frost (22.51uC befall at 94 days after planting) on the four studied varieties was assessed.
Plots were homogeneously managed to assure non-limiting factors for achieving potential production. Based on previous fertilization trials in the region [23], N was applied at 60-80 kg/ha and P 2 O 5 at 100-120 kg/ha, in addition to 5 t/ha of bovine manure. Potassium was not applied due to the high levels of this element in these soils. Crops were rain-fed, with supplemental irrigation to bring the soil to field capacity when soil moisture dropped below 75% of field capacity, as determined by weekly soil samplings. Ridomil MZ (metalaxyl-M and mancozeb) and Bravo 500 (chlorothalonil) were used against potato late blight (caused by Phytophthora infestans Mont. de Bary). Synthetic pyrethroids, Karate (lambda-cyhalothrine) and Lorsban (chlorpyrifos) were used to control Andean potato weevils (Premnotrypes latithorax Piercei, P. solaniperda Kuschel, and Rigopsidius tucumanus Heller) and potato tuber moths (Phthorimaea operculella Zeller, Symmetrischema tangolias Geyen and Paraschema detectendum Povolny).  Precipitation, minimum and maximum temperatures were measured in the field using Hobo data loggers (ONSET Pro Series H08-032-08, USA). Global solar radiation was measured at meteorological stations representative of the experimental sites. Photosynthetically active radiation (PAR) was calculated as half of the global solar radiation [24].
During the crop growing period, plants were sampled three to five times from the plots. Intermediate harvests consisted of 4 plants per plot, while the final harvest consisted of 20 plants. Biomass was divided into leaves, stems, and tubers. Roots were not measured since they are difficult to collect. For each biomass group, total fresh weight was determined and their dry matter fraction estimated using samples of 150 g fresh weight subsequently dried at 75uC until constant weight which occurred after 48 to 72 hours.
The fraction of green canopy cover was estimated in triplicate in each experimental unit, every first or second week, from emergence to the final harvest. A wooden frame of 70 cm690 cm divided into 100 cells of 7 cm69 cm was used. The frame covered the area of three plants and cells were marked as covered if leaf surfaces occupied more than 50% of the cell area.

Data Processing
The analysis of variance for the comparative growth analysis of all four genotypes was conducted by pooling the results from the four locations, once the homogeneity of the variance across experiments was demonstrated with the Bartlett test. The crop growth responses tested included: Maximum Canopy Cover (MCC), Light Use Efficiency (LUE), Dry Tuber Yield (DTY) and Harvest Index (HI). A summary of the variables and parameters of this study is described in Table 2.   The linear additive model used was: Growth analysis. Field measurements of plant growth components such as leaf area and the weight of plant parts provided for the parameters required for growth analysis. Maximum canopy cover MCC was calculated from the observed data through regression analysis. Growth parameters such as the rate of relative increase of light interception (Ro), the initial fraction of light interception at plant emergence (Fo), and the time when the fraction of light intercepted was reduced to 50% (t 0.5 ) were estimated using curve fitting. The tuber initiation onset (TIO), the maximum tuber growth or bulking rate (TG max ), and cessation of tuberization (TM) were calculated by fitting a logistic function (Sigma Plot V.9, USA) to the measured tuber weight as a function of thermal time (uCd). The average LUE was determined as the slope of the curve between cumulative total dry matter and intercepted PAR [24].
Model description and calibration. Results from growth analysis can be integrated into mathematical models of the growth process, which constitute a robust tool for yield forecasting. This basic approach was followed in the present study to assess the response of the four potato genotypes and the results were used to calibrate the mathematical equations included in a slightly modified version of the LINTUL-potato model [24,25,4]. The model's routines were originally developed in Visual Basic of Microsoft Excel and, once validated, programmed in C++. The model operates on a daily time step and its equations and parameters are described in [21]. The results of two experiments (Patacamaya 2 and Puchuni) were selected to calibrate the model and the results of the remaining two experiments were used for validation.
For modeling purposes, single plot variety mixtures were constructed, based on common practice by farmers. Our simulated mixtures included 65% Gendarme, 25% Luki, 7% Ajanhuiri and 3% Alpha [11,13]. The yields from the mixtures in single-plots were compared against monoculture yield for each variety.
Model validation. Data sets from trials conducted in the other two independent sites (Patacamaya 3 and Laurani) were used to validate the calibrated model. Simulated data for MCC and DTY were compared with field results. Mean Bias Error (MBE), the Root Mean Square Error (RMSE), the Nash coefficient (Nash), and R 2 [26,27] were used to test the suitability of the calibrated model to simulate the growth of different genotypes under the experimental conditions.
Regression diagnostics of modeling. Regression diagnostics is the general class of techniques for detecting problems -with either the model or the data set -in regression analysis [28]. We adapted the residual analysis to evaluate simulated canopy cover data in time. This analysis was portrayed by plotting the phenological time series and residual values for detecting problem areas with outliers. The Durbin-Watson test was used to detect the presence of serial correlation in the residuals and influence statistics of points was measured through Cook's D, DFFITS, DFBETAS and COVRATIO statistics [28]. The time series plotting was made in Sigma Plot V.9 software and the statistical analysis in SAS V.8.
Scenario analysis. What if type scenarios were constructed based on the validated model into which a frost routine, based on the data coming from the Wichukollu experiment, was incorporated to estimate potato yield losses for every genotype [29,4]. Two issues were addressed with the scenario analysis: 1) How would different genotypes behave when exposed to progressive frosts (increasing by 21uC) during any of the four critical (30, 60, 90 and 120 DAE) phenological phases? 2) How would mixtures of varieties behave under extreme frosts, of frequent occurrence in the high Andes (23uC at 60 and 90 DAE), compared to the monocultures? The model was set to simulate expected impact of progressive and extreme frosts on yield and it was compared to the production without frost stress.

Growth Analysis
Harvest index was the only response variable affected (P,0.05) by location (Table 3). Average LUE, MCC and DTY did not differ (P.0.05) among locations. The highest HI was recorded in Puchuni whereas the lowest HI occurred in Patacamaya 2. The average yield across all locations was six to sevenfold higher than the national average. MCC, DTY and HI differed (P,0.001) among genotypes ( The calculated values of the light use efficiency, (around 2.7 g DM MJ 21 as shown in Table 3) measured as the slope of the line describing the relationship between biomass production and absorbed PAR, were within the range of values reported in the literature for the tuberosum genotypes under temperate conditions [24,30,31].
Harvest index is generally regarded as the most important variable to estimate the productivity efficiency of the crops [32]. Generally the potato crop has higher HI values (about 0.80) than the other important crops (e.g. oil seeds 0.3 and cereals 0.5). In our case, HI was $0.8. However, in spite of the genetic diversity included in the study, the correlation with tuber yield was low, contrasting with literature findings [10]. This was probably due to MCC values that did not reach the level of a full coverage (MCC,1), which is likely due to the low F 0 and slow R 0 values in Andean genotypes, that might be explained by the low temperatures to which they are usually exposed during their growth in the Andean highlands. This poor value of the coverage level in the high Andes associated to its short duration contributes to a low biomass production [30,33].

Genotype by Environment Interaction
DTY (P,0.01) and MCC (P,0.05) were significant for genotype x environment interactions. Figure 1, shows the performance of each genotype in the four experiments in which frost was not considered, against the average performance of all genotypes in each contrasting environment defined by cumulated rainfall as dry, intermediate and moist condition. Under all conditions, Gendarme outperformed the average yield whereas Ajanhuiri outperformed the average yield in dry and intermediate conditions only but was negatively affected by moist conditions. Alpha yielded much less than the average in all environments, a behavior explained by the overall low MCC. Luki was not influenced by environments and its yield was similar to the average.    (Table 4). Alpha showed the highest (5610 23 ) light interception at emergence F 0 value whereas the lowest value was shown by Luki (2610 23 ). Rate of increase of light interception R 0 was highest for Gendarme (7.3610 23 uCd 21 ), and lowest for Luki (5.5610 23 uCd 21 ). The duration of senescence (d) was longest for Alpha (approximately 38 days after the attainment of MCC) and shortest for Ajanhuiri (23 days following maximum canopy cover MCC). In the Bolivian Andes, Alpha is clearly an early maturing type [time of reduction to 50% of light interception t 0.5 = 106 Days after Planting (DAE) and tuber initiation onset TIO = 24 DAE] while Luki is a very late maturing with a t 0.5 of 120 DAE and TIO of 60 DAE ( Table 4).
The growth and developmental differences between Alpha and Luki genotypes is further evidenced by the time when maximum bulking rate TG max is attained; 60 DAE and 90 DAE for Alpha and Luki, respectively. The intercept and the slope were equal (P.0.05) to 0 and 1, respectively. The tuberization dynamics and parameter values for all genotypes are shown in Table 5.
The R 2 values (Table 6) show that the model explained more than 81% of the variations in DTY obtained experimentally for all varieties. For Luki and Ajanhuiri, the model explained more than 81 and 88% of the variation in the measured variable, respectively. For Alpha it explained more that 89% of the variation whereas for Gendarme, the model performed quite well, with a R 2 value of 96% for DTY. All other statistical metrics used to test the validity of the model (MB, RMSE and Nash) to simulate potato yield and yield components under the conditions encountered in the high Andes also showed the model's robustness ( Table 6).
The S. tuberosum ssp. andigenum, S. ajanhuiri and S. juzepczukii seem to be well fitted for the harsh and variable environmental conditions prevalent in the high Andes, as suggested by all varieties. Concerning frost, differences in frost tolerance among the tested germplasm were evidenced by the differences in canopy cover retention following a frost episode. Figure 2 shows that Luki was not affected by frost; the effect on Ajanjuiri was minimal. On the other hand, Alpha was severely affected by the same frost episode. Nonetheless, genotypes from the subspecies tuberosum might play a role for specific objectives and climatic conditions, as suggested by Alpha performance under less stressful conditions. The genetic diversity of potatoes in the high Andes and their different responses to stresses suggests that potato production there can be highly increased with an adequate choice of germplasm, crop management and seed quality [34,33].
Tuberization dynamics provides a clear differentiation among the genotypes studied [21]. Earliness defined as tuber initiation TIO is one of the growth parameters that express differences among genotypes. The earliness of Alpha was probably prompted by environmental conditions in the experimental sites where mean temperature was around 12uC with around 12 h of photoperiod as it has been reported that short photoperiods and low mean temperatures trigger an early onset of tuber growth and bulking in S. tuberosum ssp. Tuberosum [35].
It seems that the simultaneous occurrence of early tuber bulking and fast bulking rate traits is of high importance in the environmental conditions found in the high Andes. In spite of its earliness, the slow bulking rate of Alpha was associated with its comparative lower yields. On the contrary, Gendarme with its intermediate bulking initiation but fast-bulking rates produced the highest yield. Intermediate tuber initiation TIO combined with intermediate bulking rates produced intermediate yields.
Alpha accumulates assimilates earlier thus relying on a mechanism for escaping frost damage provided it happens late in the growing season. On the other hand, Luki presents a delayed TIO and a slow tuberization rate and thus extreme frost events (temp,25uC) can disrupts the accumulation of assimilates and might cause losses.
The juzepczukii (var. Luki) and ajanhuiri genotypes outperformed the yield obtained with the mixtures under severe frosts. These results highlight the importance of the role played by local varieties tolerant to frost and drought stresses. Since the simulations were based on models parameterized with experimental results, yields obtained are feasible if good quality seeds and adequate management are implemented.

Modeling Diagnostics
The residual analysis showed a good spread of the values throughout the growth cycle. Luki presented a few points out of the accepted confidence envelop during the MCC and senescent phases. This was probably due to the rosette-type morphology which makes this variety susceptible to flattening (Figure 3).
For the emergence, maximum growth rate and MCC phases, the residual values were distributed within two standard deviations from the mean. Nonetheless, the regression diagnostics (Cook's D, DFFITS, DFBETAS and COVRATIO), showed flagging points during senescence. This might be due to the lack of accuracy of the visual estimation of canopy cover, as the human eye cannot differentiate photosynthetically active leaves from senescent ones, particularly at the early stages of senescence. The use of digital cameras and appropriate segmentation techniques can reduce these problems (unpublished results from Production Systems and Environment subprogram at CIP).

Scenario Analysis
The validated model was used to run what if type scenarios, to assess the likely impact of frost events -at specific phenologic stages -on yields. A threshold for absolute growth cessation with no  (Figure 4). It is noteworthy that when the frost event occurs at 60 DAE, there is a drastic yield reduction in all genotypes. This seems to be associated to the fact that this is the time when canopy maximum growth rate and tuber initiation are also occurring. Alpha was the most affected with losses increasing exponentially as the severity of the frost augmented. Total loss in Alpha was caused by a 25uC frost. Gendarme presents a minimal damage at 21uC (0.2 t/ha), important losses at 23uC (9.5 t/ha) and a severe loss of tuber production at 25uC (14 t/ha). Ajanhuiri, in turn, was not affected at 21uC, presented a minor reduction at 23uC (2.7 t/ha) and an important reduction of tuber production at 25uC (10.3 t/ha). Luki was not affected at 21 and 23uC, but at 25uC was almost devastated ( Figure 5).
We also simulated total yield comparing monoculture and single plot mixtures containing fractions of total acreage with four varieties distributed in the field (Alpha at 0.03, Gendarme at 0.65, Ajanhuiri at 0.07, and Luki at 0.25) the proportions mimicking the actual genotype distributions encountered in farmer fields. Scenarios without and with frost occurrence (23uC at 60 and 90 DAE) were compared (Table 7). Frost free seasons resulted in average yields of 13 t/ha for the mixture. Gendarme in monoculture out yielded the simulated mixture. On the other hand, the mixture yielded 5.3 t/ha when affected by frost simulations. This yield was higher than the resulting from the commercial varieties Alpha and Gendarme but lower than those produced by Ajanhuiri and Luki as monocultures. If the tendency to a reduced planting of frost tolerant genotypes (Ajanhuiri and Luki) is confirmed, the production based on commercial potatoes can attain only 3.5 t/ha (Gendarme) when frosts occur (Table 7), thus local frost tolerant varieties provide an insurance under uncertain climate.
The varietal single plot mixtures used by the high Andes farmers seem to be a robust strategy to cope with climatic risks. This is particularly true under uncertain climate when a portfolio of options is warranted. The inclusion of native potatoes especially bitter varieties (S. juzepczukii) in the mix is a must under highly variable climatic conditions. The scenarios tested were purposely limited to fixed frosts at specific times during the growing period, just as examples of the modeling possibilities. The reality is that frosts are likely to occur randomly at different times thus highlighting the possible role of mixing varieties -as practiced for millennia by local farmers -to cope with those extreme events. Although the projection of future climate lean towards warmer nights, the increased variability of both night temperature and precipitation together with increased radiation cooling [17] warrant the inclusion of the drought and frost tolerance traits of juzepczukii and ajanhuiri genotypes. Therefore, the genetic diversity must be maintained since temperature and rainfall variability is expected to increase.

Conclusions
The correspondence between the simulations and the results of the independent set of experiments conducted under various agroecological conditions demonstrated the adequacy of the model. This robustness applies not only to yield results under non-limiting factors but also to the forecast of the effects of frost events on canopy cover and correlated yield, which is a plus of the model. All statistical metrics used to test the validity of the model to simulate potato yield and yield components under the conditions encountered in the high Andes showed the model's robustness.
The results indicate that Ajanhuiri and Gendarme varieties seem to be well suited for conditions with intermediate moisture levels, whereas, Ajanhuiri showed a better adaptation to drier areas. Ajanhuiri is a non-bitter variety tolerant to frost. Reintroducing good quality seeds into high Andean farming systems could enhance yields given its hardiness to cope with frequent frosts. For moister -rainfed or with strategic irrigation -areas the clear recommendation seems to be Gendarme and the like genotypes. Luki, on the contrary seems to be insensitive to environmental variations within the Altiplano.
Current conditions in the high Andes include a prolonged growing period of around 180 days. Small farmers diversify their risk through an assortment of genotypes in single plots. The observed diversity of tuberization dynamics of different genotypes supports this risk management strategy given the climatic and altitudinal variability and the low opportunity cost of the land. As the access to land decreases and direct as well as indirect effects of climate change impinge on the crop, new strategies must be sought. One such a strategy could be to develop genotypes with an early TIO and a very fast early bulking rate, able to escape from early biotic or abiotic shocks while sustaining a sizable yield.