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

Function analysis of a cotton R2R3 MYB transcription factor GhMYB3 in regulating plant trichome development

Cotton is an important fibre-producing crop. Cotton fibres consist of highly elongated trichomes derived from the ovule. To improve the quality of cotton, it is necessary to identify the genes regulating fibre development. GhMYB3 was identified through bioinfomatic analysis and introduced to Arabidopsis and cotton to observe the phenotype. Protein inteaction and promoter bingding assays were conducted to explore the role of GhMYB3 in trichome fibre growth. Cotton fibre development might share a similar regulatory mechanism to Arabidopsis leaf trichomes, which is determined by the essential regulatory complex, MYB-bHLH-WD40. The GL1-like R2R3 MYB transcription factor GhMYB3 interacts with the AtGL3 protein involved in Arabidopsis trichome development. Ectopic expression of GhMYB3 could rescue the glabrous phenotype of the Arabidopsis gl1 mutant and produced more ectopic trichomes on inflorescence stems and floral organs, confirming its orthologous function in plant trichome development. The expression of GhMYB3 increased in response to exogenous gibberellin (GA₃ ), auxin (IAA) and methyl jasmonate (MeJA). Overexpression of this gene in cotton leads to a slight increase in fibre length and lint percentage, possibly by activating the transcription of its downstream gene GhRDL1 or other fibre-related genes. The results increase our understanding of the key role of GhMYB3 in positively controlling plant trichome development, and this gene could be a potential target for molecular breeding in cotton.

Identification and profiling of upland cotton microRNAs at fiber initiation stage under exogenous IAA application

Background

Cotton is the most essential textile crop worldwide, and phytohormones are critical for cotton fiber development. One example is the role of auxin in fiber initiation, but we know little molecular basis. MicroRNAs (miRNAs) have a significant function in cotton development; nevertheless their role in fiber initiation remains unclear. Here, exogenous IAA was applied to cotton plant before anthesis. Utilizing small RNA sequencing, the mechanism underlying miRNA-mediated regulation of fiber initiation under exogenous IAA treatment was investigated.

Results

With exogenous IAA application, the endogenous IAA and GA contents of IAA treated (IT) ovules were higher than control (CK) ovules at the fiber initiation stage, while endogenous ABA content was lower in IT than CK. Using scanning electron microscopy, we found the fiber number and size were significantly promoted in IT at 0 DPA. Fiber quality analysis showed that fiber length, uniformity, strength, elongation, and micronaire of IT were higher than CK, though not statistically significant, while lint percent was significantly higher in IT. We generated six small RNA libraries using − 3, 0, and 3 DPA ovules of IT and CK, and identified 58 known miRNAs and 83 novel miRNAs together with the target genes. The differential expressed miRNAs number between IT and CK at − 3, 0, 3 DPA was 34, 16 and 24, respectively. Gene ontology and KEGG pathway enrichment analyses for the target genes of the miRNAs expressed in a differential manner showed that they were significantly enriched in 30 terms and 8 pathways. QRT-PCR for those identified miRNAs and the target genes related to phytohormones and fiber development was performed, and results suggested a potential role of these miRNAs in fiber initiation.

Conclusions

The exogenous IAA application affected the relative phytohormone contents in ovule and promoted fiber initiation in cotton. Identification and profiling of miRNAs and their targets at the fiber initiation stage provided insights for miRNAs’ regulation function of fiber initiation. These findings not only shed light on the regulatory network of fiber growth but also offer clues for cotton fiber amelioration strategies in cotton.

Electronic supplementary material

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

Characterization of bHLH/HLH genes that are involved in brassinosteroid (BR) signaling in fiber development of cotton (Gossypium hirsutum)

BACKGROUND

Basic helix-loop-helix/helix-loop-helix (bHLH/HLH) transcription factors play important roles in plant development. Many reports have suggested that bHLH/HLH proteins participate in brassinosteroid (BR) hormone signaling pathways to promote cell elongation. Cotton fibers are single-cells and derived from seed surface. To explore the roles of bHLH/HLH proteins in cotton fiber development progress by modulating BR signaling pathway, we performed a systematic analysis of the bHLH/HLH gene family in upland cotton (Gossypium hirsutum) genome.

RESULTS

