Evolution and expression of MYB genes in diploid and polyploid cotton

Genome-wide identification of xylan glucuronosyltransferase family in cotton and function characterization of GhGUX5 in regulating Verticillium wilt resistance

Xylan glucuronosyltransferase (GUX) is widely involved in a variety of physiological processes in plants, including plant development, growth and the defense response to pathogens. However, the function of GUX regulators in Verticillium dahliae (V. dahliae) infection has not been considered previously in cotton. Overall, 119 GUX genes were identified from multiple species and were phylogenetically categorized into seven classes. Duplication event analysis indicated that GUXs in Gossypium hirsutum primarily originated from segmental duplication. GhGUXs promoter analysis indicated cis-regulatory elements capable of reacting to several different stresses. RNA-Seq data and qRT-PCR analysis both indicated that most GhGUXs were associated with V. dahliae infection. Gene interaction network analysis showed that GhGUX5 interacted with 11 proteins, and the relative expression of these 11 proteins changed significantly following V. dahliae infection. In addition, silencing and overexpression of GhGUX5 results to enhance and reduce plant's susceptibility to V. dahliae. Further study showed that TRV: GhGUX5 silenced cotton plants exhibited a decrease in the degree of lignification, total lignin content, gene expression levels involved in lignin biosynthesis, and enzyme activity compared with TRV: 00. The above results indicate that GhGUX5 enhances Verticillium wilt resistance through the lignin biosynthesis pathway.

Genome-wide Characterization of the JmjC Domain-Containing Histone Demethylase Gene Family Reveals GhJMJ24 and GhJMJ49 Involving in Somatic Embryogenesis Process in Cotton

The Jumonji C (JmjC) domain-containing protein family, an important family of histone demethylase in plants, can directly reverse histone methylation and play important roles in various growth and development processes. In the present study, 51 JmjC genes (GhJMJs) were identified by genome-wide analysis in upland cotton (Gossypium hirsutum), which can be categorized into six distinct groups by phylogenetic analysis. Extensive syntenic relationship events were found between G. hirsutum and Theobroma cacao. We have further explored the putative molecular regulatory mechanisms of the JmjC gene family in cotton. GhJMJ24 and GhJMJ49 were both preferentially expressed in embryogenic callus compared to nonembryogenic callus in cotton tissue culture, which might be regulated by transcription factors and microRNAs to some extent. Further experiments indicated that GhJMJ24 and GhJMJ49 might interact with SUVH4, SUVH6, DDM1, CMT3, and CMT1 in the nucleus, potentially in association with demethylation of H3K9me2. Taken together, our results provide a foundation for future research on the biological functions of GhJMJ genes in cotton, especially in somatic embryogenesis in cotton tissue culture, which is crucial for the regeneration of transgenic plants.

TaMYB29: A Novel R2R3-MYB Transcription Factor Involved in Wheat Defense Against Stripe Rust

Members of the R2R3-MYB transcription factor superfamily have been implicated in plant development, improved disease resistance, and defense responses to several types of stresses. To study the function of TaMYB29 transcription factor-a member of the R2R3-MYB superfamily-in response to an avirulent race of stripe rust pathogen, Puccinia striiformis f. sp. tritici (Pst), we identified and cloned the TaMYB29 gene from wheat cultivar (cv.) AvS+Yr10 following infection with Pst. The TaMYB29 protein, comprising 261 amino acids, contains two highly conserved MYB domains. We first showed that TaMYB29 is a transcription factor, whose transcriptional levels are significantly induced by salicylic acid (SA), abscisic acid (ABA), jasmonic acid (JA), ethylene (ET), and Pst. The results showed that TaMYB29 is involved in the wheat response to stipe rust. The overexpression of the TaMYB29 gene resulted in the accumulation of reactive oxygen species (ROS) and pathogen-independent cell death in Nicotiana benthamiana leaves. The silencing of TaMYB29 gene in wheat cv. AvS+Yr10, containing the stripe rust resistance gene Yr10, promoted hyphae growth, significantly downregulated the expression of pathogenesis-related (PR) genes, and substantially reduced the wheat resistance to Pst compared with the non-silenced control. In addition, the accumulation of hydrogen peroxide (H2O2) significantly decreased, and the activity of catalase, an enzyme required for H2O2 scavenging, was elevated. Altogether, TaMYB29 positively regulates the defense response against stripe rust in wheat AvS+Yr10 by enhancing H2O2 accumulation, PR gene expression, and SA signaling pathway-induced cell death. These results provide new insights into the contribution of TaMYB29 to the defense response against rust pathogens in wheat.

Genome-wide identification and evolution of Dof transcription factor family in cultivated and ancestral cotton species

The DNA-binding with one finger (Dof) proteins are transcription factors involved in many biological processes in plants. To predict the evolutionary pattern, a genome-wide in-silico analysis of Dof TFs family in diploid (Gossypium arboreum and Gossypium raimondii) and allotetraploid (Gossypium hirsutum and Gossypium barbadense) cotton species were carried out. In G. arboreum, we have identified 58 non-redundant genes encoding Dof proteins renamed as GaDof (G. arboreum Dof), 55 Dof genes were identified in G. raimondii (GrDof), 89 were predicted ffrom G. hirsutum (GhDof) and the highest, 110 Dof genes were identified in G. barbadense (GbDof). The phylogenetic analysis, physical location, gene structure, conserved domain analyses were also investigated for G. arboreum, G. raimondii, and G. hirsutum. The gene expression pattern in G. hirsutum, at different growth stages, revealing the probable involvement of some GhDof genes in growth and development. These genes may improve seed germination and growth in cotton.

Evolutionary comparison of competitive protein-complex formation of MYB, bHLH, and WDR proteins in plants

A protein complex consisting of a MYB, basic Helix-Loop-Helix, and a WDR protein, the MBW complex, regulates five traits, namely the production of anthocyanidin, proanthocyanidin, and seed-coat mucilage, and the development of trichomes and root hairs. For complexes involved in trichome and root hair development it has been shown that the interaction of two MBW proteins can be counteracted by the respective third protein (called competitive complex formation). We examined competitive complex formation for selected MBW proteins from Arabidopsis thaliana, Arabis alpina, Gossypium hirsutum, Petunia hybrida, and Zea mays. Quantitative analyses of the competitive binding of MYBs and WDRs to bHLHs were done by pull-down assays using ProtA- and luciferase-tagged proteins expressed in human HEC cells. We found that some bHLHs show competitive complex formation whilst others do not. Competitive complex formation strongly correlated with a phylogenetic tree constructed with the bHLH proteins under investigation, suggesting a functional relevance. We demonstrate that this different behavior can be explained by changes in one amino acid and that this position is functionally relevant in trichome development but not in anthocyanidin regulation.

