Comprehensive identification and expression analysis of CRY gene family in Gossypium.

BACKGROUND: The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance. RESULTS: In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments. CONCLUSION: We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

Comparative Genomics and Functional Studies of Putative m6A Methyltransferase (METTL) Genes in Cotton.

N6-methyladenosine (m6A) RNA modification plays important regulatory roles in plant development and adapting to the environment, which requires methyltransferases to achieve the methylation process. However, there has been no research regarding m6A RNA methyltransferases in cotton. Here, a systematic analysis of the m6A methyltransferase (METTL) gene family was performed on twelve cotton species, resulting in six METTLs identified in five allotetraploid cottons, respectively, and three to four METTLs in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that METTL genes from cottons, Arabidopsis thaliana, and Homo sapiens could be classified into three clades (METTL3, METTL14, and METTL-like clades). Cis-element analysis predicated the possible functions of METTL genes in G. hirsutum. RNA-seq data revealed that GhMETTL14 (GH_A07G0817/GH_D07G0819) and GhMETTL3 (GH_A12G2586/GH_D12G2605) had high expressions in root, stem, leaf, torus, petal, stamen, pistil, and calycle tissues. GhMETTL14 also had the highest expression in 20 and 25 dpa fiber cells, implying a potential role at the cell wall thickening stage. Suppressing GhMETTL3 and GhMETTL14 by VIGS caused growth arrest and even death in G. hirsutum, along with decreased m6A abundance from the leaf tissues of VIGS plants. Overexpression of GhMETTL3 and GhMETTL14 produced distinct differentially expressed genes (DEGs) in A. thaliana, indicating their possible divergent functions after gene duplication. Overall, GhMETTLs play indispensable but divergent roles during the growth of cotton plants, which provides the basis for the systematic investigation of m6A in subsequent studies to improve the agronomic traits in cotton.

Characterization of cotton ARF factors and the role of GhARF2b in fiber development.

BACKGROUND: Cotton fiber is a model system for studying plant cell development. At present, the functions of many transcription factors in cotton fiber development have been elucidated, however, the roles of auxin response factor (ARF) genes in cotton fiber development need be further explored. RESULTS: Here, we identify auxin response factor (ARF) genes in three cotton species: the tetraploid upland cotton G. hirsutum, which has 73 ARF genes, and its putative extent parental diploids G. arboreum and G. raimondii, which have 36 and 35 ARFs, respectively. Ka and Ks analyses revealed that in G. hirsutum ARF genes have undergone asymmetric evolution in the two subgenomes. The cotton ARFs can be classified into four phylogenetic clades and are actively expressed in young tissues. We demonstrate that GhARF2b, a homolog of the Arabidopsis AtARF2, was preferentially expressed in developing ovules and fibers. Overexpression of GhARF2b by a fiber specific promoter inhibited fiber cell elongation but promoted initiation and, conversely, its downregulation by RNAi resulted in fewer but longer fiber. We show that GhARF2b directly interacts with GhHOX3 and represses the transcriptional activity of GhHOX3 on target genes. CONCLUSION: Our results uncover an important role of the ARF factor in modulating cotton fiber development at the early stage.

Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium.

BACKGROUND: Cotton (Gossypium spp.) is the most important world-wide fiber crop but salt stress limits cotton production in coastal and other areas. Growth regulation factors (GRFs) play regulatory roles in response to salt stress, but their roles have not been studied in cotton under salt stress. RESULTS: We identified 19 GRF genes in G. raimondii, 18 in G. arboreum, 34 in G. hirsutum and 45 in G. barbadense, respectively. These GRF genes were phylogenetically analyzed leading to the recognition of seven GRF clades. GRF genes from diploid cottons (G. raimondii and G. arboreum) were largely retained in allopolyploid cotton, with subsequent gene expansion in G. barbadense relative to G. hirsutum. Most G. hirsutum GRF (GhGRF) genes are preferentially expressed in young and growing tissues. To explore their possible role in salt stress, we used qRT-PCR to study expression responses to NaCl treatment, showing that five GhGRF genes were down-regulated in leaves. RNA-seq experiments showed that seven GhGRF genes exhibited decreased expression in leaves under NaCl treatment, three of which (GhGRF3, GhGRF4, and GhGRF16) were identified by both RNA-seq and qRT-PCR. We also identified six and three GRF genes that exhibit decreased expression under salt stress in G. arboreum and G. barbadense, respectively. Consistent with its lack of leaf withering or yellowing under the salt treatment conditions, G. arboreum had better salt tolerance than G. hirsutum and G. barbadense. Our results suggest that GRF genes are involved in salt stress responses in Gossypium. CONCLUSION: In summary, we identified candidate GRF genes that were involved in salt stress responses in cotton.