In this study, we identified 437 bHLH/HLH genes in upland cotton (G. hirsutum) genome. Phylogenetic analysis revealed that GhbHLH/HLH proteins were split into twenty six clades in the tree. These GhbHLH/HLH genes are distributed unevenly in different chromosomes of cotton genome. Segmental duplication is the predominant gene duplication event and the major contributor for amplification of GhbHLH/HLH gene family. The GhbHLH/HLHs within the same group have conserved exon/intron pattern and their encoding proteins show conserved motif composition. Based on transcriptome data, we identified 77 GhbHLH/HLH candidates that are expressed at relatively high levels in cotton fibers. As adding exogenous BR (brassinolide, BL) or brassinazole (Brz, a BR biosynthesis inhibitor), expressions of these GhbHLH/HLH genes were up-regulated or down-regulated in cotton fibers. Furthermore, overexpression of GhbHLH282 (one of the BR-response genes) in Arabidopsis not only promoted the plant growth, but also changed plant response to BR signaling.

CONCLUSION

Collectively, these data suggested that these GhbHLH/HLH genes may participate in BR signaling transduction during cotton fiber development. Thus, our results may provide a valuable reference data as the basis for further studying the roles of these bHLH/HLH genes in cotton fiber development.

Transcriptomic profiling of developing fiber in levant cotton (Gossypium herbaceum L.)

Cotton (Gossypium spp.) is an imperative economic crop of the globe due to its natural textile fiber. Molecular mechanisms of fiber development have been greatly revealed in allotetraploid cotton but remained unexplored in Gossypium herbaceum. G. herbaceum can withstand the rigors of nature like drought and pests but produce coarse lint. This undesirable characteristic strongly needs the knowledge of fiber development at molecular basis. The present study reported the transcriptome sequence of the developing fiber of G. herbaceum on pyrosequencing and its analysis. About 1.38 million raw and 1.12 million quality trimmed reads were obtained followed by de novo assembly-generated 20,125 unigenes containing 14,882 coding sequences (CDs). BLASTx-based test of homology indicated that A1-derived transcripts shared a high similarity with Gossypium arboreum (A2). Functional annotation of the CDs using the UniProt categorized them into biological processes, cellular components, and molecular function, COG classification showed that a large number of CDs have significant homology in COG database (6215 CDs), and mapping of CDs with Kyoto Encyclopedia of Genes and Genomes (KEGG) database generated 200 pathways ultimately showing predominant engagement in the fiber development process. Transcription factors were predicted by comparison with Plant Transcription Factor Database, and their differential expression between stages exposed their important regulatory role in fiber development. Differential expression analysis based on reads per kilobase of transcript per million mapped reads (RPKM) value revealed activities of specific gene related to carbohydrate and lipid synthesis, carbon metabolism, energy metabolism, signal transduction, etc., at four stages of fiber development, and was validated by qPCR. Overall, this study will help as a valuable foundation for diploid cotton fiber improvement.

Identification of candidate genes for fiber length quantitative trait loci through RNA-Seq and linkage and physical mapping in cotton

BACKGROUND

Cotton (Gossypium spp.) fibers are single-celled elongated trichomes, the molecular aspects of genetic variation in fiber length (FL) among genotypes are currently unknown. In this study, two backcross inbred lines (BILs), i.e., NMGA-062 ("Long") and NMGA-105 ("Short") with 32.1 vs. 27.2 mm in FL, respectively, were chosen to perform RNA-Seq on developing fibers at 10 days post anthesis (DPA). The two BILs differed in 4 quantitative trait loci (QTL) for FL and were developed from backcrosses between G. hirsutum as the recurrent parent and G. barbadense.

RESULTS

In total, 51.7 and 54.3 million reads were obtained and assembled to 49,508 and 49,448 transcripts in the two genotypes, respectively. Of 1551 differentially expressed genes (DEGs) between the two BILs, 678 were up-regulated and 873 down-regulated in "Long"; and 703 SNPs were identified in 339 DEGs. Further physical mapping showed that 8 DEGs were co-localized with the 4 FL QTL identified in the BIL population containing the two BILs. Four SNP markers in 3 DEGs that showed significant correlations with FL were developed. Among the three candidate genes encoding for proline-rich protein, D-cysteine desulfhydrase, and thaumatin-like protein, a SNP of thaumatin-like protein gene showed consistent correlations with FL across all testing environments.