Genetics and evolution of MIXTA genes regulating cotton lint fiber development

Cotton, with cellulose-enriched mature fibers, is the largest source of natural textiles. Through a map-based cloning strategy, we isolated an industrially important lint fiber development gene (Li₃ ) that encodes an MYB-MIXTA-like transcription factor (MML) on chromosome D12 (GhMML4_D12). Virus-induced gene silencing or decreasing the expression of the GhMML4_D12 gene in n₂ NSM plants resulted in a significant reduction in epidermal cell prominence and lint fiber production. GhMML4_D12 is arranged in tandem with GhMML3, another MIXTA gene responsible for fuzz fiber development. These two very closely related MIXTA genes direct fiber initiation production in two specialized cell forms: lint and fuzz fibers. They may control the same metabolic pathways in different cell types. The MIXTAs expanded in Malvaceae during their evolution and produced a Malvaceae-specific family that regulates epidermal cell differentiation, different from the gene family that regulates leaf hair trichome development. Cotton has developed a unique transcriptional regulatory network for fiber development. Characterization of target genes regulating fiber production has provided insights into the molecular mechanisms underlying cotton fiber development and has allowed the use of genetic engineering to increase lint yield by inducing more epidermal cells to develop into lint rather than fuzz fibers.

Identification and Characterization of the Diverse Stress-Responsive R2R3-RMYB Transcription Factor from Hibiscus sabdariffa L

Various regulatory proteins play a fundamental role to manage the healthy plant growth under stress conditions. Differential display reverse transcriptase PCR and random amplification of cDNA ends (RACE) was used to explore the osmotic stress-responsive transcripts. We identified and characterized the salt stress-responsive R2R3 type RMYB transcription factor from Hibiscus sabdariffa which has an open reading frame of 690 bp, encoding 229 long chain amino acids. In silico analysis confirmed the conserved R2 and R3 domain as well as an NLS-1 localization site. The deduced amino acids of RMYB shared 83, 81, 80, 79, 72, 71, and 66% homology with Arabidopsis thaliana, Glycine max, Oryza sativa, Zea maize, Malus domestica, Populus tremula × Populus alba, and Medicago sativa specific MYB family, respectively. We observed the gene upregulation in stem, leaf, and root tissue in response to abiotic stress. Furthermore, RMYB gene was cloned into plant expression vector under CaMV35S promoter and transformed to Gossypium hirsutum: a local cotton cultivar. Overexpression of RMYB was observed in transgenic plants under abiotic stresses which further suggests its regulatory role in response to stressful conditions. The RMYB transcription factor-overexpressing in transgenic cotton plants may be used as potential agent for the development of stress tolerant crop cultivars.

Genome-wide characterization and expression analysis enables identification of abiotic stress-responsive MYB transcription factors in cassava (Manihot esculenta)

The myeloblastosis (MYB) transcription factor superfamily is the largest transcription factor family in plants, playing different roles during stress response. However, abiotic stress-responsive MYB transcription factors have not been systematically studied in cassava (Manihot esculenta), an important tropical tuber root crop. In this study, we used a genome-wide transcriptome analysis to predict 299 putative MeMYB genes in the cassava genome. Under drought and cold stresses, many MeMYB genes exhibited different expression patterns in cassava leaves, indicating that these genes might play a role in abiotic stress responses. We found that several stress-responsive MeMYB genes responded to abscisic acid (ABA) in cassava leaves. We characterize four MeMYBs, namely MeMYB1, MeMYB2, MeMYB4, and MeMYB9, as R2R3-MYB transcription factors. Furthermore, RNAi-driven repression of MeMYB2 resulted in drought and cold tolerance in transgenic cassava. Gene expression assays in wild-type and MeMYB2-RNAi cassava plants revealed that MeMYB2 may affect other MeMYBs as well as MeWRKYs under drought and cold stress, suggesting crosstalk between MYB and WRKY family genes under stress conditions in cassava.

Genome-wide identification and characterization of phospholipase C gene family in cotton (Gossypium spp.)

Phospholipase C (PLC) are important regulatory enzymes involved in several lipid and Ca²⁺-dependent signaling pathways. Previous studies have elucidated the versatile roles of PLC genes in growth, development and stress responses of many plants, however, the systematic analyses of PLC genes in the important fiber-producing plant, cotton, are still deficient. In this study, through genome-wide survey, we identified twelve phosphatidylinositol-specific PLC (PI-PLC) and nine non-specific PLC (NPC) genes in the allotetraploid upland cotton Gossypium hirsutum and nine PI-PLC and six NPC genes in two diploid cotton G. arboretum and G.raimondii, respectively. The PI-PLC and NPC genes of G. hirsutum showed close phylogenetic relationship with their homologous genes in the diploid cottons and Arabidopsis. Segmental and tandem duplication contributed greatly to the formation of the gene family. Expression profiling indicated that few of the PLC genes are constitutely expressed, whereas most of the PLC genes are preferentially expressed in specific tissues and abiotic stress conditions. Promoter analyses further implied that the expression of these PLC genes might be regulated by MYB transcription factors and different phytohormones. These results not only suggest an important role of phospholipase C members in cotton plant development and abiotic stress response but also provide good candidate targets for future molecular breeding of superior cotton cultivars.

Expression of the MYB transcription factor gene BplMYB46 affects abiotic stress tolerance and secondary cell wall deposition in Betula platyphylla

