Transcriptome profiling, molecular biological, and physiological studies reveal a major role for ethylene in cotton fiber cell elongation

AKR2A participates in the regulation of cotton fibre development by modulating biosynthesis of very-long-chain fatty acids

The biosynthesis of very-long-chain fatty acids (VLCFAs) and their transport are required for fibre development. However, whether other regulatory factors are involved in this process is unknown. We report here that overexpression of an Arabidopsis gene ankyrin repeat-containing protein 2A (AKR2A) in cotton promotes fibre elongation. RNA-Seq analysis was employed to elucidate the mechanisms of AKR2A in regulating cotton fibre development. The VLCFA content and the ratio of VLCFAs to short-chain fatty acids increased in AKR2A transgenic lines. In addition, AKR2A promotes fibre elongation by regulating ethylene and synergizing with the accumulation of auxin and hydrogen peroxide. Analysis of RNA-Seq data indicates that AKR2A up-regulates transcript levels of genes involved in VLCFAs' biosynthesis, ethylene biosynthesis, auxin and hydrogen peroxide signalling, cell wall and cytoskeletal organization. Furthermore, AKR2A interacted with KCS1 in Arabidopsis both in vitro and in vivo. Moreover, the VLCFA content and the ratio of VLCFAs to short-chain fatty acids increased significantly in seeds of AKR2A-overexpressing lines and AKR2A/KCS1 co-overexpressing lines, while AKR2A mutants are the opposite trend. Our results uncover a novel cotton fibre growth mechanism by which the critical regulator AKR2A promotes fibre development via activating hormone signalling cascade by mediating VLCFA biosynthesis. This study provides a potential candidate gene for improving fibre yield and quality through genetic engineering.

Insight into the relationship between S-lignin and fiber quality based on multiple research methods

Cotton (Gossypium hirsutum) is an important cash crop, providing people with high quality natural fiber. Lignin is the main component of cotton fiber, second only to cellulose. As a main substance filled in the cellulose framework during the secondary wall thickening process, lignin plays a key role in the formation of cotton fiber quality. However, the mechanism behind it is still unclear. In this research, we screened candidate genes involved in lignin biosynthesis based on analysis of cotton genome and transcriptome sequence data. The authenticity of the transcriptome data was verified by qRT-PCR assay. Total 62 genes were identified from nine gene families. In the process, we found the key gene GhCAD7 that affects the biosynthesis of S-lignin and the ratio of syringyl/guaiacyl (S/G). In addition, in combination with the metabolites and transcriptome profiles of the line 0-153 with high fiber quality and the line sGK9708 with low fiber quality during cotton fiber development, we speculate that the ratio of syringyl/guaiacyl (S/G) is inseparable from the quality of cotton fiber. Finally, the S-type lignin synthesis branch may play a more important role in the formation of high-quality fiber. This work provides insights into the synthesis of lignin in cotton and lays the foundation for future research into improving fiber quality.

Comprehensive analysis of the Gossypium hirsutum L. respiratory burst oxidase homolog (Ghrboh) gene family

BACKGROUND

Plant NADPH oxidase (NOX), also known as respiratory burst oxidase homolog (rboh), encoded by the rboh gene, is a key enzyme in the reactive oxygen species (ROS) metabolic network. It catalyzes the formation of the superoxide anion (O2•-), a type of ROS. In recent years, various studies had shown that members of the plant rboh gene family were involved in plant growth and developmental processes as well as in biotic and abiotic stress responses, but little is known about its functional role in upland cotton.

RESULTS

In the present study, 26 putative Ghrboh genes were identified and characterized. They were phylogenetically classified into six subfamilies and distributed at different densities across 18 of the 26 chromosomes or scaffolds. Their exon-intron structures, conserved domains, synteny and collinearity, gene family evolution, regulation mediated by cis-acting elements and microRNAs (miRNAs) were predicted and analyzed. Additionally, expression profiles of Ghrboh gene family were analyzed in different tissues/organs and at different developmental stages and under different abiotic stresses, using RNA-Seq data and real-time PCR. These profiling studies indicated that the Ghrboh genes exhibited temporal and spatial specificity with respect to expression, and might play important roles in cotton development and in stress tolerance through modulating NOX-dependent ROS induction and other signaling pathways.

CONCLUSIONS

This comprehensive analysis of the characteristics of the Ghrboh gene family determined features such as sequence, synteny and collinearity, phylogenetic and evolutionary relationship, expression patterns, and cis-element- and miRNA-mediated regulation of gene expression. Our results will provide valuable information to help with further gene cloning, evolutionary analysis, and biological function analysis of cotton rbohs.

Investigation of the EIL/EIN3 Transcription Factor Gene Family Members and Their Expression Levels in the Early Stage of Cotton Fiber Development

The ethylene-insensitive3-like/ethylene-insensitive3 (EIL/EIN3) protein family can serve as a crucial factor for plant growth and development under diverse environmental conditions. EIL/EIN3 protein is a form of a localized nuclear protein with DNA-binding activity that potentially contributes to the intricate network of primary and secondary metabolic pathways of plants. In light of recent research advances, next-generation sequencing (NGS) and novel bioinformatics tools have provided significant breakthroughs in the study of the EIL/EIN3 protein family in cotton. In turn, this paved the way to identifying and characterizing the EIL/EIN3 protein family. Hence, the high-throughput, rapid, and cost-effective meta sequence analyses have led to a remarkable understanding of protein families in addition to the discovery of novel genes, enzymes, metabolites, and other biomolecules of the higher plants. Therefore, this work highlights the recent advance in the genomic-sequencing analysis of higher plants, which has provided a plethora of function profiles of the EIL/EIN3 protein family. The regulatory role and crosstalk of different metabolic pathways, which are apparently affected by these transcription factor proteins in one way or another, are also discussed. The ethylene hormone plays an important role in the regulation of reactive oxygen species in plants under various environmental stress circumstances. EIL/EIN3 proteins are the key ethylene-signaling regulators and play important roles in promoting cotton fiber developmental stages. However, the function of EIL/EIN3 during initiation and early elongation stages of cotton fiber development has not yet been fully understood. The results provided valuable information on cotton EIL/EIN3 proteins, as well as a new vision into the evolutionary relationships of this gene family in cotton species.

Comparative transcriptomics of stem bark reveals genes associated with bast fiber development in Boehmeria nivea L. gaud (ramie)

BACKGROUND

Boehmeria nivea L. Gaud (Ramie) produces one of the longest natural fibers in nature. The bark of ramie mainly comprises of the phloem tissue of stem and is the raw material for fiber. Therefore, identifying the molecular regulation of phloem development is important for understanding of bast fiber biosynthesis and improvement of fiber quality in ramie.

RESULTS

In this study, we collected top bud (TB), bark from internode elongating region (ER) and bark from internode fully elongated region (FER) from the ramie variety Zhongzhu No. 1. Histological study indicated that these samples contain phloem tissues at different developmental and maturation stages, with a higher degree of maturation of phloem tissue in FER. RNA sequencing (RNA-seq) was performed and de novo transcriptome was assembled. Unigenes and differentially expressed genes (DEGs) in these three samples were identified. The analysis of DEGs by using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed clear differences in gene expression between ER and FER. Some unigenes involved in secondary cell wall biosynthesis were up-regulated in both ER and FER, while unigenes for some cell wall components or cell wall modifications showed differential expression between ER and FER. In addition, the ethylene respond factors (ERFs) in the ethylene signaling pathway were up-regulated in FER, and ent-kaurenoic acid oxidase (KAO) and GA 20-oxidase (GA20ox) for gibberellins biosynthesis were up-regulated while GA 2-oxidase (GA2ox) for gibberellin inactivation was down-regulated in FER.

CONCLUSIONS

Both morphological study and gene expression analysis supported a burst of phloem and vascular developmental processes during the fiber maturation in the ramie stem, and ethylene and gibberellin are likely to be involved in this process. Our findings provide novel insights into the phloem development and fiber maturation in ramie, which could be useful for fiber improvement in ramie and other fiber crops.

Histone Deacetylase (HDAC) Gene Family in Allotetraploid Cotton and Its Diploid Progenitors: In Silico Identification, Molecular Characterization, and Gene Expression Analysis under Multiple Abiotic Stresses, DNA Damage and Phytohormone Treatments

