A genetic and metabolic analysis revealed that cotton fiber cell development was retarded by flavonoid naringenin

Laccase GhLac1 Modulates Broad-Spectrum Biotic Stress Tolerance via Manipulating Phenylpropanoid Pathway and Jasmonic Acid Synthesis

Plants are constantly challenged by a multitude of pathogens and pests, which causes massive yield and quality losses annually. A promising approach to reduce such losses is to enhance the immune system of plants through genetic engineering. Previous work has shown that laccases (p-diphenol:dioxygen oxidoreductase, EC 1.10.3.2) function as lignin polymerization enzymes. Here we demonstrate that transgenic manipulation of the expression of the laccase gene GhLac1 in cotton (Gossypium hirsutum) can confer an enhanced defense response to both pathogens and pests. Overexpression of GhLac1 leads to increased lignification, associated with increased tolerance to the fungal pathogen Verticillium dahliae and to the insect pests cotton bollworm (Helicoverpa armigera) and cotton aphid (Aphis gosypii). Suppression of GhLac1 expression leads to a redirection of metabolic flux in the phenylpropanoid pathway, causing the accumulation of JA and secondary metabolites that confer resistance to V. dahliae and cotton bollworm; it also leads to increased susceptibility to cotton aphid. Plant laccases therefore provide a new molecular tool to engineer pest and pathogen resistance in crops.

SSR-based association mapping of fiber quality in upland cotton using an eight-way MAGIC population

The quality of fiber is significant in the upland cotton industry. As complex quantitative traits, fiber quality traits are worth studying at a genetic level. To investigate the genetic architecture of fiber quality traits, we conducted an association analysis using a multi-parent advanced generation inter-cross (MAGIC) population developed from eight parents and comprised of 960 lines. The reliable phenotypic data for six major fiber traits of the MAGIC population were collected from five environments in three locations. Phenotypic analysis showed that the MAGIC lines have a wider variation amplitude and coefficient than the founders. A total of 284 polymorphic SSR markers among eight parents screened from a high-density genetic map were used to genotype the MAGIC population. The MAGIC population showed abundant genetic variation and fast linkage disequilibrium (LD) decay (0.76 cM, r² > 0.1), which revealed the advantages of high efficiency and power in QTL exploration. Association mapping via a mixed linear model identified 52 significant loci associated with six fiber quality traits; 14 of them were mapped in reported QTL regions with fiber-related or other agronomic traits. Nine markers demonstrated the pleiotropism that controls more than two fiber traits. Furthermore, two SSR markers, BNL1231 and BNL3452, were authenticated as hotspots that were mapped with multi-traits. In addition, we provided candidate regions and screened six candidate genes for identified loci according to the LD decay distance. Our results provide valuable QTL for further genetic mapping and will facilitate marker-based breeding for fiber quality in cotton.

Simultaneous Editing of Two Copies of Gh14-3-3d Confers Enhanced Transgene-Clean Plant Defense Against Verticillium dahliae in Allotetraploid Upland Cotton

Gossypium hirsutum is an allotetraploid species, meaning that mutants that are difficult to be generated by classical approaches due to gene redundancy. The CRISPR/Cas9 genome editing system is a robust and highly efficient tool for generating target gene mutants, by which the genes of interest may be functionally dissected and applied through genotype-to-phenotype approaches. In this study, the CRISPR/Cas9 genome editing system was developed in G. hirsutum through editing the Gh14-3-3d gene. In T0 transgenic plants, lots of insertions and deletions (indels) in Gh14-3-3d at the expected target site were detected in the allotetraploid cotton At or Dt subgenomes. The results of the PCR, T7EI digestion and sequencing analyses showed that the indels in Gh14-3-3d gene can be stably transmitted to the next generation. Additionally, the indels in the At and Dt subgenomes were segregated in the T1 transgenic plants following Mendelian law, independing on the T-DNA segregation. Two homozygous Gh14-3-3d-edited plants free of T-DNA were chosen by PCR and sequencing assays in the T1 plants, which were called transgene-clean editing plants and were designated ce1 and ce2 in the T2 lines showed higher resistance to Verticillium dahliae infestation compared to the wild-type plants. Thus, the two transgene-clean edited lines can be used as a germplasm to breed disease-resistant cotton cultivars, possibly avoiding complex and expensive safety assessments of the transgenic plants.

A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation

