Transcript profiling by microarray and marker analysis of the short cotton (Gossypium hirsutum L.) fiber mutant Ligon lintless-1 (Li 1 )
© Gilbert et al.; licensee BioMed Central Ltd. 2013
Received: 27 March 2013
Accepted: 12 June 2013
Published: 17 June 2013
Cotton fiber length is very important to the quality of textiles. Understanding the genetics and physiology of cotton fiber elongation can provide valuable tools to the cotton industry by targeting genes or other molecules responsible for fiber elongation. Ligon Lintless-1 (Li 1 ) is a monogenic mutant in Upland cotton (Gossypium hirsutum) which exhibits an early cessation of fiber elongation resulting in very short fibers (< 6 mm) at maturity. This presents an excellent model system for studying the underlying molecular and cellular processes involved with cotton fiber elongation. Previous reports have characterized Li 1 at early cell wall elongation and during later secondary cell wall synthesis, however there has been very limited analysis of the transition period between these developmental time points.
Physical and morphological measurements of the Li 1 mutant fibers were conducted, including measurement of the cellulose content during development. Affymetrix microarrays were used to analyze transcript profiles at the critical developmental time points of 3 days post anthesis (DPA), the late elongation stage of 12 DPA and the early secondary cell wall synthesis stage of 16 DPA. The results indicated severe disruption to key hormonal and other pathways related to fiber development, especially pertaining to the transition stage from elongation to secondary cell wall synthesis. Gene Ontology enrichment analysis identified several key pathways at the transition stage that exhibited altered regulation. Genes involved in ethylene biosynthesis and primary cell wall rearrangement were affected, and a primary cell wall-related cellulose synthase was transcriptionally repressed. Linkage mapping using a population of 2,553 F2 individuals identified SSR markers associated with the Li 1 genetic locus on chromosome 22. Linkage mapping in combination with utilizing the diploid G. raimondii genome sequences permitted additional analysis of the region containing the Li 1 gene.
The early termination of fiber elongation in the Li 1 mutant is likely controlled by an early upstream regulatory factor resulting in the altered regulation of hundreds of downstream genes. Several elongation-related genes that exhibited altered expression profiles in the Li 1 mutant were identified. Molecular markers closely associated with the Li 1 locus were developed. Results presented here will lay the foundation for further investigation of the genetic and molecular mechanisms of fiber elongation.
KeywordsCotton Ligon Lintless-1 Microarray Cotton fiber elongation Microsatellite markers
Cotton seed fibers are single-celled trichomes that initiate from the ovule epidermal cells on or about the day of anthesis (DOA) . Approximately 25% of the ovule epidermal cells differentiate into fiber cells during the initiation stage of cotton fiber development and subsequently undergo a period of rapid elongation known as the elongation stage [2, 3]. The rate of fiber elongation peaks at approximately 6 to 12 days post-anthesis (DPA) and nears cessation around 22 DPA . During peak elongation fiber cells can increase in length at rates of 2 mm / day or more depending on environmental factors and genotypes [5–7]. Beginning at 12–16 DPA and overlapping with the elongation phase is the secondary cell wall (SCW) biosynthesis stage. During this stage cellulose is synthesized and deposited between the primary cell wall and the plasmalemma [8, 9]. The period of overlap between the elongation stage and the initial stage of SCW biosynthesis is referred as the transition period. Elongation and SCW biosynthesis continue until the fibers reach full length [25–30 mm in Upland cotton (Gossypium hirsutum L.) cultivars] , after which the cotton bolls open and the fibers desiccate under exposure to the environment. The environmental and genetic factors that influence the timing of these processes have been shown to also influence the development of desirable fiber traits such as lint yield and fiber quality [7, 11–13].
Several naturally occurred cotton mutations affecting a range of fiber phenotypes have been genetically and functionally characterized in cotton. Examples include the completely glabrous seeds (lintless and fiberless) observed in MD17 , the fuzzless/lintless (fl) mutant of XZ-142 [15, 16], and lines with seeds containing only lint and no fuzz, such as the naked seed lines N 1 and n 2 . Mutant lines exhibiting very short seed fibers include the Ligon Lintless-1 and −2 lines (Li 1 and Li 2 ) [18, 19]. Recently, Cai et al.  analyzed a man-made mutant Lix that showed similar phenotype to Li 1 . The understanding that initiation, elongation, and secondary cell wall synthesis are distinct developmental processes often leads to the utilization of the applicable mutant to study the specific process of interest. For example, fl mutant seeds lacking any fiber emergence have served as models for studying initiation processes where enrichment of the homeodomain–leucine zipper transcription factor (GhHD1) and GhMyb25 were identified as important for initiation [21, 22]. Likewise, N 1 , with its lack of fuzz fiber and sparsely-distributed lint fibers has been used to characterize fiber elongation processes .
In a near-isogenic state with the cotton line Texas Marker-1 (TM-1), both the Li 1 and Li 2 mutants have seed fibers that are extremely short (< 6 mm) compared to wild type (WT) fibers that are typically greater than 20 mm in length [19, 23, 24]. As a monogenic dominant trait, the short-fiber phenotypes of Li 1 and Li 2 are identical in either a homozygous dominant or heterozygous state. Unlike the Li2 mutant which appears healthy and morphologically identical to the homozygous recessive wild-type plants with the exception of shorter seed fibers, the Li 1 mutant exhibits pleiotropy in the form of severely stunted and deformed plants in both the homozygous dominant and heterozygous state .
Since the seed fibers of Li 1 and Li 2 fibers are shortened lint and fuzz fibers, these cotton mutants represent excellent candidates to study the molecular mechanisms of fiber elongation. Previously, our laboratory conducted extensive analysis of the Li 2 mutant using microarray technology, molecular mapping and metabolomic analysis [25, 26]. We developed microsatellite markers associated with the Li 2 genetic locus, and identified transcripts or genes and metabolites that were affected by the Li 2 mutation. In order to gain more comprehensive knowledge about cotton fiber development, and especially fiber elongation, we included the Li 1 mutant as a subject of our investigation.
The Li 1 mutant has been used as a model to study both primary and secondary cell wall processes [27–30]. However, previous microarray experiments with the Li 1 mutant conducted during either very early elongation or later SCW stage failed to identify significant numbers of differentially expressed transcripts. For example, the microarray experiments conducted by Bolten et al. using 24 DPA fibers only identified ~100 differentially expressed transcripts, notable among them SuSy, Expansins, and Myb transcription factors. However, apparent phenotypic differences in the Li 1 as early as 3 DPA  indicating that altered gene expression may exist at or before this stage. Noting this, a microarray experiment conducted by Liu et al.  analyzed the Li 1 mutant at the initiation and elongation stages of 0, 3 and 6 DPA. Their findings concurred with several earlier studies on the relevance of auxin, gibberellins, brassinosteroid and ethylene-related pathways in fiber development. Elongation stage (6 DPA) fibers from Li 1 demonstrated a significant alteration in transcript profiles, with 1,398 target sequences showing altered expression in the mutant. Despite this, a crucial gap remains in our understanding of how the Li 1 mutation affects the transcript profile at the transition period (later elongation stages and early SCW stages). This paper is the first attempt to analyze gene expression patterns in the Li 1 mutant using microarray technology at these critical developmental stages. Here we provide a more complete picture of the molecular events directly controlling fiber elongation. Further, it will better define the mutation in terms of its effects on primary cell wall elongation and early secondary cell wall synthesis.
Prior research has determined that the Li 1 gene is located on chromosome 22 using both SSR  and RFLP markers . Karaca et al.  identified the SSR marker MP4030 that was 12.83 cM away from the Li 1 locus. Rong et al.  provided the highest resolution to date, determining Li 1 was flanked by RFLP markers Gate4CA09 and Coau1J04 at 2.7 and 1.3 cM away, respectively, based on 151 F2 progeny derived from an interspecific cross of G. barbadense Pima S-7 × Li 1 mutant.
In order to conduct a comprehensive study of the Li 1 mutant, we first created near-isogenic lines (NIL) in DP5690 genetic background by implementing an extensive backcross scheme. The use of DP5690, a modern variety, exhibits stronger growth characteristics than TM-1 in many climates and permits for additional analysis with the previously characterized Ligon Lintless-2 (Li 2 ), which is also in the DP5690 background . Using the two NILs as parents, we made a very large F2 population comprising 2,553 progeny which was used to identify molecular markers closely associated with the Li 1 locus. To understand the molecular events that control fiber elongation and identify regulatory elements involved in this process, we obtained transcript profiles at 3 DPA (beginning of elongation), 12 DPA (late elongation stage) and 16 DPA (early SCW) using Affymetrix microarrays, and analyzed morphological characteristics of fibers from the Li 1 NILs at different developmental stages. The objectives of this research were to determine the location of the Li 1 locus on the chromosome, and identify genes that were differentially expressed during the development of WT and mutant Li 1 fibers. This information will, in turn, be useful to identify the Li 1 gene, and help to elucidate the molecular mechanisms of this gene on fiber elongation.
Plant materials for microarray and quantitative PCR (qPCR) experiments
For the microarray and qPCR experiments, a total of 100 homozygous (Li 1 Li 1 ) Li 1 mutant and 100 WT (li 1 li 1 ) plants were planted in a field at the Southern Regional Research Center, New Orleans, LA in the summers of 2011 and 2012. The soil type in New Orleans was Aquent dredged over alluvium in an elevated location to provide adequate drainage. Flowers were tagged and sample collections were made before 10:00 am and immediately placed on ice. All samples collected within each developmental stage were tagged and collected on the same day. Fruits were randomly grouped into 3 individual replicates with 5–10 fruits per replicate. Fruits were then dissected, frozen in liquid nitrogen and stored at −80°C until further processing.
A WT DP5690 (li 1 li 1 ) was used as the female in a cross with its near-isogenic mutant line (Li 1 Li 1 homozygous plant). Two thousand five hundred fifty-three F2 plants derived from approximately 20 F1 plants were planted along with their parents in a field in Stoneville, MS in 2012. The Li 1 trait of each F2 progeny plant was evaluated after boll maturation and opening (about 60 DPA). Standard conventional field practices were applied during the growing season. The soil type in Stoneville, MS was Bosket very fine sandy loam.
Fiber length and cellulose content measurements
Fiber length was measured using the method described by Schubert et al. . Two replicate samples with 10 ovules each were measured. For mass determination, air-dried fibers were gently removed from all ovules of each sample, and weighted on an analytical balance. Cellulose content from each fiber sample was measured using the method described by Updegraff  with minor modification. Dried fiber samples were cut into small pieces. Ten mg of the blended fibers were placed in 5 mL Reacti-Vials™ (Thermo Fisher Scientific, Waltham, MA). Non-cellulosic materials in fibers were hydrolyzed with acetic-nitric reagent. The remaining cellulose was hydrolyzed with sulfuric acid and measured by a colorimetric assay with anthrone using Avicel PH-101 (FMC, Rockland, ME) as a cellulose standard. The average cellulose content for each fiber sample was obtained from two biological and three technical replications.
Imaging analysis of fiber cross-sections
Fibers from 28 and 40 DPA were manually separated from the seed. After bundling the fibers together, a new razor blade was used to cut 2.5 – 3 mm of fiber from the end to be sectioned. This ensured a majority of fibers were cross-sectioned near the middle of the fiber. The fiber samples were embedded, thin-section cut, and photographed using the method previously described . The images were taken using a Nikon Cambridge Quantimet 900 microscope at 200× magnification with a Hitachi KP-DSO camera.
RNA isolation from cotton fibers
RNA was isolated as previously described . Briefly, material was obtained from developing ovules using a glass bead shearing technique . To separate the fibers from the ovules the samples were shaken vigorously enough to break fibers without damaging the ovules. Isolation of RNA was conducted using the Sigma Spectrum™ Plant Total RNA Kit (Sigma-Aldrich, St. Louis, MO) with on-column DNaseI digestion according to the manufacturer’s instructions. RNA quantity was determined by using a Nanodrop 2000 spectrophotometer (NanoDrop Technologies Inc., Wilmington, DE). A RNA integrity number (RIN) was determined for each sample using an Agilent Bioanalyzer 2100 and the RNA 6000 Nano Kit Chip (Agilent Technologies Inc., Santa Clara, CA). Only samples with RIN values of 7.0 or higher were used for further analysis.
Microarray hybridizations and data analysis
The minimum information about microarray experiments (MIAME) guidelines were followed for all microarray experiments conducted in this study . The microarray chips used for this study were the commercially available Affymetrix GeneChip® Cotton Genome Microarray (Affymetrix Inc., Santa Clara, CA), comprising 239,777 probe sets representing 21,854 cotton transcripts from a variety of EST databases. Labeling of the RNA was conducted using the Affymetrix GeneChip® 3′ IVT Express Kit and hybridizations were conducted according to the manufacturer’s protocols. Hybridizations were conducted on 3, 12, and 16 DPA samples with two biological replicates from each developmental stage. Data normalization and the determination of statistically relevant deviations in expression patterns were performed as described . To assist in analysis of biological processes represented in the data, Gene-Ontology Enrichment Analysis (GOEA) was performed using the agriGO Singular Enrichment Analysis tool . The statistical test method used was the Fisher’s Exact test (significance level 0.05). Annotation of the probe sets was accomplished with Blast2Go , and analysis of the cellulose synthase probes was conducted by translating all 6 reading frames of the probe sets and subjecting them to blastp analysis.
To investigate the activity of known cell elongation and cell wall-related genes, the microarray probe sets were compared with published lists of genes [27, 40]. Each probe set reference sequence was aligned to the G. raimondii reference genome  with blastn and the best hit with e value of 1×10-30 or smaller was used to establish an annotation. The G. raimondii gene annotations specify an arabidopsis homolog for each gene, which was used to classify the functions of each Affymetrix probe set. Seventy-four probe sets were classified as elongation genes based on an earlier report . To classify probe sets as primary or secondary cell wall genes, co-expression with arabidopsis microarray data with known primary or secondary cell wall cellulose synthase genes  was obtained from ATTEDII . This strategy produced lists of 81 primary cell wall and 43 secondary cell wall gene probe sets.
Reverse transcription and quantitative PCR
The cDNA reactions were performed using the iScript™ cDNA Synthesis kit (Bio-Rad Laboratories, Hercules, CA) per the manufacturer’s instructions. The reaction without reverse transcriptase served as negative control for testing genomic DNA contamination of the RNA samples. This reaction was then used as template in a qPCR reaction to verify that no amplification occurred. After cDNA synthesis, the qPCR reaction was conducted using iTaq™ SYBR® Green Supermix (Bio-Rad Laboratories) in a Bio-Rad CFX96 real time PCR detection system. PCR conditions and the protocols for determining primer efficiencies were as previously described . Ubiquitin-conjugated protein (UCP) (Genebank AI730710) was used as the endogenous reference gene. Primer sequences are listed in Additional file 2.
SSR marker analysis and genetic mapping
Young leaves were collected from each individual F2 plant and parents, and stored at −80°C. Total DNA was extracted from frozen leaves according to Fang et al. . The Li 1 gene was previously determined to reside on chromosome 22 (Chr. 22) [31, 32]. To rapidly identify molecular markers closely linked to the Li 1 locus, all simple sequence repeat (SSR) markers that were mapped on both Chr.22 and its homeologous Chr. 4 based on the high density consensus genetic map  were selected for analysis. RFLP markers reported by Rong et al.  were not screened due to unavailability of probes, and technical difficulties for RFLP marker analysis. The SSR markers mapped on Chr. 5 and its homeologous Chr.19 were also included because of a known translocation between Chr. 4 and Chr. 5 [45, 46]. All together, a total of 921 SSR markers mapped in these four chromosomes were screened for polymorphism between DNA bulks. For the WT bulk, DNA of 10 F2 plants with WT phenotype were pooled at equal ratio and diluted to 50ng/μL. The mutant type bulk consisted of DNA from 10 F2 plants with short seed fiber. Four DNA bulks were made, two for each type. The polymorphic markers were then analyzed using 96 F2 progeny to identify markers closely associated with the Li 1 locus. Only SSR markers that revealed less than 10 recombination events were analyzed among the total 2,553 F2 progeny plants.The PCR amplification conditions and marker data acquisition were according to Fang et al. . All SSR primer sequences can be obtained from Cotton Marker database (http://www.cottonmarker.org). Segregation data for the Li 1 trait and SSR markers were mapped using program JoinMap4.0  with logarithm of odds score =25.
Functional analysis of the Li 1 region
To obtain sequences in the region determined to contain the Li 1 mutant gene, the identified SSR marker sequences were aligned to the diploid G. raimondii D5 genome . We also blasted other SSR marker sequences from the Li 1 interval of Chr. 22 based on the high density consensus map . This permitted comparison of the gene annotations for this interval with the annotations provided by Affymetrix to identify candidate genes. qPCR on select genes was conducted as described above.
Sequencing introns to develop additional polymorphic markers
To identify additional polymorphisms that could facilitate finer mapping of the Li 1 locus, we sequenced introns from the annotated genes in the 17 Mb interval between the markers TMB2500 and DPL0489 according to the reference G. raimondii genome sequences. Our rational to sequence introns was that introns tend to have higher sequence variations than exons. We designed primers that flanked 23 introns but were anchored in protein coding sequences. Three genes that resided in this interval and showed differential expression between WT and mutant were included. Amplicons were generated from total genomic DNAs of the parental lines and two F2 individuals (one WT type and one mutant type) and were Sanger sequenced.
Fiber structure analysis
Early literature indicated that the Li 1 mutant in a TM-1 genetic background demonstrated a reduced rate of crystalline cellulose deposition during primary cell wall synthesis and an increased rate during SCW synthesis, resulting in a “thickened” appearance of the cell wall in the mutant [18, 30]. However, it remains unclear if this effect is due primarily to the mutation causing inhibition of fiber elongation processes or due to the mutation affecting SCW synthesis processes such as cellulose deposition, or both. To better characterize phenotypic changes in the mutant during late elongation and secondary cell wall deposition, changes in fiber length, mass, density, and cellulose content over developmental times were measured. Grown in standard field conditions in New Orleans, LA, the homozygous Li 1 mutant demonstrated a characteristic short fiber phenotype and other previously characterized pleiotropic characteristics (Figure 1A and Additional file 3). Heterozygous individuals derived from DP5690 x Li 1 Li 1 crosses demonstrated a dosage effect resulting in an intermediate plant size (Additional file 3). Measurements of fiber physical properties indicated that Li 1 mutants had a dramatic difference in both length and fiber dry weight (mg/seed) at the developmental stages measured (Figure 1B) as compared to the WT. Elongation ceased or remained static in the Li 1 mutant by 12 DPA, however fiber dry weight continued to increase during SCW synthesis and through to maturity, likely due to the continued deposition of cellulose. Calculating fiber density further illustrated the continued increase in biomass observed over the developmental stages and a much higher density in comparison to the WT. However, the overall cellulose content per unit mass was similar between WT and the Li 1 mutant in 12, 16, 28 DPA and mature fibers. Image analysis indicated the Li 1 mutant fibers at 28 DPA, (the youngest stage that was technically obtainable), were in general “thicker” than WT fibers (Figure 1C), which corroborated our observations on the increasing mass of the fibers. Combined, these results suggested that while thickening of the secondary cell wall occurred through developmental stages, it was not due to a relative increase in overall cellulose content or rate of production per fiber cell in the Li 1 mutant. Rather, due to early cessation of fiber elongation, the stunted fiber did not distribute the cellulose along its longitudinal axis, thus resulting in a thickened appearance, and an increased mass per unit length.
Gene ontology analysis results of target sequences that are up or down regulated in the Li 1 mutant at 12 DPA or 16 DPA
Gene ontology categorization*
Representative annotations of probed sequences (Li1/Wt expression ratio)
Unique to 12 DPA
Response to hormone stimulus GO:0009725 (0.00023)
eg. indole-3-acetic acid-amido synthetase (0.48), xyloglucan endotransglucosylase (2.2), shaggy-related kinase (2.1), auxin-responsive protein (.34) and auxin response (0.42).
Response to cytokinin stimulus GO:0009735 (0.00014)
eg. ap2 erf domain-containing transcription factor (2.1), homeobox protein knotted-1-like 3-like (2.1)
Transmembrane receptor protein tyr kinase signaling pathway GO:0007169 (1.9e-06)
e.g. strubbelig-receptor family 7 protein (0.45), lrr receptor-like serine threonine-protein kinase (0.35),
Protein amino acid phosphorylation GO:0006468 (1.5e-05)
e.g. cyclin-dependent kinase f-4-like (0.48), shaggy-related protein kinase eta-like (2.07), serine threonine-protein kinase aurora-1 (.43), feronia receptor-like kinase (0.43)
Hydrolase activity, hydrolyzing O-glycosyl compounds GO:0004553 (3e-07)
e.g. cobra-like 4 protein (3.74), beta galactosidase 1 (0.46), xyloglucan endotransglucosylase hydrolase (0.44), acidic chitinase (3.32), xyloglucan endotransglucosylase hydrolase protein 2 (2.33), beta-xylosidase alpha-l-arabinofuranosidase 2-like (0.43)
Xyloglucan:xyloglucosyl transferase activity GO:0016762 (0.00016)
e.g. xyloglucan endotransglucosylase hydrolase (0.41), xyloglucan endotransglucosylase hydrolase protein 2 (2.33), probable xyloglucan endotransglucosylase hydrolase protein 32-like (2.56)
Protein kinase activity GO:0004672 (0.00075)
serine threonine-protein kinase (3.12), probable receptor-like serine threonine-protein kinase at5g57670-like (11.02)
Unique to 16 DPA
Cellular nitrogen compound metabolic process GO:0034641 (1.1e-05)
glutamine synthetase (2.05), phenylalanine ammonia-lyase (0.26), asparagine synthetase (2.06), serine threonine protein kinase 2 (3.23)
A significant number of processes that were common to 12 DPA and 16 DPA, but not to 3 DPA showed altered expression in the mutant (Additional file 5). GOEA analysis of this category of genes indicated a large decrease in probe sets categorized in nucleosome assembly (GO:0006334) and lipid transport (GO:0006869). There was an increase in mitochondrial electron transport (GO:0006120), which includes NADH-dehydrogenase genes and NADH-plastoquinon oxidoreductase subunits, both members of Complex I of the electron transport chain, known to be high producers of reactive oxygen species (ROS) .
Select elongation-related probe sets showing altered regulation in the Li 1 mutant*
Genes up-regulated in theLi 1 mutant relative to its near isogenic wild type line
GhACO3, ethylene biosynthesis, elongation related
1-aminocyclopropane-1-carboxylate oxidase homolog 4-like
no homology known
brassinosteroid-regulated protein bru1
homology to xyloglucan endotransglcosylase (AT4G14130.1), expression correlates with elongation inhibition
fasciclin-like arabinogalactan protein
EST from fiber, homology to G. hirs. fasciclin-like arabinogalactan protein 3 (FLA3). AtFLA3 overexpression leads to defective elongation in stamen.
At4g37260, MYB73. Highly responsive to ethylene, ABA. Downregulated in fl mutant elonagtion phase.
ethylene-responsive transcription factor wri1
wri1-1 mutants demonstrate defective elongation of hypocotyl.
sucrose synthase sus1
GhSus1 isoform C, targeted to cell wall during secondary cell wall synthesis.
xyloglucan endotransglucosylase hydrolase
TCH4, dwarf A. thaliana mutants have reduced expression, elongation related.
pollen ole e 1 allergen and extensin family protein
Genes down-regulated in the Li 1 mutant relative to its near isogenic wild type line
ACO4, ethylene biosynthesis, induced in ovule culture by fatty acid synthesis, increasing ethylene production and root elongation.
EST, significant homology to ACO proteins, but not identical to ACO1-4
Highest homology with GhACS 3/6
gibberellin synthesis, stimulated by Auxin treatment in A. thaliana, which also causes hypocotyl cell elongation. Increased expression in transgenic cotton causes increased fiber length.
myb-like transcription factor 3
GhMYB3, contains gibberellin responsive GLABROUS1, which promotes trichome formation in A. thaliana.
xyloglucan endotransglucosylase hydrolase
GhXTH2, high homology to coding region of GhXTH1 which produces longer fibers when over-expressed in transgenic cotton.
rho gtpase activation protein
High homology to ROP1, tip-localized GTPase responsible for cell elongation and polarity.
Twelve probe sets representing an unknown number of AP2/ERF (ethylene response factor) domain-containing transcription factors were significantly differentially expressed, with the majority demonstrating down-regulation (e.g. Ghi.7874) (Table 2 and Additional file 6). Eight probe sets measuring levels of the pathogenesis-related 10 (PR10) proteins, a family of defense and stress-related genes regulated by jasmonic acid, ethylene and other effectors  exhibited expression patterns 8 to 96 fold higher in the Li 1 mutant. Sucrose Synthase 1 (Sus1) was highly up-regulated (56-fold) at 3 DPA in the Li 1 mutant, then not differentially expressed at 12 DPA and 16 DPA.
Transcription analysis of cellulose synthase activity
Corroboration of microarray data
Microarray expression ratios of Li 1 /WT and corroboration by RT-qPCR analysis*
cellulose synthase catalytic subunit
cellulose synthase-like protein e6-like
pathogenesis-related protein 10
Tubulin alpha 4 (Tua 4)
Mapping the Li 1 locus region with SSR markers
A total of 13 kb sequences from the introns of 23 genes residing in the 17 Mb interval between markers TMB2500 and DPL0489 were compared between the WT and Li 1 mutant. However, we were unable to identify additional sequence polymorphisms to facilitate finer genetic mapping of the Li 1 locus.
Functional analysis of the Li 1 region
A previous microarray experiment examining 0 and 3 DPA fibers from the Li 1 mutant found little evidence of global or significant alterations in gene expression patterns . Conversely, microarray data of 6 DPA fibers and proteomic analysis of 12 DPA from Li 1 fibers demonstrated significant perturbation of expression profiles in the mutant, indicating that processes related to primary cell wall elongation are affected in the Li 1 mutant [27, 29]. In addition, earlier studies of the Li 1 mutant focusing on secondary cell wall synthesis failed to report findings of upstream processes being significantly affected . Thus this analysis on late elongation and early SCW stages would complement and extend these earlier findings aimed to better characterize the effects of the mutation at both primary and secondary cell wall synthesis, and provide data for the fiber development during the transition period. The morphological and molecular data presented here supports the model that the Li 1 mutation is an upstream factor primarily targeting elongation processes. Several ontological categories of genes and individual genes that have previously been identified as having a role in fiber elongation were identified, as were new targets for investigation.
Previous proteomic analysis of Li 1 at 12 DPA  identified a limited number of proteins that could be corroborated by our microarray data. A comprehensive analysis revealed that both data sets include the down regulation of cytoskelatal tubulin proteins and metabolism-related proteins (ie. glycolipid transfer protein and a pyrophosphatase). The stress response proteins nucleoredoxin and germin-like proteins, a flavanone-3-hydroxylase and translation factor 6 were increased in both data sets. However, the remainder of proteins identified as showing altered regulation were not corroborated by our data. A detail list is shown in the Additional file 8. This could be due to the difference in genetic background, or due to the technical limitations of 2-dimensional gel electrophoresis. Analysis of more recent Affymetrix microarray data obtained at 6 DPA  is limited by the fact that the authors did not provide probe set information, thus preventing a more rigorous or comprehensive analysis from being conducted. However, they also identified several actin and tubulin genes that were down-regulated in the mutant. At 6 DPA there was an increase in xyloglucan endotransglucosylase, fascilin-like arabinogalactan protein, and the ethylene synthesis-related 1-aminocyclopropane −1-carboxylate oxidase (GhACO3) which our data indicated remained true at 12 DPA. As expected several differences also exist. Analysis of 6 DPA fibers showed decrease expression levels of auxin-IAA-related genes and the ethylene responsive AP2/ERF family genes, which then became up-regulated in the mutant at 12 DPA. Conversely, another gene involved in ethylene production, GhACO4, was up-regulated at 6 DPA but down-regulated in the mutant at 12 DPA. Genes that were down-regulated in Li 1 at 6 DPA, but were at WT levels in 12 DPA were enriched in plasma membrane associated proteins (GO:0005886, p = 0.0008), such as filament-like plant protein 4, myosin heavy chain, and perk1-like protein kinase.
The use of the Li 1 mutant as a model system for both elongation and secondary cell wall synthesis was based on earlier studies utilizing Li 1 in a TM-1 genetic background that measured the ratio of fiber weight to length and showed an increase in fiber mass throughout development, implicating continued or increased cellulose deposition . Despite its very short fibers, the authors claimed the dry mass of the Li 1 fiber was approximately 83% of its NIL. A second study using the same genetic line that measured [14C] glucose deposition in the secondary cell wall indicated a 5-fold higher rate of cellulose deposition per mm of fiber in the Li 1 mutant . However, much of our data using Li 1 in a DP5690 genetic background suggested that molecular events involved in cellulose deposition and secondary cell wall synthesis were not affected to the degree that primary cell wall processes were. Measurements indicated that the dry mass of the fiber was approximately 20% of its NIL and that the actual cellulose content per unit mass remained unchanged (Figure 1B). Further, qPCR of cellulose synthase genes showed that secondary cell wall genes remained largely unaffected. These results support the model that thickening of the cell wall is due to inhibited elongation processes, but not due to increased cellulose production. The “increased rate” of glucose incorporation observed by Kohel et al.  was based on measurements of cellulose production per unit length instead of per unit mass. However, since elongation ceases in the Li 1 mutant but its cellulose production continues at WT levels, we speculate that Li 1 fibers expand outwardly, thus resulting in an increased cellulose production per unit length. It is worthy to mention that the use of a different NIL in our study as compared with previous studies may also account for differences in fiber measurements.
The ontological categories of the 434 unique probes differentially expressed at 12 DPA included Xyloglucan:xyloglucosyl transferase activity (GO:0016762) (Table 1). Xyloglucan endotransglycosylase/hydrolase (XTH) enzymes are proposed to disrupt xyloglucan-cellulose crosslinks, thus permitting cell wall rearrangement and fiber elongation . Six probe sets identified as XTH showed differential regulation in the Li 1 mutant. One probe set (Ghi.7911.1) was down-regulated at 3 DPA and significantly up-regulated (9 fold) at 12 DPA and 16 DPA (Additional file 6). This gene [GenBank:AY476737] remains uncharacterized in Gossypium, however its homeologous gene in A. thaliana TCH4/XTH22 (TAIR;AT5G57560) demonstrated very similar results in a study conducted on petiole elongation that indicated XTH22 is under different genetic regulation than other characterized XTH’s . Of the remaining XTHs identified in our microarray data, three probe sets exhibited decreased expression at 12 DPA and one slightly increased. While it is clearly plausible that a decreased expression of XTHs could lead to the inhibition of elongation, the role of XTHs that increased in the Li 1 mutant warrants further investigation. Other ontological categories affected in the mutant included “responses to hormone synthesis” (GO:0009725) (ie. IAA synthetase and auxin-response genes) and “hydrolase activity, hydrolyzing O-glycosyl compounds” (GO:0004553), which contain 19 probe sets likely correlated with hormonal activation to signal entry into the primary to secondary cell wall transition stage . The role of IAA, auxin response genes, and hormonal regulatory enzymes such as glucosyltranferases in fiber development processes has been widely documented . It is of interest to note here that the only category identifying this set of related genes was the list of probes that were up or down-regulated only at 12 DPA. Probe sets that were also enriched or decreased in 16 DPA mutant tissues (Additional file 5) were not enriched in these ontological categories, implying that altered hormonal regulation in the mutant was occurring primarily before and at 12 DPA. Analysis of hormone contents by Chen et al.  identified altered levels of Abscisic Acid, dihydrozeatin, and others in the Li 1 mutant, but only analyzed up to 8 DPA. Microarray data on 6 DPA fibers reported altered expression of several hormone related peptides in the Li 1 mutant, such as auxin-related genes, gibberellins, brassinosteroid, abscisic acid and jasmonic acid-related genes . These data suggest, when taken together with molecular and morphological data, that the early cessation of elongation in Li 1 , is in part due to the culmination of altered hormonal factors that occurs during early to mid elongation.
Consistent with what has been reported at earlier developmental stages, multiple genes related to ethylene biosynthesis and ethylene response were differentially expressed in the Li 1 fibers . A key enzyme in ethylene biosynthesis, 1-aminocyclopropane-1-carboxylate synthase (ACS), hybridized by probe set Ghi.5451and demonstrated a 2-fold down-regulation at 12 DPA (Table 2). Previous studies have correlated increased ACS activity with cotton fiber elongation . 1-aminocyclopropane-1-carboxylate oxidases (ACO) 1–4 are also involved in ethylene biosynthesis . ACO3 (Ghi.798.1) showed a significant up-regulation at all stages analyzed, and ACO4 (Ghi.8025.1) exhibited down-regulation at 12 DPA and 16 DPA. ACO1 and ACO2 (Ghi.6953.1 and Gra.2141.1, respectively) did not show altered expression in the Li 1 mutant. Consistent with our data, previous studies have implicated ACO3 as having peak expression at the late elongation stage near 12 DPA. ACO1 and 2 have demonstrated different expression patterns with peak expression in early secondary cell wall synthesis stage [54, 55], implying that they may be involved in later developmental stages. The fact that ACO3 was related to elongation in previous studies and showed altered regulation in our microarray (in contrast to ACO1 and 2) further supports Li 1 being a key elongation-related mutation.
Pathogenesis-related protein 10 family of proteins (PR10) exhibited a significant alteration in their expression patterns in the Li 1 mutant, ranging from 8 to 96 fold difference between Li 1 and WT. This family of proteins consists of a large and functionally diverse group of proteins, ranging in function from antimicrobial/antiviral activity, hormone/ligand binding, secondary metabolism, and abiotic stress response (reviewed in ). Further implying a role in growth and development, individual PR10 proteins have been found to be regulated by multiple phytohormone-related cis-regulatory sequences including ethylene response elements  and brassinosteroids . However it remains unclear if the stress response seen in the Li 1 is related to previously mentioned ROS levels or a response to altered hormone expression.
The cellulose synthase activity reported here is of particular interest. Previous reports have speculated that an increase in secondary cell wall cellulose synthase activity may account for the thickened cell wall of the Li 1 , however, our data failed to confirm this. Rather, another GhCes, probed by Ghi.58712, demonstrated altered expression at the elongation stage of development and showed an expression pattern consistent with elongation-related activity in the WT. Translation of the probe target’s consensus sequence and TAIR blast search showed the Ghi.58712 target sequence shared highest homology, although not identity with AtCesA9,-2,-5, -6 and −3 (E = 0.0). AtcesA9 is only expressed during embryogenesis, however the remaining AtCesA2, -5, -6, and −3 are members of a primary cell wall associated cellulose synthase complex (CSC) [72, 80]. Substantial evidence exists that these closely related genes are cell elongation related. Atces2 null mutations showed a severe dwarf phenotype in A. thaliana, and functional studies have demonstrated that AtCesA2 and −6 were partially functionally redundant during elongation . A. thaliana mutants at the PROCUST1 locus, which encodes cesA6, exhibited cell elongation effects in a pleiotropic manner . cesA6 promoter-GUS fusion experiments have demonstrated that its expression occurs throughout the hypocotyl and root, peaking in the cell elongation zone of the expanding root .
The Li 1 genetic locus was previously identified as residing on chromosome 22 [31, 32]. Our results confirmed this chromosomal assignment. Based on the high density consensus genetic map constructed by Blenda et al. , it could be seen that the Li 1 locus might be close to the centromeric region as original indicated by Rong et al. . It is worthy of mention that the genetic distances observed in our experiment were larger than those reported by Rong et al. . This greater recombination observed in our research might be due to our much larger (151 vs 2,553) population size and different population structure (interspecific vs intraspecific cross). There is a large gap (9.6 cM) between the Li 1 locus and marker DPL0489 (Figure 4). We have screened all the SSR markers mapped in this interval based on the high density consensus map  and a newly published map , and could not further close the gap. Additionally, we were unable to identify sequence polymorphisms among the 23 gene introns that were located in the region harboring the Li 1 locus even though some of these genes demonstrated differential expression between WT and mutant. This also indicates that altered gene expression may not be necessary due to gene sequence change. This result may also imply that the genomic region harboring Li 1 locus is highly monomorphic. Recently, Cai et al.  reported a similar phenomenon when mapping the Lix locus. One of their flanking markers, NAU3469, was 24.5 cM away from the Lix locus.
The mapping data coupled with the recently released G. raimondii sequences provided the opportunity for additional analysis of sequences in the vicinity of the Li 1 locus. Probe sets from the microarray with high homology to a glycosyl hydrolase family 38, a xyloglucan endotrans-glucosylase hydrolase 32, and AtTCP-20 are near the Li 1 locus, demonstrated differential expression in the Li 1 mutant and have apparent associations with elongation processes. In addition to the already discussed XTH enzymes, glycosyl hydrolase family 38, a family of related mannosidases, affect cell wall phenotypes when mutated in A. thaliana. Interference with AtTCP-20 in planta in A. thaliana by fusion with a repressor domain resulted in severe developmental phenotypes characterized by reduced cellular elongation . Additionally, a transcription factor identified in Gossypium barbadense with a highly homologous TCP domain, GbTCP, produced a short fiber phenotype when silenced by RNAi .
Determining which, if any of these is the Li 1 mutation is currently under investigation. Next-generation sequencing of transcripts by RNA-seq is currently under way in our laboratory, and may reveal SNPs or splice variants responsible for either altered regulation of an elongation-specific gene, or a mutation resulting in a nonfunctional protein. Identification of these variations will facilitate developing closer markers, and help to eventual cloning of the Li 1 gene.
The identification of the gene responsible for the Li 1 phenotype would provide an invaluable tool in the quest to understand fiber elongation processes. In the meantime, the data generated here, in combination with published data from other developmental time points, has the potential to provide a sound basis for the examination of key hormonal, structural, and other pathways involved in cotton fiber elongation.
Measurements of fiber characteristics and microarray analysis of the Li 1 mutant and its WT were conducted at 3, 12, and 16 DPA with the goal of enhancing our understanding of cotton fiber elongation. Both methodologies supported the notion that the early cessation of elongation in Li 1 was due to disruption of primary cell wall elongation-related processes. Further, we identified and discussed several elongation-related genes that exhibited altered expression profiles in the Li 1 mutant, including a putative primary cell-wall related cellulose synthase. We conducted SSR marker analysis on a large population, and using the G. raimondii reference sequence identified elongation-related genes near the Li 1 locus with altered expression levels. The data here will contribute to developing a comprehensive understanding of cotton fiber elongation.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its Additional files.
Cellulose synthase complex
Day of anthesis
Expressed sequence tag
Gene ontology enrichment analysis
Primary cell wall
Restriction fragment length polymorphism
Reverse transcription quantitative polymerase chain reaction
Secondary cell wall
Single nucleotide polymorphism
Simple sequence repeat
Texas Marker 1
This research was funded by United States Department of Agriculture-Agricultural Research Service CRIS project 6435-21000-016-00D. We greatly thank Mr. Chris Florane and Mrs. Tracy Condon for their assistance with field work, sample collection, and fiber structure analysis. We thank Mrs. Holly King for her expertise in fiber sectioning and imaging. Our great appreciation goes to Mrs. Sheron Simpson and Dr. Brian Scheffler at Genomics and Bioinformatics Research Unit at Stoneville, MS for their excellent support in SSR marker analysis. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U. S. Department of Agriculture that is an equal opportunity provider and employer.
- Lee JJ, Woodward AW, Chen ZJ: Gene expression changes and early events in cotton fibre development. Ann Bot. 2007, 100 (7): 1391-1401.PubMed CentralView ArticlePubMedGoogle Scholar
- Basra AS, Malik CP: Development of the cotton fiber. Int Rev Cytol. 1984, 89: 65-113.View ArticleGoogle Scholar
- Tiwari SC, Wilkins TA: Cotton (Gossypium hirsutum) seed trichomes expand via diffuse growing mechanism. Can J Bot. 1995, 73 (5): 746-757.View ArticleGoogle Scholar
- Meinert MC, Delmer DP: Changes in biochemical composition of the cell wall of the cotton fiber during development. Plant Physiol. 1977, 59 (6): 1088-1097.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim HJ, Triplett BA: Cotton fiber growth in planta and in vitro. Models for plant cell elongation and cell wall biogenesis. Plant Physiol. 2001, 127 (4): 1361-1366.PubMed CentralView ArticlePubMedGoogle Scholar
- Schubert AM, Benedict CR, Berlin JD, Kohel RJ: Cotton fiber development - kinetics of cell elongation and secondary wall thickening. Crop Sci. 1973, 13 (6): 704-709.View ArticleGoogle Scholar
- Hinchliffe DJ, Meredith WR, Delhom CD, Thibodeaux DP, Fang DD: Elevated growing degree days influence transition stage timing during cotton fiber development resulting in increased fiber-bundle strength. Crop Sci. 2011, 51: 1683-1692.View ArticleGoogle Scholar
- Willison JH, Brown RM: An examination of the developing cotton fiber: wall and plasmalemma. Protoplasma. 1977, 92 (1–2): 21-42.View ArticleGoogle Scholar
- Seagull RW: Changes in microtubule organization and wall microfibril orientation during in vitro cotton fiber development: an immunofluorescent study. Can J Bot. 1986, 64 (7): 1373-1381.View ArticleGoogle Scholar
- Ruan YL: Recent advances in understanding cotton fibre and seed development. Seed Science Research. 2005, 15 (4): 269-280.View ArticleGoogle Scholar
- Liakatas A, Roussopoulos D, Whittington WJ: Controlled-temperature effects on cotton yield and fibre properties. J Agric Sci. 1998, 130 (4): 463-471.View ArticleGoogle Scholar
- Roussopoulos D, Liakatas A, Whittington WJ: Controlled-temperature effects on cotton growth and development. J Agric Sci. 1998, 130 (4): 451-462.View ArticleGoogle Scholar
- Davidonis GH, Johnson AS, Landivar JA, Fernandez CJ: Cotton fiber quality is related to boll location and planting date. Agron J. 2004, 96 (1): 42-47.View ArticleGoogle Scholar
- Turley RB, Kloth RH: Identification of a third fuzzless seed locus in upland cotton (Gossypium hirsutum L.). J Hered. 2002, 93 (5): 359-364.View ArticlePubMedGoogle Scholar
- Ji SJ, Lu YC, Feng JX, Wei G, Li J, Shi YH, Fu Q, Liu D, Luo JC, Zhu YX: Isolation and analyses of genes preferentially expressed during early cotton fiber development by subtractive PCR and cDNA array. Nucleic Acids Res. 2003, 31 (10): 2534-2543.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang QQ, Liu F, Chen XS, Ma XJ, Zeng HQ, Yang ZM: Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant. Genomics. 2010, 96 (6): 369-376.View ArticlePubMedGoogle Scholar
- Lee JJ, Hassan OS, Gao W, Wei NE, Kohel RJ, Chen XY, Payton P, Sze SH, Stelly DM, Chen ZJ: Developmental and gene expression analyses of a cotton naked seed mutant. Planta. 2006, 223 (3): 418-432.View ArticlePubMedGoogle Scholar
- Kohel R, Quisenberry JE, Benedict CR: Fiber elongation and dry weight changes in mutant lines of cotton. Crop Sci. 1974, 14: 471-475.View ArticleGoogle Scholar
- Kohel RJ, Narbuth EV, Benedict CR: Fiber development of Ligon lintless-2 mutant of cotton. Crop Sci. 1992, 32 (3): 733-735.View ArticleGoogle Scholar
- Cai C, Tong X, Liu F, Lv F, Wang H, Zhang T, Guo W: Discovery and identification of a novel Ligon lintless-like mutant (Lix) similar to the Ligon lintless (Li1) in allotetraploid cotton. Theor Appl Genet. 2013, 126 (4): 963-970.View ArticlePubMedGoogle Scholar
- Wu Y, Machado AC, White RG, Llewellyn DJ, Dennis ES: Expression profiling identifies genes expressed early during lint fibre initiation in cotton. Plant and Cell Physiology. 2006, 47 (1): 107-127.View ArticlePubMedGoogle Scholar
- Wu Y, Llewellyn DJ, White R, Ruggiero K, Al-Ghazi Y, Dennis ES: Laser capture microdissection and cDNA microarrays used to generate gene expression profiles of the rapidly expanding fibre initial cells on the surface of cotton ovules. Planta. 2007, 226 (6): 1475-1490.View ArticlePubMedGoogle Scholar
- Kohel RJ: Linkage tests in Upland cotton, Gossypium hirsutum L. II. Crop Science. 1972, 12: 66-69.View ArticleGoogle Scholar
- Narbuth EV, Kohel RJ: Inheritance and linkage analysis of a new fiber mutant in cotton. J Hered. 1990, 81: 131-133.Google Scholar
- Hinchliffe DJ, Turley RB, Naoumkina M, Kim HJ, Tang Y, Yeater KM, Li P, Fang DD: A combined functional and structural genomics approach identified an EST-SSR marker with complete linkage to the Ligon lintless-2 genetic locus in cotton (Gossypium hirsutum L.). BMC Genomics. 2011, 12: 445-PubMed CentralView ArticlePubMedGoogle Scholar
- Naoumkina M, Hinchliffe DJ, Turley RB, Bland JM, Fang DD: Integrated metabolomics and genomics analysis provides new insights into the fiber elongation process in Ligon lintless-2 mutant cotton (Gossypium hirsutum L.). BMC Genomics. 2013, 14: 155-PubMed CentralView ArticlePubMedGoogle Scholar
- Liu K, Sun J, Yao L, Yuan Y: Transcriptome analysis reveals critical genes and key pathways for early cotton fiber elongation in Ligon lintless-1 mutant. Genomics. 2012, 100 (1): 42-50.View ArticlePubMedGoogle Scholar
- Bolton JJ, Soliman KM, Wilkins TA, Jenkins JN: Aberrant expression of critical genes during secondary cell wall biogenesis in a cotton mutant, Ligon Lintless-1 (Li-1). Comp Funct Genomics. 2009, Article ID 659301-Google Scholar
- Zhao PM, Wang LL, Han LB, Wang J, Yao Y, Wang HY, Du XM, Luo YM, Xia GX: Proteomic identification of differentially expressed proteins in the Ligon lintless mutant of upland cotton (Gossypium hirsutum L.). J Proteome Res. 2010, 9 (2): 1076-1087.View ArticlePubMedGoogle Scholar
- Kohel RJ, Benedict CR, Jividen GM: Incorporation of (14C) glucose into crystalline cellulose in aberrant fibers of a cotton mutant. Crop Sci. 1993, 33 (5): 1036-1040.View ArticleGoogle Scholar
- Karaca M, Saha S, Jenkins JN, Zipf A, Kohel R, Stelly DM: Simple sequence repeat (SSR) markers linked to the Ligon lintless (Li1) mutant in cotton. J Hered. 2002, 93 (3): 221-224.View ArticlePubMedGoogle Scholar
- Rong J, Pierce GJ, Waghmare VN, Rogers CJ, Desai A, Chee PW, May OL, Gannaway JR, Wendel JF, Wilkins TA, et al: Genetic mapping and comparative analysis of seven mutants related to seed fiber development in cotton. Theor Appl Genet. 2005, 111 (6): 1137-1146.View ArticlePubMedGoogle Scholar
- Updegraff DM: Semi-micro determination of cellulose in biological materials. Anal Biochem. 1969, 32 (3): 420-424.View ArticlePubMedGoogle Scholar
- Boylston EK, Thibodeaux DP, Evans JP: Applying microscopy to the development of a reference method for cotton fiber maturity. Textile Res J. 1993, 63 (2): 80-87.View ArticleGoogle Scholar
- Taliercio EW, Boykin D: Analysis of gene expression in cotton fiber initials. BMC Plant Biol. 2007, 7: 22-PubMed CentralView ArticlePubMedGoogle Scholar
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, et al: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29 (4): 365-371.View ArticlePubMedGoogle Scholar
- Benedito VA, Torres-Jerez I, Murray JD, Andriankaja A, Allen S, Kakar K, Wandrey M, Verdier J, Zuber H, Ott T, et al: A gene expression atlas of the model legume Medicago truncatula. The Plant Journal. 2008, 55 (3): 504-513.View ArticlePubMedGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38: W64-W70.PubMed CentralView ArticlePubMedGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676.View ArticlePubMedGoogle Scholar
- Persson S, Wei H, Milne J, Page GP, Somerville CR: Identification of genes required for cellulose synthesis by regression analysis of public microarray data sets. Proc Natl Acad Sci U S A. 2005, 102 (24): 8633-8638.PubMed CentralView ArticlePubMedGoogle Scholar
- Paterson AH, Wendel JF, Gundlach H, Guo H, Jenkins J, Jin D, Llewellyn D, Showmaker KC, Shu S, Udall J, et al: Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature. 2012, 492 (7429): 423-427.View ArticlePubMedGoogle Scholar
- Obayashi T, Nishida K, Kasahara K, Kinoshita K: ATTED-II Updates: condition-specific gene coexpression to extend coexpression analyses and applications to a broad range of flowering plants. Plant and Cell Physiology. 2011, 52 (2): 213-219.PubMed CentralView ArticlePubMedGoogle Scholar
- Fang DD, Xiao J, Canci PC, Cantrell RG: A new SNP haplotype associated with blue disease resistance gene in cotton (Gossypium hirsutum L.). Theor Appl Genet. 2010, 120 (5): 943-953.View ArticlePubMedGoogle Scholar
- Blenda A, Fang DD, Rami J-F, Garsmeur O, Luo F, Lacape J-M: A high density consensus genetic map of tetraploid cotton that integrates multiple component maps through molecular marker redundancy check. PLoS One. 2012, 7 (9): e45739-PubMed CentralView ArticlePubMedGoogle Scholar
- Rong J, Abbey C, Bowers JE, Brubaker CL, Chang C, Chee PW, Delmonte TA, Ding X, Garza JJ, Marler BS, et al: A 3347-locus genetic recombination map of sequence-tagged sites reveals features of genome organization, transmission and evolution of cotton (Gossypium). Genetics. 2004, 166 (1): 389-417.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu JZ, Kohel RJ, Fang DD, Cho J, Van Deynze A, Ulloa M, Hoffman SM, Pepper AE, Stelly DM, Jenkins JN, et al: A high-density simple sequence repeat and single nucleotide polymorphism genetic map of the tetraploid cotton genome. G3, Gene|Genome|Genetics. 2012, 2 (1): 43-58.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Ooijen JW: JoinMap 4.0: Software for the calculation of genetic linkage maps in experimental populations. 2006, Wageningen, The Netherlands: Kyazma B.VGoogle Scholar
- Shao MY, Wang XD, Ni M, Bibi N, Yuan SN, Malik W, Zhang HP, Liu YX, Hua SJ: Regulation of cotton fiber elongation by xyloglucan endotransglycosylase/hydrolase genes. Genet Mol Res. 2011, 10 (4): 3771-3782.View ArticlePubMedGoogle Scholar
- Thaker VS, Saroop S, Singh YD: Physiological and biochemical changes associated with cotton fibre development. IV. Glycosidases and β-1,3-glucanase activities. Ann Bot. 1987, 60 (5): 579-585.Google Scholar
- Li J, Nam KH, Vafeados D, Chory J: BIN2, a new brassinosteroid-insensitive locus in Arabidopsis. Plant Physiol. 2001, 127 (1): 14-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Evans ML: The action of auxin on plant cell elongation. CRC Crit Rev Plant Sci. 1985, 2 (4): 317-365.View ArticlePubMedGoogle Scholar
- Murphy MP: How mitochondria produce reactive oxygen species. Biochem J. 2009, 417 (1): 1-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Zhang Y, Zhang J, Cheng C, Guo X: Molecular characterization of a transient expression gene encoding for 1-aminocyclopropane-1-carboxylate synthase in cotton (Gossypium hirsutum L.). J Biochem Mol Biol. 2007, 40 (5): 791-800.View ArticlePubMedGoogle Scholar
- Wang H, Mei W, Qin Y, Zhu Y: 1-Aminocyclopropane-1-carboxylic acid synthase 2 is phosphorylated by calcium-dependent protein kinase 1 during cotton fiber elongation. Acta Biochem Biophys Sin (Shanghai). 2011, 43 (8): 654-661.View ArticleGoogle Scholar
- Shi YH, Zhu SW, Mao XZ, Feng JX, Qin YM, Zhang L, Cheng J, Wei LP, Wang ZY, Zhu YX: Transcriptome profiling, molecular biological, and physiological studies reveal a major role for ethylene in cotton fiber cell elongation. Plant Cell. 2006, 18 (3): 651-664.PubMed CentralView ArticlePubMedGoogle Scholar
- Hare PD, Moller SG, Huang LF, Chua NH: LAF3, a novel factor required for normal phytochrome A signaling. Plant Physiol. 2003, 133 (4): 1592-1604.PubMed CentralView ArticlePubMedGoogle Scholar
- Li J, Yu M, Geng L-L, Zhao J: The fasciclin-like arabinogalactan protein gene, FLA3, is involved in microspore development of Arabidopsis. The Plant Journal. 2010, 64: 482-497.View ArticlePubMedGoogle Scholar
- Chen YH, Yang XY, He K, Liu M, Li J, Gao Z, Lin ZQ, Zhang Y, Wang X, Qiu X, et al: The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol. 2006, 60 (1): 107-124.View ArticleGoogle Scholar
- Padmalatha KR, Patil DP, Kumar K, Dhandapani G, Kanakachari M, Phanindra ML, Kumar S, Mohan TC, Jain N, Prakash AH, et al: Functional genomics of fuzzless-lintless mutant of Gossypium hirsutum L. cv. MCU5 reveal key genes and pathways involved in cotton fibre initiation and elongation. BMC Genomics. 2012, 13: 624-PubMed CentralView ArticlePubMedGoogle Scholar
- Cernac A, Andre C, Hoffmann-Benning S, Benning C: WRI1 is required for seed germination and seedling establishment. Plant Physiol. 2006, 141 (2): 745-757.PubMed CentralView ArticlePubMedGoogle Scholar
- Brill E, van Thournout M, White RG, Llewellyn D, Campbell PM, Engelen S, Ruan Y, Arioli T, Furbank RT: A novel isoform of sucrose synthase is targeted to the cell wall during secondary cell wall synthesis in cotton fiber. Plant Physiol. 2011, 157 (1): 40-54.PubMed CentralView ArticlePubMedGoogle Scholar
- Takahashi N, Nakazawa M, Shibata K, Yokota T, Ishikawa A, Suzuki K, Kawashima M, Ichikawa T, Shimada H, Matsui M: shk1-D, a dwarf Arabidopsis mutant caused by activation of the CYP72C1 gene, has altered brassinosteroid levels. The Plant Journal. 2005, 42 (1): 13-22.View ArticlePubMedGoogle Scholar
- Qin YM, Hu CY, Pang Y, Kastaniotis AJ, Hiltunen JK, Zhu YX: Saturated very-long-chain fatty acids promote cotton fiber and Arabidopsis cell elongation by activating ethylene biosynthesis. Plant Cell. 2007, 19 (11): 3692-3704.PubMed CentralView ArticlePubMedGoogle Scholar
- Frigerio M, Alabadi D, Perez-Gomez J, Garcia-Carcel L, Phillips AL, Hedden P, Blazquez MA: Transcriptional regulation of gibberellin metabolism genes by auxin signaling in Arabidopsis. Plant Physiol. 2006, 142 (2): 553-563.PubMed CentralView ArticlePubMedGoogle Scholar
- Cowling RJ, Harberd NP: Gibberellins control Arabidopsis hypocotyl growth via regulation of cellular elongation. J Exp Bot. 1999, 50 (337): 1351-1357.View ArticleGoogle Scholar
- Xiao YH, Li DM, Yin MH, Li XB, Zhang M, Wang YJ, Dong J, Zhao J, Luo M, Luo XY, et al: Gibberellin 20-oxidase promotes initiation and elongation of cotton fibers by regulating gibberellin synthesis. J Plant Physiol. 2010, 167 (10): 829-837.View ArticlePubMedGoogle Scholar
- Loguerico LL, Zhang JQ, Wilkins TA: Differential regulation of six novel MYB-domain genes defines two distinct expression patterns in allotetraploid cotton (Gossypium hirsutum L.). Mol Gen Genet. 1999, 261 (4–5): 660-671.PubMedGoogle Scholar
- Perazza D, Vachon G, Herzog M: Gibberellins promote trichome formation by Up-regulating GLABROUS1 in arabidopsis. Plant Physiol. 1998, 117 (2): 375-383.PubMed CentralView 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.View ArticlePubMedGoogle Scholar
- Hwang JU, Vernoud V, Szumlanski A, Nielsen E, Yang Z: A tip-localized RhoGAP controls cell polarity by globally inhibiting Rho GTPase at the cell apex. Curr Biol. 2008, 18 (24): 1907-1916.PubMed CentralView ArticlePubMedGoogle Scholar
- Seo PJ, Lee AK, Xiang F, Park CM: Molecular and functional profiling of Arabidopsis pathogenesis-related genes: insights into their roles in salt response of seed germination. Plant Cell Physiol. 2008, 49 (3): 334-344.View ArticlePubMedGoogle Scholar
- Desprez T, Juraniec M, Crowell EF, Jouy H, Pochylova Z, Parcy F, Höfte H, Gonneau M, Vernhettes S: Organization of cellulose synthase complexes involved in primary cell wall synthesis in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2007, 104 (39): 15572-15577.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu FJ, Liang WH, Zhang TZ: Genetic analysis of homozygous-dominant-surviving ligon lintless recombinants in cotton (Gossypium hirsutum L.). Acta Agronomica Sinica. 2010, 36 (12): 2020-2027.View ArticleGoogle Scholar
- Sasidharan R, Chinnappa CC, Staal M, Elzenga JT, Yokoyama R, Nishitani K, Voesenek LA, Pierik R: Light quality-mediated petiole elongation in Arabidopsis during shade avoidance involves cell wall modification by xyloglucan endotransglucosylase/hydrolases. Plant Physiol. 2010, 154 (2): 978-990.PubMed CentralView ArticlePubMedGoogle Scholar
- Sembdner G, Atzorn R, Schneider G: Plant hormone conjugation. Plant Mol Biol. 1994, 26 (5): 1459-1481.View ArticlePubMedGoogle Scholar
- Chen J, Du X, Zhou X, Zhao H: Levels of cytokinens in the ovules of cotton mutants with altered fiber development. J Plant Growth Regul. 1997, 16: 181-185.View ArticleGoogle Scholar
- Liu J, Ekramoddoullah AK: The family 10 of plant pathogenesis-related proteins: their structure, regulation, and function in response to biotic and abiotic stresses. Physiol Mol Plant Pathol. 2006, 68 (1–3): 3-13.View ArticleGoogle Scholar
- Liu J, Ekramoddoullah AK, Piggott N, Zamani A: Molecular cloning of a pathogen/wound-inducible PR10 promoter from Pinus monticola and characterization in transgenic Arabidopsis plants. Planta. 2005, 221 (2): 159-169.View ArticlePubMedGoogle Scholar
- Markovic-Housley Z, Degano M, Lamba D, von Roepenack-Lahaye E, Clemens S, Susani M, Ferreira F, Scheiner O, Breiteneder H: Crystal structure of a hypoallergenic isoform of the major birch pollen allergen Bet v 1 and its likely biological function as a plant steroid carrier. J Mol Biol. 2003, 325 (1): 123-133.View ArticlePubMedGoogle Scholar
- Persson S, Paredez A, Carroll A, Palsdottir H, Doblin M, Poindexter P, Khitrov N, Auer M, Somerville CR: Genetic evidence for three unique components in primary cell-wall cellulose synthase complexes in Arabidopsis. Proc Natl Acad Sci USA. 2007, 104 (39): 15566-15571.PubMed CentralView ArticlePubMedGoogle Scholar
- Fagard M, Desnos T, Desprez T, Goubet F, Refregier G, Mouille G, McCann M, Rayon C, Vernhettes S, Hofte H: PROCUSTE1 encodes a cellulose synthase required for normal cell elongation specifically in roots and dark-grown hypocotyls of Arabidopsis. Plant Cell. 2000, 12 (12): 2409-2423.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao L, Lv Y, Cai C, Tong X, Chen X, Zhang W, Du H, Guo X, Guo W: Toward allotetraploid cotton genome assembly: integration of a high-density molecular genetic linkage map with DNA sequence information. BMC Genomics. 2012, 13: 539-PubMed CentralView ArticlePubMedGoogle Scholar
- Liebminger E, Huttner S, Vavra U, Fischl R, Schoberer J, Grass J, Blaukopf C, Seifert GJ, Altmann F, Mach L, et al: Class I alpha-mannosidases are required for N-glycan processing and root development in Arabidopsis thaliana. Plant Cell. 2009, 21 (12): 3850-3867.PubMed CentralView ArticlePubMedGoogle Scholar
- Herve C, Dabos P, Bardet C, Jauneau A, Auriac MC, Ramboer A, Lacout F, Tremousaygue D: In vivo interference with AtTCP20 function induces severe plant growth alterations and deregulates the expression of many genes important for development. Plant Physiol. 2009, 149 (3): 1462-1477.PubMed CentralView ArticlePubMedGoogle Scholar
- Hao J, Tu L, Hu H, Tan J, Deng F, Tang W, Nie Y, Zhang X: GbTCP, a cotton TCP transcription factor, confers fibre elongation and root hair development by a complex regulating system. J Exp Bot. 2012, 63 (17): 6267-6281.PubMed CentralView 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.