Histone deacetylases (HDACs) play a significant role in a plant’s development and response to various environmental stimuli by regulating the gene transcription. However, HDACs remain unidentified in cotton. In this study, a total of 29 HDACs were identified in allotetraploid Gossypium hirsutum, while 15 and 13 HDACs were identified in Gossypium arboretum and Gossypium raimondii, respectively. Gossypium HDACs were classified into three groups (reduced potassium dependency 3 (RPD3)/HDA1, HD2-like, and Sir2-like (SRT) based on their sequences, and Gossypium HDACs within each subgroup shared a similar gene structure, conserved catalytic domains and motifs. Further analysis revealed that Gossypium HDACs were under a strong purifying selection and were unevenly distributed on their chromosomes. Gene expression data revealed that G. hirsutumHDACs were differentially expressed in various vegetative and reproductive tissues, as well as at different developmental stages of cotton fiber. Furthermore, some G. hirsutum HDACs were co-localized with quantitative trait loci (QTLs) and single-nucleotide polymorphism (SNPs) of fiber-related traits, indicating their function in fiber-related traits. We also showed that G. hirsutum HDACs were differentially regulated in response to plant hormones (abscisic acid (ABA) and auxin), DNA damage agent (methyl methanesulfonate (MMS)), and abiotic stresses (cold, salt, heavy metals and drought), indicating the functional diversity and specification of HDACs in response to developmental and environmental cues. In brief, our results provide fundamental information regarding G.hirsutumHDACs and highlight their potential functions in cotton growth, fiber development and stress adaptations, which will be helpful for devising innovative strategies for the improvement of cotton fiber and stress tolerance.

Designed Manipulation of the Brassinosteroid Signal to Enhance Crop Yield

Brassinosteroid (BR), a plant steroid hormone, plays crucial role in modulating plant growth and development, which affect crop architecture and yield. However, BR application cannot highly benefit to agricultural production as expectation, because it regulates multiple processes in different tissues and leads to side effect. In addition, accurately modifying BR signal at transcriptional level is difficult. Effective manipulation of the BR signal and avoidance of side effects are required to enhance yield in different crops. Application of BR by spraying at specific developmental stages can enhance crop yield, but this method is impractical for use on a large scale. The accurate molecular design of crops would be much more helpful to manipulate the BR signal in specific organs and/or at particular developmental stages to enhance crop yield. This minireview summarizes the BR regulation of yield in different crops, especially horticultural crops, and the strategies used to regulate the BR signal to enhance crop yield. One popular strategy is to directly modulate the BR signal through modifying the functions of important components in the BR signal transduction pathway. Another strategy is to identify and modulate regulators downstream of, or in crosstalk with, the BR signal to manipulate its role in specific processes and increase crop yield. Efforts to accurately design a BR manipulation strategy will ultimately lead to effective control of the BR signal to avoid side effects and enhance crop yield.

Analysis of the MIR160 gene family and the role of MIR160a_A05 in regulating fiber length in cotton

The MIR160 family in Gossypium hirsutum and G. barbadense was characterized, and miR160a_A05 was found to increase cotton-fiber length by downregulating its target gene (ARF17) and several GH3 genes. Cotton fiber is the most important raw material for the textile industry. MicroRNAs are involved in regulating cotton-fiber development, but a role in fiber elongation has not been demonstrated. In this study, miR160a was found to be differentially expressed in elongating fibers between two interspecific (between Gossypium hirsutum and G. barbadense) backcross inbred lines (BILs) with different fiber lengths. The gene MIR160 colocalized with a previously mapped fiber-length quantitative trait locus. Its target gene ARF17 was differentially expressed between the two BILs during fiber elongation, but in the inverse fashion. Bioinformatics was used to analyze the MIR160 family in both G. hirsutum and G. barbadense. Moreover, qRT-PCR analysis identified MIR160a as the functional MIR160 gene encoding the miR160a precursor during fiber elongation. Using virus-induced gene silencing and overexpression, overexpressed MIR160a_A05 resulted in significantly longer fibers compared with wild type, whereas suppression of miR160 resulted in significantly shorter fibers. Expression levels of the target gene auxin-response factor 17 (ARF17) and related genes GH3 in the two BILs and/or the virus-infected plants demonstrated similar changes in response to modulation of miR160a level. Finally, overexpression or suppression of miR160 increased or decreased, respectively, the cellular level of indole-3-acetic acid, which is involved in fiber elongation. These results describe a specific regulatory mechanism for fiber elongation in cotton that can be utilized for future crop improvement.

High-density linkage map construction and QTL analyses for fiber quality, yield and morphological traits using CottonSNP63K array in upland cotton (Gossypium hirsutum L.)

Background

Improving fiber quality and yield are the primary research objectives in cotton breeding for enhancing the economic viability and sustainability of Upland cotton production. Identifying the quantitative trait loci (QTL) for fiber quality and yield traits using the high-density SNP-based genetic maps allows for bridging genomics with cotton breeding through marker assisted and genomic selection. In this study, a recombinant inbred line (RIL) population, derived from cross between two parental accessions, which represent broad allele diversity in Upland cotton, was used to construct high-density SNP-based linkage maps and to map the QTLs controlling important cotton traits.

Results

Molecular genetic mapping using RIL population produced a genetic map of 3129 SNPs, mapped at a density of 1.41 cM. Genetic maps of the individual chromosomes showed good collinearity with the sequence based physical map. A total of 106 QTLs were identified which included 59 QTLs for six fiber quality traits, 38 QTLs for four yield traits and 9 QTLs for two morphological traits. Sub-genome wide, 57 QTLs were mapped in A sub-genome and 49 were mapped in D sub-genome. More than 75% of the QTLs with favorable alleles were contributed by the parental accession NC05AZ06. Forty-six mapped QTLs each explained more than 10% of the phenotypic variation. Further, we identified 21 QTL clusters where 12 QTL clusters were mapped in the A sub-genome and 9 were mapped in the D sub-genome. Candidate gene analyses of the 11 stable QTL harboring genomic regions identified 19 putative genes which had functional role in cotton fiber development.

Conclusion

We constructed a high-density genetic map of SNPs in Upland cotton. Collinearity between genetic and physical maps indicated no major structural changes in the genetic mapping populations. Most traits showed high broad-sense heritability. One hundred and six QTLs were identified for the fiber quality, yield and morphological traits. Majority of the QTLs with favorable alleles were contributed by improved parental accession. More than 70% of the mapped QTLs shared the similar map position with previously reported QTLs which suggest the genetic relatedness of Upland cotton germplasm. Identification of QTL clusters could explain the correlation among some fiber quality traits in cotton. Stable and major QTLs and QTL clusters of traits identified in the current study could be the targets for map-based cloning and marker assisted selection (MAS) in cotton breeding. The genomic region on D12 containing the major stable QTLs for micronaire, fiber strength and lint percentage could be potential targets for MAS and gene cloning of fiber quality traits in cotton.

Extensin arabinoside chain length is modulated in elongating cotton fibre

Cotton fibre provides a unicellular model system for studying cell expansion and secondary cell wall deposition. Mature cotton fibres are mainly composed of cellulose while the walls of developing fibre cells contain a variety of polysaccharides and proteoglycans required for cell expansion. This includes hydroxyproline-rich glycoproteins (HRGPs) comprising the subgroup, extensins. In this study, extensin occurrence in cotton fibres was assessed using carbohydrate immunomicroarrays, mass spectrometry and monosaccharide profiling. Extensin amounts in three species appeared to correlate with fibre quality. Fibre cell expression profiling of the four cotton cultivars, combined with extensin arabinoside chain length measurements during fibre development, demonstrated that arabinoside side-chain length is modulated during development. Implications and mechanisms of extensin side-chain length dynamics during development are discussed.

Comparative Genome-wide Analysis and Expression Profiling of Histone Acetyltransferase (HAT) Gene Family in Response to Hormonal Applications, Metal and Abiotic Stresses in Cotton

Post-translational modifications are involved in regulating diverse developmental processes. Histone acetyltransferases (HATs) play vital roles in the regulation of chromation structure and activate the gene transcription implicated in various cellular processes. However, HATs in cotton, as well as their regulation in response to developmental and environmental cues, remain unidentified. In this study, 9 HATs were identified from Gossypium raimondi and Gossypium arboretum, while 18 HATs were identified from Gossypium hirsutum. Based on their amino acid sequences, Gossypium HATs were divided into three groups: CPB, GNAT, and TAF_II250. Almost all the HATs within each subgroup share similar gene structure and conserved motifs. Gossypium HATs are unevenly distributed on the chromosomes, and duplication analysis suggests that Gossypium HATs are under strong purifying selection. Gene expression analysis showed that Gossypium HATs were differentially expressed in various vegetative tissues and at different stages of fiber development. Furthermore, all the HATs were differentially regulated in response to various stresses (salt, drought, cold, heavy metal and DNA damage) and hormones (abscisic acid (ABA) and auxin (NAA)). Finally, co-localization of HAT genes with reported quantitative trait loci (QTL) of fiber development were reported. Altogether, these results highlight the functional diversification of HATs in cotton growth and fiber development, as well as in response to different environmental cues. This study enhances our understanding of function of histone acetylation in cotton growth, fiber development, and stress adaptation, which will eventually lead to the long-term improvement of stress tolerance and fiber quality in cotton.

Comparative Proteomic Analysis of Molecular Differences between Leaves of Wild-Type Upland Cotton and Its Fuzzless-Lintless Mutant

Fuzzless-lintless mutant (fl) ovules of upland cotton have been used to investigate cotton fiber development for decades. However, the molecular differences of green tissues between fl and wild-type (WT) cotton were barely reported. Here, we found that gossypol content, the most important secondary metabolite of cotton leaves, was higher in Gossypium hirsutum L. cv Xuzhou-142 (Xu142) WT than in fl. Then, we performed comparative proteomic analysis of the leaves from Xu142 WT and its fl. A total of 4506 proteins were identified, of which 103 and 164 appeared to be WT- and fl-specific, respectively. In the 4239 common-expressed proteins, 80 and 74 were preferentially accumulated in WT and fl, respectively. Pathway enrichment analysis and protein–protein interaction network analysis of both variety-specific and differential abundant proteins showed that secondary metabolism and chloroplast-related pathways were significantly enriched. Quantitative real-time PCR confirmed that the expression levels of 12 out of 16 selected genes from representative pathways were consistent with their protein accumulation patterns. Further analyses showed that the content of chlorophyll a in WT, but not chlorophyll b, was significantly increased compared to fl. This work provides the leaf proteome profiles of Xu142 and its fl mutant, indicating the necessity of further investigation of molecular differences between WT and fl leaves.

Molecular identification and functional characterization of GhAMT1.3 in ammonium transport with a high affinity from cotton (Gossypium hirsutum L.)

Ammonium (NH₄ ⁺ ) represents a primary nitrogen source for many plants, its effective transport into and between tissues and further assimilation in cells determine greatly plant nitrogen use efficiency. However, biological components involved in NH₄ ⁺ movement in woody plants are unclear. Here, we report kinetic evidence for cotton NH₄ ⁺ uptake and molecular identification of certain NH₄ ⁺ transporters (AMTs) from cotton (Gossypium hirustum). A substrate-influx assay using ¹⁵ N-isotope revealed that cotton possessed a high-affinity transport system with a K_m of 58 μM for NH₄ ⁺ . Sequence analysis showed that GhAMT1.1-1.3 encoded respectively a membrane protein containing 485, 509 or 499 amino acids. Heterologous functionality test demonstrated that GhAMT1.1-1.3 expression mediated NH₄ ⁺ permeation across the plasma membrane (PM) of yeast and/or Arabidopsis qko-mutant cells, allowing a growth restoration of both mutants on NH₄ ⁺ . Quantitative PCR measurement showed that GhAMT1.3 was expressed in roots and leaves and markedly up-regulated by N-starvation, repressed by NH₄ ⁺ resupply and regulated diurnally and age-dependently, suggesting that GhAMT1.3 should be a N-responsive gene. Importantly, GhAMT1.3 expression in Arabidopsis improved plant growth on NH₄ ⁺ and enhanced total nitrogen accumulation (∼50% more), conforming with the observation of 2-fold more NH₄ ⁺ absorption by GhAMT1.3-transformed qko plant roots during a 1-h root influx period. Together with its targeting to the PM and saturated transport kinetics with a K_m of 72 μM for NH₄ ⁺ , GhAMT1.3 is suggested to be a high-affinity NH₄ ⁺ permease that may play a significant role in cotton NH₄ ⁺ acquisition and utilization, adding a new member in the plant AMT family.

Genome-Wide Analysis of Cotton Auxin Early Response Gene Families and Their Roles in Somatic Embryogenesis

Auxin is well known to regulate growth and development processes. Auxin early response genes serve as a critical component of auxin signaling and mediate auxin regulation of diverse physiological processes. In the present study, a genome-wide identification and comprehensive analysis of auxin early response genes were conducted in upland cotton. A total of 71 auxin response factor (ARF), 86 Auxin/Indole-3-Acetic Acid (Aux/IAA), 63 Gretchen Hagen3 (GH3), and 194 small auxin upregulated RNA (SAUR) genes were identified in upland cotton, respectively. Phylogenetic analysis revealed that the ARF, GH3, and SAUR families were likely subject to extensive evolutionary divergence between Arabidopsis and upland cotton, while the Aux/IAA family was evolutionary conserved. Expression profiles showed that the ARF, Aux/IAA, GH3, and SAUR family genes were extensively involved in embryogenic competence acquisition of upland cotton callus. The Aux/IAA family genes generally showed a higher expression level in the non-embryogenic callus (NEC) of highly embryogenic cultivar CCRI24 than that of recalcitrant cultivar CCRI12, which may be conducive to initializing the embryogenic transformation. Auxin early response genes were tightly co-expressed with most of the known somatic embryogenesis (SE) related genes, indicating that these genes may regulate upland cotton SE by interacting with auxin early response genes.

GhVTC1, the Key Gene for Ascorbate Biosynthesis in Gossypium hirsutum, Involves in Cell Elongation Under Control of Ethylene

L-Ascorbate (Asc) plays important roles in cell growth and plant development, and its de novo biosynthesis was catalyzed by the first rate-limiting enzyme VTC1. However, the function and regulatory mechanism of VTC1 involved in cell development is obscure in Gossypium hirsutum. Herein, the Asc content and AsA/DHA ratio were accumulated and closely linked with fiber development. The GhVTC1 encoded a typical VTC1 protein with functional conserved domains and expressed preferentially during fiber fast elongation stages. Functional complementary analysis of GhVTC1 in the loss-of-function Arabidopsis vtc1-1 mutants indicated that GhVTC1 is genetically functional to rescue the defects of mutants to normal or wild type (WT). The significant shortened primary root in vtc1-1 mutants was promoted to the regular length of WT by the ectopic expression of GhVTC1 in the mutants. Additionally, GhVTC1 expression was induced by ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC), and the GhVTC1 promoter showed high activity and included two ethylene-responsive elements (ERE). Moreover, the 5'-truncted promoters containing the ERE exhibited increased activity by ACC treatment. Our results firstly report the cotton GhVTC1 function in promoting cell elongation at the cellular level, and serve as a foundation for further understanding the regulatory mechanism of Asc-mediated cell growth via the ethylene signaling pathway.

