Functional genomics of fuzzless-lintless mutant of Gossypium hirsutum L. cv. MCU5 reveal key genes and pathways involved in cotton fibre initiation and elongation

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Genome-wide identification and expression analysis of GL2-interacting-repressor (GIR) genes during cotton fiber and fuzz development

GL2-interacting-repressor (GIR) family members may contribute to fiber/fuzz formation via a newly discovered unique pathway in Gossypium arboreum. There are similarities between cotton fiber development and the formation of trichomes and root hairs. The GL2-interacting-repressors (GIRs) are crucial regulators of root hair and trichome formation. The GaFzl gene, annotated as GaGIR1, is negatively associated with trichome development and fuzz initiation. However, there is relatively little available information regarding the other GIR genes in cotton, especially regarding their effects on cotton fiber development. In this study, 21 GIR family genes were identified in the diploid cotton species Gossypium arboreum; these genes were divided into three groups. The GIR genes were characterized in terms of their phylogenetic relationships, structures, chromosomal distribution and evolutionary dynamics. These GIR genes were revealed to be unequally distributed on 12 chromosomes in the diploid cotton genome, with no GIR gene detected on Ga06. The cis-acting elements in the promoter regions were predicted to be responsive to light, phytohormones, defense activities and stress. The transcriptomic data and qRT-PCR results revealed that most GIR genes were not differentially expressed between the wild-type control and the fuzzless mutant line. Moreover, 14 of 21 family genes were expressed at high levels, indicating these genes may play important roles during fiber development and fuzz formation. Furthermore, Ga01G0231 was predominantly expressed in root samples, suggestive of a role in root hair formation rather than in fuzz initiation and development. The results of this study have enhanced our understanding of the GIR genes and their potential utility for improving cotton fiber through breeding.

Comparative Metabolomics Analysis Reveals Sterols and Sphingolipids Play a Role in Cotton Fiber Cell Initiation

Cotton fiber is a seed trichome that protrudes from the outer epidermis of cotton ovule on the day of anthesis (0 day past anthesis, 0 DPA). The initial number and timing of fiber cells are closely related to fiber yield and quality. However, the mechanism underlying fiber initiation is still unclear. Here, we detected and compared the contents and compositions of sphingolipids and sterols in 0 DPA ovules of Xuzhou142 lintless-fuzzless mutants (Xufl) and Xinxiangxiaoji lintless-fuzzless mutants (Xinfl) and upland cotton wild-type Xuzhou142 (XuFL). Nine classes of sphingolipids and sixty-six sphingolipid molecular species were detected in wild-type and mutants. Compared with the wild type, the contents of Sphingosine-1-phosphate (S1P), Sphingosine (Sph), Glucosylceramide (GluCer), and Glycosyl-inositol-phospho-ceramides (GIPC) were decreased in the mutants, while the contents of Ceramide (Cer) were increased. Detail, the contents of two Cer molecular species, d18:1/22:0 and d18:1/24:0, and two Phyto-Cer molecular species, t18:0/22:0 and t18:0/h22:1 were significantly increased, while the contents of all GluCer and GIPC molecular species were decreased. Consistent with this result, the expression levels of seven genes involved in GluCer and GIPC synthesis were decreased in the mutants. Furthermore, exogenous application of a specific inhibitor of GluCer synthase, PDMP (1-phenyl-2-decanoylamino-3-morpholino-1-propanol), in ovule culture system, significantly inhibited the initiation of cotton fiber cells. In addition, five sterols and four sterol esters were detected in wild-type and mutant ovules. Compared with the wild type, the contents of total sterol were not significantly changed. While the contents of stigmasterol and campesterol were significantly increased, the contents of cholesterol were significantly decreased, and the contents of total sterol esters were significantly increased. In particular, the contents of campesterol esters and stigmasterol esters increased significantly in the two mutants. Consistently, the expression levels of some sterol synthase genes and sterol ester synthase genes were also changed in the two mutants. These results suggested that sphingolipids and sterols might have some roles in the initiation of fiber cells. Our results provided a novel insight into the regulatory mechanism of fiber cell initiation.

Mitochondria: Key Organelles Accelerating Cell Wall Material Accumulation in Juice Sacs of Pummelo (Citrus grandis L. Osbeck) Fruits during Postharvest Storage

Granulation is a physiological disorder of juice sacs in citrus fruits, which develops through secondary cell wall formation. However, the synergistic changes in the cytoplasm of juice sac cells remain largely unknown. This study investigated the dynamic ultrastructure of juice sacs of “Guanxi” pummelo fruits by transmission electron microscopy and determined their cell wall material, soluble sugar, and organic acid contents. The results showed that lignin and hemicellulose are accumulated in juice sacs isolated from dorsal vascular bundles, while lignin and cellulose contribute to the granulation of juice sacs isolated from septal vascular bundles. The significant differences in lignin, cellulose, and hemicellulose contents between the two types of juice sacs began to be observed at 30 days of storage. Fructose levels were elevated in juice sacs isolated from the dorsal vascular bundles from 10 to 60 days. Sucrose contents significantly decreased in juice sacs isolated from the septal vascular bundles from 30 to 60 days. Meanwhile glucose, citric acid, and malic acid contents exhibited no apparent changes in both types of juice sacs. Based on the comprehensive analysis of the ultrastructure of both types of juice sacs, it was clearly found that plasma membrane ruptures induce cell wall material synthesis in intracellular spaces; however, cell wall substance contents did not significantly increase until the number of mitochondria sharply increased. In particular, sucrose contents began to decrease significantly just after the mitochondria amount largely increased in juice sacs isolated from the septal vascular bundles, indicating that mitochondria play a key role in regulating carbon source sugar partitioning for cell wall component synthesis.

Genetic mapping and transcriptomic characterization of a new fuzzless-tufted cottonseed mutant

Fiber mutants are unique and valuable resources for understanding the genetic and molecular mechanisms controlling initiation and development of cotton fibers that are extremely elongated single epidermal cells protruding from the seed coat of cottonseeds. In this study, we reported a new fuzzless-tufted cotton mutant (Gossypium hirsutum) and showed that fuzzless-tufted near-isogenic lines (NILs) had similar agronomic traits and a higher ginning efficiency compared to their recurrent parents with normal fuzzy seeds. Genetic analysis revealed that the mutant phenotype is determined by a single incomplete dominant locus, designated N5. The mutation was fine mapped to an approximately 250-kb interval containing 33 annotated genes using a combination of bulked segregant sequencing, SNP chip genotyping, and fine mapping. Comparative transcriptomic analysis using 0-6 days post-anthesis (dpa) ovules from NILs segregating for the phenotypes of fuzzless-tufted (mutant) and normal fuzzy cottonseeds (wild-type) uncovered candidate genes responsible for the mutant phenotype. It also revealed that the flanking region of the N5 locus is enriched with differentially expressed genes (DEGs) between the mutant and wild-type. Several of those DEGs are members of the gene families with demonstrated roles in cell initiation and elongation, such as calcium-dependent protein kinase and expansin. The transcriptome landscape of the mutant was significantly reprogrammed in the 6 dpa ovules and, to a less extent, in the 0 dpa ovules, but not in the 2 and 4 dpa ovules. At both 0 and 6 dpa, the reprogrammed mutant transcriptome was mainly associated with cell wall modifications and transmembrane transportation, while transcription factor activity was significantly altered in the 6 dpa mutant ovules. These results imply a similar molecular basis for initiation of lint and fuzz fibers despite certain differences.

