Integrative transcriptomics reveals association of abscisic acid and lignin pathways with cassava whitefly resistance

BACKGROUND

Whiteflies are a global threat to crop yields, including the African subsistence crop cassava (Manihot esculenta). Outbreaks of superabundant whitefly populations throughout Eastern and Central Africa in recent years have dramatically increased the pressures of whitefly feeding and virus transmission on cassava. Whitefly-transmitted viral diseases threaten the food security of hundreds of millions of African farmers, highlighting the need for developing and deploying whitefly-resistant cassava. However, plant resistance to whiteflies remains largely poorly characterized at the genetic and molecular levels. Knowledge of cassava-defense programs also remains incomplete, limiting characterization of whitefly-resistance mechanisms. To better understand the genetic basis of whitefly resistance in cassava, we define the defense hormone- and Aleurotrachelus socialis (whitefly)-responsive transcriptome of whitefly-susceptible (COL2246) and whitefly-resistant (ECU72) cassava using RNA-seq. For broader comparison, hormone-responsive transcriptomes of Arabidopsis thaliana were also generated.

RESULTS

Whitefly infestation, salicylic acid (SA), jasmonic acid (JA), ethylene (ET), and abscisic acid (ABA) transcriptome responses of ECU72 and COL2246 were defined and analyzed. Strikingly, SA responses were largely reciprocal between the two cassava genotypes and we suggest candidate regulators. While susceptibility was associated with SA in COL2246, resistance to whitefly in ECU72 was associated with ABA, with SA-ABA antagonism observed. This was evidenced by expression of genes within the SA and ABA pathways and hormone levels during A. socialis infestation. Gene-enrichment analyses of whitefly- and hormone-responsive genes suggest the importance of fast-acting cell wall defenses (e.g., elicitor recognition, lignin biosynthesis) during early infestation stages in whitefly-resistant ECU72. A surge of ineffective immune and SA responses characterized the whitefly-susceptible COL2246's response to late-stage nymphs. Lastly, in comparison with the model plant Arabidopsis, cassava's hormone-responsive genes showed striking divergence in expression.

CONCLUSIONS

This study provides the first characterization of cassava's global transcriptome responses to whitefly infestation and defense hormone treatment. Our analyses of ECU72 and COL2246 uncovered possible whitefly resistance/susceptibility mechanisms in cassava. Comparative analysis of cassava and Arabidopsis demonstrated that defense programs in Arabidopsis may not always mirror those in crop species. More broadly, our hormone-responsive transcriptomes will also provide a baseline for the cassava community to better understand global responses to other yield-limiting pests/pathogens.

Genome-wide analyses of cassava Pathogenesis-related (PR) gene families reveal core transcriptome responses to whitefly infestation, salicylic acid and jasmonic acid

Background

Whiteflies are a threat to cassava (Manihot esculenta), an important staple food in many tropical/subtropical regions. Understanding the molecular mechanisms regulating cassava’s responses against this pest is crucial for developing control strategies. Pathogenesis-related (PR) protein families are an integral part of plant immunity. With the availability of whole genome sequences, the annotation and expression programs of the full complement of PR genes in an organism can now be achieved. An understanding of the responses of the entire complement of PR genes during biotic stress and to the defense hormones, salicylic acid (SA) and jasmonic acid (JA), is lacking. Here, we analyze the responses of cassava PR genes to whiteflies, SA, JA, and other biotic aggressors.

Results

The cassava genome possesses 14 of the 17 plant PR families, with a total of 447 PR genes. A cassava PR gene nomenclature is proposed. Phylogenetic relatedness of cassava PR proteins to each other and to homologs in poplar, rice and Arabidopsis identified cassava-specific PR gene family expansions. The temporal programs of PR gene expression in response to the whitefly (Aleurotrachelus socialis) in four whitefly-susceptible cassava genotypes showed that 167 of the 447 PR genes were regulated after whitefly infestation. While the timing of PR gene expression varied, over 37% of whitefly-regulated PR genes were downregulated in all four genotypes. Notably, whitefly-responsive PR genes were largely coordinately regulated by SA and JA. The analysis of cassava PR gene expression in response to five other biotic stresses revealed a strong positive correlation between whitefly and Xanthomonas axonopodis and Cassava Brown Streak Virus responses and negative correlations between whitefly and Cassava Mosaic Virus responses. Finally, certain associations between PR genes in cassava expansions and response to biotic stresses were observed among PR families.

Conclusions

This study represents the first genome-wide characterization of PR genes in cassava. PR gene responses to six biotic stresses and to SA and JA are demonstrably different to other angiosperms. We propose that our approach could be applied in other species to fully understand PR gene regulation by pathogens, pests and the canonical defense hormones SA and JA.

A metabolomics characterisation of natural variation in the resistance of cassava to whitefly

BACKGROUND

