Lili Tu

Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry

The incompatible pathosystem between resistant cotton (Gossypium barbadense cv. 7124) and Verticillium dahliae strain V991 was used to study the cotton transcriptome changes after pathogen inoculation by RNA-Seq. Of 32 774 genes detected by mapping the tags to assembly cotton contigs, 3442 defence-responsive genes were identified. Gene cluster analyses and functional assignments of differentially expressed genes indicated a significant transcriptional complexity. Quantitative real-time PCR (qPCR) was performed on selected genes with different expression levels and functional assignments to demonstrate the utility of RNA-Seq for gene expression profiles during the cotton defence response. Detailed elucidation of responses of leucine-rich repeat receptor-like kinases (LRR-RLKs), phytohormone signalling-related genes, and transcription factors described the interplay of signals that allowed the plant to fine-tune defence responses. On the basis of global gene regulation of phenylpropanoid metabolism-related genes, phenylpropanoid metabolism was deduced to be involved in the cotton defence response. A closer look at the expression of these genes, enzyme activity, and lignin levels revealed differences between resistant and susceptible cotton plants. Both types of plants showed an increased level of expression of lignin synthesis-related genes and increased phenylalanine-ammonia lyase (PAL) and peroxidase (POD) enzyme activity after inoculation with V. dahliae, but the increase was greater and faster in the resistant line. Histochemical analysis of lignin revealed that the resistant cotton not only retains its vascular structure, but also accumulates high levels of lignin. Furthermore, quantitative analysis demonstrated increased lignification and cross-linking of lignin in resistant cotton stems. Overall, a critical role for lignin was believed to contribute to the resistance of cotton to disease.

Isolation and Characterization of Genes Associated to Cotton Somatic Embryogenesis by Suppression Subtractive Hybridization and Macroarray

Somatic embryogenesis (SE) is the developmental reprogramming of somatic cells toward the embryogenesis pathway and is a notable illustration of cell totipotency. To identify genes involved in SE, subtractive polymerase chain reaction (PCR) was performed to generate transcripts highly enriched for SE-related genes, using cDNA prepared from a mixture of embryogenic callus and preglobular somatic embryos, as the tester, and cDNA from nonembryogenic callus, as the driver. After differential screening and subsequent confirmation by reverse Northern blot analysis, a total of 671 differentially expressed cDNA fragments were identified, and 242 unigenes significantly up-regulated during cotton SE were recovered, as confirmed by Northern blot and reverse-transcription PCR analysis of representative cases, including most previously published SE-related genes in plants. In total, more than half had not been identified previously as SE-related genes, including dominant crucial genes involved in transcription, posttranscription, and transportation, and about one-third had not been reported previously to GenBank or were expected to be unknown, or newly identified genes. We used cDNA arrays to further investigate the expression patterns of these genes in differentiating gradient culture, ranging from proembryogenic masses to somatic embryos at every stage. The cDNA collection is composed of a broad repertoire of SE genes which is an important resource for understanding the genetic interactions underlying SE signaling and regulation. Our results suggested that a complicated and concerted mechanism involving multiple cellular pathways is responsible for cotton SE. This report represents a systematic and comprehensive analysis of genes involved in the process of somatic embryogenesis.

The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres

Gossypium hirsutum contributes the most production of cotton fibre, but G. barbadense is valued for its better comprehensive resistance and superior fibre properties. However, the allotetraploid genome of G. barbadense has not been comprehensively analysed. Here we present a high-quality assembly of the 2.57 gigabase genome of G. barbadense, including 80,876 protein-coding genes. The double-sized genome of the A (or At) (1.50 Gb) against D (or Dt) (853 Mb) primarily resulted from the expansion of Gypsy elements, including Peabody and Retrosat2 subclades in the Del clade, and the Athila subclade in the Athila/Tat clade. Substantial gene expansion and contraction were observed and rich homoeologous gene pairs with biased expression patterns were identified, suggesting abundant gene sub-functionalization occurred by allopolyploidization. More specifically, the CesA gene family has adapted differentially temporal expression patterns, suggesting an integrated regulatory mechanism of CesA genes from At and Dt subgenomes for the primary and secondary cellulose biosynthesis of cotton fibre in a “relay race”-like fashion. We anticipate that the G. barbadense genome sequence will advance our understanding the mechanism of genome polyploidization and underpin genome-wide comparison research in this genus.