Genetic Identification and Transcriptome Analysis of Lintless and Fuzzless Traits in Gossypium arboreum L

Cotton fibres, as single cells arising from the seed coat, can be classified as lint and fuzz according to their final length. Gossypium arboreum is a cultivated diploid cotton species and a potential donor of the A subgenome of the more widely grown tetraploid cottons. In this study, we performed genetic studies on one lintless and seven fuzzless G. arboreum accessions. Through association and genetic linkage analyses, a recessive locus on Chr06 containing GaHD-1 was found to be the likely gene underlying the lintless trait. GaHD-1 carried a mutation at a splicing acceptor site that resulted in alternative splicing and a deletion of 247 amino acid from the protein. The regions containing GaGIR1 and GaMYB25-like were found to be associated with fuzz development in G. arboreum, with the former being the major contributor. Comparative transcriptome analyses using 0-5 days post-anthesis (dpa) ovules from lintless, fuzzless, and normal fuzzy seed G. arboreum accessions revealed gene modules and hub genes potentially important for lint and fuzz initiation and development. Three significant modules and 26 hub genes associated with lint fibre initiation were detected by weighted gene co-expression network analysis. Similar analyses identified three vital modules and 10 hub genes to be associated with fuzz development. The findings in this study contribute to understanding the complex molecular mechanism(s) regulating fibre initiation and development and indicate that G. arboreum may have fibre developmental pathways different from tetraploid cotton. It also provides candidate genes for further investigation into modifying fibre development in G. arboreum.

Map-based cloning of qBWT-c12 discovered brassinosteroid-mediated control of organ size in cotton

Assuring fiber yield stability is the primary objective for cotton breeders since the world population is on the rise, and the demand for cotton fiber is increasing every year. Thus, enhancing average cotton boll weight (BWT) could improve seed cotton production, and ultimately to increase cotton fiber yield. This study accomplished the map-based cloning of a novel boll weight regulating locus, qBWT-c12, in cotton. Bulk segregation analysis detected linked markers, aided in the detection of a stable BWT regulating locus, qBWT-c12, on Chr12 in a novel boll size mutant, BS41. Progeny evaluation confined the qBWT-c12 to a 0.89 cM interval between the AD-A12_07 and AD-FM_44 markers in recombinant derived F₃ and F₄ populations. Homology mapping detected a 40 bp insertion-deletion (InDel) site in the AD-FM_44 clone sequence situated +341 downstream of GhBRH1_A12, which showed complete linkage to the BWT phenotype. The suppressed expression of GhBRH1_A12 suggested its putative involvement during early boll development events in BS41. Although brassinosteroid (BR) biosynthesis and signaling pathway genes were up regulated in different tissues, but the organ growth was suppressed leading to dwarf plants, smaller leaves, and de-morphed smaller bolls in BS41. Thus, a disruption in the BR signal cascade is anticipated and could be related to lower GhBRH1_A12 expression in BS41.This study firstly reported the genetic dissection of boll size regulation of G. barbadense in G. hirsutum background using map-based cloning of a BWT regulating locus, qBWT-c12. Moreover, it also emphasized the putative role GhBRH1_A12 in regulating BR homeostasis and its potential to modulate plant growth and boll development in cotton.

Comparative Proteomic Analysis of Molecular Differences between Leaves of Wild-Type Upland Cotton and Its Fuzzless-Lintless Mutant

Fuzzless-lintless mutant (fl) ovules of upland cotton have been used to investigate cotton fiber development for decades. However, the molecular differences of green tissues between fl and wild-type (WT) cotton were barely reported. Here, we found that gossypol content, the most important secondary metabolite of cotton leaves, was higher in Gossypium hirsutum L. cv Xuzhou-142 (Xu142) WT than in fl. Then, we performed comparative proteomic analysis of the leaves from Xu142 WT and its fl. A total of 4506 proteins were identified, of which 103 and 164 appeared to be WT- and fl-specific, respectively. In the 4239 common-expressed proteins, 80 and 74 were preferentially accumulated in WT and fl, respectively. Pathway enrichment analysis and protein–protein interaction network analysis of both variety-specific and differential abundant proteins showed that secondary metabolism and chloroplast-related pathways were significantly enriched. Quantitative real-time PCR confirmed that the expression levels of 12 out of 16 selected genes from representative pathways were consistent with their protein accumulation patterns. Further analyses showed that the content of chlorophyll a in WT, but not chlorophyll b, was significantly increased compared to fl. This work provides the leaf proteome profiles of Xu142 and its fl mutant, indicating the necessity of further investigation of molecular differences between WT and fl leaves.

Genome-wide analysis of cotton C2H2-zinc finger transcription factor family and their expression analysis during fiber development

BACKGROUND

C2H2-zinc finger protein family is commonly found in the plant, and it is known as the key actors in the regulation of transcription and vital component of chromatin structure. A large number of the C2H2-zinc finger gene members have not been well characterized based on their functions and structure in cotton. However, in other plants, only a few C2H2-zinc finger genes have been studied.

RESULTS

In this work, we performed a comprehensive analysis and identified 386, 196 and 195 C2H2-zinc finger genes in Gossypium hirsutum (upland cotton), Gossypium arboreum and Gossypium raimondii, respectively. Phylogenetic tree analysis of the C2H2-zinc finger proteins encoding the C2H2-zinc finger genes were classified into seven (7) subgroups. Moreover, the C2H2-zinc finger gene members were distributed in all cotton chromosomes though with asymmetrical distribution patterns. All the orthologous genes were detected between tetraploid and the diploid cotton, with 154 orthologous genes pair detected between upland cotton and Gossypium arboreum while 165 orthologous genes were found between upland cotton and Gossypium raimondii. Synonymous (Ks) and non-synonymous (Ka) nucleotide substitution rates (Ka/Ks) analysis indicated that the cotton C2H2-zinc finger genes were highly influenced mainly by negative selection, which maintained their protein levels after the duplication events. RNA-seq data and RT-qPCR validation of the RNA seq result revealed differential expression pattern of some the C2H2-zinc finger genes at different stages of cotton fiber development, an indication that the C2H2-zinc finger genes play an important role in initiating and regulating fiber development in cotton.

CONCLUSIONS

This study provides a strong foundation for future practical genome research on C2H2-zinc finger genes in upland cotton. The expression levels of C2H2-zinc finger genes family is a pointer of their involvement in various biochemical and physiological functions which are directly related to cotton fiber development during initiation and elongation stages. This work not only provides a basis for determining the nominal role of the C2H2-zinc finger genes in fiber development but also provide valuable information for characterization of potential candidate genes involved in regulation of cotton fiber development.

