Deep sequencing reveals differences in the transcriptional landscapes of fibers from two cultivated species of cotton.

Cotton (Gossypium) fiber is the most prevalent natural product used in the textile industry. The two major cultivated species, G. hirsutum (Gh) and G. barbadense (Gb), are allotetraploids with contrasting fiber quality properties. To better understand the molecular basis for their fiber differences, EST pyrosequencing was used to document the fiber transcriptomes at two key development stages, 10 days post anthesis (dpa), representing the peak of fiber elongation, and 22 dpa, representing the transition to secondary cell wall synthesis. The 617,000 high quality reads (89% of the total 692,000 reads) from 4 libraries were assembled into 46,072 unigenes, comprising 38,297 contigs and 7,775 singletons. Functional annotation of the unigenes together with comparative digital gene expression (DGE) revealed a diverse set of functions and processes that were partly linked to specific fiber stages. Globally, 2,770 contigs (7%) showed differential expression (>2-fold) between 10 and 22 dpa (irrespective of genotype), with 70% more highly expressed at 10 dpa, while 2,248 (6%) were differentially expressed between the genotypes (irrespective of stage). The most significant genes with differential DGE at 10 dpa included expansins and lipid transfer proteins (higher in Gb), while at 22 dpa tubulins, cellulose, and sucrose synthases showed higher expression in Gb. DGE was compared with expression data of 10 dpa-old fibers from Affymetrix microarrays. Among 543 contigs showing differential expression on both platforms, 74% were consistent in being either over-expressed in Gh (242 genes) or in Gb (161 genes). Furthermore, the unigene set served to identify 339 new SSRs and close to 21,000 inter-genotypic SNPs. Subsets of 88 SSRs and 48 SNPs were validated through mapping and added 65 new loci to a RIL genetic map. The new set of fiber ESTs and the gene-based markers complement existing available resources useful in basic and applied research for crop improvement in cotton.

Meta-analysis of cotton fiber quality QTLs across diverse environments in a Gossypium hirsutum x G. barbadense RIL population.

BACKGROUND: Cotton fibers (produced by Gossypium species) are the premier natural fibers for textile production. The two tetraploid species, G. barbadense (Gb) and G. hirsutum (Gh), differ significantly in their fiber properties, the former having much longer, finer and stronger fibers that are highly prized. A better understanding of the genetics and underlying biological causes of these differences will aid further improvement of cotton quality through breeding and biotechnology. We evaluated an inter-specific Gh x Gb recombinant inbred line (RIL) population for fiber characteristics in 11 independent experiments under field and glasshouse conditions. Sites were located on 4 continents and 5 countries and some locations were analyzed over multiple years. RESULTS: The RIL population displayed a large variability for all major fiber traits. QTL analyses were performed on a per-site basis by composite interval mapping. Among the 651 putative QTLs (LOD > 2), 167 had a LOD exceeding permutation based thresholds. Coincidence in QTL location across data sets was assessed for the fiber trait categories strength, elongation, length, length uniformity, fineness/maturity, and color. A meta-analysis of more than a thousand putative QTLs was conducted with MetaQTL software to integrate QTL data from the RIL and 3 backcross populations (from the same parents) and to compare them with the literature. Although the global level of congruence across experiments and populations was generally moderate, the QTL clustering was possible for 30 trait x chromosome combinations (5 traits in 19 different chromosomes) where an effective co-localization of unidirectional (similar sign of additivity) QTLs from at least 5 different data sets was observed. Most consistent meta-clusters were identified for fiber color on chromosomes c6, c8 and c25, fineness on c15, and fiber length on c3. CONCLUSIONS: Meta-analysis provided a reliable means of integrating phenotypic and genetic mapping data across multiple populations and environments for complex fiber traits. The consistent chromosomal regions contributing to fiber quality traits constitute good candidates for the further dissection of the genetic and genomic factors underlying important fiber characteristics, and for marker-assisted selection.

Transcript profiling during fiber development identifies pathways in secondary metabolism and cell wall structure that may contribute to cotton fiber quality.

