Genome-Wide Identification of R2R3-MYB Transcription Factors Regulating Secondary Cell Wall Thickening in Cotton Fiber Development

Cytochemical localization and synthesis mechanism of the glucomannan in pseudobulbs of Bletilla striata Reichb. f

The dried pseudobulbs of Bletilla striata, an important traditional Chinese medicine named BaiJi, have an extraordinary polysaccharide content and excellent prospects for medicinal effects. However, the distribution and molecular mechanism underlying biosynthesis are poorly understood. In this study, chemical and immunologic analyses were performed in representative tissues of B. striata, and the results showed that what are conventionally termed Bletilla striata polysaccharides (BSPs) are water-soluble polysaccharides deposited only in pseudobulbs. The structural component of BSPs is glucomannan, with a mannose:glucose mass ratio of ~3:2. BSPs are present in the parenchyma of the pseudobulbs in cells known as glucomannan idioblasts and distributed in the cytoplasm within cellular membranes, but are not contained in the vacuole. Comparative transcriptomics and bioinformatics analyses mapped the pathway from sucrose to BSP and identified BsGPI, BsmanA, and BsCSLAs as the key genes of BSP biosynthesis, suggesting that the functional differentiation of the cellulose synthase-like family A (CSLA) may be critical for the flow of glucomannan to the BSP or cell wall. Subsequently, virus-mediated gene silencing showed that silencing of two CSLAs (Bs03G11846 and Bs03G11849) led to a decrease in BSP content, and yeast two-hybrid and luciferase complementation experiments confirmed that four CSLAs (Bs03G11846, Bs03G11847, Bs03G11848, and Bs03G11849) can form homo- or heterodimers, suggesting that multiple CSLAs may form a large complex that functions in BSP synthesis. Our results provide cytological evidence of BSP and describe the isolation and characterization of candidate genes involved in BSP synthesis, laying a solid foundation for further research on its regulation mechanisms and the genetic engineering breeding of B. striata.

An aldehyde dehydrogenase gene, GhALDH7B4_A06, positively regulates fiber strength in upland cotton (Gossypium hirsutum L.)

High fiber strength (FS) premium cotton has significant market demand. Consequently, enhancing FS is a major objective in breeding quality cotton. However, there is a notable lack of known functionally applicable genes that can be targeted for breeding. To address this issue, our study used specific length-amplified fragment sequencing combined with bulk segregant analysis to study FS trait in an F2 population. Subsequently, we integrated these results with previous quantitative trait locus mapping results regarding fiber quality, which used simple sequence repeat markers in F2, F2:3, and recombinant inbred line populations. We identified a stable quantitative trait locus qFSA06 associated with FS located on chromosome A06 (90.74-90.83 Mb). Within this interval, we cloned a gene, GhALDH7B4_A06, which harbored a critical mutation site in coding sequences that is distinct in the two parents of the tested cotton line. In the paternal parent Ji228, the gene is normal and referred to as GhALDH7B4_A06O; however, there is a nonsense mutation in the maternal parent Ji567 that results in premature termination of protein translation, and this gene is designated as truncated GhALDH7B4_A06S. Validation using recombinant inbred lines and gene expression analysis revealed that this mutation site is correlated with cotton FS. Virus-induced gene silencing of GhALDH7B4 in cotton caused significant decreases in FS and fiber micronaire. Conversely, GhALDH7B4_A06O overexpression in Arabidopsis boosted cell wall component contents in the stem. The findings of our study provide a candidate gene for improving cotton fiber quality through molecular breeding.

Fingerprint Finder: Identifying Genomic Fingerprint Sites in Cotton Cohorts for Genetic Analysis and Breeding Advancement

Genomic data in Gossypium provide numerous data resources for the cotton genomics community. However, to fill the gap between genomic analysis and breeding field work, detecting the featured genomic items of a subset cohort is essential for geneticists. We developed FPFinder v1.0 software to identify a subset of the cohort's fingerprint genomic sites. The FPFinder was developed based on the term frequency-inverse document frequency algorithm. With the short-read sequencing of an elite cotton pedigree, we identified 453 pedigree fingerprint genomic sites and found that these pedigree-featured sites had a role in cotton development. In addition, we applied FPFinder to evaluate the geographical bias of fiber-length-related genomic sites from a modern cotton cohort consisting of 410 accessions. Enriching elite sites in cultivars from the Yangtze River region resulted in the longer fiber length of Yangze River-sourced accessions. Apart from characterizing functional sites, we also identified 12,536 region-specific genomic sites. Combining the transcriptome data of multiple tissues and samples under various abiotic stresses, we found that several region-specific sites contributed to environmental adaptation. In this research, FPFinder revealed the role of the cotton pedigree fingerprint and region-specific sites in cotton development and environmental adaptation, respectively. The FPFinder can be applied broadly in other crops and contribute to genetic breeding in the future.