Genome-wide classification, evolutionary analysis and gene expression patterns of the kinome in Gossypium

The protein kinase (PK, kinome) family is one of the largest families in plants and regulates almost all aspects of plant processes, including plant development and stress responses. Despite their important functions, comprehensive functional classification, evolutionary analysis and expression patterns of the cotton PK gene family has yet to be performed on PK genes. In this study, we identified the cotton kinomes in the Gossypium raimondii, Gossypium arboretum, Gossypium hirsutum and Gossypium barbadense genomes and classified them into 7 groups and 122-24 subfamilies using software HMMER v3.0 scanning and neighbor-joining (NJ) phylogenetic analysis. Some conserved exon-intron structures were identified not only in cotton species but also in primitive plants, ferns and moss, suggesting the significant function and ancient origination of these PK genes. Collinearity analysis revealed that 16.6 million years ago (Mya) cotton-specific whole genome duplication (WGD) events may have played a partial role in the expansion of the cotton kinomes, whereas tandem duplication (TD) events mainly contributed to the expansion of the cotton RLK group. Synteny analysis revealed that tetraploidization of G. hirsutum and G. barbadense contributed to the expansion of G. hirsutum and G. barbadense PKs. Global expression analysis of cotton PKs revealed stress-specific and fiber development-related expression patterns, suggesting that many cotton PKs might be involved in the regulation of the stress response and fiber development processes. This study provides foundational information for further studies on the evolution and molecular function of cotton PKs.

Inhibition of Heat Shock proteins HSP90 and HSP70 induce oxidative stress, suppressing cotton fiber development

Cotton fiber is a specialized unicellular structure useful for the study of cellular differentiation and development. Heat shock proteins (HSPs) have been shown to be involved in various developmental processes. Microarray data analysis of five Gossypium hirsutum genotypes revealed high transcript levels of GhHSP90 and GhHSP70 genes at different stages of fiber development, indicating their importance in the process. Further, we identified 26 and 55 members of HSP90 and HSP70 gene families in G. hirsutum. The treatment of specific inhibitors novobiocin (Nov; HSP90) and pifithrin/2-phenylethynesulfonamide (Pif; HSP70) in in-vitro cultured ovules resulted in a fewer number of fiber initials and retardation in fiber elongation. The molecular chaperone assay using bacterially expressed recombinant GhHSP90-7 and GhHSP70-8 proteins further confirmed the specificity of inhibitors. HSP inhibition disturbs the H₂O₂ balance that leads to the generation of oxidative stress, which consequently results in autophagy in the epidermal layer of the cotton ovule. Transmission electron microscopy (TEM) of inhibitor-treated ovule also corroborates autophagosome formation along with disrupted mitochondrial cristae. The perturbations in transcript profile of HSP inhibited ovules show differential regulation of different stress and fiber development-related genes and pathways. Altogether, our results indicate that HSP90 and HSP70 families play a crucial role in cotton fiber differentiation and development by maintaining cellular homeostasis.

Genetic dissection of the fuzzless seed trait in Gossypium barbadense

Cotton fibres are single-celled trichomes arising from the epidermal cells of the seed coat and may be either long (lint) or very short (fuzz). The dominant fuzzless N1 of Gossypium hirsutum is a defective allele of the At-subgenome homoeolog of MYB25-like, but the genetic components underlying the recessive fuzzless trait from G. barbadense (Gb) are unknown. We have identified five genetic loci, including a major contributing locus containing MYB25-like_Dt, associated with Gb fuzzless seeds based on genotyping of fuzzy and fuzzless near isogenic lines (NILs) from an interspecies cross (G. barbadense × G. hirsutum). At 3 d post-anthesis when fuzz fibres are initiating, expression of MYB25-like_Dt was significantly lower in fuzzless NILs than in fuzzy seeded NILs, while higher MYB25-like_Dt expression was associated with more seed fuzz across different cotton genotypes. Phenotypic and genotypic analysis of MYB25-like homoeoalleles in cottons showing different fibre phenotypes and their crossing progeny indicated that both MYB25-like_At and MYB25-like_Dt are associated with lint development, and that fuzz development is mainly determined by the expression level of MYB25-like_Dt at ~3 d post-anthesis. Expression of Gb fuzzless seeds depends on genetic background and interactions amongst the multiple loci identified. MYB25-like_Dt is one of the best candidates for N2.

Transcriptomic profiling of developing fiber in levant cotton (Gossypium herbaceum L.)

Cotton (Gossypium spp.) is an imperative economic crop of the globe due to its natural textile fiber. Molecular mechanisms of fiber development have been greatly revealed in allotetraploid cotton but remained unexplored in Gossypium herbaceum. G. herbaceum can withstand the rigors of nature like drought and pests but produce coarse lint. This undesirable characteristic strongly needs the knowledge of fiber development at molecular basis. The present study reported the transcriptome sequence of the developing fiber of G. herbaceum on pyrosequencing and its analysis. About 1.38 million raw and 1.12 million quality trimmed reads were obtained followed by de novo assembly-generated 20,125 unigenes containing 14,882 coding sequences (CDs). BLASTx-based test of homology indicated that A1-derived transcripts shared a high similarity with Gossypium arboreum (A2). Functional annotation of the CDs using the UniProt categorized them into biological processes, cellular components, and molecular function, COG classification showed that a large number of CDs have significant homology in COG database (6215 CDs), and mapping of CDs with Kyoto Encyclopedia of Genes and Genomes (KEGG) database generated 200 pathways ultimately showing predominant engagement in the fiber development process. Transcription factors were predicted by comparison with Plant Transcription Factor Database, and their differential expression between stages exposed their important regulatory role in fiber development. Differential expression analysis based on reads per kilobase of transcript per million mapped reads (RPKM) value revealed activities of specific gene related to carbohydrate and lipid synthesis, carbon metabolism, energy metabolism, signal transduction, etc., at four stages of fiber development, and was validated by qPCR. Overall, this study will help as a valuable foundation for diploid cotton fiber improvement.

Transcript profiling of genes expressed during fibre development in diploid cotton (Gossypium arboreum L.)

Background

Cotton fibre is a single cell and it is one of the best platforms for unraveling the genes express during various stages of fibre development. There are reports devoted to comparative transcriptome study on fiber cell initiation and elongation in tetraploid cultivated cotton. However, in the present investigation, comparative transcriptome study was made in diploid cultivated cotton using isogenic fuzzy-lintless (Fl) and normal fuzzy linted (FL) lines belong to Gossypium arboreum, diploid species at two stages, 0 and 10 dpa (days post anthesis), using Affymetrix cotton GeneChip genome array.

Result