Updates on molecular mechanisms in the development of branched trichome in Arabidopsis and nonbranched in cotton

Trichomes are specialized epidermal cells and a vital plant organ that protect plants from various harms and provide valuable resources for plant development and use. Some key genes related to trichomes have been identified in the model plant Arabidopsis thaliana through glabrous mutants and gene cloning, and the hub MYB-bHLH-WD40, consisting of several factors including GLABRA1 (GL1), GL3, TRANSPARENT TESTA GLABRA1 (TTG1), and ENHANCER OF GLABRA3 (EGL3), has been established. Subsequently, some upstream transcription factors, phytohormones and epigenetic modification factors have also been studied in depth. In cotton, a very important fibre and oil crop globally, in addition to the key MYB-like factors, more important regulators and potential molecular mechanisms (e.g. epigenetic modifiers, distinct metabolic pathways) are being exploited during different fibre developmental stages. This occurs due to increased cotton research, resulting in the discovery of more complex regulation mechanisms from the allotetraploid genome of cotton. In addition, some conservative as well as specific mediators are involved in trichome development in other species. This study summarizes molecular mechanisms in trichome development and provides a detailed comparison of the similarities and differences between Arabidopsis and cotton, analyses the possible reasons for the discrepancy in identification of regulators, and raises future questions and foci for understanding trichome development in more detail.

Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

BACKGROUND

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

RESULTS

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

CONCLUSION

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

