Generation and analysis of a large-scale expressed sequence Tag database from a full-length enriched cDNA library of developing leaves of Gossypium hirsutum L

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.)

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole-genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome-wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.

LSD 3.0: a comprehensive resource for the leaf senescence research community

The leaf senescence database (LSD) is a comprehensive resource of senescence-associated genes (SAGs) and their corresponding mutants. Through manual curation and extensive annotation, we updated the LSD to a new version LSD 3.0, which contains 5853 genes and 617 mutants from 68 species. To provide sustainable and reliable services for the plant research community, LSD 3.0 (https://bigd.big.ac.cn/lsd/) has been moved to and maintained by the National Genomics Data Center at Beijing Institute of Genomics, Chinese Academy of Sciences. In the current release, we added some new features: (i) Transcriptome data of leaf senescence in poplar were integrated; (ii) Leaf senescence-associated transcriptome data information in Arabidopsis, rice and soybean were included; (iii) Senescence-differentially expressed small RNAs (Sen-smRNA) in Arabidopsis were identified; (iv) Interaction pairs between Sen-smRNAs and senescence-associated transcription factors (Sen-TF) were established; (v) Senescence phenotypes of 90 natural accessions (ecotypes) and 42 images of ecotypes in Arabidopsis were incorporated; (vi) Mutant seed information of SAGs in rice obtained from Kitbase was integrated; (vii) New options of search engines for ecotypes and transcriptome data were implemented. Together, the updated database bears great utility to continue to provide users with useful resources for studies of leaf senescence.

Establishing a System for Functional Characterization of Full-Length cDNAs of Camellia sinensis

Tea (Camellia sinensis) is enriched with bioactive secondary metabolites, and is one of the most popular nonalcoholic beverages globally. Two tea reference genomes have been reported; however, the functional analysis of tea genes has lagged, mainly due to tea’s recalcitrance to genetic transformation and the absence of alternative high throughput heterologous expression systems. A full-length cDNA collection with a streamlined cloning system is needed in this economically important woody crop species. RNAs were isolated from nine different vegetative tea tissues, pooled, then used to construct a normalized full-length cDNA library. The titer of unamplified and amplified cDNA library was 6.89 × 10⁶ and 1.8 × 10¹⁰ cfu/mL, respectively; the library recombinant rate was 87.2%. Preliminary characterization demonstrated that this collection can complement existing tea reference genomes and facilitate rare gene discovery. In addition, to streamline tea cDNA cloning and functional analysis, a binary vector (pBIG2113SF) was reengineered, seven tea cDNAs isolated from this library were successfully cloned into this vector, then transformed into Arabidopsis. One FL-cDNA, which encodes a putative P_1B-type ATPase 5 (CsHMA5), was characterized further as a proof of concept. We demonstrated that overexpression of CsHMA5 in Arabidopsis resulted in copper hyposensitivity. Thus, our data demonstrated that this represents an efficient system for rare gene discovery and functional characterization of tea genes. The integration of a tea FL-cDNA collection with efficient cloning and a heterologous expression system would facilitate functional annotation and characterization of tea genes.

Global analysis of DNA methylation in young (J1) and senescent (J2) Gossypium hirsutum L. cotyledons by MeDIP-Seq

DNA methylation is an important epigenetic modification regulating gene expression, genomic imprinting, transposon silencing and chromatin structure in plants and plays an important role in leaf senescence. However, the DNA methylation pattern during Gossypium hirsutum L. cotyledon senescence is poorly understood. In this study, global DNA methylation patterns were compared between two cotyledon development stages, young (J1) and senescence (J2), using methylated DNA immunoprecipitation (MeDIP-Seq). Methylated cytosine occurred mostly in repeat elements, especially LTR/Gypsy in both J1 and J2. When comparing J1 against J2, there were 1222 down-methylated genes and 623 up-methylated genes. Methylated genes were significantly enriched in carbohydrate metabolism, biosynthesis of other secondary metabolites and amino acid metabolism pathways. The global DNA methylation level decreased from J1 to J2, especially in gene promoters, transcriptional termination regions and regions around CpG islands. We further investigated the expression patterns of 9 DNA methyltransferase-associated genes and 2 DNA demethyltransferase-associated genes from young to senescent cotyledons, which were down-regulated during cotyledon development. In this paper, we first reported that senescent cotton cotyledons exhibited lower DNA methylation levels, primarily due to decreased DNA methyltransferase activity and which also play important role in regulating secondary metabolite process.

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.

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.

Genomic organization, differential expression, and functional analysis of the SPL gene family in Gossypium hirsutum

SQUAMOSA promoter binding protein-like (SPL) genes encode plant-specific transcription factors that are involved in many fundamental developmental processes. Certain SPL genes contain sequences complementary to miR156, a microRNA (miRNA) that plays a role in modulating plant gene expression. In this study, 30 SPL genes were identified in the reference genome of Gossypium raimondii and 24 GhSPLs were cloned from Gossypium hirsutum. G. raimondii is regarded as the putative contributor of the D-subgenome of G. hirsutum. Comparative analysis demonstrated sequence conservation between GhSPLs and other plant species. GhSPL genes could be classified into seven subclades based on phylogenetic analysis, diverse intron-exon structure, and motif prediction. Within each subclade, genes shared a similar structure. Sequence and experimental analysis predicted that 18 GhSPL genes are putative targets of GhmiR156. Additionally, tissue-specific expression analysis of GhSPL genes showed that their spatiotemporal expression patterns during development progressed differently, with most genes having high transcript levels in leaves, stems, and flowers. Finally, overexpression of GhSPL3 and GhSPL18 in Arabidopsis plants demonstrated that these two genes are involved in the development of leaves and second shoots and play an integral role in promoting flowering. The flowering integrator GhSOC1 may bind to the promoter of GhSPL3 but not GhSPL18 to regulate flowering. In conclusion, our analysis of GhSPL genes will provide some gene resources and a further understanding of GhSPL3 and GhSPL18 function in flowering promotion. Furthermore, the comparative genomics and functional analysis deepened our understanding of GhSPL genes during upland cotton vegetative and reproductive growth.