A global gene expression profiling study at different stages of fiber development was undertaken on two cotton species cultivated for fiber, Gossypium hirsutum (L.) and G. barbadense (L.). A large proportion of the genome was expressed during both fiber elongation and subsequent secondary cell wall thickening. There was a major shift in abundance of transcripts for gene regulation, cell organization and metabolism between fiber elongation and fiber thickening that was fundamentally similar in both species. Each stage had its own distinctive features represented by specific metabolic and regulatory genes, a number of which have been noted previously. Many of the genes expressed in the fibers were of a similar type and developmental expression to those seen in other fiber-producing plants, indicating a conservation of mechanisms of cell elongation and wall thickening across diverse plant genera. Secondary metabolism and pectin synthesis and modification genes were amongst the most statistically significant differentially expressed categories between the two species during fiber elongation. The gene profiles of the fiber thickening stage, however, were almost identical between the two species, suggesting that their different final fiber quality properties may be established at earlier stages of fiber development. Expression levels of representative phenylpropanoid and pectin modification genes showed high correlations with specific fiber properties in an inter-specific cotton recombinant inbred line (RIL) population, supporting a role in determining fiber quality.

PRD, an Arabidopsis AINTEGUMENTA-like gene, is involved in root architectural changes in response to phosphate starvation.

Changes in root architecture are one of the adaptive strategies used by plants to compensate for local phosphate (Pi) deficiency in soils. Root architecture variables triggered by Pi availability are well documented in Arabidopsis (Arabidopsis thaliana), but the molecular mechanisms behind these adaptive responses remain to be elucidated. By the use of transcriptomic and quantitative RT-PCR analysis, we observed that an AINTEGUMENTA-like gene, named PRD for Phosphate Root Development, was rapidly repressed in roots under low Pi conditions. The physiological function of the PRD gene was analyzed through the null allele mutant prd, which displayed less development of primary and lateral roots under Pi-starvation conditions than wild-type plants. Complementation of the prd mutant with the wild-type gene led to a similar response to Pi starvation as wild-type plants, indicating the complete rescue of the mutant phenotype. These results suggest that PRD gene is involved in the regulation of root architectural responses to Pi starvation by controlling primary and lateral root elongation. This model is in agreement with the tissue-specific pattern of PRD gene expression, which was observed to occur specifically in the apex in both the primary and lateral roots. However, Pi influx, anionic profiles and root expression of genes typically induced by Pi starvation, such as high affinity Pi transporters (PHT1;1 and PHT1;4) and an acid phosphatase (AtACP5), were similar in wild type and prd plants in response to Pi starvation. These results support the hypothesis that the PRD gene is not a checkpoint for Pi-starvation responses, but acts specifically as a regulator of root architectural responses to Pi starvation.

Laser capture microdissection and cDNA microarrays used to generate gene expression profiles of the rapidly expanding fibre initial cells on the surface of cotton ovules.

Cotton (Gossypium hirsutum L.) fibre initial cells undergo a rapid cellular re-programming around anthesis to form the long cellulose fibres prized for textile manufacture. On the day of anthesis the cotton fibre initial cells balloon out from the ovule surface and so are clearly distinguished from adjacent epidermal pavement cells. To enhance our understanding of the molecular processes that determine which cells become fibres and why adjacent epidermal cells remain in a different developmental state we studied the expression profiles of the two respective cell types. Using laser-capture microdissection, coupled with an in vitro RNA amplification system, we used cDNA microarray slides to profile the gene expression in expanding fibre initials compared to the non-expanding epidermal cells at an early stage just after the fibre initials are discernable. Except for a few regulatory genes, the genes that are up-regulated in the cotton fibre initials relative to epidermal cells predominantly encode proteins involved in generating the components for the extra cell membrane and primary cell wall needed for the rapid cell expansion of the initials. This includes synthesis of enzymes and cell wall proteins, carbohydrates, and lipids. An analysis of single channel fluorescence levels confirmed that these classes of genes were also the most highly expressed genes in fibre initials. Genes involved in DNA metabolism were also well represented in the expanding fibre cell, consistent with the limited endoreduplication we previously reported to occur in fibre initial cells.