Cassava whitefly outbreaks were initially reported in East and Central Africa cassava (Manihot esculenta Crantz) growing regions in the 1990's and have now spread to other geographical locations, becoming a global pest severely affecting farmers and smallholder income. Whiteflies impact plant yield via feeding and vectoring cassava mosaic and brown streak viruses, making roots unsuitable for food or trading. Deployment of virus resistant varieties has had little impact on whitefly populations and therefore development of whitefly resistant varieties is also necessary as part of integrated pest management strategies. Suitable sources of whitefly resistance exist in germplasm collections that require further characterization to facilitate and assist breeding programs.

RESULTS

In the present work, a hierarchical metabolomics approach has been employed to investigate the underlying biochemical mechanisms associated with whitefly resistance by comparing two naturally occurring accessions of cassava, one susceptible and one resistant to whitefly. Quantitative differences between genotypes detected at pre-infestation stages were consistently observed at each time point throughout the course of the whitefly infestation. This prevalent differential feature suggests that inherent genotypic differences override the response induced by the presence of whitefly and that they are directly linked with the phenotype observed. The most significant quantitative changes relating to whitefly susceptibility were linked to the phenylpropanoid super-pathway and its linked sub-pathways: monolignol, flavonoid and lignan biosynthesis. These findings suggest that the lignification process in the susceptible variety is less active, as the susceptible accession deposits less lignin and accumulates monolignol intermediates and derivatives thereof, differences that are maintained during the time-course of the infestation.

CONCLUSIONS

Resistance mechanism associated to the cassava whitefly-resistant accession ECU72 is an antixenosis strategy based on reinforcement of cell walls. Both resistant and susceptible accessions respond differently to whitefly attack at biochemical level, but the inherent metabolic differences are directly linked to the resistance phenotype rather than an induced response in the plant.

Antagonistic role of 9-lipoxygenase-derived oxylipins and ethylene in the control of oxidative stress, lipid peroxidation and plant defence

9-lipoxygenases (9-LOXs) initiate fatty acid oxygenation in plant tissues, with formation of 9-hydroxy-10,12,15-octadecatrienoic acid (9-HOT) from linolenic acid. A lox1 lox5 mutant, which is deficient in 9-LOX activity, and two mutants noxy6 and noxy22 (non-responding to oxylipins), which are insensitive to 9-HOT, have been used to investigate 9-HOT signalling. Map-based cloning indicated that the noxy6 and noxy22 mutations are located at the CTR1 (CONSTITUTIVE ETHYLENE RESPONSE1) and ETO1 (ETHYLENE-OVERPRODUCER1) loci, respectively. In agreement, the noxy6 and noxy22 mutants, renamed as ctr1-15 and eto1-14, respectively, showed enhanced ethylene (ET) production. The correlation between increased ET production and reduced 9-HOT sensitivity indicated by these results was supported by experiments in which exogenously added ethylene precursor ACC (1-aminocyclopropane-1-carboxylic acid) impaired the responses to 9-HOT. Moreover, a reciprocal interaction between ET and 9-HOT signalling was indicated by results showing that the effect of ACC was reduced in the presence of 9-HOT. We found that the 9-LOX and ET pathways regulate the response to the lipid peroxidation-inducer singlet oxygen. Thus, the massive transcriptional changes seen in wild-type plants in response to singlet oxygen were greatly affected in the lox1 lox5 and eto1-14 mutants. Furthermore, these mutants displayed enhanced susceptibility to both singlet oxygen and Pseudomonas syringae pv. tomato, in the latter case leading to increased accumulation of the lipid peroxidation product malondialdehyde. These findings demonstrate an antagonistic relationship between products of the 9-LOX and ET pathways, and suggest a role for the 9-LOX pathway in modulating oxidative stress, lipid peroxidation and plant defence.

A new chromosome nomenclature system for oat (Avena sativa L. and A. byzantina C. Koch) based on FISH analysis of monosomic lines

Fluorescent in situ hybridization (FISH) with multiple probes was used to analyze mitotic and meiotic chromosome spreads of Avena sativa cv 'Sun II' monosomic lines, and of A. byzantina cv 'Kanota' monosomic lines from spontaneous haploids. The probes used were A. strigosa pAs120a (a repetitive sequence abundant in A-genome chromatin), A. murphyi pAm1 (a repetitive sequence abundant in C-genome chromatin), A. strigosa pITS (internal transcribed spacer of rDNA) and the wheat rDNA probes pTa71 (nucleolus organizer region or NOR) and pTa794 (5S). Simultaneous and sequential FISH employing pairs of these probes allowed the identification and genome assignation of all chromosomes. FISH mapping using mitotic and meiotic metaphases facilitated the genomic and chromosomal identification of the monosome in each line. Of the 17 'Sun II' lines analyzed, 13 distinct monosomic lines were found, corresponding to four monosomes of the A-genome, five of the C-genome and four of the D-genome. In addition, 12 distinct monosomic lines were detected among the 20 'Kanota' lines examined, corresponding to six monosomes of the A-genome, three of the C-genome and three of the D-genome. The results show that 19 chromosomes out of 21 of the complement are represented by monosomes between the two genetic backgrounds. The identity of the remaining chromosomes can be deduced either from one intergenomic translocation detected on both 'Sun II' and 'Kanota' lines, or from the single reciprocal, intergenomic translocation detected among the 'Sun II' lines. These results permit a new system to be proposed for numbering the 21 chromosome pairs of the hexaploid oat complement. Accordingly, the A-genome contains chromosomes 8A, 11A, 13A, 15A, 16A, 17A and 19A; the C-genome contains chromosomes 1C, 2C, 3C, 4C, 5C, 6C and 7C; and the D-genome consists of chromosomes 9D, 10D, 12D, 14D, 18D, 20D and 21D. Moreover, the FISH patterns of 16 chromosomes in 'Sun II' and 15 in 'Kanota' suggest that these chromosomes could be involved in intergenomic translocations. By comparing the identities of individually translocated chromosomes in the two hexaploid species with those of other hexaploids, we detected different types of intergenomic translocations.

