Parallel expression evolution of oxidative stress-related genes in fiber from wild and domesticated diploid and polyploid cotton (Gossypium)

Domestication over Speciation in Allopolyploid Cotton Species: A Stronger Transcriptomic Pull

Cotton has been domesticated independently four times for its fiber, but the genomic targets of selection during each domestication event are mostly unknown. Comparative analysis of the transcriptome during cotton fiber development in wild and cultivated materials holds promise for revealing how independent domestications led to the superficially similar modern cotton fiber phenotype in upland (G. hirsutum) and Pima (G. barbadense) cotton cultivars. Here we examined the fiber transcriptomes of both wild and domesticated G. hirsutum and G. barbadense to compare the effects of speciation versus domestication, performing differential gene expression analysis and coexpression network analysis at four developmental timepoints (5, 10, 15, or 20 days after flowering) spanning primary and secondary wall synthesis. These analyses revealed extensive differential expression between species, timepoints, domestication states, and particularly the intersection of domestication and species. Differential expression was higher when comparing domesticated accessions of the two species than between the wild, indicating that domestication had a greater impact on the transcriptome than speciation. Network analysis showed significant interspecific differences in coexpression network topology, module membership, and connectivity. Despite these differences, some modules or module functions were subject to parallel domestication in both species. Taken together, these results indicate that independent domestication led G. hirsutum and G. barbadense down unique pathways but that it also leveraged similar modules of coexpression to arrive at similar domesticated phenotypes.

Cotton fiber as a model for understanding shifts in cell development under domestication

Cotton fiber provides the predominant plant textile in the world, and it is also a model for plant cell wall biosynthesis. The development of the single-celled cotton fiber takes place across several overlapping but discrete stages, including fiber initiation, elongation, the transition from elongation to secondary cell wall formation, cell wall thickening, and maturation and cell death. During each stage, the developing fiber undergoes a complex restructuring of genome-wide gene expression change and physiological/biosynthetic processes, which ultimately generate a strikingly elongated and nearly pure cellulose product that forms the basis of the global cotton industry. Here, we provide an overview of this developmental process focusing both on its temporal as well as evolutionary dimensions. We suggest potential avenues for further improvement of cotton as a crop plant.

Transcriptional loss of domestication-driven cytoskeletal GhPRF1 gene causes defective floral and fiber development in cotton (Gossypium)

Constitutive- and fiber-specific RNAi of GhPRF1 gene illustrated strong correlation between domestication-driven profilin genes and floral/fiber architecture in cotton. During morpho-transformation of short-fuzz of wild cotton into the elongating spinnable fibers under the millennia of human selection, actin-polymerizing cytoskeletal profilin genes had undergone significant sequence alterations and spatiotemporal shift in their transcription levels. To comprehend the expression dynamics of profilin genes with their phenotypic implications, transgenic expression modulation of cotton profilin 1 (GhPRF1) gene was performed in the constitutive- and fiber-specific manner in Coker 310FR cotton cultivar. The constitutive GhPRF1-RNAi lines (35S:GhPRF1-RNAi) exhibited distorted 'monadelphous' staminal-tube, reduced pollen-viability and poorly developed fibers, whereas floral and fiber development of fiber-specific GhPRF1-RNAi lines showed no abnormalities. Moreover, the fiber-specific GhPRF1 overexpression lines (FBP7:GhPRF1-Ox) showed increased emergence of fiber-initials on the ovule surface, on the contrary to no fiber-initials in fiber-specific RNAi lines (FBP7:GhPRF1-RNAi). Interestingly, the average seed weight and fiber weight of FBP7:GhPRF1-Ox lines increased > 60% and > 38%, respectively, compared with FBP7:GhPRF1-RNAi lines and untransformed control seeds. On a molecular basis, the aberrant floral and fiber development of 35S:GhPRF1-RNAi lines was largely associated with sugar metabolism and hormone-signaling mechanisms. These observations illustrated the strong correlation between domestication-driven GhPRF genes, and floral/fiber development in cotton. Also, the enhanced agronomic traits in GhPRF1-Ox lines of cotton empowered us to recognize their imperative roles, and their future deployment for the sustainable cotton crop improvement.

Global expression dynamics and miRNA evolution profile govern floral/fiber architecture in the modern cotton (Gossypium)

Majority of differentially expressed miRNAs with functional attributes have been recruited independently and parallelly during allopolyploidy followed by the millennia of human selection of both domesticated G. hirsutum and G. barbadense. The genus Gossypium is a marvelous evolutionary model for studying allopolyploidy and morpho-evolution of long-spinnable fibers from the ancestral wild-fuzz. Many genes, transcription factors, and notably, the regulatory miRNAs essentially govern such remarkable modern fiber phenotypes. To comprehend the impact of allopolyploidy on the evolutionary selection of transcriptional dynamicity of key miRNAs, comparative transcriptome profiling of vegetative and fiber tissues of domesticated diploid G. arboreum (A₂) and allopolyploid cotton species G. hirsutum (AD₁), and G. barbadense (AD₂) identified > 300 differentially expressed miRNAs (DEmiRs) within or between corresponding tissues of A₂, AD₁ and AD₂ species. Up to 49% and 32% DEmiRs were up- and down-regulated at fiber initiation stage of AD₁ and AD₂ species, respectively, whereas 50% and 18% DEmiRs were up- and down-regulated at fiber elongation stage of both the allopolyploid species. Interestingly, A-subgenome-specific DEmiRs exhibit expression dominance in the allopolyploid genetic backgrounds. Comparative spatio-temporal expression analyses of AD₁ and AD₂ species discovered that a majority of DEmiRs were recruited independently under millennia of human selection during domestication. Functional annotations of these DEmiRs revealed selection of associated molecular functions such as hormone-signaling, calcium-signaling and reactive oxygen species (ROS) signaling during fiber initiation and elongation. To validate the functional attributes of annotated DEmiRs, we demonstrated for the first time that the target-mimicry-based constitutive diminution of auxin-signaling associated miR167 directly affected the differentiation of floral and fiber tissues of transgenic cotton. These results strongly suggested that the evolutionarily favored DEmiRs including miR167 were involved in the transcriptional regulation of numerous genes during cotton evolution for enhanced fiber-associated agronomic traits.

Evolution of AITR family genes in cotton and their functions in abiotic stress tolerance