Transcriptome analysis reveals differences in the mechanisms of fiber initiation and elongation between long- and short-fiber cotton (Gossypium hirsutum L.) lines

Background

Improving the yield and fiber quality of upland cotton is a goal of plant breeders. However, increasing the yield and quality of cotton fibers is becoming more urgent. While the growing human population needs more cotton fiber, climate change is reducing the amount of land on which cotton can be planted, or making it difficult to ensure that water and other resources will be available in optimal quantities. The most logical means of improving yield and quality is understanding and manipulating the genes involved. Here, we used comparative transcriptomics to explore differences in gene expression between long- and short-fiber cotton lines to identify candidate genes useful for cotton improvement.

Results

Light and electron microscopy revealed that the initial fiber density was significantly greater in our short-fiber group (SFG) than in our long-fiber group (LFG). Compared with the SFG fibers, the LFG fibers were longer at all developmental stages. Comparison of the LFG and SFG transcriptomes revealed a total of 3538 differentially expressed genes (DEGs). Notably, at all three developmental stages examined, two expression patterns, consistently downregulated (profile 0) and consistently upregulated (profile 7), were identified, and both were significantly enriched in the SFG and LFG. Twenty-two DEGs known to be involved in fiber initiation were detected in profile 0, while 31 DEGs involved in fiber elongation were detected in profile 7. Functional annotation suggested that these DEGs, which included ERF1, TUA2, TUB1, and PER64, affect fiber elongation by participating in the ethylene response, microtubule synthesis, and/or the peroxidase (POD) catalytic pathway. qRT-PCR was used to confirm the RNA sequencing results for select genes.

Conclusions

A comparison of SFG and LFG transcription profiles revealed modest but important differences in gene expression between the groups. Notably, our results confirm those of previous studies suggesting that genes involved in ethylene, tubulin, and POD pathways play important roles in fiber development. The 22 consistently downregulated DEGs involved in fiber initiation and the 31 consistently upregulated genes involved in fiber elongation are seemingly good candidate genes for improving fiber initiation and elongation in cotton.

Electronic supplementary material

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

Comprehensive analyses of ZFP gene family and characterization of expression profiles during plant hormone response in cotton

BACKGROUND

Zinc finger proteins (ZFPs) containing only a single zinc finger domain play important roles in the regulation of plant growth and development, as well as in biotic and abiotic stress responses. To date, the evolutionary history and functions of the ZFP gene family have not been identified in cotton.

RESULTS

In this paper, we identified 29 ZFP genes in Gossypium hirsutum. This gene family was divided into seven subfamilies, 22 of which were distributed over 17 chromosomes. Bioinformatic analysis revealed that 20 GhZFP genes originated from whole genome duplications and two originated from dispersed duplication events, indicating that whole genome duplication is the main force in the expansion of the GhZFP gene family. Most GhZFP8 subfamily genes, except for GhZFP8-3, were highly expressed during fiber cell growth, and were induced by brassinosteroids in vitro. Furthermore, we found that a large number of GhZFP genes contained gibberellic acid responsive elements, auxin responsive elements, and E-box elements in their promoter regions. Exogenous application of these hormones significantly stimulated the expression of these genes.

CONCLUSIONS

Our findings reveal that GhZFP8 genes are involved in cotton fiber development and widely induced by auxin, gibberellin and BR, which provides a foundation for the identification of more downstream genes with potential roles in phytohormone stimuli, and a basis for breeding better cotton varieties in the future.

Comprehensive analysis of WOX genes uncovers that WOX13 is involved in phytohormone-mediated fiber development in cotton

BACKGROUND

The WOX (WUSCHEL-RELATED HOMEOBOX) gene family encodes a class of transcription factors that are unique to green plants, where they are involved in regulating the development of plant tissues and organs by determining cell fate. Although the importance of the WOX gene is well known, there are few studies describing their functions in cotton.

RESULTS

In this study, 32 WOX genes were found in Gossypium hirsutum. Phylogenetic analysis showed that WOX proteins of cotton can be divided into three clades: the ancient, intermediate, and WUS clades. The number of WOX proteins in the WUS clade was greater than the sum of the proteins in the other two clades. Our analysis revealed that 20 GhWOX genes are distributed on 16 cotton chromosomes and that duplication events are likely to have contributed to the expansion of the GhWOX family. All GhWOX genes have introns, and each GhWOX protein contains multiple motifs. RNA-seq data and real-time PCR showed that GhWOX13 gene subfamily is specifically expressed at a high level in cotton fibers. We also identified putative GA, NAA, and BR response elements in the promoter regions of the GhWOX13 genes and GhWOX13 transcripts were significantly induced by GA, NAA, and BR.

CONCLUSIONS

Our data provides a useful resource for future studies on the functional roles of cotton WOX genes and shows that the GhWOX13 genes may influence cotton fiber development. Our results also provide an approach for identifying and characterizing WOX protein genes in other species.

1-Aminocyclopropane-1-Carboxylic Acid Oxidase (ACO): The Enzyme That Makes the Plant Hormone Ethylene

The volatile plant hormone ethylene regulates many plant developmental processes and stress responses. It is therefore crucial that plants can precisely control their ethylene production levels in space and time. The ethylene biosynthesis pathway consists of two dedicated steps. In a first reaction, S-adenosyl-L-methionine (SAM) is converted into 1-aminocyclopropane-1-carboxylic acid (ACC) by ACC-synthase (ACS). In a second reaction, ACC is converted into ethylene by ACC-oxidase (ACO). Initially, it was postulated that ACS is the rate-limiting enzyme of this pathway, directing many studies to unravel the regulation of ACS protein activity, and stability. However, an increasing amount of evidence has been gathered over the years, which shows that ACO is the rate-limiting step in ethylene production during certain dedicated processes. This implies that also the ACO protein family is subjected to a stringent regulation. In this review, we give an overview about the state-of-the-art regarding ACO evolution, functionality and regulation, with an emphasis on the transcriptional, post-transcriptional, and post-translational control. We also highlight the importance of ACO being a prime target for genetic engineering and precision breeding, in order to control plant ethylene production levels.

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.

Modification of phytosterol composition influences cotton fiber cell elongation and secondary cell wall deposition

BACKGROUND

Cotton fiber is a single cell that arises from the epidermis of ovule. It is not only a main economic product of cotton, but an ideal material for studying on the growth and development of plant cell. Our previous study indicated that phytosterol content and the ratio of campesterol to sitosterol fluctuated regularly in cotton fiber development. However, what effects of modified phytosterol content and composition on the growth and development of cotton fiber cell is unknown. In this study, we overexpressed the GhSMT2-1, a cotton homologue of sterol C-24 methyltransferase 2 gene in transgenic upland cotton plants to modify phytosterol content and composition in fiber cells and investigated the changes on fiber elongation and secondary cell wall deposition.

RESULTS

GhSMT2-1 overexpression led to changes of phytosterol content and the ratio of campesterol to sitosterol in fiber cell. At the rapid elongation stage of fiber cell, total phytosterol and sitosterol contents were increased while campesterol content was decreased in transgenic fibers when compared to control fibers. Accordingly, the ratio of campesterol to sitosterol declined strikingly. Simultaneously, the transgenic fibers were shorter and thicker than control fibers. Exogenous application of sitosterol or campesterol separately inhibited control fiber cell elongation in cotton ovule culture system in vitro. In addition, campesterol treatment partially rescued transgenic fiber elongation.

CONCLUSION

These results elucidated that modification of phytosterol content and composition influenced fiber cell elongation and secondary cell wall formation. High sitosterol or low ratio of campesterol to sitosterol suppresses fiber elongation and/or promote secondary cell wall deposition. The roles of sitosterol and campesterol were discussed in fiber cell development. There might be a specific ratio of campesterol to sitosterol in different developmental stage of cotton fibers, in which GhSMT2-1 play an important role. Our study, at a certain degree, provides novel insights into the regulatory mechanisms of fiber cell development.

Role of Cytokinin, Strigolactone, and Auxin Export on Outgrowth of Axillary Buds in Apple

