Analysis of sea-island cotton and upland cotton in response to Verticillium dahliaeinfection by RNA sequencing
© Sun et al.; licensee BioMed Central Ltd. 2013
Received: 28 July 2013
Accepted: 26 November 2013
Published: 5 December 2013
Cotton Verticillium wilt is a serious soil-borne vascular disease that causes great economic loss each year. However, due to the lack of resistant varieties of upland cotton, the molecular mechanisms of resistance to this disease, especially to the pathogen Verticillium dahliae, remain unclear.
We used the RNA-seq method to research the molecular mechanisms of cotton defence responses to different races of Verticillium dahliae by comparing infected sea-island cotton and upland cotton. A total of 77,212 unigenes were obtained, and the unigenes were subjected to BLAST searching and annotated using the GO and KO databases. Six sets of digital gene expression data were mapped to the reference transcriptome. The gene expression profiles of cotton infected with Verticillium dahliae were compared to those of uninfected cotton; 44 differentially expressed genes were identified. Regarding genes involved in the phenylalanine metabolism pathway, the hydroxycinnamoyl transferase gene (HCT) was upregulated in upland cotton whereas PAL, 4CL, CAD, CCoAOMT, and COMT were upregulated in sea-island cotton. Almost no differentially expressed genes in this pathway were identified in sea-island cotton and upland cotton when they were infected with V. dahliae V991 and V. dahliae D07038, respectively.
Our comprehensive gene expression data at the transcription level will help elucidate the molecular mechanisms of the cotton defence response to V. dahliae. By identifying the genes involved in the defence response of each type of cotton to V. dahliae, our data not only provide novel molecular information for researchers, but also help accelerate research on genes involved in defences in cotton.
KeywordsCotton Verticillium dahliae RNA sequencing Transcriptome analysis Mechanism of resistance
Cotton is the nature fibre crop with the greatest economic importance. Cotton Verticillium wilt (CVW) is a serious soil-borne vascular disease caused by Verticillium dahliae. The disease can cause cotton yellowing, wilt, defoliating, and finally death . CVW was first reported in upland cotton in the United States, and spread to China on imported American cotton in 1935[3–5]. At present, more than 200 million hectares of cotton are subject to CVW, and the economic loss is tremendous every year, especially defoliated CVW . Currently, no fungicides are available to cure commercial upland cotton (Gossypium hirsutum L., the main cultivated cultivar) once they become infected, although the sea-land cotton species (Gossypium barbadense L., not the main cultivated cultivar) is immune to V. dahliae. Because of the difficulty obtaining a new highly resistant upland cotton breed by cross-breeding of upland cotton and sea-island cotton, no completely suitable high-resistance upland cotton variety exists. The breeding of disease-resistant cultivars remains the primary control method for cotton. At present, the primary method of breeding for disease-resistant cotton is crossbreeding of resistant cotton, which has been used to create a series of varieties with resistance to CVW [7, 8].
In recent years, researchers have investigated the molecular mechanisms of resistance to V. dahliae with the aim of ultimately using genetic engineering to breed cultivars resistant to CVW. Resistance to CVW is controlled by a single dominant gene in sea-island cotton, whereas that in upland cotton may be controlled by either a single dominant gene or two dominant genes [5, 8–13]. When cotton is infected by the CVW pathogen, the defence system is activated, generating a series of cascade reactions. Three aspects of this system, tissue resistance, physiological and biochemical resistance, and micro-organism resistance, are involved in the mechanism of resistance [14, 15]. Several V. dahliae resistance-related genes have been isolated and studied, including the following four groups: cinnamyl alcohol dehydrogenase (CAD), farnesyl diphosphate synthase, isopentenyl diphosphate isomerase, and 3-hydroxy-3-methylglutaryl coenzyme A reductase, the key regulatory enzymes in the synthesis of phytoalexin, which can kill pathogens through increased local concentrations in the plant [16, 17] (Group 1); lipid transfer proteins, chitinase, β-1,3 glucanase, nonexpressor of pathogenesis-related genes 1, thionin, gastrodia antifungal protein, and lectin-like protein, which may play important roles in the response to pathogen infections in cotton plants [18, 19] (Group 2); phenylalanine ammonia lyase (PAL), cinnamate-4-hydroxylase, peroxisome, polygalacturonase-inhibiting proteins, and laccase, which are defence proteins that can slow down the spread of pathogens [20, 21] (Group 3); and resistance gene analogues and defence gene analogues that have been cloned and further transformed by genetic engineering for functional research, such as lipid transfer protein, nonrace-specific disease resistance 1, MAP kinase kinase 2, and others [9, 22] (Group 4).
The phenylpropanoid pathway plays a key role in cotton resistance to V. dahliae[8, 23–25]. Even so, the molecular mechanisms of V. dahliae resistance are unclear. In recent years, novel, high-throughput, deep-sequencing transcriptome analysis, termed RNA-seq, has made it possible to efficiently generate large-scale expressed sequence tag (EST) libraries and improved the speed of gene discovery . Although RNA-seq technologies have been used for researching gene expression profiles in the resistance response of sea-island cotton to V. dahliae, the phenylpropanoid metabolism pathway and genes associated with lignin metabolism has been found to play a central role . Whether a different response mechanism to the wilt fungus V. dahliae exists in cotton of different races and whether different response mechanisms exist in different cotton species (G. hirsutum L. and G. barbadense L.) are questions for further research. Therefore, in this study, four samples, including ZhongZhiMian KV-1 (G. hirsutum L.) and XinHai 15 (G. barbadense L.) infected with V. dahliae strain V991 (highly toxic) and D07038 (intermediately toxic), respectively, and two uninfected samples include ZhongZhiMian KV-1 (G. hirsutum L.) and XinHai 15 (G. barbadense L.) were sequenced in order to addressed the above issues using RNA-seq technologies.
Illumina sequencing and sequence assembly
Functional annotation and classification
All of the unigenes were compared to the sequences in public databases, including the NCBI nonredundant protein (nr) database, the Clusters of Orthologous Groups (COG) database, the Swiss-Prot protein database, and the KEGG database, using BLASTX with a cutoff e-value of 10-5. A total of 66,208 unigenes (78% of all unigenes) returned a significant BLAST result.
Pathway annotation to KEGG
No. of unigenes (28395)
Biosynthesis of secondary metabolites
Plant hormone signal transduction
Ribosome biogenesis in eukaryotes
Protein processing in endoplasmic reticulum
mRNA surveillance pathway
Starch and sucrose metabolism
Ubiquitin mediated proteolysis
Ether lipid metabolism
Circadian rhythm - plant
Amino sugar and nucleotide sugar metabolism
Stilbenoid, diarylheptanoid and gingerol biosynthesis
Limonene and pinene degradation
Phosphatidylinositol signaling system
DGE library sequencing
Statistics of DGE sequencing
Total mapped reads
<=2 bp mismatch
Total unmapped reads
Different components of the raw reads in each sample
Clean reads only
Gene expression variation among the different samples
Genes related to resistance to CVW resistance
Based on the above comparisons, the DEGs were functionally analysed using MapMan software. The results showed that a total of 1975 genes occurred simultaneously in the four comparisons above (a–d), including 1535 upregulated genes and 440 downregulated genes. These genes were principally associated with cell wall, lipids, and secondary metabolism and so on (Additional file 1). To identify cotton genes most associated with resistance to CVW, an absolute value of the fold change times > =5 was used as a threshold, resulting in 44 DEGs being screened out (including 41 upregulated and 3 downregulated genes, Table 4).
Validation of RNA-Seq-based gene expression
Cotton is an economically important crop, and cotton CVW causes significant economic losses. Although some cotton cultivars resistant to V. dahliae have been generated using conventional breeding techniques, the molecular mechanisms of disease resistance in cotton are still unclear. Xu et al.  found that lignin metabolism is a key pathway in the resistance of cotton, but the authors mainly analysed the RNA-seq results by analysing the resistance mechanism of the resistant Gossypium barbadense in contrast to the susceptible Gossypium hirsutum. Are other key genes involved in the resistance of cotton? In this study, we used the commercial upland cotton (Gossypium hirsutum) variety ZhongZhiMian KV-1 and the sea-island cotton Xinhai 15, both highly resistant to V. dahliae, to study the resistance relationships between different cotton species and different V. dahliae strains.
Using transcriptome sequence analysis, 56,836,100 reads were obtained, corresponding to about 4.7 Gb of raw sequence data. Then the raw sequences were further assembled using the short-reads assembly program Trinity and used to predict unigenes. In spite of the genome sequence of a diploid cotton Gossypium raimondi had been published, whose progenitor is the putative contributor of the D-subgenome to cotton species Gossypium barbadense and Gossypium hirsutum[28, 29], due to the lack of A-genome sequence information, the D-genome sequence were not used in our bioinformation analysis in this study, and the quality of these assembled unigenes was confirmed by qRT-PCR. The results showed that signals for 19 out of 20 tested genes were detected, demonstrating the high quality of the assembled results, and also indicating that relatively short reads from Illumina sequencing can be effectively assembled and used for novel gene discovery. The 66208 predicted unigenes were subjected to BLAST annotation; 78% of the unigenes returned a significant BLAST result. About 20% of the predicted unigenes could not be assigned annotation information, which may be caused by the limited information about the genomes or transcriptomes of cotton and its related species. Some unigenes were too short to allow statistically meaningful matches. From the KEGG annotation results, 2238 unigenes (7.88% of all KEGG annotation unigenes) were found to be involved in the group 'plant-pathogen interaction’, which may be associated with the defence response of cotton to V. dahliae. This group was the fourth largest group in our results, which suggests that it is feasible through the immense capacity of Illumina sequencing to discover genes related to cotton resistance to V. dahliae in metabolic pathways.
To understand the regulatory pathways and molecular mechanisms of resistance to V. dahliae, we created six DGE libraries to analyse gene expression patterns following inoculation with V. dahliae. About 80% of the reads were mapped to our transcriptome reference database. The high proportion of mapped reads supports the validity of our analysis. The quality of the DGE libraries was also confirmed by qRT-PCR analysis. Because the cotton disease resistance-related genes may begin to be expressed after inoculation with V. dahlia, cotton seedlings of different species inoculated with different races of V. dahliae were compared to controls. Different genes may be involved in disease resistance, and the common genes in the four comparison groups reduced the number of DEGs that may be related to cotton disease resistance. There were no replicated data using the Audic-Claverie method in DGE data analysing, but the results are credible because we applied a series of correction method, such as: a false discovery rate of <0.001 and an absolute value of the log 2 ratio >1 as the threshold were used for judging the significance of the gene expression differences, and combined with the quality results of the DGE libraries been further confirmed by qRT-PCR.
The inoculated V. dahliae cotton samples compared with uninoculated samples exhibited much more upregulated genes and downregulated genes, this possibly indicates that V. dahliae activates cotton resistance genes during the infection process. Furthermore, the number of upregulated and downregulated genes of G. barbadense was much greater than that in G. hirsutum (both inoculated with V. dahliae races V991 and D07038) in the four compared groups, suggesting that the resistance genes to V.dahliae of the sea-island cotton G. barbadense is more than that of the upland cotton G. hirsutum.
Different expression genes screening out
cell wall.cellulose synthesis.cellulose synthase
cellulose synthase-like D5
TPS14 (TERPENE SYNTHASE 14)
early-responsive to dehydration protein-related
early-responsive to dehydration protein-related
misc.UDP glucosyl and glucoronyl transferases
XT1 (XYLOSYLTRANSFERASE 1)
fatty acid (omega-1)-hydroxylase/ oxygen binding
LTL1 (LI-TOLERANT LIPASE 1)
RNA.regulation of transcription.AP2/EREBP, APETALA2/Ethylene-responsive element binding protein family
ABR1 (ABA REPRESSOR1)
RNA.regulation of transcription.AP2/EREBP, APETALA2/Ethylene-responsive element binding protein family
ap2/EREBP transcription factor
RNA.regulation of transcription.bHLH,Basic Helix-Loop-Helix family
basic helix-loop-helix (bHLH) family protein
RNA.regulation of transcription.C2C2(Zn) CO-like, Constans-like zinc finger family
JAZ10 (JASMONATE-ZIM-DOMAIN PROTEIN 10)
RNA.regulation of transcription.MYB-related transcription factor family
ATRL5 (ARABIDOPSIS RAD-LIKE 5)
calmodulin-binding family protein
NAC domain-containing protein 68
amino acid/polyamine transporter II
BOR4 (REQUIRES HIGH BORON 4)
Although the molecular functions of the cotton defence response to V. dahliae and the associated signal transduction pathways remain largely unknown, the present transcriptome analysis provides valuable information regarding cotton resistance, which may facilitate future investigations of the detailed regulatory mechanisms and pathways. Using Illumina sequencing technology, we obtained about 5 Gb of raw sequence data and predicted 77,212 assembled unigenes, 66,208 of which were able to be annotated. Finally 44 differentially expressed genes were identified. The results suggest that HCT in upland cotton and CCoAOMT, PAL, 4CL, CAD, and COMT in sea-island cotton are involved in cotton defence response to V. dahliae. We believe that our data not only provide novel molecular information, but will also help to accelerate research into gene functions of cotton defence responses.
Preparation of material
One highly aggressive strain of the defoliating fungus V. dahliae, V991, and one intermediately aggressive strain, D07038, from the Institute of Plant Protection (IPPCAAS) and the Cotton Institute of the Chinese Academy of Agriculture Sciences(CICAAS), respectively, were used for inoculation . The cotton varieties ZhongZhiMian KV-1 (G. hirsutum L., from IPPCAAS) and XinHai 15 (G. barbadense L.,from the Economic Crop Institute of the Xinjiang Academy of Agriculture Sciences,ECIXAAS, and CICAAS) resistant to V. dahliae were grown in sterilized soil (a mix of peat and sawdust) in a natural environment . Each seedling was inoculated with 10 mL V. dahliae spore suspension of 2 × 107 spores per mL by watering injured roots from plants at the two-true-leaf growth stage. Control plants were not inoculated but were otherwise treated and sampled with distilled water in the same way. Plants from six individual seedlings were collected for each treatment at each sampling time point (including 24, 48, and 96 h after inoculation) after washed by 75% alcohol and sterile water.
cDNA preparation for Illumina sequencing
Total RNA was extracted using a Plant RNA EASYspin Plus Kit (aidlab, Peking, China) according to the manufacturer’s instructions. Six RNA samples (for each sample, all time points of each individual treatment were mixed), including four samples of ZhongZhiMian KV-1 (G. hirsutum) and XinHai 15 (G. barbadense) infected with V. dahliae strain V991 (highly toxic) and D07038 (intermediately toxic), respectively. Two plant samples treated with sterile water were further verified for the quantity of RNA using an ultraviolet spectrometer and electrophoresis on a denaturing formaldehyde agarose gel. The samples were as follows: 1, ZhongZhiMian KV-1 infected with D07038 (named: KV1/D); 2, XinHai 15 infected with D07038(named: XH/D); 3, ZhongZhiMian KV-1 infected with V991(named: KV1/V); 4, XinHai 15 infected with V991(named: XH/V); 5, uninfected ZhongZhiMian KV-1(named: KV1/U); and 6, uninfected XinHai 15(named: XH/U). All samples were mixed equally for transcriptome sequencing, but subjected individually to DGE sequencing. Oligo(dT) beads were used to isolate poly(A) + mRNA from total RNA. Fragmentation buffer was added to disrupt the mRNA into short fragments. These short fragments were used as templates with random hexamer primer to synthesise first-strand cDNA. Second-strand cDNA was synthesised by adding buffer, dNTPs, RNase, and DNA polymerase I. The resulting short fragments were purified with a QiaQuick PCR extraction kit and resolved with EB buffer for end reparation and addition of a poly(A) tail. Next, the short fragments were connected with sequencing adapters. Following agarose gel electrophoresis, suitable fragments were selected as templates for PCR. Finally, the library was sequenced using an Illumina HiSeq™ 2000.
Analysis of transcriptome sequencing data
The raw data from the images were collected using Illumina GA Pipeline 1.6 by removing low-quality reads (reads with unknown sequences 'N’), adaptor sequence fragments, and empty reads. Next, de novo assembly of the transcriptome into unigenes was carried out with Trinity, a short-read assembly program . In brief, Trinity firstly combines reads with certain length of overlap to form longer fragments, which are called contigs. Then the reads are mapped back to contigs; with paired-end reads it is able to detect contigs from the same transcript as well as the distances between these contigs. Finally, Trinity connects the contigs, and get sequences that cannot be extended on either end. Such sequences are defined as Unigenes. Subsequently, a Basic Local Alignment Search Tool (BLAST) BLASTx alignment (e-value < 0.00001) was performed for the unigenes with protein databases, including the non-redundant (nr), Swiss-Prot, Kyoto Encyclopaedia of Genes and Genomes (KEGG), and COG databases, and the best alignments were used to decide the sequence direction of the unigenes. If the results from the different databases conflicted with each other, a priority order of nr, Swiss-Prot, KEGG, and COG was followed to decide the sequence direction of the unigenes. When a unigene was not aligned in any of the above databases, ESTScan was used to predict its coding regions and decide its sequence direction . In the final step, using nr annotation, the Blast2GO program was used to obtain the Gene Ontology (GO) and KEGG annotations for the unigenes . After the GO annotation was obtained for each unigene, WEGO software was used to classify the unigenes by function and to determine the distribution of gene functions in the species at the macro level . All raw transcriptome data have been deposited in SRA (NCBI).
Analysis and mapping of DGE reads
The raw image data were transformed by base calling into sequence data. To map the DGE reads, the sequenced raw data were filtered to remove adaptor sequences, low-quality sequences (reads with unknown sequences 'N’), empty reads. For reads annotation, clean reads sequences were mapped to our transcriptome reference database, allowing no more than two nucleotide mismatch [43, 47]. Gene coverage is the percentage of a gene covered by reads. This value is equal to the ratio of the base number in a gene covered by unique mapping reads to the total bases number of that gene. For gene expression analysis, the number of expressed sequence was calculated and normalised to the number of Reads per Kb per Million reads (RPKM).
Identification of differentially expressed genes
To compare the differences in gene expression, the reads frequency in each DGE library was statistically analysed according to the method of Audic and Claverie . We used a false discovery rate of <0.001 and an absolute value of the log 2 ratio >1 as the threshold for judging the significance of the gene expression differences . To determine the genes related to disease resistance, the inoculated V. dahliae cotton samples (samples 1 to 4) were compared to water-inoculated cotton samples (samples 5 and 6). Next, the differentially expressed genes were subjected to GO and KEGG Ontology (KO) enrichment analysis . At the same time, the differentially expressed genes sequences of each group were uploaded to the Mercator web application to assign MapMan bins, and then differentially expressed genes between the four compared groups of related pathways were screened and analysed using MapMan software .
Data validation by qPCR
primers for qRT-PCR
Celullose synthase-like D protein
Terpene synthase 3
Early-responsive to dehydration protein-related
Early-responsive to dehydration protein-related
AP2/ERF domain-containing transcription factor
AP2/ERF domain-containing transcription factor
Symbiotic ammonium transporter
Protein TIFY 9
MYB transcription factor
Mitogen-activated protein kinase kinase kinase 17
F-box family protein-like
NAC domain protein
Amino acid transporter
Anion exchanger family protein
We acknowledge the Beijing Genomics institute at Shenzhen for its assistance in original data processing and related bioinformatics analysis and thank the National Medium-term Gene Bank of Cotton in China, Professor Guiliang Jian in IPP, Professor Xueyuan Li in ECIXAAS, Dr.Heqin Zhu in CICAAS for the resource materials supply. This work was supported by the Natural Science Foundation of China (NSFC grant No. 31071461) and Chongqing Natural Science Foundation (cstc2012jjA80037;cstc2011AB1095) and Chongqing Municipal Commission of Education (KJ130502, KJ110506).
- Sal'kova EG, Guseva NN: The role of pectolyitc enzymes of the Verticillium dahliae fungus in the development of cotton wilt. Dokl Akad Nauk SSSR. 1965, 163: 515-522.PubMedGoogle Scholar
- Sink KC, Grey WE: A root-injection method to assess Verticillium wilt resistance of peppermint (Mentha × piperita L.) and its use in identifying resistant somaclones of cv. Black Mitcham. Euphytica. 1999, 106 (3): 223-230. 10.1023/A:1003591908308.View ArticleGoogle Scholar
- Du W-S, Du X, Ma Z: Progress of Inheritance and Molecular Biology of Verticillium wilt resistance in Cotton. Cotton Sci. 2002, 14 (5): 55-61.Google Scholar
- Bugbee W, Sappenfield W: Effect of Verticillium wilt on cotton yield, fiber properties, and seed quality. Crop Sci. 1970, 10 (6): 649-652. 10.2135/cropsci1970.0011183X001000060011x.View ArticleGoogle Scholar
- Cai Y, Xiaohong H, Mo J, Sun Q, Yang J, Liu J: Molecular research and genetic engineering of resistance to Verticillium wilt in cotton: a review. Afr J Biotechnol. 2009, 8 (25): 7363-7372.Google Scholar
- Jing YL, Liu YB, Fan WF, Xiao BC: Advance in the Study on Verticillium Wilt of Cotton and It’s Breeding for Resistance. Acta Agric Boreali-occidentalis Sin. 1999, 8 (3): 106-110.Google Scholar
- Jian G, Lu M: Study on the method of breeding cotton for resistance to Verticillium dahliae. Acta Phytopathologica Sin. 2004, 34 (4): 70-74.Google Scholar
- Ma C, Jian G, Zheng C: The advances in cotton breeding resistance to Fusarium and Verticillium wilts in China during past fifty years. Agric Sci China. 2002, 35 (5): 508-513.Google Scholar
- Qi J, Ma C, Zhao L, Liu S: Study on Heredity of Verticillium wilt resistance of G. barbadense L. Acta Gossypii Sin. 2000, 12 (4): 169-171.Google Scholar
- Barrow JR: Heterozygosity in inheritance of Verticillium wilt tolerance in cotton. Phytopathology. 1970, 60 (2): 301-303. 10.1094/Phyto-60-301.View ArticleGoogle Scholar
- Verhalen LM, Brinkerhoff L, Fun K-C, Morrison WC: A quantitative genetic study of Verticillium wilt resistance among selected lines of upland cotton. Crop Sci. 1971, 11 (3): 407-412. 10.2135/cropsci1971.0011183X001100030029x.View ArticleGoogle Scholar
- Jiang F, Zhao J, Zhou L, Guo W, Zhang T: Molecular mapping of Verticillium wilt resistance QTL clustered on chromosomes D7 and D9 in upland cotton. Sci China Ser C, Life Sci / Chin Acad Sci. 2009, 52 (9): 872-884. 10.1007/s11427-009-0110-8.View ArticleGoogle Scholar
- Wang L, Dai X: Progress on Molecular Research of Cotton Verticillium wilt resistances. Mol Plant Breed. 2003, 1 (1): 97-102.Google Scholar
- Smit F, Dubery IA: Cell wall reinforcement in cotton hypocotyls in response to a Verticillium dahliae elicitor. Phytochemistry. 1997, 44 (5): 811-815. 10.1016/S0031-9422(96)00595-X.View ArticleGoogle Scholar
- Zhou Y, Du H, Yuan H, Zhang Y, Zhu B: Isolation and purification of antifungal protein from Paenibacillus to Verticillium dahliae. Cotton Sci. 2007, 19 (2): 98-101.Google Scholar
- Gaffe J, Bru J-P, Causse M, Vidal A, Stamitti-Bert L, Carde J-P, Gallusci P: LEFPS1, a tomato farnesyl pyrophosphate gene highly expressed during early fruit development. Plant Physiol. 2000, 123 (4): 1351-1362. 10.1104/pp.123.4.1351.PubMed CentralView ArticlePubMedGoogle Scholar
- Mannanov RN: The use of natural bio-agents for the control of cotton phytopathogens. Meded Rijksuniv Gent Fak Landbouwkd Toegep Biol Wet. 2001, 66 (3a): 183-186.PubMedGoogle Scholar
- Cheng H, Jian G, Ni W, Yang H, Z-x WANG, W-j SUN, B-l ZHANG, WANG X-f MAC, S-r JIA: Increase of Fusarium-and Verticillium-resistance by transferring chitinase and glucanase gene into cotton. Sci Agric Sin. 2005, 38 (6): 1160-1166.Google Scholar
- Zhang Y, Wang X, Cheng C, Gao Q, Liu J, Guo X: Molecular cloning and characterization of GhNPR1, a gene implicated in pathogen responses from cotton (Gossypium hirsutum L.). Biosci Rep. 2008, 28: 7-14. 10.1042/BSR20070028.View ArticlePubMedGoogle Scholar
- Smith CG, Rodgers MW, Zimmerlin A, Ferdinando D, Bolwell GP: Tissue and subcellular immunolocalisation of enzymes of lignin synthesis in differentiating and wounded hypocotyl tissue of French bean (Phaseolus vulgaris L.). Planta. 1994, 192 (2): 155-164. 10.1007/BF01089030.View ArticleGoogle Scholar
- H-k HUANG, Y-y WU, W-x XU, Q-y ZHOU: Cloning, Sequencing and Expression Profile Analysis of A Polygalacturonase-inhibiting Protein (PGIP) Gene in Cotton. J Henan Agric Sci. 2008, 6: 012-Google Scholar
- Wang HM, Lin ZX, Zhang XL, Chen W, Guo XP, Nie YC, Li YH: Mapping and quantitative trait loci analysis of Verticillium wilt resistance genes in cotton. J Integr Plant Biol. 2008, 50 (2): 174-182. 10.1111/j.1744-7909.2007.00612.x.View ArticlePubMedGoogle Scholar
- Pomar F, Novo M, Bernal MA, Merino F, Barceló AR: Changes in stem lignins (monomer composition and crosslinking) and peroxidase are related with the maintenance of leaf photosynthetic integrity during Verticillium wilt in Capsicum annuum. New Phytol. 2004, 163 (1): 111-123. 10.1111/j.1469-8137.2004.01092.x.View ArticleGoogle Scholar
- Gayoso C, Pomar F, Novo-Uzal E, Merino F, de Ilárduya ÓM: The Ve-mediated resistance response of the tomato to Verticillium dahliae involves H2O2, peroxidase and lignins and drives PAL gene expression. BMC Plant Biol. 2010, 10 (1): 232-10.1186/1471-2229-10-232.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu L, Zhu L, Tu L, Liu L, Yuan D, Jin L, Long L, Zhang X: Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. J Exp Bot. 2011, 62 (15): 5607-5621. 10.1093/jxb/err245.PubMed CentralView ArticlePubMedGoogle Scholar
- Ansorge WJ: Next-generation DNA sequencing techniques. New Biotechnol. 2009, 25 (4): 195-203. 10.1016/j.nbt.2008.12.009.View ArticleGoogle Scholar
- Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32 (suppl 1): D277-D280.PubMed CentralView ArticlePubMedGoogle Scholar
- Andrew H, Paterson JFW, Heidrun G, Hui G, Jerry J, Dianchuan J, Danny L, Kurtis C, Showmaker S, Shengqiang S, Joshua U, Mi-jeong Y, Robert B, Wei C, Adi D-F, Mary V, Duke LG, Jane G, Corrinne G, Kara G, Guanjing H, Tae-ho L, Jingping L, Lifeng L, Tao L, Barry S, Marler J, Justin , Page , Alison , et al: Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature. 2012, 492 (7429): 5-Google Scholar
- Wang K, Wang Z, Li F, Ye W, Wang J, Song G, Yue Z, Cong L, Shang H, Zhu S, et al: The draft genome of a diploid cotton Gossypium raimondii. Nat Genet. 2012, 44 (10): 1098-1103. 10.1038/ng.2371.View ArticlePubMedGoogle Scholar
- Burton RA, Wilson SM, Hrmova M, Harvey AJ, Shirley NJ, Medhurst A, Stone BA, Newbigin EJ, Bacic A, Fincher GB: Cellulose synthase-like CslF genes mediate the synthesis of cell wall (1, 3; 1, 4)-ß-D-glucans. Science. 2006, 311 (5769): 1940-1942. 10.1126/science.1122975.View ArticlePubMedGoogle Scholar
- Doblin MS, Pettolino FA, Wilson SM, Campbell R, Burton RA, Fincher GB, Newbigin E, Bacic A: A barley cellulose synthase-like CSLH gene mediates (1,3;1,4)-beta-D-glucan synthesis in transgenic Arabidopsis. Proc Natl Acad Sci U S A. 2009, 106 (14): 5996-6001. 10.1073/pnas.0902019106.PubMed CentralView ArticlePubMedGoogle Scholar
- Li M, Xiong G, Li R, Cui J, Tang D, Zhang B, Pauly M, Cheng Z, Zhou Y: Rice cellulose synthase-like D4 is essential for normal cell‒wall biosynthesis and plant growth. Plant J. 2009, 60 (6): 1055-1069. 10.1111/j.1365-313X.2009.04022.x.View ArticlePubMedGoogle Scholar
- Vega-Sanchez ME, Verhertbruggen Y, Christensen U, Chen X, Sharma V, Varanasi P, Jobling SA, Talbot M, White RG, Joo M, et al: Loss of Cellulose synthase-like F6 function affects mixed-linkage glucan deposition, cell wall mechanical properties, and defense responses in vegetative tissues of rice. Plant Physiol. 2012, 159 (1): 56-69. 10.1104/pp.112.195495.PubMed CentralView ArticlePubMedGoogle Scholar
- Rose JK, Braam J, Fry SC, Nishitani K: The XTH family of enzymes involved in xyloglucan endotransglucosylation and endohydrolysis: current perspectives and a new unifying nomenclature. Plant Cell Physiol. 2002, 43 (12): 1421-1435. 10.1093/pcp/pcf171.View ArticlePubMedGoogle Scholar
- Shao M, Wang X, Ni M, Bibi N, Yuan S, Malik W, Zhang H, Liu Y, Hua S: Regulation of cotton fiber elongation by xyloglucan endotransglycosylase/hydrolase genes. Genet Mol Res. 2011, 10: 3771-3782. 10.4238/2011.October.27.1.View ArticlePubMedGoogle Scholar
- Lee J, Burns TH, Light G, Sun Y, Fokar M, Kasukabe Y, Fujisawa K, Maekawa Y, Allen RD: Xyloglucan endotransglycosylase/hydrolase genes in cotton and their role in fiber elongation. Planta. 2010, 232 (5): 1191-1205. 10.1007/s00425-010-1246-2.View ArticlePubMedGoogle Scholar
- Michailidis G, Argiriou A, Darzentas N, Tsaftaris A: Analysis of xyloglucan endotransglycosylase/hydrolase (XTH) genes from allotetraploid (Gossypium hirsutum) cotton and its diploid progenitors expressed during fiber elongation. J Plant Physiol. 2009, 166 (4): 403-416. 10.1016/j.jplph.2008.06.013.View ArticlePubMedGoogle Scholar
- Danner H, Boeckler GA, Irmisch S, Yuan JS, Chen F, Gershenzon J, Unsicker SB, Köllner TG: Four terpene synthases produce major compounds of the gypsy moth feeding-induced volatile blend of Populus trichocarpa. Phytochemistry. 2011, 72 (9): 897-908. 10.1016/j.phytochem.2011.03.014.View ArticlePubMedGoogle Scholar
- Chen XY, Chen Y, Heinstein P, Davisson VJ: Cloning, expression, and characterization of (+)-delta-cadinene synthase: a catalyst for cotton phytoalexin biosynthesis. Arch Biochem Biophys. 1995, 324 (2): 255-266. 10.1006/abbi.1995.0038.View ArticlePubMedGoogle Scholar
- Chen XY, Wang M, Chen Y, Davisson VJ, Heinstein P: Cloning and heterologous expression of a second (+)-delta-cadinene synthase from Gossypium arboreum. J Nat Prod. 1996, 59 (10): 944-951. 10.1021/np960344w.View ArticlePubMedGoogle Scholar
- Shadle G, Chen F, Srinivasa Reddy M, Jackson L, Nakashima J, Dixon RA: Down-regulation of hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase in transgenic alfalfa affects lignification, development and forage quality. Phytochemistry. 2007, 68 (11): 1521-1529. 10.1016/j.phytochem.2007.03.022.View ArticlePubMedGoogle Scholar
- ZHAO L, S-z WANG, F-j QI, W-w ZHANG, G-l JIAN: Isolation, Characterization and Phylogenetic Analysis of the Resistance Gene Analogues (RGAs) in Upland Cotton (Gossypium hirsutum L.) Zhongzhimian KV-1 [J]. Cotton Sci. 2009, 6: 007-Google Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-652. 10.1038/nbt.1883.PubMed CentralView ArticlePubMedGoogle Scholar
- Iseli C, Jongeneel CV, Bucher P: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. ISMB. 1999, 1999: 138-148.Google Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (suppl 2): W293-W297.PubMed CentralView ArticlePubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.View ArticlePubMedGoogle Scholar
- Audic S, Claverie J-M: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.PubMedGoogle Scholar
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29 (4): 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M: mapman: a user‒driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37 (6): 914-939. 10.1111/j.1365-313X.2004.02016.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.