Abiotic stresses are major environmental factors inhibiting plant growth and development. AITRs (ABA-induced transcription repressors) are a novel family of transcription factors regulating ABA (abscisic acid) signalling and plant responses to abiotic stresses in Arabidopsis. However, the composition and evolution history of AITRs and their roles in the cotton genus are largely unknown. A total of 12 putative AITRs genes were identified in cultivated tetraploid cotton, Gossypium hirsutum. Phylogenetic analysis of GhAITRs in these tetraploid cottons and their closely related species implicate ancient genome-wide duplication occurring after speciation of Gossypium, and Theobroma could generate duplicates of GhAITRs. Duplicated GhAITRs were stably inherited following diploid speciation and further allotetraploidy in Gossypium. Homologous GhAITRs shared common expression patterns in response to ABA, drought and salinity treatments, and drought tolerance induced in transgenic Arabidopsis plants expressing GhAITR-A1. Together, our findings reveal that duplicates in the GhAITRs gene family were achieved by whole genome duplication rather than three individual duplication events, and that GhAITRs function as transcription repressors and are involved in the regulation of plant responses to ABA and drought stress. These results provide insights towards the improvement of abiotic stress tolerance in cotton using GhAITRs.

QTL controlling fiber quality traits under salt stress in upland cotton (Gossypium hirsutum L.)

QTL for fiber quality traits under salt stress discerned candidate genes controlling fatty acid metabolism. Salinity stress seriously affects plant growth and limits agricultural productivity of crop plants. To dissect the genetic basis of response to salinity stress, a recombinant inbred line population was developed to compare fiber quality in upland cotton (Gossypium hirsutum L.) under salt stress and normal conditions. Based on three datasets of (1) salt stress, (2) normal growth, and (3) the difference value between salt stress and normal conditions, 51, 70, and 53 QTL were mapped, respectively. Three QTL for fiber length (FL) (qFL-Chr1-1, qFL-Chr5-5, and qFL-Chr24-4) were detected under both salt and normal conditions and explained 4.26%, 9.38%, and 3.87% of average phenotypic variation, respectively. Seven genes within intervals of two stable QTL (qFL-Chr1-1 and qFL-Chr5-5) were highly expressed in lines with extreme long fiber. A total of 35 QTL clusters comprised of 107 QTL were located on 18 chromosomes and exhibited pleiotropic effects. Thereinto, two clusters were responsible for improving five fiber quality traits, and 6 influenced FL and fiber strength (FS). The QTL with positive effect for fiber length exhibited active effects on fatty acid synthesis and elongation, but the ones with negative effect played passive roles on fatty acid degradation under salt stress.

Overexpression of plant ferredoxin-like protein promotes salinity tolerance in rice (Oryza sativa)

High-salinity stress is one of the major limiting factors on crop productivity. Physiological strategies against high-salinity stress include generation of reactive oxygen species (ROS), induction of stress-related genes expression, accumulation of abscisic acid (ABA) and up-regulation of antiporters. ROS are metabolism by-products and involved in signal transduction pathway. Constitutive expression of plant ferrodoxin-like protein (PFLP) gene enhances pathogen-resistance activities and root-hair growth through promoting ROS generation. However, the function of PFLP in abiotic stress responses is unclear. In this study, PFLP-1 and PFLP-2-transgenic rice plants were generated to elucidate the role of PFLP under salinity stress. PFLP overexpression significantly increased salt tolerance in PFLP-transgenic rice plants compared with non-transgenic plants (Oryza sativa japonica cv. Tainung 67, designated as TNG67). In high-salinity conditions, PFLP-transgenic plants exhibited earlier ROS production, higher antioxidant enzyme activities, higher ABA accumulation, up-regulated expression of stress-related genes (OsRBOHa, Cu/Zn SOD, OsAPX, OsNCED2, OsSOS1, OsCIPK24, OsCBL4, and OsNHX2), and leaf sodium ion content was lower compared with TNG67 plant. In addition, transgenic lines maintained electron transport rates and contained lower malondialdhyde (MDA) content than TNG67 plant did under salt-stress conditions. Overall results indicated salinity tolerance was improved by PFLP overexpression in transgenic rice plant. The PFLP gene is a potential candidate for improving salinity tolerance for valuable agricultural crops.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton

Extreme elongation distinguishes about one-fourth of cotton (Gossypium sp.) seed epidermal cells as "lint" fibers, useful for the textile industry, from "fuzz" fibers (<5 mm). Ligon lintless-2 (Li ₂ ), a dominant mutation that results in no lint fiber but normal fuzz fiber, offers insight into pathways and mechanisms that differentiate spinnable cotton from its progenitors. A genetic map developed using 1,545 F2 plants showed that marker CISP15 was 0.4 cM from Li ₂ , and "dominant" simple sequence repeat (SSR) markers (i.e. with null alleles in the Li ₂ genotype) SSR7 and SSR18 showed complete linkage with Li ₂ Nonrandom distribution of markers with null alleles suggests that the Li ₂ phenotype results from a 176- to 221-kb deletion of the terminal region of chromosome 18 that may have been masked in prior pooled-sample mapping strategies. The deletion includes 10 genes with putative roles in fiber development. Two Glycosyltransferase Family 1 genes showed striking expression differences during elongation of wild-type versus Li ₂ fiber, and virus-induced silencing of these genes in the wild type induced Li ₂ -like phenotypes. Further, at least 7 of the 10 putative fiber development genes in the deletion region showed higher expression in the wild type than in Li ₂ mutants during fiber development stages, suggesting coordinated regulation of processes in cell wall development and cell elongation, consistent with the hypothesis that some fiber-related quantitative trait loci comprise closely spaced groups of functionally diverse but coordinately regulated genes.

Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F₂ individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.

Eriodictyol can modulate cellular auxin gradients to efficiently promote in vitro cotton fibre development

BACKGROUND

Flavonoids have essential roles in flower pigmentation, fibre development and disease resistance in cotton. Previous studies show that accumulation of naringenin in developing cotton fibres significantly affects fibre growth. This study focused on determining the effects of the flavonoids naringenin, dihydrokaempferol, dihydroquerectin and eriodictyol on fibre development in an in vitro system.

RESULTS