Plant MYB transcription factors control diverse biological processes, such as differentiation, development and abiotic stress responses. In this study, we characterized BplMYB46, an MYB gene from Betula platyphylla (birch) that is involved in both abiotic stress tolerance and secondary wall biosynthesis. BplMYB46 can act as a transcriptional activator in yeast and tobacco. We generated transgenic birch plants with overexpressing or silencing of BplMYB46 and subjected them to gain- or loss-of-function analysis. The results suggest that BplMYB46 improves salt and osmotic tolerance by affecting the expression of genes including SOD, POD and P5CS to increase both reactive oxygen species scavenging and proline levels. In addition, BplMYB46 appears to be involved in controlling stomatal aperture to reduce water loss. Overexpression of BplMYB46 increases lignin deposition, secondary cell wall thickness and the expression of genes in secondary cell wall formation. Further analysis indicated that BplMYB46 binds to MYBCORE and AC-box motifs and may directly activate the expression of genes involved in abiotic stress responses and secondary cell wall biosynthesis whose promoters contain these motifs. The transgenic BplMYB46-overexpressing birch plants, which have improved salt and osmotic stress tolerance, higher lignin and cellulose content and lower hemicellulose content than the control, have potential applications in the forestry industry.

Cotton GhMYB7 is predominantly expressed in developing fibers and regulates secondary cell wall biosynthesis in transgenic Arabidopsis

The secondary cell wall in mature cotton fibers contains over 90% cellulose with low quantities of xylan and lignin. However, little is known regarding the regulation of secondary cell wall biosynthesis in cotton fibers. In this study, we characterized an R2R3-MYB transcription factor, GhMYB7, in cotton. GhMYB7 is expressed at a high level in developing fibers and encodes a MYB protein that is targeted to the cell nucleus and has transcriptional activation activity. Ectopic expression of GhMYB7 in Arabidopsis resulted in small, curled, dark green leaves and also led to shorter inflorescence stems. A cross-sectional assay of basal stems revealed that cell wall thickness of vessels and interfascicular fibers was higher in transgenic lines overexpressing GhMYB7 than in the wild type. Constitutive expression of GhMYB7 in Arabidopsis activated the expression of a suite of secondary cell wall biosynthesis-related genes (including some secondary cell wall-associated transcription factors), leading to the ectopic deposition of cellulose and lignin. The ectopic deposition of secondary cell walls may have been initiated before the cessation of cell expansion. Moreover, GhMYB7 was capable of binding to the promoter regions of AtSND1 and AtCesA4, suggesting that GhMYB7 may function upstream of NAC transcription factors. Collectively, these findings suggest that GhMYB7 is a potential transcriptional activator, which may participate in regulating secondary cell wall biosynthesis of cotton fibers.

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

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

The cotton transcription factor TCP14 functions in auxin-mediated epidermal cell differentiation and elongation

Plant-specific TEOSINTE-BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factors play crucial roles in development, but their functional mechanisms remain largely unknown. Here, we characterized the cellular functions of the class I TCP transcription factor GhTCP14 from upland cotton (Gossypium hirsutum). GhTCP14 is expressed predominantly in fiber cells, especially at the initiation and elongation stages of development, and its expression increased in response to exogenous auxin. Induced heterologous overexpression of GhTCP14 in Arabidopsis (Arabidopsis thaliana) enhanced initiation and elongation of trichomes and root hairs. In addition, root gravitropism was severely affected, similar to mutant of the auxin efflux carrier PIN-FORMED2 (PIN2) gene. Examination of auxin distribution in GhTCP14-expressing Arabidopsis by observation of auxin-responsive reporters revealed substantial alterations in auxin distribution in sepal trichomes and root cortical regions. Consistent with these changes, expression of the auxin uptake carrier AUXIN1 (AUX1) was up-regulated and PIN2 expression was down-regulated in the GhTCP14-expressing plants. The association of GhTCP14 with auxin responses was also evidenced by the enhanced expression of auxin response gene IAA3, a gene in the AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) family. Electrophoretic mobility shift assays showed that GhTCP14 bound the promoters of PIN2, IAA3, and AUX1, and transactivation assays indicated that GhTCP14 had transcription activation activity. Taken together, these results demonstrate that GhTCP14 is a dual-function transcription factor able to positively or negatively regulate expression of auxin response and transporter genes, thus potentially acting as a crucial regulator in auxin-mediated differentiation and elongation of cotton fiber cells.

Transgenic wheat expressing Thinopyrum intermedium MYB transcription factor TiMYB2R-1 shows enhanced resistance to the take-all disease

The disease take-all, caused by the fungus Gaeumannomyces graminis, is one of the most destructive root diseases of wheat worldwide. Breeding resistant cultivars is an effective way to protect wheat from take-all. However, little progress has been made in improving the disease resistance level in commercial wheat cultivars. MYB transcription factors play important roles in plant responses to environmental stresses. In this study, an R2R3-MYB gene in Thinopyrum intermedium, TiMYB2R-1, was cloned and characterized. The gene sequence includes two exons and an intron. The expression of TiMYB2R-1 was significantly induced following G. graminis infection. An in vitro DNA binding assay proved that TiMYB2R-1 protein could bind to the MYB-binding site cis-element ACI. Subcellular localization assays revealed that TiMYB2R-1 was localized in the nucleus. TiMYB2R-1 transgenic wheat plants were generated, characterized molecularly, and evaluated for take-all resistance. PCR and Southern blot analyses confirmed that TiMYB2R-1 was integrated into the genomes of three independent transgenic wheat lines by distinct patterns and the transgene was heritable. Reverse transcription-PCR and western blot analyses revealed that TiMYB2R-1 was highly expressed in the transgenic wheat lines. Based on disease response assessments for three successive generations, the significantly enhanced resistance to take-all was observed in the three TiMYB2R-1-overexpressing transgenic wheat lines. Furthermore, the transcript levels of at least six wheat defence-related genes were significantly elevated in the TiMYB2R-1 transgenic wheat lines. These results suggest that engineering and overexpression of TiMYB2R-1 may be used for improving take-all resistance of wheat and other cereal crops.

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

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

Expression of the R2R3-MYB transcription factor TaMYB14 from Trifolium arvense activates proanthocyanidin biosynthesis in the legumes Trifolium repens and Medicago sativa