Scanning electron microscopy (SEM) analysis uncovered the occurrence of few fibre cell initials in the Fl line as compared to many in Normal FL at −2 and 0 dpa. However, at 10 dpa there were no fibre cells found elongated in Fl but many elongated cells were found in FL line. Up-regulation of transcription factors, AP2-EREBP, C2H2, C3H, HB and WRKY was observed at 0 dpa whereas in 10 dpa transcription factors, AP2-EREBP, AUX/IAA, bHLH, C2H2, C3H, HB, MYB, NAC, Orphans, PLATZ and WRKY were found down regulated in Fl line. These transcription factors were mainly involved in metabolic pathways such as phytohormone signaling, energy metabolism of cell, fatty acid metabolism, secondary metabolism and other signaling pathways and are related directly or indirectly in fiber development. Quantitative real-time PCR was performed to check fold up or down-regulation of these genes and transcription factors (TFs) down regulated in mutants as compared to normal at 0 and 10 dpa.

Conclusion

This study elucidates that the up-regulation of transcription factors like AP2-EREBP, C2H2, C3H, HB, WRKY and phytohormone signaling genes at 0 dpa and their down-regulation at the 10 dpa might have constrain the fibre elongation in fuzzy-lintless line. Along with this the down-regulation of genes involved in synthesis of VLCFA chain, transcripts necessary for energy and cell wall metabolism, EXPANSINs, arabinogalactan proteins (AGPs), tubulin might also be the probable reason for reduced growth of fibres in the Fl. Plant receptor-like kinases (RLKs), Leucine Rich Repeats) LRR- family protein and signal transduction coding for mitogen-activated protein kinase (MAPK) cascade, have been engaged in coordination of cell elongation and SCW biosynthesis, down-regulation of these might loss the function leads to reduced fibre growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4066-y) contains supplementary material, which is available to authorized users.

In vivo imaging of Ca2+ accumulation during cotton fiber initiation using fluorescent indicator YC3.60

Non-tip-focused Ca ²⁺ gradient indicated by genetically expressing a FRET-based calcium sensor YC3.60 was established in spherical expanding cotton fibers, which is vital for cotton fiber initiation. Cotton fiber is a single cell elongated from ovule epidermis. It is not only the most important natural fiber used in the textile industry but also an ideal model for studying cell differentiation and elongation. Before linear cell growth, cotton fibers undergo spherical expansion at the beginning of initiation. Ca²⁺, as an important secondary messenger, plays a central role in polarized cell growth including cotton fiber elongation. However, the role of Ca²⁺ in fiber initiation is far from well understood. In this paper, through ovule culture we demonstrate that Ca²⁺ is crucial for fiber initiation. Using transgenic cotton expressing the fluorescent Ca²⁺ indicator YC3.60, we show cellular and intracellular distribution of Ca²⁺ in cotton ovule epidermis and fiber cells. In the initiating fiber cell, Ca²⁺ accumulated mainly at the base of the cell, while in the fast elongating cell, the Ca²⁺ was enriched in the tip region. This cellular distribution of Ca²⁺ reported by YC3.60 was confirmed by the staining with a Ca²⁺-sensitive dye fluo-3/AM. Compared to the fluorescent dye staining, the YC3.60 system can reveal more detailed information on the intracellular distribution without photobleaching. Taken together, our data suggest that Ca²⁺ plays an important role in spherical expansion of cotton fiber initials.

Identification of candidate genes for fiber length quantitative trait loci through RNA-Seq and linkage and physical mapping in cotton

BACKGROUND

Cotton (Gossypium spp.) fibers are single-celled elongated trichomes, the molecular aspects of genetic variation in fiber length (FL) among genotypes are currently unknown. In this study, two backcross inbred lines (BILs), i.e., NMGA-062 ("Long") and NMGA-105 ("Short") with 32.1 vs. 27.2 mm in FL, respectively, were chosen to perform RNA-Seq on developing fibers at 10 days post anthesis (DPA). The two BILs differed in 4 quantitative trait loci (QTL) for FL and were developed from backcrosses between G. hirsutum as the recurrent parent and G. barbadense.

RESULTS

In total, 51.7 and 54.3 million reads were obtained and assembled to 49,508 and 49,448 transcripts in the two genotypes, respectively. Of 1551 differentially expressed genes (DEGs) between the two BILs, 678 were up-regulated and 873 down-regulated in "Long"; and 703 SNPs were identified in 339 DEGs. Further physical mapping showed that 8 DEGs were co-localized with the 4 FL QTL identified in the BIL population containing the two BILs. Four SNP markers in 3 DEGs that showed significant correlations with FL were developed. Among the three candidate genes encoding for proline-rich protein, D-cysteine desulfhydrase, and thaumatin-like protein, a SNP of thaumatin-like protein gene showed consistent correlations with FL across all testing environments.

CONCLUSIONS

This study represents one of the first investigations of positional candidate gene approach of QTL in cotton in integrating transcriptome and SNP identification based on RNA-Seq with linkage and physical mapping of QTL and genes, which will facilitate eventual cloning and identification of genes responsible for FL QTL. The candidate genes may serve as the foundation for further in-depth studies of the molecular mechanism of natural variation in fiber elongation.

GhCaM7-like, a calcium sensor gene, influences cotton fiber elongation and biomass production

Calcium signaling regulates many developmental processes in plants. Calmodulin (CaM) is one of the most conserved calcium sensors and has a flexible conformation in eukaryotes. The molecular functions of CaM are unknown in cotton, which is a major source of natural fiber. In this study, a Gossypium hirsutum L.CaM7-like gene was isolated from upland cotton. Bioinformatics analysis indicated that the GhCaM7-like gene was highly conserved as compared with Arabidopsis AtCaM7. The GhCaM7-like gene showed a high expression level in elongating fibers. Expression of β-glucuronidase was observed in trichomes on the stem, leaf and root in transgenic Arabidopsis plants of a PRO_GhCaM7-like:GUS fusion. Silencing of the GhCaM7-like gene resulted in decreased fiber length, but also caused reduction in stem height, leaf dimensions, seed length and 100-seed weight, in comparison with those of the control. Reduced expression of the GhCaM7-like gene caused decreased Ca²⁺ influx in cells of the leaf hypodermis and stem apex, and down-regulation of GhIQD1 (IQ67-domain containing protein), GhAnn2 (Annexins) and GhEXP2 (Expansin). These results indicate that the GhCaM7-like gene plays a vital role in calcium signaling pathways, and may regulate cotton fiber elongation and biomass production by affecting Ca²⁺ signatures and downstream signaling pathways of CaM.

Comparative Proteomic Analysis of Cotton Fiber Development and Protein Extraction Method Comparison in Late Stage Fibers

The distinct stages of cotton fiber development and maturation serve as a single-celled model for studying the molecular mechanisms of plant cell elongation, cell wall development and cellulose biosynthesis. However, this model system of plant cell development is compromised for proteomic studies due to a lack of an efficient protein extraction method during the later stages of fiber development, because of a recalcitrant cell wall and the presence of abundant phenolic compounds. Here, we compared the quality and quantities of proteins extracted from 25 dpa (days post anthesis) fiber with multiple protein extraction methods and present a comprehensive quantitative proteomic study of fiber development from 10 dpa to 25 dpa. Comparative analysis using a label-free quantification method revealed 287 differentially-expressed proteins in the 10 dpa to 25 dpa fiber developmental period. Proteins involved in cell wall metabolism and regulation, cytoskeleton development and carbohydrate metabolism among other functional categories in four fiber developmental stages were identified. Our studies provide protocols for protein extraction from maturing fiber tissues for mass spectrometry analysis and expand knowledge of the proteomic profile of cotton fiber development.