20 μM eriodictyol treatment produced a maximum fibre growth, in terms of fibre length and total fibre units. To gain insight into the associated transcriptional regulatory networks, RNA-seq analysis was performed on eriodictyol-treated elongated fibres, and computational analysis of differentially expressed genes revealed that carbohydrate metabolism and phytohormone signaling pathways were differentially modulated. Eriodictyol treatment also promoted the biosynthesis of quercetin and dihydroquerectin in ovules and elongating fibres through enhanced expression of genes encoding chalcone isomerase, chalcone synthase and flavanone 3-hydroxylase. In addition, auxin biosynthesis and signaling pathway genes were differentially expressed in eriodictyol-driven in vitro fibre elongation. In absence of auxin, eriodictyol predominantly enhanced fibre growth when the localized auxin gradient was disrupted by the auxin transport inhibitor, triiodobenzoic acid.

CONCLUSION

Eriodictyol was found to significantly enhance fibre development through accumulating and maintaining the temporal auxin gradient in developing unicellular cotton fibres.

Genetic Basis of Fiber Improvement and Decreased Stress Tolerance in Cultivated Versus Semi-Domesticated Upland Cotton

Crop domestication from wild ancestors has resulted in the wide adaptation coupled with improved yield and quality traits. However, the genetic basis of many domesticated characteristics remains to be explored. Upland cotton (Gossypium hirsutum) is the most important tetraploid cotton species, accounting for about 90% of world cotton commerce. Here, we reveal the effects of domestication on fiber and stress traits through comprehensive analyses of semi-domesticated races and cultivated cotton accessions. A total of 416 cotton accessions were genotyped, and a decrease in genetic diversity from races to landraces and modern cultivars was detected. Furthermore, 71 domestication selective sweeps (DSS) and 14 improvement selective sweeps (ISS) were identified, with the Dt sub-genome experiencing stronger selection than the At sub-genome during the both selection types. The more expressed genes and a delay in the expression peak of genes related to secondary cell wall (SCW) development in modern cultivars compared to semi-domesticated cotton races, may have contributed to long fibers in these plants. However, down-regulation of genes related to stress response was responsible for decreasing stress tolerance in modern cultivars. We further experimentally confirmed that silencing of PR1 and WRKY20, genes that showed higher expression in the semi-domesticated races, drastically compromised cotton resistance to V. dahliae. Our results reveal fiber improvement and decreased stress tolerance as a result of the domestication of modern upland cotton cultivars.

Genome-wide analysis of UDP-glycosyltransferase super family in Brassica rapa and Brassica oleracea reveals its evolutionary history and functional characterization

Background

Glycosyltransferases comprise a highly divergent and polyphyletic multigene family that is involved in widespread modification of plant secondary metabolites in a process called glycosylation. According to conserved domains identified in their amino acid sequences, these glycosyltransferases can be classified into a single UDP-glycosyltransferase (UGT) 1 superfamily.

Results

We performed genome-wide comparative analysis of UGT genes to trace evolutionary history in algae, bryophytes, pteridophytes, and angiosperms; then, we further investigated the expansion mechanisms and function characterization of UGT gene families in Brassica rapa and Brassica oleracea. Using Hidden Markov Model search, we identified 3, 21, 140, 200, 115, 147, and 147 UGTs in Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Arabidopsis thaliana, B. rapa, and B. oleracea, respectively. Phylogenetic analysis revealed that UGT80 gene family is an ancient gene family, which is shared by all plants and UGT74 gene family is shared by ferns and angiosperms, but the remaining UGT gene families were shared by angiosperms. In dicot lineage, UGTs among three species were classified into three subgroups containing 3, 6, and 12 UGT gene families. Analysis of chromosomal distribution indicates that 98.6 and 71.4% of UGTs were located on B. rapa and B. oleracea pseudo-molecules, respectively. Expansion mechanism analyses uncovered that whole genome duplication event exerted larger influence than tandem duplication on expansion of UGT gene families in B. rapa, and B. oleracea. Analysis of selection forces of UGT orthologous gene pairs in B. rapa, and B. oleracea compared to A. thaliana suggested that orthologous genes in B. rapa, and B. oleracea have undergone negative selection, but there were no significant differences between A. thaliana –B. rapa and A. thaliana –B. oleracea lineages. Our comparisons of expression profiling illustrated that UGTs in B. rapa performed more discrete expression patterns than these in B. oleracea indicating stronger function divergence. Combing with phylogeny and expression analysis, the UGTs in B. rapa and B. oleracea experienced parallel evolution after they diverged from a common ancestor.

Conclusion

We first traced the evolutionary history of UGT gene families in plants and revealed its evolutionary and functional characterization of UGTs in B. rapa, and B. oleracea. This study provides novel insights into the evolutionary history and functional divergence of important traits or phenotype-related gene families in plants.

Electronic supplementary material

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

GhABF2, a bZIP transcription factor, confers drought and salinity tolerance in cotton (Gossypium hirsutum L.)

The bZIP transcription factor (TF) act as an important regulator for the abscisic acid (ABA) mediated abiotic stresses signaling pathways in plants. Here, we reported the cloning and characterization of GhABF2, encoding for typical cotton bZIP TF. Overexpression of GhABF2 significantly improved drought and salt stress tolerance both in Arabidopsis and cotton. However, silencing of GhABF2 made transgenic cotton sensitive to PEG osmotic and salt stress. Expression of GhABF2 was induced by drought and ABA treatments but repressed by high salinity. Transcriptome analysis indicated that GhABF2 increases drought and salt tolerance by regulating genes related to ABA, drought and salt response. The proline contents, activity of superoxide dismutase (SOD) and catalase (CAT) were also significantly increased in GhABF2-overexpression cottons in comparison to wild type after drought and salt treatment. Further, an increase in fiber yield under drought and saline-alkali wetland exhibited the important role of GhABF2 in enhancing the drought and salt tolerance in transgenic lines. In conclusion, manipulation of GhABF2 by biotechnological tools could be a sustainable strategy to deploy drought and salt tolerance in cotton.

Mapping of fiber quality QTLs reveals useful variation and footprints of cotton domestication using introgression lines