Proanthocyanidins (PAs) are oligomeric flavonoids and one group of end products of the phenylpropanoid pathway. PAs have been reported to be beneficial for human and animal health and are particularly important in pastoral agricultural systems for improved animal production and reduced greenhouse gas emissions. However, the main forage legumes grown in these systems, such as Trifolium repens and Medicago sativa, do not contain any substantial amounts of PAs in leaves. We have identified from the foliar PA-accumulating legume Trifolium arvense an R2R3-MYB transcription factor, TaMYB14, and provide evidence that this transcription factor is involved in the regulation of PA biosynthesis in legumes. TaMYB14 expression is necessary and sufficient to up-regulate late steps of the phenylpropanoid pathway and to induce PA biosynthesis. RNA interference silencing of TaMYB14 resulted in almost complete cessation of PA biosynthesis in T. arvense, whereas Nicotiana tabacum, M. sativa, and T. repens plants constitutively expressing TaMYB14 synthesized and accumulated PAs in leaves up to 1.8% dry matter. Targeted liquid chromatography-multistage tandem mass spectrometry analysis identified foliar PAs up to degree of polymerization 6 in leaf extracts. Hence, genetically modified M. sativa and T. repens plants expressing TaMYB14 provide a viable option for improving animal health and mitigating the negative environmental impacts of pastoral animal production systems.

Expression of a wheat MYB gene in transgenic tobacco enhances resistance to Ralstonia solanacearum, and to drought and salt stresses

MYB transcription factors play diverse roles in plant growth, developmental processes and stress responses. A full-length cDNA sequence of a MYB gene, namely TaPIMP1, was isolated from wheat (Triticum aestivum L.). The TaPIMP1 transcript level was significantly up-regulated by inoculation with a fungal pathogen Bipolaris sorokiniana and by drought treatment. TaPIMP1 encodes the MYB protein TaPIMP1 consisting of 323 amino acids. TaPIMP1 contains two MYB DNA binding domains (R2, R3), two putative nuclear localization sites and two putative transcription activation domains. TaPIMP1 is a new member of the R2R3-MYB transcription factor subfamily. Transient expression in onion epidermal cells of GFP fused with TaPIMP1 proved that subcellular localization of TaPIMP1 occurred in the nucleus. The TaPIMP1 gene was transferred into tobacco (Nicotiana tabacum L.) cultivar W38 by Agrobacterium-mediated transformation. After screening through PCR and RT-PCR analyses, transgenic tobacco lines expressing TaPIMP1 were identified and evaluated for pathogen resistance, and drought and salt tolerance. Compared to untransformed tobacco host plants, TaPIMP1 expressing plants displayed significantly enhanced resistance to Ralstonia solanacearum and exhibited improved tolerances to drought and salt stresses. In these transgenic lines, the activities of phenylalanine ammonia-lyase (PAL) and superoxide dismutase (SOD) were significantly increased relative to wild-type tobacco plants. The results suggested that the wheat R2R3-MYB transcription factor plays an important role in modulating responses to biotic and abiotic stresses.

Marker-assisted breeding as next-generation strategy for genetic improvement of productivity and quality: can it be realized in cotton?

The dawdling development in genetic improvement of cotton with conventional breeding program is chiefly due to lack of complete knowledge on and precise manipulation of fiber productivity and quality. Naturally available cotton continues to be a resource for the upcoming breeding program, and contemporary technologies to exploit the available natural variation are outlined in this paper for further improvement of fiber. Particularly emphasis is given to application, obstacles, and perspectives of marker-assisted breeding since it appears to be more promising in manipulating novel genes that are available in the cotton germplasm. Deployment of system quantitative genetics in marker-assisted breeding program would be essential to realize its role in cotton. At the same time, role of genetic engineering and in vitro mutagenesis cannot be ruled out in genetic improvement of cotton.

GhMYB25-like: a key factor in early cotton fibre development

MYB transcription factors have been implicated in regulation of the development of ovule epidermal cells into the elongated seed fibres of cotton. An R2R3 MYB, GhMYB25-like, identified from its reduced expression in a fibreless mutant of cotton (Xu142 fl), is here shown to play a key role in the very early stages of fibre cell differentiation. A GhMYB25-like promoter-GUS construct was expressed predominantly in the epidermal layers of cotton ovules before anthesis (-3days post-anthesis, dpa), increasing in expression in 0-dpa ovules, primarily in those epidermal cells expanding into fibres, and then in elongating fibres at +3dpa, declining thereafter. This was consistent with GhMYB25-like transcript abundance during fibre development. RNA interference suppression of GhMYB25-like resulted in cotton plants with fibreless seeds, but normal trichomes elsewhere, phenocopying the Xu142 fl mutant. Like Xu142 fl these plants had reduced expression of the fibre-expressed MYBs, GhMYB25 and GhMYB109, indicating that GhMYB25-like is upstream from those MYBs. This hierarchy was supported by the absence of any change in transcript level of GhMYB25-like in GhMYB25- and GhMYB109-silenced transgenic lines. Transgenic cotton with an additional copy of the native gene had elevated expression of GhMYB25-like in ovules, but no obvious increase in fibre initials, suggesting that there may be other factors that interact with GhMYB25-like to differentiate epidermal cells into fibre cells.

Polyploidization altered gene functions in cotton (Gossypium spp.)

Cotton (Gossypium spp.) is an important crop plant that is widely grown to produce both natural textile fibers and cottonseed oil. Cotton fibers, the economically more important product of the cotton plant, are seed trichomes derived from individual cells of the epidermal layer of the seed coat. It has been known for a long time that large numbers of genes determine the development of cotton fiber, and more recently it has been determined that these genes are distributed across At and Dt subgenomes of tetraploid AD cottons. In the present study, the organization and evolution of the fiber development genes were investigated through the construction of an integrated genetic and physical map of fiber development genes whose functions have been verified and confirmed. A total of 535 cotton fiber development genes, including 103 fiber transcription factors, 259 fiber development genes, and 173 SSR-contained fiber ESTs, were analyzed at the subgenome level. A total of 499 fiber related contigs were selected and assembled. Together these contigs covered about 151 Mb in physical length, or about 6.7% of the tetraploid cotton genome. Among the 499 contigs, 397 were anchored onto individual chromosomes. Results from our studies on the distribution patterns of the fiber development genes and transcription factors between the At and Dt subgenomes showed that more transcription factors were from Dt subgenome than At, whereas more fiber development genes were from At subgenome than Dt. Combining our mapping results with previous reports that more fiber QTLs were mapped in Dt subgenome than At subgenome, the results suggested a new functional hypothesis for tetraploid cotton. After the merging of the two diploid Gossypium genomes, the At subgenome has provided most of the genes for fiber development, because it continues to function similar to its fiber producing diploid A genome ancestor. On the other hand, the Dt subgenome, with its non-fiber producing D genome ancestor, provides more transcription factors that regulate the expression of the fiber genes in the At subgenome. This hypothesis would explain previously published mapping results. At the same time, this integrated map of fiber development genes would provide a framework to clone individual full-length fiber genes, to elucidate the physiological mechanisms of the fiber differentiation, elongation, and maturation, and to systematically study the functional network of these genes that interact during the process of fiber development in the tetraploid cottons.