Shoot branching is regulated by phytohormones, including cytokinin (CK), strigolactone (SL), and auxin in axillary buds. The correlative importance of these phytohormones in the outgrowth of apple axillary buds remains unclear. In this study, the outgrowth dynamics of axillary buds of a more-branching mutant (MB) and its wild-type (WT) of Malus spectabilis were assessed using exogenous chemical treatments, transcriptome analysis, paraffin section, and reverse transcription-quantitative PCR analysis (RT-qPCR). High contents of CK and abscisic acid coincided in MB axillary buds. Exogenous CK promoted axillary bud outgrowth in the WT but not in MB, whereas exogenous gibberellic had no significant effect on bud outgrowth in the WT. Functional analysis of transcriptome data and RT-qPCR analysis of gene transcripts revealed that MB branching were associated with CK signaling, auxin transport, and SL signaling. Transcription of the SL-related genes MsMAX1, MsD14, and MsMAX2 in the axillary buds of MB was generally upregulated during bud outgrowth, whereas MsBRC1/2 were generally downregulated both in WT and MB. Exogenous SL inhibited outgrowth of axillary buds in the WT and the apple varieties T337, M26, and Nagafu 2, whereas axillary buds of the MB were insensitive to SL treatment. Treatment with N-1-naphthylphalamic acid (NPA; an auxin transport inhibitor) inhibited bud outgrowth in plants of the WT and MB. The transcript abundance of MsPIN1 was generally decreased in response to NPA and SL treatments, and increased in CK and decapitation treatments, whereas no consistent pattern was observed for MsD14 and MsMAX2. Collectively, the present results suggest that in apple auxin transport from the axillary bud to the stem may be essential for the outgrowth of axillary buds, and at least, is involved in the process of bud outgrowth.

Cotton fiber elongation requires the transcription factor GhMYB212 to regulate sucrose transportation into expanding fibers

Cotton is white gold across the globe and composed of fiber cells derived from the outer integument of cotton ovules. Fiber elongation uses sucrose as a direct carbon source. The molecular mechanism transcriptionally controlling sucrose transport from ovules into the elongating fibers remains elusive. In this study the involvement of GhMYB212 in the regulation of sucrose transportion into expanding fibers was investigated. GhMYB212 RNAi plants (GhMYB212i) accumulated less sucrose and glucose in developing fibers, and had shorter fibers and a lower lint index. RNA-seq and protein-DNA binding assays revealed that GhMYB212 was closely linked to the pathways of sucrose and starch transportation and metabolism, directly controling the expression of a sucrose transporter gene GhSWEET12. GhSWEET12 RNAi plants (GhSWEET12i) possessed similar fiber phenotypes to those of GhMYB212i. Exogenous sucrose supplementation in ovule cultures did not rescue the shorter fiber phenotype of GhMYB212i and GhSWEET12i. This finding supported the idea that the attenuated rate of sucrose transport from the outer seed coat into the fibers is responsible for the retardation of fiber elongation. Current investigations support the idea that GhMYB212 functions as the main regulator of fiber elongation by controlling the expression of GhSWEET12, and therefore it is important to study cell expansion and sugar transportation during seed development.

CtACO1 Overexpression Resulted in the Alteration of the Flavonoids Profile of Safflower

Background: Flavonoids with various structures play a vital role in plant acclimatization to varying environments as well as in plant growth, development, and reproduction. Exogenous applications of ethylene and 1-aminocyclopropane carboxylic acid (ACC), could affect the accumulation of flavonoids. Very few attempts have been made to investigate the effect of 1-aminocyclopropane carboxylic acid oxidase (ACO), a unique enzyme that catalyzes ACC to ethylene, on genes and metabolites in the flavonoid biosynthetic pathway. In this study, two ACOs in safflower (CtACOs) were cloned, and then transgenic safflower with overexpressed CtACO1 was generated through the Agrobacterium-mediated floral dipping method. Results: CtACO1 and CtACO2 were both characterized by the 2-oxoglutarate binding domain RxS and the ferrous iron binding site HxDxnH as ACOs from other plants. However, the transcript levels of CtACO1 in flowers at stages I, II, III, and IV were all higher than those of CtACO2. At the cellular level, by using electroporation transformation, CtACO1 was found to be localized at the cytomembrane in onion epidermal cells. CtACO1 overexpression had varying effects on genes involved in the ethylene and flavonoid biosynthetic pathways. The metabolites analysis showed that CtACO1 overexpression lines had a higher accumulation of quercetin and its glycosylated derivatives (quercetin 3-β-d-glucoside and rutin). In contrast, the accumulation of quinochalcones (hydroxysafflor yellow A and carthamin), kaempferol glycosylated derivatives (kaempferol-3-O-β-rutinoside and kaempferol-3-O-β-d-glucoside), apigenin, and luteolin in CtACO1 overexpression lines were decreased. Conclusion: This study confirmed the feasibility of applying the floral dipping method to safflower and showed a novel regulatory effect of CtACO1 in the flavonoid biosynthetic pathway. It provides hypothetical and practical groundwork for further research on regulating the overall metabolic flux of flavonoids in safflower, particularly hydroxysafflor yellow A and other quinochalcones, by using appropriate genetic engineering strategies.

Feasibility and Performance of Loop-Mediated Isothermal Amplification Assay in the Diagnosis of Pulmonary Tuberculosis in Decentralized Settings in Eastern China

Early diagnosis is essential for the control and prevention of tuberculosis (TB). The objective of this study was to investigate the feasibility and performance of loop-mediated isothermal amplification (LAMP) in the diagnosis of pulmonary TB in county-level microscopy centers in Qingdao, Eastern China. A total of 523 presumptive TB patients were consecutively recruited between July 2017 and April 2018, and 22 patients were excluded from the analysis. Of 102 culture-positive cases, TB-LAMP identified 91 cases, demonstrating a sensitivity of 89.2%. In comparison, the sensitivity of routine smear microscopy was 69.6% (71/102), which was significantly lower than that of TB-LAMP (P=0.001). In addition, TB-LAMP sensitivities in smear-positive and smear-negative samples were 98.6% and 67.7%, respectively. In conclusion, our data demonstrate that TB-LAMP outperforms conventional smear microscopy in TB diagnosis, which could be used as an alternative method for smear microscopy in resource-limited settings in China.

Ethylene Treatment in Studying Leaf Senescence in Arabidopsis

Plant leaf senescence is the final process of leaf development that involves the mobilization of nutrients from old leaves to newly growing tissues. Leaf senescence involves a coordinated action at the cellular and organism levels under the control of a highly regulated genetic program. However, leaf senescence is also influenced by multiple internal and environmental signals that are integrated into the age information. Among these internal factors, the simple gaseous phytohormone ethylene is well-known to be an important endogenous modulator of plant leaf senescence and fruit ripening. Ethylene and its biosynthetic precursor ACC (1-aminocyclopropane-l-carboxylic acid) as well as inhibitors of ethylene biosynthesis aminoethoxyvinylglycine (AVG) and action (AgNO₃) or the inhibitor of ethylene perception 1-MCP (1-methylcyclopropene) have been widely used by the research community of plant leaf senescence. However, until now, no systemically experimental method about the usage of ethylene in studying the molecular mechanisms of leaf senescence in Arabidopsis is available. Here, we provide detailed methods for exogenous application of ethylene and its inhibitors in studying Arabidopsis leaf senescence, which has been successfully used in our laboratory.

Fiber Quality Improvement in Upland Cotton (Gossypium hirsutum L.): Quantitative Trait Loci Mapping and Marker Assisted Selection Application

Genetic improvement in fiber quality is one of the main challenges for cotton breeders. Fiber quality traits are controlled by multiple genes and are classified as complex quantitative traits, with a negative relationship with yield potential, so the genetic gain is low in traditional genetic improvement by phenotypic selection. The availability of Gossypium genomic sequences facilitates the development of high-throughput molecular markers, quantitative trait loci (QTL) fine mapping and gene identification, which helps us to validate candidate genes and to use marker assisted selection (MAS) on fiber quality in breeding programs. Based on developments of high density linkage maps, QTLs fine mapping, marker selection and omics, we have performed trait dissection on fiber quality traits in diverse populations of upland cotton. QTL mapping combined with multi-omics approaches such as, RNA sequencing datasets to identify differentially expressed genes have benefited the improvement of fiber quality. In this review, we discuss the application of molecular markers, QTL mapping and MAS for fiber quality improvement in upland cotton.

A Pivotal Role of Hormones in Regulating Cotton Fiber Development