Fiber quality improvement is a driving force for further cotton domestication and breeding. Here, QTLs for fiber quality were mapped in 115 introgression lines (ILs) first developed from two intraspecific populations of cultivated and feral cotton landraces. A total of 60 QTLs were found, which explained 2.03–16.85% of the phenotypic variance found in fiber quality traits. A total of 36 markers were associated with five fiber traits, 33 of which were found to be associated with QTLs in multiple environments. In addition, nine pairs of common QTLs were identified; namely, one pair of QTLs for fiber elongation, three pairs for fiber length, three pairs for fiber strength and two pairs for micronaire (qMICs). All common QTLs had additive effects in the same direction in both IL populations. We also found five QTL clusters, allowing cotton breeders to focus their efforts on regions of QTLs with the highest percentages of phenotypic variance. Our results also reveal footprints of domestication; for example, fourteen QTLs with positive effects were found to have remained in modern cultivars during domestication, and two negative qMICs that had never been reported before were found, suggesting that the qMICs regions may be eliminated during artificial selection.

Domestication-driven Gossypium profilin 1 (GhPRF1) gene transduces early flowering phenotype in tobacco by spatial alteration of apical/floral-meristem related gene expression

BACKGROUND

Plant profilin genes encode core cell-wall structural proteins and are evidenced for their up-regulation under cotton domestication. Notwithstanding striking discoveries in the genetics of cell-wall organization in plants, little is explicit about the manner in which profilin-mediated molecular interplay and corresponding networks are altered, especially during cellular signalling of apical meristem determinacy and flower development.

RESULTS

Here we show that the ectopic expression of GhPRF1 gene in tobacco resulted in the hyperactivation of apical meristem and early flowering phenotype with increased flower number in comparison to the control plants. Spatial expression alteration in CLV1, a key meristem-determinacy gene, is induced by the GhPRF1 overexpression in a WUS-dependent manner and mediates cell signalling to promote flowering. But no such expression alterations are recorded in the GhPRF1-RNAi lines. The GhPRF1 transduces key positive flowering regulator AP1 gene via coordinated expression of FT4, SOC1, FLC1 and FT1 genes involved in the apical-to-floral meristem signalling cascade which is consistent with our in silico profilin interaction data. Remarkably, these positive and negative flowering regulators are spatially controlled by the Actin-Related Protein (ARP) genes, specifically ARP4 and ARP6 in proximate association with profilins. This study provides a novel and systematic link between GhPRF1 gene expression and the flower primordium initiation via up-regulation of the ARP genes, and an insight into the functional characterization of GhPRF1 gene acting upstream to the flowering mechanism. Also, the transgenic plants expressing GhPRF1 gene show an increase in the plant height, internode length, leaf size and plant vigor.

CONCLUSIONS

Overexpression of GhPRF1 gene induced early and increased flowering in tobacco with enhanced plant vigor. During apical meristem determinacy and flower development, the GhPRF1 gene directly influences key flowering regulators through ARP-genes, indicating for its role upstream in the apical-to-floral meristem signalling cascade.

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.

Analyzing serial cDNA libraries revealed reactive oxygen species and gibberellins signaling pathways in the salt response of Upland cotton (Gossypium hirsutum L.)

By comparing series full-length cDNA libraries stressed and control, the dynamic process of salt stress response in Upland cotton was studied, and reactive oxygen species and gibberellins signaling pathways were proposed. The Upland cotton is the most important fiber plant with highly salt tolerance. However, the molecular mechanism underlying salt tolerance in domesticated cotton was unclear. Here, seven full-length cDNA libraries were constructed for seedling roots of Upland cotton 'Zhong G 5' at 0, 3, 12 and 48 h after the treatment of control or 150 mM NaCl stress. About 3300 colonies in each library were selected robotically for 5'-end pyrosequencing, resulting in 20,358 expressed sequence tags (ESTs) totally. And 8516 uniESTs were then assembled, including 2914 contigs and 5602 singletons, and explored for Gene Ontology (GO) function. GO comparison between serial stress libraries and control reflected the growth regulation, stimulus response, signal transduction and biology regulation processes were conducted dynamically in response to salt stress. MYB, MYB-related, WRKY, bHLH, GRAS and ERF families of transcription factors were significantly enriched in the early response. 65 differentially expressed genes (DEGs), mainly associated with reactive oxygen species (ROS) scavenging, gibberellins (GAs) metabolism, signal transduction, transcription regulation, stress response and transmembrane transport, were identified and confirmed by quantitative real-time PCR. Overexpression of selected DEGs increased tolerance against salt stress in transgenic yeast. Results in this study supported that a ROS-GAs interacting signaling pathway of salt stress response was activated in Upland cotton. Our results provided valuable gene resources for further investigation of the molecular mechanism of salinity tolerance.

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.

Deciphering the molecular bases for drought tolerance in Arabidopsis autotetraploids

Whole genome duplication (autopolyploidy) is common in many plant species and often leads to better adaptation to adverse environmental conditions. However, little is known about the physiological and molecular mechanisms underlying these adaptations. Drought is one of the major environmental conditions limiting plant growth and development. Here, we report that, in Arabidopsis thaliana, tetraploidy promotes alterations in cell proliferation and organ size in a tissue-dependent manner. Furthermore, it potentiates plant tolerance to salt and drought stresses and decreases transpiration rate, likely through controlling stomata density and closure, abscisic acid (ABA) signalling and reactive oxygen species (ROS) homeostasis. Our transcriptomic analyses revealed that tetraploidy mainly regulates the expression of genes involved in redox homeostasis and ABA and stress response. Taken together, our data have shed light on the molecular basis associated with stress tolerance in autopolyploid plants.

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.

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.

Next generation genetic mapping of the Ligon-lintless-2 (Li₂) locus in upland cotton (Gossypium hirsutum L.)

Mapping-by-sequencing and novel subgenome-specific SNP markers were used to fine map the Ligon-lintless 2 ( Li 2 ) short-fiber gene in tetraploid cotton. These methodologies will accelerate gene identification in polyploid species. Next generation sequencing offers new ways to identify the genetic mechanisms that underlie mutant phenotypes. The release of a reference diploid Gossypium raimondii (D5) genome and bioinformatics tools to sort tetraploid reads into subgenomes has brought cotton genetic mapping into the genomics era. We used multiple high-throughput sequencing approaches to identify the relevant region of reference sequence and identify single nucleotide polymorphisms (SNPs) near the short-fiber mutant Ligon-lintless 2 (Li 2) gene locus. First, we performed RNAseq on 8-day post-anthesis (DPA) fiber cells from the Li 2 mutant and its wild type near isogenic line (NIL) Gossypium hirsutum cv. DP5690. We aligned sequence reads to the D5 genome, sorted the reads into A and D subgenomes with PolyCat and called SNPs with InterSNP. We then identified SNPs that would result in non-synonymous substitutions to amino acid sequences of annotated genes. This step allowed us to identify a 1-Mb region with 24 non-synonymous SNPs, representing the introgressed region that differentiates Li 2 from its NIL. Next, we sequenced total DNA from pools of F2 plants, using a super bulked segregant analysis sequencing (sBSAseq) approach. The sBSAseq predicted 82 non-synonymous SNPs among 3,494 SNPs in a 3-Mb region that includes the region identified by RNAseq. We designed subgenome-specific SNP markers and tested them in an F2 population of 1,733 individuals to construct a genetic map. Our resulting genetic interval contains only one gene, an aquaporin, which is highly expressed in wild-type fibers and is significantly under-expressed in elongating Li 2 fiber cells.