The MYB transcription factor GhMYB25 regulates early fibre and trichome development

Little is still known about the developmental control of the long seed coat trichomes of cotton (Gossypium hirsutum L.). In Arabidopsis, leaf trichome initiation is regulated by a group of well-defined transcription factors that includes MYB and homeodomain types. Many MYBs are expressed in fibres, but their roles in fibre development remain unclear. We analysed the function of one MYB transcription factor, GhMYB25, identified from transcriptome comparisons between wild-type and fibreless cotton mutants. A GhMYB25 promoter-GUS construct in transgenic cotton was expressed in the epidermis of ovules, developing fibre initials and fibres, in the trichomes of a number of tissues including leaves, stems and petals, as well as in the anthers, pollen and the epidermal layers of roots and root initials, but not in root hairs. Cotton is an allotetraploid with two very similar GhMYB25 genes that were silenced by a single RNAi construct. GhMYB25-silenced cotton showed alterations in the timing of rapid fibre elongation, resulting in short fibres, dramatic reductions in trichomes on other parts of the plant, and reductions in seed production. Reciprocal crosses between transgenic and non-transgenic plants indicated that pollen and ovule viability per se were not disrupted. Ectopic over-expression of GhMYB25 had more subtle impacts, with increases in cotton fibre initiation and leaf trichome number. High expression appeared to adversely affect fertility. Our results provide convincing evidence for a role of GhMYB25, like other MIXTA-like MYBS, in regulating specialized outgrowths of epidermal cells, including, in this case, cotton fibres.

The HD-Zip IV gene GaHOX1 from cotton is a functional homologue of the Arabidopsis GLABRA2

Most of the plant homeodomain-containing proteins play important roles in organ patterning and development, and Arabidopsis GLABRA2 (GL2), a member of the class IV homeodomain-leucine zipper (HD-ZIP) proteins, is a trichome and non-root hair cell regulator. Here we report the analysis of two cotton homeodomain-containing proteins, GaHOX1 and GaHOX2, isolated from the diploid cotton Gossypium arboreum. Both GaHOX1 and GaHOX2 belong to the class IV HD-ZIP family. When expressed under the control of the GL2 promoter, GaHOX1 rescued trichome development of an Arabidopsis glabrous mutant of gl2-2 (SALK_130213), whereas GaHOX2 did not. On the other hand, expression of GaHOX1 with a Cauliflower mosaic virus (CaMV) 35S promoter in the wild-type Arabidopsis plants suppressed the trichome development just as the GL2 ectopic expression. Expression analysis by Northern, RT-PCR and in situ hybridization indicated that GaHOX1 is predominantly expressed in cotton fiber cells at early developmental stages, consistent with its putative role in regulating cotton fiber development, while GaHOX2 is expressed in both fiber and other ovular tissues, including outer and inner integuments. Our results suggest that GaHOX1 is a functional homolog of GL2 in plant trichome development.

Gene-rich islands for fiber development in the cotton genome

Cotton fiber is an economically important seed trichome and the world's leading natural fiber used in the manufacture of textiles. As a step toward elucidating the genomic organization and distribution of gene networks responsible for cotton fiber development, we investigated the distribution of fiber genes in the cotton genome. Results revealed the presence of gene-rich islands for fiber genes with a biased distribution in the tetraploid cotton (Gossypium hirsutum L.) genome that was also linked to discrete fiber developmental stages based on expression profiles. There were 3 fiber gene-rich islands associated with fiber initiation on chromosome 5, 3 islands for the early to middle elongation stage on chromosome 10, 3 islands for the middle to late elongation stage on chromosome 14, and 1 island on chromosome 15 for secondary cell wall deposition, for a total of 10 fiber gene-rich islands. Clustering of functionally related gene clusters in the cotton genome displaying similar transcriptional regulation indicates an organizational hierarchy with significant implications for the genetic enhancement of particular fiber quality traits. The relationship between gene-island distribution and functional expression profiling suggests for the first time the existence of functional coupling gene clusters in the cotton genome.

Cotton (Gossypium spp.) R2R3-MYB transcription factors SNP identification, phylogenomic characterization, chromosome localization, and linkage mapping

R2R3-MYB transcription factors of plants are involved in the regulation of trichome length and density. Several of them are differentially expressed during initiation and elongation of cotton fibers. We report sequence phylogenomic characterization of the six MYB genes, their chromosomal localization, and linkage mapping via SNP marker in AD-genome cotton (2n = 52). Phylogenetic grouping and comparison to At- and Dt-genome putative ancestral diploid species of allotetraploid cotton facilitated differentiation between genome-specific polymorphisms (GSPs) and marker-suitable locus-specific polymorphisms (LSPs). The SNP frequency averaged one per 77 bases overall, and one per 106 and 30 bases in coding and non-coding regions, respectively. SNP-based multivariate relationships conformed to independent evolution of the six MYB homoeologous loci in the four tetraploid species. Nucleotide diversity analysis indicated that the six MYB loci evolved more quickly in the Dt- than At-genome. The greater variation in the Dt-D genome comparisons than that in At-A genome comparisons showed no significant bias among synonymous substitution, non-synonymous substitution, and nucleotide change in non-coding regions. SNPs were concordantly mapped by deletion analysis and linkage mapping, which confirmed their value as candidate gene markers and indicated the reliability of the SNP discovery strategy in tetraploid cotton species. We consider that these SNPs may be useful for genetic dissection of economically important fiber and yield traits because of the role of these genes in fiber development.

Identification of target genes for a MYB-type anthocyanin regulator in Gerbera hybrida