Cotton is the main source of renewable fiber in the world and is primarily used for textile production. Cotton fibers are single cells differentiated from the ovule epidermis and are an excellent model system for studying cell elongation, polyploidization, and cell wall biosynthesis. Plant hormones, which are present in relatively low concentrations, play important roles in various developmental processes, and recently, multiple reports have revealed the pivotal roles of hormones in regulating cotton fiber development. For example, exogenous application of hormones has been shown to promote the initiation and growth of fiber cells. However, a comprehensive understanding about phytohormone regulating fiber development is still unknown. Here, we focus on recent advances in elucidating the roles of multiple phytohormones in the control of fiber development, namely auxin, gibberellin, brassinosteroid, ethylene, cytokinin, abscisic acid, and strigolactones. We not only review the identification of genes involved in hormone biosynthetic and signaling pathways but also discuss the mechanisms of these phytohormones in regulating the initiation and elongation of fiber cells in cotton. Auxin, gibberellin, brassinosteroid, ethylene, jasmonic acid, and strigolactones play positive roles in fiber development, whereas cytokinin and abscisic acid inhibit fiber growth. Our aim is to provide a comprehensive review of the role of phytohormones in cotton fiber development that will serve as the basis for further elucidation of the mechanisms by which plant hormones regulate fiber growth.

Transcriptome Profiling of the Elongating Internode of Cotton (Gossypium hirsutum L.) Seedlings in Response to Mepiquat Chloride

The plant growth retardant mepiquat chloride (MC) has been extensively used to produce compact plant canopies and increase yield in cotton (Gossypium hirsutum L.). Previous studies showed that MC reduced plant height and internode length by inhibiting GA biosynthesis and cell elongation. However, whether there are other molecular mechanisms underlying MC-induced growth retardation has remained largely unknown. In the present study, we conducted histological, transcriptomic, and phytohormone analyses of the second elongating internodes of cotton seedlings treated with MC. Histological analysis revealed that the MC shortened the internodes through suppressing both cell division and cell elongation. Consistent with the observed phenotype, many genes related to cell growth were significantly downregulated by MC. Transcriptome profiling showed that the expression of genes related not only to GA, but also to auxin, brassinosteroid (BR), and ethylene metabolism and signaling was remarkably suppressed, whereas that of genes related to cytokinin (CK) and abscisic acid (ABA) metabolism was induced by MC. Consistent with the expression pattern, significant decrease of endogenous GA, auxin, and BR content, but an increase in CK content was observed after MC treatment. Most of these hormone related genes displayed opposite regulation pattern by exogenous GA₃ treatment compared to MC; simultaneous application of MC and GA₃ could alleviate the genes expression changes induced by MC treatment, indicating MC does not directly affect other plant hormones, but through the inhibition of the GA biosynthesis. In addition, the expression of genes related to secondary metabolism and many transcription factors (TFs) were differentially regulated by MC. In summary, we confirmed the important role of GA in MC-induced growth inhibition of cotton, and further found that other hormones were also involved in this process in a GA-dependent manner. This study provides novel insights into the molecular mechanism underlying the MC-mediated inhibition of internode elongation in cotton seedlings.

Genome-wide characterization and phylogenetic analysis of GSK gene family in three species of cotton: evidence for a role of some GSKs in fiber development and responses to stress

BACKGROUND

The glycogen synthase kinase 3/shaggy kinase (GSK3) is a serine/threonine kinase with important roles in animals. Although GSK3 genes have been studied for more than 30 years, plant GSK genes have been studied only since the last decade. Previous research has confirmed that plant GSK genes are involved in diverse processes, including floral development, brassinosteroid signaling, and responses to abiotic stresses.

RESULT

In this study, 20, 15 (including 5 different transcripts) and 10 GSK genes were identified in G. hirsutum, G. raimondii and G. arboreum, respectively. A total of 65 genes from Arabidopsis, rice, and cotton were classified into 4 clades. High similarities were found in GSK3 protein sequences, conserved motifs, and gene structures, as well as good concordance in gene pairwise comparisons (G. hirsutum vs. G. arboreum, G. hirsutum vs. G. raimondii, and G. arboreum vs. G. raimondii) were observed. Whole genome duplication (WGD) within At and Dt sub-genomes has been central to the expansion of the GSK gene family. Furthermore, GhSK genes showed diverse expression patterns in various tissues. Additionally, the expression profiles of GhSKs under different stress treatments demonstrated that many are stress-responsive genes. However, none were induced by brassinolide treatment. Finally, nine co-expression sub-networks were observed for GhSKs and the functional annotations of these genes suggested that some GhSKs might be involved in cotton fiber development.

CONCLUSION

In this present work, we identified 45 GSK genes from three cotton species, which were divided into four clades. The gene features, muti-alignment, conversed motifs, and syntenic blocks indicate that they have been highly conserved during evolution. Whole genome duplication was determined to be the dominant factor for GSK gene family expansion. The analysis of co-expressed sub-networks and tissue-specific expression profiles suggested functions of GhSKs during fiber development. Moreover, their different responses to various abiotic stresses indicated great functional diversity amongst the GhSKs. Briefly, data presented herein may serve as the basis for future functional studies of GhSKs.

Mitochondrial alternative oxidase-dependent autophagy involved in ethylene-mediated drought tolerance in Solanum lycopersicum

Mitochondrial alternative oxidase (AOX) is involved in a large number of plant physiological processes, such as growth, development and stress responses; however, the exact role of AOX in response to drought remains unclear. In our study, we provide solid evidences that the activated AOX capacity positively involved in ethylene-induced drought tolerance, in tomato (Solanum lycopersicum), accompanied by the changing level of hydrogen peroxide (H2 O2 ) and autophagy. In AOX1a-RNAi plants, the ethylene-induced drought tolerance was aggravated and associated with decreasing level of autophagy. The H2 O2 level was relatively higher in AOX1a-RNAi plants, whereas it was lower in AOX1a-overexpressing (35S-AOX1a-OE) plants after 1-(aminocarbonyl)-1-cyclopropanecarboxylic acid (ACC) pretreatment in the 14th day under drought stress. Interestingly, the accumulation of autophagosome was accompanied by the changing level of reactive oxygen species (ROS) in AOX transgenic tomato under drought stress whether or not pretreated with ACC. Pharmacological scavenging of H2 O2 accumulation in AOX1a-RNAi (aox19) stimulated autophagy acceleration under drought stress, and it seems that AOX-dependent ROS signalling is critical in triggering autophagy. Lower levels of ROS signalling positively induce autophagy activity, whereas higher ROS level would lead to rapid programmed cell death (PCD), especially in ethylene-mediated drought tolerance. Moreover, ethylene-induced autophagy during drought stress also can be through ERF5 binding to the promoters of ATG8d and ATG18h. These results demonstrated that AOX plays an essential role in ethylene-induced drought tolerance and also played important roles in mediating autophagy generation via balancing ROS level.

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.

Genome-wide analysis of BES1 genes in Gossypium revealed their evolutionary conserved roles in brassinosteroid signaling

Brassinosteroids (BRs), which are essential phytohormones for plant growth and development, are important for cotton fiber development. Additionally, BES1 transcription factors are critical for BR signal transduction. However, cotton BES1 family genes have not been comprehensively characterized. In this study, we identified 11 BES1 genes in G. arboreum, 11 in G. raimondii, 16 in G. barbadense, and 22 in G. hirsutum. The BES1 sequences were significantly conserved in the Arabidopsis thaliana, rice, and upland cotton genomes. A total of 94 BES1 genes from 10 different plant species were divided into three clades according to the neighbor-joining and minimum-evolution methods. Moreover, the exon/intron patterns and motif distributions were highly conserved among the A. thaliana and cotton BES1 genes. The collinearity among the orthologs from the At and Dt subgenomes was estimated. Segmental duplications in the At and Dt subgenomes were primarily responsible for the expansion of the cotton BES1 gene family. Of the GhBES1 genes, GhBES1.4_At/Dt exhibited BL-induced expression and was predominantly expressed in fibers. Furthermore, Col-0/mGhBES1.4_At plants produced curled leaves with long and bent petioles. These transgenic plants also exhibited decreased hypocotyl sensitivity to brassinazole and constitutive BR induced/repressed gene expression patterns. The constitutive BR responses of the plants overexpressing mGhBES1.4_At were similar to those of the bes1-D mutant.

In Vitro Activity of PBTZ169 against Multiple Mycobacterium Species

In this study, we demonstrate that PBTZ169 exhibits significant differences in in vitro activity against multiple Mycobacterium species. The amino acid polymorphism at codon 387 of decaprenylphosphoryl-beta-d-ribose oxidase (DprE1) can be used as a surrogate marker for in vitro susceptibility to PBTZ169 in mycobacteria. In addition, the amino acid substitution at codon 154 in DprE1 may be associated with acquired resistance to PBTZ169 in the Mycobacterium fortuitum mutants.