Phytosterol content and the campesterol:sitosterol ratio influence cotton fiber development: role of phytosterols in cell elongation

Phytosterols play an important role in plant growth and development, including cell division, cell elongation, embryogenesis, cellulose biosynthesis, and cell wall formation. Cotton fiber, which undergoes synchronous cell elongation and a large amount of cellulose synthesis, is an ideal model for the study of plant cell elongation and cell wall biogenesis. The role of phytosterols in fiber growth was investigated by treating the fibers with tridemorph, a sterol biosynthetic inhibitor. The inhibition of phytosterol biosynthesis resulted in an apparent suppression of fiber elongation in vitro or in planta. The determination of phytosterol quantity indicated that sitosterol and campesterol were the major phytosterols in cotton fibers; moreover, higher concentrations of these phytosterols were observed during the period of rapid elongation of fibers. Furthermore, the decrease and increase in campesterol:sitosterol ratio was associated with the increase and decease in speed of elongation, respectively, during the elongation stage. The increase in the ratio was associated with the transition from cell elongation to secondary cell wall synthesis. In addition, a number of phytosterol biosynthetic genes were down-regulated in the short fibers of ligon lintless-1 mutant, compared to its near-isogenic wild-type TM-1. These results demonstrated that phytosterols play a crucial role in cotton fiber development, and particularly in fiber elongation.

Antisense expression of Gossypium barbadense UGD6 in Arabidopsis thaliana significantly alters cell wall composition

