Molecular evolution and phylogenetic analysis of genes related to cotton fibers development from wild and domesticated cotton species in Gossypium

The domestication of both diploid and tetraploid cotton species was carried out for fiber utilization. To understand the origin and domestication of fibers, 18 genes related to fiber development were individually cloned and sequenced from 22 different cotton species. Their structures, phylogenetic relationship and molecular evolution were further studied. In the orthologous and homeologous loci of the 18 genes, the sequence and structure of 72.22% were conserved and 27.78% were diverse. Tree topologies constructed based on the combined sequences showed that all 13 D-genome species were congruent with Fryxell's subsection taxonomy, the A- and D-subgenomes independently evolved in the allopolyploid after polyploid formation, and Gossypium raimondii had the closest relationship with all allotetraploids of D-subgenomes. The molecular evolutionary rates revealed approximately equivalent rates among different D-genome species, and purifying selection acted on all genes in the wild D-genome species. Among orthologs and homeologs, the D-subgenomes had higher evolutionary rates than the A-subgenomes in tetraploid cotton species, and the cultivars had higher evolutionary rates than either the semi-domesticated or wild species. Our study revealed that human domestication altered the molecular evolutionary pattern of genes related to fiber development, and Gossypium hirsutum endured greater selective pressures than Gossypium barbadense during the domestication process.     

Analysis of xyloglucan endotransglycosylase/hydrolase (XTH) genes from allotetraploid (Gossypium hirsutum) cotton and its diploid progenitors expressed during fiber elongation

Multiple cellular pathways have been shown to be involved during fiber initiation and elongation stages in the cultivated allotetraploid cotton (Gossypium hirsutum). The cell wall enzymes xyloglucan endotransglycosylase/hydrolases (XTH) have been reported to be associated with the biosynthesis of the cell wall and the growth of cotton fibers, probably regulating the plasticity of the primary cell wall. Among various cotton fiber cDNAs found to be preferentially expressed in cotton fibers, a xyloglucan endotransglycosylase (XTH) cDNA was significantly up-regulated during the elongation stage of cotton fiber development. In the present study, we isolated and characterized genomic clones encoding cotton XTH from cultivated cotton (Gossypium hirsutum) and its diploid progenitors (Gossypium arboreum and Gossypium raimondii), designated GhXTH1-1, GhXTH1-2, GaXTH1 and GrXTH, respectively. In addition, we isolated and characterized, by in silico methods, the putative promoter of XTH1 from Gossypium hirsutum. Sequence analysis revealed more than 50% homology to XTH's at the protein level. DNA gel blot hybridization indicated that at least two copies of GhXTH1 are present in Gossypium hirsutum whereas the diploid progenitor species Gossypium arboreum and Gossypium raimondii has only a single copy. Quantitative real-time PCR and high-resolution melting experiments indicated that in Gossypium hirsutum cultivars, in cotton fibers during early stages of fiber elongation specifically expressing only the GhXTH1-1 gene and expression levels of GhXTH1-1 in fibers varies among cultivars differing in fiber percentage and fiber length.     

The evolution of spinnable cotton fiber entailed prolonged development and a novel metabolism

A central question in evolutionary biology concerns the developmental processes by which new phenotypes arise. An exceptional example of evolutionary innovation is the single-celled seed trichome in Gossypium (“cotton fiber”). We have used fiber development in Gossypium as a system to understand how morphology can rapidly evolve. Fiber has undergone considerable morphological changes between the short, tightly adherent fibers of G. longicalyx and the derived long, spinnable fibers of its closest relative, G. herbaceum, which facilitated cotton domestication. We conducted comparative gene expression profiling across a developmental time-course of fibers from G. longicalyx and G. herbaceum using microarrays with ~22,000 genes. Expression changes between stages were temporally protracted in G. herbaceum relative to G. longicalyx, reflecting a prolongation of the ancestral developmental program. Gene expression and GO analyses showed that many genes involved with stress responses were upregulated early in G. longicalyx fiber development. Several candidate genes upregulated in G. herbaceum have been implicated in regulating redox levels and cell elongation processes. Three genes previously shown to modulate hydrogen peroxide levels were consistently expressed in domesticated and wild cotton species with long fibers, but expression was not detected by quantitative real time-PCR in wild species with short fibers. Hydrogen peroxide is important for cell elongation, but at high concentrations it becomes toxic, activating stress processes that may lead to early onset of secondary cell wall synthesis and the end of cell elongation. These observations suggest that the evolution of long spinnable fibers in cotton was accompanied by novel expression of genes assisting in the regulation of reactive oxygen species levels. Our data suggest a model for the evolutionary origin of a novel morphology through differential gene regulation causing prolongation of an ancestral developmental program.     

Correspondence of trichome mutations in diploid and tetraploid cottons