Genome-Wide Analysis of Cotton MYB Transcription Factors and the Functional Validation of GhMYB in Response to Drought Stress

MYB transcription factors play important roles during abiotic stress responses in plants. However, little is known about the accurate systematic analysis of MYB genes in the four cotton species, Gossypium hirsutum, G. barbadense, G. arboreum and G. raimondii. Herein, we performed phylogenetic analysis and showed that cotton MYBs and Arabidopsis MYBs were clustered in the same subfamilies for each species. The identified cotton MYBs were distributed unevenly on chromosomes in various densities for each species, wherein genome-wide tandem and segment duplications were the main driving force of MYB family expansion. Synteny analysis suggested that the abundant collinearity pairs of MYBs were identified between G. hirsutum and the other three species, and that they might have undergone strong purification selection. Characteristics of conserved motifs, along with their consensus sequence, promoter cis elements and gene structure, revealed that MYB proteins might be highly conserved in the same subgroups for each species. Subsequent analysis of differentially expressed genes and expression patterns indicated that most GhMYBs might be involved in response to drought (especially) and salt stress, which was supported by the expression levels of nine GhMYBs using real-time quantitative PCR. Finally, we performed a workflow that combined virus-induced gene silencing and the heterologous transformation of Arabidopsis, which confirmed the positive roles of GhMYBs under drought conditions, as validated by determining the drought-tolerant phenotypes, damage index and/or water loss rate. Collectively, our findings not only expand our understanding of the relationships between evolution and function of MYB genes, but they also provide candidate genes for cotton breeding.

The transcription factor ERF108 interacts with AUXIN RESPONSE FACTORs to mediate cotton fiber secondary cell wall biosynthesis

Phytohormones play indispensable roles in plant growth and development. However, the molecular mechanisms underlying phytohormone-mediated regulation of fiber secondary cell wall (SCW) formation in cotton (Gossypium hirsutum) remain largely underexplored. Here, we provide mechanistic evidence for functional interplay between the APETALA2/ethylene response factor (AP2/ERF) transcription factor GhERF108 and auxin response factors GhARF7-1 and GhARF7-2 in dictating the ethylene-auxin signaling crosstalk that regulates fiber SCW biosynthesis. Specifically, in vitro cotton ovule culture revealed that ethylene and auxin promote fiber SCW deposition. GhERF108 RNA interference (RNAi) cotton displayed remarkably reduced cell wall thickness compared with controls. GhERF108 interacted with GhARF7-1 and GhARF7-2 to enhance the activation of the MYB transcription factor gene GhMYBL1 (MYB domain-like protein 1) in fibers. GhARF7-1 and GhARF7-2 respond to auxin signals that promote fiber SCW thickening. GhMYBL1 RNAi and GhARF7-1 and GhARF7-2 virus-induced gene silencing (VIGS) cotton displayed similar defects in fiber SCW formation as GhERF108 RNAi cotton. Moreover, the ethylene and auxin responses were reduced in GhMYBL1 RNAi plants. GhMYBL1 directly binds to the promoters of GhCesA4-1, GhCesA4-2, and GhCesA8-1 and activates their expression to promote cellulose biosynthesis, thereby boosting fiber SCW formation. Collectively, our findings demonstrate that the collaboration between GhERF108 and GhARF7-1 or GhARF7-2 establishes ethylene-auxin signaling crosstalk to activate GhMYBL1, ultimately leading to the activation of fiber SCW biosynthesis.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Advances in understanding epigenetic regulation of plant trichome development: a comprehensive review

Plant growth and development are controlled by a complex gene regulatory network, which is currently a focal point of research. It has been established that epigenetic factors play a crucial role in plant growth. Trichomes, specialized appendages that arise from epidermal cells, are of great significance in plant growth and development. As a model system for studying plant development, trichomes possess both commercial and research value. Epigenetic regulation has only recently been implicated in the development of trichomes in a limited number of studies, and microRNA-mediated post-transcriptional regulation appears to dominate in this context. In light of this, we have conducted a review that explores the interplay between epigenetic regulations and the formation of plant trichomes, building upon existing knowledge of hormones and transcription factors in trichome development. Through this review, we aim to deepen our understanding of the regulatory mechanisms underlying trichome formation and shed light on future avenues of research in the field of epigenetics as it pertains to epidermal hair growth.

Identification of MYB gene family and functional analysis of GhMYB4 in cotton (Gossypium spp.)