Uridine diphosphate-glucose dehydrogenase (UGD, EC1.1.1.22 oxidizes UDP-Glc (UDP-D-glucose) to UDP-GlcA (UDP-D-glucuronate), a critical precursor of cell wall polysaccharides. GbUGD6 from Gossypium barbadense is more highly expressed late in the elongation of cotton fibers (15 d post-anthesis (DPA)) and during the stage of secondary cell wall thickening (30 DPA). Subcellular localization analysis in onion epidermis revealed that fluorescently labeled GbUGD6 protein was distributed throughout the cell membrane, as well as the nucleus and vacuoles. Examination of UGD function in Arabidopsis revealed that the antisense GbUGD6 lines had shorter roots, deferred blossoming, compared to wild-type plants. Activities of associated enzymes were also affected by UGD reduction, and biochemical analysis of cell wall samples showed an increase in cellulose levels and a decrease in UGP-GlcA contents. The results of the present study as well as previous studies on UGD support the conclusion that UGD plays a major role in synthesizing polysaccharides synthesis in the cell wall.

Ethylene Biosynthesis Is Promoted by Very-Long-Chain Fatty Acids during Lysigenous Aerenchyma Formation in Rice Roots

In rice (Oryza sativa) roots, lysigenous aerenchyma, which is created by programmed cell death and lysis of cortical cells, is constitutively formed under aerobic conditions, and its formation is further induced under oxygen-deficient conditions. Ethylene is involved in the induction of aerenchyma formation. reduced culm number1 (rcn1) is a rice mutant in which the gene encoding the ATP-binding cassette transporter RCN1/OsABCG5 is defective. Here, we report that the induction of aerenchyma formation was reduced in roots of rcn1 grown in stagnant deoxygenated nutrient solution (i.e. under stagnant conditions, which mimic oxygen-deficient conditions in waterlogged soils). 1-Aminocyclopropane-1-carboxylic acid synthase (ACS) is a key enzyme in ethylene biosynthesis. Stagnant conditions hardly induced the expression of ACS1 in rcn1 roots, resulting in low ethylene production in the roots. Accumulation of saturated very-long-chain fatty acids (VLCFAs) of 24, 26, and 28 carbons was reduced in rcn1 roots. Exogenously supplied VLCFA (26 carbons) increased the expression level of ACS1 and induced aerenchyma formation in rcn1 roots. Moreover, in rice lines in which the gene encoding a fatty acid elongase, CUT1-LIKE (CUT1L; a homolog of the gene encoding Arabidopsis CUT1, which is required for cuticular wax production), was silenced, both ACS1 expression and aerenchyma formation were reduced. Interestingly, the expression of ACS1, CUT1L, and RCN1/OsABCG5 was induced predominantly in the outer part of roots under stagnant conditions. These results suggest that, in rice under oxygen-deficient conditions, VLCFAs increase ethylene production by promoting 1-aminocyclopropane-1-carboxylic acid biosynthesis in the outer part of roots, which, in turn, induces aerenchyma formation in the root cortex.

Water Stress and Foliar Boron Application Altered Cell Wall Boron and Seed Nutrition in Near-Isogenic Cotton Lines Expressing Fuzzy and Fuzzless Seed Phenotypes

Our previous research, conducted under well-watered conditions without fertilizer application, showed that fuzziness cottonseed trait resulted in cottonseed nutrition differences between fuzzy (F) and fuzzless (N) cottonseed. Under water stress conditions, B mobility is further limited, inhibiting B movement within the plant, affecting seed nutrition (quality). Therefore, we hypothesized that both foliar B and water stress can affect B mobility, altering cottonseed protein, oil, and mineral nutrition. The objective of the current research was to evaluate the effects of the fuzziness seed trait on boron (B) and seed nutrition under water stress and foliar B application using near-isogenic cotton lines (NILs) grown in a repeated greenhouse experiment. Plants were grown under-well watered conditions (The soil water potential was kept between -15 to -20 kPa, considered field capacity) and water stress conditions (soil water potential between -100 and -150 kPa, stressed conditions). Foliar B was applied at a rate of 1.8 kg B ha(-1) as H3BO3. Under well-watered conditions without B the concentrations of seed oil in N lines were higher than in F lines, and seed K and N levels were lower in N lines than in F lines. Concentrations of K, N, and B in leaves were higher in N lines than in F lines, opposing the trend in seeds. Water-stress resulted in higher seed protein concentrations, and the contribution of cell wall (structural) B to the total B exceeded 90%, supporting the structural role of B in plants. Foliar B application under well-watered conditions resulted in higher seed protein, oil, C, N, and B in only some lines. This research showed that cottonseed nutrition differences can occur due to seed fuzziness trait, and water stress and foliar B application can alter cottonseed nutrition.

Comparative Transcriptomics Reveals Jasmonic Acid-Associated Metabolism Related to Cotton Fiber Initiation

Analysis of mutants and gene expression patterns provides a powerful approach for investigating genes involved in key stages of plant fiber development. In this study, lintless-fuzzless XinWX and linted-fuzzless XinFLM with a single genetic locus difference for lint were used to identify differentially expressed genes. Scanning electron microscopy showed fiber initiation in XinFLM at 0 days post anthesis (DPA). Fiber transcriptional profiling of the lines at three initiation developmental stages (-1, 0, 1 DPA) was performed using an oligonucleotide microarray. Loop comparisons of the differentially expressed genes within and between the lines was carried out, and functional classification and enrichment analysis showed that gene expression patterns during fiber initiation were heavily associated with hormone metabolism, transcription factor regulation, lipid transport, and asparagine biosynthetic processes, as previously reported. Further, four members of the allene-oxide cyclase (AOC) family that function in jasmonate biosynthesis were parallel up-regulation in fiber initiation, especially at -1 DPA, compared to other tissues and organs in linted-fuzzed TM-1. Real time-quantitative PCR (RT-qPCR) analysis in different fiber mutant lines revealed that AOCs were up-regulated higher at -1 DPA in lintless-fuzzless than that in linted-fuzzless and linted-fuzzed materials, and transcription of the AOCs was increased under jasmonic acid (JA) treatment. Expression analysis of JA biosynthesis-associated genes between XinWX and XinFLM showed that they were up-regulated during fiber initiation in the fuzzless-lintless mutant. Taken together, jasmonic acid-associated metabolism was related to cotton fiber initiation. Parallel up-regulation of AOCs expression may be important for normal fiber initiation development, while overproduction of AOCs might disrupt normal fiber development.

Cottonseed protein, oil, and mineral status in near-isogenic Gossypium hirsutum cotton lines expressing fuzzy/linted and fuzzless/linted seed phenotypes under field conditions

Cotton is an important crop in the world and is a major source of oil for human consumption and cotton meal for livestock. Cottonseed nutrition (seed composition: protein, oil, and minerals) determines the quality of seeds. Therefore, maintaining optimum levels of cottonseed nutrition is critical. Physiological and genetic mechanisms controlling the levels of these constituents in cottonseed are still largely unknown. Our previous research conducted under greenhouse conditions showed that seed and leaf nutrition differed between fuzzless and fuzzy seed isolines. Therefore, the objective of this research was to investigate the seed fuzz phenotype (trait) effects on seed protein, oil, N, C, S, and minerals in five sets of near-isogenic mutant cotton lines for seed fuzz in a 2-year experiment under field condition to evaluate the stability of the effect of the trait on seed nutrition. The isolines (genotypes) in each set differ for the seed fuzz trait (fuzzless/linted seed line, N lines, and fuzzy/linted seed line, F lines). Results showed that seed protein was higher in the fuzzy genotype in all sets, but seed oil was higher in fuzzless genotype in all sets. The concentrations of seed Ca and C were higher in all fuzzless genotypes, but N, S, B, Fe, and Zn were higher in most of the fuzzy genotypes. Generally, minerals were higher in leaves of F lines, suggesting the translocation of minerals from leaves to seeds was limited. The research demonstrated that fiber development could be involved in cottonseed composition. This may be due to the involvement of fiber development in carbon and nitrogen metabolism, and the mobility of nutrients from leaves (source) to seed (sink). This information is beneficial to breeders to consider fuzzless cottonseed for potential protein and oil use and select for higher oil or higher protein content, and to physiologists to further understand the mobility of minerals to increase the quality of cottonseed nutrition for food and feed.

Quantitative proteomics and transcriptomics reveal key metabolic processes associated with cotton fiber initiation

An iTRAQ-based proteomics of ovules from the upland cotton species Gossypium hirsutum and its fuzzless-lintless mutant was performed, and finally 2729 proteins that preferentially accumulated at anthesis in wild-type ovules were identified. We confirmed that the gene expression levels of 2005 among these proteins also increased by performing an RNA sequencing transcriptomics. Expression of proteins involved in carboxylic acid metabolism, small-molecule metabolic processes, hormone regulation, and lipid metabolism was significantly enhanced in wild-type ovules. Quantitative real-time PCR verified the increased expression of 26 genes involved in these processes. Cotton 3-hydroxyacyl-CoA dehydratase (GhPAS2) catalyzing the third reaction of very long-chain fatty acid (VLCFA) biosynthesis, accumulated at anthesis in wild-type ovules. Heterogeneous expression of GhPAS2 restored viability to the Saccharomyces cerevisiae haploid psh1-deletion strain deficient in PAS2 activity. Application of VLCFA biosynthesis inhibitor acetochlor (2-chloro-N-[ethoxymethyl]-N-[2-ethyl-6-methyl-phenyl]-acetamide; ACE) and gibberellic acid to the unfertilized cotton ovules significantly suppressed fiber cell protrusion. In this study, the profiling of gene expression at both transcriptome and proteome levels provides new insights into cotton fiber cell initiation.

BIOLOGICAL SIGNIFICANCE

Cotton fiber initiation determines the ultimate number of fibers per ovule, thereby determining fiber yield. In total, 2729 proteins were preferentially accumulated in wild-type ovules at anthesis. The most up-regulated proteins were assigned to carboxylic acid metabolism, small-molecule metabolic processes, hormone regulation, and lipid metabolism. In consistence with these findings, we characterized GhPAS2 gene coding for the enzyme that catalyzes VLCFA production. VLCFA biosynthesis inhibitor, acetochlor, was shown to significantly suppress fiber initiation. This study provides a genome-scale transcriptomic and proteomic characterization of fiber initial cells, laying a solid basis for further investigation of the molecular processes governing fiber cell development.

Cotton fiber elongation network revealed by expression profiling of longer fiber lines introgressed with different Gossypium barbadense chromosome segments

Background

Cotton fiber, a highly elongated, thickened single cell of the seed epidermis, is a powerful cell wall research model. Fiber length, largely determined during the elongation stage, is a key property of fiber quality. Several studies using expressed sequence tags and microarray analysis have identified transcripts that accumulate preferentially during fiber elongation. To further show the mechanism of fiber elongation, we used Digital Gene Expression Tag Profiling to compare transcriptome data from longer fiber chromosome introgressed lines (CSILs) containing segments of various Gossypium barbadense chromosomes with data from its recurrent parent TM-1 during fiber elongation (from 5 DPA to 20 DPA).

Results

A large number of differentially expressed genes (DEGs) involved in carbohydrate, fatty acid and secondary metabolism, particularly cell wall biosynthesis, were highly upregulated during the fiber elongation stage, as determined by functional enrichment and pathway analysis. Furthermore, DEGs related to hormone responses and transcription factors showed upregulated expression levels in the CSILs. Moreover, metabolic and regulatory network analysis indicated that the same pathways were differentially altered, and distinct pathways exhibited altered gene expression, in the CSILs. Interestingly, mining of upregulated DEGs in the introgressed segments of these CSILs based on D-genome sequence data showed that these lines were enriched in glucuronosyltransferase, inositol-1, 4, 5-trisphosphate 3-kinase and desulfoglucosinolate sulfotransferase activity. These results were similar to the results of transcriptome analysis.

Conclusions

This report provides an integrative network about the molecular mechanisms controlling fiber length, which are mainly tied to carbohydrate metabolism, cell wall biosynthesis, fatty acid metabolism, secondary metabolism, hormone responses and Transcription factors. The results of this study provide new insights into the critical factors associated with cell elongation and will facilitate further research aimed at understanding the mechanisms underlying cotton fiber elongation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-838) contains supplementary material, which is available to authorized users.