CONCLUSIONS

This study represents one of the first investigations of positional candidate gene approach of QTL in cotton in integrating transcriptome and SNP identification based on RNA-Seq with linkage and physical mapping of QTL and genes, which will facilitate eventual cloning and identification of genes responsible for FL QTL. The candidate genes may serve as the foundation for further in-depth studies of the molecular mechanism of natural variation in fiber elongation.

The cotton β-galactosyltransferase 1 (GalT1) that galactosylates arabinogalactan proteins participates in controlling fiber development

Arabinogalactan proteins (AGPs) are highly glycosylated proteins that play pivotal roles in diverse developmental processes in plants. Type-II AG glycans, mostly O-linked to the hydroxyproline residues of the protein backbone, account for up to 95% w/w of the AGP, but their functions are still largely unclear. Cotton fibers are extremely elongated single-cell trichomes on the seed epidermis; however, little is known of the molecular basis governing the regulation of fiber cell development. Here, we characterized the role of a CAZy glycosyltransferase 31 (GT31) family member, GhGalT1, in cotton fiber development. The fiber length of the transgenic cotton overexpressing GhGalT1 was shorter than that of the wild type, whereas in the GhGalT1-silenced lines there was a notable increase in fiber length compared with wild type. The carbohydrate moieties of AGPs were altered in fibers of GhGalT1 transgenic cotton. The galactose: arabinose ratio of AG glycans was higher in GhGalT1 overexpression fibers, but was lower in GhGalT1-silenced lines, compared with that in the wild type. Overexpression of GhGalT1 upregulates transcript levels of a broad range of cell wall-related genes, especially the fasciclin-like AGP (FLA) backbone genes. An enzyme activity assay demonstrated that GhGalT1 is a β-1,3-galactosyltransferase (β-1,3-GalT) involved in biosynthesis of the β-1,3-galactan backbone of the type-II AG glycans of AGPs. We also show that GhGalT1 can form homo- and heterodimers with other cotton GT31 family members to facilitate AG glycan assembly of AGPs. Thus, our data demonstrate that GhGalT1 influences cotton fiber development via controlling the glycosylation of AGPs, especially FLAs.

Functional characterization of a basic helix-loop-helix (bHLH) transcription factor GhDEL65 from cotton (Gossypium hirsutum)

Cotton fiber is proposed to share some similarity with the Arabidopsis thaliana leaf trichome, which is regulated by the MYB-bHLH-WD40 transcription complex. Although several MYB transcription factors and WD40 family proteins in cotton have been characterized, little is known about the role of bHLH family proteins in cotton. Here, we report that GhDEL65, a bHLH protein from cotton (Gossypium hirsutum), is a functional homologue of Arabidopsis GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) in regulating trichome development. Transcripts of GhDEL65 were detected in 0 ∼ 1 days post-anthesis (DPA) ovules and abundant in 3-DPA fibers, implying that GhDEL65 may act in early fiber development. Ectopic expression of GhDEL65 in Arabidopsis gl3 egl3 double mutant partly rescued the trichome development, and constitutive expression of GhDEL65 in wild-type plants led to increased trichome density on rosette leaves and stems, mainly by activating the transcription of two key positive regulators of trichome development, GLABRA1 (GL1) and GLABRA2 (GL2), and suppressed the expression of a R3 single-repeat MYB factor TRIPTYCHON (TRY). GhDEL65 could interact with cotton R2R3 MYB transcription factors GhMYB2 and GhMYB3, as well as the WD40 protein GhTTG3, suggesting that the MYB-bHLH-WD40 protein complex also exists in cotton fiber cell, though its function in cotton fiber development awaits further investigation.

Identification of Multiple Stress Responsive Genes by Sequencing a Normalized cDNA Library from Sea-Land Cotton (Gossypium barbadense L.)

Background

Plants often face multiple stresses including drought, extreme temperature, salinity, nutrition deficiency and biotic stresses during growth and development. All the stresses result in a series of physiological and metabolic reactions and then generate reversible inhibition of metabolism and growth and can cause seriously irreversible damage, even death. At each stage of cotton growth, environmental stress conditions pose devastating threats to plant growth and development, especially yield and quality. Due to the complex stress conditions and unclear molecular mechanisms of stress response, there is an urgent need to explore the mechanisms of cotton response against abiotic stresses.