Myeloblastosis (MYB) transcription factors (TFs) form a large gene family involved in a variety of biological processes in plants. Little is known about their roles in the development of cotton pigment glands. In this study, 646 MYB members were identified in Gossypium hirsutum genome and phylogenetic classification was analyzed. Evolution analysis revealed assymetric evolution of GhMYBs during polyploidization and sequence divergence of MYBs in G. hirustum was preferentially happend in D sub-genome. WGCNA (weighted gene co-expression network analysis) showed that four modules had potential relationship with gland development or gossypol biosynthesis in cotton. Eight differentially expressed GhMYB genes were identified by screening transcriptome data of three pairs of glanded and glandless cotton lines. Of these, four were selected as candidate genes for cotton pigment gland formation or gossypol biosynthesis by qRT-PCR assay. Silencing of GH_A11G1361 (GhMYB4) downregulated expression of multiple genes in gossypol biosynthesis pathway, indicating it could be involved in gossypol biosynthesis. The potential protein interaction network suggests that several MYBs may have indirect interaction with GhMYC2-like, a key regulator of pigment gland formation. Our study was the systematic analysis of MYB genes in cotton pigment gland development, providing candidate genes for further study on the roles of cotton MYB genes in pigment gland formation, gossypol biosynthesis and future crop plant improvement.

Genome-wide analysis elucidates the roles of GhHMA genes in different abiotic stresses and fiber development in upland cotton

The heavy metal-binding domain is involved in heavy metal transporting and plays a significant role in plant detoxification. However, the functions of HMAs are less well known in cotton. In this study, a total of 143 GhHMAs (heavy metal-binding domain) were detected by genome-wide identification in G. hirsutum L. All the GhHMAs were classified into four groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GhHMA genes was highly conserved. 212 paralogous GhHMA gene pairs were identified, and the segmental duplications were the main role to the expansion of GhHMAs. The Ka/Ks values suggested that the GhHMA gene family has undergone purifying selection during the long-term evolutionary process. GhHMA3 and GhHMA75 were located in the plasma membrane, while GhHMA26, GhHMA117 and GhHMA121 were located in the nucleus, respectively. Transcriptomic data and qRT-PCR showed that GhHMA26 exhibited different expression patterns in each tissue and during fiber development or under different abiotic stresses. Overexpressing GhHMA26 significantly promoted the elongation of leaf trichomes and also improved the tolerance to salt stress. Therefore, GhHMA26 may positively regulate fiber elongation and abiotic stress. Yeast two-hybrid assays indicated that GhHMA26 and GhHMA75 participated in multiple biological functions. Our results suggest some genes in the GhHMAs might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GhHMA genes in cotton.

Systematic Analysis and Functional Characterization of R2R3-MYB Genes in Scutellaria baicalensis Georgi

R2R3-MYB transcription factors participate in multiple critical biological processes, particularly as relates to the regulation of secondary metabolites. The dried root of Scutellaria baicalensis Georgi is a traditional Chinese medicine and possesses various bioactive attributes including anti-inflammation, anti-HIV, and anti-COVID-19 properties due to its flavonoids. In the current study, a total of 95 R2R3-MYB genes were identified in S. baicalensis and classified into 34 subgroups, as supported by similar exon–intron structures and conserved motifs. Among them, 93 R2R3-SbMYBs were mapped onto nine chromosomes. Collinear analysis revealed that segmental duplications were primarily responsible for driving the evolution and expansion of the R2R3-SbMYB gene family. Synteny analyses showed that the ortholog numbers of the R2R3-MYB genes between S. baicalensis and other dicotyledons had a higher proportion compared to that which is found from the monocotyledons. RNA-seq data indicated that the expression patterns of R2R3-SbMYBs in different tissues were different. Quantitative reverse transcriptase-PCR (qRT-PCR) analysis showed that 36 R2R3-SbMYBs from different subgroups exhibited specific expression profiles under various conditions, including hormone stimuli treatments (methyl jasmonate and abscisic acid) and abiotic stresses (drought and cold shock treatments). Further investigation revealed that SbMYB18/32/46/60/70/74 localized in the nucleus, and SbMYB18/32/60/70 possessed transcriptional activation activity, implying their potential roles in the regulatory mechanisms of various biological processes. This study provides a comprehensive understanding of the R2R3-SbMYBs gene family and lays the foundation for further investigation of their biological function.

Molecular Mechanisms of Plant Trichome Development