Genetic modification of the flavonoid pathway has been used to produce novel colours and colour patterns in ornamental plants as well as to modify the nutritional and pharmaceutical properties of food crops. It has been suggested that co-ordinate control of multiple steps of the pathway with the help of regulatory genes would lead to a more predictable control of metabolic flux. Regulation of anthocyanin biosynthesis has been studied in a common ornamental plant, Gerbera hybrida (Asteraceae). An R2R3-type MYB factor, GMYB10, shares high sequence similarity and is phylogenetically grouped together with previously characterized regulators of anthocyanin pigmentation. Ectopic expression of GMYB10 leads to strongly enhanced accumulation of anthocyanin pigments as well as to an altered pigmentation pattern in transgenic gerbera plants. Anthocyanin analysis indicates that GMYB10 specifically induces cyanidin biosynthesis in undifferentiated callus and in vegetative tissues. Furthermore, in floral tissues enhanced pelargonidin production is detected. Microarray analysis using the gerbera 9K cDNA array revealed a highly predicted set of putative target genes for GMYB10 including new gene family members of both early and late biosynthetic genes of the flavonoid pathway. However, completely new candidate targets, such as a serine carboxypeptidase-like gene as well, as two new MYB domain factors, GMYB11 and GMYB12, whose exact function in phenylpropanoid biosynthesis is not clear yet, were also identified.

HlMyb3, a putative regulatory factor in hop (Humulus lupulus L.), shows diverse biological effects in heterologous transgenotes

A hop-specific cDNA library from glandular tissue-enriched hop cones was screened for Myb transcription factors. cDNA encoding for R2R3 Myb, designated HlMyb3, was cloned and characterized. According to the amino acid (aa) sequence, HlMyb3 shows the highest homology to GhMyb5 from cotton and is unrelated to the previously characterized HlMyb1 from the hop. Southern blot analyses indicated that HlMyb3 is a unique gene, which was detected in various Humulus lupulus cultivars, but not in Humulus japonicus. Reverse transcription and real-time PCR revealed the highest levels of HlMyb3 mRNA in hop cones at a late stage of maturation and in colored petiole epidermis, while the lowest levels were observed in hop flowers. Two alternative open reading frames starting in the N-terminal domain of HlMyb3, encoding for proteins having 269 and 265 amino acids with apparent molecular masses of 30.3 and 29.9 kDa, respectively, were analyzed as transgenes that were overexpressed in Arabidopsis thaliana, Nicotiana benthamiana, and Petunia hybrida plants. Transformation with the longer 269 aa variant designated l-HlMyb3 led to a flowering delay and to a strong inhibition of seed germination in A. thaliana. Nearly complete flower sterility, dwarfing, and leaf curling of P. hybrida and N. benthamiana l-HlMyb3 transgenotes were noted. On the contrary, the shorter 265-aa-encoding s-HlMyb3 transgene led in A. thaliana to the stimulation of initial seed germination, to fast initiation of the lateral roots, and to quite specific branching phenotypes with many long lateral stems formed at angles near 90 degrees . Limited plant sterility but growth stimulation and rather branched phenotypes were evident for s-HlMyb3 transgenotes of P. hybrida and N. benthamiana. It was found that both HlMyb3 transgenes interfere in the accumulation and composition of flavonol glycosides and phenolic acids in transformed plants. These effects on heterologous transgenotes suggest that the HlMyb3 gene may influence hop morphogenesis, as well as metabolome composition during lupulin gland maturation.

Transcriptome profiling, sequence characterization, and SNP-based chromosomal assignment of the EXPANSIN genes in cotton

The knowledge of biological significance associated with DNA markers is very limited in cotton. SNPs are potential functional marker to tag genes of biological importance. Plant expansins are a group of extracellular proteins that directly modify the mechanical properties of cell walls, enable turgor-driven cell extension, and likely affect length and quality of cotton fibers. Here, we report the expression profiles of EXPANSIN transcripts during fiber elongation and the discovery of SNP markers, assess the SNP characteristics, and localize six EXPANSIN A genes to chromosomes. Transcriptome profiling of cotton fiber oligonucleotide microarrays revealed that seven EXPANSIN transcripts were differentially expressed when there was parallel polar elongation during morphogenesis at early stage of fiber development, suggesting that major and minor isoforms perform discrete functions during polar elongation and lateral expansion. Ancestral and homoeologous relationships of the six EXPANSIN A genes were revealed by phylogenetic grouping and comparison to extant A- and D-genome relatives of contemporary AD-genome cottons. The average rate of SNP per nucleotide was 2.35% (one SNP per 43 bp), with 1.74 and 3.99% occurring in coding and noncoding regions, respectively, in the selected genotypes. An unequal evolutionary rate of the EXPANSIN A genes at the subgenome level of tetraploid cotton was recorded. Chromosomal locations for each of the six EXPANSIN A genes were established by gene-specific SNP markers. Results revealed a strategy for discovering SNP markers in a polyploidy species like cotton. These markers could be useful to associate candidate genes with the complex fiber traits in MAS.

Evolutionary dynamics of Waxy and the origin of hexaploid Spartina species (Poaceae)

We investigated the evolutionary dynamics of duplicated copies of the granule-bound starch synthase I gene (GBSSI or Waxy) within polyploid Spartina species. Molecular cloning, sequencing, and phylogenetic analyses revealed incongruences between the expected species phylogeny and the inferred gene trees. Some genes within species were more divergent than expected from ploidy level alone, suggesting the existence of paralogous sets of Waxy loci in Spartina. Phylogenetic analyses indicate that this paralogy originated from a duplication that occurred prior to the divergence of Spartina from other Chloridoideae. Gene tree topologies revealed three divergent homoeologous sequences in the hexaploid S. alterniflora that are consistent with the proposal of an allopolyploid origin of the hexaploid clade. Waxy sequences differ in insertion-deletion events in introns, which may be used to diagnose gene copies. Both paralogous and homoeologous coding regions appear to evolving under selective constraints.

Innovation in anti-herbivore defense systems during neopolypoloidy - the functional consequences of instantaneous speciation