Core cis-element variation confers subgenome-biased expression of a transcription factor that functions in cotton fiber elongation.

Cotton cultivars have evolved to produce extensive, long, seed-born fibers important for the textile industry, but we know little about the molecular mechanism underlying spinnable fiber formation. Here, we report how PACLOBUTRAZOL RESISTANCE 1 (PRE1) in cotton, which encodes a basic helix-loop-helix (bHLH) transcription factor, is a target gene of spinnable fiber evolution. Differential expression of homoeologous genes in polyploids is thought to be important to plant adaptation and novel phenotypes. PRE1 expression is specific to cotton fiber cells, upregulated during their rapid elongation stage and A-homoeologous biased in allotetraploid cultivars. Transgenic studies demonstrated that PRE1 is a positive regulator of fiber elongation. We determined that the natural variation of the canonical TATA-box, a regulatory element commonly found in many eukaryotic core promoters, is necessary for subgenome-biased PRE1 expression, representing a mechanism underlying the selection of homoeologous genes. Thus, variations in the promoter of the cell elongation regulator gene PRE1 have contributed to spinnable fiber formation in cotton. Overexpression of GhPRE1 in transgenic cotton yields longer fibers with improved quality parameters, indicating that this bHLH gene is useful for improving cotton fiber quality.

Targeting insect mitochondrial complex I for plant protection.

Plant engineered to express double-stranded RNA (dsRNA) can target the herbivorous insect gene for silencing. Although mounting evidence has emerged to support feasibility of this new pest control technology, field application is slow largely due to lack of potent targets. Here, we show that suppression of the gene encoding NDUFV2, a subunit of mitochondrial complex I that catalyses NADH dehydrogenation in respiratory chain, was highly lethal to insects. Feeding cotton bollworm (Helicoverpa armigera) larvae with transgenic cotton tissues expressing NDUFV2 dsRNA led to mortality up to 80% within 5 days, and almost no larvae survived after 7 days of feeding, due to the altered mitochondrial structure and activity. Transcriptome comparisons showed a drastic repression of dopa decarboxylase genes. Reciprocal assays with Asian corn borer (Ostrinia furnacalis), another lepidopteran species, revealed the sequence-specific effect of NDUFV2 suppression. Furthermore, the hemipteran lugus Apolygus lucorum was also liable to NDUFV2 repression. These data demonstrate that the mitochondrial complex I is a promising target with both sequence specificity and wide applicability for the development of new-generation insect-proof crops.

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.

GhHDZ76, a cotton HD-Zip transcription factor, involved in regulating the initiation and early elongation of cotton fiber development in G. hirsutum.