Methodology and Principal Findings

A normalized cDNA library was constructed using Gossypium barbadense Hai-7124 treated with different stress conditions (heat, cold, salt, drought, potassium and phosphorus deficit and Verticillium dahliae infection). Random sequencing of this library generated 6,047 high-quality expressed sequence tags (ESTs). The ESTs were clustered and assembled into 3,135 uniESTs, composed of 2,497 contigs and 638 singletons. The blastx results demonstrated 2,746 unigenes showing significant similarity to known genes, 74 uniESTs displaying significant similarity to genes of predicted proteins, and 315 uniESTs remain uncharacterized. Functional classification unveiled the abundance of uniESTs in binding, catalytic activity, and structural molecule activity. Annotations of the uniESTs by the plant transcription factor database (PlantTFDB) and Plant Stress Protein Database (PSPDB) disclosed that transcription factors and stress-related genes were enriched in the current library. The expression of some transcription factors and specific stress-related genes were verified by RT-PCR under various stress conditions.

Conclusions/Significance

Annotation results showed that a huge number of genes respond to stress in our study, such as MYB-related, C2H2, FAR1, bHLH, bZIP, MADS, and mTERF. These results will improve our knowledge of stress tolerance in cotton. In addition, they are also helpful in discovering candidate genes related to stress tolerance. The publicly available ESTs from G. barbadense are a valuable genomic resource that will facilitate further molecular study and breeding of stress-tolerant cotton.

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.

Global analysis of the Gossypium hirsutum L. Transcriptome during leaf senescence by RNA-Seq

BACKGROUND

Leaf senescence is an important developmental programmed degeneration process that dramatically affects crop quality and yield. The regulation of senescence is highly complex. Although senescence regulatory genes have been well characterized in model species such as Arabidopsis and rice, there is little information on the control of this process in cotton. Here, the senescence process in cotton (Gossypium hirsutum L.) leaves was investigated over a time course including young leaf, mature leaf and leaf samples from different senescence stages using RNA-Seq.

RESULTS

Of 24,846 genes detected by mapping the tags to Gossypium genomes, 3,624 genes were identified as differentially expressed during leaf senescence. There was some overlap between the genes identified here and senescence-associated genes previously identified in other species. Most of the genes related to photosynthesis, chlorophyll metabolism and carbon fixation were downregulated; whereas those for plant hormone signal transduction were upregulated. Quantitative real-time PCR was used to evaluate the results of RNA-Seq for gene expression profiles. Furthermore, 519 differentially expressed transcription factors were identified, notably WRKY, bHLH and C3H. In addition, 960 genes involved in the metabolism and regulation of eight hormones were identified, of which many genes involved in the abscisic acid, brassinosteroid, jasmonic acid, salicylic acid and ethylene pathways were upregulated, indicating that these hormone-related genes might play crucial roles in cotton leaf development and senescence. However, most auxin, cytokinin and gibberellin pathway-related genes were downregulated, suggesting that these three hormones may act as negative regulators of senescence.

CONCLUSIONS

This is the first high-resolution, multiple time-course, genome-wide comprehensive analysis of gene expression in cotton. These data are the most comprehensive dataset currently available for cotton leaf senescence, and will serve as a useful resource for unraveling the functions of many specific genes involved in cotton leaf development and senescence.

Two cotton fiber-associated glycosyltransferases, GhGT43A1 and GhGT43C1, function in hemicellulose glucuronoxylan biosynthesis during plant development