Plant trichomes, protrusions formed from specialized aboveground epidermal cells, provide protection against various biotic and abiotic stresses. Trichomes can be unicellular, bicellular or multicellular, with multiple branches or no branches at all. Unicellular trichomes are generally not secretory, whereas multicellular trichomes include both secretory and non-secretory hairs. The secretory trichomes release secondary metabolites such as artemisinin, which is valuable as an antimalarial agent. Cotton trichomes, also known as cotton fibers, are an important natural product for the textile industry. In recent years, much progress has been made in unraveling the molecular mechanisms of trichome formation in Arabidopsis thaliana, Gossypium hirsutum, Oryza sativa, Cucumis sativus, Solanum lycopersicum, Nicotiana tabacum, and Artemisia annua. Here, we review current knowledge of the molecular mechanisms underlying fate determination and initiation, elongation, and maturation of unicellular, bicellular and multicellular trichomes in several representative plants. We emphasize the regulatory roles of plant hormones, transcription factors, the cell cycle and epigenetic modifications in different stages of trichome development. Finally, we identify the obstacles and key points for future research on plant trichome development, and speculated the development relationship between the salt glands of halophytes and the trichomes of non-halophytes, which provides a reference for future studying the development of plant epidermal cells.

Transcriptome Time-Course Analysis in the Whole Period of Cotton Fiber Development

Gossypium hirsutum and Gossypium barbadense are the widely cultivated tetraploid cottons around the world, which evolved great differences in the fiber yield and quality due to the independent domestication process. To reveal the genetic basis of the difference, we integrated 90 samples from ten time points during the fiber developmental period for investigating the dynamics of gene expression changes associated with fiber in G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124 and acc. 3-79. Globally, 44,484 genes expressed in all three cultivars account for 61.14% of the total genes. About 61.39% (N = 3,412) of the cotton transcription factors were involved in fiber development, which consisted of 58 cotton TF families. The differential analysis of intra- and interspecies showed that 3 DPA had more expression changes. To discover the genes with temporally changed expression profiles during the whole fiber development, 1,850 genes predominantly expressed in G. hirsutum and 1,050 in G. barbadense were identified, respectively. Based on the weighted gene co-expression network and time-course analysis, several candidate genes, mainly involved in the secondary cell wall synthesis and phytohormones, were identified in this study, underlying possibly the transcriptional regulation and molecular mechanisms of the fiber quality differences between G. barbadense and G. hirsutum. The quantitative real-time PCR validation of the candidate genes was consistent with the RNA-seq data. Our study provides a strong rationale for the analysis of gene function and breeding of high-quality cotton.

Two R2R3-MYB genes cooperatively control trichome development and cuticular wax biosynthesis in Prunus persica

The fruit surface has an enormous impact on the external appearance and postharvest shelf-life of fruit. Here, we report two functionally redundant genes, PpMYB25 and PpMYB26, involved in regulation of fruit skin texture in peach. PpMYB25 can activate transcription of PpMYB26 and they both induce trichome development and cuticular wax accumulation, resulting in peach fruit with a fuzzy and dull appearance. By contrast, nonfunctional mutation of PpMYB25 caused by an insertional retrotransposon in the last exon in nectarine fails to activate transcription of PpMYB26, resulting in nectarine fruit with a smooth and shiny appearance due to loss of trichome initiation and decreased cuticular wax accumulation. Secondary cell wall biosynthesis in peach fruit pubescence is controlled by a transcriptional regulatory network, including the master regulator PpNAC43 and its downstream MYB transcription factors such as PpMYB42, PpMYB46 and PpMYB83. Our results show that PpMYB25 and PpMYB26 coordinately regulate fruit pubescence and cuticular wax accumulation and their simultaneous perturbation results in the origin of nectarine, which is botanically classified as a subspecies of peach.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

Genome-Wide Identification of Cotton (Gossypium spp.) Glycerol-3-Phosphate Dehydrogenase (GPDH) Family Members and the Role of GhGPDH5 in Response to Drought Stress

Glycerol-3-phosphate dehydrogenase (GPDH) is a key enzyme in plant glycerol synthesis and metabolism, and plays an important role in plant resistance to abiotic stress. Here, we identified 6, 7, 14 and 14 GPDH genes derived from Gossypium arboreum, Gossypium raimondii, Gossypium barbadense and Gossypium hirsutum, respectively. Phylogenetic analysis assigned these genes into three classes, and most of the genes within the family were expanded by whole-genome duplication (WGD) and segmental duplications. Moreover, determination of the nonsynonymous substitution rate/synonymous substitution rate (Ka/Ks) ratio showed that the GPDH had an evolutionary preference for purifying selection. Transcriptome data revealed that GPDH genes were more active in the early stages of fiber development. Additionally, numerous stress-related cis-elements were identified in the potential promoter region. Then, a protein-protein-interaction (PPI) network of GPDH5 in G. hirsutum was constructed. In addition, we predicted 30 underlying miRNAs in G. hirsutum. Functional validation results indicated that silencing GhGPDH5 diminished drought tolerance in the upland cotton TM-1 line. In summary, this study provides a fundamental understanding of the GPDH gene family in cotton, GhGPDH5 exerts a positive effect during drought stress and is potentially involved in stomatal closure movements.