In this study, the whole HD-Zip family members of G. hirsutum were identified, and GhHDZ76 was classified into the HD-Zip IV subgroup. GhHDZ76 was predominantly expressed in the 0-5 DPA of fiber development stage and localized in the nucleus. Overexpression of GhHDZ76 significantly increased the length and density of trichomes in Arabidopsis thaliana. The fiber length of GhHDZ76 knockout lines by CRISPR/Cas9 was significantly shorter than WT at the early elongation and mature stage, indicating that GhHDZ76 positively regulate the fiber elongation. Scanning electron microscopy showed that the number of ovule surface protrusion of 0 DPA of GhHDZ76 knockout lines was significantly lower than WT, suggesting that GhHDZ76 can also promote the initiation of fiber development. The transcript level of GhWRKY16, GhRDL1, GhEXPA1 and GhMYB25 genes related to fiber initiation and elongation in GhHDZ76 knockout lines were significantly decreased. Yeast two-hybrid and Luciferase complementation imaging (LCI) assays showed that GhHDZ76 can interact with GhWRKY16 directly. As a transcription factor, GhHDZ76 has transcriptional activation activity, which could bind to L1-box elements of the promoters of GhRDL1 and GhEXPA1. Double luciferase reporter assay showed that the GhWRKY16 could enhance the transcriptional activity of GhHDZ76 to pGhRDL1, but it did not promote the transcriptional activity of GhHDZ76 to pGhEXPA1. GhHDZ76 protein may also promote the transcriptional activity of GhWRKY16 to the downstream target gene GhMYB25. Our results provided a new gene resource for fiber development and a theoretical basis for the genetic improvement of cotton fiber quality.

TRY intron2 determined its expression in inflorescence activated by SPL9 and MADS-box genes in Arabidopsis.

Plant trichomes are specialized epidermal cells that protect plants from insects and pathogens. In Arabidopsis, epidermal hairs decrease as internodes increase in height, with only few epidermal hairs produced on the sepals abaxial surface of the early flowers. TRIPTYCHON (TRY) is known to be a negative regulator of epidermal hair development in Arabidopsis, suppressing the formation of ectopic epidermal hairs in the inflorescence. Here, we reported that the second intron of TRY gene plays a critical role in trichome spatial distribution in Arabidopsis. The expression of TRY rises with the increasing stem nodes and reaches the peak in the inflorescence, while the trichomes distribution decrease. The transgenic plants showed that TRY promoter could only drive the genomic instead of coding sequences combined with GUS reporter gene, which indicates that the regulatory elements of TRY expression in inflorescence could be located in the intron regions. Multiple SPLs and MADS-box binding sites were found in the TRY intron2 sequence. Further genetic and biochemistry assays revealed that the flowering-related genes such as SPL9 could bind to these cis-elements directly, contributing to the TRY spatial expression. Since cotton fiber and Arabidopsis trichomes share similar regulatory mechanism, extended analysis showed that the intron2 of cotton TRY genes also contain the cis-elements. Thus, the introns harboring the transcription element may be the general way to regulate the gene expression in different plants and provides molecular clues for the related crops' traits design.

RES transformation for biosynthesis and detoxification.

The reactive electrophilic species (RES), typically the molecules bearing α,ß-unsaturated carbonyl group, are widespread in living organisms and notoriously known for their damaging effects. Many of the mycotoxins released from phytopathogenic fungi are RES and their contamination to cereals threatens food safety worldwide. However, due to their high reactivity, RES are also used by host organisms to synthesize specific metabolites. The evolutionary conserved glyoxalase (GLX) system scavenges the cytotoxic α-oxoaldehydes that bear RES groups, which cause host disorders and diseases. In cotton, a specialized enzyme derived from glyoxalase I (GLXI) through gene duplications and named as specialized GLXI (SPG), acts as a distinct type of aromatase in the gossypol pathway to transform the RES intermediates into the phenolic products. In this review, we briefly introduce the research progress in understanding the RES, especially the RES-type mycotoxins, the GLX system and SPG, and discuss their application potential in detoxification and synthetic biology.

The miR319-Targeted GhTCP4 Promotes the Transition from Cell Elongation to Wall Thickening in Cotton Fiber.