Xylan is the major hemicellulosic constituent in dicot secondary cell walls. Cell wall composition of cotton fiber changes dynamically throughout development. Not only the amounts but also the molecular sizes of the hemicellulosic polysaccharides show substantial changes during cotton fiber development. However, none of the genes encoding glycosyltransferases (GTs) responsible for synthesizing xylan have been isolated and characterized in cotton fiber. In this study, we applied a bioinformatics approach and identified two putative GTs from cotton, designated GhGT43A1 and GhGT43C1, which belong to the CAZy GT43 family and are closely related to Arabidopsis IRX9 and IRX14, respectively. We show that GhGT43A1 is highly and preferentially expressed in 15 and 20 days post-anthesis (dpa) cotton fiber, whereas GhGT43C1 is ubiquitously expressed in most organs, with especially high expression in 15 dpa fiber and hypocotyl. Complementation analysis demonstrates that GhG43A1 and GhGT43C1 are orthologs of Arabidopsis IRX9 and IRX14, respectively. Furthermore, we show that overexpression of GhGT43A1 or GhGT43C1 in Arabidopsis results in increased xylan content. We also show that overexpression of GhGT43A1 or GhGT43C1 leads to more cellulose deposition. These findings suggest that GhGT43A1 and GhGT43C1 likely participate in xylan synthesis during fiber development.

Elucidation of Nuclear and Organellar Genomes of Gossypium hirsutum: Furthering Studies of Species Evolution and Applications for Crop Improvement

Plant genomes are larger and more complex than other eukaryotic organisms, due to small and large duplication events, recombination and subsequent reorganization of the genetic material. Commercially important cotton is the result of a polyploidization event between Old and New World cottons that occurred over one million years ago. Allotetraploid cotton has properties that are dramatically different from its progenitors-most notably, the presence of long, spinnable fibers. Recently, the complete genome of a New World cotton ancestral species, Gossypium raimondii, was completed. Future genome sequencing efforts are focusing on an Old World progenitor, G. arboreum. This sequence information will enable us to gain insights into the evolution of the cotton genome that may be used to understand the evolution of other plant species. The chloroplast genomes of multiple cotton species and races have been determined. This information has also been used to gain insight into the evolutionary history of cotton. Analysis of the database of nuclear and organellar sequences will facilitate the identification of potential genes of interest and subsequent development of strategies for improving cotton.

Transcriptome profiling of Gossypium barbadense inoculated with Verticillium dahliae provides a resource for cotton improvement

BACKGROUND

Verticillium wilt, caused by the fungal pathogen Verticillium dahliae, is the most severe disease in cotton (Gossypium spp.), causing great lint losses worldwide. Disease management could be achieved in the field if genetically improved, resistant plants were used. However, the interaction between V. dahliae and cotton is a complicated process, and its molecular mechanism remains obscure. To understand better the defense response to this pathogen as a means for obtaining more tolerant cultivars, we monitored the transcriptome profiles of roots from resistant plants of G. barbadense cv. Pima90-53 that were challenged with V. dahliae.

RESULTS

In all, 46,192 high-quality expressed sequence tags (ESTs) were generated from a full-length cDNA library of G. barbadense. They were clustered and assembled into 23126 unigenes that comprised 2661 contigs and 20465 singletons. Those unigenes were assigned Gene Ontology terms and mapped to 289 KEGG pathways. A total of 3027 unigenes were found to be homologous to known defense-related genes in other plants. They were assigned to the functional classification of plant-pathogen interactions, including disease defenses and signal transduction. The branch of "SA→NPR1→TGA→PR-1→Disease resistance" was first discovered in the interaction of cotton-V. dahliae, indicating that this wilt process includes both biotrophic and necrotrophic stages. In all, 4936 genes coding for putative transcription factors (TF) were identified in our library. The most abundant TF family was the NAC group (527), followed by G2-like (440), MYB (372), BHLH (331), bZIP (271) ERF, C3H, and WRKY. We also analyzed the expression of genes involved in pathogen-associated molecular pattern (PAMP) recognition, the activation of effector-triggered immunity, TFs, and hormone biosynthesis, as well as genes that are pathogenesis-related, or have roles in signaling/regulatory functions and cell wall modification. Their differential expression patterns were compared among mock-/inoculated- and resistant/susceptible cotton. Our results suggest that the cotton defense response has significant transcriptional complexity and that large accumulations of defense-related transcripts may contribute to V. dahliae resistance in cotton. Therefore, these data provide a resource for cotton improvement through molecular breeding approaches.

CONCLUSIONS