A thaumatin-like protein gene involved in cotton fiber secondary cell wall development enhances resistance against Verticillium dahliae and other stresses in transgenic tobacco.

For the first time, a sea-island cotton (Gossypium barbadense L.) thaumatin-like protein gene (GbTLP1) with a potential role in secondary cell wall development has been overexpressed in tobacco to elucidate its function. The presence of the transgene was verified by Southern blotting and higher expression levels of GbTLP1 in transgenic tobacco plants were revealed by reverse-transcription and quantitative real-time polymerase chain reaction analyses. Transgenic plants with constitutively higher expression of the GbTLP1 showed enhanced resistance against different stress agents, particularly, its performance against Verticillium dahliae was exceptional. Transgenic tobacco plants also exhibited considerable resistance against Fusarium oxysporum and some abiotic stresses including salinity and drought. In this experiment, transgenic plants without GbTLP1 expression were also used as controls, which behaved similar to non-transgenic control plants. Overexpression of GbTLP1 had no significant deleterious effect on plant growth except that flowering was delayed for 3-5 weeks. The apparent pleiotropic effect of this novel gene has given us insight to the plant defense mechanism.

Detection of somaclonal variation of cotton (Gossypium hirsutum) using cytogenetics, flow cytometry and molecular markers

Two protocols of plant regeneration for cotton were adopted in this study, namely, 2, 4-D and kinetin hormone combination and IBA and kinetin hormone combination. Twenty-eight embryogenic cell lines via somatic embryogenesis and 67 regenerated plants from these embryogenic calli were selected and used for random amplified polymorphic DNA (RAPD), simple sequence repeat (SSR), chromosomal number counting, and flow cytometric analysis. The roles of RAPD and SSR markers in detecting somaclonal variation of cotton (Gossypium hirsutum L.) were evaluated. Two cluster analyses were performed to express, in the form of dendrograms, the relationships among the hormone combinations and the genetic variability. Both DNA-based techniques were able to amplify all of the cell clones and regenerated plantlets genomes and relative higher genetic variation could be detected in the culture type with 2, 4-D and kinetin hormone combination. The result suggested that 2, 4-D and kinetin hormone combination could induce relative high somaclonal variation and RAPD and SSR markers are useful in detecting somaclonal variation of regenerated cotton plants via somatic embryogenesis. Chromosome number counting and flow cytometry analysis revealed that the number of chromosomes and ploidy levels were nearly stable in all regenerated plants except two regenerated plantlets (lost 4 and 5 chromosomes, respectively) which meant that cytological changes were not correlated with the frequency of RAPD and SSR polymorphisms. This result also might mean that the cell lines with variation of chromosome numbers were difficult to regenerate plants.

GbTCP, a cotton TCP transcription factor, confers fibre elongation and root hair development by a complex regulating system

As the most important natural raw material for textile industry, cotton fibres are an excellent model for studying single-cell development. Although expression profiling and functional genomics have provided some data, the mechanism of fibre development is still not well known. A class I TCP transcription factor (designated GbTCP), encoding 344 amino acids, was isolated from the normalized cDNA library of sea-island cotton fibre (from –2 to 25 days post anthesis). GbTCP was preferentially expressed in the elongating cotton fibre from 5 to 15 days post anthesis. Some expression was also observed in stems, apical buds, and petals. RNAi silencing of GbTCP produced shorter fibre, a reduced lint percentage, and a lower fibre quality than the wild-type plants. Overexpression of GbTCP enhanced root hair initiation and elongation in Arabidopsis and regulated branching. Solexa sequencing and Affymetrix GeneChip analysis indicated that GbTCP positively regulates the level of jasmonic acid (JA) and, as a result, activates downstream genes (reactive oxygen species, calcium signalling, ethylene biosynthesis and response, and several NAC and WRKY transcription factors) necessary for elongation of fibres and root hairs. JA content analysis in cotton also confirmed that GbTCP has a profound effect on JA biosynthesis. In vitro ovule culture showed that an appropriate concentration of JA promoted fibre elongation. The results suggest that GbTCP is an important transcription factor for fibre and root hair development by regulating JA biosynthesis and response and other pathways, including reactive oxygen species, calcium channel and ethylene signalling.