TaWRKY51 promotes lateral root formation through negative regulation of ethylene biosynthesis in wheat (Triticum aestivum L.)

Common wheat (Triticum aestivum L.) is an important staple food crop worldwide. Lateral roots (LRs), as the major component of root architecture, affect water and nutrient uptake in wheat. The phytohormone ethylene is known to affect LR formation; however, the factor(s) modulating ethylene during this process have not yet been elucidated in wheat. Here we identified wheat TaWRKY51 as a key factor that functions in LR formation by modulating ethylene biosynthesis. Wheat TaWRKY51RNA interference lines (TaWRKY51-RNAi) and the homozygous mutants tawrky51-2a and tawrky51-2b all produced fewer LRs than the wild type and negative transgenic plants, whereas the TaWRKY51 overexpression lines (TaWRKY51-OE) had the opposite phenotype. Transcription analysis revealed that 1-aminocyclopropane-1-carboxylic acid synthase (ACS) genes (TaACS2, TaACS7 and TaACS8) involved in ethylene biosynthesis were downregulated in TaWRKY51-OE lines but upregulated in TaWRKY51-RNAi lines. The rate of ethylene production also decreased in TaWRKY51-OE lines but increased in TaWRKY51-RNAi lines compared with their respective negative transgenic controls. Electrophoretic mobility shift and transient expression assays revealed that TaWRKY51 inhibits the expression of ACS genes by binding to the W-box cis-element present in their promoter region. Moreover, overexpression of ACS2 or exogenous application of 1-aminocyclopropane-1-carboxylic acid reversed the phenotype of enhanced LR number in TaWRKY51-OE Arabidopsis lines, and overexpression of TaWRKY51 in the ethylene-overproducing mutant eto1-1 rescued its LR defect phenotype. In addition, genetic evidence demonstrates that TaWRKY51-regulated LR formation is also dependent on ethylene and auxin signaling pathways. Our findings reveal a molecular genetic mechanism by which a WRKY gene coordinates ethylene production and LR formation in wheat.

Role of GhHDA5 in H3K9 deacetylation and fiber initiation in Gossypium hirsutum

Cotton fibers are single-celled trichomes that initiate from the epidermal cells of the ovules at or before anthesis. Here, we identified that the histone deacetylase (HDAC) activity is essential for proper cotton fiber initiation. We further identified 15 HDACs homoeologs in each of the A- and D-subgenomes of Gossypium hirsutum. Few of these HDAC homoeologs expressed preferentially during the early stages of fiber development [-1, 0 and 6 days post-anthesis (DPA)]. Among them, GhHDA5 expressed significantly at the time of fiber initiation (-1 and 0 DPA). The in vitro assay for HDAC activity indicated that GhHDA5 primarily deacetylates H3K9 acetylation marks. Moreover, the reduced expression of GhHDA5 also suppresses fiber initiation and lint yield in the RNA interference (RNAi) lines. The 0 DPA ovules of GhHDA5RNAi lines also showed alterations in reactive oxygen species homeostasis and elevated autophagic cell death in the developing fibers. The differentially expressed genes (DEGs) identified through RNA-seq of RNAi line (DEP12) and their pathway analysis showed that GhHDA5 modulates expression of many stress and development-related genes involved in fiber development. The reduced expression of GhHDA5 in the RNAi lines also resulted in H3K9 hyper-acetylation on the promoter region of few DEGs assessed by chromatin immunoprecipitation assay. The positively co-expressed genes with GhHDA5 showed cumulative higher expression during fiber initiation, and gene ontology annotation suggests their involvement in fiber development. Furthermore, the predicted protein interaction network in the positively co-expressed genes indicates HDA5 modulates fiber initiation-specific gene expression through a complex involving reported repressors.

Genome-wide identification of the GhARF gene family reveals that GhARF2 and GhARF18 are involved in cotton fibre cell initiation

Auxin signalling plays an essential role in regulating plant development. Auxin response factors (ARFs), which are critical components of auxin signalling, modulate the expression of early auxin-responsive genes by binding to auxin response factor elements (AuxREs). However, there has been no comprehensive characterization of this gene family in cotton. Here, we identified 56 GhARF genes in the assembled Gossypium hirsutum genome. This gene family was divided into 17 subfamilies, and 44 members of them were distributed across 21 chromosomes. GhARF6 and GhARF11 subfamily genes were predominantly expressed in vegetative tissues, whereas GhARF2 and GhARF18 subfamily genes were highly expressed during seed fibre cell initiation. GhARF2-1 and GhARF18-1 were exclusively expressed in trichomes, organs similar to cotton seed fibre cells, and overexpression of these genes in Arabidopsis enhances trichome initiation. Comparative transcriptome analysis combined with AuxRE prediction revealed 11 transcription factors as potential target genes of GhARF2 and GhARF18. Six of these genes were significantly expressed during seed fibre cell initiation and were bound by GhARF2-1 and GhARF18-1 in yeast one-hybrid assays. Our results suggest that GhARF2 and GhARF18 genes may be key regulators of cotton seed fibre initiation by regulating the expression of several transcription factor genes. This study deepens our understanding of auxin-mediated initiation of cotton seed fibre cells and helps us in breeding better cotton varieties in the future.

Resequencing a core collection of upland cotton identifies genomic variation and loci influencing fiber quality and yield

Upland cotton is the most important natural-fiber crop. The genomic variation of diverse germplasms and alleles underpinning fiber quality and yield should be extensively explored. Here, we resequenced a core collection comprising 419 accessions with 6.55-fold coverage depth and identified approximately 3.66 million SNPs for evaluating the genomic variation. We performed phenotyping across 12 environments and conducted genome-wide association study of 13 fiber-related traits. 7,383 unique SNPs were significantly associated with these traits and were located within or near 4,820 genes; more associated loci were detected for fiber quality than fiber yield, and more fiber genes were detected in the D than the A subgenome. Several previously undescribed causal genes for days to flowering, fiber length, and fiber strength were identified. Phenotypic selection for these traits increased the frequency of elite alleles during domestication and breeding. These results provide targets for molecular selection and genetic manipulation in cotton improvement.

Flavonoid biosynthesis controls fiber color in naturally colored cotton

The existence of only natural brown and green cotton fibers (BCF and GCF, respectively), as well as poor fiber quality, limits the use of naturally colored cotton (Gossypium hirsutum L.). A better understanding of fiber pigment regulation is needed to surmount these obstacles. In this work, transcriptome analysis and quantitative reverse transcription PCR revealed that 13 and 9 phenylpropanoid (metabolic) pathway genes were enriched during pigment synthesis, while the differential expression of phenylpropanoid (metabolic) and flavonoid metabolic pathway genes occurred among BCF, GCF, and white cotton fibers (WCF). Silencing the chalcone flavanone isomerase gene in a BCF line resulted in three fiber phenotypes among offspring of the RNAi lines: BCF, almost WCF, and GCF. The lines with almost WCF suppressed chalcone flavanone isomerase, while the lines with GCF highly expressed the glucosyl transferase (3GT) gene. Overexpression of the Gh3GT or Arabidopsis thaliana 3GT gene in BCF lines resulted in GCF. Additionally, the phenylpropanoid and flavonoid metabolites of BCF and GCF were significantly higher than those of WCF as assessed by a metabolomics analysis. Thus, the flavonoid biosynthetic pathway controls both brown and green pigmentation processes. Like natural colored fibers, the transgenic colored fibers were weaker and shorter than WCF. This study shows the potential of flavonoid pathway modifications to alter cotton fibers’ color and quality.

A Gossypium hirsutum GDSL lipase/hydrolase gene (GhGLIP) appears to be involved in promoting seed growth in Arabidopsis