This study generated a substantial amount of cotton transcript sequences that are related to defense responses against V. dahliae. These genomics resources and knowledge of important related genes contribute to our understanding of host-pathogen interactions and the defense mechanisms utilized by G. barbadense, a non-model plant system. These tools can be applied in establishing a modern breeding program that uses marker-assisted selections and oligonucleotide arrays to identify candidate genes that can be linked to valuable agronomic traits in cotton, including disease resistance.

Cotton PRP5 gene encoding a proline-rich protein is involved in fiber development

Proline-rich proteins contribute to cell wall structure of specific cell types and are involved in plant growth and development. In this study, a fiber-specific gene, GhPRP5, encoding a proline-rich protein was functionally characterized in cotton. GhPRP5 promoter directed GUS expression only in trichomes of both transgenic Arabidopsis and tobacco plants. The transgenic Arabidopsis plants with overexpressing GhPRP5 displayed reduced cell growth, resulting in smaller cell size and consequently plant dwarfs, in comparison with wild type plants. In contrast, knock-down of GhPRP5 expression by RNA interference in cotton enhanced fiber development. The fiber length of transgenic cotton plants was longer than that of wild type. In addition, some genes involved in fiber elongation and wall biosynthesis of cotton were up-regulated or down-regulated in the transgenic cotton plants owing to suppression of GhPRP5. Collectively, these data suggested that GhPRP5 protein as a negative regulator participates in modulating fiber development of cotton.

Cotton GalT1 encoding a putative glycosyltransferase is involved in regulation of cell wall pectin biosynthesis during plant development

Arabinogalactan proteins (AGPs), are a group of highly glycosylated proteins that are found throughout the plant kingdom. To date, glycosyltransferases that glycosylate AGP backbone have remained largely unknown. In this study, a gene (GhGalT1) encoding a putative β-1,3-galactosyltransferase (GalT) was identified in cotton. GhGalT1, belonging to CAZy GT31 family, is the type II membrane protein that contains an N-terminal transmembrane domain and a C-terminal galactosyltransferase functional domain. A subcellular localization assay demonstrated that GhGalT1 was localized in the Golgi apparatus. RT-PCR analysis revealed that GhGalT1 was expressed at relatively high levels in hypocotyls, roots, fibers and ovules. Overexpression of GhGalT1 in Arabidopsis promoted plant growth and metabolism. The transgenic seedlings had much longer primary roots, higher chlorophyll content, higher photosynthetic efficiency, the increased biomass, and the enhanced tolerance to exogenous D-arabinose and D-galactose. In addition, gas chromatography (GC) analysis of monosaccharide composition of cell wall fractions showed that pectin was changed in the transgenic plants, compared with that of wild type. Three genes (GAUT8, GAUT9 and xgd1) involved in pectin biosynthesis were dramatically up-regulated in the transgenic lines. These data suggested that GhGalT1 may be involved in regulation of pectin biosynthesis required for plant development.

Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species

Background

Cotton is the leading fiber crop worldwide. Gossypium barbadense is an important species of cotton because of its extra-long staple fibers with superior luster and silkiness. However, a systematic analysis and utilization of cDNA sequences from G. barbadense fiber development remains understudied.

Results

A total of 21,079 high quality sequences were generated from two non-normalized cDNA libraries prepared by using a mixture of G. barbadense Hai7124 fibers and ovules. After assembly processing, a set of 8,653 unigenes were obtained. Of those, 7,786 were matched to known proteins and 7,316 were assigned to functional categories. The molecular functions of these unigenes were mostly related to binding and catalytic activity, and carbohydrate, amino acid, and energy metabolisms were major contributors among the subsets of metabolism. Sequences comparison between G. barbadense and G. hirsutum revealed that 8,245 unigenes from G. barbadense were detected the similarity with those released publicly in G. hirsutum, however, the remaining 408 sequences had no hits against G. hirsutum unigenes database. Furthermore, 13,275 putative ESTs InDels loci involved in the orthologous and/or homoeologous differences between/within G. barbadense and G. hirsutum were discovered by in silico analyses, and 2,160 InDel markers were developed by ESTs with more than five insertions or deletions. By gel electrophoresis combined with sequencing verification, 71.11% candidate InDel loci were reconfirmed orthologous and/or homoeologous loci polymorphisms using G. hirsutum acc TM-1 and G. barbadense cv Hai7124. Blastx result showed among 2,160 InDel loci, 81 with significant function similarity with known genes associated with secondary wall synthesis process, indicating the important roles in fiber quality in tetraploid cultivated cotton species.