Down-regulating annexin gene GhAnn2 inhibits cotton fiber elongation and decreases Ca2+ influx at the cell apex

Cotton fiber is a single cell that differentiates from the ovule epidermis and undergoes synchronous elongation with high secretion and growth rate. Apart from economic importance, cotton fiber provides an excellent single-celled model for studying mechanisms of cell-growth. Annexins are Ca(2+)- and phospholipid-binding proteins that have been reported to be localized in multiple cellular compartments and involved in control of vesicle secretions. Although several annexins have been found to be highly expressed in elongating cotton fibers, their functional roles in fiber development remain unknown. Here, 14 annexin family members were identified from the fully sequenced diploid G. raimondii (D5 genome), half of which were expressed in fibers of the cultivated tetraploid species G. hirsutum (cv. YZ1). Among them, GhAnn2 from the D genome of the tetraploid species displayed high expression level in elongating fiber. The expression of GhAnn2 could be induced by some phytohormones that play important roles in fiber elongation, such as IAA and GA3. RNAi-mediated down-regulation of GhAnn2 inhibited fiber elongation and secondary cell wall synthesis, resulting in shorter and thinner mature fibers in the transgenic plants. Measurement with non-invasive scanning ion-selective electrode revealed that the rate of Ca(2+) influx from extracellular to intracellular was decreased at the fiber cell apex of GhAnn2 silencing lines, in comparison to that in the wild type. These results indicate that GhAnn2 may regulate fiber development through modulating Ca(2+) fluxes and signaling.

Related to ABA-Insensitive3(ABI3)/Viviparous1 and AtABI5 transcription factor coexpression in cotton enhances drought stress adaptation

Drought tolerance is an important trait being pursued by the agbiotech industry. Abscisic acid (ABA) is a stress hormone that mediates a multitude of processes in growth and development, water use efficiency (WUE) and gene expression during seed development and in response to environmental stresses. Arabidopsis B3-domain transcription factor Related to ABA-Insensitive3 (ABI3)/Viviparous1 (namely AtRAV2) and basic leucine zipper (bZIPs) AtABI5 or AtABF3 transactivated ABA-inducible promoter:GUS reporter expression in a maize mesophyll protoplast transient assay and showed synergies in reporter transactivation when coexpressed. Transgenic cotton (Gossypium hirsutum) expressing AtRAV1/2 and/or AtABI5 showed resistance to imposed drought stress under field and greenhouse conditions and exhibited improved photosynthesis and WUEs associated with absorption through larger root system and greater leaf area. We observed synergy for root biomass accumulation in the greenhouse, intrinsic WUE in the field and drought tolerance in stacked AtRAV and AtABI5 double-transgenic cotton. We assessed AtABI5 and AtRAV1/2 involvement in drought stress adaptations through reactive oxygen species scavenging and osmotic adjustment by marker gene expression in cotton. Deficit irrigation-grown AtRAV1/2 and AtABI5 transgenics had 'less-stressed' molecular and physiological phenotypes under drought, likely due to improved photoassimilation and root and shoot sink strengths and enhanced expression of endogenous GhRAV and genes for antioxidant and osmolyte biosynthesis. Overexpression of bZIP and RAV TFs could impact sustainable cotton agriculture and potentially other crops under limited irrigation conditions.

The calcium sensor GhCaM7 promotes cotton fiber elongation by modulating reactive oxygen species (ROS) production

Fiber elongation is the key determinant of fiber quality and output in cotton (Gossypium hirsutum). Although expression profiling and functional genomics provide some data, the mechanism of fiber development is still not well understood. Here, a gene encoding a calcium sensor, GhCaM7, was isolated based on its high expression level relative to other GhCaMs in fiber cells at the fast elongation stage. The level of expression of GhCaM7 in the wild-type and the fuzzless/lintless mutant correspond to the presence and absence, respectively, of fiber initials. Overexpressing GhCaM7 promotes early fiber elongation, whereas GhCaM7 suppression by RNAi delays fiber initiation and inhibits fiber elongation. Reactive oxygen species (ROS) play important roles in early fiber development. ROS induced by exogenous hydrogen peroxide (H2 O2 ) and Ca(2+) starvation promotes early fiber elongation. GhCaM7 overexpression fiber cells show increased ROS concentrations compared with the wild-type, while GhCaM7 RNAi fiber cells have reduced concentrations. Furthermore, we show that H2 O2 enhances Ca(2+) influx into the fiber and feedback-regulates the expression of GhCaM7. We conclude that GhCaM7, Ca(2+) and ROS are three important regulators involved in early fiber elongation. GhCaM7 might modulate ROS production and act as a molecular link between Ca(2+) and ROS signal pathways in early fiber development.

Comparative evolutionary and developmental dynamics of the cotton (Gossypium hirsutum) fiber transcriptome

The single-celled cotton (Gossypium hirsutum) fiber provides an excellent model to investigate how human selection affects phenotypic evolution. To gain insight into the evolutionary genomics of cotton domestication, we conducted comparative transcriptome profiling of developing cotton fibers using RNA-Seq. Analysis of single-celled fiber transcriptomes from four wild and five domesticated accessions from two developmental time points revealed that at least one-third and likely one-half of the genes in the genome are expressed at any one stage during cotton fiber development. Among these, ∼5,000 genes are differentially expressed during primary and secondary cell wall synthesis between wild and domesticated cottons, with a biased distribution among chromosomes. Transcriptome data implicate a number of biological processes affected by human selection, and suggest that the domestication process has prolonged the duration of fiber elongation in modern cultivated forms. Functional analysis suggested that wild cottons allocate greater resources to stress response pathways, while domestication led to reprogrammed resource allocation toward increased fiber growth, possibly through modulating stress-response networks. This first global transcriptomic analysis using multiple accessions of wild and domesticated cottons is an important step toward a more comprehensive systems perspective on cotton fiber evolution. The understanding that human selection over the past 5,000+ years has dramatically re-wired the cotton fiber transcriptome sets the stage for a deeper understanding of the genetic architecture underlying cotton fiber synthesis and phenotypic evolution.

Ever since Darwin biologists have recognized that comparative study of crop plants and their wild relatives offers a powerful framework for generating insights into the mechanisms that underlie evolutionary change. Here, we study the domestication process in cotton, Gossypium hirsutum, an allopolyploid species (containing two different genomes) which initially was domesticated approximately 5000 years ago, and which primarily is grown for its single-celled seed fibers. Strong directional selection over the millennia was accompanied by transformation of the short, coarse, and brown fibers of wild plants into the long, strong, and fine white fibers of the modern cotton crop plant. To explore the evolutionary genetics of cotton domestication, we conducted transcriptome profiling of developing cotton fibers from multiple accessions of wild and domesticated cottons. Comparative analysis revealed that the domestication process dramatically rewired the transcriptome, affecting more than 5,000 genes, and with a more evenly balanced usage of the duplicated copies arising from genome doubling. We identify many different biological processes that were involved in this transformation, including those leading to a prolongation of fiber elongation and a reallocation of resources toward increased fiber growth in modern forms. The data provide a rich resource for future functional analyses targeting crop improvement and evolutionary objectives.