Comparative transcriptome analysis of short fiber mutants Ligon-lintless 1 and 2 reveals common mechanisms pertinent to fiber elongation in cotton (Gossypium hirsutum L.)

Understanding the molecular processes affecting cotton (Gossypium hirsutum) fiber development is important for developing tools aimed at improving fiber quality. Short fiber cotton mutants Ligon-lintless 1 (Li₁) and Ligon-lintless 2 (Li₂) are naturally occurring, monogenic mutations residing on different chromosomes. Both mutations cause early cessation in fiber elongation. These two mutants serve as excellent model systems to elucidate molecular mechanisms relevant to fiber length development. Previous studies of these mutants using transcriptome analysis by our laboratory and others had been limited by the fact that very large numbers of genes showed altered expression patterns in the mutants, making a targeted analysis difficult or impossible. In this research, a comparative microarray analysis was conducted using these two short fiber mutants and their near isogenic wild type (WT) grown under both field and greenhouse environments in order to identify key genes or metabolic pathways common to fiber elongation. Analyses of three transcriptome profiles obtained from different growth conditions and mutant types showed that most differentially expressed genes (DEGs) were affected by growth conditions. Under field conditions, short fiber mutants commanded higher expression of genes related to energy production, manifested by the increasing of mitochondrial electron transport activity or responding to reactive oxygen species when compared to the WT. Eighty-eight DEGs were identified to have altered expression patterns common to both short fiber mutants regardless of growth conditions. Enrichment, pathway and expression analyses suggested that these 88 genes were likely involved in fiber elongation without being affected by growth conditions.

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.

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

The cotton (Gossypium spp.) fiber is a unique elongated cell that is useful for investigating cell differentiation. Previous studies have demonstrated the importance of factors such as sugar metabolism, the cytoskeleton, and hormones, which are commonly known to be involved in plant cell development, while the secondary metabolites have been less regarded. By mining public data and comparing analyses of fiber from two cotton species (Gossypium hirsutum and Gossypium barbadense), we found that the flavonoid metabolism is active in early fiber cell development. Different flavonoids exhibited distinct effects on fiber development during ovule culture; among them, naringenin (NAR) could significantly retard fiber development. NAR is a substrate of flavanone 3-hydroxylase (F3H), and silencing the F3H gene significantly increased the NAR content of fiber cells. Fiber development was suppressed following F3H silencing, but the overexpression of F3H caused no obvious effects. Significant retardation of fiber growth was observed after the introduction of the F3H-RNA interference segment into the high-flavonoid brown fiber G. hirsutum T586 line by cross. A greater accumulation of NAR as well as much shorter fibers were also observed in the BC1 generation plants. These results suggest that NAR is negatively associated with fiber development and that the metabolism mediated by F3H is important in fiber development, thus highlighting that flavonoid metabolism represents a novel pathway with the potential for cotton fiber improvement.

Crop plants as models for understanding plant adaptation and diversification

Since the time of Darwin, biologists have understood the promise of crop plants and their wild relatives for providing insight into the mechanisms of phenotypic evolution. The intense selection imposed by our ancestors during plant domestication and subsequent crop improvement has generated remarkable transformations of plant phenotypes. Unlike evolution in natural settings, descendent and antecedent conditions for crop plants are often both extant, providing opportunities for direct comparisons through crossing and other experimental approaches. Moreover, since domestication has repeatedly generated a suite of "domestication syndrome" traits that are shared among crops, opportunities exist for gaining insight into the genetic and developmental mechanisms that underlie parallel adaptive evolution. Advances in our understanding of the genetic architecture of domestication-related traits have emerged from combining powerful molecular technologies with advanced experimental designs, including nested association mapping, genome-wide association studies, population genetic screens for signatures of selection, and candidate gene approaches. These studies may be combined with high-throughput evaluations of the various "omics" involved in trait transformation, revealing a diversity of underlying causative mutations affecting phenotypes and their downstream propagation through biological networks. We summarize the state of our knowledge of the mutational spectrum that generates phenotypic novelty in domesticated plant species, and our current understanding of how domestication can reshape gene expression networks and emergent phenotypes. An exploration of traits that have been subject to similar selective pressures across crops (e.g., flowering time) suggests that a diversity of targeted genes and causative mutational changes can underlie parallel adaptation in the context of crop evolution.

A bountiful harvest: genomic insights into crop domestication phenotypes

Human selection during crop domestication has resulted in remarkable transformations of plant phenotypes, providing a window into the genetic basis of morphological evolution. Recent progress in our understanding of the genetic architecture of novel plant traits has emerged from combining advanced molecular technologies with improved experimental designs, including nested association mapping, genome-wide association studies, population genetic screens for signatures of selection, and candidate gene approaches. These studies reveal a diversity of underlying causative mutations affecting phenotypes important in plant domestication and crop improvement, including coding sequence substitutions, presence/absence and copy number variation, transposon activation leading to novel gene structures and expression patterns, diversification following gene duplication, and polyploidy leading to altered combinatorial capabilities. The genomic regions unknowingly targeted by human selection include both structural and regulatory genes, often with results that propagate through the transcriptome as well as to other levels in the biosynthetic and morphogenetic networks.

Molecular tools for exploring polyploid genomes in plants