Allopolyploid hybridization instantly merges two differentially adapted genomes into one individual. Allopolyploids are often evolutionarily successful, undergoing adaptive radiations despite the genetic and physiological challenges of merging genomes. We examine a suite of induced herbivore resistance traits in three independent lines of the synthetic allopolyploid Nicotianaxmierata (Nma) and its parent species, N. miersii (Nmi) and N. attenuata (Na), to determine how a complex polygenetic adaptation fares during the early stages of neoallopolyploid formation. All species responded to Manduca sexta oral secretions (OS) with a temporally prolonged jasmonate (JA) burst. In one parent (Na), the JA burst was additionally amplified and associated with the elicitation of direct and indirect defenses. In the other parent (Nmi), OS neither amplified the JA burst nor elicited defense responses, although applied MeJA confirmed the inducibility of the defense responses. All lines of Nma retained enough aspects of Na's JA signaling to recognize OS and to accumulate sufficient direct defenses to impair the growth of Manduca larvae. Most defense-related metabolites were retained in Nma even if inherited from only one parent; however, OS-elicited volatiles, trypsin protease inhibitors (TPIs) and chlorogenic acid were lost in some lines, even though MeJA treatment elicited similar responses in all lines. Herbivore defense systems are flexibly inherited in allopolyploids, causing individuals to diverge over only a few generations; for example, line 1 of Nma could not produce TPIs after OS elicitation, whereas lines 2 and 3 could. This flexible integration of defense signaling systems with a diversity of elicited responses may explain why adaptive radiations are commonly found in allopolyploid lineages.

A global assembly of cotton ESTs

Approximately 185,000 Gossypium EST sequences comprising >94,800,000 nucleotides were amassed from 30 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including drought stress and pathogen challenges. These libraries were derived from allopolyploid cotton (Gossypium hirsutum; A(T) and D(T) genomes) as well as its two diploid progenitors, Gossypium arboreum (A genome) and Gossypium raimondii (D genome). ESTs were assembled using the Program for Assembling and Viewing ESTs (PAVE), resulting in 22,030 contigs and 29,077 singletons (51,107 unigenes). Further comparisons among the singletons and contigs led to recognition of 33,665 exemplar sequences that represent a nonredundant set of putative Gossypium genes containing partial or full-length coding regions and usually one or two UTRs. The assembly, along with their UniProt BLASTX hits, GO annotation, and Pfam analysis results, are freely accessible as a public resource for cotton genomics. Because ESTs from diploid and allotetraploid Gossypium were combined in a single assembly, we were in many cases able to bioinformatically distinguish duplicated genes in allotetraploid cotton and assign them to either the A or D genome. The assembly and associated information provide a framework for future investigation of cotton functional and evolutionary genomics.

Developmental and gene expression analyses of a cotton naked seed mutant

Cotton fiber development is a fundamental biological phenomenon, yet the molecular basis of fiber cell initiation is poorly understood. We examined molecular and cellular events of fiber cell development in the naked seed mutant (N1N1) and its isogenic line of cotton (Gossypium hirsutum L. cv. Texas Marker-1, TM-1). The dominant mutation not only delayed the process of fiber cell formation and elongation but also reduced the total number of fiber cells, resulting in sparsely distributed short fibers. Gene expression changes in TM-1 and N1N1 mutant lines among four tissues were analyzed using spotted cotton oligo-gene microarrays. Using the Arabidopsis genes, we selected and designed approximately 1,334 70-mer oligos from a subset of cotton fiber ESTs. Statistical analysis of the microarray data indicates that the number of significantly differentially expressed genes was 856 in the leaves compared to the ovules (3 days post-anthesis, DPA), 632 in the petals relative to the ovules (3 DPA), and 91 in the ovules at 0 DPA compared to 3 DPA, all in TM-1. Moreover, 117 and 30 genes were expressed significantly different in the ovules at three and 0 DPA, respectively, between TM-1 and N1N1. Quantitative RT-PCR analysis of 23 fiber-associated genes in seven tissues including ovules, fiber-bearing ovules, fibers, and non-fiber tissues in TM-1 and N1N1 indicates a mode of temporal regulation of the genes involved in transcriptional and translational regulation, signal transduction, and cell differentiation during early stages of fiber development. Suppression of the fiber-associated genes in the mutant may suggest that the N1N1 mutation disrupts temporal regulation of gene expression, leading to a defective process of fiber cell elongation and development.

Genome evolution of allopolyploids: a process of cytological and genetic diploidization

Allopolyploidy is a prominent mode of speciation in higher plants. Due to the coexistence of closely related genomes, a successful allopolyploid must have the ability to invoke and maintain diploid-like behavior, both cytologically and genetically. Recent studies on natural and synthetic allopolyploids have raised many discrepancies. Most species have displayed non-Mendelian behavior in the allopolyploids, but others have not. Some species have demonstrated rapid genome changes following allopolyploid formation, while others have conserved progenitor genomes. Some have displayed directed, non-random genome changes, whereas others have shown random changes. Some of the genomic changes have appeared in the F1 hybrids, which have been attributed to the union of gametes from different progenitors, while other changes have occurred during or after genome doubling. Although these observations provide significant novel insights into the evolution of allopolyploids, the overall mechanisms of the event are still elusive. It appears that both genetic and epigenetic operations are involved in the diploidization process of allopolyploids. Overall, genetic and epigenetic variations are often associated with the activities of repetitive sequences and transposon elements. Specifically, genomic sequence elimination and chromosome rearrangement are probably the major forces guiding cytological diploidization. Gene non-functionalization, sub-functionalization, neo-functionalization, as well as other kinds of epigenetic modifications, are likely the leading factors promoting genetic diploidization.

Control of plant trichome development by a cotton fiber MYB gene

Cotton (Gossypium spp) plants produce seed trichomes (cotton fibers) that are an important commodity worldwide; however, genes controlling cotton fiber development have not been characterized. In Arabidopsis thaliana the MYB gene GLABRA1 (GL1) is a central regulator of trichome development. Here, we show that promoter of a cotton fiber gene, RD22-like1 (RDL1), contains a homeodomain binding L1 box and a MYB binding motif that confer trichome-specific expression in Arabidopsis. A cotton MYB protein GaMYB2/Fiber Factor 1 transactivated the RDL1 promoter both in yeast and in planta. Real-time PCR and in situ analysis showed that GaMYB2 is predominantly expressed early in developing cotton fibers. After transferring into Arabidopsis, GL1::GaMYB2 rescued trichome formation of a gl1 mutant, and interestingly, 35S::GaMYB2 induced seed-trichome production. We further demonstrate that the first intron of both GL1 and GaMYB2 plays a role in patterning trichomes: it acts as an enhancer in trichome and a repressor in nontrichome cells, generating a trichome-specific pattern of MYB gene expression. Disruption of a MYB motif conserved in intron 1 of GL1, WEREWOLF, and GaMYB2 genes affected trichome production. These results suggest that cotton and Arabidopsis use similar transcription factors for regulating trichomes and that GaMYB2 may be a key regulator of cotton fiber development.