Comparative transcriptome analysis reveals the regulation network for fiber strength in cotton

OBJECTIVE

Determine the effect of secondary cell wall (SCW) thickness and microcrystalline cellulose content (MCC) on mature fiber strength (FS) and reveal through comparative transcriptome analysis the molecular regulation network governing FS in cotton.

RESULTS

Transmission electron microscope (TEM) analysis of two parent varieties, Prema with elite FS and 86-1 with weak fiber, revealed significant difference in the SCW but not in MCC. Transcriptome analysis revealed that genes differentially expressed during SCW thickening (20 DPA) are highly related to FS; in particular, up-regulated genes such as UDPG, CESA2, and NAC83 were important in SCW thickening, likely contributing to higher FS. GO and KEGG enrichment analysis revealed the common up-regulated genes to be enriched in carbon metabolism and terms relating to the cell wall.

CONCLUSIONS

We developed two recombinant inbred lines with elite FS, selected from the filial generation of Prema and 86-1. By comparing transcriptomic data, we revealed the gene expression network governing SCW thickness in mature fiber. Our results provide solid insights into the relationship of the SCW and FS.

Genome-Wide Identification of the HD-ZIP III Subfamily in Upland Cotton Reveals the Involvement of GhHB8-5D in the Biosynthesis of Secondary Wall in Fiber and Drought Resistance

A subfamily of transcription factors known as HD-ZIP III plays distinct roles in the secondary cell wall biosynthesis, which could be attributed to the quality of cotton fiber and adaptation to drought stress. In this study, 18 HD-ZIP III genes were identified as genome wide from the upland cotton (Gossypium hirsutum). These genes are distributed on 14 different chromosomes, and all of them have undergone segmental duplications. Numerous cis-elements were identified in the promoter regions, which are related to phytohormone responses and abiotic stresses. Expression profiling of these genes by quantitative real-time (qRT)-PCR illustrated their differential spatial expression, with preferential expression in cotton fiber. Among these genes, GhHB8-5D was predicted to encode a protein that is targeted to the cell nucleus and having self-activation ability. In addition, the ectopic expression of GhHB8-5D or its synonymous mutant GhHB8-5Dm in Arabidopsis resulted in stunted plant growth, curly leaves, and twisted inflorescence stems. Microscopy examination revealed that the morphology of vascular bundles and deposition of secondary wall had substantially altered in stems, which is concomitant with the significant alteration in the transcription levels of secondary wall-related genes in these transgenic Arabidopsis. Further, ectopic expression of GhHB8-5D or GhHB8-5Dm in Arabidopsis also led to significant increase in green seedling rate and reduction in root length relative to wild type when the plants were grown under mimicked drought stress conditions. Taken together, our results may shed new light on the functional roles of GhHB8-5D that is attributable for secondary cell wall thickening in response to drought stress. Such a finding may facilitate a novel strategy for improving plant adaptations to environmental changes via regulating the biosynthesis of secondary cell wall.

Genome-Wide Identification of R2R3-MYB Transcription Factors: Discovery of a "Dual-Function" Regulator of Gypenoside and Flavonol Biosynthesis in Gynostemma pentaphyllum

The R2R3-MYB gene family participates in several plant physiological processes, especially the regulation of the biosynthesis of secondary metabolites. However, little is known about the functions of R2R3-MYB genes in Gynostemma pentaphyllum (G. pentaphyllum), a traditional Chinese medicinal herb that is an excellent source of gypenosides (a class of triterpenoid saponins) and flavonoids. In this study, a systematic genome-wide analysis of the R2R3-MYB gene family was performed using the recently sequenced G. pentaphyllum genome. In total, 87 R2R3-GpMYB genes were identified and subsequently divided into 32 subgroups based on phylogenetic analysis. The analysis was based on conserved exon-intron structures and motif compositions within the same subgroup. Collinearity analysis demonstrated that segmental duplication events were majorly responsible for the expansion of the R2R3-GpMYB gene family, and Ka/Ks analysis indicated that the majority of the duplicated R2R3-GpMYB genes underwent purifying selection. A combination of transcriptome analysis and quantitative reverse transcriptase-PCR (qRT-PCR) confirmed that Gynostemma pentaphyllum myeloblastosis 81 (GpMYB81) along with genes encoding gypenoside and flavonol biosynthetic enzymes exhibited similar expression patterns in different tissues and responses to methyl jasmonate (MeJA). Moreover, GpMYB81 could bind to the promoters of Gynostemma pentaphyllum farnesyl pyrophosphate synthase 1 (GpFPS1) and Gynostemma pentaphyllum chalcone synthase (GpCHS), the key structural genes of gypenoside and flavonol biosynthesis, respectively, and activate their expression. Altogether, this study highlights a novel transcriptional regulatory mechanism that suggests that GpMYB81 acts as a "dual-function" regulator of gypenoside and flavonol biosynthesis in G. pentaphyllum.