Polyploidy is a very common phenomenon in the plant kingdom, where even diploid species are often described as paleopolyploids. The polyploid condition may bring about several advantages compared to the diploid state. Polyploids often show phenotypes that are not present in their diploid progenitors or exceed the range of the contributing species. Some of these traits may play a role in heterosis or could favor adaptation to new ecological niches. Advances in genomics and sequencing technology may create unprecedented opportunities for discovering and monitoring the molecular effects of polyploidization. Through this review, we provide an overview of technologies and strategies that may allow an in-depth analysis of polyploid genomes. After introducing some basic aspects on the origin and genetics of polyploids, we highlight the main tools available for genome and gene expression analysis and summarize major findings. In the last part of this review, the implications of next generation sequencing are briefly discussed. The accumulation of knowledge on polyploid formation, maintenance, and divergence at whole-genome and subgenome levels will not only help plant biologists to understand how plants have evolved and diversified, but also assist plant breeders in designing new strategies for crop improvement.

Genomes for jeans: cotton genomics for engineering superior fiber

Twenty years ago, scientists predicted that better understanding of fiber development would lead to novel ways to engineer superior cotton fiber. Advances in genetic resources, DNA markers, DNA sequence information, and gene expression data have indeed provided new insights into fiber initiation, elongation and maturation. Many exciting applications of this knowledge offer the potential to select better cotton genotypes more effectively in mainstream breeding programs or engineer genotypes with improved agronomic and/or quality traits. Here, we discuss recent progress in understanding genes involved in fiber development, and their regulation and manipulation to engineer improved fibers. Better understanding of quantitative trait loci/gene interactions that influence fiber quality and yield may help to tailor superior cotton genotypes to diverse environments.

Duplicate gene evolution, homoeologous recombination, and transcriptome characterization in allopolyploid cotton

BACKGROUND

Modern allotetraploid cotton contains an "A" and "D" genome from an ancestral polyploidy event that occurred approximately 1-2 million years ago. Diploid A- and D-genome species can be compared to the A- and D-genomes found within these allotetraploids to make evolutionary inferences about polyploidy. In this paper we present a comprehensive EST assembly derived from diploid and model allotetraploid cottons and demonstrate several evolutionary inferences regarding genic evolution that can be drawn from these data.

RESULTS

We generated a set of cotton expressed sequence tags (ESTs), comprising approximately 4.4 million Sanger and next-generation (454) transcripts supplemented by approximately 152 million Illumina reads from diploid and allotetraploid cottons. From the EST alignments we inferred 259,192 genome-specific single nucleotide polymorphisms (SNPs). Molecular evolutionary analyses of protein-coding regions demonstrate that the rate of nucleotide substitution has increased among both allotetraploid genomes relative to the diploids, and that the ratio of nonsynonymous to synonymous substitutions has increased in one of the two polyploid lineages we sampled. We also use these SNPs to show that a surprisingly high percentage of duplicate genes (~7 %) show a signature of non-independent evolution in the allotetraploid nucleus, having experienced one or more episodes of nonreciprocal homoeologous recombination (NRHR).

CONCLUSIONS

In this study we characterize the functional and mutational properties of the cotton transcriptome, produce a large genome-specific SNP database, and detect illegitimate genetic exchanges between duplicate genomes sharing a common allotetraploid nucleus. Our findings have important implications for our understanding of the consequences of polyploidy and duplicate gene evolution. We demonstrate that cotton genes have experienced an increased rate of molecular evolution following duplication by polyploidy, and that polyploidy has enabled considerable levels of nonreciprocal exchange between homoeologous genes.

Duplicate gene evolution, homoeologous recombination, and transcriptome characterization in allopolyploid cotton

BACKGROUND

Modern allotetraploid cotton contains an "A" and "D" genome from an ancestral polyploidy event that occurred approximately 1-2 million years ago. Diploid A- and D-genome species can be compared to the A- and D-genomes found within these allotetraploids to make evolutionary inferences about polyploidy. In this paper we present a comprehensive EST assembly derived from diploid and model allotetraploid cottons and demonstrate several evolutionary inferences regarding genic evolution that can be drawn from these data.

RESULTS

We generated a set of cotton expressed sequence tags (ESTs), comprising approximately 4.4 million Sanger and next-generation (454) transcripts supplemented by approximately 152 million Illumina reads from diploid and allotetraploid cottons. From the EST alignments we inferred 259,192 genome-specific single nucleotide polymorphisms (SNPs). Molecular evolutionary analyses of protein-coding regions demonstrate that the rate of nucleotide substitution has increased among both allotetraploid genomes relative to the diploids, and that the ratio of nonsynonymous to synonymous substitutions has increased in one of the two polyploid lineages we sampled. We also use these SNPs to show that a surprisingly high percentage of duplicate genes (~7 %) show a signature of non-independent evolution in the allotetraploid nucleus, having experienced one or more episodes of nonreciprocal homoeologous recombination (NRHR).

CONCLUSIONS

In this study we characterize the functional and mutational properties of the cotton transcriptome, produce a large genome-specific SNP database, and detect illegitimate genetic exchanges between duplicate genomes sharing a common allotetraploid nucleus. Our findings have important implications for our understanding of the consequences of polyploidy and duplicate gene evolution. We demonstrate that cotton genes have experienced an increased rate of molecular evolution following duplication by polyploidy, and that polyploidy has enabled considerable levels of nonreciprocal exchange between homoeologous genes.

Parallel up-regulation of the profilin gene family following independent domestication of diploid and allopolyploid cotton (Gossypium)

Cotton is remarkable among our major crops in that four species were independently domesticated, two allopolyploids and two diploids. In each case thousands of years of human selection transformed sparsely flowering, perennial shrubs into highly productive crops with seeds bearing the vastly elongated and abundant single-celled hairs that comprise modern cotton fiber. The genetic underpinnings of these transformations are largely unknown, but comparative gene expression profiling experiments have demonstrated up-regulation of profilin accompanying domestication in all three species for which wild forms are known. Profilins are actin monomer binding proteins that are important in cytoskeletal dynamics and in cotton fiber elongation. We show that Gossypium diploids contain six profilin genes (GPRF1-GPRF6), located on four different chromosomes (eight chromosomes in the allopolyploid). All but one profilin (GPRF6) are expressed during cotton fiber development, and both homeologs of GPRF1-GPRF5 are expressed in fibers of the allopolyploids. Remarkably, quantitative RT-PCR and RNAseq data demonstrate that GPRF1-GPRF5 are all up-regulated, in parallel, in the three independently domesticated cottons in comparison with their wild counterparts. This result was additionally supported by iTRAQ proteomic data. In the allopolyploids, there This usage of novel should be fine, since it refers to a novel evolutionary process, not a novel discovery has been novel recruitment of the sixth profilin gene (GPRF6) as a result of domestication. This parallel up-regulation of an entire gene family in multiple species in response to strong directional selection is without precedent and suggests unwitting selection on one or more upstream transcription factors or other proteins that coordinately exercise control over profilin expression.

