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 volume 12, Article number: 445 (2011)
Cotton fiber length is an important quality attribute to the textile industry and longer fibers can be more efficiently spun into yarns to produce superior fabrics. There is typically a negative correlation between yield and fiber quality traits such as length. An understanding of the regulatory mechanisms controlling fiber length can potentially provide a valuable tool for cotton breeders to improve fiber length while maintaining high yields. The cotton (Gossypium hirsutum L.) fiber mutation Ligon lintless-2 is controlled by a single dominant gene (Li 2 ) that results in significantly shorter fibers than a wild-type. In a near-isogenic state with a wild-type cotton line, Li 2 is a model system with which to study fiber elongation.
Two near-isogenic lines of Ligon lintless-2 (Li 2 ) cotton, one mutant and one wild-type, were developed through five generations of backcrosses (BC5). An F2 population was developed from a cross between the two Li 2 near-isogenic lines and used to develop a linkage map of the Li 2 locus on chromosome 18. Five simple sequence repeat (SSR) markers were closely mapped around the Li 2 locus region with two of the markers flanking the Li 2 locus at 0.87 and 0.52 centimorgan. No apparent differences in fiber initiation and early fiber elongation were observed between the mutant ovules and the wild-type ones. Gene expression profiling using microarrays suggested roles of reactive oxygen species (ROS) homeostasis and cytokinin regulation in the Li 2 mutant phenotype. Microarray gene expression data led to successful identification of an EST-SSR marker (NAU3991) that displayed complete linkage to the Li 2 locus.
In the field of cotton genomics, we report the first successful conversion of gene expression data into an SSR marker that is associated with a genomic region harboring a gene responsible for a fiber trait. The EST-derived SSR marker NAU3991 displayed complete linkage to the Li 2 locus on chromosome 18 and resided in a gene with similarity to a putative plectin-related protein. The complete linkage suggests that this expressed sequence may be the Li 2 gene.
Cotton seed fibers are initially ovule epidermal cells that terminally differentiate into fiber cells typically on the 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 [1, 2]. The rate of fiber elongation peaks at approximately 6 to 12 DPA and nears cessation at 22 DPA . During peak elongation fiber cells can increase in length at rates of 2 mm/day or more depending on environment and genotype [1, 4]. The length of fibers is mostly variety specific, but can also be affected by environmental conditions such as temperature during the elongation stage of development . The elongation stage is followed by a brief period known as the transition stage that usually begins from 12 to 16 DPA in field conditions and depending on environmental factors such as lower temperatures that are shown to delay the onset of the transition stage . The SCW stage immediately follows transition and is characterized by a dramatic increase in SCW-related gene transcripts like cellulose synthases and changes in cell wall composition as large amounts of cellulose are deposited in the SCW . The SCW stage persists until about 32 DPA at which time the fiber cell is composed of approximately 95% cellulose with the remaining 5% of non-cellulosic materials comprised of proteins, polysaccharides, pectins, and waxes that reside mostly in the PCW and cuticle [3, 8]. The final stage of fiber development is maturation that ceases from 40 to 60 DPA depending on environment and genotype . At this time the cotton bolls crack and open, exposing the seed fibers to external ambient conditions causing them to desiccate and take on the fluffy appearance normally associated with cotton fibers.
Varieties of cultivated Upland cotton (Gossypium hirsutum L.) that display fiber mutation phenotypes including lintless and fuzzless seeds were first described in the early twentieth century [10, 11]. Currently, numerous naturally occurring cotton fiber mutations have been identified globally and characterized at the genetic, and more recently, gene expression levels [12–16]. These fiber mutations include, among others, the glabrous seeds in the fiberless mutant lines MD17, SL1-7-1, and XZ142w [17, 18]; seeds with only lint fibers and no fuzz fibers in the Naked seed lines N 1  and n 2 ; and seeds that are described as extremely short lint fibers in the Ligon lintless-1 (Li 1 )  and Ligon lintless-2 (Li 2 )  mutant lines. The fiber mutations of N 1 , Li 1 , and Li 2 are single gene dominant traits [20, 21] while the n 2 fiber mutation is a single gene recessive trait . Recent genetic studies on the Li 2 mutation also indicate that it may have incomplete penetrance as evidenced by mutant and WT fibers at different boll locations on the same plants , or possibly phenotype variation due to epigenetics. These naturally occurring mutants and their wild-type fiber NILs provide a unique and powerful model system to study cotton fibers at various stages of development including initiation, elongation, and secondary cell wall biosynthesis.
The Li 1 gene was mapped to chromosome 22 using SSR  and RFLP  markers while the Li 2 gene was mapped to chromosome 18 by phenotype association with cotton aneuploid stocks , and linkage analysis by RFLP markers . More recently, a draft of the physical map of the diploid cotton D-genome progenitor G. raimondii was released and used along with tetraploid cotton A- and D- subgenome genetic maps to generate a consensus genetic-physical map of the cotton genome that included flanking markers for the Li 2 gene on chromosome 18 . These were RFLP markers designated A1552 and Gate4BC11 that flanked the Li 2 gene in an 8.9 cM region and were mapped using an interspecific F2 mapping population composed of 158 individuals from a cross between G. hirsutum Li 2 and G. barbadense cv. Pima S-7. Previously, another RFLP marker designated Gate4BF10 was reported to flank the Li 2 gene along with A1552 in a 1.5 cM region of chromosome 18 . However, the more recently released cotton consensus genetic-physical map developed by the same laboratory indicated that Gate4BF10 was mapped at two locations on chromosome 18 , which leave the mapping accuracy of this marker in doubt. In a separate study, two Li 2 F2 segregating populations were developed and used to screen SSR markers for linkage to the Li 2 genetic locus. The closest SSR marker was mapped to chromosome 18 and located 6.051 cM from the Li 2 gene in an interspecific segregating population, and 9.266 cM from the Li 2 gene in an intraspecific segregating population .
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 WT fibers that are typically greater than 20 mm in length [20, 21, 28]. 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 Li 1 mutant, which exhibits pleiotropy in the form of severely stunted and deformed plants in both the homozygous dominant and heterozygous state , the Li 2 mutant plants appear healthy and morphologically identical to the homozygous recessive wild-type plants with the exception of shorter seed fibers .
Cytological evidence suggests that the seed fibers of Li 1 mutants undergo initiation in the same manner as WT fibers, but begin to show some distorted morphological features during the early elongation stage of development . Since the seed fibers of Li 1 and Li 2 fibers are shortened lint fibers, these cotton mutants represent excellent candidates to study the molecular mechanisms of fiber elongation. A recent gene expression study using microarrays on Li 1 mutant and WT cotton NILs that focused on the SCW stage of fiber development identified genes potentially responsible for the phenotypic differences observed in mutant Li 1 fibers compared to WT fibers . Several genes in particular were differentially expressed during SCW biosynthesis that could potentially be involved with the Li 1 phenotype including EXPANSINS, tubulin genes, sucrose synthase (SuSy), and genes encoding MYB transcription factors .
Since extensive research has been ongoing with the Li 1 mutant, our laboratory selected the Li 2 mutant in an established near-isogenic state with the Upland cotton variety DP5690 as a model system to study fiber elongation events using a combined functional and structural genomics approach of microarray gene expression and molecular marker analysis. The Li 2 mutant was also selected due to concerns over the pleiotropic effects of the Li 1 mutation on development timing and the desire to harvest fiber samples from mutant and WT plants simultaneously. Understanding the molecular events that control fiber elongation and identifying regulatory elements involved in this process can provide cotton researchers and breeders with means of improving fiber length while maintaining yield either through marker-assisted selection or a transgenic approach. The main objective of this research was to identify genes that were differentially expressed during the development of WT and mutant Li 2 fibers and convert the gene expression data into portable molecular markers for use in association mapping to identify the Li 2 locus.
Here we report: 1) the development of two Li 2 NILs of cotton (G. hirsutum) in the backcross five (BC5) generation; 2) no apparent phenotypic differences in seed fibers of mutant Li 2 Li 2 plants compared to WT li 2 li 2 plants during the initiation and early elongation stages of fiber development; 3) mapping the Li 2 locus with SSR markers; 4) the identification of genes differentially expressed in fibers of mutant Li 2 plants compared to WT plants using microarray gene expression analysis on selected developmental time-points; 5) confirmation of microarray gene expression profiles by RT-qPCR; 6) successful conversion of differentially expressed genes into EST-SSR markers; and 7) the identification of an EST-SSR marker representing a putative plectin-related protein or regulatory element that has complete linkage to the Li 2 genetic locus suggesting it may be the Li 2 gene.
Plant materials and greenhouse experimental design
Two NILs of Li 2 Upland cottons that were homozygous dominant (Li 2 Li 2 ) and homozygous recessive (li 2 li 2 ) for the Li 2 locus were developed in a backcross program at Stoneville, MS in field and greenhouse environments. Mutant Texas marker-1 (TM-1) cotton plants containing the Li 2 gene were crossed with the Upland cotton variety DP5690 and F1 progeny were backcrossed for five generations (BC5) by SSD to DP5690 which served as the recurrent parent in each backcross. The DP5690 recurrent parent was a pure inbred line that was self-pollinated for nine generations via SSD. Progeny in each backcross were selected based on phenotype for the Li 2 short-fiber mutation. The pedigree of the two Li 2 NILs is detailed in Additional file 1.
A total of 102 Li 2 cotton plants and 80 WT li 2 li 2 cotton plants were planted in the greenhouse on six tables. The mutant Li 2 plants were six BC5F3 lines that originated via SSD. A total of 17 individual plants were grown for each line to confirm that they were homozygous Li 2 Li 2 . Once the genotypes were confirmed the populations were culled to 72 mutant Li 2 Li 2 plants and 72 WT li 2 li 2 plants that were placed in the same greenhouse on six tables in a randomized complete block design. The individual mutant and WT cotton plants were labeled into three pools representing three biological replicates. Cotton bolls were harvested at the following time-points during development: -3, -1, 0, 1, 3, 5, 8, 12, 16, and 20 DPA. Bolls from the same cotton line, biological replicate, and developmental time-point were harvested from all six tables to account for environmental variability within the greenhouse and bulked for subsequent analyses. The number of bolls per bulked sample varied according to developmental time-point, with a greater number of bolls required for the earliest time-point to ensure sufficient biological material and progressively fewer bolls required for each successive time-point. For example, at opposite ends of the developmental time-course, ovules from approximately 20 - 30 bolls were bulked for each -3 DPA sample, and ovules with fibers attached from approximately 8 -10 bolls were bulked for each 20 DPA sample. Harvested bolls were placed immediately on ice and transported to the laboratory where they were dissected on ice and the majority of the ovules frozen in liquid nitrogen and stored at -80°C. A small number of ovules from each sample from the -1 to 5 DPA time-points were used for SEM as described herein.
A mutant Li 2 Li 2 homozygous plant was used as the female in a cross with its near-isoline WT li 2 li 2 DP5690. One hundred and thirty-six F2 plants were planted in the field in Stoneville, MS in 2009. The Li 2 trait of each F2 progeny plant was evaluated twice at approximately 30 DPA and after boll maturation and opening.
SSR marker analysis and genetic mapping
Young leaves were collected from each one of the F2 plants in the described mapping population. Total DNA was extracted from fresh leaves using 2.0% hexadecyltrimethylammonium bromide . DNA was purified using Omega EZNA® DNA isolation column (Omega Bio-Tek, Norcross, GA). As previously reported, the Li 2 locus resides on chromosome 18 [24, 25, 27]. To rapidly identify SSR markers closely linked to the Li 2 locus, we first selected 86 SSR markers that were previously mapped on either chromosome 18 or its homeologous chromosome 13 based on several published maps [31–35]. The probes of the RFLP markers reported by Rong et al. (2005)  were not available to us, and thus were not evaluated in our population. Bulked segregant analysis  was then used to identify potentially linked markers. For the Li 2 bulk, DNA of 10 F2 plants that had the Li 2 phenotype were pooled at equal ratio and diluted to 10 ng/μL. The WT bulk consisted of pooled DNA from 10 F2 progeny that had normal lint phenotype. SSR primers that generated polymorphic patterns between bulks were tested using the 20 individual DNA samples that were included in the bulks. The markers linked to the Li 2 locus were analyzed on 136 individual F2 progeny plants as previously described . All SSR primer sequences can be obtained from Cotton Marker database (http://www.cottonmarker.org) except DPL0547. The primer sequences of the SSR markers associated with Li 2 locus are listed in Additional file 2. Segregation data for the Li 2 trait and SSR markers were mapped using program JoinMap3.0 [38, 39] with logarithm of odds score = 25.
Scanning electron microscopy
To prepare the samples for SEM analysis, cotton ovules were placed in tissue fixative consisting of 3% (v/v) glutaraldehyde in 0.1 M sodium phosphate, pH 7.0 and stored at 4°C. The time-points utilized for SEM were -1 to 5 DPA. After fixation, the cotton ovules were dehydrated in a graded ethanol series starting from 20% (v/v) up to 100% (v/v) ethanol. After three changes of 100% ethanol, the ovules were placed in American Optical microporous specimen capsules under 100% ethanol  and critical point dried from liquid carbon dioxide by standard methodology in a Ladd Critical Point Dryer model 28,000 (Ladd Research, Williston, VT). The ovules were mounted on standard Cambridge SEM stubs using double-stick Avery photo tabs, #06001. The SEM mounts were coated with 60/40 gold/palladium using a Hummer™ II Sputter Coater (Ladd Research, Williston, VT) to a thickness of 200 nm. The specimens were examined in a XL30 Environmental Scanning Electron Microscope (FEI Company, Hillsboro, OR) at an accelerating voltage from 10-15 kV under high vacuum conditions.
Cotton fiber total RNA isolation
Cotton fibers were isolated from developing ovules using a glass bead shearing technique to separate fibers from the ovules . Total RNA was isolated from detached fibers using the Sigma Spectrum™ Plant Total RNA Kit (Sigma-Aldrich, St. Louis, MO) with the optional on-column DNase1 digestion according to the manufacturer's protocol. The concentration of each RNA sample was determined using a NanoDrop 2000 spectrophotometer (NanoDrop Technologies Inc., Wilmington, DE). The RNA quality for each sample was determined by RNA integrity number (RIN) using an Agilent Bioanalyzer 2100 and the RNA 6000 Nano Kit Chip (Agilent Technologies Inc., Santa Clara, CA) with 250 ng of total RNA per sample.
Reverse transcription quantitative real-time PCR
The experimental procedures and data analysis related to RT-qPCR were performed according to the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines . The cDNA synthesis reactions were performed using the iScript™ cDNA Synthesis Kit (Bio-Rad Laboratories, Hercules, CA) according to the manufacturer's instructions with 1 μg of total RNA per reaction used as template. Control cDNA synthesis reactions to check for genomic DNA contamination during RT-qPCR consisted of the same template and components as the experimental reactions without the reverse transcriptase enzyme. The RT-qPCR reactions were performed with iTaq™ SYBR® Green Supermix (Bio-Rad Laboratories) in a Bio-Rad CFX96 real time PCR detection system. Thermal cycler parameters for RT-qPCR were as follows: 95°C 3 minutes, 50 cycles of 95°C 15 seconds, 60°C 30 seconds. A dissociation curve was generated and used to validate that a single amplicon was present for each RT-qPCR reaction. The calculations for amplification efficiencies of the target and reference genes, RNA inhibition assays, and the relative quantifications of the different target gene transcript abundances were performed using the comparative Cq method as described in the ABI Guide to Performing Relative Quantitation of Gene Expression Using Real-Time Quantitative PCR (Applied Biosystems, Foster City, CA) with the following modification: the average of three reference gene Cq values was determined by taking the geometric mean which was used to calculate the ΔCq values for the individual target genes . The endogenous reference genes used in the RT-qPCR reactions were the 18S rRNA (Genbank accession U42827), ubiquitin-conjugating protein (Genebank AI730710), and alpha-tubulin 4 (TubA4, Genbank AF106570) . The reference and target gene primer sequences, and target gene descriptions including in silico specificity screens using BLASTx  are shown in Additional file 2.
Initial analyses of specific transcripts in fibers of Li 2 mutant and WT plants were performed to determine the most informative time-points to use for microarray analysis. The genes selected were based on previous studies that indicated developmental regulation of these genes during different stages of fiber development, specifically elongation and SCW biosynthesis. The elongation stage-related genes α-expansin1 (GhExp1)  and Cu/Zn superoxide dismutase (GhCSD1)  and the SCW-related genes cellulose synthase2 (GhCesA2)  and a β-1,3-glucanase-like gene  that was previously shown to be up-regulated in the fiber SCW stage  were selected for preliminary RT-qPCR. The nucleotide primer sequences and gene accessions are shown in Additional file 2. The RT-qPCR reactions were performed with cDNA from fibers at initiation, elongation, and SCW thickening stages at the following time-points: DOA, 1, 3, 5, 8, 12, 16, and 20 DPA. Three biological replications were used for each time-point sample.
Microarray hybridizations and data analysis
The microarray experiments conducted in this study followed the minimum information about a microarray experiment (MIAME) guidelines . The microarray utilized was the commercially available Affymetrix GeneChip® Cotton Genome Array (Affymetrix Inc., Santa Clara, CA), that represents 21,854 transcripts gathered from four species of cotton (G. arboreum, G. barbadense, G. hirsutum, and G. raimondii) and a variety of tissue types including fibers at the initiation, elongation, and SCW biosynthesis stages of development. For each sample 500 ng of cotton fiber RNA was utilized for labeling using the Affymetrix GeneChip® 3' IVT Express Kit and Cotton Genome Array hybridizations were performed according to standardized Affymetrix protocols. The developmental time-points from Li 2 mutant and WT fibers used for microarray hybridizations were 0, 8, and 12 DPA with two biological replicates used for each time-point sample. Procedures for data normalization and assessment of statistically and biologically significant genes were performed as described by Benedito et al. . The Affymetrix microarray dataset was deposited in the ArrayExpress database with the expression number E-MEXP-3306.
To gain insight into the biological processes represented by the differentially expressed genes in the WT and mutant Li 2 fibers, an over-representation analysis (ORA) was performed on the genes that were more abundantly transcribed and statistically significant in fibers of each cotton NIL at each time-point using the program Blast2GO . The Blast2GO program utilizes a Fisher's exact test for ORA and corrects the p-value to control the increased false discovery rate (FDR) associated with multiple hypotheses testing. For this experiment the FDR-corrected p-value cutoff for significance was 0.05.
Confirmation of homozygosity in the mutant Li 2 lines
Segregation analysis of the mutant phenotype was accomplished by visual inspection of the ovules from each one of the 102 plants in the six mutant lines. An example of the Li 2 mutant and WT fiber development stages from the DOA to maturity are shown in Figure 1 with slight differences in fiber length evident as early as 5 DPA. Definitive confirmation of the mutant Li 2 phenotype was possible by 16 DPA and indicated that five of the six lines were homozygous Li 2 Li 2 , while one line was heterozygous Li 2 li 2 . The heterozygous line was removed from the greenhouse and both the mutant and WT population were reduced to 72 individual plants each to achieve the randomized complete block design of the experiment.
Morphology of mutant Li 2 fibers during initiation and early elongation
SEM analyses revealed no discernable differences in the appearances of ovules and fibers from Li 2 mutant and WT (Figure 2). The WT and mutant fibers appeared to undergo normal initiation as visualized on the DOA with no apparent differences in the distribution or density of fiber initials on the ovule surfaces (Figure 2, C & 2I). Likewise, WT and mutant fiber morphology and length appeared the same during the early elongation stage of fiber development up to 5 DPA (Figure 2, D-F; J-L). While no differences between WT and mutant fibers were evident in the SEM (Figure 2), there does appear to be some slight length differences at 5 DPA as shown in Figure 1, so it can be stated that the mutant fiber phenotype becomes evident by approximately 5 DPA. These data are consistent with comparable fiber length measurements obtained from developing fibers of TM-1, Li 2 , and F3 progeny derived from a cross of TM-1 and Li 2 .
Preliminary RT-qPCR using genes differentially expressed during cotton fiber development
The expression of the elongation stage-related gene GhExp1 followed a similar pattern of expression in WT and Li 2 mutant fibers, with transcript abundance decreasing at the beginning of the SCW stage from 16 - 20 DPA. However, the transcript abundance of GhExp1 was a significantly decreased (2.24-fold) in Li 2 mutant fibers compared to WT fibers during the presumed peak of fiber elongation at 8 DPA. The transcript abundance of the second elongation stage-related gene, GhCSD1, followed the previously shown pattern of expression during fiber development with transcript abundance higher during fiber elongation compared to the SCW biosynthesis stage . In fibers of the Li 2 mutant, the expression levels of GhCSD1 remained significantly lower compared to WT fibers during elongation stage time-points (5 - 12 DPA) and remained at a seemingly basal level of expression with no change in transcript abundance in mutant fibers over time-points normally associated with elongation and SCW biosynthesis from 5 - 20 DPA (Figure 3A).
The SCW-related genes GhCesA2 and β-1,3-glucanase-like followed the expected expression profiles in WT and Li 2 mutant fibers with up-regulation occurring concurrent with onset of the transition stage from 12 - 16 DPA and SCW biosynthesis 16 - 20 DPA (Figure 3B). While the pattern of expression for the SCW-related genes was the same in WT and Li 2 mutant fibers, the transcript abundance of both genes was significantly higher in the WT fibers at 16 and 20 DPA for GhCesA2 (3.31- and 4.62-fold, respectively) and 20 DPA for the β-1,3-glucanase-like gene (6.64-fold) as shown in Figure 3B.
Microarray gene expression analysis
The initial RT-qPCR analysis indicated significantly different expression levels of the elongation stage-related genes GhExp1 and GhCSD1 during time-points associated with the elongation stage of fiber development (Figure 3A). Since there were no apparent morphological differences between WT and Li 2 mutant fibers during initiation or early elongation (Figure 2), this microarray study focused primarily on time-points related to fiber elongation, more specifically, 8 and 12 DPA. These time points typically represent peak rates of elongation and the beginning of the transition stage when fibers are still elongating at lower rates. An additional time-point, DOA, was added to the microarray experiment to serve as a reference and also to confirm that there were no significant differences in fiber initiation at the gene expression level. Global gene expression analysis by microarray revealed a significant number of differentially expressed genes in fibers of the mutant and WT NILs at 8 and 12 DPA. At 8 and 12 DPA, there were 1079 and 1106 genes, respectively, more abundantly transcribed (≥ 2-fold; < Bonferroni-corrected p-value threshold 2.07194E-06) in fibers of the WT NIL compared to the mutant fibers. In fibers of the Li 2 mutant NIL 1064 and 1473 genes were significantly and more abundantly transcribed at 8 and 12 DPA, respectively, compared to WT fibers. Far fewer genes were differentially expressed in fibers of the WT and mutant lines in the DOA time-point with 17 genes more abundantly transcribed in WT fibers and 141 genes more abundantly transcribed in mutant fibers.
On the DOA, there were no significantly enriched biological processes represented by the differentially expressed genes in WT or mutant fibers. This result corresponded with the lack of any observable differences in fiber morphology between WT and mutant Li 2 fibers during initiation and early elongation. In 8 DPA fibers of the WT line, over-represented biological processes included, among others, receptor activity, signaling, cytoskeleton, and cell growth, while genes more abundantly transcribed at 8 DPA in mutant Li 2 fibers were enriched in stress and stimulus responses. The genes involved in signaling, cell growth, and cytoskeleton biological processes at 8 DPA in the WT fibers included genes whose products are known to regulate cell elongation and expansion in Arabidopsis thaliana. Among these were genes known to encode protein kinases that regulate cell elongation in A. thaliana through brassinosteroid (BR) signaling pathways such as orthologues of BRASSINOSTEROID INSENSITIVE1 (BRI1), and HERCULES1 (HERK1) and HERCULES2 (HERK2). Some of the genes that encode proteins with the potential to alter or regulate cell elongation are shown in Table 1 along with the functionally characterized phenotypes resulting from inhibition or over-expression in A. thaliana. Genes of interest with higher transcript abundance in Li 2 mutant fibers compared to WT fibers at 8 and 12 DPA, particularly those related to stress and stimulus responses, involved phytohormone signaling pathways for ethylene and cytokinin biosynthesis and regulation of phytohormone homeostasis. The potential implications of genes and biological processes over-represented in WT and Li 2 mutant fibers during elongation are discussed later and summarized in Table 1. The complete results of the ORA analysis are presented in Additional file 3 which includes the individual microarray probeset IDs for the over-represented biological processes for each cotton NIL and time-point.
Corroboration of the microarray data
A total of twelve genes were selected for corroboration of the microarray data, in addition to the cotton Cu/Zn superoxide dismutase (GhCSD1) gene that was utilized for preliminary RT-qPCR data to select fiber developmental time-points for microarray analysis (Figure 3). Of the four genes selected for the preliminary RT-qPCR, only the GhCSD1 gene is accurately represented in the Affymetrix microarray data with the probeset ID Ghi.2036.1.S1_s_at (Table 2, Additional file 2). The GhCesA2 gene shares 100% homology with the probesets nucleotide sequence, but is below the limit of detection during the initiation and elongation stages selected for microarray analysis. The GhExp1 and β-1,3-glucanase-like genes used to design primers for RT-qPCR (Additional file 2) do not share enough homology with the Affymetrix nucleotide sequences in the regions corresponding to the probsets to allow for accurate comparison of transcript abundance between microarray and RT-qPCR.
The genes selected for microarray corroboration by RT-qPCR were chosen primarily from differentially expressed genes (≥ 2-fold; < Bonferroni-corrected p-value threshold 2.07194E-06) in WT and mutant fibers at 8 and 12 DPA. Some selected genes were more abundantly expressed at 8 and/or 12 DPA in WT fibers, while others were more abundantly expressed in mutant fibers at the same time-points. These genes encode regulatory elements such as transcription factors including the G. hirsutum fiber initiation-related MYB transcription factor GhMYB109; gene products involved in carbohydrate metabolism such as beta-galactosidase, glycosyltransferases; and a cellulose-synthase-like cell wall-related protein. Genes that were not differentially expressed during elongation time-points were also selected to confirm the validity of the microarray data. These included genes encoding a R2R3-MYB transcription factor (probeset ID Ghi.9209.1.S1_at) and sucrose synthase (probeset ID Ghi.8665.1.S1_s_at). The microarray data were successfully corroborated for most time-points with the exception of the SANT/MYB transcription factor gene (probeset ID Ghi.1711.1.S1_s_at) at 8 DPA. The complete results of the microarray data and the corresponding RT-qPCR data for all twelve genes and the GhCSD1 gene are summarized in Table 2 and Additional file 4.
Fine mapping the Li 2 locus region with codominant SSR markers
Of the 86 SSR markers screened, 6 (6.97%) were polymorphic between two DNA bulks from Li 2 mutant and WT plants. Except for the marker TMB2295 that had 5 recombinations with the Li 2 trait when the 20 individuals comprising the bulks were tested, the markers CIR216, DC30513, DPL0547, DPL0922, and MUSB1135 co-segregated with the Li 2 trait. These 5 markers were further evaluated in the whole F2 population. A map was constructed around the Li 2 region (Figure 4). The Li 2 locus was flanked by markers DPL0547 and DPL0922 at genetic distances of 0.87 cM and 0.52 cM, respectively.
Identification of an EST-SSR marker with complete linkage to the Li 2 locus
In an effort of using gene expression data to identify candidate molecular markers for association mapping of the Li 2 gene, we conducted a BLAST search  of the Cotton Marker Database (CMD) (http://www.cottonmarker.org) using consensus sequences and singletons representing differentially expressed genes from the microarray data. The search criteria were set at a stringent level (e-value > E-50) in order to select sequences with a relatively high sequence alignment length and homology to the microarray query sequences. Query sequences included genes more abundantly expressed (≥ 2-fold) in fibers of both mutant and WT cottons at 8 and 12 DPA. The numbers of query sequences from WT fibers were 1031 and the numbers of query sequences from mutant fibers were 982. A total of 692 SSR markers were identified based on sequence homology between the microarray query nucleotide sequences and the SSR source nucleotide sequences. Of these 692 SSR markers, 179 were previously mapped by different research groups [31, 33–35, 55], and 14 of them were mapped on chromosome 18 or 13. Except markers DPL0249 and TMB1767 that were developed from genomic DNA libraries, the other 12 were EST-SSRs (http://www.cottonmarker.org). We analyzed these 14 SSR markers in the F2 population, and only one (NAU3991) segregated among progeny plants. The marker NAU3991 was developed from G. raimondii EST  and displayed complete linkage to the Li 2 genetic locus in our intraspecific F2 population composed of 136 individuals (Figure 4). The Blastx hits for the full-length cDNA with highest similarity in the NCBI non-redundant database included a putative transcription factor isolated from Ricinus communis; plectin-related proteins from A. thaliana and Vitis vinifera; and unknown proteins from Populus trichocarpa and Glycine max. The full-length predicted protein sequences are highly conserved among cotton and the other plant species for all of the above mentioned protein sequences. The NCBI conserved domain database (CDD)  indicated a putative DNA binding domain present in the predicted protein sequence.
The microarray probeset that led to the discovery of the linkage between NAU3991 and Li 2 locus, Ghi.8501.1.A1_at was represented as TC276355 in the Dana-Farber Cancer Institute Cotton Gene Index (DFCI) 11.0 database http://compbio.dfci.harvard.edu/tgi/. The RT-qPCR primers used to corroborate the microarray data for probeset Ghi.8501.1.A1_at was designed from TC276355 (Additional file 2). The expression profile of TC276355 matched the microarray data for probesets Ghi.8501.1.A1_at and indicated that the gene was up-regulated in WT fibers during the fiber elongation and transition stages (8 - 16 DPA) with significantly higher transcript abundance in WT fibers compared to mutant fibers at 8 and 12 DPA (Figure 5).
Naturally occurring cotton fiber mutations present an excellent opportunity to study the molecular events controlling fiber development. When a monogenic fiber mutation in the homozygous state is coupled with its homozygous WT counterpart, especially as NILs, an ideal model system is created for a comparative study of fiber development. Here we have presented such a model system in the Li 2 short-fiber mutation as an advanced (BC5) NIL in the G. hirsutum DP5690 background. In addition to the Li 2 short-fiber mutant phenotype being caused by a single dominant allele, no pleiotropic effects are observed in the development, morphology, or fecundity of mutant cotton plants as either homozygotes or heterozygotes. Also, no pleiotropic effects of the Li 2 mutation are observed in mutant fiber development before the elongation stage. In this regard, the Li 2 mutant is useful to study molecular events specific to and controlling fiber elongation.
It was previously demonstrated that the fibers of another dominant, monogenic short-fiber mutation similar to Li 2 , designated Li 1 , have a thicker SCW than WT fibers and Li 2 mutant fibers. This was demonstrated by measuring the ratio of fiber weight-to-length to infer the degree of SCW thickening in fibers of Li 1 , Li 2 , TM-1, and normal homozygous recessive Li 2 segregates (F3) from a cross between Li 2 and TM-1 . The thickened SCW of Li 1 fibers was further confirmed in another study that measured the incorporation of [14C] glucose in the SCW and indicated a five-fold increase in crystalline cellulose synthesis per millimeter of fiber in the Li 1 mutant compared to WT fibers . This suggests that the pleiotropic effect of the short lint fiber Li 1 mutation also affects multiple stages of fiber development, while the Li 2 mutation remains specific to fiber length. It was also speculated that the thickened SCW walls of Li 1 fibers may be due to increased SCW cellulose synthase activity or an increased density of cellulose synthase subunits per unit length of fiber causing an increase in cellulose deposition in the SCW . The later speculation would likely result in an increase in the transcript abundances of SCW CesA genes such as GhCesA2 in Li 1 mutant fibers compared to WT fibers during the SCW stage. A recent comparative gene expression study on the SCW stage of Li 1 mutant and WT TM-1 fibers would have provided insight into this theory, however the SCW CesA genes were not discussed and the microarray dataset was not made publically available . The authors did report 8.8-fold higher transcript levels of a sucrose synthase (SuSy) gene (Genbank accession U73588) in Li 1 mutant fibers compared to WT TM-1 fibers that could support the observed increased cellulose deposition in the SCW of Li 1 mutant fibers [12, 57]. Our RT-qPCR data for the SCW-related genes GhCesA2 and β-1,3-glucanase-like indicated significantly higher levels of expression at the beginning of SCW synthesis in WT fiber compared to Li 2 mutant fibers (Figure 3B). This would be the converse of the Li 1 fibers if the density of CesA subunits per unit length of fiber were higher in Li 1 mutant fibers compared to WT fibers and supports the fiber elongation-specific, non-pleiotropic nature of the Li 2 mutation. The SuSy gene that was selected as a non-differentially expressed gene for microarray corroboration also supports this hypothesis (Table 2; Additional file 4). The RT-qPCR on all developmental time-points in our study using primers specific for the G. hirsutum SuSy gene (Genbank accession U73588) matched the GhCesA2 expression at 20 DPA (Figure 3B) with statistically significant 2.54-fold higher SuSy transcript abundance in WT fibers compared to mutant fibers (Additional file 4).
The significantly higher expression levels of the GhCSD1 gene in WT fibers compared to mutant fibers coincided well with a previously postulated model that suggested a role for the modulation of ROS in fiber cell elongation. The model suggested that short fiber cotton species such as the diploid G. longicalyx are unable to regulate ROS accumulation during elongation, leading to stress conditions, an increase in the expression of stress-related genes, and a decrease in the rate of fiber cell elongation . Comparison of the fiber gene expression profiles of G. longicalyx and diploid G. arboreum during elongation (5 DPA) revealed over-representation of biological processes involved in response to stress, presumably due to over-accumulation of ROS, in the shorter fibers of G. longicalyx. Results indicating the same differences in ROS homeostasis were obtained by the same research laboratory in two more microarray gene expression studies that compared short and long fiber cotton species during elongation [59, 60]. A similar over-representation of stress response-related genes (GO:0006950) was observed in Li 2 mutant fibers compared to WT fibers at 8 DPA in our study, including genes responsive to oxidative stress and involved in ROS homeostasis (Additional file 3). A total of 105 genes were included in the over-representation of stress responses in Li 2 mutant fibers at 8 DPA and included ROS-related genes such as an NADPH-dependent aldo-keto reductase (AKR) with high similarity to the A. thaliana AKR4C9 protein involved in the reduction of reactive α,β-unsaturated carbonyls produced through lipid peroxidation ; an alternative oxidase highly similar to the stress-responsive A. thaliana ALTERNATIVE OXIDASE 1A (AOX1A); a G. hirsutum ascorbate peroxidase (GhAPX1) involved in H2O2 homeostasis ; a G. hirsutum class III peroxidase (GhPOX1) suggested to play a role in ROS production ; a transcription factor highly similar to the A. thaliana MULTIPROTEIN BRIDGING FACTOR 1C (AtMBF1C) that is over-expressed in response to multiple stimuli including elevated H2O2 levels ; and several putative APETELA2 (AP2)-domain transcription factors that are also up-regulated in transgenic A. thaliana over-expressing AtMBF1C . Given that some of these gene products are involved in the production of ROS required for normal fiber cell elongation, while others are involved in the reduction oxidative stress and repair of oxidative damage, it is possible that the mechanisms of ROS modulation and homeostasis are compromised in Li 2 mutant fibers.
Numerous genes involved in biological processes such as cell expansion and elongation that are up-regulated in the WT fibers compared to mutant fibers during the time-points coinciding with peak rates of fiber elongation. Genes that were more abundantly transcribed (≥ 2-fold) and statistically significant in either Li 2 mutant or WT fibers during elongation were closely examined to discern the possible effects of their up- or down-regulation based on functional analysis in cotton and heterologous systems. Some examples of these genes with high similarity to their A. thaliana orthologues are briefly discussed. Detailed descriptions of these genes are also shown in Table 1. The A. thaliana knockout mutants for many of these genes developed by either T-DNA insertion, ethyl methanesulfonate (EMS) mutagenesis, or x-ray mutagenesis suggest roles in cell expansion and elongation. Some examples of the genes more abundantly transcribed in WT fibers during elongation included putative A. thaliana orthologues that when silenced by mutagenesis caused shortened petiole cells in herk1 theseus1 (the1) double mutants  and herk1 herk2 the1 triple mutants ; inhibition of root hair tip growth in mutants for the potassium transporter TINY ROOT HAIR 1 (TRH1) ; a 50% reduction in root hair length in mutants for a type XI myosin protein, XIK ; shortened leaf trichomes, reduced polarized expansion of hypocotyls, and shorter hypocotyl cells in mutants for the actin polymerization factor CROOKED (ARPC5) ; reduced length of root cells in mutants for the glycosylphosphatidylinositol (GPI)-anchored protein COBRA ; and reduced cell sizes resulting in dwarfism in the mutant brassinosteroid insensitive 1 (bri1) . The functionality of a cotton orthologue of BRI1 was previously demonstrated by complementation to rescue the WT phenotype of an A. thaliana bri1 mutant. Furthermore, chemical inhibition of a BR mediated response was demonstrated to inhibit fiber initiation and elongation in in vitro cotton ovule cultures and in planta [73, 74]. Two more G. hirsutum genes encoded the transcription factor GhMYB109 and the microtubule protein Beta tubulin1 (GhTub1) more abundantly in WT fibers compared to mutant fibers. Antisense suppression of GhMYB109 in transgenic cotton caused smaller and shrunken fibers initials and inhibition of fiber elongation resulting in a short-fiber phenotype [75, 76]. The GhTub1 gene was transcribed 46-fold and 7-fold higher in WT fibers at 8 and 12 DPA, respectively, and resulted in longitudinal cell growth when over-expressed in fission yeast cells under the control of an inducible promoter .
It is well-documented that increased levels of natural and/or synthetic cytokinins inhibit cotton fiber cell elongation in the in vitro cotton ovule culture system [78, 79]. Genes that were more abundantly transcribed in Li 2 mutant fibers compared to WT fibers encode proteins that are involved in cytokinin signaling and included putative orthologues of the A. thaliana HISTIDINE-CONTAINING PHOSPHOTRANSMITTER 1 (AtAHP1) that functions as a positive regulator of cytokinin signaling ; the A. thaliana homeobox domain transcription factor HAT22 that is up-regulated in response to exogenous application of cytokinins ; and orthologues of the A. thaliana UDP-glycosyltransferases UGT73C2 and UGT73C5/DON-GLUCOSYLTRANSFERASE 1 (AtDOGT1), glycosylate cytokinins which are implicated in cytokinin homeostasis . Over-expression of UGT73C5/DOGT1 in transgenic A. thaliana results in a dwarf phenotype with reduced elongation of hypocotyls similar to BR-deficient mutants . Increased cytokinin levels as determined indirectly by enzyme-linked immunosorbent assay measurements were previously reported in the ovules and fibers of four cotton fiber mutants compared to WT . This study included Ligon lintless cotton, but did not specify Li 1 or Li 2 . While immunoassays that measure plant hormone are generally semi-quantitative, the results do coincide with our gene expression results and suggest a possible role for cytokinin signaling and homeostasis in the Li 2 mutant phenotype (Table 1). Any one of these genes or gene families could plausibly explain the Li 2 short-fiber phenotype based on extensive functional characterization in cotton and heterologous systems such as A. thaliana. However, the fact that Li 2 is a monogenic dominant trait suggests that the differential expression of these genes in Li 2 mutant and WT fibers are downstream effects of a defect in a regulatory element controlling fiber development, specifically elongation in the case of the Li 2 mutation.
Little information is available regarding the cotton gene harboring the NAU3991 EST-SSR marker with complete linkage to the Li 2 locus. The function of this gene in plants is speculative and could be structural or regulatory in nature based on the putative DNA binding motif and homology to the mammalian plectin protein that acts as a microfilament crosslinker in mammalian systems . Presently, no functional information is available for this gene in any plant species and no phenotypes listed for mutants of A. thaliana or other plant species.
One of the major limitations of a gene expression analysis study is the lack of empirical evidence for functionality of any of the identified candidate genes that may cause an aberrant phenotype or enhance a trait of interest. This leads to a great deal of speculation as we have done here and in other cotton gene expression studies. This is especially true for a quantitative trait such as cotton fiber quality. The main goal of this gene expression study focused on practical concerns such as conversion of gene expression data into quantifiable, portable, and useable systems such as molecular markers. In regards to the regulation of a complex and multi-tiered event such as cotton fiber development by a single dominant gene, it should be expected that many downstream events from the developmental point of global regulation will be dramatically altered resulting in a high false discovery rate of candidate genes that truly are differentially expressed between, for example, mutant and WT fibers, but may be acting synergistically to cause the mutant phenotype. Since Li 2 is a monogenic dominant trait, it is possible that not one of these genes alone is the actual mutant allele causing the phenotype. The discovery of an EST-SSR marker based on gene expression data and with complete linkage to the Li 2 locus presents a novel approach for data mining a genetic marker database. Presently, we are cloning the full-length cDNAs and genomic DNAs for the NAU3991 marker gene from the Li 2 mutant and WT plants and have begun screening a G. hirsutum bacterial artificial chromosome (BAC) library. Identification of BAC clones containing the NAU3991 marker gene(s) will allow analysis of the genomic regions flanking the Li 2 locus and facilitate functional analysis of the Li 2 gene in transgenic cotton.
The cotton short fiber mutation Ligon lintless-2 (Li 2 ) appears specific to the elongation stage of fiber development with no apparent morphological differences between WT and mutant fibers until approximately 5 DPA. Comparative microarray gene expression profiling identified genes differentially expressed between WT and mutant fibers and suggested a role for ROS and cytokinin regulation that could result in the short fiber phenotype observed in the Li 2 mutant fibers. The microarray gene expression data was successfully converted into an EST-derived SSR marker, NAU3991, which displayed complete linkage to the Li 2 locus on chromosome 18. The complete linkage suggested that the gene harboring the NAU3991, or another unknown gene closely linked to the EST-SSR marker, may be the Li 2 gene. The function of the gene harboring the NAU3991 marker is unknown in plant species, but shares homology with a gene encoding a plectin protein that acts as a microfilament crosslinker in mammalian systems.
day of anthesis
expressed sequence tag
primary cell wall
restriction fragment length polymorphism
reactive oxygen species
reverse transcription quantitative polymerase chain reaction
secondary cell wall
scanning electron microscopy
single nucleotide polymorphism
single seed decent
simple sequence repeat
Basra AS, Malik CP: Development of the cotton fiber. Int Rev Cytol. 1984, 89: 65-113.
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: 1361-1366. 10.1104/pp.010724.
Meinert MC, Delmer DP: Changes in the biochemical composition of the cell wall of the cotton fiber during development. Plant Physiol. 1977, 59: 1088-1097. 10.1104/pp.59.6.1088.
Schubert AM, Benedict CR, Gates CE, Kohel RJ: Growth and development of the lint fibers of Pima S-4 cotton. Crop Sci. 1976, 16: 539-543. 10.2135/cropsci1976.0011183X001600040026x.
Thaker VS, Saroop S, Vaishnav PP, Singh YD: Genotypic variations and influence of diurnal temperature on cotton fiber development. Field Crop Res. 1989, 22: 129-141. 10.1016/0378-4290(89)90063-4.
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. 10.2135/cropsci2010.10.0569.
Haigler CH, Singh B, Wang G, Zhang D: From genomics of cotton fiber secondary cell wall deposition and cellulose biogenesis. Genetics and Genomics of Cotton. Edited by: Paterson AH. 2009, New York: Springer Science+Business Media, LLC, 3: 385-417. Jorgensen RA (Series Editor) Plant Genetics/Genomics
Gamble GR: Variation in surface chemical constituents of cotton (Gossypium hirsutum) fiber as a function of maturity. J Agr Food Chem. 2003, 51: 7995-7998. 10.1021/jf030178f.
Wilkins TA, Jernstedt JA: From molecular genetics of developing cotton fibers. Cotton Fibers: Developmental Biology, Quality Improvement, and Textile Processing. Edited by: Basra AS. 1999, New York: The Hawthorn Press, Inc., 231-269.
Griffee F, Ligon LL: Occurrence of "lintless" cotton plants and the inheritance of the character "lintless". J Am Soc Agron. 1929, 21: 711-717. 10.2134/agronj1929.00021962002100070001x.
Kearney TH, Harrison GJ: Inheritance of smooth seeds in cotton. J Agric Res. 1927, 35: 193-217.
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
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: 418-432. 10.1007/s00425-005-0098-7.
Percy RG, Kohel RJ: Quantitative genetics. Cotton: origin, history, technology, and production. Edited by: Smith CW. 1999, New York: John Wiley and Sons, Inc., 319-360.
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: 369-376. 10.1016/j.ygeno.2010.08.009.
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: 1475-1490. 10.1007/s00425-007-0580-5.
Turley RB, Kloth RH: The inheritance model for the fiberless trait in upland cotton (Gossypium hirsutum L.) line SL1-7-1: variation on a theme. Euphytica. 2008, 164: 123-132. 10.1007/s10681-008-9670-0.
Zhang TZ, Pan JJ: Genetic analysis of a fuzzless-lintless mutant in Gossypium hirsutum L. Jiangsu J Agr Sci. 1991, 7: 13-16.
Ware JO, Benedict LI, Rolfe WH: A recessive naked seed character in Upland cotton. J Hered. 1947, 38: 313-319.
Kohel RJ: Linkage tests in Upland cotton, Gossypium hirsutum L. II. Crop Sci. 1972, 12: 66-69. 10.2135/cropsci1972.0011183X001200010022x.
Narbuth EV, Kohel RJ: Inheritance and linkage analysis of a new fiber mutant in cotton. J Hered. 1990, 81: 131-133.
An C, Jenkins JN, McCarty JC, Saha S: Atypical Ligon lintless-2 phenotype in cotton. J Cotton Sci. 2010, 14: 13-16.
Karaca M, Saha S, Jenkins JN, Zipf A, Kohel R, Stelly DM: Simple sequence repeat (SSR) markers linked to the Ligon lintless (Li 1 ) mutant in cotton. J Hered. 2002, 93: 221-224. 10.1093/jhered/93.3.221.
Rong J, Pierce GJ, Waghmare VN, Rogers CJ, Desai A, Chee PW, May OL, Gannaway JR, Wendel JF, Wilkins TA, Paterson AH: Genetic mapping and comparative analysis of seven mutants related to seed fiber development in cotton. Theor Appl Genet. 2005, 111: 1137-1146. 10.1007/s00122-005-0041-0.
Kohel RJ, Stelly DM, Yu J: Tests of six cotton (Gossypium hirsutum L.) mutants for association with aneuploids. J Hered. 2002, 93: 130-132. 10.1093/jhered/93.2.130.
Lin L, Pierce GJ, Bowers JE, Estill JC, Compton RO, Rainville LK, Kim C, Lemke C, Rong J, Tang H, Wang X, Braidotti M, Chen AH, Chicola K, Collura K, Epps E, Golser W, Grover C, Ingles J, Karunakaran S, Kudrna D, Olive J, Tabassum N, Um E, Wissotski M, Yu Y, Zuccolo A, ur Rahman M, Peterson DG, Wing RA, Wendel JF, Paterson AH: A draft physical map of a D-genome cotton species (Gossypium raimondii). BMC Genomics. 2010, 11: 395-10.1186/1471-2164-11-395.
Liu F, Liang W, Guo W, Zhang T: Genetic mapping of short fiber mutant gene Li 2 in upland cotton. Cotton Sci. 2010, 22: 398-402.
Kohel RJ, Narbuth EV, Benedict CR: Fiber development of Ligon lintless-2 mutant of cotton. Crop Sci. 1992, 32: 733-735. 10.2135/cropsci1992.0011183X003200030031x.
Sedcole JR: Number of plants necessary to recover a trait. Crop Sci. 1977, 17: 677-678.
Paterson AH, Brubaker CL, Wendel JF: A rapid method for extraction of Cotton (Gossypium spp) genomic DNA suitable for RFLP or PCR analysis. Plant Mol Biol Rep. 1993, 11: 122-127. 10.1007/BF02670470.
Guo W, Cai C, Wang C, Han Z, Song X, Wang K, Niu X, Wang C, Lu K, Shi B, Zhang T: A microsatellite-based, gene-rich linkage map reveals genome structure, function and evolution in Gossypium. Genetics. 2007, 176: 527-541. 10.1534/genetics.107.070375.
Guo W, Cai C, Wang C, Zhao L, Wang L, Zhang T: A preliminary analysis of genome structure and composition in Gossypium hirsutum. BMC Genomics. 2008, 9: 314-10.1186/1471-2164-9-314.
Lacape JM, Nguyen TB, Thibivilliers S, Bojinov B, Courtois B, Cantrell RG, Burr B, Hau B: A combined RFLP-SSR-AFLP map of tetraploid cotton based on a Gossypium hirsutum x Gossypium barbadense backcross population. Genome. 2003, 46: 612-626. 10.1139/g03-050.
Nguyen TB, Giband M, Brottier P, Risterucci AM, Lacape JM: Wide coverage of the tetraploid cotton genome using newly developed microsatellite markers. Theor Appl Genet. 2004, 109: 167-175. 10.1007/s00122-004-1612-1.
Xiao J, Wu K, Fang DD, Stelly DM, Yu J, Cantrell RG: New SSR markers for use in cotton (Gossypium spp.) improvement. J Cotton Sci. 2009, 13: 75-157.
Michelmore RW, Paran I, Kesseli RV: Identification of markers linked to disease-resistance genes by bulked segregant analysis: a rapid method to detect markers in specific genomic regions by using segregating populations. Proc Natl Acad Sci USA. 1991, 88: 9828-9832. 10.1073/pnas.88.21.9828.
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: 943-953. 10.1007/s00122-009-1223-y.
Stam P: Construction of integrated genetic linkage maps by means of a new computer package: JoinMap. Plant J. 1993, 3: 739-744. 10.1111/j.1365-313X.1993.00739.x.
Van Ooijen JW, Voorrips RE: JoinMap 3.0, Software for the calculation of genetic linkage maps. Plant Research International. 2001, Wageningen, The Netherlands
Robinson G: The use of porous electron microscope capsules to prepare small lymph nodes for histology. Stain Technol. 1978, 53: 177-178.
Taliercio EW, Boykin D: Analysis of gene expression in cotton fiber initials. BMC Plant Biol. 2007, 7: 22-10.1186/1471-2229-7-22.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-
Whittaker DJ, Triplett BA: Gene-specific changes in alpha-tubulin transcript accumulation in developing cotton fibers. Plant Physiol. 1999, 121: 181-8. 10.1104/pp.121.1.181.
Gish W, States DJ: Identification of protein coding regions by database similarity search. Nature Genet. 1993, 3: 266-272. 10.1038/ng0393-266.
Harmer SE, Orford SJ, Timmis JN: Characterisation of six alpha-expansin genes in Gossypium hirsutum (upland cotton). Mol Genet Genomics. 2002, 268: 1-9. 10.1007/s00438-002-0721-2.
Kim HJ, Kato N, Kim S, Triplett B: CU/ZN superoxide dismutase in developing cotton fibers: evidence for an extracellular form. Planta. 2008, 228: 281-292. 10.1007/s00425-008-0734-0.
Pear JR, Kawagoe Y, Schreckengost WE, Delmer DP, Stalker DM: Higher plants contain homologs of the bacterial celA genes encoding the catalytic subunit of cellulose synthase. Proc Natl Acad Sci USA. 1996, 93: 12637-12642. 10.1073/pnas.93.22.12637.
Yang SS, Cheung F, Lee JJ, Ha M, Wei NE, Sze SH, Stelly DM, Thaxton P, Triplett B, Town CD, Jeffrey Chen Z: Accumulation of genome-specific transcripts, transcription factors and phytohormonal regulators during early stages of fiber cell development in allotetraploid cotton. Plant J. 2006, 47: 761-775. 10.1111/j.1365-313X.2006.02829.x.
Hinchliffe DJ, Meredith WR, Yeater KM, Kim HJ, Woodward AW, Chen ZJ, Triplett BA: Near-isogenic cotton germplasm lines that differ in fiber-bundle strength have temporal differences in fiber gene expression patterns as revealed by comparative high-throughput profiling. Theor Appl Genet. 2010, 120: 1347-1366. 10.1007/s00122-010-1260-6.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.
Benedito V, Torez-Jerez I, Murray J, Andriankaja A, Allen S, Kakar K, Wandrey M, Verdier J, Zuber H, Ott T, Moreau S, Niebel A, Frickey T, Weiller G, He J, Dai X, Zhao P, Tang Y, Udvardi M: A gene expression atlas of the model legume Medicago truncatula. Plant J. 2008, 55: 504-513. 10.1111/j.1365-313X.2008.03519.x.
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: 3674-3676. 10.1093/bioinformatics/bti610.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Yu Y, Yuan D, Liang S, Li X, Wang X, Lin Z, Zhang X: Genome structure of cotton revealed by a genome-wide SSR genetic map constructed from a BC1 population between Gossypium hirsutum and G. barbadense. BMC Genomics. 2011, 12: 15-10.1186/1471-2164-12-15.
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: D225-229. 10.1093/nar/gkq1189.
Kohel RJ, Benedict CR, Jividen GM: Incorporation of [14C] glucose into crystalline cellulose in aberrant fibers of a cotton mutant. Crop Sci. 1993, 33: 1036-1040. 10.2135/cropsci1993.0011183X003300050032x.
Hovav R, Udall JA, Chaudhary B, Hovav E, Flagel L, Hu G, Wendel JF: The evolution of spinnable cotton fiber entailed prolonged development and a novel metabolism. PLoS Genet. 2008, 4: e25-10.1371/journal.pgen.0040025.
Chaudhary B, Hovav R, Flagel L, Mittler R, Wendel JF: Parallel expression evolution of oxidative stress-related genes in fiber from wild and domesticated diploid and polyploid cotton (Gossypium). BMC Genomics. 2009, 10: 378-10.1186/1471-2164-10-378.
Chaudhary B, Hovav R, Rapp R, Verma N, Udall JA, Wendel JF: Global analysis of gene expression in cotton fibers from wild and domesticated Gossypium barbadense. Evol Dev. 2008, 10: 567-82. 10.1111/j.1525-142X.2008.00272.x.
Yamauchi Y, Hasegawa A, Taninaka A, Mizutani M, Sugimoto Y: NADPH-dependent reductases involved in the detoxification of reactive carbonyls in plants. J Biol Chem. 2011, 286: 6999-7009. 10.1074/jbc.M110.202226.
Li HB, Qin YM, Pang Y, Song WQ, Mei WQ, Zhu YX: A cotton ascorbate peroxidase is involved in hydrogen peroxide homeostasis during fibre cell development. New Phytol. 2007, 175: 462-471. 10.1111/j.1469-8137.2007.02120.x.
Mei W, Qin Y, Song W, Li J, Zhu Y: Cotton GhPOX1 encoding plant class III peroxidase may be responsible for the high level of reactive oxygen species production that is related to cotton fiber elongation. J Genet Genomics. 2009, 36: 141-150. 10.1016/S1673-8527(08)60101-0.
Arce DP, Godoy AV, Tsuda K, Yamazaki K, Valle EM, Iglesias MJ, Di Mauro MF, Casalongué CA: The analysis of an Arabidopsis triple knock-down mutant reveals functions for MBF1 genes under oxidative stress conditions. Plant Physiol. 2010, 167: 194-200. 10.1016/j.jplph.2009.09.003.
Suzuki N, Rizhsky L, Liang H, Shuman J, Shulaev V, Mittler R: Enhanced tolerance to environmental stress in transgenic plants expressing the transcriptional coactivator multiprotein bridging factor 1c. Plant Physiol. 2005, 139: 1313-1322. 10.1104/pp.105.070110.
Guo H, Li L, Ye H, Yu X, Algreen A, Yin Y: Three related receptor-like kinases are required for optimal cell elongation in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2009, 106: 7648-7653. 10.1073/pnas.0812346106.
Guo H, Ye H, Li L, Yin Y: A family of receptor-like kinases are regulated by BES1 and involved in plant growth in Arabidopsis thaliana. Plant Signal Behav. 2009, 4: 784-786. 10.4161/psb.4.8.9231.
Rigas S, Debrosses G, Haralampidis K, Vicente-Agullo F, Feldmann KA, Grabov A, Dolan L, Hatzopoulos P: TRH1 encodes a potassium transporter required for tip growth in Arabidopsis root hairs. Plant Cell. 2001, 13: 139-151.
Ojangu EL, Järve K, Paves H, Truve E: Arabidopsis thaliana myosin XIK is involved in root hair as well as trichome morphogenesis on stems and leaves. Protoplasma. 2007, 230: 193-202. 10.1007/s00709-006-0233-8.
Mathur J, Mathur N, Kirik V, Kernebeck B, Srinivas BP, Hülskamp M: Arabidopsis CROOKED encodes for the smallest subunit of the ARP2/3 complex and controls cell shape by region specific fine F-actin formation. Development. 2003, 130: 3137-3146. 10.1242/dev.00549.
Hauser MT, Morikami A, Benfey PN: Conditional root expansion mutants of Arabidopsis. Development. 1995, 121: 1237-1252.
Clouse SD, Langford M, McMorris TC: A brassinosteroid-insensitive mutant in Arabidopsis thaliana exhibits multiple defects in growth and development. Plant Physiol. 1996, 111: 671-678. 10.1104/pp.111.3.671.
Sun Y, Fokar M, Asami T, Yoshida S, Allen RD: Characterization of the brassinosteroid insensitive 1 genes of cotton. Plant Mol Biol. 2004, 54: 221-232.
Sun Y, Veerabomma S, Abdel-Mageed HA, Fokar M, Asami T, Yoshida S, Allen RD: Brassinosteroid regulates fiber development on cultured cotton ovules. Plant Cell Physiol. 2005, 46: 1384-1391. 10.1093/pcp/pci150.
Pu L, Li Q, Fan X, Yang W, Xue Y: The R2R3 MYB transcription factor GhMYB109 is required for cotton fiber development. Genetics. 2008, 180: 811-820. 10.1534/genetics.108.093070.
Suo J, Liang X, Pu L, Zhang Y, Xue Y: Identification of GhMYB109 encoding a R2R3 MYB transcription factor that expressed specifically in fiber initials and elongating fibers of cotton (Gossypium hirsutum L.). Biochim Biophys Acta. 2003, 1630: 25-34.
Ji S, Lu Y, Li J, Wei G, Liang X, Zhu Y: A beta-tubulin-like cDNA expressed specifically in elongating cotton fibers induces longitudinal growth of fission yeast. Biochem Biophys Res Commun. 2002, 296: 1245-50. 10.1016/S0006-291X(02)02069-7.
Beasley CA, Ting IP: Effects of plant growth substances on in vitro fiber development from fertilized cotton ovules. Amer J Bot. 1973, 60: 130-139. 10.2307/2441099.
Beasley CA, Ting IP: Effects of plant growth substances on in vitro fiber development from unfertilized cotton ovules. Amer J Bot. 1974, 61: 188-194. 10.2307/2441189.
Hutchison CE, Li J, Argueso C, Gonzalez M, Lee E, Lewis MW, Maxwell BB, Perdue TD, Schaller GE, Alonso JM, Ecker JR, Kieber JJ: The Arabidopsis histidine phosphotransfer proteins are redundant positive regulators of cytokinin signaling. Plant Cell. 2006, 18: 3073-3087. 10.1105/tpc.106.045674.
Brenner WG, Romanov GA, Köllmer I, Bürkle L, Schmülling T: Immediate-early and delayed cytokinin response genes of Arabidopsis thaliana identified by genome-wide expression profiling reveal novel cytokinin-sensitive processes and suggest cytokinin action through transcriptional cascades. Plant J. 2005, 44: 314-333. 10.1111/j.1365-313X.2005.02530.x.
Hou B, Lim EK, Higgins GS, Bowles DJ: N-glucosylation of cytokinins by glycosyltransferases of Arabidopsis thaliana. J Biol Chem. 2004, 279: 47822-47832. 10.1074/jbc.M409569200.
Poppenberger B, Fujioka S, Soeno K, George GL, Vaistij FE, Hiranuma S, Seto H, Takatsuto S, Adam G, Yoshida S, Bowles D: The UGT73C5 of Arabidopsis thaliana glucosylates brassinosteroids. Proc Natl Acad Sci USA. 2005, 102: 15253-15258. 10.1073/pnas.0504279102.
Chen J-G, Du X-M, Zhou X, Zhao H-Y: Levels of cytokinins in the ovules of cotton mutants with altered fiber development. J Plant Growth Regul. 1997, 16: 181-185. 10.1007/PL00006994.
Wiche G: Role of plectin in cytoskeleton organization and dynamics. J Cell Sci. 1998, 111: 2477-2486.
Li J, Chory J: A putative leucine-rich repeat receptor kinase involved in brassinosteroid signal transduction. Cell. 1997, 90: 929-938. 10.1016/S0092-8674(00)80357-8.
Belkhadir Y, Durbak A, Wierzba M, Schmitz RJ, Aguirre A, Michel R, Rowe S, Fujioka S, Tax FE: Intragenic suppression of a trafficking-defective brassinosteroid receptor mutant in Arabidopsis. Genetics. 2010, 185: 1283-1296. 10.1534/genetics.109.111898.
Benfey PN, Linstead PJ, Roberts K, Schiefelbein JW, Hauser MT, Aeschbacher RA: Root development in Arabidopsis: four mutants with dramatically altered root morphogenesis. Development. 1993, 119: 57-70.
Gao M, Wang X, Wang D, Xu F, Ding X, Zhang Z, Bi D, Cheng YT, Chen S, Li X, Zhang Y: Regulation of cell death and innate immunity by two receptor-like kinases in Arabidopsis. Cell Host Microbe. 2009, 6: 34-44. 10.1016/j.chom.2009.05.019.
This research was funded by United States Department of Agriculture-Agricultural Research Service project 6435-21000-015-00D. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.
DJH had the main responsibility for the study including the greenhouse experiment; tagging and sample harvest; RNA isolations; RT-qPCR and statistical analysis of the RT-qPCR data; gene selection from the microarray data for SSR data mining; and nucleotide sequence analysis. RBT developed the Li 2 mutant and WT NILs and the F2 mapping population used as plant material. MN performed RNA isolations and assessment of RNA quality; gene selection for corroboration of the microarray results; RT-qPCR and statistical analysis of the RT-qPCR data; and nucleotide sequence analysis. HJK prepared samples for SEM analysis. YT performed statistical analysis on the microarray data. KMY designed the greenhouse experiment for proper experimental design and statistical analysis of the data. PL assisted with the greenhouse experiment; molecular marker analysis; and conducted cDNA cloning for sequence analysis. DDF conceived the experiment, coordinated and supervised the research, conducted molecular marker analysis, association mapping, and linkage map construction. DJH and DDF wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Pedigree of the Li 2 mutant and WT NILs. The NILs utilized in this study were in the G. hirsutum cv. DP5690 genetic background and in the BC5 generation. (TIFF 1 MB)
Additional file 2: Primer pair sequences designed for marker analysis and RT-qPCR. The forward and reverse nucleotide primer sequences are shown along with the microarray probesets IDs, BLASTx sequence descriptions, and relevant accession numbers. (XLS 26 KB)
Additional file 3: Over-representation analysis (ORA) of genes differentially expressed in Li 2 WT and mutant fibers ≥ 2-fold. The original p-values for the ORA are shown along with the FDR-corrected and FWER-corrected p-values. The statistical significance cutoff for the ORA was an FDR-corrected p-value of 0.05. The TestSeqs are the microarray probesets IDs over-represented for the indicated biological process. (XLS 364 KB)
Additional file 4: Corroboration of the microarray gene expression data by RT-qPCR on Li 2 mutant and WT fibers. Gene expression profiles of twelve genes selected to verify the microarray gene expression data. Affymetrix probeset IDs and predicted gene products are shown on the graph titles for each gene. The DPA time-points that revealed a significant (≥ 2-fold; p-value < 0.05) difference in transcript abundance are indicated by an asterisk and the fold-change in transcript abundance is shown on the graphs above each indicated time-point. Error bars indicate standard deviation from 3 biological replicates. (TIFF 6 MB)
About this article
Cite this article
Hinchliffe, D.J., Turley, R.B., Naoumkina, M. et al. 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 12, 445 (2011). https://doi.org/10.1186/1471-2164-12-445