Novel insights into water-deficit-responsive mRNAs and lncRNAs during fiber development in Gossypium hirsutum

BACKGROUND

The fiber yield and quality of cotton are greatly and periodically affected by water deficit. However, the molecular mechanism of the water deficit response in cotton fiber cells has not been fully elucidated.

RESULTS

In this study, water deficit caused a significant reduction in fiber length, strength, and elongation rate but a dramatic increase in micronaire value. To explore genome-wide transcriptional changes, fibers from cotton plants subjected to water deficit (WD) and normal irrigation (NI) during fiber development were analyzed by transcriptome sequencing. Analysis showed that 3427 mRNAs and 1021 long noncoding RNAs (lncRNAs) from fibers were differentially expressed between WD and NI plants. The maximum number of differentially expressed genes (DEGs) and lncRNAs (DERs) was identified in fibers at the secondary cell wall biosynthesis stage, suggesting that this is a critical period in response to water deficit. Twelve genes in cotton fiber were differentially and persistently expressed at ≥ five time points, suggesting that these genes are involved in both fiber development and the water-deficit response and could potentially be used in breeding to improve cotton resistance to drought stress. A total of 540 DEGs were predicted to be potentially regulated by DERs by analysis of coexpression and genomic colocation, accounting for approximately 15.76% of all DEGs. Four DERs, potentially acting as target mimics for microRNAs (miRNAs), indirectly regulated their corresponding DEGs in response to water deficit.

CONCLUSIONS

This work provides a comprehensive transcriptome analysis of fiber cells and a set of protein-coding genes and lncRNAs implicated in the cotton response to water deficit, significantly affecting fiber quality during the fiber development stage.

Comparative transcriptome analysis of fiber and nonfiber tissues to identify the genes preferentially expressed in fiber development in Gossypium hirsutum

Cotton is an important natural fiber crop and economic crop worldwide. The quality of cotton fiber directly determines the quality of cotton textiles. Identifying cotton fiber development-related genes and exploring their biological functions will not only help to better understand the elongation and development mechanisms of cotton fibers but also provide a theoretical basis for the cultivation of new cotton varieties with excellent fiber quality. In this study, RNA sequencing technology was used to construct transcriptome databases for different nonfiber tissues (root, leaf, anther and stigma) and fiber developmental stages (7 days post-anthesis (DPA), 14 DPA, and 26 DPA) of upland cotton Coker 312. The sizes of the seven transcriptome databases constructed ranged from 4.43 to 5.20 Gb, corresponding to approximately twice the genome size of Gossypium hirsutum (2.5 Gb). Among the obtained clean reads, 83.32% to 88.22% could be compared to the upland cotton TM-1 reference genome. By analyzing the differential gene expression profiles of the transcriptome libraries of fiber and nonfiber tissues, we obtained 1205, 1135 and 937 genes with significantly upregulated expression at 7 DPA, 14 DPA and 26 DPA, respectively, and 124, 179 and 213 genes with significantly downregulated expression. Subsequently, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analyses were performed, which revealed that these genes were mainly involved in catalytic activity, carbohydrate metabolism, the cell membrane and organelles, signal transduction and other functions and metabolic pathways. Through gene annotation analysis, many transcription factors and genes related to fiber development were screened. Thirty-six genes were randomly selected from the significantly upregulated genes in fiber, and expression profile analysis was performed using qRT-PCR. The results were highly consistent with the gene expression profile analyzed by RNA-seq, and all of the genes were specifically or predominantly expressed in fiber. Therefore, our RNA sequencing-based comparative transcriptome analysis will lay a foundation for future research to provide new genetic resources for the genetic engineering of improved cotton fiber quality and for cultivating new transgenic cotton germplasms for fiber quality improvement.

The Basic Helix-Loop-Helix Transcription Factor SmbHLH1 Represses Anthocyanin Biosynthesis in Eggplant