Plant cell growth involves a complex interplay among cell-wall expansion, biosynthesis, and, in specific tissues, secondary cell wall (SCW) deposition, yet the coordination of these processes remains elusive. Cotton fiber cells are developmentally synchronous, highly elongated, and contain nearly pure cellulose when mature. Here, we report that the transcription factor GhTCP4 plays an important role in balancing cotton fiber cell elongation and wall synthesis. During fiber development the expression of miR319 declines while GhTCP4 transcript levels increase, with high levels of the latter promoting SCW deposition. GhTCP4 interacts with a homeobox-containing factor, GhHOX3, and repressing its transcriptional activity. GhTCP4 and GhHOX3 function antagonistically to regulate cell elongation, thereby establishing temporal control of fiber cell transition to the SCW stage. We found that overexpression of GhTCP4A upregulated and accelerated activation of the SCW biosynthetic pathway in fiber cells, as revealed by transcriptome and promoter activity analyses, resulting in shorter fibers with varied lengths and thicker walls. In contrast, GhTCP4 downregulation led to slightly longer fibers and thinner cell walls. The GhHOX3-GhTCP4 complex may represent a general mechanism of cellular development in plants since both are conserved factors in many species, thus providing us a potential molecular tool for the design of fiber traits.

Promoter of a cotton fibre MYB gene functional in trichomes of Arabidopsis and glandular trichomes of tobacco.

Cotton fibres are unicellular seed trichomes. Our previous study suggested that the cotton R2R3 MYB transcript factor GaMYB2 is a functional homologue of the Arabidopsis trichome regulator GLABRA1 (GL1). Here, the GaMYB2 promoter activity is reported in cotton (Gossypium hirsutum), tobacco (Nicotiana tabacum), and Arabidopsis plants. A 2062 bp promoter of GaMYB2 was isolated from G. arboreum, and fused to a beta-glucuronidase (GUS) reporter gene. In cotton, the GaMYB2 promoter exhibited activities in developing fibre cells and trichomes of other aerial organs, including leaves, stems and bracts. In Arabidopsis the promoter was specific to trichomes. Different from Arabidopsis and cotton that have unicellular non-glandular simple trichomes, tobacco plants contain more than one type of trichome, including multicellular simple and glandular secreting trichomes (GSTs). Interestingly, in tobacco plants the GaMYB2 promoter directed GUS expression exclusively in glandular cells of GSTs. A series of 5'-deletions revealed that a 360 bp fragment upstream to the translation initiation codon was sufficient to drive gene expression. A putative cis-element of the T/G-box was located at -233 to -214; a yeast one-hybrid assay showed that Arabidopsis bHLH protein GLABRA3 (GL3), also a trichome regulator, and GhDEL65, a GL3-like cotton protein, had high binding activities to the T/G-box motif. Overexpression of GL3 or GhDEL65 enhanced the GaMYB2 promoter activity in transgenic Arabidopsis plants. A comparison of GaMYB2 promoter specificities in trichomes of different plant species with different types of trichomes provides a tool for further dissection of plant trichome structure and development.

Gossypium barbadense genome sequence provides insight into the evolution of extra-long staple fiber and specialized metabolites.

Of the two cultivated species of allopolyploid cotton, Gossypium barbadense produces extra-long fibers for the production of superior textiles. We sequenced its genome (AD)2 and performed a comparative analysis. We identified three bursts of retrotransposons from 20 million years ago (Mya) and a genome-wide uneven pseudogenization peak at 11-20 Mya, which likely contributed to genomic divergences. Among the 2,483 genes preferentially expressed in fiber, a cell elongation regulator, PRE1, is strikingly At biased and fiber specific, echoing the A-genome origin of spinnable fiber. The expansion of the PRE members implies a genetic factor that underlies fiber elongation. Mature cotton fiber consists of nearly pure cellulose. G. barbadense and G. hirsutum contain 29 and 30 cellulose synthase (CesA) genes, respectively; whereas most of these genes (>25) are expressed in fiber, genes for secondary cell wall biosynthesis exhibited a delayed and higher degree of up-regulation in G. barbadense compared with G. hirsutum, conferring an extended elongation stage and highly active secondary wall deposition during extra-long fiber development. The rapid diversification of sesquiterpene synthase genes in the gossypol pathway exemplifies the chemical diversity of lineage-specific secondary metabolites. The G. barbadense genome advances our understanding of allopolyploidy, which will help improve cotton fiber quality.