A combined functional and structural genomics approach identified an EST-SSR marker with complete linkage to the Ligon lintless-2 genetic locus in cotton (Gossypium hirsutum L.)

BACKGROUND

Cotton fiber length is an important quality attribute to the textile industry and longer fibers can be more efficiently spun into yarns to produce superior fabrics. There is typically a negative correlation between yield and fiber quality traits such as length. An understanding of the regulatory mechanisms controlling fiber length can potentially provide a valuable tool for cotton breeders to improve fiber length while maintaining high yields. The cotton (Gossypium hirsutum L.) fiber mutation Ligon lintless-2 is controlled by a single dominant gene (Li2) that results in significantly shorter fibers than a wild-type. In a near-isogenic state with a wild-type cotton line, Li2 is a model system with which to study fiber elongation.

RESULTS

Two near-isogenic lines of Ligon lintless-2 (Li2) cotton, one mutant and one wild-type, were developed through five generations of backcrosses (BC5). An F2 population was developed from a cross between the two Li2 near-isogenic lines and used to develop a linkage map of the Li2 locus on chromosome 18. Five simple sequence repeat (SSR) markers were closely mapped around the Li2 locus region with two of the markers flanking the Li2 locus at 0.87 and 0.52 centimorgan. No apparent differences in fiber initiation and early fiber elongation were observed between the mutant ovules and the wild-type ones. Gene expression profiling using microarrays suggested roles of reactive oxygen species (ROS) homeostasis and cytokinin regulation in the Li2 mutant phenotype. Microarray gene expression data led to successful identification of an EST-SSR marker (NAU3991) that displayed complete linkage to the Li2 locus.

CONCLUSIONS

In the field of cotton genomics, we report the first successful conversion of gene expression data into an SSR marker that is associated with a genomic region harboring a gene responsible for a fiber trait. The EST-derived SSR marker NAU3991 displayed complete linkage to the Li2 locus on chromosome 18 and resided in a gene with similarity to a putative plectin-related protein. The complete linkage suggests that this expressed sequence may be the Li2 gene.

Generation, annotation and analysis of first large-scale expressed sequence tags from developing fiber of Gossypium barbadense L

Background

Cotton fiber is the world's leading natural fiber used in the manufacture of textiles. Gossypium is also the model plant in the study of polyploidization, evolution, cell elongation, cell wall development, and cellulose biosynthesis. G. barbadense L. is an ideal candidate for providing new genetic variations useful to improve fiber quality for its superior properties. However, little is known about fiber development mechanisms of G. barbadense and only a few molecular resources are available in GenBank.

Methodology and Principal Findings

In total, 10,979 high-quality expressed sequence tags (ESTs) were generated from a normalized fiber cDNA library of G. barbadense. The ESTs were clustered and assembled into 5852 unigenes, consisting of 1492 contigs and 4360 singletons. The blastx result showed 2165 unigenes with significant similarity to known genes and 2687 unigenes with significant similarity to genes of predicted proteins. Functional classification revealed that unigenes were abundant in the functions of binding, catalytic activity, and metabolic pathways of carbohydrate, amino acid, energy, and lipids. The function motif/domain-related cytoskeleton and redox homeostasis were enriched. Among the 5852 unigenes, 282 and 736 unigenes were identified as potential cell wall biosynthesis and transcription factors, respectively. Furthermore, the relationships among cotton species or between cotton and other model plant systems were analyzed. Some putative species-specific unigenes of G. barbadense were highlighted.

Conclusions/Significance

The ESTs generated in this study are from the first large-scale EST project for G. barbadense and significantly enhance the number of G. barbadense ESTs in public databases. This knowledge will contribute to cotton improvements by studying fiber development mechanisms of G. barbadense, establishing a breeding program using marker-assisted selection, and discovering candidate genes related to important agronomic traits of cotton through oligonucleotide array. Our work will also provide important resources for comparative genomics, polyploidization, and genome evolution among Gossypium species.

SHINE transcription factors act redundantly to pattern the archetypal surface of Arabidopsis flower organs

Floral organs display tremendous variation in their exterior that is essential for organogenesis and the interaction with the environment. This diversity in surface characteristics is largely dependent on the composition and structure of their coating cuticular layer. To date, mechanisms of flower organ initiation and identity have been studied extensively, while little is known regarding the regulation of flower organs surface formation, cuticle composition, and its developmental significance. Using a synthetic microRNA approach to simultaneously silence the three SHINE (SHN) clade members, we revealed that these transcription factors act redundantly to shape the surface and morphology of Arabidopsis flowers. It appears that SHNs regulate floral organs' epidermal cell elongation and decoration with nanoridges, particularly in petals. Reduced activity of SHN transcription factors results in floral organs' fusion and earlier abscission that is accompanied by a decrease in cutin load and modified cell wall properties. SHN transcription factors possess target genes within four cutin- and suberin-associated protein families including, CYP86A cytochrome P450s, fatty acyl-CoA reductases, GSDL-motif lipases, and BODYGUARD1-like proteins. The results suggest that alongside controlling cuticular lipids metabolism, SHNs act to modify the epidermis cell wall through altering pectin metabolism and structural proteins. We also provide evidence that surface formation in petals and other floral organs during their growth and elongation or in abscission and dehiscence through SHNs is partially mediated by gibberellin and the DELLA signaling cascade. This study therefore demonstrates the need for a defined composition and structure of the cuticle and cell wall in order to form the archetypal features of floral organs surfaces and control their cell-to-cell separation processes. Furthermore, it will promote future investigation into the relation between the regulation of organ surface patterning and the broader control of flower development and biological functions.

Structure, expression differentiation and evolution of duplicated fiber developmental genes in Gossypium barbadense and G. hirsutum

BACKGROUND

Both Gossypium hirsutum and G. barbadense probably originated from a common ancestor, but they have very different agronomic and fiber quality characters. Here we selected 17 fiber development-related genes to study their structures, tree topologies, chromosomal location and expression patterns to better understand the interspecific divergence of fiber development genes in the two cultivated tetraploid species.