Effects of fuzzless cottonseed phenotype on cottonseed nutrient composition in near isogenic cotton (Gossypium hirsutum L.) mutant lines under well-watered and water stress conditions

There is no information available on the effect of fuzzless seed trait on cottonseed nutrient composition (minerals, N, S, protein, and oil) under drought stress. The objective of this research was to investigate the effect of the fuzzless seed trait on cottonseed nutrients using five sets of near-isogenic lines (NILs). Each set consists of two lines that share the same genetic background, but differ in seed fuzziness (fuzzy, F; fuzzless, N). The near isogenic lines will enable us to compare the effect of the trait without confounding the genotypic background effects. We hypothesized that since the fuzzless trait involved in fiber initiation development, and was reported to be involved in biochemical, molecular, and genetic processes, this trait may also alter cottonseed nutrient composition. Results showed that NIL sets accumulated different levels of minerals in seeds and leaves, and the fuzzless trait (N) in most of the lines altered seed and leaf mineral accumulations when compared with fuzzy lines (F) or the control line. For example, K, P, Mg, Cu, and Na concentrations in seeds were higher in MD N and STV N than in their equivalent MD F and STV F lines. Leaf concentrations of Ca, K, Mg, S, B, Cu, and Fe in MD N lines were higher than MD F line. Lower levels of nutrients in seeds and leaves were observed under water stress conditions, especially Ca, Mg, N, and B in seeds.Generally and with few exceptions, seed protein was higher in fuzzy lines than in fuzzless lines; however, seed oil was higher in fuzzless lines than in fuzzy lines. Our research demonstrated that fuzzless trait altered the composition and level of nutrients in seed and leaves in well watered and water stressed plants. Differences in protein and oil between fuzzy and fuzzless seeds may indicate alteration in nitrogen and carbon fixation and metabolism. The differential accumulation of seed nutrients in this germplasm could be used by cotton breeders to select for higher cottonseed quality.

Glycoproteome of elongating cotton fiber cells

Cotton ovule epidermal cell differentiation into long fibers primarily depends on wall-oriented processes such as loosening, elongation, remodeling, and maturation. Such processes are governed by cell wall bound structural proteins and interacting carbohydrate active enzymes. Glycosylation plays a major role in the structural, functional, and localization aspects of the cell wall and extracellular destined proteins. Elucidating the glycoproteome of fiber cells would reflect its wall composition as well as compartmental requirement, which must be system specific. Following complementary proteomic approaches, we have identified 334 unique proteins comprising structural and regulatory families. Glycopeptide-based enrichment followed by deglycosylation with PNGase F and A revealed 92 unique peptides containing 106 formerly N-linked glycosylated sites from 67 unique proteins. Our results showed that structural proteins like arabinogalactans and carbohydrate active enzymes were relatively more abundant and showed stage- and isoform-specific expression patterns in the differentiating fiber cell. Furthermore, our data also revealed the presence of heterogeneous and novel forms of structural and regulatory glycoproteins. Comparative analysis with other plant glycoproteomes highlighted the unique composition of the fiber glycoproteome. The present study provides the first insight into the identity, abundance, diversity, and composition of the glycoproteome within single celled cotton fibers. The elucidated composition also indirectly provides clues about unicellular compartmental requirements underlying single cell differentiation.

Gene expression profile analysis of Ligon lintless-1 (Li1) mutant reveals important genes and pathways in cotton leaf and fiber development

Ligon lintless-1 (Li1) is a monogenic dominant mutant of Gossypium hirsutum (upland cotton) with a phenotype of impaired vegetative growth and short lint fibers. Despite years of research involving genetic mapping and gene expression profile analysis of Li1 mutant ovule tissues, the gene remains uncloned and the underlying pathway of cotton fiber elongation is still unclear. In this study, we report the whole genome-level deep-sequencing analysis of leaf tissues of the Li1 mutant. Differentially expressed genes in leaf tissues of mutant versus wild-type (WT) plants are identified, and the underlying pathways and potential genes that control leaf and fiber development are inferred. The results show that transcription factors AS2, YABBY5, and KANDI-like are significantly differentially expressed in mutant tissues compared with WT ones. Interestingly, several fiber development-related genes are found in the downregulated gene list of the mutant leaf transcriptome. These genes include heat shock protein family, cytoskeleton arrangement, cell wall synthesis, energy, H2O2 metabolism-related genes, and WRKY transcription factors. This finding suggests that the genes are involved in leaf morphology determination and fiber elongation. The expression data are also compared with the previously published microarray data of Li1 ovule tissues. Comparative analysis of the ovule transcriptomes of Li1 and WT reveals that a number of pathways important for fiber elongation are enriched in the downregulated gene list at different fiber development stages (0, 6, 9, 12, 15, 18dpa). Differentially expressed genes identified in both leaf and fiber samples are aligned with cotton whole genome sequences and combined with the genetic fine mapping results to identify a list of candidate genes for Li1.

Elucidation of Nuclear and Organellar Genomes of Gossypium hirsutum: Furthering Studies of Species Evolution and Applications for Crop Improvement

Plant genomes are larger and more complex than other eukaryotic organisms, due to small and large duplication events, recombination and subsequent reorganization of the genetic material. Commercially important cotton is the result of a polyploidization event between Old and New World cottons that occurred over one million years ago. Allotetraploid cotton has properties that are dramatically different from its progenitors-most notably, the presence of long, spinnable fibers. Recently, the complete genome of a New World cotton ancestral species, Gossypium raimondii, was completed. Future genome sequencing efforts are focusing on an Old World progenitor, G. arboreum. This sequence information will enable us to gain insights into the evolution of the cotton genome that may be used to understand the evolution of other plant species. The chloroplast genomes of multiple cotton species and races have been determined. This information has also been used to gain insight into the evolutionary history of cotton. Analysis of the database of nuclear and organellar sequences will facilitate the identification of potential genes of interest and subsequent development of strategies for improving cotton.

Transcript profiling by microarray and marker analysis of the short cotton (Gossypium hirsutum L.) fiber mutant Ligon lintless-1 (Li1)

Background

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₁) 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₁ 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.

Results

Physical and morphological measurements of the Li₁ 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 F₂ individuals identified SSR markers associated with the Li₁ 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₁ gene.

Conclusions

The early termination of fiber elongation in the Li₁ 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₁ mutant were identified. Molecular markers closely associated with the Li₁ locus were developed. Results presented here will lay the foundation for further investigation of the genetic and molecular mechanisms of fiber elongation.