Control of cotton fibre elongation by a homeodomain transcription factor GhHOX3.

Cotton fibres are unusually long, single-celled epidermal seed trichomes and a model for plant cell growth, but little is known about the regulation of fibre cell elongation. Here we report that a homeodomain-leucine zipper (HD-ZIP) transcription factor, GhHOX3, controls cotton fibre elongation. GhHOX3 genes are localized to the 12th homoeologous chromosome set of allotetraploid cotton cultivars, associated with quantitative trait loci (QTLs) for fibre length. Silencing of GhHOX3 greatly reduces (>80%) fibre length, whereas its overexpression leads to longer fibre. Combined transcriptomic and biochemical analyses identify target genes of GhHOX3 that also contain the L1-box cis-element, including two cell wall loosening protein genes GhRDL1 and GhEXPA1. GhHOX3 interacts with GhHD1, another homeodomain protein, resulting in enhanced transcriptional activity, and with cotton DELLA, GhSLR1, repressor of the growth hormone gibberellin (GA). GhSLR1 interferes with the GhHOX3-GhHD1 interaction and represses target gene transcription. Our results uncover a novel mechanism whereby a homeodomain protein transduces GA signal to promote fibre cell elongation.

A cotton BURP domain protein interacts with α-expansin and their co-expression promotes plant growth and fruit production.

Plant growth requires cell wall extension. The cotton AtRD22-Like 1 gene GhRDL1, predominately expressed in elongating fiber cells, encodes a BURP domain-containing protein. Here, we show that GhRDL1 is localized in cell wall and interacts with GhEXPA1, an α-expansin functioning in wall loosening. Transgenic cotton overexpressing GhRDL1 showed an increase in fiber length and seed mass, and an enlargement of endopleura cells of ovules. Expression of either GhRDL1 or GhEXPA1 alone in Arabidopsis led to a substantial increase in seed size; interestingly, their co-expression resulted in the increased number of siliques, the nearly doubled seed mass, and the enhanced biomass production. Cotton plants overexpressing GhRDL1 and GhEXPA1 proteins produced strikingly more fruits (bolls), leading to up to 40% higher fiber yield per plant without adverse effects on fiber quality and vegetative growth. We demonstrate that engineering cell wall protein partners has a great potential in promoting plant growth and crop yield.

Down-regulation of S-adenosyl-L: -homocysteine hydrolase reveals a role of cytokinin in promoting transmethylation reactions.

S-adenosyl-L: -homocysteine hydrolase (SAHH) is a key enzyme for maintenance of cellular transmethylation potential. Although a cytokinin-binding activity had been hypothesized for SAHH, the relation between cytokinin and transmethylation reactions has not been elucidated. Here we show that, of the two Arabidopsis thaliana SAHH genes, AtSAHH1 has a much higher expression level than AtSAHH2. A T-DNA insertion mutant of AtSAHH1 (sahh1-1) and the RNA interference (RNAi) plants (dsAtSAHH2) accumulated a higher level of cytokinins, exhibited phenotypic changes similar to those of cytokinin-overproducers, and their global DNA methylation status was reduced. On the other hand, cytokinins positively regulate the transmethylation pathway genes, including AtSAHH1, AtADK1 (for adenosine kinase), and this regulation involves the cytokinin activity. Furthermore, expression of three cytosine DNA methyltransferase genes examined was inducible by cytokinin treatment. Unlike adenine and adenosine which are SAHH inhibitors, the adenine-type cytokinins have no effect on SAHH activity at protein level. Changing of endogenous cytokinin levels by transgene expression resulted in alterations of DNA methylation status in the sahh1-1 background, suggesting that cytokinins promote DNA methylation, at least under transmethylation stringent conditions. These data demonstrate that the phytohormone cytokinin plays a role in promoting transmethylation reactions, including DNA methylation.