Cloning and characterization of resistance gene analogs from Avena species

Sequences analogous to plant resistance genes of the NBS–LRR class were cloned from the genomic DNA of 11 Avena species with different genomes and levels of ploidy. Three pairs of degenerate primers were used, based on conserved DNA sequence motifs belonging to the NBS domain, and 33 sequences were identified. These were subdivided into 7 classes depending on nucleotide sequence identity. Despite the high level of degeneracy, the primers behaved in a highly selective way; the majority of sequences from the different species obtained with every primer combin ation showed strong identity and were considered homologous. For most species, only one sequence of each class was identified in each genome, suggesting that duplicated sequences are fairly divergent. The strong identity among specific NBS sequences precludes any conclusions being made on the evolution of these species. The genomic organization of the RGA sequences was explored using those of A. strigosa as probes in Southern blots involving digested DNA from 15 Avena species. The hybridization patterns showed wide diversity both among sequences within a species and among species for each sequence. However, the dendrogram generated using the RFLPs showed relationships among species to be in good agreement with those previously established using other molecular markers.Key words: resistance gene analog (RGA), disease resistance genes, diversity, Avena, oats.

Isolation and mapping of resistance gene analogs from the Avena strigosa genome

Degenerate primers based on conserved regions of the nucleotide binding site (NBS) domain (encoded by the largest group of cloned plant disease resistance genes) were used to isolate a set of 15 resistance gene analogs (RGA) from the diploid species Avena strigosa Schreb. These were grouped into seven classes on the basis of 60% or greater nucleic acid sequence identity. Representative clones were used for genetic mapping in diploid and hexaploid oats. Two RGAs were mapped at two loci of the linkage group AswBF belonging to the A. strigosa x A. wiestii Steud map, and ten RGAs were mapped at 15 loci in eight linkage groups belonging to the A. byzantina C. Koch cv. Kanota x A. sativa L. cv. Ogle map. A similar approach was used for targeting genes encoding receptor-like kinases. Three different sequences were obtained and mapped to two linkage groups of the hexaploid oat map. Associations were explored between already known disease resistance loci mapped in different populations and the RGAs. Molecular markers previously linked to crown rust and barley yellow dwarf resistance genes or quantitative trait loci were found in the Kanota x Ogle map linked to RGAs at a distance ranging from 0 cM to 20 cM. Homoeologous RGAs were found to be linked to loci either conferring resistance to different isolates of the same pathogen or to different pathogens. This suggests that these RGAs identify genome regions containing resistance gene clusters.

Fluorescence in situ hybridization mapping of Avena sativa L. cv. SunII and its monosomic lines using cloned repetitive DNA sequences

Fluorescent in situ hybridization (FISH) employing multiple probes was used with mitotic or meiotic chromosome spreads of Avena sativa L. cv. SunII and its monosomic lines to produce physical chromosome maps. The probes used were Avena strigosa pAs120a (which hybridizes exclusively to A-genome chromosomes), Avena murphyi pAm1 (which hybridizes exclusively to C-genome chromosomes), A. strigosa pAs121 (which hybridizes exclusively to A- and D-genome chromosomes), and the wheat rDNA probes pTa71 and pTa794. Simultaneous and sequential FISH employing two-by-two combinations of these probes allowed the unequivocal identification and genome assignation of all chromosomes. Ten pairs were found carrying intergenomic translocations: (i) between the A and C genomes (chromosome pair 5A); (ii) between the C and D genomes (pairs 1C, 2C, 4C, 10C, and 16C); and (iii) between the D and C genomes (pairs 9D, 11D, 13D, and 14D). The existence of a reciprocal intergenomic translocation (10C-14D) is also proposed. Comparing these results with those of other hexaploids, three intergenomic translocations (10C, 9D, and 14D) were found to be unique to A. sativa cv. SunII, supporting the view that 'SunII' is genetically distinct from other hexaploid Avena species and from cultivars of the A. sativa species. FISH mapping using meiotic and mitotic metaphases facilitated the genomic and chromosomal identification of the aneuploid chromosome in each monosomic line. Of the 18 analyzed, only 11 distinct monosomic lines were actually found, corresponding to 5 lines of the A genome, 2 lines of the C genome, and 4 lines of the D genome. The presence or absence of the 10C-14D interchange was also monitored in these lines.