Alternative splicing (AS) is a crucial regulatory mechanism in eukaryotes, which acts by greatly increasing transcriptome diversity. The extent and complexity of AS has been revealed in model plants using high-throughput next-generation sequencing. However, this technique is less effective in accurately identifying transcript isoforms in polyploid species because of the high sequence similarity between coexisting subgenomes. Here we characterize AS in the polyploid species cotton. Using Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq), we developed an integrated pipeline for Iso-Seq transcriptome data analysis (https://github.com/Nextomics/pipeline-for-isoseq). We identified 176 849 full-length transcript isoforms from 44 968 gene models and updated gene annotation. These data led us to identify 15 102 fibre-specific AS events and estimate that c. 51.4% of homoeologous genes produce divergent isoforms in each subgenome. We reveal that AS allows differential regulation of the same gene by miRNAs at the isoform level. We also show that nucleosome occupancy and DNA methylation play a role in defining exons at the chromatin level. This study provides new insights into the complexity and regulation of AS, and will enhance our understanding of AS in polyploid species. Our methodology for Iso-Seq data analysis will be a useful reference for the study of AS in other species.

Interaction between calcium and potassium modulates elongation rate in cotton fiber cells

Calcium (Ca2+) is necessary for fiber cell development in cotton (Gossypium hirsutum), both as a cell wall structural component and for environmental signaling responses. It is also known that potassium (K+) plays a critical role in cotton fiber cell elongation. However, it is unclear whether Ca2+ integrates its activities with K+ to regulate fiber elongation. Here, we report the novel discovery that Ca2+ deficiency, when integrated with K+ signaling, promotes fiber elongation. Using inductively coupled plasma-mass spectrometry (ICP-MS), we determined dynamic profiles of the ionome in ovules and fibers at different developmental stages, and found that a high accumulation of macro-elements, but not Ca2+, was associated with longer fibers. Using an in vitro ovule culture system, we found that under Ca2+-deficient conditions, sufficient K+ (52 mM) rapidly induced ovule and fiber browning, while reduced K+ (2 or 27 mM) not only suppressed tissue browning but also altered fiber elongation. Reduced K+ also enhanced reactive oxygen species scavenging ability and maintained abscisic acid and jasmonic acid levels, which in turn compensated for Ca2+ deficiency. Ca2+ deficiency combined with reduced K+ (0 mM Ca2+ and 27 mM K+) produced longer fibers in cultured ovules, due to cell wall loosening by phytosulfokine (PSK), expansin (EXP), and xyloglucan endotransglycosylase/hydrolase (XTH), and an increase of the K+ content of fiber cells. Using transgenic cotton, we showed that the CBL-INTERACTING PROTEIN KINASE 6 (GhCIPK6) gene mediates the uptake of K+ under Ca2+-deficient conditions. This study establishes a new link between Ca2+, K+, and fiber elongation.

Flavonoid Profile of the Cotton Plant, Gossypium hirsutum: A Review

Cotton, Gossypium hirsutum L., is a plant fibre of significant economic importance, with seeds providing an additional source of protein in human and animal nutrition. Flavonoids play a vital role in maintaining plant health and function and much research has investigated the role of flavonoids in plant defence and plant vigour and the influence these have on cotton production. As part of ongoing research into host plant/invertebrate pest interactions, we investigated the flavonoid profile of cotton reported in published, peer-reviewed literature. Here we report 52 flavonoids representing seven classes and their reported distribution within the cotton plant. We briefly discuss the historical research of flavonoids in cotton production and propose research areas that warrant further investigation.

Comparative transcriptome analysis of cotton fiber development of Upland cotton (Gossypium hirsutum) and Chromosome Segment Substitution Lines from G. hirsutum × G. barbadense

Background

How to develop new cotton varieties possessing high yield traits of Upland cotton and superior fiber quality traits of Sea Island cotton remains a key task for cotton breeders and researchers. While multiple attempts bring in little significant progresses, the development of Chromosome Segment Substitution Lines (CSSLs) from Gossypium barbadense in G. hirsutum background provided ideal materials for aforementioned breeding purposes in upland cotton improvement. Based on the excellent fiber performance and relatively clear chromosome substitution segments information identified by Simple Sequence Repeat (SSR) markers, two CSSLs, MBI9915 and MBI9749, together with the recurrent parent CCRI36 were chosen to conduct transcriptome sequencing during the development stages of fiber elongation and Secondary Cell Wall (SCW) synthesis (from 10DPA and 28DPA), aiming at revealing the mechanism of fiber development and the potential contribution of chromosome substitution segments from Sea Island cotton to fiber development of Upland cotton.

Results

In total, 15 RNA-seq libraries were constructed and sequenced separately, generating 705.433 million clean reads with mean GC content of 45.13% and average Q30 of 90.26%. Through multiple comparisons between libraries, 1801 differentially expressed genes (DEGs) were identified, of which the 902 up-regulated DEGs were mainly involved in cell wall organization and response to oxidative stress and auxin, while the 898 down-regulated ones participated in translation, regulation of transcription, DNA-templated and cytoplasmic translation based on GO annotation and KEGG enrichment analysis. Subsequently, STEM software was performed to explicate the temporal expression pattern of DEGs. Two peroxidases and four flavonoid pathway-related genes were identified in the “oxidation-reduction process”, which could play a role in fiber development and quality formation. Finally, the reliability of RNA-seq data was validated by quantitative real-time PCR of randomly selected 20 genes.

Conclusions

The present report focuses on the similarities and differences of transcriptome profiles between the two CSSLs and the recurrent parent CCRI36 and provides novel insights into the molecular mechanism of fiber development, and into further exploration of the feasible contribution of G. barbadense substitution segments to fiber quality formation, which will lay solid foundation for simultaneously improving fiber yield and quality of upland cotton through CSSLs.

Electronic supplementary material

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

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.

Proteomic Analysis of Differences in Fiber Development between Wild and Cultivated Gossypium hirsutum L

Upland cotton (Gossypium hirsutum L.) is one of the world's most important fiber crops, accounting for more than 90% of all cotton production. While their wild progenitors have relatively short and coarse, often tan-colored fibers, modern cotton cultivars possess longer, finer, stronger, and whiter fiber. In this study, the wild and cultivated cottons (YU-3 and TM-1) selected show significant differences on fibers at 10 days postanthesis (DPA), 20 DPA, and mature stages at the morphological level. To explore the effects of domestication, reveal molecular mechanisms underlying these phenotypic differences, and better inform our efforts to further enhance cotton fiber quality, isobaric tags for relative and absolute protein quantification-facilitated proteomic methods were performed on developing fibers. There were 6990 proteins identified; among them, 336 were defined as differentially expressed proteins between fibers of wild versus domesticated cotton. The down- or up-regulated proteins in wild cotton were involved in phenylpropanoid biosynthesis, zeatin biosynthesis, fatty acid elongation, and other processes. Association analysis between transcriptome and proteome showed positive correlations between transcripts and proteins at both 10 DPA and 20 DPA. Differences in proteomics have been verified at the mRNA level by quantitative real-time polymerase chain reaction and have been validated at the physiological and biochemical levels by POD (peroxidase) activity assays and ZA (zeatin) content estimates. This work corroborates the major pathways involved in cotton fiber development and demonstrates that POD activity and zeatin content have a great potential related to fiber elongation and thickening.

Tissue and cell-specific transcriptomes in cotton reveal the subtleties of gene regulation underlying the diversity of plant secondary cell walls

BACKGROUND

Knowledge of plant secondary cell wall (SCW) regulation and deposition is mainly based on the Arabidopsis model of a 'typical' lignocellulosic SCW. However, SCWs in other plants can vary from this. The SCW of mature cotton seed fibres is highly cellulosic and lacks lignification whereas xylem SCWs are lignocellulosic. We used cotton as a model to study different SCWs and the expression of the genes involved in their formation via RNA deep sequencing and chemical analysis of stem and seed fibre.

RESULTS

Transcriptome comparisons from cotton xylem and pith as well as from a developmental series of seed fibres revealed tissue-specific and developmentally regulated expression of several NAC transcription factors some of which are likely to be important as top tier regulators of SCW formation in xylem and/or seed fibre. A so far undescribed hierarchy was identified between the top tier NAC transcription factors SND1-like and NST1/2 in cotton. Key SCW MYB transcription factors, homologs of Arabidopsis MYB46/83, were practically absent in cotton stem xylem. Lack of expression of other lignin-specific MYBs in seed fibre relative to xylem could account for the lack of lignin deposition in seed fibre. Expression of a MYB103 homolog correlated with temporal expression of SCW CesAs and cellulose synthesis in seed fibres. FLAs were highly expressed and may be important structural components of seed fibre SCWs. Finally, we made the unexpected observation that cell walls in the pith of cotton stems contained lignin and had a higher S:G ratio than in xylem, despite that tissue's lacking many of the gene transcripts normally associated with lignin biosynthesis.

CONCLUSIONS

Our study in cotton confirmed some features of the currently accepted gene regulatory cascade for 'typical' plant SCWs, but also revealed substantial differences, especially with key downstream NACs and MYBs. The lignocellulosic SCW of cotton xylem appears to be achieved differently from that in Arabidopsis. Pith cell walls in cotton stems are compositionally very different from that reported for other plant species, including Arabidopsis. The current definition of a 'typical' primary or secondary cell wall might not be applicable to all cell types in all plant species.

Asymmetric subgenome selection and cis-regulatory divergence during cotton domestication

Comparative population genomics offers an excellent opportunity for unraveling the genetic history of crop domestication. Upland cotton (Gossypium hirsutum) has long been an important economic crop, but a genome-wide and evolutionary understanding of the effects of human selection is lacking. Here, we describe a variation map for 352 wild and domesticated cotton accessions. We scanned 93 domestication sweeps occupying 74 Mb of the A subgenome and 104 Mb of the D subgenome, and identified 19 candidate loci for fiber-quality-related traits through a genome-wide association study. We provide evidence showing asymmetric subgenome domestication for directional selection of long fibers. Global analyses of DNase I-hypersensitive sites and 3D genome architecture, linking functional variants to gene transcription, demonstrate the effects of domestication on cis-regulatory divergence. This study provides new insights into the evolution of gene organization, regulation and adaptation in a major crop, and should serve as a rich resource for genome-based cotton improvement.

A high-efficiency CRISPR/Cas9 system for targeted mutagenesis in Cotton (Gossypium hirsutum L.)

The complex allotetraploid genome is one of major challenges in cotton for repressing gene expression. Developing site-specific DNA mutation is the long-term dream for cotton breeding scientists. The clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system is emerging as a robust biotechnology for targeted-DNA mutation. In this study, two sgRNAs, GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2, were designed in the identical genomic regions of GhMYB25-like A and GhMYB25-like D, which were encoded by cotton A subgenome and the D subgenome, respectively, was assembled to direct Cas9-mediated allotetraploid cotton genome editing. High proportion (14.2–21.4%) CRISPR/Cas9-induced specific truncation events, either from GhMYB25-like A DNA site or from GhMYB25-like D DNA site, were detected in 50% examined transgenic cotton through PCR amplification assay and sequencing analyses. Sequencing results also demonstrated that 100% and 98.8% mutation frequency were occurred on GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2 target site respectively. The off-target effect was evaluated by sequencing two putative off-target sites, which have 3 and 1 mismatched nucleotides with GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2, respectively; all the examined samples were not detected any off-target-caused mutation events. Thus, these results demonstrated that CRISPR/Cas9 is qualified for generating DNA level mutations on allotetraploid cotton genome with high-efficiency and high-specificity.

GhJAZ2 negatively regulates cotton fiber initiation by interacting with the R2R3-MYB transcription factor GhMYB25-like

Jasmonic acid (JA) signaling has been well studied in Arabidopsis. Most reports focus on the role of JA in biological pathways, such as stress resistance, trichome initiation and anthocyanin accumulation. The JASMONATE ZIM-DOMAIN (JAZ) protein is one of the important repressors in the JA signaling pathway. Previous studies showed that JA functions in fiber initiation and elongation, but little is known about how JAZ genes function in fiber development. In this study, a cotton JAZ protein (GhJAZ2) containing a highly conserved TIFY motif and a C-terminal Jas domain was identified, and its function during cotton fiber development was analysed. Gene expression analysis showed that GhJAZ2 was preferentially expressed in the root, hypocotyl, flower and ovule 1 day before anthesis. Overexpression of GhJAZ2 inhibited both lint and fuzz fiber initiation, and reduced the fiber length. Yeast two-hybrid assays showed that GhJAZ2 interacted with the R2R3-MYB transcription factors GhMYB25-like and GhGL1, the bHLH transcription factor GhMYC2, the WD repeat protein GhWD40 and the unknown protein GhJI1. Among these transcription factors, previous studies showed that downregulation of GhMYB25-like leads to a fiberless phenotype in cotton seeds. Molecular and genetic evidence showed that the GhJAZ2 protein suppresses fiber initiation in the overexpressing lines by interacting with GhMYB25-like and suppressing GhMYB25-like activity. Our results suggested that GhJAZ2 functions as a primary transcription repressor during lint and fuzz fiber initiation by interacting with GhMYB25-like, GhGL1, GhMYC2, GhWD40 and GhJI1 to regulate the JA signaling pathway.

Molecular characterization of flavanone 3-hydroxylase gene and flavonoid accumulation in two chemotyped safflower lines in response to methyl jasmonate stimulation

BACKGROUND

Among secondary metabolites, flavonoids are particularly crucial for plant growth, development, and reproduction, as well as beneficial for maintenance of human health. As a flowering plant, safflower has synthesized a striking variety of flavonoids with various pharmacologic properties. However, far less research has been carried out on the genes involved in the biosynthetic pathways that generate these amazing flavonoids, especially characterized quinochalcones. In this study, we first cloned and investigated the participation of a presumed flavanone 3-hydroxylase gene (F3H) from safflower (CtF3H) in a flavonoid biosynthetic pathway.

RESULTS

Bioinformation analysis showed that CtF3H shared high conserved residues and confidence with F3H from other plants. Subcellular localization uncovered the nuclear and cytosol localization of CtF3H in onion epidermal cells. The functional expressions of CtF3H in Escherichia coli BL21(DE3)pLysS cells in the pMAL-C5x vector led to the production of dihydrokaempferol when naringenin was the substrate. Furthermore, the transcriptome expression of CtF3H showed a diametrically opposed expression pattern in a quinochalcone-type safflower line (with orange-yellow flowers) and a flavonol-type safflower line (with white flowers) under external stimulation by methyl jasmonate (MeJA), which has been identified as an elicitor of flavonoid metabolites. Further metabolite analysis showed the increasing tendency of quinochalcones and flavonols, such as hydroxysafflor yellow A, kaempferol-3-O-β-D-glucoside, kaempferol-3-O-β-rutinoside, rutin, carthamin, and luteolin, in the quinochalcone-type safflower line. Also, the accumulation of kaempferol-3-O-β-rutinoside and kaempferol-3-O-β-D-glucoside in flavonols-typed safflower line showed enhanced accumulation pattern after MeJA treatment. However, other flavonols, such as kaempferol, dihydrokaempferol and quercetin-3-O-β-D-glucoside, in flavonols-typed safflower line presented down accumulation respond to MeJA stimulus.

CONCLUSIONS

Our results showed that the high expression of CtF3H in quinochalcone-type safflower line was associated with the accumulation of both quinochalcones and flavonols, whereas its low expression did not affect the increased accumulation of glycosylated derivatives (kaempferol-3-O-β-rutinoside and rutin) in flavonols-typed safflower line but affect the upstream precursors (D-phenylalanine, dihydrokaempferol, kaempferol), which partly revealed the function of CtF3H in different phenotypes and chemotypes of safflower lines.

Fibre elongation requires normal redox homeostasis modulated by cytosolic ascorbate peroxidase in cotton (Gossypium hirsutum)

High-quality cotton fibre equates to a more comfortable textile. Fibre length is an important index of fibre quality. Hydrogen peroxide (H2O2) acts as a signalling molecule in the regulation of fibre elongation. Results from in vitro ovule culture suggest that the alteration of fibre cell H2O2 levels affects fibre development. Ascorbate peroxidase (APX) is an important reactive oxygen species (ROS) scavenging enzyme, and we found that GhAPX1AT/DT encoded one member of the previously unrealized group of cytosolic APXs (cAPXs) that were preferentially expressed during the fibre elongation stage. Transgenic cottons with up- and down-regulation of GhAPX1AT/DT were generated to control fibre endogenous levels of H2O2 Suppression of all cAPX (IAO) resulted in a 3.5-fold increase in H2O2 level in fibres and oxidative stress, which significantly suppressed fibre elongation. The fibre length of transgenic lines with over-expression or specific down-regulation of GhAPX1AT/DT did not show any obvious change. However, the fibres in the over-expression lines exhibited higher tolerance to oxidative stress. Differentially expressed genes (DEGs) in fibres at 10 days post-anthesis (DPA) of IAO lines identified by RNA-seq were related to redox homeostasis, signalling pathways, stress responses and cell wall synthesis, and the DEGs that were up-regulated in IAO lines were also up-regulated in the 10 DPA and 20 DPA fibres of wild cotton compared with domesticated cotton. These results suggest that optimal H2O2 levels and redox state regulated by cytosolic APX are key mechanisms regulating fibre elongation, and dysregulation of the increase in H2O2 induces oxidative stress and results in shorter fibres by initiating secondary cell wall-related gene expression.

Integrative transcriptome, proteome, phosphoproteome and genetic mapping reveals new aspects in a fiberless mutant of cotton

To investigate the molecular mechanisms of fiber initiation in cotton (Gossypium spp.), an integrated approach combining transcriptome, iTRAQ-based proteome and genetic mapping was taken to compare the ovules of the Xuzhou 142 wild type (WT) with its fuzzless-lintless (fl) mutant at −3 and 0 day post-anthesis. A total of 1,953 mRNAs, 187 proteins, and 131 phosphoproteins were differentially expressed (DE) between WT and fl, and the levels of transcripts and their encoded proteins and phosphoproteins were highly congruent. A functional analysis suggested that the abundance of proteins were mainly involved in amino sugar, nucleotide sugar and fatty acid metabolism, one carbon pool for folate metabolism and flavonoid biosynthesis. qRT-PCR, Western blotting, and enzymatic assays were performed to confirm the regulation of these transcripts and proteins. A molecular mapping located the lintless gene li3 in the fl mutant on chromosome 26 for the first time. A further in-silico physical mapping of DE genes with sequence variations between fl and WT identified one and four candidate genes in the li3 and n2 regions, respectively. Taken together, the transcript abundance, phosphorylation status of proteins at the fiber initiation stage and candidate genes have provided insights into regulatory processes underlying cotton fiber initiation.

Multi-omics maps of cotton fibre reveal epigenetic basis for staged single-cell differentiation

Epigenetic modifications are highlighted for their great importance in regulating plant development, but their function associated with single-cell differentiation remains undetermined. Here, we used the cotton fibre, which is the epidermal hair on the cotton ovule, as a model to investigate the regulatory role of DNA methylation in cell differentiation. The level of CHH (H = A, T, or C) DNA methylation level was found to increase during fibre development, accompanied by a decrease in RNA-directed DNA methylation (RdDM). Examination of nucleosome positioning revealed a gradual transition from euchromatin to heterochromatin for chromatin dynamics in developing fibres, which could shape the DNA methylation landscape. The observed increase in DNA methylation in fibres, compared with other ovule tissue, was demonstrated to be mediated predominantly by an active H3K9me2-dependent pathway rather than the RdDM pathway, which was inactive. Furthermore, integrated multi-omics analyses revealed that dynamic DNA methylation played a role in the regulation of lipid biosynthesis and spatio-temporal modulation of reactive oxygen species during fibre differentiation. Our study illustrates two divergent pathways mediating a continuous increase of DNA methylation and also sheds further light on the epigenetic basis for single-cell differentiation in plants. These data and analyses are made available to the wider research community through a comprehensive web portal.

RNA Interference for Functional Genomics and Improvement of Cotton (Gossypium sp.)

RNA interference (RNAi), is a powerful new technology in the discovery of genetic sequence functions, and has become a valuable tool for functional genomics of cotton (Gossypium sp.). The rapid adoption of RNAi has replaced previous antisense technology. RNAi has aided in the discovery of function and biological roles of many key cotton genes involved in fiber development, fertility and somatic embryogenesis, resistance to important biotic and abiotic stresses, and oil and seed quality improvements as well as the key agronomic traits including yield and maturity. Here, we have comparatively reviewed seminal research efforts in previously used antisense approaches and currently applied breakthrough RNAi studies in cotton, analyzing developed RNAi methodologies, achievements, limitations, and future needs in functional characterizations of cotton genes. We also highlighted needed efforts in the development of RNAi-based cotton cultivars, and their safety and risk assessment, small and large-scale field trials, and commercialization.

AtRAV1 and AtRAV2 overexpression in cotton increases fiber length differentially under drought stress and delays flowering

There is a longstanding problem of an inverse relationship between cotton fiber qualities versus high yields. To better understand drought stress signaling and adaptation in cotton (Gossypium hirsutum) fiber development, we expressed the Arabidopsis transcription factors RELATED_TO_ABA-INSENSITIVE3/VIVIPAROUS1/(RAV1) and AtRAV2, which encode APETALA2-Basic3 domain proteins shown to repress transcription of FLOWERING_LOCUS_T (FT) and to promote stomatal opening cell-autonomously. In three years of field trials, we show that AtRAV1 and AtRAV2-overexpressing cotton had ∼5% significantly longer fibers with only marginal decreases in yields under well-watered or drought stress conditions that resulted in 40-60% yield penalties and 3-7% fiber length penalties in control plants. The longer transgenic fibers from drought-stressed transgenics could be spun into yarn which was measurably stronger and more uniform than that from well-watered control fibers. The transgenic AtRAV1 and AtRAV2 lines flowered later and retained bolls at higher nodes, which correlated with repression of endogenous GhFT-Like (FTL) transcript accumulation. Elevated expression early in development of ovules was observed for GhRAV2L, GhMYB25-Like (MYB25L) involved in fiber initiation, and GhMYB2 and GhMYB25 involved in fiber elongation. Altered expression of RAVs controlling critical nodes in developmental and environmental signaling hierarchies has the potential for phenotypic modification of crops.

Long noncoding RNAs and their proposed functions in fibre development of cotton (Gossypium spp.)

Long noncoding RNAs (lncRNAs) are transcripts of at least 200 bp in length, possess no apparent coding capacity and are involved in various biological regulatory processes. Until now, no systematic identification of lncRNAs has been reported in cotton (Gossypium spp.). Here, we describe the identification of 30 550 long intergenic noncoding RNA (lincRNA) loci (50 566 transcripts) and 4718 long noncoding natural antisense transcript (lncNAT) loci (5826 transcripts). LncRNAs are rich in repetitive sequences and preferentially expressed in a tissue-specific manner. The detection of abundant genome-specific and/or lineage-specific lncRNAs indicated their weak evolutionary conservation. Approximately 76% of homoeologous lncRNAs exhibit biased expression patterns towards the At or Dt subgenomes. Compared with protein-coding genes, lncRNAs showed overall higher methylation levels and their expression was less affected by gene body methylation. Expression validation in different cotton accessions and coexpression network construction helped to identify several functional lncRNA candidates involved in cotton fibre initiation and elongation. Analysis of integrated expression from the subgenomes of lncRNAs generating miR397 and its targets as a result of genome polyploidization indicated their pivotal functions in regulating lignin metabolism in domesticated tetraploid cotton fibres. This study provides the first comprehensive identification of lncRNAs in Gossypium.

Metabolomic and transcriptomic insights into how cotton fiber transitions to secondary wall synthesis, represses lignification, and prolongs elongation

Background

The morphogenesis of single-celled cotton fiber includes extreme elongation and staged cell wall differentiation. Designing strategies for improving cotton fiber for textiles and other uses relies on uncovering the related regulatory mechanisms. In this research we compared the transcriptomes and metabolomes of two Gossypium genotypes, Gossypium barbadense cv Phytogen 800 and G. hirsutum cv Deltapine 90. When grown in parallel, the two types of fiber developed similarly except for prolonged fiber elongation in the G. barbadense cultivar. The data were collected from isolated fibers between 10 to 28 days post anthesis (DPA) representing: primary wall synthesis to support elongation; transitional cell wall remodeling; and secondary wall cellulose synthesis, which was accompanied by continuing elongation only in G. barbadense fiber.

Results

Of 206 identified fiber metabolites, 205 were held in common between the two genotypes. Approximately 38,000 transcripts were expressed in the fiber of each genotype, and these were mapped to the reference set and interpreted by homology to known genes. The developmental changes in the transcriptomes and the metabolomes were compared within and across genotypes with several novel implications. Transitional cell wall remodeling is a distinct stable developmental stage lasting at least four days (18 to 21 DPA). Expression of selected cell wall related transcripts was similar between genotypes, but cellulose synthase gene expression patterns were more complex than expected. Lignification was transcriptionally repressed in both genotypes. Oxidative stress was lower in the fiber of G. barbadense cv Phytogen 800 as compared to G. hirsutum cv Deltapine 90. Correspondingly, the G. barbadense cultivar had enhanced capacity for management of reactive oxygen species during its prolonged elongation period, as indicated by a 138-fold increase in ascorbate concentration at 28 DPA.

Conclusions

The parallel data on deep-sequencing transcriptomics and non-targeted metabolomics for two genotypes of single-celled cotton fiber showed that a discrete developmental stage of transitional cell wall remodeling occurs before secondary wall cellulose synthesis begins. The data showed how lignification can be transcriptionally repressed during secondary cell wall synthesis, and they implicated enhanced capacity to manage reactive oxygen species through the ascorbate-glutathione cycle as a positive contributor to fiber length.

Electronic supplementary material

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

The Basic/Helix-Loop-Helix Protein Family in Gossypium: Reference Genes and Their Evolution during Tetraploidization

Basic/helix-loop-helix (bHLH) proteins comprise one of the largest transcription factor families and play important roles in diverse cellular and molecular processes. Comprehensive analyses of the composition and evolution of the bHLH family in cotton are essential to elucidate their functions and the molecular basis of cotton development. By searching bHLH homologous genes in sequenced diploid cotton genomes (Gossypium raimondii and G. arboreum), a set of cotton bHLH reference genes containing 289 paralogs were identified and named as GobHLH001-289. Based on their phylogenetic relationships, these cotton bHLH proteins were clustered into 27 subfamilies. Compared to those in Arabidopsis and cacao, cotton bHLH proteins generally increased in number, but unevenly in different subfamilies. To further uncover evolutionary changes of bHLH genes during tetraploidization of cotton, all genes of S5a and S5b subfamilies in upland cotton and its diploid progenitors were cloned and compared, and their transcript profiles were determined in upland cotton. A total of 10 genes of S5a and S5b subfamilies (doubled from A- and D-genome progenitors) maintained in tetraploid cottons. The major sequence changes in upland cotton included a 15-bp in-frame deletion in GhbHLH130D and a long terminal repeat retrotransposon inserted in GhbHLH062A, which eliminated GhbHLH062A expression in various tissues. The S5a and S5b bHLH genes of A and D genomes (except GobHLH062) showed similar transcription patterns in various tissues including roots, stems, leaves, petals, ovules, and fibers, while the A- and D-genome genes of GobHLH110 and GobHLH130 displayed clearly different transcript profiles during fiber development. In total, this study represented a genome-wide analysis of cotton bHLH family, and revealed significant changes in sequence and expression of these genes in tetraploid cottons, which paved the way for further functional analyses of bHLH genes in the cotton genus.

CHS silencing suggests a negative cross-talk between wax and flavonoid pathways in tomato fruit cuticle

Tomato fruits (Solanum lycopersicum L.) accumulate flavonoids in their cuticle and epidermal cells during ripening. These flavonoids come from de novo biosynthesis due to a significant increase in chalcone synthase (CHS) activity during ripening. Virus-induced gene silencing (VIGS) of tomato fruits have been used to down-regulate SlCHS expression during ripening and analyze the effects at the epidermal and cuticle level. Besides the expected change in fruit color due to a lack of flavonoids incorporated to the cuticle, several other modifications such as a decrease in the amount of cutin and polysaccharides were observed. These indicate a role for either flavonoids or CHS in the alteration of the expression levels of some genes involved in cuticle biosynthesis. Moreover, a negative interaction between the 2 cuticle components, flavonoids and waxes, suggests a relationship between these 2 metabolic pathways.

Stable transformation and expression of GhEXPA8 fiber expansin gene to improve fiber length and micronaire value in cotton

Cotton fiber is multigenic trait controlled by number of genes. Previous studies suggest that one of these genes may be responsible for switching cotton fiber growth on and off to influence the fiber quality produced from a cotton seed. In the present study, the Gossypium hirsutum GhEXPA8 fiber expansin gene was introduced into local cotton variety NIAB 846 by using an Agrobacterium-mediated gene transformation. The neomycin phosphotransferase (NPTII) gene was used as a selection marker for screening of putative transgenic cotton plants. Integration and expression of the fiber expansin gene in cotton plants was confirmed with molecular techniques including Southern blot analyses, real-time PCR. Cellulose assay was used for measurement of cellulose contents of transgenic cotton fiber. The data collected from 3 years of field performance of the transgenic cotton plants expressing GhEXPA8 showed that significant improvement has been made in fiber lengths and micronaire values as compared to control G. hirsutum variety NIAB 846 cotton plants. Statistical techniques were also used for analysis of fiber and agronomic characteristics. The results of this study support improvement of cotton fiber through genetic modification.

Proteomics profiling of fiber development and domestication in upland cotton (Gossypium hirsutum L.)

Comparative proteomic analyses were performed to detail the evolutionary consequences of strong directional selection for enhanced fiber traits in modern upland cotton (Gossypium hirsutum L.). Using two complementary proteomic approaches, 2-DE and iTRAQ LC-MS/MS, fiber proteomes were examined for four representative stages of fiber development. Approximately 1,000 protein features were characterized using each strategy, collectively resulting in the identification and functional categorization of 1,223 proteins. Unequal contributions of homoeologous proteins were detected for over a third of the fiber proteome, but overall expression was balanced with respect to the genome-of-origin in the allopolyploid G. hirsutum. About 30% of the proteins were differentially expressed during fiber development within wild and domesticated cotton. Notably, domestication was accompanied by a doubling of protein developmental dynamics for the period between 10 and 20 days following pollination. Expression levels of 240 iTRAQ proteins and 293 2-DE spots were altered by domestication, collectively representing multiple cellular and metabolic processes, including metabolism, energy, protein synthesis and destination, defense and stress response. Analyses of homoeolog-specific expression indicate that duplicated gene products in cotton fibers can be differently regulated in response to selection. These results demonstrate the power of proteomics for the analysis of crop domestication and phenotypic evolution.

Cotton WRKY1 mediates the plant defense-to-development transition during infection of cotton by Verticillium dahliae by activating JASMONATE ZIM-DOMAIN1 expression

Plants have evolved an elaborate signaling network to ensure an appropriate level of immune response to meet the differing demands of developmental processes. Previous research has demonstrated that DELLA proteins physically interact with JASMONATE ZIM-DOMAIN1 (JAZ1) and dynamically regulate the interaction of the gibberellin (GA) and jasmonate (JA) signaling pathways. However, whether and how the JAZ1-DELLA regulatory node is regulated at the transcriptional level in plants under normal growth conditions or during pathogen infection is not known. Here, we demonstrate multiple functions of cotton (Gossypium barbadense) GbWRKY1 in the plant defense response and during development. Although GbWRKY1 expression is induced rapidly by methyl jasmonate and infection by Verticillium dahliae, our results show that GbWRKY1 is a negative regulator of the JA-mediated defense response and plant resistance to the pathogens Botrytis cinerea and V. dahliae. Under normal growth conditions, GbWRKY1-overexpressing lines displayed GA-associated phenotypes, including organ elongation and early flowering, coupled with the down-regulation of the putative targets of DELLA. We show that the GA-related phenotypes of GbWRKY1-overexpressing plants depend on the constitutive expression of Gossypium hirsutum GhJAZ1. We also show that GhJAZ1 can be transactivated by GbWRKY1 through TGAC core sequences, and the adjacent sequences of this binding site are essential for binding specificity and affinity to GbWRKY1, as revealed by dual-luciferase reporter assays and electrophoretic mobility shift assays. In summary, our data suggest that GbWRKY1 is a critical regulator mediating the plant defense-to-development transition during V. dahliae infection by activating JAZ1 expression.

Comparison of the transcriptome between two cotton lines of different fiber color and quality

To understand the mechanism of fiber development and pigmentation formation, the mRNAs of two cotton lines were sequenced: line Z128 (light brown fiber) was a selected mutant from line Z263 (dark brown fiber). The primary walls of the fiber cell in both Z263 and Z128 contain pigments; more pigments were laid in the lumen of the fiber cell in Z263 compared with that in Z128. However, Z263 contained less cellulose than Z128. A total of 71,895 unigenes were generated: 13,278 (20.26%) unigenes were defined as differentially expressed genes (DEGs) by comparing the library of Z128 with that of Z263; 5,345 (8.16%) unigenes were up-regulated and 7,933 (12.10%) unigenes were down-regulated. qRT-PCR and comparative transcriptional analysis demonstrated that the pigmentation formation in brown cotton fiber was possibly the consequence of an interaction between oxidized tannins and glycosylated anthocyanins. Furthermore, our results showed the pigmentation related genes not only regulated the fiber color but also influenced the fiber quality at the fiber elongation stage (10 DPA). The highly expressed flavonoid gene in the fiber elongation stage could be related to the fiber quality. DEGs analyses also revealed that transcript levels of some fiber development genes (Ca²⁺/CaM, reactive oxygen, ethylene and sucrose phosphate synthase) varied dramatically between these two cotton lines.

Small RNA and degradome profiling reveals a role for miRNAs and their targets in the developing fibers of Gossypium barbadense

microRNAs (miRNAs) are 20-24 nucleotide non-coding small RNAs that play important roles in plant development. The stages of cotton fiber development include initiation, elongation, secondary wall thickening (SWT) and maturation. We constructed seven fiber RNA libraries representing the initiation, elongation and SWT stages. In total, 47 conserved miRNA families and seven candidate miRNAs were profiled using small RNA sequencing. Northern blotting and real-time polymerase chain reaction (PCR) analyses revealed the dynamic expression of miRNAs during fiber development. In addition, 140 targets of 30 conserved miRNAs and 38 targets of five candidate miRNAs were identified through degradome sequencing. Analysis of correlated expression between miRNAs and their targets demonstrated that specific miRNAs suppressed the expression of transcription factors, SBP and MYB, a leucine-rich receptor-like protein kinase, a pectate lyase, α-tubulin, a UDP-glucuronic acid decarboxylase and cytochrome C oxidase subunit 1 to affect fiber development. Histochemical analyses detected the biological activity of miRNA156/157 in ovule and fiber development. Suppressing miRNA156/157 function resulted in the reduction of mature fiber length, illustrating that miRNA156/157 plays an essential role in fiber elongation.

Transcriptome dynamics during fibre development in contrasting genotypes of Gossypium hirsutum L

Understanding the contribution of genetic background in fibre quality traits is important for the development of future cotton varieties with superior fibre quality. We used Affymetrix microarray (Santa Clara, CA) and Roche 454 GSFLX (Branford, CT) for comparative transcriptome analysis between two superior and three inferior genotypes at six fibre developmental stages. Microarray-based analysis of variance (ANOVA) for 89 microarrays encompassing five contrasting genotypes and six developmental stages suggests that the stages of the fibre development have a more pronounced effect on the differentially expressed genes (DEGs) than the genetic background of genotypes. Superior genotypes showed enriched activity of cell wall enzymes, such as pectin methyl esterase, at early elongation stage, enriched metabolic activities such as lipid, amino acid and ribosomal protein subunits at peak elongation, and prolonged combinatorial regulation of brassinosteroid and auxin at later stages. Our efforts on transcriptome sequencing were focused on changes in gene expression at 25 DPA. Transcriptome sequencing resulted in the generation of 475 658 and 429 408 high-quality reads from superior and inferior genotypes, respectively. A total of 24 609 novel transcripts were identified manually for Gossypium hirsutum with no hits in NCBI 'nr' database. Gene ontology analyses showed that the genes for ribosome biogenesis, protein transport and fatty acid biosynthesis were over-represented in superior genotype, whereas salt stress, abscisic acid stimuli and water deprivation leading to the increased proteolytic activity were more pronounced in inferior genotype.

Transcriptome and biochemical analyses revealed a detailed proanthocyanidin biosynthesis pathway in brown cotton fiber

Brown cotton fiber is the major raw material for colored cotton industry. Previous studies have showed that the brown pigments in cotton fiber belong to proanthocyanidins (PAs). To clarify the details of PA biosynthesis pathway in brown cotton fiber, gene expression profiles in developing brown and white fibers were compared via digital gene expression profiling and qRT-PCR. Compared to white cotton fiber, all steps from phenylalanine to PA monomers (flavan-3-ols) were significantly up-regulated in brown fiber. Liquid chromatography mass spectrometry analyses showed that most of free flavan-3-ols in brown fiber were in 2, 3-trans form (gallocatechin and catechin), and the main units of polymeric PAs were trihydroxylated on B ring. Consistent with monomeric composition, the transcript levels of flavonoid 3′, 5′-hydroxylase and leucoanthocyanidin reductase in cotton fiber were much higher than their competing enzymes acting on the same substrates (dihydroflavonol 4-reductase and anthocyanidin synthase, respectively). Taken together, our data revealed a detailed PA biosynthesis pathway wholly activated in brown cotton fiber, and demonstrated that flavonoid 3′, 5′-hydroxylase and leucoanthocyanidin reductase represented the primary flow of PA biosynthesis in cotton fiber.

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.

Transcriptome analysis reveals crosstalk of responsive genes to multiple abiotic stresses in cotton (Gossypium hirsutum L.)

Abiotic stress is a major environmental factor that limits cotton growth and yield, moreover, this problem has become more and more serious recently, as multiple stresses often occur simultaneously due to the global climate change and environmental pollution. In this study, we sought to identify genes involved in diverse stresses including abscisic acid (ABA), cold, drought, salinity and alkalinity by comparative microarray analysis. Our result showed that 5790, 3067, 5608, 778 and 6148 transcripts, were differentially expressed in cotton seedlings under treatment of ABA (1 μM ABA), cold (4°C), drought (200 mM mannitol), salinity (200 mM NaCl) and alkalinity (pH=11) respectively. Among the induced or suppressed genes, 126 transcripts were shared by all of the five kinds of abiotic stresses, with 64 up-regulated and 62 down-regulated. These common members are grouped as stress signal transduction, transcription factors (TFs), stress response/defense proteins, metabolism, transport facilitation, as well as cell wall/structure, according to the function annotation. We also noticed that large proportion of significant differentially expressed genes specifically regulated in response to different stress. Nine of the common transcripts of multiple stresses were selected for further validation with quantitative real time RT-PCR (qRT-PCR). Furthermore, several well characterized TF families, for example, WRKY, MYB, NAC, AP2/ERF and zinc finger were shown to be involved in different stresses. As an original report using comparative microarray to analyze transcriptome of cotton under five abiotic stresses, valuable information about functional genes and related pathways of anti-stress, and/or stress tolerance in cotton seedlings was unveiled in our result. Besides this, some important common factors were focused for detailed identification and characterization. According to our analysis, it suggested that there was crosstalk of responsive genes or pathways to multiple abiotic or even biotic stresses, in cotton. These candidate genes will be worthy of functional study under diverse stresses.

Proteomic profiling of developing cotton fibers from wild and domesticated Gossypium barbadense

Pima cotton (Gossypium barbadense) is widely cultivated because of its long, strong seed trichomes ('fibers') used for premium textiles. These agronomically advanced fibers were derived following domestication and thousands of years of human-mediated crop improvement. To gain an insight into fiber development and evolution, we conducted comparative proteomic and transcriptomic profiling of developing fiber from an elite cultivar and a wild accession. Analyses using isobaric tag for relative and absolute quantification (iTRAQ) LC-MS/MS technology identified 1317 proteins in fiber. Of these, 205 were differentially expressed across developmental stages, and 190 showed differential expression between wild and cultivated forms, 14.4% of the proteome sampled. Human selection may have shifted the timing of developmental modules, such that some occur earlier in domesticated than in wild cotton. A novel approach was used to detect possible biased expression of homoeologous copies of proteins. Results indicate a significant partitioning of duplicate gene expression at the protein level, but an approximately equal degree of bias for each of the two constituent genomes of allopolyploid cotton. Our results demonstrate the power of complementary transcriptomic and proteomic approaches for the study of the domestication process. They also provide a rich database for mining for functional analyses of cotton improvement or evolution.

The flavonoid pathway regulates the petal colors of cotton flower

Although biochemists and geneticists have studied the cotton flower for more than one century, little is known about the molecular mechanisms underlying the dramatic color change that occurs during its short developmental life following blooming. Through the analysis of world cotton germplasms, we found that all of the flowers underwent color changes post-anthesis, but there is a diverse array of petal colors among cotton species, with cream, yellow and red colors dominating the color scheme. Genetic and biochemical analyses indicated that both the original cream and red colors and the color changes post-anthesis were related to flavonoid content. The anthocyanin content and the expression of biosynthesis genes were both increased from blooming to one day post-anthesis (DPA) when the flower was withering and undergoing abscission. Our results indicated that the color changes and flavonoid biosynthesis of cotton flowers were precisely controlled and genetically regulated. In addition, flavonol synthase (FLS) genes involved in flavonol biosynthesis showed specific expression at 11 am when the flowers were fully opened. The anthocyanidin reductase (ANR) genes, which are responsible for proanthocyanidins biosynthesis, showed the highest expression at 6 pm on 0 DPA, when the flowers were withered. Light showed primary, moderate and little effects on flavonol, anthocyanin and proanthocyanidin biosynthesis, respectively. Flavonol biosynthesis was in response to light exposure, while anthocyanin biosynthesis was involved in flower color changes. Further expression analysis of flavonoid genes in flowers of wild type and a flavanone 3-hydroxylase (F3H) silenced line showed that the development of cotton flower color was controlled by a complex interaction between genes and light. These results present novel information regarding flavonoids metabolism and flower development.