RESULTS

The sequence and structure of 70.59% genes were conserved with the same exon length and numbers in different species, while 29.41% genes showed diversity. There were 15 genes showing independent evolution between the A- and D-subgenomes after polyploid formation, while two evolved via different degrees of colonization. Chromosomal location showed that 22 duplicate genes were located in which at least one fiber quality QTL was detected. The molecular evolutionary rates suggested that the D-subgenome of the allotetraploid underwent rapid evolutionary differentiation, and selection had acted at the tetraploid level. Expression profiles at fiber initiation and early elongation showed that the transcripts levels of most genes were higher in Hai7124 than in TM-1. During the primary-secondary transition period, expression of most genes peaked earlier in TM-1 than in Hai7124. Homeolog expression profile showed that A-subgenome, or the combination of A- and D-subgenomes, played critical roles in fiber quality divergence of G. hirsutum and G. barbadense. However, the expression of D-subgenome alone also played an important role.

CONCLUSION

Integrating analysis of the structure and expression to fiber development genes, suggests selective breeding for certain desirable fiber qualities played an important role in divergence of G. hirsutum and G. barbadense.

Gene expression in developing fibres of Upland cotton (Gossypium hirsutum L.) was massively altered by domestication

Background

Understanding the evolutionary genetics of modern crop phenotypes has a dual relevance to evolutionary biology and crop improvement. Modern upland cotton (Gossypium hirsutum L.) was developed following thousands of years of artificial selection from a wild form, G. hirsutum var. yucatanense, which bears a shorter, sparser, layer of single-celled, ovular trichomes ('fibre'). In order to gain an insight into the nature of the developmental genetic transformations that accompanied domestication and crop improvement, we studied the transcriptomes of cotton fibres from wild and domesticated accessions over a developmental time course.

Results

Fibre cells were harvested between 2 and 25 days post-anthesis and encompassed the primary and secondary wall synthesis stages. Using amplified messenger RNA and a custom microarray platform designed to interrogate expression for 40,430 genes, we determined global patterns of expression during fibre development. The fibre transcriptome of domesticated cotton is far more dynamic than that of wild cotton, with over twice as many genes being differentially expressed during development (12,626 versus 5273). Remarkably, a total of 9465 genes were diagnosed as differentially expressed between wild and domesticated fibres when summed across five key developmental time points. Human selection during the initial domestication and subsequent crop improvement has resulted in a biased upregulation of components of the transcriptional network that are important for agronomically advanced fibre, especially in the early stages of development. About 15% of the differentially expressed genes in wild versus domesticated cotton fibre have no homology to the genes in databases.

Conclusions

We show that artificial selection during crop domestication can radically alter the transcriptional developmental network of even a single-celled structure, affecting nearly a quarter of the genes in the genome. Gene expression during fibre development within accessions and expression alteration arising from evolutionary change appears to be 'modular' - complex genic networks have been simultaneously and similarly transformed, in a coordinated fashion, as a consequence of human-mediated selection. These results highlight the complex alteration of the global gene expression machinery that resulted from human selection for a longer, stronger and finer fibre, as well as other aspects of fibre physiology that were not consciously selected. We illustrate how the data can be mined for genes that were unwittingly targeted by aboriginal and/or modern domesticators during crop improvement and/or which potentially control the improved qualities of domesticated cotton fibre.

See Commentary: http://www.biomedcentral.com/1741-7007/8/137

Transcriptome divergence between the hexaploid salt-marsh sister species Spartina maritima and Spartina alterniflora (Poaceae)

Invasive species are ideal model systems to investigate the evolutionary processes associated with their ecological success by comparison with closely related species. In this article, we explore transcriptome evolution following divergence between two closely related salt-marsh species, the invasive Spartina alterniflora (native to the East-American Atlantic coast, introduced in several continents) and the declining Spartina maritima (native to the Euro-African Atlantic coast). We have explored the utility of cross-species hybridization microarrays using rice (Oryza sativa) oligo-microarrays to compare leaf expression patterns between these species. Coding sequence comparisons from 10 nuclear genes (2256 bp) revealed that nucleotide divergence between Spartina and Oryza range from 8% to 14%. More than 70% of the 60-mer oligonucleotide sequences spotted on the rice microarray exhibited stable and repeatable patterns when hybridized against Spartina RNA. In total, 9353 (44.5%) genes on the array hybridized with both species S. maritima and S. alterniflora. Among these genes, 1247 genes were found to be differentially expressed between the two Spartina species, most of them (957) being up-regulated in S. alterniflora. In particular, developmental and cellular growth genes (gene ontology, biological process) were highly up-regulated in S. alterniflora and down-regulated in S. maritima, whereas genes involved in stress response were up-regulated in S. maritima. Our findings indicate the suitability of cross-species microarray hybridization between Spartina and O. sativa and reveal the extent of leaf transcriptome evolution that took place during the divergence between S. alterniflora and S. maritima. Expression patterns are consistent with the morphological differentiation and differential expansion of the two species.

Evolutionary rate variation, genomic dominance and duplicate gene expression evolution during allotetraploid cotton speciation

Here, we describe the evolution of gene expression among a diversified cohort of five allopolyploid species in the cotton genus (Gossypium). Using this phylogenetic framework and comparisons with expression changes accompanying F(1) hybridization, we provide a temporal perspective on expression diversification following a shared genome duplication. Global patterns of gene expression were studied by the hybridization of petal RNAs to a custom microarray. This platform measures total expression for c. 42 000 duplicated genes, and genome-specific expression for c. 1400 homoeologs (genes duplicated by polyploidy). We report homoeolog expression bias favoring the allopolyploid D genome over the A genome in all species (among five polyploid species, D biases ranging from c. 54 to 60%), in addition to conservation of biases among genes. Furthermore, we find surprising levels of transgressive up- and down-regulation in the allopolyploids, a diminution of the level of bias in genomic expression dominance but not in its magnitude, and high levels of rate variation among allotetraploid species. We illustrate how phylogenetic and temporal components of expression evolution may be partitioned and revealed following allopolyploidy. Overall patterns of expression evolution are similar among the Gossypium allotetraploids, notwithstanding a high level of interspecific rate variation, but differ strikingly from the direction of genomic expression dominance patterns in the synthetic F(1) hybrid.