A Genetic and Metabolic Analysis Revealed that Cotton Fiber Cell Development Was Retarded by Flavonoid Naringenin

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

GbPDF1 Is Involved in Cotton Fiber Initiation via the Core cis-Element HDZIP2ATATHB2

Cotton (Gossypium spp.) fiber cells are seed trichomes derived from the epidermal layer of the cotton seed coat. The molecular components responsible for regulating fiber cell differentiation have not been fully elucidated. A cotton PROTODERMAL FACTOR1 gene (GbPDF1) was found to be expressed preferentially during fiber initiation and early elongation, with highest accumulation in fiber cells 5 d post anthesis. PDF1 silencing caused retardation of fiber initiation and produced shorter fibers and lower lint percentage compared with the wild type, indicating that the gene is required for cotton fiber development. Further analysis showed that a higher accumulation of hydrogen peroxide occurred in the RNA interference transgenic cotton lines. Meanwhile, the expression of several genes related to ethylene and pectin synthesis or sugar transport during cotton fiber growth was found to be significantly reduced in the PDF1-suppressed cotton. Three proteins interacting with GbPDF1 in yeast and in planta might involve cellular signaling or metabolism. GbPDF1 promoter::GUS constructs in transgenic cotton were predominantly expressed in the epidermis of ovules and developing fibers. Progressive deletions of the GbPDF1 promoter showed that a 236-bp promoter fragment was sufficient for basal GbPDF1 transcription in cotton. Mutation of putative regulatory sequences showed that HDZIP2ATATHB2, an element within the fragment, was essential for PGbPDF1-1 expression. The binding activity between this cis-element and nuclear extracts from fiber-bearing cotton ovules at 5 d post anthesis was specific. We conclude that GbPDF1 plays a critical role together with interaction partners in hydrogen peroxide homeostasis and steady biosynthesis of ethylene and pectin during fiber development via the core cis-element HDZIP2ATATHB2.

Long noncoding RNAs and their proposed functions in fibre development of cotton (Gossypium spp.)

Long noncoding RNAs (lncRNAs) are transcripts of at least 200 bp in length, possess no apparent coding capacity and are involved in various biological regulatory processes. Until now, no systematic identification of lncRNAs has been reported in cotton (Gossypium spp.). Here, we describe the identification of 30 550 long intergenic noncoding RNA (lincRNA) loci (50 566 transcripts) and 4718 long noncoding natural antisense transcript (lncNAT) loci (5826 transcripts). LncRNAs are rich in repetitive sequences and preferentially expressed in a tissue-specific manner. The detection of abundant genome-specific and/or lineage-specific lncRNAs indicated their weak evolutionary conservation. Approximately 76% of homoeologous lncRNAs exhibit biased expression patterns towards the At or Dt subgenomes. Compared with protein-coding genes, lncRNAs showed overall higher methylation levels and their expression was less affected by gene body methylation. Expression validation in different cotton accessions and coexpression network construction helped to identify several functional lncRNA candidates involved in cotton fibre initiation and elongation. Analysis of integrated expression from the subgenomes of lncRNAs generating miR397 and its targets as a result of genome polyploidization indicated their pivotal functions in regulating lignin metabolism in domesticated tetraploid cotton fibres. This study provides the first comprehensive identification of lncRNAs in Gossypium.

Asymmetric subgenome selection and cis -regulatory divergence during cotton domestication

Comparative population genomics offers an excellent opportunity for unraveling the genetic history of crop domestication. Upland cotton (Gossypium hirsutum) has long been an important economic crop, but a genome-wide and evolutionary understanding of the effects of human selection is lacking. Here, we describe a variation map for 352 wild and domesticated cotton accessions. We scanned 93 domestication sweeps occupying 74 Mb of the A subgenome and 104 Mb of the D subgenome, and identified 19 candidate loci for fiber-quality-related traits through a genome-wide association study. We provide evidence showing asymmetric subgenome domestication for directional selection of long fibers. Global analyses of DNase I–hypersensitive sites and 3D genome architecture, linking functional variants to gene transcription, demonstrate the effects of domestication on cis-regulatory divergence. This study provides new insights into the evolution of gene organization, regulation and adaptation in a major crop, and should serve as a rich resource for genome-based cotton improvement.