Conclusion

Sequence comparisons and InDel markers development will lay the groundwork for promoting the identification of genes related to superior agronomic traits, genetic differentiation and comparative genomic studies between G. hirsutum and G. barbadense.

Isolation, characterization and mapping of genes differentially expressed during fibre development between Gossypium hirsutum and G. barbadense by cDNA-SRAP

Gossypium hirsutum and G. barbadense are two cultivated tetraploid cotton species with differences in fibre quality. The fibre of G. barbadense is longer, stronger and finer than that of G. hirsutum. To isolate genes expressed differently between the two species during fibre development, cDNA-SRAP (sequence-related amplified polymorphism) was applied. This technique was used to analyse genes at different stages of fibre development in G. hirsutum cv. Emian22 and G. barbadense acc. 3-79, the parents of our interspecific mapping population. A total of 4096 SRAP primer combinations were used to screen polymorphism between the DNA of the parents, and 275 highly polymorphic primers were picked out to analyse DNA and RNA from leaves and fibres at different developmental stages of the parents. A total of 168 DNA fragments were isolated from gels and sequenced: 54, 30, 38 and 41 from fibres of 5, 10, 15 and 20 days post-anthesis, respectively, and five from multi stages. To genetically map these sequences, 104 sequence-specific primers were developed and were used to screened polymorphism between the mapping parents. Finally, six markers were mapped on six chromosomes of our backbone interspecific genetic map. This work can give us a primary knowledge of differences in mechanism of fibre development between G. hirsutum and G. barbadense.

Characterization of head transcriptome and analysis of gene expression involved in caste differentiation and aggression in Odontotermes formosanus (Shiraki)

BACKGROUND

The subterranean termite Odontotermes formosanus (Shiraki) is a serious insect pest of trees and dams in China. To date, very little is known about genomic or transcriptomic data for caste differentiation and aggression in O. formosanus. Hence, studies on transcriptome and gene expression profiling are helpful to better understand molecular basis underlying caste differentiation and aggressive behavior in O. formosanus.

METHODOLOGY AND PRINCIPAL FINDINGS

Using the Illumina sequencing, we obtained more than 57 million sequencing reads derived from the heads of O. formosanus. These reads were assembled into 116,885 unique sequences (mean size  =  536 bp). Of the unigenes, 30,646 (26.22%) had significant similarity with proteins in the NCBI nonredundant protein database and Swiss-Prot database (E-value<10(-5)). Of these annotated unigenes, 10,409 and 9,009 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. In total, 19,611 (25.52%) unigenes were mapped onto 242 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). A total of 11,661 simple sequence repeats (SSRs) were predicted from the current transcriptome database. Moreover, we detected seven putative genes involved in caste differentiation and six putative genes involved in aggression. The qPCR analysis showed that there were significant differences in the expression levels of the three putative genes hexamerin 2, β-glycosidase and bicaudal D involved in caste differentiation and one putative gene Cyp6a20 involved in aggression among workers, soldiers and larvae of O. formosanus.

CONCLUSIONS

To our knowledge, this is the first study to characterize the complete head transcriptome of a higher fungus-cultivating termite using high-throughput sequencing. Our study has provided the comprehensive sequence resources available for elucidating molecular basis underlying caste differentiation and aggressive behavior in O. formosanus.

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.

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.

Bioinformatics opportunities for identification and study of medicinal plants

Plants have been used as a source of medicine since historic times and several commercially important drugs are of plant-based origin. The traditional approach towards discovery of plant-based drugs often times involves significant amount of time and expenditure. These labor-intensive approaches have struggled to keep pace with the rapid development of high-throughput technologies. In the era of high volume, high-throughput data generation across the biosciences, bioinformatics plays a crucial role. This has generally been the case in the context of drug designing and discovery. However, there has been limited attention to date to the potential application of bioinformatics approaches that can leverage plant-based knowledge. Here, we review bioinformatics studies that have contributed to medicinal plants research. In particular, we highlight areas in medicinal plant research where the application of bioinformatics methodologies may result in quicker and potentially cost-effective leads toward finding plant-based remedies.