Many basic helix-loop-helix transcription factors (TFs) have been reported to promote anthocyanin biosynthesis in numerous plant species, but little is known about bHLH TFs that inhibit anthocyanin accumulation. In this study, SmbHLH1 from Solanum melongena was identified as a negative regulator of anthocyanin biosynthesis. However, SmbHLH1 showed high identity with SmTT8, which acts as a SmMYB113-dependent positive regulator of anthocyanin-biosynthesis in plants. Overexpression of SmbHLH1 in eggplant caused a dramatic decrease in anthocyanin accumulation. Only the amino acid sequences at the N and C termini of SmbHLH1 differed from the SmTT8 sequence. Expression analysis revealed that the expression pattern of SmbHLH1 was opposite to that of anthocyanin accumulation. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that SmbHLH1 could not interact with SmMYB113. Dual-luciferase assay demonstrated that SmbHLH1 directly repressed the expression of SmDFR and SmANS. Our results demonstrate that the biological function of bHLHs in anthocyanin biosynthesis may have evolved and provide new insight into the molecular functions of orthologous genes from different plant species.

GhMYB7 promotes secondary wall cellulose deposition in cotton fibres by regulating GhCesA gene expression through three distinct cis-elements

Cotton fibre is the most important source for natural textiles. The secondary cell walls (SCWs) of mature cotton fibres contain the highest proportion of cellulose content (> 90%) in any plant. The onset and progression of SCW cellulose synthesis need to be tightly controlled to balance fibre elongation and cell wall deposition. However, regulatory mechanisms that control cellulose synthesis during cotton fibre growth remain elusive. Here, we conducted genetic and functional analyses demonstrating that the R2R3-MYB GhMYB7 controls cotton fibre cellulose synthesis. Overexpression of GhMYB7 in cotton sped up SCW cellulose biosynthesis in fibre cells, and led to shorter fibres with thicker walls. By contrast, RNA interference (RNAi) silencing of GhMYB7 delayed fibre SCW cellulose synthesis and resulted in elongated fibres with thinner walls. Furthermore, we demonstrated that GhMYB7 regulated cotton fibre SCW cellulose synthases by directly binding to three distinct cis-elements in the respective GhCesA4, GhCesA7 and GhCesA8 promoters. We found that this regulatory mechanism of cellulose synthesis was 'hi-jacked' also by other GhMYBs. Together, our findings uncover a hitherto-unknown mechanism that cotton fibre employs to regulate SCW cellulose synthesis. Our results also provide a strategy for genetic improvement of SCW thickness of cotton fibre.

GhUBX controlling helical growth results in production of stronger cotton fiber

Cotton fiber is an excellent model for studying plant cell elongation and cell wall biogenesis as well because they are highly polarized and use conserved polarized diffuse growth mechanism. Fiber strength is an important trait among cotton fiber qualities due to ongoing changes in spinning technology. However, the molecular mechanism of fiber strength forming is obscure. Through map-based cloning, we identified the fiber strength gene GhUBX. Increasing its expression, the fiber strength of the transgenic cotton was significantly enhanced compared to the receptor W0 and the helices number of the transgenic fiber was remarkably increased. Additionally, we proved that GhUBX regulates the fiber helical growth by degrading the GhSPL1 via the ubiquitin 26S–proteasome pathway. Taken together, we revealed the internal relationship between fiber helices and fiber stronger. It will be useful for improving the fiber quality in cotton breeding and illustrating the molecular mechanism for plant twisted growth.

•

Isolation of the first fiber strength gene GhUBX using map-based cloning strategy

•

Verification of the function of GhUBX experimentally in transgenic cotton

•

Link helices to the cotton fiber strength, that more helices make fiber stronger

•

An ubiquitin–proteasome system regulating the development of cotton fiber

Alternative splicing of MaMYB16L regulates starch degradation in banana fruit during ripening

The alternative splicing of select genes is an important mechanism to regulate responses to endogenous and environmental signals in plants. However, the role of alternative splicing in regulating fruit ripening remains unclear. Here, we discovered that MaMYB16L, an R1-type MYB transcription factor, undergoes alternative splicing and generates two transcripts, the full-length isoform MaMYB16L and a truncated form MaMYB16S, in banana fruit. During banana fruit ripening, the alternative splicing process intensifies with downregulated MaMYB16L and upregulated MaMYB16S. Moreover, MaMYB16L is a transcriptional repressor that directly binds with the promoters of many genes associated with starch degradation and MaDREB2, a positive ripening regulator, and represses their expression. In contrast, MaMBY16S lacks a DNA-binding domain but competitively combines and forms non-functional heterodimers with functional MaMYB16L. MaMYB16L-MaMYB16S heterodimers decrease the binding capacity and transrepression activity of MaMYB16L. The downregulation of MaMYB16L and the upregulation of MaMYB16S, that is, a decreased ratio of active to non-active isoforms, facilitates the activation of ripening-related genes and thereby promotes fruit ripening. Furthermore, the transient overexpression of MaMYB16S promotes banana fruit ripening, whereas the overexpression of MaMYB16L delays this process. Therefore, the alternative splicing of MaMYB16L might generate a self-controlled regulatory loop to regulate banana fruit ripening.

