Global changes in gene expression during compatible and incompatible interactions of cowpea (Vigna unguiculata L.) with the root parasitic angiosperm Striga gesnerioides
BMC Genomics volume 13, Article number: 402 (2012)
Cowpea, Vigna unguiculata L. Walp., is one of the most important food and forage legumes in the semi-arid tropics. While most domesticated forms of cowpea are susceptible to the root parasitic weed Striga gesnerioides, several cultivars have been identified that show race-specific resistance. Cowpea cultivar B301 contains the RSG3-301 gene for resistance to S. gesnerioides race SG3, but is susceptible to race SG4z. When challenged by SG3, roots of cultivar B301 develop a strong resistance response characterized by a hypersensitive reaction and cell death at the site of parasite attachment. In contrast, no visible response occurs in B301 roots parasitized by SG4z.
Gene expression in the roots of the cowpea cultivar B301 during compatible (susceptible) and incompatible (resistant) interactions with S. gesnerioides races SG4z and SG3, respectively, were investigated at the early (6 days post-inoculation (dpi)) and late (13 dpi) stages of the resistance response using a Nimblegen custom design cowpea microarray. A total of 111 genes were differentially expressed in B301 roots at 6 dpi; this number increased to 2102 genes at 13 dpi. At 13 dpi, a total of 1944 genes were differentially expressed during compatible (susceptible) interactions of B301 with SG4z. Genes and pathways involved in signal transduction, programmed cell death and apoptosis, and defense response to biotic and abiotic stress were differentially expressed in the early resistance response; at the later time point, enrichment was primarily for defense-related gene expression, and genes encoding components of lignifications and secondary wall formation. In compatible interactions (B301 – SG4z), multiple defense pathways were repressed, including those involved in lignin biosynthesis and secondary cell wall modifications, while cellular transport processes for nitrogen and sulfur were increased.
Distinct changes in global gene expression profiles occur in host roots following successful and unsuccessful attempted parasitism by Striga. Induction of specific defense related genes and pathways defines components of a unique resistance mechanism. Some genes and pathways up-regulated in the host resistance response to SG3 are repressed in the susceptible interactions, suggesting that the parasite is targeting specific components of the host’s defense. These results add to our understanding of plant-parasite interactions and the evolution of resistance to parasitic weeds.
Cowpea (Vigna unguiculata L. Walp.) is the most important grain legume grown in sub-Saharan Africa [1, 2]. Approximately 12.5 million tons of cowpea grains are produced worldwide each year, with the majority of the production (over 64%) taking place on low-input, subsistence farms in west and central Africa . In these regions, cowpea is often referred to as “poor man’s meat” because of its high protein content (20-25%) and good nutritional value . The fruits are consumed at all stages of growth (e.g., green pods, fresh or dry seeds) and the young leaves are often used for soups and stews . In addition to its value as human food, cowpea hay is an important source of animal fodder . Two characteristics add to the plant’s agronomic importance: it is generally drought tolerant and it fixes nitrogen symbiotically, thereby enhancing soil fertility, especially when used in rotation with cereals .
Like most crops, cowpea growth and grain yields are greatly reduced by a variety of biotic pests (e.g., bacterial, fungal and viral diseases; insects; nematodes; herbivores) and biotic stresses (severe drought, salinity, and heat) . Among the major biotic constraints is parasitism by Striga gesnerioides L Walp. (Orobanchaceae), commonly referred to as witch weed. Witch weed is a noxious and persistent pest in farm fields: yield losses due to S. gesnerioides parasitism are extensive in the Sudano-Sahelian belt of west and central Africa . Control of the parasite is difficult because it produces thousands of seeds per generation that remain in the seed bank for years and also because most of the damage to its host plant occurs prior to its emergence from the ground . The damaging effects of Striga in this region are further compounded by poor soils and drought.
While most cowpea plants are susceptible to Striga parasitism, some local landraces and wild accessions have been identified that are resistant and in most reports, resistance is a dominant characteristic, inherited in a monogenic manner [2, 10]. Complicating the identification of Striga-resistant germplasm is the variable nature of the parasite with at least seven distinct races of S. gesnerioides [designated SG1 (Burkina Faso), SG2 (Mali), SG3 (Nigeria and Niger), SG4 (Benin), SG4z (Zakpota region of Benin), SG5 (Cameroon), and SG6 (Senegal)] now identified throughout West Africa [11–13]. Analysis of several advanced populations -- segregating for resistance to one or more of the different races of S. gesnerioides has resulted in the genetic mapping of several race-specific resistance (R) genes within the cowpea genome and the development of molecular markers linked to these genes . Using a positional cloning approach, Li and Timko  isolated and characterized a gene (designated RSG3-301) capable of conferring resistance to S. gesnerioides race 3 (SG3). RSG3-301 encodes an R protein homolog containing a coiled-coil (CC) protein-protein interaction domain at the N-terminus, a nucleotide binding site (NBS), and a leucine-rich repeat domain at the C-terminus. Silencing of RSG3-301 expression in the resistant cultivar B301 leads to susceptibility to race SG3, but does not affect resistance to other races of the parasite, underscoring the specificity of the resistance response .
Resistant cowpea genotypes exhibit two different response mechanisms to Striga attack. When challenged by a known race, cultivars carrying the appropriate race-specific resistance gene exhibit a rapid and robust hypersensitive response (HR) typified by browning and necrosis at the site of parasite attachment, and subsequent rapid death of the parasite within 3-4 days [13, 15, 16] (Figure 1). In host plants lacking the appropriate resistance gene, the parasite rapidly penetrates the host root cortex, forms connections to the host vascular system, swells to form a tubercle, and expands its cotyledons leading to subsequent above ground growth and flowering.
To date, only limited information is available on transcription changes in hosts in response to Striga parasitism [17–19]. To address this gap in our knowledge and to better understand the mechanics of host resistance to parasitic weeds, we have investigated the nature of global gene expression changes that occur during compatible and incompatible S. gesnerioides -cowpea interactions, using a custom designed cowpea Nimblegen microarray. These studies have identified candidate genes and potential pathways involved in resistant and susceptible responses to attempted Striga parasitism which can be tested further to gain additional information on this unique plant-plant interaction.
Design of cowpea microarray platform
Whole genome microarrays provide a means to access gene expression involved in a particular cellular or developmental event on a global scale. Since a microarray platform was not commercially available for cowpea at the time of this study, we designed a 385,000 feature long-nucleotide custom microarray to interrogate 43,253 cowpea unigenes of known and unknown function, previously identified in a gene-space sequencing study [20, 21]. On the fabricated microarray, each predicted gene coding region is represented by six to eight 60-mer oligonucleotide probes.
To validate our microarray design, we examined the distribution of mean probe set intensity and probe set variance over all assayed samples ( Additional file 1) to identify problematic probe sets that failed to hybridize (extremely low mean and variance) or hybridized nonspecifically (extremely high mean and low variance). The results of this analysis showed that of 49,287 probe sets on the microarray, 20,224 probe sets (41%) exhibited intensity profiles consistent with target-specific RNA hybridization. High intra-probe set covariance between individual probe intensities indicate that 75% of probe sets were composed of functional probes whose variation in intensities robustly represented abundance of the intended target RNA species ( Additional file 2).
Investigating global changes in gene expression during compatible and incompatible cowpea-Striga interactions
To examine the response of cowpea roots to attempted parasitism by S. gesnerioides, we used cowpea cultivar B301 as our host plant. Cowpea cultivar B301 is a land-race from Botswana that was originally identified as resistant to Striga in field studies and subsequently shown to be resistant to all races of S. gesnerioides known at the time . However, in the early 1990s, a new hyper virulent race of the parasite (SG4z), identified in the Zakpota region of Benin, was discovered to overcome the resistance present in B301 . B301 still maintains its resistance to all other races outside of this area. The differential response of B301 parasite races provides a good model for examining global changes in gene expression in the same genetic background during compatible/susceptible (B301-SG4z) and incompatible (B301-SG3) interactions (Figure 1). For our analysis, samples from three biological replicates were collected at two time points (6 dpi and 13 dpi). These time points were chosen on the basis of previous histological studies indicating that the earliest evidence of a hypersensitive response can be found at 6 dpi, with complete host cortical root browning and necrosis of the parasite evident by 13 dpi . Conversely, by 13 dpi of compatible interaction, host-parasite vascular connections have been established, and tubercle swelling and parasite cotyledon expansion have initiated.
Differential expression analysis was conducted as described in the Materials and Methods section, and genes were identified as differentially expressed if the adjusted P value was less than or equal to 0.05 (i.e. a 5% false discovery rate [FDR] threshold). Based on these criteria, at 6 dpi, the expression of 111 genes was significantly altered in B301 roots infected with S. gesnerioides race SG3 (resistant interaction) compared to uninfected control roots (Table 1). At 13 dpi, the number of significantly altered genes increases dramatically to 2102, among which 52 are genes whose expression was also altered at the earlier 6 dpi time point ( Additional file 3). A greater proportion (54%) of the differentially regulated genes were upregulated at both the early and late stages in the incompatible interactions. The large increase in the number of differentially expressed genes for the late time point of the resistant response is not unexpected since while activation of a defense cascade leading to the resistant phenotype and its associated hypersensitive response (i.e., browning, apoptosis) is rapid, propagation of the response and subsequent biosynthetic and physiological changes are broader and more persistent.
In compatible interactions of B301 with SG4z, 1944 genes had their expression altered compared to control roots at 13 dpi (Table 1). The expression of the majority (1089/1944) of these genes was down-regulated compared to uninfected control roots. Of the differentially regulated genes found, 923 genes were common in both resistant and susceptible interactions at 13 dpi (Additional file 3). There were 1179 and 1021 genes that were unique in the resistant and susceptible interactions, respectively.
The fold-change ratios, adjusted P values, and most significant GenBank annotations for selected genes from the various treatments are presented in Tables 2, 3, and 4. The genes were chosen based on the assumption that genes with the highest fold-change ratio are most likely to have a direct role in the plant-parasite interaction.
In the incompatible interaction between B301 and SG3 (Table 2), among the most highly up-regulated genes at 6 dpi are those encoding various transcription factors and signal transduction components, defense-related proteins (e.g., chitinase and peroxidase homologs, narbonin), and proteins involved in cell wall biogenesis. Among the most highly down-regulated genes at 6 dpi are those encoding cell cycle regulators and cellular transporters. At the later stages (13 dpi) of incompatible interactions (Table 3), up-regulated gene expression is dominated by genes encoding various transcription factors, signal transduction components, proteins involved in cellular energy metabolism and proteins associated with developmental regulation.
In contrast, in B301 roots being successfully parasitized by SG4z (Table 4), among the most highly down-regulated genes are those encoding defense-related proteins (e.g., peroxidase, chitinases), plant growth regulators (including multiple auxin transporters and auxin responsive regulators), and proteins involved in cell wall biogenesis (e.g., expansins, cellulose synthase, phenylalanine ammonium lyase and glucosyltransferase). Of the various genes up-regulated in this compatible interaction, of particular note are those encoding various nitrogen and sulfur transporters and permeases.
In order to verify that the changes in gene expression detected in our microarray analysis reflect actual changes in transcript levels, quantitative reverse-transcriptase PCR (qRT-PCR) was carried out using total RNAs samples isolated from plant tissue samples from identical treatments harvested independently. In most cases for the various genes tested by qRT-PCR, the levels of induction or repression (fold-change in expression) observed were comparable (Additional file 4).
Co-regulated genes and pathways in incompatible cowpea B301-S. gesnerioides race SG3 interactions
To identify genes and pathways that are co-regulated within the various treatments, differentially expressed genes were clustered based on their intensity values profile across treatments and expressed as a colored grid (heatmap), with genes ordered by a dendrogram representing covariance across treatments (Figure 2, Additional file 5). Several clusters of co-regulated genes can be identified within replicates of early (6 dpi) and late (13 dpi) incompatible interactions, and within compatible B301-SG4z interactions. Representative genes from these various clusters on the heatmap that show a full expression profile at the 0.1% FDR threshold are depicted in Additional file 6.
We used GOrilla (Gene Ontology enRIchment anaLysis and visuaLizAtion tool) to identify pathways enriched (P-values less than 10-3) with differentially expressed genes in the early (6 dpi) incompatible interaction between B301 and SG3 (Figure 3, Additional file 7). Among the biological processes showing the most significant changes in gene expression were pathways associated with cell death, programmed cell death, and apoptosis, as well as pathways associated with response to stress and biotic stress. We also see significant enrichment of differentially expressed genes associated with cell wall biogenesis and organization, and genes categorized under molecular processes, such as hydrolase activity and catalytic activity. Enrichment in these categories is not unexpected. One of the key features of incompatible interactions between cowpea and S. gesnerioides is the elicitation of a hypersensitive response when host roots expressing the appropriate race-specific resistance gene are challenged by a parasite with a recognized avirulence factor, such as in the case of B301 parasitized by S. gesnerioides race SG3 .
During the late stages of the B301-SG3 incompatible interaction (13 dpi), defense pathways associated with response to a variety of biotic stresses are among the biological processes showing the most significant changes in gene expression (Figure 4, Additional file 8). Significant alteration to gene expression was also seen in pathways involved in response to oxidative stress, signaling pathways for secondary metabolism, and pathways involved in multiple plants’ biosynthetic and chemical detoxification processes.
Co-regulated genes and pathways in compatible cowpea B301-S. gesnerioides race SG4z interactions
In striking contrast to the response of B301 roots to SG3, B301 roots parasitized by SG4z show no phenotypic manifestation during parasite ingress through the cortex. Nonetheless, dramatic alterations in gene expression are elicited with the greatest enrichment of differential expression found in pathways associated with cellular differentiation and growth, cell signaling and metabolism, and defense signaling (Figure 5, Additional file 9). These observed changes in cellular function are consistent with the dramatic down-regulation of gene expression involved in auxin transport and signaling, both critical for cellular growth and proliferation, and down-regulation of the expression of genes involved in cell wall growth (e.g., expansins) and reinforcement (e.g., enzymes of cellulose, lignin, and callus formation).
As it might be expected, B301 roots that were actively parasitized by SG4z also had gene expression differences enriched in pathways associated with defense signaling and response to biotic and abiotic stress. The active suppression of defense response cascades by pest- and pathogen-derived effectors that attack host cellular targets is well characterized; the global alterations in signaling and stress response cascades observed in our studies indicate that similar activities are likely in play in these host-parasite associations.
Gene expression in incompatible Striga-host interactions
The mechanics of plant host response to attempted parasitism and the nature of the gene expression changes occurring during both compatible and incompatible interactions between S. gesnerioides and cowpea have been investigated using oligonucleotide-based microarrays. Up until now, only limited information was available on global changes in transcription taking place in hosts in response to Striga parasitism. Our data show that alterations, both positive and negative, in a wide range of processes occur as a consequence of Striga attack. In cultivars harboring resistance genes capable of interfacing with a defense signal transduction pathway, as in the case of cowpea cultivar B301 which contains multirace resistance to S. gesnerioides, a rapid and robust defense response is elicited in the form of hypersensitivity at the site of parasite attachment, under conditions where the invading parasite is recognized. We found that accompanying the visible hypersensitive response is the up-regulation of genes involved in signal transduction and biosynthetic processes associated with formation barriers to prevent parasite ingress. These include activation of gene expression associated with cell wall biogenesis and strengthening (e.g., lignifications), as well as processes leading to programmed cell death at the host-parasite interface.
In a study of the interaction of S. hermonthica with its host plant sorghum, Hiraoka and Sugimoto  observed that attempted parasitism induced JA-responsive genes and suppressed SA-responsive genes in the roots of a highly susceptible cultivar, and suggested that susceptible hosts recognize Striga parasitism as wounding stress rather than microbial stress. In contrast, they found that in the roots of a moderately resistant sorghum cultivar, both SA- and JA-responsive genes were induced suggesting that resistance involved pathways associated with both wounding and pathogen challenge. Pretreatment of the roots of the moderately resistant sorghum cultivar with SA led to an inhibition of Striga development suggesting that the SA-responsive genes are directly involved in host resistance mechanism. In a subsequent study, subtractive hybridization was used to examine the interaction of Lotus japonicus with P. aegyptiaca (a compatible interaction) and S. hermonthica (incompatible) . These investigations found little or no overlap among the Phelipanche and Striga-induced transcripts, suggesting that L. japonicus roots are able to distinguish the compatible parasite from the incompatible one. Among the genes specifically induced by P. aegyptiaca were those encoding components of JA biosynthesis, whereas S. hermonthica parasitism induced genes in phytoalexin biosynthesis.
Examination of the differential expression of genes encoding candidate disease resistance and signal transduction components during the interactions of S. gesnerioides with various cowpea cultivars revealed that PR-5 transcript levels were dramatically elevated in the roots of cowpea genotypes resistant to Striga compared to uninfected roots and roots challenged with a race of Striga to which it was susceptible or adapted to another host species . In contrast, transcript levels of COI1 and EDS1 increased during susceptible and non-host interactions but were unchanged during resistance response. The COI1 gene product plays a pivotal role in activation of JA-mediated response cascades and, in some cases, serves as an inhibitor of SA signaling. EDS1 encodes lipase-like proteins that control defense activation and programmed cell death in plants. Induction of COI1 gene expression in the compatible Striga-cowpea interaction suggests that COI1 may down-regulate or suppress the resistance response and block SA signaling pathway in cowpea plants, thus allowing Striga attachment and further development.
Swarbrick et al.  characterized gene expression changes in the roots of the S. hermonthica resistant and susceptible rice cultivars Nipponbare and IAC 165 using a rice microarray. These investigators found that the levels of a large number of transcripts in rice roots were either positively (up-regulated) or negatively (down-regulated) affected by parasitism. Among the genes up-regulated in the Striga-resistant cultivar Nipponbare were ones encoding HR protein homologs; PR-proteins associated with microbial pathogenesis including endochitinases (PR-3), glucanases (PR-2), and thaumatin-like proteins (PR-5); pleiotropic drug resistance ABC transporters; and enzymes in phenylpropanoid metabolism. In addition, transcripts encoding several WRKY transcription factors (TFs), previously implicated in other SA-dependent resistance responses , were observed to be more abundant in parasitized roots. Large-scale down-regulation of gene expression was observed in the susceptible cultivar IAC 165, particularly transcripts whose encoded products annotate to Gene Ontology (GO) functional categories of plant growth regulator signaling and metabolism, biogenesis of cellular components, and cell division. Interestingly, a majority of the genes down-regulated following attempted Striga parasitism in both IAC165 and Nipponbare roots encodes products that annotate as proteins of unknown function.
A recent study that used whole genome oligonucleotide arrays to examine the transcriptome of A. thaliana roots after inoculation with pre-germinated S. hermonthica seeds found that large numbers of genes involved in cell wall synthesis, defense signaling, regulation of transcription and protein synthesis, oxidative stress, and primary and secondary metabolism were up-regulated at the earliest stages of parasite infection . This included up-regulation of many genes (EDS1, EDS5, PAD3, NPR1, NIMIN1, PR2) involved in the SA signaling pathway, as well as the up-regulation of a key WRKY transcription factor (AtWRKY70) that regulates the expression of genes involved in the SA signaling pathway and is thought to have a role in determining the balance between SA and JA signaling. In addition, there was evidence for the activation of genes involved in the JA and ethylene biosynthetic pathways.
In the present study of cowpea-S. gesnerioides interactions, among the most highly induced genes in the early resistance response (6 dpi) of B301 to SG3 are those involved in response to biotic stimuli, abiotic stimuli, wounding, oxidative stress and HR, and components of the JA and ETH signaling pathways. Among the genes most highly induced were a number of chitinases and chitinase homologs including narbonin, a protein with glycosyl hydrolase activity  whose biological function has yet to be elucidated . Narbonin has been previously suggested as interfering with parasite growth in the root by altering cell wall extensibility .
Lignifications and callused deposition have been previously observed at the interface of host plant cells with invading Orobanche haustorium  and is thought to be correlated with levels of host resistance to Orobanche. We similarly observed an increased expression of genes involved in cell wall biogenesis suggesting that lignifications may be part the mechanism of defense employed by cowpea against S. gesnerioides ingress.
A number of cytochrome P450s, previously implicated in various cellular detoxification pathways in plants, were also found to be induced during the resistance response of B301. Increased cytochrome P450 gene expression has been previously reported during the resistance reaction of rice infected by S. hermonthica.
The only other global study of cowpea transcriptional profiling published thus far is the work by Das et al.  who investigated differences in gene expression in resistant and susceptible cowpea genotypes elicited by feeding of the root-nematode, Meloidogyne incognita, using a soybean Affymetrix Gene Chip expression array. At 3 dpi in both compatible and incompatible interactions, more genes were down-regulated than up-regulated. When expression between infected resistant and susceptible genotypes was compared, a greater number and proportion of genes were down-regulated in the resistant than in the susceptible genotype, whereas more genes were up-regulated in the susceptible than in the resistant genotype in response to nematode infection. Gene ontology based functional categorization revealed that the typical defense response was partially suppressed in resistant roots, allowing nematode juvenile development. Das et al.  suggests that suppression of genes in ROS formation and other defense related responses might be important negative resistance mechanism. Our findings contrast with this, since we see activation of a wide range of processes as a consequence of Striga attack. Direct comparisons of gene lists suggest that in both early stages (3 dpi) of incompatible and compatible interactions there are many induced genes common in both host-parasite interactions suggesting that these may be part of a basal defense response. Among these were genes involved in lignifications and cross-linking. This response leads to heavy lignifications of cell walls and creates a mechanical barrier for the pathogen. In the incompatible cowpea-nematode interaction, expansin, thought to be important for maintaining the specialized feeding structures in hosts, was highly down-regulated. This was not observed in the cowpea-Striga interaction but down-regulation of expansins were observed in the response of a resistant rice cultivar to attack by S. hermonthica.
Gene expression in compatible Striga-host interactions
We know that root parasitic angiosperms secrete pectin methylesterases, polygalacturonase, and other cell wall softening enzymes during attempted penetration of the host root cortex, and that such secretions may in fact assist in overcoming the cell wall reinforcement that occurs in the host root in response to parasite parasitism [30–33]. It is now evident that many different phytopathogens have also evolved specific effectors and virulence factors that are capable of entering the host cell and suppressing the host resistance machinery or bypassing surveillance [34–36]. Among the most highly down-regulated genes in B301 roots during a susceptible interaction with SG4z were genes in the phenylpropanoid and lignin biosynthesis pathways, primary and secondary cell wall biogenesis, and components of the SA and JA signal transduction pathways.
S. gesnerioides-infected cowpea often demonstrated severe reduction in the growth rate and reduction of the host biomass . This was due to two major factors: direct carbon transfer from the host to S. gesnerioides, and reduction of photosynthetic capacity . One of the noticeable changes in gene expression during the susceptible response in B301 was the suppression of plant growth regulators, particularly auxin and gibberellins. Some genes, such as cellulose synthase, known to be involved in cell wall development  were also observed to be down-regulated.
One of the keys to the success of a parasitic plant-host interaction is the translocation of water and nutrients (including nitrates, amino acids, carbohydrates and minerals) from host to parasite through the haustoria . The identification of differentially up-regulated transcripts for transporters responsible for nitrogen, sulfur and amino acids in B301 roots during susceptible interactions with invading SG4z suggests that, in addition to suppressing some host functions to facilitate entry, the parasite is also modifying others to provide a source of nutrition. This is clearly an area that warrants additional study.
Distinct changes in global gene expression profiles occur in host roots following successful and unsuccessful attempted parasitism by Striga. Induction of specific defense-related genes and pathways defines components of a unique resistance mechanism. Some genes and pathways up-regulated in the host resistance response to SG3 are repressed in the susceptible interactions, suggesting that the parasite is targeting specific components of host defense. These results add to our understanding of plant-parasite interactions and the evolution of resistance to parasitic weeds.
Plant growth and infection with Striga gesnerioides
Cowpea cultivar B301 seeds were surface-sterilized with 1% hypochlorite for 5 min. They were then placed between two sheets of moist glass fiber filter paper (GF/A Whatman, Piscataway, NJ), held between two blocks of moistened rockwool (Grodan Inc., Milton, ON) for 5 days and then transferred to a growth chamber, which consists of a 24 cm x 24 cm x 3 cm Petri dish containing rockwool with a 100 μm mesh separating the cowpea roots from the rockwool . The chamber was covered with aluminum foil to keep light away from the roots. Cowpea seedlings were grown for another 7 days in a controlled environment growth room under a 12 H light–dark photoperiod at 30oC. Seeds of S. gesnerioides races SG3 and SG4z were surface sterilized and pre-conditioned for 9 days as previously described . The germination of Striga seeds was triggered using root exudates from cowpea cultivar B301. The roots of each cowpea seedling were inoculated with 15 mg of pre-germinated S. gesnerioides using a paintbrush. Root material was collected 6 and 13 days post-inoculation (dpi). Unattached Striga seeds and seedlings were removed by rinsing the roots with distilled water and the infected areas on the roots were collected by cutting at 0.5 cm at either site of the infection site. Sections (1 cm long) of control roots were collected from similar positions along the root where infections occurred in Striga treated samples. Three biological replicates were generated for each treatment with each replicate consisting of three independently inoculated cowpea plants. Tissue samples were frozen using liquid nitrogen and stored in a -80oC freezer.
RNA isolation, quantification, and hybridization to the cowpea microarray
Total RNA was isolated from frozen, finely ground root tissues as described by  with minor modification. RNA samples were treated with DNase (Roche Applied Science, Indianapolis, IN) and repurified using Qiagen RNeasy kit (Qiagen, Valencia, CA) according to the manufacturers’ instructions. RNA concentration and quality was determined using an Experion (Bio-Rad, Hercules, CA).
A 385,000 feature microarray was fabricated (Nimblegen Inc., Madison, WI) using 60-nucleotide long oligonucleotide synthesized based upon the nucleotide sequence of 43,253 cowpea unigenes of known and unknown function previously identified in a gene-space sequencing study [20, 21]. On the fabricated microarray, each predicted gene coding region is represented by six to eight 60-mer oligonucleotide probes. For hybridization, cDNA was prepared according to the protocol provided by Roche (Roche Applied Science, Indianapolis, IN). cDNA quality control and hybridization to the 385 K arrays were performed at Nimblegen.
Analysis of microarray data
Statistical analysis of Nimblegen microarray data was performed using software from the Bioconductor project  and implemented in the R environment for statistical computing and graphics (R Development Core Team, 2009), version 2.10.1.
Prior to reading the data from the raw expression values from the XYS files, an annotation package was build for the Nimblegen custom design file “2007-05-01_UVir_CowPea_expr.ndf” using the pdInfoBuilder package protocol  version 1.10.1. The array data was read into R using the oligo package protocol  version 1.10.4. RMA background subtraction and normalization were applied to the microarray set. The experiment was performed using three biological replicates.
Differential expression analysis was conducted using the Limma package, version 3.2.3. Limma analysis applies linear models to each gene and an empirical Bayesian method to moderate standard errors for estimation of fold changes. A multiple testing correction was performed using the using Benjamin-Hochberg error correction model. Genes were identified as differentially expressed if the adjusted P-value, false discovery rate (FDR) was less than 0.1.
The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus  and are accessible through GEO Series accession number (GSE39348, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE39348).
Assembled GSR sequences were mapped against their Arabidopsis protein sequences from TAIR10  using FASTX  version 36.06 with a P-value threshold of 10E-3. GSR sequences ordered by FDR values were analyzed for gene enrichment using GOrilla  with a P-value threshold of 10E-3. Assembled GSR sequences were annotated against their best hit against plant refseq  using FASTX with a P-value threshold of 10E-3.
Annotated and assembled GSRs were used to generate a heatmap of normalized intensity values at 0.1% FDR threshold using enhanced heatmap function in R. GSRs with significant P-values of 10E-3 have been grouped based on their intensity profile and were used to generate GO Trees using Amigo . Enriched GO terms were included on the side and were matched to corresponding genes on the heatmap.
Quantitative reverse transcriptase PCR verification of transcript levels
Transcript levels in various tissue samples were validated by qRT-PCR using new total RNA samples prepared from B301 root tissues before and following attempted Striga parasitism as described above. cDNAs were synthesized using the Invitrogen Thermo script Kit (Invitrogen, Caarlsbad, CA). qRT-PCR reactions were carried out using the iQ SYBR Green Super mix Kit (Bio-Rad, Hercules, CA) and gene-specific primers (Additional file 10) against selected differentially expressed target genes. Reactions were carried out on an cycler Optical Module PCR according to the manufacturer’s instructions and gene expression fold changes were calculated according to Schmittgen and Livak . Expression was normalized using an actin reference gene (XM_003521168; sequence id: 33645964, GI: 356504867) giving a 132 bp amplicon. Previous studies  have shown that expression of this actin reference gene is constant throughout development and during interactions with S. gesnerioides.
Ehlers JD, Hall AE: Cowpea (Vigna unguiculata L Walp). Field Crop Res. 1997, 53: 187-204. 10.1016/S0378-4290(97)00031-2.
Timko MP, Singh BB: Cowpea, a multifunctional legume. Genomics of Tropical Crop Plants. Edited by: Moore PH, Ming R. 2008, New York: Springer Science + Business Media, LLC, 227-258.
Langyintuo AS, Lowenberg-DeBoer J, Faye M, Lambert D, Ibro G, Moussa B, Kergna A, Kushwaha S, Musa S, Ntoukam G: Cowpea supply and demand in West and Central Africa. Field Crop Res. 2003, 82: 215-231. 10.1016/S0378-4290(03)00039-X.
Diouf D, Hilu K: Microsatellites and RAPD markers to study genetic relationships among cowpea breeding lines and local varieties in Senegal. Genet Resour Crop Ev. 2005, 52: 1057-1067. 10.1007/s10722-004-6107-z.
Quaye W, Adofo K, Madode Y, Abdul-Razak A: Exploratory and multidisciplinary survey of the cowpea network in the Tolon-Kumbungu district of Ghana: A food sovereignty perspective. Afr J Agric Res. 2009, 4: 311-320.
Tarawali SA, Singh BB, Gupta SC, Tabo R, Harris F, Nokoe S, Fernández-Rivera S, Bationo A, Manyong VM, Makinde K, Odion EC: Cowpea as a key factor for a new approach to intergrated crop-livestock system research in the dry savannas of West Africa. Challenges and opportunities for enhancing sustainable cowpea production. Edited by: Fatokun CA, Tarawali SA, Singh BB, Kormawa PM, Tamò M. 2002, Ibadan, Nigeria: Proceedings of the World Cowpea Conference III held at the International Institute of Tropical Agriculture (IITA), 233-251. 4–8 September 2000
Sanginga N, Dashiell KE, Diels J, Vanlauwe B, Lyasse O, Carsky RJ, Tarawali S, Asafo-Adjei B, Menkir A, Schulz S, Singh BB, Chikoye D, Keatinge D, Ortiz R: Sustainable resource management coupled to resilient germplasm to provide new intensive cereal-grain-legume-livestock systems in the dry savanna. Agr Ecosyst Environ. 2003, 100: 305-314. 10.1016/S0167-8809(03)00188-9.
Parker C: Observations on the current status of Orobanche and Striga problems worldwide. Pest Manag Sci. 2009, 65: 453-459. 10.1002/ps.1713.
Ejeta G: Breeding for Striga resistance in sorghum: Exploitation of an intricate host-parasite biology. Crop Sci. 2007, 47: S216-S227.
Timko M, Ehlers JD, Roberts PA: Cowpea. Genome Mapping and Molecular Breeding in Plants. Volume 3 Pulses, Sugar and Tuber Crops. Edited by: Kole C. 2007, Berlin Heidelberg: Springer Verlag, 49-67.
Lane JA, Moore THM, Child DV, Cardwell KF, Singh BB, Bailey JA: Virulence characteristics of a new race of the parasitic angiosperm, Striga gesnerioides, from southern Benin on cowpea (Vigna-unguiculata). Euphytica. 1993, 72: 183-188.
Lane JA: Resitance of cowpea [Vigna unguiculata (L.) Walp.] to Striga gesnerioides (Willd.) Vatke, a parasitic angiosperm. New Phytol. 1993, 125: 405-412. 10.1111/j.1469-8137.1993.tb03893.x.
Botanga CJ, Timko MP: Phenetic relationships among different races of Striga gesnerioides (Willd.) Vatke from West Africa. Genome. 2006, 49: 1351-1365. 10.1139/g06-086.
Li J, Timko MP: Gene-for-gene resistance in Striga-cowpea associations. Science. 2009, 325: 1094-10.1126/science.1174754.
Botanga CJ, Timko MP: Genetic structure and analysis of host and nonhost interactions of Striga gesnerioides (witchweed) from central Florida. Phytopathology. 2005, 95: 1166-1173. 10.1094/PHYTO-95-1166.
Bailey JA, Reiss GC: Striga gesnerioides parasitising cowpea: Development of infection structures and mechanisms of penetration. Annals of Bot. 1998, 81: 431-440. 10.1006/anbo.1997.0577.
Hiraoka Y, Sugimoto Y: Molecular responses of sorghum to purple witchweed (Striga hermonthica) parasitism. Weed Sci. 2008, 56: 356-363. 10.1614/WS-07-136.1.
Li J, Lis KE, Timko MP: Molecular genetics of race-specific resistance of cowpea to Striga gesnerioides (Willd.). Pest Manag Sci. 2009, 65: 520-527. 10.1002/ps.1722.
Swarbrick PJ: Global patterns of gene expression in rice cultivars undergoing a susceptible or resistant interaction with the parasitic plant Striga hermonthica. New Phytol. 2008, 179: 515-529. 10.1111/j.1469-8137.2008.02484.x.
Chen X, Laudeman TW, Rushton PJ, Spraggins TA, Timko MP: CGKB: an annotation knowledge base for cowpea (Vigna unguiculata L.) methylation filtered genomic genespace sequences. BMC Bioinformatics. 2007, 8: 129-10.1186/1471-2105-8-129.
Timko MP, Rushton PJ, Laudeman TW, Bokowiec MT, Chipumuro E, Cheung F, Town CD, Chen X: Sequencing and analysis of the gene-rich space of cowpea. BMC Genomics. 2008, 9: 103-10.1186/1471-2164-9-103.
Singh BB, Emechebe AM: Inheritance of Striga resistance in cowpea genotype B301. Crop Sci. 1990, 30: 879-881. 10.2135/cropsci1990.0011183X003000040023x.
Hiraoka Y, Ueda H, Sugimoto Y: Molecular responses of Lotus japonicus to parasitism by the compatible species Orobanche aegyptiaca and the incompatible species Striga hermonthica. J Exp Bot. 2009, 60: 641-650. 10.1093/jxb/ern316.
Ryu HS, Han M, Lee SK, Cho JI, Ryoo N, Heu S, Lee YH, Bhoo S, Wang GL, Hahn TR, Jeon JS: A comprehensive expression analysis of the WRKY gene superfamily in rice plants during defense response. Plant Cell Rep. 2006, 25: 836-847. 10.1007/s00299-006-0138-1.
Vasey R: The identification of novel sources of resistance to Striga species. 2005, Sheffield: PhD thesis. Sheffield University
Coulson AF: A proposed structure for family-18 chitinases - a possible function for narbonin. FEBS Lett. 1994, 354: 41-44. 10.1016/0014-5793(94)01084-6.
Steffens P, Nong VH, Hillmer S, Saalbach I, Muntz K: Subcellular localization of the 2S globulin narbonin in seeds of Vicia narbonensis. Planta. 1997, 203: 44-50.
Letousey P, de Zelicourt A, Dos Santos CV, Thoiron S, Monteau F, Simier P, Thalouarn P, Delavault P: Molecular analysis of resistance mechanisms to Orobanche cumana in sunflower. Plant Pathol. 2007, 56: 536-546. 10.1111/j.1365-3059.2007.01575.x.
Das S, Ehlers JD, Close TJ, Roberts PA: Transcriptional profiling of root nematode induced feeding sites in cowpea (Vigna unguiculata L. Walp) using a soybean genome array. BMC Genomics. 2010, 11: 480-10.1186/1471-2164-11-480.
Borevitz JO, Xia YJ, Blount J, Dixon RA, Lamb C: Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell. 2000, 12: 2383-2393.
Scholes JD, Press MC: Striga infestation of cereal crops - an unsolved problem in resource limited agriculture. Curr Opin Plant Biol. 2008, 11: 180-186. 10.1016/j.pbi.2008.02.004.
Irving LJ, Cameron DD: You are what you eat: interactions between root parasitic plants and their hosts. Adv Bot Res. 2009, 50: 87-138.
Scholes JD, Timko MP: Host reaction to attack by root parasitic plant. Parasitic Orobanchaceae - Parasitic mechanisms and control strategies. Edited by: Joel DM, Musselman LJ, Gressel J. 2013, Springer Verlag, In Press
Boller T, He SY: Innate immunity in plants: an arms race between pattern recognition receptors in plants and effectors in microbial pathogens. Science. 2009, 324: 742-744. 10.1126/science.1171647.
Torto-Alalibo T, Collmer CW, Lindeberg M, Bird D, Collmer A, Tyler BM: Common and contrasting themes in host cell-targeted effectors from bacterial, fungal, oomycete and nematode plant symbionts described using the Gene Ontology. BMC Microbiol. 2009, 9:
Torto-Alalibo T, Meng S, Dean RA: Infection strategies of filamentous microbes described with the Gene Ontology. Trends Microbiol. 2009, 17: 320-327. 10.1016/j.tim.2009.05.003.
Parker C, Riches CR: Parasitic weeds of the world: biology and control. 1993, Wallingford: CAB International
Kalluri UC, Joshi CP: Differential expression patterns of two cellulose synthase genes are associated with primary and secondary cell wall development in aspen trees. Planta. 2004, 220: 47-55. 10.1007/s00425-004-1329-z.
Wang X, Tian W, Li Y: Development of an efficient protocol of RNA isolation from recalcitrant tree tissues. Mol Biotechnol. 2008, 38: 57-64. 10.1007/s12033-007-0073-6.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge YC, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang JH: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.
The Arabidopsis Database:http://www.ftp.arabidopsis.org/home/tair/Sequences/blast_datasets/TAIR10_blastsets/TAIR10_pep_20101028,
Pearson WR: Flexible sequence similarity searching with the FASTA3 program package. Methods Mol Biol. 2000, 132: 185-
Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009, 10: 48-10.1186/1471-2105-10-48.
The Plant Refseq:http://ftp.ncbi.nih.gov/refseq/release/plant/,
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C-T method. Nat Protoc. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.
Mellor KE, Hoffman AM, Timko MP: Use of ex vitro composite plants to study the interaction of cowpea (Vigna unguiculata L.) with the root parasitic angiosperm Striga gesnerioides. Plant Methods. 2012, 8: 22-10.1186/1746-4811-8-22.
We wish to thank Stephen Turner for his help in submission of the microarray data. We also thank Tom Laudemann, Jim Roberts, and Shuran Li for their help with various aspects of these studies. This work was supported in part by grants from the National Science Foundation (DBI-0701748 and IBN-0322420) and Kirkhouse Trust to MPT.
Authors declare that they have no competing interests.
MPT conceived of the project and was responsible for directing all of the research activities. KEM and KH carried out the experimental aspects of the work and assisted SNP, MJL, and AJM in computational analysis. SNP, MJL, and AJM carried out the statistical analysis of the microarray data, and worked with KEM and KH in refining gene annotations and expression analysis. All authors have assisted in the writing of the manuscript and have read and approved the final submitted version of the manuscript.
Kan Huang, Karolina E Mellor contributed equally to this work.
Electronic supplementary material
Additional file 3: Venn diagrams for genes up-regulated and down-regulated during resistant and susceptible interaction of cowpea with A) Striga race SG3 and B) Striga race SG4z.(PDF 251 KB)
Additional file 6: Representative genes from the heatmap showing a full expression profile at 0.1% FDR threshold.(PDF 218 KB)
About this article
Cite this article
Huang, K., Mellor, K.E., Paul, S.N. et al. Global changes in gene expression during compatible and incompatible interactions of cowpea (Vigna unguiculata L.) with the root parasitic angiosperm Striga gesnerioides. BMC Genomics 13, 402 (2012). https://doi.org/10.1186/1471-2164-13-402
- Striga gesnerioides
- Vigna unguiculata