Quantitative variation for leaf trichome number is observed within and among Gossypium species, varying from glabrous to densely pubescent phenotypes. Moreover, economically important cotton lint fibers are modified trichomes. Earlier studies have mapped quantitative trait loci (QTLs) affecting leaf pubescence in Gossypium using allotetraploids. In this study, we mapped genes responsible for leaf trichome density in a diploid A genome cross. We were able to map 3 QTLs affecting leaf pubescence based on trichome counts obtained from young leaves (YL) and mature leaves (ML). When the F(2) progeny were classified as pubescent versus glabrous, their ratio did not deviate significantly from a 3:1 model, suggesting that glabrousness is inherited in a simple Mendelian fashion. The glabrous mutation mapped to linkage group A3 at the position of major QTL YL1 and ML1 and appeared orthologous to the t1 locus of the allotetraploids. Interestingly, a fiber mutation, sma-4(ha), observed in the same F(2) population cosegregated with the glabrous marker, which indicates either close linkage or common genetic control of lint fiber and leaf trichomes. Studies of A genome diploids may help to clarify the genetic control of trichomes and fiber in both diploid and tetraploid cottons.     

Genotypic and developmental evidence for the role of plasmodesmatal regulation in cotton fiber elongation mediated by callose turnover

Cotton fibers are single-celled hairs that elongate to several centimeters long from the seed coat epidermis of the tetraploid species (Gossypium hirsutum and Gossypium barbadense). Thus, cotton fiber is a unique system to study the mechanisms of rapid cell expansion. Previous work has shown a transient closure of plasmodesmata during fiber elongation (Y.-L. Ruan, D.J. Llewellyn, R.T. Furbank [2001] Plant Cell 13: 47-60). To examine the importance of this closure in fiber elongation, we compared the duration of the plasmodesmata closure among different cotton genotypes differing in fiber length. Confocal imaging of the membrane-impermeant fluorescent molecule carboxyfluorescein revealed a genotypic difference in the duration of the plasmodesmata closure that positively correlates with fiber length among three tetraploid genotypes and two diploid progenitors. In all cases, the closure occurred at the rapid phase of elongation. Aniline blue staining and immunolocalization studies showed that callose deposition and degradation at the fiber base correlates with the timing of plasmodesmata closure and reopening, respectively. Northern analyses showed that the expression of a fiber-specific beta-1,3-glucanase gene, GhGluc1, was undetectable when callose was deposited at the fiber base but became evident at the time of callose degradation. Genotypically, the level of GhGluc1 expression was high in the short fiber genotype and weak in the intermediate and long fiber genotypes. The data provide genotypic and developmental evidence that (1) plasmodesmata closure appears to play an important role in elongating cotton fibers, (2) callose deposition and degradation may be involved in the plasmodesmata closure and reopening, respectively, and (3) the expression of GhGluc1 could play a role in this process by degrading callose, thus opening the plasmodesmata.     

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

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

Sequence Divergence of Microsatellites and Phylogeny Analysis in Tetraploid Cotton Species and Their Putative Diploid Ancestors

To determine the level of microsatellite sequence differences and to use the information to construct a phylogenetic relationship for cultivated tetraploid cotton （Gossypium spp.） species and their putative diploid ancestors, 10 genome-derived microsatellite primer pairs were used to amplify eight species, including two tetraploid and six diploid species, in Gossypium. A total of 92 unique amplicons were resolved using polyacrylamide gel electrophoresis. Each amplicon was cloned, sequenced, and analyzed using standard phylogenetic software. Allelic diversities were caused mostly by changes in the number of simple sequence repeat （SSR） motif repeats and only a small proportion resulted from interruption of the SSR motif within the locus for the same genome. The frequency of base substitutions was 0.5%-1.0% in different genomes, with only few indels found. Based on the combined 10 SSR flanking sequence data, the homology of A-genome diploid species averaged 98.9%, even though most of the amplicons were of the same size, and the sequence homology between G. gossypioides （Ulbr.） Standl. and three other D-genome species （G. raimondii Ulbr., G. davidsonii Kell., and G. thurberi Tod.） was 98.5%, 98.6%, and 98.5%, respectively. Phylogenetic trees of the two allotetraploid species and their putative diploid progenitors showed that homoelogous sequences from the A- and D-subgenome were still present in the polyploid subgenomes and they evolved independently. Meanwhile, homoelogous sequence interaction that duplicated loci in the polyploid subgenomes became phylogenetic sisters was also found in the evolutionary history of tetraploid cotton species. The results of the present study suggest that evaluation of SSR variation at the sequence level can be effective in exploring the evolutionary relationships among Gossypuim species.     

Genome-Wide Analysis of the Sus Gene Family in Cotton