Genome-Wide Analysis of the R2R3-MYB Gene Family in Fragaria × ananassa and Its Function Identification During Anthocyanins Biosynthesis in Pink-Flowered Strawberry

The strawberry (Fragaria × ananassa) is an economically important fruit throughout the world. The large R2R3-MYB gene family participates in a variety of plant functions, including anthocyanin biosynthesis. The present study is the first genome-wide analysis of the MYB gene family in the octoploid strawberry and describes the identification and characterization of the family members using the recently sequenced F. × ananassa genome. Specifically, we aimed to identify the key MYBs involved in petal coloration in the pink-flowered strawberry, which increases its ornamental value. A comprehensive, genome-wide analysis of F. × ananassa R2R3-FaMYBs was performed, investigating gene structures, phylogenic relationships, promoter regions, chromosomal locations, and collinearity. A total of 393 R2R3-FaMYB genes were identified in the F. × ananassa genome and divided into 36 subgroups based on phylogenetic analysis. Most genes with similar functions in the same subgroup exhibited similar exon-intron structures and motif compositions. These R2R3-FaMYBs were unevenly distributed over 28 chromosomes. The expansion of the R2R3-FaMYB gene family in the F. × ananassa genome was found to be caused mainly by segmental duplication. The Ka/Ks analysis indicated that duplicated R2R3-FaMYBs mostly experienced purifying selection and showed limited functional divergence after the duplication events. To elucidate which R2R3-FaMYB genes were associated with anthocyanin biosynthesis in the petals of the pink-flowered strawberry, we compared transcriptional changes in different flower developmental stages using RNA-seq. There were 131 differentially expressed R2R3-FaMYB genes identified in the petals, of which three genes, FaMYB28, FaMYB54, and FaMYB576, appeared likely, based on the phylogenetic analysis, to regulate anthocyanin biosynthesis. The qRT-PCR showed that these three genes were more highly expressed in petals than in other tissues (fruit, leaf, petiole and stolon) and their expressions were higher in red compared to pink and white petals. These results facilitate the clarification on the roles of the R2R3-FaMYB genes in petal coloration in the pink-flowered strawberry. This work provides useful information for further functional analysis on the R2R3-FaMYB gene family in F. × ananassa.

Genome-wide association analysis reveals loci and candidate genes involved in fiber quality traits in sea island cotton (Gossypium barbadense)

BACKGROUND

Sea island cotton (Gossypium barbadense) has markedly superior high quality fibers, which plays an important role in the textile industry and acts as a donor for upland cotton (G. hirsutum) fiber quality improvement. The genetic characteristics analysis and the identification of key genes will be helpful to understand the mechanism of fiber development and breeding utilization in sea island cotton.

RESULTS

In this study, 279 sea island cotton accessions were collected from different origins for genotyping and phenotyping fiber quality traits. A set of 6303 high quality single nucleotide polymorphisms (SNPs) were obtained by high-density CottonSNP80K array. The population characteristics showed that the sea island cotton accessions had wide genetic diversity and were clustered into three groups, with Group1 closely related to Menoufi, an original sea island cotton landrace, and Group2 and Group3 related to widely introduced accessions from Egypt, USA and Former Soviet Union. Further, we used 249 accessions and evaluated five fiber quality traits under normal and salt environments over 2 years. Except for fiber uniformity (FU), fiber length (FL) and fiber elongation (FE) were significantly decreased in salt conditions, while fiber strength (FS) and fiber micronaire (MIC) were increased. Based on 6303 SNPs and genome-wide association study (GWAS) analysis, a total of 34 stable quantitative trait loci (QTLs) were identified for the five fiber quality traits with 25 detected simultaneously under normal and salt environments. Gene Ontology (GO) analysis indicated that candidate genes in the 25 overlapped QTLs were enriched mostly in "cellular and biological process". In addition, "xylem development" and "response to hormone" pathways were also found. Haplotype analyses found that GB_A03G0335 encoding an E3 ubiquitin-protein ligase in QTL TM6004 had SNP variation (A/C) in gene region, was significantly correlated with FL, FS, FU, and FE, implying a crucial role in fiber quality.

CONCLUSIONS

The present study provides a foundation for genetic diversity of sea island cotton accessions and will contribute to fiber quality improvement in breeding practice.