Optimizing imputation of marker data from genotyping-by-sequencing (GBS) for genomic selection in non-model species: Rubber tree (Hevea brasiliensis) as a case study.

Genotyping-by-sequencing (GBS) provides the marker density required for genomic predictions (GP). However, GBS gives a high proportion of missing SNP data which, for species without a chromosome-level genome assembly, must be imputed without knowing the SNP physical positions. Here, we compared GP accuracy with seven map-independent and two map-dependent imputation approaches, and when using all SNPs against the subset of genetically mapped SNPs. We used two rubber tree (Hevea brasiliensis) datasets with three traits. The results showed that the best imputation approaches were LinkImputeR, Beagle and FImpute. Using the genetically mapped SNPs increased GP accuracy by 4.3%. Using LinkImputeR on all the markers allowed avoiding genetic mapping, with a slight decrease in GP accuracy. LinkImputeR gave the highest level of correctly imputed genotypes and its performances were further improved by its ability to define a subset of SNPs imputed optimally. These results will contribute to the efficient implementation of genomic selection with GBS. For Hevea, GBS is promising for rubber yield improvement, with GP accuracies reaching 0.52.

Gene deletion of Corynespora cassiicola cassiicolin Cas1 suppresses virulence in the rubber tree.

Corynespora cassiicola is an ascomycete fungus causing important damages in a wide range of plant hosts, including rubber tree. The small secreted protein cassiicolin is suspected to play a role in the onset of the disease in rubber tree, based on toxicity and gene expression profiles. However, its exact contribution to virulence, compared to other putative effectors, remains unclear. We created a deletion mutant targeting the cassiicolin gene Cas1 from the highly aggressive isolate CCP. Wild-type CCP and mutant ccpΔcas1 did not differ in terms of mycelium growth, sporulation, and germination rate in vitro. Cas1 gene deletion induced a complete loss of virulence on the susceptible clones PB260 and IRCA631, as revealed by inoculation experiments on intact (non-detached) leaves. However, residual symptoms persisted when inoculations were conducted on detached leaves, notably with longer incubation times. Complementation with exogenous cassiicolin restored the mutant capacity to colonize the leaf tissues. We also compared the toxicity of CCP and ccpΔcas1 culture filtrates, through electrolyte leakage measurements on abraded detached leaves, over a range of clones as well as an F1 population derived from the cross between the clones PB260 (susceptible) and RRIM600 (tolerant). On average, filtrate toxicity was lower but not fully suppressed in ccpΔcas1 compared to CCP, with clone-dependent variations. The two QTL, previously found associated with sensitivity to CPP filtrate or to the purified cassiicolin, were no longer detected with the mutant filtrate, while new QTL were revealed. Our results demonstrate that: (1) cassiicolin is a necrotrophic effector conferring virulence to the CCP isolate in susceptible rubber clones and (2) other effectors produced by CCP contribute to residual filtrate toxicity and virulence in senescing/wounded tissues. These other effectors may be involved in saprotrophy rather than necrotrophy.

Genetic Determinism of Sensitivity to Corynespora cassiicola Exudates in Rubber Tree (Hevea brasiliensis).

An indirect phenotyping method was developed in order to estimate the susceptibility of rubber tree clonal varieties to Corynespora Leaf Fall (CLF) disease caused by the ascomycete Corynespora cassiicola. This method consists in quantifying the impact of fungal exudates on detached leaves by measuring the induced electrolyte leakage (EL%). The tested exudates were either crude culture filtrates from diverse C. cassiicola isolates or the purified cassiicolin (Cas1), a small secreted effector protein produced by the aggressive isolate CCP. The test was found to be quantitative, with the EL% response proportional to toxin concentration. For eight clones tested with two aggressive isolates, the EL% response to the filtrates positively correlated to the response induced by conidial inoculation. The toxicity test applied to 18 clones using 13 toxinic treatments evidenced an important variability among clones and treatments, with a significant additional clone x treatment interaction effect. A genetic linkage map was built using 306 microsatellite markers, from the F1 population of the PB260 x RRIM600 family. Phenotyping of the population for sensitivity to the purified Cas1 effector and to culture filtrates from seven C. cassiicola isolates revealed a polygenic determinism, with six QTL detected on five chromosomes and percentages of explained phenotypic variance varying from 11 to 17%. Two common QTL were identified for the CCP filtrate and the purified cassiicolin, suggesting that Cas1 may be the main effector of CCP filtrate toxicity. The CCP filtrate clearly contrasted with all other filtrates. The toxicity test based on Electrolyte Leakage Measurement offers the opportunity to assess the sensitivity of rubber genotypes to C. cassiicola exudates or purified effectors for genetic investigations and early selection, without risk of spreading the fungus in plantations. However, the power of this test for predicting field susceptibility of rubber clones to CLF will have to be further investigated.

Characterization of branching in two Hevea brasiliensis clones.

Rubber has been grown for several decades in Africa, notably in Ivory Coast. Although the yields obtained with selected clones are good, some problems with adaptation to edaphic and climatic conditions have yet to be solved. Of these problems, tree breakage due to violent winds during frequent storms is a major handicap. Some clones are more resistant to wind damage. However, this resistance trait remains difficult to measure, and detecting it at an early stage appeared important. As no differences have been found in the physical properties of the wood of different clones, the search for differences between clones displaying resistance and susceptibility to wind damage was switched to architectural structures. Architectural traits were thus studied on two rubber clones with very different resistances to wind. Observations were focused both on young trees and adult trees. Of all the architectural traits observed, a description of branching provided information likely to explain differences in the way the clones being studied developed over time. The clone GT1 which is resistant to wind, developed more lateral branches but they did not grow much in length. On the other hand, fewer axillary branches appeared on the susceptible clone PB235 but they grew longer. Moreover, as some of these traits proved to be similar in both young and adult trees, they could be detected at a sufficiently early stage for use as early prediction traits, if their involvement in susceptibility to wind damage is confirmed.