Sucrose synthase (Sus) is a key enzyme in plant sucrose metabolism. In cotton, Sus (EC 2.4.1.13) is the main enzyme that degrades sucrose imported into cotton fibers from the phloem of the seed coat. This study demonstrated that the genomes of Gossypium arboreum L., G. raimondii Ulbr., and G. hirsutum L., contained 8, 8, and 15 Sus genes, respectively. Their structural organizations, phylogenetic relationships, and expression profiles were characterized. Comparisons of genomic and coding sequences identified multiple introns, the number and positions of which were highly conserved between diploid and allotetraploid cotton species. Most of the phylogenetic clades contained sequences from all three species, suggesting that the Sus genes of tetraploid G. hirsutum derived from those of its diploid ancestors. One Sus group (Sus I) underwent expansion during cotton evolution. Expression analyses indicated that most Sus genes were differentially expressed in various tissues and had development-dependent expression profiles in cotton fiber cells. Members of the same orthologous group had very similar expression patterns in all three species. These results provide new insights into the evolution of the cotton Sus gene family, and insight into its members' physiological functions during fiber growth and development.     

利用SSR标记技术研究棉属A、D染色体组的进化

利用SSR分子标记技术,对棉属A、D染色体二倍体及四倍体代表棉种进行了遗传多样性分析。供试的10个二倍体代表棉种间遗传多态性丰富,分子聚类结果与Fryxell棉属分类结果相同。分子水平上进一步揭示出属于D染色体组的拟似棉与其他D染色体组棉种的相似系数最低,A,D染色体组间相似系数很高,该结果支持拟全民族似棉是D染色体组最原始棉种,棉属不同染色体组是共同起源,单元进化的理论,利用栽培的异源四倍体棉种不太适于研究棉属A、D染色体组的进化。     

The cotton kinesin-like calmodulin-binding protein associates with cortical microtubules in cotton fibers

Microtubules in interphase plant cells form a cortical array, which is critical for plant cell morphogenesis. Genetic studies imply that the minus end-directed microtubule motor kinesin-like calmodulin-binding protein (KCBP) plays a role in trichome morphogenesis in Arabidopsis. However, it was not clear whether this motor interacted with interphase microtubules. In cotton (Gossypium hirsutum) fibers, cortical microtubules undergo dramatic reorganization during fiber development. In this study, cDNA clones of the cotton KCBP homolog GhKCBP were isolated from a cotton fiber-specific cDNA library. During cotton fiber development from 10 to 21 DPA, the GhKCBP protein level gradually decreases. By immunofluorescence, GhKCBP was detected as puncta along cortical microtubules in fiber cells of different developmental stages. Thus our results provide evidence that GhKCBP plays a role in interphase cell growth likely by interacting with cortical microtubules. In contrast to fibers, in dividing cells of cotton, GhKCBP localized to the nucleus, the microtubule preprophase band, mitotic spindle, and the phragmoplast. Therefore KCBP likely exerts multiple roles in cell division and cell growth in flowering plants.     

A New Synthetic Allotetraploid (A1A1G2G2) between Gossypium herbaceum and G. australe: Bridging for Simultaneously Transferring Favorable Genes from These Two Diploid Species into Upland Cotton

Gossypium herbaceum, a cultivated diploid cotton species (2n = 2x = 26, A₁A₁), has favorable traits such as excellent drought tolerance and resistance to sucking insects and leaf curl virus. G. australe, a wild diploid cotton species (2n = 2x = 26, G₂G₂), possesses numerous economically valuable characteristics such as delayed pigment gland morphogenesis (which is conducive to the production of seeds with very low levels of gossypol as a potential food source for humans and animals) and resistance to insects, wilt diseases and abiotic stress. Creating synthetic allotetraploid cotton from these two species would lay the foundation for simultaneously transferring favorable genes into cultivated tetraploid cotton. Here, we crossed G. herbaceum (as the maternal parent) with G. australe to produce an F₁ interspecific hybrid and doubled its chromosome complement with colchicine, successfully generating a synthetic tetraploid. The obtained tetraploid was confirmed by morphology, cytology and molecular markers and then self-pollinated. The S₁ seedlings derived from this tetraploid gradually became flavescent after emergence of the fifth true leaf, but they were rescued by grafting and produced S₂ seeds. The rescued S₁ plants were partially fertile due to the existence of univalents at Metaphase I of meiosis, leading to the formation of unbalanced, nonviable gametes lacking complete sets of chromosomes. The S₂ plants grew well and no flavescence was observed, implying that interspecific incompatibility, to some extent, had been alleviated in the S₂ generation. The synthetic allotetraploid will be quite useful for polyploidy evolutionary studies and as a bridge for transferring favorable genes from these two diploid species into Upland cotton through hybridization.