Electronic Journal of Biotechnology
versión ISSN 0717-3458
Electron. J. Biotechnol. vol.15 no.5 Valparaíso set. 2012
Isolation and analysis of differentially expressed genes from peanut in response to challenge with Ralstonia solanacearum
Yu Fei Ding1 · Chuan Tang Wang*1 · Yue Yi Tang1 · Xiu Zhen Wang1 · Qi Wu1 · Dong Qing Hu2 · Hong Tao Yu1 · Jian Cheng Zhang1 · Feng Gao Cui1 · Guo Sheng Song1 · Hua Yuan Gao3 · Shan Lin Yu1
1Shandong Peanut Research Institute, Qingdao, China
*Corresponding author: firstname.lastname@example.org
Financial support: The authors are most grateful to the following funding sources, which made the research possible: China Agricultural Research System (CARS-14), Qingdao Science & Technology Support Program (10-3-3-20-nsh, 09-1-3-67-jch), Shandong Natural Science Foundation (Y2008D11), and Shandong Key Project of Science & Technology (2009GG10009008).
Keywords: Arachis hypogaea, ARF, cyclophilin, differentially expressed genes, genefishing, Ralstonia solanacearum.
Background: Bacterial wilt caused by Ralstonia solanacearum is the most devastating disease in peanut. Planting resistant peanut cultivars is deemed as the sole economically viable means for effective control of the disease. To understand the molecular mechanism underlying resistance and facilitate breeding process, differences in gene expression between seeds of Rihua 1 (a Virginia type peanut variety resistant to bacterial wilt) inoculated with the bacterial pathogen suspension (109 cfu ml-1) and seeds of the same cultivar treated with water (control), were studied using the GenefishingTM technology. Results: A total of 25 differentially expressed genes were isolated. Expression of genes encoding cyclophilin and ADP-ribosylation factor, respectively, were further studied by real time RT-PCR, and full length cDNAs of both genes were obtained by rapid amplification of cDNA ends. Conclusions: The study provided candidate genes potentially useful for breeding peanut cultivars with both high yield and bacterial wilt resistance, although confirmation of their functions through transgenic studies is still needed.
Being a rich source of cooking oil and dietary protein, the cultivated peanut (Arachis hypogaea L.) is regarded as one of the most important cash crops in the world. Bacterial wilt (BW) disease caused by Ralstonia solanacearum (Yabuuchi et al. 1995) poses a great threat to peanut production in China and Southeast Asia. More than 10% of the area under peanut is affected in China (Yu et al. 2011). Yield reduction generally ranges from 10% to 20%; however, in heavily infested field, over 50% yield losses are not uncommon. In extreme cases, the disease may even cause total crop failure (Yu et al. 2011).
Like any other plant BW diseases, peanut BW is difficult to control. No desirable chemical management measure is currently available. Planting resistant peanut cultivars is deemed as the sole economically viable means for effective control (Yu et al. 2011). Unfortunately, most of the resistant germplasm lines identified are small-seeded genotypes with low yield potential; transferring BW resistance to high yielding adapted peanut cultivars has therefore become an urgent task (Yu et al. 2011).
Understanding the mechanism underlying BW resistance at molecular level may hasten the breeding process. Thus far, in peanut, there have been several reports regarding identification of DNA markers related to BW resistance (Yu et al. 2011), and transcripts involved in response to biotic and abiotic stresses other than BW have been identified. These included, transcripts related to resistance to late spot disease (Luo et al. 2005; Nobile et al. 2008; Kumar and Kirti, 2011), Aspergillus flavus and A. parasiticus infection (Guo et al. 2008; Guo et al. 2011), and root-knot nematode (Meloidogyne arenaria) parasitization (Guimãraes et al. 2010; Tirumalaraju et al. 2011), and desiccation (Jain et al. 2001; Gopalakrishna et al. 2001) and chilling responsive genes (Tang et al. 2011); in contrast, there is only one report on differentially expressed genes (DEGs) between BW resistant and susceptible peanut genotypes. Peng et al. (2011) identified 119 transcription-derived fragments (TDFs) after root wounding inoculation with R. solanacearum, from Yuanza 9102 (a Spanish type peanut cultivar with BW resistance) and Zhonghua 12 (a susceptible Spanish type peanut cultivar) using cDNA-AFLP and further studied their expression patterns. 98 TDFs were cloned and sequenced, 40 of which were found to have homology to sequences deposited in non-redundant (nr) database of NCBI with known function, while 15 had homology to sequences with unknown function and 43 had no homology.
In the present communication, we reported isolation of DEGs in response to challenge with Ralstonia solanacearum from a large-seeded peanut cultivar Rihua 1 after treatment with the bacterial pathogen suspension using GenefishingTM technology and cloning of full length cDNAs of genes coding for cyclophilin (CyP) and ADP-ribosylation factor (ARF), respectively, by rapid amplification of cDNA ends (RACE).
Peanut variety and bacterial strain. Rihua 1, a Virginia type peanut cultivar with verified high resistance to BW both in field and at laboratory, was kindly provided by Mr. Dian Wen Zhang, a peanut breeder from San Zhuang Town (119º8E, 35º30N), Rizhao City, Shandong Province, China. The causal pathogen, R. solanacearum RZ strain was isolated from diseased peanut plants collected in Rizhao field by the first author and maintained at Laboratory of Biotech Division, Shandong Peanut Research Institute.
Inoculum preparation and inoculation. The bacterial strain was streaked onto TZC Agar (nutrient agar supplemented with 0.05% tetrazolium chloride) (Kelman, 1954) and incubated at 28ºC for 48 to 72 hrs. The bacterial strain was then flushed off the surface of the culture media with sterile double distilled water to prepare a suspension of 109 cfu per millilitre for inoculation.
The peanut seeds were surface sterilized with 75% (v/v) ethanol, and soaked in the bacterial suspension in a Petri dish in a growth chamber (28ºC). No light was provided. Sterile water treatment was used as control.
Isolation of total RNA. Total RNA was isolated from the peanut seeds treated with sterile double distilled water or bacterial suspension for 2, 3, 4 and 5 days, respectively, using RNAprep pure Plant Kit (Tiangen, Beijing, China). Concentration and integrity of the RNA were determined by spectrophotometry and relative intensity of brightness of GelRed (Biotium, CA, USA) stained bands resolved on a 1.2% agarose gel.
Cloning and sequencing of DEGs. DEGs from seeds of Rihua 1 treated with bacterial suspension or water (control) were identified using GenefishingTM DEG Premix Kit (Seegene, Korea) following manufacturers instructions. The reaction mixture for reverse transcription (RT) (20 µl total volume) consisted of 3 µg of total RNA isolated from control or pathogen infected samples, 1 µl of 10 µM dT-ACP1, 4 µl of 5 x RT buffer, 2 µl of 10 mM dNTP mix, 0.5 µl of RNase inhibitor (40 U l-1) (Tiangen, Beijing, China) along with 200 U of M-MLV reverse transcriptase (TaKaRa, Japan). RT was conducted at 42ºC for 90 min, followed by incubation at 70ºC for 15 min to terminate the reaction. First strand cDNA products were then diluted with 80 µl of DNase-free water and directly used in subsequent Genefishing PCR. PCR mixture (20 µl) contained 50 ng of first strand cDNA, 0.5 µM arbitrary ACP, 0.5 µM dT-ACP2 and 2 x SeeAmp ACP Master-mix. PCR program was 94ºC for 5 min, 50ºC for 3 min and 72ºC for 1 min, followed by 40 cycles of 94ºC for 40 sec, 65ºC for 40 sec and 72ºC for 40 sec, and a final extension of 72ºC for 5 min. PCR products were separated on a 2% agarose gel, stained with Gelred and visualized under UV light. Amplicons of interest from treated samples were cloned into a pGM-T vector (Tiangen, Beijing, China), and sequenced by Genscript Inc., Nanjing, China.
Sequence analysis and annotation of DEGs. Cloned DEGs were analyzed using Blastn in search of homologous sequences in nr and EST databases of NCBI. Annotation was performed based on the best match identified by Blastn against nucleonic acid databases at NCBI.
Real time RT-PCR. First strand cDNA was synthesized using total RNA, M-MLV reverse transcriptase (TaKaRa, Japan), oligo-dT, dNTPs and reverse transcriptase buffer. PCR primer pairs were designed with Beacon Designer 7.91 (Premier Biosoft International, Palo Alto, CA, USA) (Table 1). Real time RT-PCR was run on a Lightcycler 2.0 PCR machine (Roche Diagnostics, Penzberg, Germany). Both the PCR program and the component of reaction mixture were the same as previously described by Tang et al. (2011). The PCR program was followed by a melting program of 65ºC to 95ºC at a transition rate of 0.1ºC s-1 with the fluorescence continuously monitored. In each run, a negative control without cDNA template was included to evaluate overall specificity. Reactions were performed in triplicate, and the averages presented. Fold changes in RNA transcripts were calculated by the 2-ΔΔCt method (Livak and Schmittgen, 2001) with β-actin gene as an internal control (Table 1).
5 RACE and isolation of full length cDNA of cyclophilin and ADP-ribosylation factor. 5 cDNA ends were amplified with gene specific primers designed with Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA, USA) using SMARTerTM RACE cDNA Amplification Kit (Clontech, California, USA) according to manufacturers instructions. The full length cDNA sequence of cyclophilin from Rihua 1, designated as CyP, was deduced by aligning 5RACE product and 3d-A13 cDNA sequence. The primer pair, DE-CYP (Table 1), was used to amplify the open reading frame (ORF). Similarly, the full length cDNA sequence of ADP-ribosylation factor from Rihua 1, designated as ARF, was obtained by aligning 5 RACE product and 4d-A9 cDNA sequence. The primer pair, DE-ARF (Table 1), was synthesized to amplify the corresponding ORF.
To obtain the ORFs of the two genes, total RNA was used for reverse transcription in the presence of Oligo-dT, dNTPs and M-MLV reverse transcriptase enzyme (TaKaRa, Japan). First strand cDNA thus synthesized was used as template for subsequent PCR using Taq Platinum Master Polymerase (Tiangen, Beijing, China). PCR products were cloned into a pGM-T vector (Tiangen, Beijing, China) and sequenced by Genscript Inc., Nanjing, China.
Analysis of full length cDNA sequences. After vector and the primer sequences were removed, the cloned sequences were subjected to further analysis. Blastn similarity searches were performed at NCBI website (http://www.ncbi.nlm.nih.gov). Open reading frame (ORF) was analyzed with BioxM, and multiple sequence alignment with DNAStar Lasergene MegAlign 7.1.0 (DNASTAR, Inc., WI, USA). Annotation was performed using BLAST2GO (http://www.blast2go.com/b2ghome).
Isolation and annotation of DEGs from peanut in response to challenge with R. solanacearum
DEGs in Rihua 1 upon infection with R. solanacearum were isolated using Seegenes Genefishing DEG Premix Kit. Partial results were shown in Figure 1. The kit utilizes Annealing Control Primer (ACP) technology (Hwang et al. 2003) suitable for high annealing temperature, thereby essentially eliminating false positives (Kumar and Kirti, 2011). Twenty-five partial cDNAs were cloned and sequenced from Rihua 1 in response to bacterial inoculation (Table 2).
Verification of DEGs
Of the 25 unique cDNA sequences, 11 (44%) had no known function, 4 (16%) were potentially related to diseases defence, 2 (8%), and 1 (4%), 2 (8%), 3 (12%) and 2 (8%) were involved in signal transduction, transcription, protein synthesis and basic metabolism, respectively.
Two of the four transcripts encoding putative diseases defence related proteins, 3d-A13 and 4d-A9 (Table 2) were verified by real time RT-PCR. The results showed that both genes expressed differentially in stressed and control peanut seeds, with relative expression of 1.34 and 0.82 in the 2nd day, 4.88 and 0.28 in 3rd day, 3.69 and 3.49 in 4th day, and 1.97 and 1.11 in 5th day, respectively (Figure 2), indicating that a marked increase in relative expression appeared in the 3rd day for 3d-A13 and 4th day for 4d-A9.
Changes in gene expression in bacterial suspension treated seeds over time were investigated. The expression of the 2 genes in the 3rd, 4th and 5th day relative to that in the 2nd day was shown in Figure 3. Both genes reached peak relative expression in the 4th day.
Isolation and analysis of full length cDNAs of cyclophilin and ADP-ribosylation factor
Considering the fact that the sequence 3d-A13 encoded a cyclophilin-like protein potentially related to diseases defence, 5 cDNA sequence of the gene was amplified with 5RACE (Lane 1 in Figure 4), and full length cDNA obtained. The full length cDNA was 907 bp in length potentially encoding a single polypeptide of 172 amino acids. An 87 bp 5UTR and a 296 bp 3UTR including the poly-A tail were present flanking the ORF. The full length cDNA of cyclophilin (CyP) was submitted to the Genbank under the accession number JN379456. The corresponding ORF amplified with primer pair DE-CYP showed high sequence identity to A. diogoi (data not shown). Multiple sequences alignment of the deduced amino acid sequence of cyclophilin from Rihua 1 in the present study and those from Arabidopsis thaliana, Glycine max and Populus trichocarpa revealed high sequence homology (Figure 5).
4d-A9, the sequence that potentially coded for ADP-ribosylation factor (ARF), was used to obtain the full length cDNA of ARF with 5RACE (Lane 2 in Figure 4). The full length cDNA was 923bp long with a 114 bp 5UTR and a 266 bp 3UTR, potentially encoding a single polypeptide of 118 amino acids. It was submitted to the Genbank under the accession number JN379456. Alignment of the deduced amino acid sequence of ARF from Rihua 1 in the present study with those from other organisms like Arabidopsis thaliana, Populus trichocarpa and Medicago sativa demonstrated that the protein was evolutionarily conserved (Figure 6).
The annotation results of both genes using BLAST2GO were shown in Table 3.
Although several DNA markers related to BW resistance were identified in previous studies, the map distances were too large to be used in peanut breeding programs (Yu et al. 2011). In contrast, differential expression analysis, which can be performed without a population in a short period of time, may result in some genes usable in peanut breeding, as long as their functions are further confirmed by transgenic studies. In the report of Peng et al. (2011), artificially wounded peanut roots were used for inoculation. In the present study, however, the seeds were soaked in the bacterial suspension instead. BW resistance used to be identified by sowing seeds directly in diseased field or by soaking seeds in the bacterial suspension prior to sowing. In either case, calculation of the survival percentage of each genotype at harvest was needed. Consequently, both methods were time-consuming. Recent years, root wounding or cuttings inoculation has been frequently used for evaluation of BW resistance in vegetable crops, and has proved to be rapid and reliable. The root wounding method is also applicable to peanut for evaluation of BW resistance, but it might be inappropriate for isolation of resistance related genes through transcriptome profiling, as roots wounding inevitably leads to the expression of wound-inducible genes unrelated to BW resistance. Up-regulated genes from the present study, however, may contain genes conferring BW resistance, as well as genes unrelated to resistance. In some cases, up-regulation may be a consequence of plant-pathogen interaction, rather than a determinant in resistance (Seevers et al. 1971). If there is a susceptible control for comparison, the unrelated genes may be excluded.
In this report, gene expression patterns of CyP and ARF were studied and full length cDNAs of both genes obtained. Huang et al. (2011) studied the dynamics of R. solanacearum population by stem injection and real time RT-PCR, and concluded that 3-5 day post inoculation was of vital importance for peanut and R. solanacearum interaction; likewise, in our study, expression of the 2 genes during this period was generally higher than that in 2 day, providing supports to Huangs observation.
CyPs are ubiquitous and constitutively expressed. However, they are also stress-responsive proteins, and up-regulated gene expression have been reported in response to abiotic/biotic stresses including heat, cold, salt, wounding, and virus infection (Dubery, 2007). For example, Dubery (2007) observed accumulation of potato CyPs mRNA in response to salicylic acid, Phytophthora infestans elicitor and P. infestans infection, and concluded that CyPs played an important role in plant stress responses. Kumar and Kirti (2011) isolated cyclophilin from A. diogoi inoculated with peanut late spot pathogen, Phaeoisariopsis personata. Constitutive heterogonous expression of the gene in transgenic tobacco enhanced resistance to R. solanacearum (Kumar and Kirti, 2011). We speculated that the CyP isolated from BW resistant cultivar Rihua 1 was also of importance to BW resistance in peanut.
ARFs have been isolated from Arabidopsis, rice, maize and wheat. Their function is mainly involved in mitosis and cell cycle control during seed development and regulation of intracellular transport (Lee et al. 2003, Table 3). Several reports have indicated that ARFs have a role in stress resistance in plants. Lee et al. (2003) reported that over-expression of rice ARF1 gene induced pathogenesis-related (PR) genes and pathogen resistance in tobacco plants, and they deduced that ARF1 might be a component of various plant defense signaling pathways in inducing the expression of a subset of PR genes. Coemans et al. (2008) identified an ARF1 involved in non-host resistance to bacteria and N gene-mediated resistance in Nicotiana benthamiana through high-throughput in planta expression screening. The ARF identified in our study might also have some implications in peanut resistance to the bacterial pathogen, R. solanacearum, whose function still deserves further confirmation.
To summarize, in this study, a total of 25 DEGs upon inoculation of seeds with R. solanacearum suspension were isolated from Rihua 1, a Virginia type BW resistant peanut cultivar by using GenefishingTM technology, and gene expression patterns of CyPand ARF were studied. However, the detailed relationship of these genes with BW resistance in peanut still remains unknown. Further research is under way to elucidate their functions in peanut by antisense/RNAi technology or over-expression, as a high-efficiency genotype-independent transgenic protocol has already been developed at our lab (Li et al. 2011). The outcome of the present study may therefore provide candidate genes potentially useful for breeding high yielding BW resistant peanut cultivars.
COEMANS, B.; TAKAHASHI, Y.; BERERICH, T.; ITO, A.; KANZAKI, H.; MATSUMURA, H.; SAITOH, H.; TSUDA, S.; KAMOUN, S.; SÁGI, L.; SWENNEN, R. and TERAUCHI, R. (2008). High-throughput in planta expression screening identifies an ADP-ribosylation factor (ARF1) involved in non-host resistance and R gene-mediated resistance. Molecular Plant Pathology, vol. 9, no. 1, p. 25-36. [CrossRef] [ Links ]
GOPALAKRISHNA, R.; KUMAR, G.; KRISHNAPRASAD, B.T.; MATHEW, M.K. and KUMAR, M.U. (2001). A stress-responsive gene from groundnut, Gdi-15, is homologous to flavonol 3-O-glucosyltransferase involved in anthocyanin biosynthesis. Biochemical and Biophysical Research Communications, vol. 284, no. 3, p. 574-579. [CrossRef] [ Links ]
GUIMARÃES, P.M.; BRASILEIRO, A.C.M.; PROITE, K.; DE ARAÚJO, A.C.G.; LEAL-BERTIOLI, S.C.M.; PIC-TAYLOR, A.; DA SILVA, F.R.; MORGANTE, C.V.; RIBEIRO, S.G. and BERTIOLI, D.J. (2010). A study of gene expression in the nematode resistant wild peanut relative, Arachis stenosperma, in response to challenge with Meloidogyne arenaria. Tropical Plant Biology, vol. 3, no. 4, p. 183-192. [CrossRef] [ Links ]
GUO, B.Z.; CHEN, X.P.; DANG, P.; SCULLY, B.T.; LIANG, X.Q.; HOLBROOK, C.C.; YU, J. and CULBREATH, A.K. (2008). Peanut gene expression profiling in developing seeds at different reproduction stages during Aspergillus parasiticus infection. BMC Development Biology, vol. 8, p. 12. [CrossRef] [ Links ]
GUO, B.Z.; FEDOROVA, N.D.; CHEN, X.P.; WAN, C.; WANG, L.; NIERMAN, W.C.; BHATNAGAR, D. and YU, J. (2011). Gene expression profiling and identification of resistance genes to Aspergillus flavus infection in peanut through EST and microarray strategies. Toxins, vol. 3, no. 7, p. 737-753. [CrossRef] [ Links ]
HUANG, J.Q.; YAN, L.Y.; YE, X.W.; LEI, Y. and LIAO, B.S. (2011). Development of Ralstonia solanacearum quantification method and its application in peanut with bacterial wilt disease. Scientia Agricultura Sinica, vol. 44, no. 1, p. 58-66. [CrossRef] [ Links ]
HWANG, I.T.; KIM, Y.J.; KIM, S.H.; KWAK, C.I.; GU, Y.Y. and CHUN, J.Y. (2003). Annealing control primer system for improving specificity of PCR amplification. BioTechniques, vol. 32, no. 6, p. 1180-1184. [ Links ]
JAIN, A.K.; BASHA, S.M. and HOLBROOK, C.C. (2001). Identification of drought-responsive transcripts in peanut (Arachis hypogaea L.). Electronic Journal of Biotechnology, vol. 4, no. 2. [CrossRef] [ Links ]
KELMAN, A. (1954). The relationship of pathogenicity in Ralstonia solanacearum to colony appearance on a tetrazolium medium. Phytopathology, vol. 44, p. 693-695. [ Links ]
KUMAR, K.R.R. and KIRTI, P.B. (2011). Differential gene expression in Arachis diogoi upon interaction with peanut late leaf spot pathogen Phaeoisariopsis personata and characterization of a pathogen induced cyclophilin. Plant Molecular Biology, vol. 75, no. 4-5, p. 497-513. [CrossRef] [ Links ]
LEE, W.Y.; HONG, J.K.; KIM, C.Y.; CHUN, H.J.; PARK, H.C.; KIM, J.C.; YUN, D-J; CHUNG, W.S.; LEE, S-H; LEE, S.Y.; CHO, M.J. and LIM, C.O. (2003). Over-expressed rice ADP-ribosylation factor 1 (RARF1) induces pathogenesis-related genes and pathogen resistance in tobacco plants. Physiologia Plantarum, vol. 119, no. 4, p. 573-581. [CrossRef] [ Links ]
LI, G.J.; WANG, X.Z; TANG, Y.Y.; WU, Q. and WANG, C.T. (2011). Isolation of peanut DNA sequences from transgenic seeds with flanking T-DNA insertion sites elevated protein content. Journal of Peanut Science, vol. 40, no. 3, p. 1-6. [ Links ]
LUO, M.; DANG, P.; BAUSHER, M.G.; HOLBROOK, C.C.; LEE, R.D.; LYNCH, R.E. and GUO, B.Z. (2005). Identification of transcripts involved in resistance responses to leaf spot disease caused by Cercosporidium personatum in peanut (Arachis hypogaea). Phytopathology, vol. 95, no. 4, p. 381-387. [CrossRef] [ Links ]
NOBILE, P.M.; LOPES, C.R.; BRASALOBRES-CAVALLARI, C.B.; QUECINI, V.; COUTINHO, L.L.; HOSHINO, A.A. and GIMENES, M.A. (2008). Peanut genes identified during initial phase of Cercosporidium personatum infection. Plant Science, vol. 174, no. 1, p. 78-87. [CrossRef] [ Links ]
PENG, W.-F.; LV, J.-W.; REN, X.P.; HUANG, L.; ZHAO, X.-Y.; WEN, Q.-G. and JIANG, H.-F. (2011). Differential expression of genes related to bacterial wilt resistance in peanut (Arachis hypogaea L.). Hereditas (Bejing), vol. 33, no. 4, p. 389-396. [CrossRef] [ Links ]
TANG, Y.Y.; WANG, C.T.; YANG, G.P.; FENG, T.; GAO, H.Y.; WANG, X.Z.; CHI, X.T.; WU, Q. and CHEN, D.X. (2011). Identification of chilling-responsive transcripts in peanut (Arachis hypogaea L.). Electronic Journal of Biotechnology, vol. 14, no. 5. [CrossRef] [ Links ]
TIRUMALARAJU, S.V.; JAIN, M. and GALLO, M. (2011). Differential gene expression in roots of nematode-resistant and - susceptible peanut (Arachis hypogaea) cultivars in response to early stages of peanut root-knot nematode (Meloidogyne arenaria) parasitization. Journal of Plant Physiology, vol. 168, no. 5, p. 481-492. [CrossRef] [ Links ]
YABUUCHI, E.; KOSAKO, Y.; YANO, I.; HOTTA, H. and NISHIUCHI, Y. (1995). Transfer of two Burkholderia and an Alcaligenes species to Ralstonia gen. nov.: Proposal of Ralstonia pickettii (Ralston, Palleroni and Doudoroff 1973) comb. nov., Ralstonia solanacearum (Smith 1896) comb. nov. and Ralstonia eutropha (Davis 1969) comb. nov. Microbiology and Immunology, vol. 39, no. 11, p. 897-904. [ Links ]
YU, S.L.; WANG, C.T.; YANG, Q.L.; ZHANG, D.X.; ZHANG, X.Y.; CAO, Y.L.; LIANG, X.Q. and LIAO, B.S. (2011). Peanut Genetics and Breeding in China. Shanghai Science and Technology Press, 565 p. ISBN 978-7-5478-0610-4. [ Links ]
Note: Electronic Journal of Biotechnology is not responsible if on-line references cited on manuscripts are not available any more after the date of publication.