GDSL lipase (GLIP) plays a pivotal role in plant cell growth as a multifunctional hydrolytic enzyme. Herein, a cotton (Gossypium hirsutum L. cv Xuzhou 142) GDSL lipase gene (GhGLIP) was obtained from developing ovules and fibers. The GhGLIP cDNA contained an open reading frame (ORF) of 1,143 base pairs (bp) and encodes a putative polypeptide of 380 amino acid residues. Sequence alignment indicated that GhGLIP includes four enzyme catalytic amino acid residue sites of Ser (S), Gly (G), Asn (N) and His (H), located in four conserved blocks. Phylogenetic tree analysis showed that GhGLIP belongs to the typical class IV lipase family with potential functions in plant secondary metabolism. Subcellular distribution analysis demonstrated that GhGLIP localized to the nucleus, cytoplasm and plasma membrane. GhGLIP was expressed predominantly at 5-15 day post anthesis (dpa) in developing ovules and elongating fibers, measured as mRNA levels and enzyme activity. Ectopic overexpression of GhGLIP in Arabidopsis plants resulted in enhanced seed development, including length and fresh weight. Meanwhile, there was increased soluble sugar and protein storage in transgenic Arabidopsis plants, coupled with the promotion of lipase activity. Moreover, the expression of cotton GhGLIP is induced by ethylene (ETH) treatment in vitro. A 1,954-bp GhGLIP promoter was isolated and expressed high activity in driving green fluorescence protein (GFP) expression in tobacco leaves. Cis-acting element analysis of the GhGLIP promoter (pGhGLIP) indicated the presence of an ethylene-responsive element (ERE), and transgenic tobacco leaves with ectopic expression of pGhGLIP::GFP-GUS showed increased GUS activity after ETH treatment. In summary, these results suggest that GhGLIP is a functional enzyme involved in ovule and fiber development and performs significant roles in seed development.

Two pivotal RNA editing sites in the mitochondrial atp1mRNA are required for ATP synthase to produce sufficient ATP for cotton fiber cell elongation

RNA editing is a post-transcriptional maturation process affecting organelle transcripts in land plants. However, the molecular functions and physiological roles of RNA editing are still poorly understood. Using high-throughput sequencing, we identified 692 RNA editing sites in the Gossypium hirsutum mitochondrial genome. A total of 422 editing sites were found in the coding regions and all the edits are cytidine (C) to uridine (U) conversions. Comparative analysis showed that two editing sites in Ghatp1, C1292 and C1415, had a prominent difference in editing efficiency between fiber and ovule. Biochemical and genetic analyses revealed that the two vital editing sites were important for the interaction between the α and β subunits of ATP synthase, which resulted in ATP accumulation and promoted cell growth in yeast. Ectopic expression of C1292, C1415, or doubly edited Ghatp1 in Arabidopsis caused a significant increase in the number of trichomes in leaves and root length. Our results indicate that editing at C1292 and C1415 sites in Ghatp1 is crucial for ATP synthase to produce sufficient ATP for cotton fiber cell elongation. This work extends our understanding of RNA editing in atp1 and ATP synthesis, and provides insights into the function of mitochondrial edited Atp1 protein in higher plants.

Recent insights into cotton functional genomics: progress and future perspectives

Functional genomics has transformed from futuristic concept to well-established scientific discipline during the last decade. Cotton functional genomics promise to enhance the understanding of fundamental plant biology to systematically exploit genetic resources for the improvement of cotton fibre quality and yield, as well as utilization of genetic information for germplasm improvement. However, determining the cotton gene functions is a much more challenging task, which has not progressed at a rapid pace. This article presents a comprehensive overview of the recent tools and resources available with the major advances in cotton functional genomics to develop elite cotton genotypes. This effort ultimately helps to filter a subset of genes that can be used to assemble a final list of candidate genes that could be employed in future novel cotton breeding programme. We argue that next stage of cotton functional genomics requires the draft genomes refinement, re-sequencing broad diversity panels with the development of high-throughput functional genomics tools and integrating multidisciplinary approaches in upcoming cotton improvement programmes.

Protecting cotton crops under elevated CO2 from waterlogging by managing ethylene

Soil waterlogging and subsequent ethylene release from cotton (Gossypium hirsutum L.) tissues has been linked with abscission of developing cotton fruits. This glasshouse study investigates the effect of a 9-day waterlogging event and CO2 enrichment (eCO2, 700 parts per million (ppm)) on a fully linted cultivar 'Empire' and a lintless cotton mutant (5B). We hypothesised that cotton performance in extreme environments such as waterlogging can be improved through mitigating ethylene action. Plants were grown at 28:20°C day:night temperature, 50-70% relative humidity and a 14:10 light:dark photoperiod under natural light and were exposed to waterlogging and eCO2 at early reproductive growth. Ethylene synthesis was inhibited by spraying aminoethoxyvinylglycine (830ppm) 1 day before waterlogging. Waterlogging significantly increased ethylene release from both cotton genotypes, although fruit production was significantly inhibited only in Empire. Aminoethoxyvinylglycine consistently reduced waterlogging-induced abscission of fruits, mainly in Empire. Limited damage to fruits in 5B, despite increased ethylene production during waterlogging, suggested that fruit abscission in 5B was inhibited by disrupting ethylene metabolism genetically. Elevated CO2 promoted fruit production in both genotypes and was more effective in 5B than in Empire plants. Hence 5B produced more fruits than Empire, providing additional sinks (existing and new fruit) that enhanced the response to CO2 enrichment.

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.

Genome-wide investigation and expression profiling of APX gene family in Gossypium hirsutum provide new insights in redox homeostasis maintenance during different fiber development stages

Ascorbate peroxidase (APX) is a member of heme-containing peroxidases which catalyze the H₂O₂-dependent oxidation of a wide range of substrates in plants and animals. As is known, H₂O₂ acts as a signaling molecule in the regulation of fiber development. Our previous work reported that ascorbate peroxidase 1 (GhAPX1) was important for cotton fiber elongation. However, knowledge about APX gene family members and their evolutionary and functional characteristics in cotton is limited. Here, we report 26 GhAPX genes by genome-wide investigation of tetraploid cotton Gossypium hirsutum. Phylogenetic and gene structure analyses classified these APX members into five clades and syntenic analysis suggested two duplication events. Expression profiling of the 26 APXs revealed that ten members are expressed in cotton fibers. Notably, GhAPX10A, GhAPX10D, GhAPX12A, and GhAPX12D showed high expression levels in 30-day fiber, while GhAPX1A/D, GhAPX3A/D, and GhAPX6A/D showed very low expression levels. The enzyme activity and H₂O₂ content assays revealed that cotton fiber kept high enzyme activity and the lowest H₂O₂ level in 30-day fibers, indicating that other than GhAPX1, the newly reported APX members are responsible for the reactive oxygen species homeostasis in the cotton fiber maturation stages. Expression profiling of ten fiber-expressed APXs after phytohormone treatments revealed their regulation patterns by different stimuli, suggesting that GhAPX1, GhAPX12A, and GhAPX12D are responsible to most phytohormone treatments. Our data provided evolutionary and functional information of GhAPX gene family members and revealed that different members are responsible to redox homeostasis during different cotton fiber development stages.

Genomewide analysis of the lateral organ boundaries domain gene family in Eucalyptus grandis reveals members that differentially impact secondary growth

Lateral Organ Boundaries Domain (LBD) proteins are plant-specific transcription factors playing crucial roles in growth and development. However, the function of LBD proteins in Eucalyptus grandis remains largely unexplored. In this study, LBD genes in E. grandis were identified and characterized using bioinformatics approaches. Gene expression patterns in various tissues and the transcriptional responses of EgLBDs to exogenous hormones were determined by qRT-PCR. Functions of the selected EgLBDs were studied by ectopically overexpressing in a hybrid poplar (Populus alba × Populus glandulosa). Expression levels of genes in the transgenic plants were investigated by RNA-seq. Our results showed that there were forty-six EgLBD members in the E. grandis genome and three EgLBDs displayed xylem- (EgLBD29) or phloem-preferential expression (EgLBD22 and EgLBD37). Confocal microscopy indicated that EgLBD22, EgLBD29 and EgLBD37 were localized to the nucleus. Furthermore, we found that EgLBD22, EgLBD29 and EgLBD37 were responsive to the treatments of indol-3-acetic acid and gibberellic acid. More importantly, we demonstrated EgLBDs exerted different influences on secondary growth. Namely, 35S::EgLBD37 led to significantly increased secondary xylem, 35S::EgLBD29 led to greatly increased phloem fibre production, and 35S::EgLBD22 showed no obvious effects. We revealed that key genes related to gibberellin, ethylene and auxin signalling pathway as well as cell expansion were significantly up- or down-regulated in transgenic plants. Our new findings suggest that LBD genes in E. grandis play important roles in secondary growth. This provides new mechanisms to increase wood or fibre production.