Functional genomics of cell elongation in developing cotton fibers

Cotton fibers are single-celled seed trichomes of major economic importance. Factors that regulate the rate and duration of cell expansion control fiber morphology and important agronomic traits. For genetic characterization of rapid cell elongation in cotton fibers, approximately 14,000 unique genes were assembled from 46,603 expressed sequence tags (ESTs) from developmentally staged fiber cDNAs of a cultivated diploid species ( Gossypium arboreum L.). Conservatively, the fiber transcriptome represents 35-40% of the genes in the cotton genome. In silico expression analysis revealed that rapidly elongating fiber cells exhibit significant metabolic activity, with the bulk of gene transcripts, represented by three major functional groups - cell wall structure and biogenesis, the cytoskeleton and energy/carbohydrate metabolism. Oligonucleotide microarrays revealed dynamic changes in gene expression between primary and secondary cell wall biogenesis showing that fiber genes in the dbEST are highly stage-specific for cell expansion - a conclusion supported by the absence of known secondary cell wall-specific genes from our fiber dbEST. During the developmental switch from primary to secondary cell wall syntheses, 2553 "expansion-associated" fiber genes are significantly down regulated. Genes (81) significantly up-regulated during secondary cell wall synthesis are involved in cell wall biogenesis and energy/carbohydrate metabolism, which is consistent with the stage of cellulose synthesis during secondary cell wall modification in developing fibers. This work provides the first in-depth view of the genetic complexity of the transcriptome of an expanding cell, and lays the groundwork for studying fundamental biological processes in plant biology with applications in agricultural biotechnology.

Ordered origin of the typical two- and three-repeat Myb genes

Myb domain proteins contain a conserved DNA-binding domain composed of one to four conserved repeat motifs. In animals, Myb proteins are encoded by a small gene family and commonly contain three repeat motifs (R1R2R3); whereas, plant Myb proteins are encoded by a very large and diverse gene family in which a motif containing two repeats (R2R3) is the most common. In contrast to the conservation in the Myb domain, other regions of Myb proteins are highly variable. To explore the evolutionary origin of Myb genes, we cloned and sequenced Myb domains from maize and sorghum, and conducted a comprehensive phylogenetic analysis of Myb genes. The results indicate that the origins of individual Myb repeats are strikingly distinct, and that the R2 repeat has evolved more slowly than the R1 and R3 repeats. However, it is not clear which repeat is the most ancient one. The evidence also suggests that R2R3 and R1R2R3 Myb genes co-existed in eukaryotes before the divergence of plants and animals. Based on our results, we propose that R1R2R3 Myb genes were derived from R2R3 Myb genes by gain of the R1 repeat through an ancient intragenic duplication; this gain model is more parsimonious than the previous proposal that R2R3 Myb genes were derived from R1R2R3 Mybs by loss of the R1 repeat. A separate group of diverse non-typical Myb proteins exhibits a polyphyletic origin and a complex evolutionary pattern. Finally, a small group of ancient Myb paralogs prior to the amplification of current Myb genes is identified. Together, these results support a new model for the ordered evolution of Myb gene family.

An ATP-binding cassette transporter GhWBC1 from elongating cotton fibers

We have isolated a cDNA (GhWBC1) from cotton (Gossypium hirsutum) that encodes an ATP-binding cassette transporter of the WBC (white/brown complex) subfamily. Members of this subfamily are half-sized transporters and are reported to mediate lipid and drug excretion in human (Homo sapiens). GhWBC1 is highly expressed in developing fiber cells, but transcripts were also detectable in other tissues except roots. The transcript level peaked in rapidly expanding fibers from 5 to 9 DPA and then decreased. The GhWBC1 expression was weak in fiber cells of an li (ligon-lintless) mutant, which is defective in fiber cell elongation. These data indicate that GhWBC1 gene expression correlates with cotton fiber elongation. Transient expression of enhanced green fluorescence protein-GhWBC1 fusion protein in onion (Allium cepa) epidermal cells revealed plasma membrane localization. The GhWBC1 cDNA driven by a constitutive 35S promoter was introduced into Arabidopsis. About 13% of the transformants produced short siliques (SSs), whereas others had normal siliques (long siliques [LSs]). In siliques of SS lines, most embryos were severely shriveled, and only several seeds per silique could be found at maturity. The transgene expression level was higher in SS lines than in LS lines. Expression of AtWBC11, the closest homolog of GhWBC1 in Arabidopsis, was not altered in either SS or LS transgenic plants examined. These data suggest that GhWBC1 interferes with substance translocation that is required for Arabidopsis seed and silique development. Characterization of Arabidopsis WBC members, particularly AtWBC11, will help to dissect the role of GhWBC1 in cotton fiber development and elongation.

Genes duplicated by polyploidy show unequal contributions to the transcriptome and organ-specific reciprocal silencing

Most eukaryotes have genomes that exhibit high levels of gene redundancy, much of which seems to have arisen from one or more cycles of genome doubling. Polyploidy has been particularly prominent during flowering plant evolution, yielding duplicated genes (homoeologs) whose expression may be retained or lost either as an immediate consequence of polyploidization or on an evolutionary timescale. Expression of 40 homoeologous gene pairs was assayed by cDNA-single-stranded conformation polymorphism in natural (1- to 2-million-yr-old) and synthetic tetraploid cotton (Gossypium) to determine whether homoeologous gene pairs are expressed at equal levels after polyploid formation. Silencing or unequal expression of one homoeolog was documented for 10 of 40 genes examined in ovules of Gossypium hirsutum. Assays of homoeolog expression in 10 organs revealed variable expression levels and silencing, depending on the gene and organ examined. Remarkably, silencing and biased expression of some gene pairs are reciprocal and developmentally regulated, with one homoeolog showing silencing in some organs and the other being silenced in other organs, suggesting rapid subfunctionalization. Duplicate gene expression was examined in additional natural polyploids to characterize the pace at which expression alteration evolves. Analysis of a synthetic tetraploid revealed homoeolog expression and silencing patterns that sometimes mirrored those of the natural tetraploid. Both long-term and immediate responses to polyploidization were implicated. Data suggest that some silencing events are epigenetically induced during the allopolyploidization process.