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

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

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

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

Co-localization and analysis of miR477b with fiber length quantitative trait loci in cotton

Cotton is an important cash crop for the textile industry. However, the understanding of natural genetic variation of fiber elongation in relation to miRNA is lacking. A miRNA gene (miR477b) was found to co-localize with a previously mapped fiber length (FL) quantitative trait locus (QTL). The miR477b was differentially expressed during fiber elongation between two backcross inbred lines (BILs) differing in FL and its precursor sequences. Bioinformatics and qRT-PCR analysis were further used to analyse the miRNA genes, which could produce mature miR477b. Cotton plants with virus-induced gene silencing (VIGS) constructs to over-express the allele of miR477b from the BIL with longer fibers had significantly longer fibers as compared with negative control plants, while the VIGS plants with suppressed miRNA expression had significantly shorter fibers. The expression level of the target gene (DELLA) and related genes (RDL1 and EXPA1 for DELLA through HOX3 protein) in the two BILs and/or the VIGS plants were generally congruent, as expected. This report represents one of the first comprehensive studies to integrate QTL linkage mapping and physical mapping of small RNAs with both small and mRNA transcriptome analysis, followed by VIGS, to identify candidate small RNA genes affecting the natural variation of fiber elongation in cotton.

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.

CsREV-CsTCP4-CsVND7 module shapes xylem patterns differentially between stem and leaf to enhance tea plant tolerance to drought

Cultivating drought-tolerant tea varieties enhances both yield and quality of tea plants in northern China. However, the mechanisms underlying their drought tolerance remain largely unknown. Here we identified a key regulator called CsREV, which differentially regulates xylem patterns between leaves and stems, thereby conferring drought tolerance in tea plants. When drought occurs, upregulation of CsREV activates the CsVND7a-dependent xylem vessel differentiation. However, when drought persists, the vessel differentiation is hindered as CsVND7a is downregulated by CsTCP4a. This, combined with the CsREV-promoted secondary-cell-wall thickness of xylem vessel, leads to the enhanced curling of leaves, a characteristic closely associated with plant drought tolerance. Notably, this inhibitory effect of CsTCP4a on CsVND7a expression is absent in stems, allowing stem xylem vessels to continuously differentiate. Overall, the CsREV-CsTCP4-CsVND7 module is differentially utilized to shape the xylem patterns in leaves and stems, potentially balancing water transportation and utilization to improve tea plant drought tolerance.

Comparative Transcriptomic Analysis of Gossypium hirsutum Fiber Development in Mutant Materials (xin w 139) Provides New Insights into Cotton Fiber Development

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Gibberellic acid promotes single-celled fiber elongation through the activation of two signaling cascades in cotton

The agricultural green revolution spectacularly enhanced crop yield through modification of gibberellin (GA) signaling. However, in cotton, the GA signaling cascades remain elusive, limiting our potential to cultivate new cotton varieties and improve yield and quality. Here, we identified that GA prominently stimulated fiber elongation through the degradation of DELLA protein GhSLR1, thereby disabling GhSLR1's physical interaction with two transcription factors, GhZFP8 and GhBLH1. Subsequently, the resultant free GhBLH1 binds to GhKCS12 promoter and activates its expression to enhance VLCFAs biosynthesis. With a similar mechanism, the free GhZFP8 binds to GhSDCP1 promoter and activates its expression. As a result, GhSDCP1 upregulates the expression of GhPIF3 gene associated with plant cell elongation. Ultimately, the two parallel signaling cascades synergistically promote cotton fiber elongation. Our findings outline the mechanistic framework that translates the GA signal into fiber cell elongation, thereby offering a roadmap to improve cotton fiber quality and yield.

Graphene enhances artemisinin production in the traditional medicinal plant Artemisia annua via dynamic physiological processes and miRNA regulation

We investigated the effects of graphene on the model herb Artemisia annua, which is renowned for producing artemisinin, a widely used pharmacological compound. Seedling growth and biomass were promoted when A. annua was cultivated with low concentrations of graphene, an effect which was attributed to a 1.4-fold increase in nitrogen uptake, a 15%-22% increase in chlorophyll fluorescence, and greater abundance of carbon cycling-related bacteria. Exposure to 10 or 20 mg/L graphene resulted in a ∼60% increase in H2O2, and graphene could act as a catalyst accelerator, leading to a 9-fold increase in catalase (CAT) activity in vitro and thereby maintaining reactive oxygen species (ROS) homeostasis. Importantly, graphene exposure led to an 80% increase in the density of glandular secreting trichomes (GSTs), in which artemisinin is biosynthesized and stored. This contributed to a 5% increase in artemisinin content in mature leaves. Interestingly, expression of miR828 was reduced by both graphene and H2O2 treatments, resulting in induction of its target gene AaMYB17, a positive regulator of GST initiation. Subsequent molecular and genetic assays showed that graphene-induced H2O2 inhibits micro-RNA (miRNA) biogenesis through Dicers and regulates the miR828-AaMYB17 module, thus affecting GST density. Our results suggest that graphene may contribute to yield improvement in A. annua via dynamic physiological processes together with miRNA regulation, and it may thus represent a new cultivation strategy for increasing yield capacity through nanobiotechnology.

Improvement of qualitative and quantitative traits in cotton under normal and stressed environments using genomics and biotechnological tools: A review

Due to the increasing demand for high-quality and high fiber-yielding cotton (Gossypium spp.), research into the development of stress-resilient cotton cultivars has acquired greater significance. Various biotic and abiotic stressors greatly affect cotton production and productivity, posing challenges to the future of the textile industry. Moreover, the content and quality of cottonseed oil can also potentially be influenced by future environmental conditions. Apart from conventional methods, genetic engineering has emerged as a potential tool to improve cotton fiber quality and productivity. Identification and modification of genome sequences and the expression levels of yield-related genes using genetic engineering approaches have enabled to increase both the quality and yields of cotton fiber and cottonseed oil. Herein, we evaluate the significance and molecular mechanisms associated with the regulation of cotton agronomic traits under both normal and stressful environmental conditions. In addition, the importance of gossypol, a toxic phenolic compound in cottonseed that can limit consumption by animals and humans, is reviewed and discussed.

GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, play a regulatory role in cotton fiber secondary cell wall biosynthesis

KEY MESSAGE

GhHB14_D10 and GhREV_D5 regulated secondary cell wall formation and played an important role in fiber development. Cotton serves as an important source of natural fiber, and the biosynthesis of the secondary cell wall plays a pivotal role in determining cotton fiber quality. Nevertheless, the intricacies of this mechanism in cotton fiber remain insufficiently elucidated. This study investigates the functional roles of GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, in secondary cell wall biosynthesis in cotton fibers. Both GhHB14_D10 and GhREV_D5 were found to be localized in the nucleus with transcriptional activation activity. Ectopic overexpression of GhHB14_D10 and GhREV_D5 in Arabidopsis resulted in changed xylem differentiation, secondary cell wall deposition, and expression of genes related to the secondary cell wall. Silencing of GhHB14_D10 and GhREV_D5 in cotton led to enhanced fiber length, reduced cell wall thickness, cellulose contents and expression of secondary cell wall-related genes. Moreover, GhHB14_D10's direct interaction with GhREV_D5, and transcriptional regulation of cellulose biosynthesis genes GhCesA4-4 and GhCesA7-2 revealed their collaborative roles in secondary cell wall during cotton fiber development. Overall, these results shed light on the roles of GhHB14_D10 and GhREV_D5 in secondary cell wall biosynthesis, offering a strategy for the genetic improvement of cotton fiber quality.

Dissecting the major genetic components underlying cotton lint development

Numerous genetic loci and several functionally characterized genes have been linked to determination of lint percentage (lint%), one of the most important cotton yield components, but we still know little about the major genetic components underlying lint%. Here, we first linked the genetic loci containing MYB25-like_At and HD1_At to the fiberless seed trait of 'SL1-7-1' and found that MYB25-like_At and HD1_At were very lowly expressed in 'SL1-7-1' ovules during fiber initiation. We then dissected the genetic components involved in determination of lint% using segregating populations derived from crosses of fuzzless mutants and intermediate segregants with different lint%, which not only confirmed the HD1_At locus but identified the HD1_Dt locus as being the major genetic components contributing to fiber initiation and lint%. The segregating populations also allowed us to evaluate the relative contributions of MYB25-like_At, MYB25-like_Dt, HD1_At, and HD1_Dt to lint%. Haplotype analysis of an Upland cotton (Gossypium hirsutum) population with 723 accessions (including 81 fuzzless seed accessions) showed that lint% of the accessions with the LP allele (higher lint%) at MYB25-like_At, MYB25-like_Dt, or HD1_At was significantly higher than that with the lp allele (lower lint%). The lint% of the Upland cotton accessions with 3 or 4 LP alleles at MYB25-like and HD1 was significantly higher than that with 2 LP alleles. The results prompted us to propose a strategy for breeding high-yielding cotton varieties, i.e. pyramiding the LP alleles of MYB25-like and HD1 with new lint% LP alleles without negative impact on seed size and fiber quality.

Single-cell RNA sequencing reveals a hierarchical transcriptional regulatory network of terpenoid biosynthesis in cotton secretory glandular cells

Plants can synthesize a wide range of terpenoids in response to various environmental cues. However, the specific regulatory mechanisms governing terpenoid biosynthesis at the cellular level remain largely elusive. In this study, we employed single-cell RNA sequencing to comprehensively characterize the transcriptome profile of cotton leaves and established a hierarchical transcriptional network regulating cell-specific terpenoid production. We observed substantial expression levels of genes associated with the biosynthesis of both volatile terpenes (such as β-caryophyllene and β-myrcene) and non-volatile gossypol-type terpenoids in secretory glandular cells. Moreover, two novel transcription factors, namely GoHSFA4a and GoNAC42, are identified to function downstream of the Gossypium PIGMENT GLAND FORMATION genes. Both transcription factors could directly regulate the expression of terpenoid biosynthetic genes in secretory glandular cells in response to developmental and environmental stimuli. For convenient retrieval of the single-cell RNA sequencing data generated in this study, we developed a user-friendly web server . Our findings not only offer valuable insights into the precise regulation of terpenoid biosynthesis genes in cotton leaves but also provide potential targets for cotton breeding endeavors.

Integrating transcriptome and physiological analyses to elucidate the molecular responses of sorghum to fluxofenim and metolachlor herbicide

The extensive use of herbicides has raised concerns about crop damage, necessitating the development of effective herbicide safeners. Fluxofenim has emerged as a promising herbicide safener; however, it's underlying mechanism remains unclear. Here, we screened two inbred lines 407B and HYZ to investigate the detoxication of fluxofenim in mitigating metolachlor damage in sorghum. Metolachlor inhibited seedling growth in both 407B and HYZ, while, fluxofenim could significantly restore the growth of 407B, but not effectively complement the growth of HYZ. Fluxofenim significantly increased the activities of glutathione-S-transferase (GST) to decrease metolachlor residue in 407B, but not in HYZ. This implys that fluxofenim may reduce metolachlor toxicity by regulating its metabolism. Furthermore, metolachlor suppressed AUX-related and JA-related genes expression, while up-regulated the expression of SA-related genes. Fluxofenim also restored the expression of AUX-related and JA-related genes inhibited by metolachlor and further increased expression of SA-related genes. Moreover, we noted a significant increase in the content of trans-zeatin O-glucoside (tZOG) and Gibberellin1 (GA1) after the fluxofenim treatment. In conclusion, fluxofenim may reduce the injury of herbicide by affecting herbicide metabolism and regulating hormone signaling pathway.

Regulatory controls of duplicated gene expression during fiber development in allotetraploid cotton

Polyploidy complicates transcriptional regulation and increases phenotypic diversity in organisms. The dynamics of genetic regulation of gene expression between coresident subgenomes in polyploids remains to be understood. Here we document the genetic regulation of fiber development in allotetraploid cotton Gossypium hirsutum by sequencing 376 genomes and 2,215 time-series transcriptomes. We characterize 1,258 genes comprising 36 genetic modules that control staged fiber development and uncover genetic components governing their partitioned expression relative to subgenomic duplicated genes (homoeologs). Only about 30% of fiber quality-related homoeologs show phenotypically favorable allele aggregation in cultivars, highlighting the potential for subgenome additivity in fiber improvement. We envision a genome-enabled breeding strategy, with particular attention to 48 favorable alleles related to fiber phenotypes that have been subjected to purifying selection during domestication. Our work delineates the dynamics of gene regulation during fiber development and highlights the potential of subgenomic coordination underpinning phenotypes in polyploid plants.

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.

The cotton miR171a-SCL6 module mediates plant resistance through regulating GhPR1 expression

Plants have developed intricate defense mechanisms in response to fluctuating environmental cues, including the use of microRNA (miRNA) as post-transcriptional regulators. However, the specific mechanisms through which miRNA contributes to disease resistance remain largely elusive. While the miR171-SCLs have been investigated in an eclectic array of plants, there has been a notable scarcity of research specifically focused on cotton (Gossypium hirsutum). In our previous miRNA-sequencing analysis, we found that ghr-miR171a displayed a differential response to infections by Verticillium dahliae. In this study, we further investigated the function of the miR171a-SCL6 module in cotton during V. dahliae infection. The ghr-miR171a was confirmed to direct the cleavage of GhSCL6 mRNA in the post-transcriptional process, as evidenced by 5' RLM-RACE, β-glucuronidase (GUS) histochemical staining and enzyme activity assay. Interestingly, we found that overexpressing ghr-miR171a reduced cotton plants' resistance to V. dahliae, while suppressing ghr-miR171a increased the plants' defense capacity. The GhSCL6 protein, when fused with green fluorescent protein (GFP), localizes in the cell nucleus, indicating its potential role in gene regulation. This was further corroborated by yeast two-hybrid assays, which verified GhSCL6's transcriptional activation ability. Through quantitative reverse transcriptase PCR (qRT-PCR), luciferase (LUC) fluorescence, and yeast one-hybrid assays, we found that GhSCL6 binds to the GT-box element of the GhPR1 promoter, activating its expression and thereby enhancing plant disease resistance. Taken together, our findings demonstrate that the cotton miR171a-SCL6 module regulates Verticillium wilt resistance in plants through the post-transcriptional process. This insight may offer new perspectives for disease resistance strategies in cotton.

The cotton MYB33 gene is a hub gene regulating the trade-off between plant growth and defense in Verticillium dahliae infection

INTRODUCTION

Sessile plants engage in trade-offs between growth and defense capacity in response to fluctuating environmental cues. MYB is an important transcription factor that plays many important roles in controlling plant growth and defense. However, the mechanism behind how it keeps a balance between these two physiological processes is still largely unknown.

OBJECTIVES

Our work focuses on the dissection of the molecular mechanism by which GhMYB33 regulates plant growth and defense.

METHODS

The CRISPR/Cas9 technique was used to generate mutants for deciphering GhMYB33 functions. Yeast two-hybrid, luciferase complementary imaging, and co-immunoprecipitation assays were used to prove that proteins interact with each other. We used the electrophoretic mobility shift assay, yeast one-hybrid, and luciferase activity assays to analyze GhMYB33 acting as a promoter. A β-glucuronidase fusion reporter and 5' RNA ligase mediated amplification of cDNA ends analysis showed that ghr-miR319c directedly cleaved the GhMYB33 mRNA.

RESULTS

Overexpressing miR319c-resistant GhMYB33 (rGhMYB33) promoted plant growth, accompanied by a significant decline in resistance against Verticillium dahliae. Conversely, its knockout mutant, ghmyb33, demonstrated growth restriction and concomitant augmentation of V. dahliae resistance. GhMYB33 was found to couple with the DELLA protein GhGAI1 and bind to the specific cis-elements of GhSPL9 and GhDFR1 promoters, thereby modulating internode elongation and plant resistance in V. dahliae infection. The ghr-miR319c was discovered to target and suppress GhMYB33 expression. The overexpression of ghr-miR319c led to enhanced plant resistance and a simultaneous reduction in plant height.

CONCLUSION

Our findings demonstrate that GhMYB33 encodes a hub protein and controls the expression of GhSPL9 and GhDFR1, implicating a pivotal role for the miR319c-MYB33 module to regulate the trade-offs between plant growth and defense.

Deciphering transcriptional mechanisms of maize internodal elongation by regulatory network analysis

The lengths of the basal internodes is an important factor for lodging resistance of maize (Zea mays). In this study, foliar application of coronatine (COR) to 10 cultivars at the V8 growth stage had different suppression effects on the length of the eighth internode, with three being categorized as strong-inhibition cultivars (SC), five as moderate (MC), and two as weak (WC). RNA-sequencing of the eighth internode of the cultivars revealed a total of 7895 internode elongation-regulating genes, including 777 transcription factors (TFs). Genes related to the hormones cytokinin, gibberellin, auxin, and ethylene in the SC group were significantly down-regulated compared to WC, and more cell-cycle regulatory factors and cell wall-related genes showed significant changes, which severely inhibited internode elongation. In addition, we used EMSAs to explore the direct regulatory relationship between two important TFs, ZmABI7 and ZmMYB117, which regulate the cell cycle and cell wall modification by directly binding to the promoters of their target genes ZmCYC1, ZmCYC3, ZmCYC7, and ZmCPP1. The transcriptome reported in this study will provide a useful resource for studying maize internode development, with potential use for targeted genetic control of internode length to improve the lodging resistance of maize.

Genome-wide analysis of the serine carboxypeptidase-like protein family reveals Ga09G1039 is involved in fiber elongation in cotton

The Gossypium is a model genus for understanding polyploidy and the evolutionary pattern of inheritance. This study aimed to investigate the characteristics of SCPLs in different cotton species and their role in fiber development. A total of 891 genes from one typical monocot and ten dicot species were naturally divided into three classes based on phylogenetic analysis. The SCPL gene family in cotton has undergone intense purifying selection with some functional variation. Segmental duplication and whole genome duplication were shown to be the two main reasons for the increase in the number of genes during cotton evolution. The identification of Gh_SCPL genes exhibiting differential expression in particular tissues or response to environmental stimuli provides a new measure for the in-depth characterization of selected genes of importance. Ga09G1039 was involved in the developmental process of fibers and ovules, and it is significantly different from proteins from other cotton species in terms of phylogenetic, gene structure, conserved protein motifs and tertiary structure. Overexpression of Ga09G1039 significantly increased the length of stem trichomes. Ga09G1039 may be a serine carboxypeptidase protein with hydrolase activity, according to functional region, prokaryotic expression, and western blotting analysis. The results provide a comprehensive overview of the genetic basis of SCPLs in Gossypium and further our knowledge in understanding the key aspects of SCPLs in cotton with their potential role in fiber development and stress resistance.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Overexpression of TCP9-like gene enhances salt tolerance in transgenic soybean

TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) transcription factors are a plant-specific family and play roles in plant growth, development, and responses to biotic and abiotic stresses. However, little is known about the functions of the TCP transcription factors in the soybean cultivars with tolerance to salt stress. In this study, TCP9-like, a TCP transcription factor, was identified in the soybean cultivars exposed to salt stress. The expression of TCP9-like gene in the roots of salt-tolerant soybean cultivars was higher than that in salt-sensitive cultivars treated with NaCl. The overexpression of TCP9-like enhanced the salt tolerance of the salt-sensitive soybean cultivar 'DN50'. In T2 generation, the plants with TCP9-like overexpression had significantly lower Na+ accumulation and higher K+ accumulation than the WT plants exposed to 200 or 250 mmol/L NaCl. The K+/Na+ ratio in the plants overexpressing TCP9-like was significantly higher than that in WT plants treated with 200 mmol/L NaCl. Meanwhile, the overexpression of TCP9-like up-regulated the expression levels of GmNHX1, GmNHX3, GmSOS1, GmSOS2-like, and GmHKT1, which were involved in the K+/Na+ homeostasis pathway. The findings indicated that TCP9-like mediated the regulation of both Na+ and K+ accumulation to improve the tolerance of soybean to salt stress.

Mapping intron retention events contributing to complex traits using splice quantitative trait locus

BACKGROUND

Alternative splicing (AS) of mRNA plays an important roles in transcriptome diversity, involving regulation of plant growth and stress response. Understanding the variation of AS events underlying GWAS loci in a crop population can provide insight into the molecular mechanisms of complex agronomic traits. To date, genome-wide association studies relating AS events to agronomic traits have rarely been conducted at the population level in crops.

RESULTS

Here, a pipeline was constructed to identify candidate AS events related to complex traits. Firstly, ovule transcriptome data were used to characterize intron retention (IR), the predominant type of AS in plants, on a genome-wide scale. This was done in a natural population consisting of 279 upland cotton lines. Secondly, splice quantitative trait locus (sQTL) analysis was carried out, which yielded a total of 2295 sQTLs involving 1607 genes. Of these, 14.25% (n = 427) were cis-regulatory loci. Integration with expression quantitative trait loci (eQTL) revealed that 53 (21.4%) cis-sGenes were regulated by both cis-sQTLs and cis-eQTLs. Finally, co-localization analysis integrated with GWAS loci in this population showed 32 cis-QTLs to be co-located with genetic regulatory loci related to fiber yield and quality traits, indicating that sQTLs are likely to participate in regulating cotton fiber yield and quality. An in-depth evaluation confirmed that differences in the IR rates of sQTL-regulated candidate genes such as GhLRRK1 and GhGC1 are associated with lint percentage (LP), which has potential in breeding applications.

CONCLUSION

This study provides a clue that AS of mRNA has an impact on crop yield, along with functional sQTLs are new genetic resources for cotton precision breeding.

A conserved brassinosteroid-mediated BES1-CERP-EXPA3 signaling cascade controls plant cell elongation

Continuous plant growth is achieved by cell division and cell elongation. Brassinosteroids control cell elongation and differentiation throughout plant life. However, signaling cascades underlying BR-mediated cell elongation are unknown. In this study, we introduce cotton fiber, one of the most representative single-celled tissues, to decipher cell-specific BR signaling. We find that gain of function of GhBES1, a key transcriptional activator in BR signaling, enhances fiber elongation. The chromatin immunoprecipitation sequencing analysis identifies a cell-elongation-related protein, GhCERP, whose transcription is directly activated by GhBES1. GhCERP, a downstream target of GhBES1, transmits the GhBES1-mediated BR signaling to its target gene, GhEXPA3-1. Ultimately, GhEXPA3-1 promotes fiber cell elongation. In addition, inter-species functional analysis of the BR-mediated BES1-CERP-EXPA3 signaling cascade also promotes Arabidopsis root and hypocotyl growth. We propose that the BES1-CERP-EXPA3 module may be a broad-spectrum pathway that is universally exploited by diverse plant species to regulate BR-promoted cell elongation.

Lint percentage and boll weight QTLs in three excellent upland cotton (Gossypium hirsutum): ZR014121, CCRI60, and EZ60

BACKGROUND

Upland cotton (Gossypium hirsutum L.) is the most economically important species in the cotton genus (Gossypium spp.). Enhancing the cotton yield is a major goal in cotton breeding programs. Lint percentage (LP) and boll weight (BW) are the two most important components of cotton lint yield. The identification of stable and effective quantitative trait loci (QTLs) will aid the molecular breeding of cotton cultivars with high yield.

RESULTS

Genotyping by target sequencing (GBTS) and genome-wide association study (GWAS) with 3VmrMLM were used to identify LP and BW related QTLs from two recombinant inbred line (RIL) populations derived from high lint yield and fiber quality lines (ZR014121, CCRI60 and EZ60). The average call rate of a single locus was 94.35%, and the average call rate of an individual was 92.10% in GBTS. A total of 100 QTLs were identified; 22 of them were overlapping with the reported QTLs, and 78 were novel QTLs. Of the 100 QTLs, 51 QTLs were for LP, and they explained 0.29-9.96% of the phenotypic variation; 49 QTLs were for BW, and they explained 0.41-6.31% of the phenotypic variation. One QTL (qBW-E-A10-1, qBW-C-A10-1) was identified in both populations. Six key QTLs were identified in multiple-environments; three were for LP, and three were for BW. A total of 108 candidate genes were identified in the regions of the six key QTLs. Several candidate genes were positively related to the developments of LP and BW, such as genes involved in gene transcription, protein synthesis, calcium signaling, carbon metabolism, and biosynthesis of secondary metabolites. Seven major candidate genes were predicted to form a co-expression network. Six significantly highly expressed candidate genes of the six QTLs after anthesis were the key genes regulating LP and BW and affecting cotton yield formation.

CONCLUSIONS

A total of 100 stable QTLs for LP and BW in upland cotton were identified in this study; these QTLs could be used in cotton molecular breeding programs. Putative candidate genes of the six key QTLs were identified; this result provided clues for future studies on the mechanisms of LP and BW developments.

Poplar coma morphogenesis and miRNA regulatory networks by combining ovary tissue sectioning and deep sequencing

Poplar coma, commonly referred to as "seed hairs", is a tuft of trichomes attached to the seed coat that helps seed dispersal. However, they can also trigger health impacts for humans, including sneezing, shortness of breath, and skin irritation. Despite efforts to study the regulatory mechanism of herbaceous trichome formation, poplar coma remains poorly understood. In this study, we showed that the epidermal cells of the funiculus and placenta are the origin of poplar coma based on observations of paraffin sections. Small RNA (sRNA) and degradome libraries were also constructed at three stages of poplar coma development, including initiation and elongation stages. Based on 7,904 miRNA-target pairs identified by small RNA and degradome sequencing, we constructed a miRNA-transcript factor and a stage-specific miRNA regulatory network. By combining paraffin section observation and deep sequencing, our research will provide greater insight into the molecular mechanisms of poplar coma development.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum

The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC₄F_3:5). A single segment substitution line (CSSL-106) screened from BC₆F₂ was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC₇F₂ individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

Cotton fiber as a model for understanding shifts in cell development under domestication

Cotton fiber provides the predominant plant textile in the world, and it is also a model for plant cell wall biosynthesis. The development of the single-celled cotton fiber takes place across several overlapping but discrete stages, including fiber initiation, elongation, the transition from elongation to secondary cell wall formation, cell wall thickening, and maturation and cell death. During each stage, the developing fiber undergoes a complex restructuring of genome-wide gene expression change and physiological/biosynthetic processes, which ultimately generate a strikingly elongated and nearly pure cellulose product that forms the basis of the global cotton industry. Here, we provide an overview of this developmental process focusing both on its temporal as well as evolutionary dimensions. We suggest potential avenues for further improvement of cotton as a crop plant.

The RNAome landscape of tomato during arbuscular mycorrhizal symbiosis reveals an evolving RNA layer symbiotic regulatory network

Arbuscular mycorrhizal symbiosis (AMS) is an ancient plant-fungus relationship that is widely distributed in terrestrial plants. The formation of symbiotic structures and bidirectional nutrient exchange requires the regulation of numerous genes. However, the landscape of RNAome during plant AMS involving different types of regulatory RNA is poorly understood. In this study, a combinatorial strategy utilizing multiple sequencing approaches was used to decipher the landscape of RNAome in tomato, an emerging AMS model. The annotation of the tomato genome was improved by a multiple-platform sequencing strategy. A total of 3,174 protein-coding genes were upregulated during AMS, 42% of which were alternatively spliced. Comparative-transcriptome analysis revealed that genes from 24 orthogroups were consistently induced by AMS in eight phylogenetically distant angiosperms. Seven additional orthogroups were specifically induced by AMS in all surveyed dicot AMS host plants. However, these orthogroups were absent or not induced in monocots and/or non-AMS hosts, suggesting a continuously evolving AMS-responsive network in addition to a conserved core regulatory module. Additionally, we detected 587 lncRNAs, ten miRNAs, and 146 circRNAs that responded to AMS, which were incorporated to establish a tomato AMS-responsive, competing RNA-responsive endogenous RNA (ceRNA) network. Finally, a tomato symbiotic transcriptome database (TSTD, https://efg.nju.edu.cn/TSTD) was constructed to serve as a resource for deep deciphering of the AMS regulatory network. These results help elucidate the reconfiguration of the tomato RNAome during AMS and suggest a sophisticated and evolving RNA layer responsive network during AMS processes.

Identification and function of miRNA-mRNA interaction pairs during lateral root development of hemi-parasitic Santalum album L. seedlings

Sandalwood (Santalum album L.) is a hemi-parasitic tree species famous for its santalol and santalene, which are extracted from its heartwood and roots. The ability to understand root functionality within its branched root system would benefit the regulation of sandalwood growth and enhance the commercial value of sandalwood. Phenotypic and anatomical evidence in this study indicated that seed germination stage 4 (SG4) seemed pivotal for lateral root (LR) morphogenesis. Small RNA (sRNA) high-throughput sequencing of root tissues at three sub-stages of SG4 (lateral root primordia initiation (LRPI), lateral root primordia development (LRPD), and lateral root primordia emergence (LRPE)) was performed to identify microRNAs (miRNAs) associated with LR development. A total of 135 miRNAs, including 70 differentially expressed miRNAs (DEMs), were screened. Ten DEMs were selected to investigate transcript abundance in different organs or developmental stages. Among 100 negative DEM-mRNA interaction pairs, four targets (Sa-miR166m_2, 408d, 858a, and novel_Sa-miR8) were selected for studying cleavage sites by 5' RLM-RACE validation. The expression mode of the four miRNA-mRNA pairs was investigated after indole-3-acetic acid (IAA) treatment. IAA enhanced the abundance of homeobox-leucine-zipper protein 32 (HOX32), laccase 12 (LAC12), myeloblastosis86 (MYB86), and pectin methylesterase inhibitor6 (PMEI6) target transcripts by reducing the expression of Sa-miR166m_2, 408d, 858a, and novel_Sa-miR8 in the first 10 min. A schematic model of miRNA-regulated LR development is proposed for this hemi-parasitic species. This novel genetic information for improving sandalwood root growth and development may allow for the cultivation of fast-growing and high-yielding plantations.

Fine-mapping and candidate gene analysis of qFL-c10-1 controlling fiber length in upland cotton (Gossypium hirsutum L.)

A fiber length QTL, qFL-c10-1, was fine-mapped to a 96.5-kb region containing one gene that has not been characterized in plants. Fiber length is an important component of cotton fiber quality, which is associated with other quality properties such as fiber strength, fiber maturity, and fineness. In our previous studies, a stable QTL qFL-c10-1 controlling fiber length had been identified on chromosome A10 in an upland cotton recombinant inbred line (RIL) population from a cross between Jimian5 and DH962. To fine-map qFL-c10-1, an F₂ population with 1081 individual plants from a cross between a recombinant line DJ61 and Jimian5 was established. Using linkage analysis and progeny recombination experiment, qFL-c10-1 was mapped into a 96.5-kb genomic region that just contained one proper transcript Ghir_A10G022020 (described as GhFL10), an undescribed gene in plants. One 214-bp deletion was identified in the promoter region of DJ61 compared with Jimian5. Quantitative real-time PCR (qRT-PCR) and comparative analysis of parental sequences suggested that GhFL10 was the most promising candidate gene for qFL-c10-1. According to RNA-seq, yeast two-hybrid assay and bimolecular fluorescence complementation (BiFC), we speculate that GhFL10 interacts with NF-YA transcription factors to negatively regulate fiber elongation.

Single-cell RNA-seq reveals fate determination control of an individual fibre cell initiation in cotton (Gossypium hirsutum)

Cotton fibre is a unicellular seed trichome, and lint fibre initials per seed as a factor determines fibre yield. However, the mechanisms controlling fibre initiation from ovule epidermis are not understood well enough. Here, with single-cell RNA sequencing (scRNA-seq), a total of 14 535 cells were identified from cotton ovule outer integument of Xu142_LF line at four developmental stages (1.5, 1, 0.5 days before anthesis and the day of anthesis). Three major cell types, fibre, non-fibre epidermis and outer pigment layer were identified and then verified by RNA in situ hybridization. A comparative analysis on scRNA-seq data between Xu142 and its fibreless mutant Xu142 fl further confirmed fibre cluster definition. The developmental trajectory of fibre cell was reconstructed, and fibre cell was identified differentiated at 1 day before anthesis. Gene regulatory networks at four stages revealed the spatiotemporal pattern of core transcription factors, and MYB25-like and HOX3 were demonstrated played key roles as commanders in fibre differentiation and tip-biased diffuse growth respectively. A model for early development of a single fibre cell was proposed here, which sheds light on further deciphering mechanism of plant trichome and the improvement of cotton fibre yield.

Strigolactones act downstream of gibberellins to regulate fiber cell elongation and cell wall thickness in cotton (Gossypium hirsutum)

Strigolactones (SLs) are a class of phytohormones that regulate plant shoot branching and adventitious root development. However, little is known regarding the role of SLs in controlling the behavior of the smallest unit of the organism, the single cell. Here, taking advantage of a classic single-cell model offered by the cotton (Gossypium hirsutum) fiber cell, we show that SLs, whose biosynthesis is fine-tuned by gibberellins (GAs), positively regulate cell elongation and cell wall thickness by promoting the biosynthesis of very long-chain fatty acids (VLCFAs) and cellulose, respectively. Furthermore, we identified two layers of transcription factors (TFs) involved in the hierarchical regulation of this GA-SL crosstalk. The top-layer TF GROWTH-REGULATING FACTOR 4 (GhGRF4) directly activates expression of the SL biosynthetic gene DWARF27 (D27) to increase SL accumulation in fiber cells and GAs induce GhGRF4 expression. SLs induce the expression of four second-layer TF genes (GhNAC100-2, GhBLH51, GhGT2, and GhB9SHZ1), which transmit SL signals downstream to two ketoacyl-CoA synthase genes (KCS) and three cellulose synthase (CesA) genes by directly activating their transcription. Finally, the KCS and CesA enzymes catalyze the biosynthesis of VLCFAs and cellulose, respectively, to regulate development of high-grade cotton fibers. In addition to providing a theoretical basis for cotton fiber improvement, our results shed light on SL signaling in plant development at the single-cell level.

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

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

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

The sea-island cotton GbTCP4 transcription factor positively regulates drought and salt stress responses

TCP transcription factors play important regulatory roles in plant growth and development; however, their function in response to salt and drought stress in sea-island cotton (Gossypium barbadense) is unknown. Here, GbTCP4 expression was induced by abscisic acid (ABA), drought, and NaCl treatments. Under drought stress, compared to wild-type (WT) Arabidopsis, transgenic GbTCP4-overexpressing Arabidopsis showed increased seed germination rate, root length and survival rate; additionally, it was ABA-insensitive at the germination stage but ABA-sensitive at the seedling stage, showing reduced stomatal opening and ABA enrichment. Under salt stress, compared to WT Arabidopsis, transgenic GbTCP4-overexpressing Arabidopsis showed greater root length, survival rate, and SPAD value and lower malondialdehyde (MDA) content. Conversely, under drought or salt stress, virus-induced gene-silenced GbTCP4 cotton showed decreased root length, area and volume and increased MDA content and sensitivity to drought and salt stress compared with control cotton. RNA-seq and quantitative real-time PCR analyses showed that GbTCP4 affected the transcription levels of genes across multiple abiotic stress-related metabolic pathways. Furthermore, GbTCP4 activated the transcription of GbUVR8 and GbbHLH130 by binding to their promoters. These results suggest that GbTCP4 positively regulates drought and salt stress responses and is a suitable candidate gene for improving plant drought and salt tolerance.

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

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

The miR319/TaGAMYB3 module regulates plant architecture and improves grain yield in common wheat (Triticum aestivum)

Plant architecture is a key determinant of crop productivity and adaptation. The highly conserved microRNA319 (miR319) family functions in various biological processes, but little is known about how miR319 regulates plant architecture in wheat (Triticum aestivum). Here, we determined that the miR319/TaGAMYB3 module controls plant architecture and grain yield in common wheat. Repressing tae-miR319 using short tandem target mimics resulted in favorable plant architecture traits, including increased plant height, reduced tiller number, enlarged spikes and flag leaves, and thicker culms, as well as enhanced grain yield in field plot tests. Overexpressing tae-miR319 had the opposite effects on plant architecture and grain yield. Although both TaPCF8 and TaGAMYB3 were identified as miR319 target genes, genetic complementation assays demonstrated that only miR319-resistant TaGAMYB3 (rTaGAMYB3) abolished tae-miR319-mediated growth inhibition of flag leaves and spikes. TaGAMYB3 functions as a transcriptional activator of downstream genes, including TaPSKR1, TaXTH23, TaMADS5 and TaMADS51, by binding to their promoters. Furthermore, TaGAMYB3 physically interacts with TaBA1, an important regulator of spike development, to additively activate the transcription of downstream genes such as TaMADS5. Our findings provide insight into how the miR319/TaGAMYB3 module regulates plant architecture and improves grain yield in common wheat.

Systematical characterization of GRF gene family in sorghum, and their potential functions in aphid resistance

Sorghum (Sorghum bicolor) is the fifth important cereal and an industrial energy crop in the world. Growth Regulation Factors (GRFs) play an important role in response to environmental stress, however, the knowledge of GRFs relating to the pest resistance is lacking. Here, we identified 8 GRF genes harboring the typical QLQ (glutamine, leucine, glutamine) and WRC (tryptophan, arginine, cysteine) domains in Sorghum, which could be classified into 4 clades through phylogenetic analysis. The SbGRF genes express in most tissues, while more than half of them express at the highest level in inflorescence. To further investigate their possible role in stress response, we analyzed the transcriptomics data. The results showed that SbGRFs could respond to the abiotic stresses including heat, salt and drought stress. Furthermore, combined the data with qRT-PCR, SbGRF1, 2, 4 and 7 were identified as dominant genes response to the aphid-induced stress. SSR markers close to these genes were also searched. Above all, we summarized the SbGRFs and provided their potential roles in aphid response.

CRISPR/Cas9-mediated P-CR domain-specific engineering of CESA4 heterodimerization capacity alters cell wall architecture and improves saccharification efficiency in poplar

Cellulose is the most abundant unique biopolymer in nature with widespread applications in bioenergy and high-value bioproducts. The large transmembrane-localized cellulose synthase (CESA) complexes (CSCs) play a pivotal role in the biosynthesis and orientation of the para-crystalline cellulose microfibrils during secondary cell wall (SCW) deposition. However, the hub CESA subunit with high potential homo/heterodimerization capacity and its functional effects on cell wall architecture, cellulose crystallinity, and saccharification efficiency remains unclear. Here, we reported the highly potent binding site containing four residues of Pro435, Trp436, Pro437, and Gly438 in the plant-conserved region (P-CR) of PalCESA4 subunit, which are involved in the CESA4-CESA8 heterodimerization. The CRISPR/Cas9-knockout mutagenesis in the predicted binding site results in physiological abnormalities, stunt growth, and deficient roots. The homozygous double substitution of W436Q and P437S and heterozygous double deletions of W436 and P437 residues potentially reduced CESA4-binding affinity resulting in normal roots, 1.5-2-fold higher plant growth and cell wall regeneration rates, 1.7-fold thinner cell wall, high hemicellulose content, 37%-67% decrease in cellulose content, high cellulose DP, 25%-37% decrease in cellulose crystallinity, and 50% increase in saccharification efficiency. The heterozygous deletion of W436 increases about 2-fold CESA4 homo/heterodimerization capacity led to the 50% decrease in plant growth and increase in cell walls thickness, cellulose content (33%), cellulose DP (20%), and CrI (8%). Our findings provide a strategy for introducing commercial CRISPR/Cas9-mediated bioengineered poplars with promising cellulose applications. We anticipate our results could create an engineering revolution in bioenergy and cellulose-based nanomaterial technologies.

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.

The Characters of Non-Coding RNAs and Their Biological Roles in Plant Development and Abiotic Stress Response

Plant growth and development are greatly affected by the environment. Many genes have been identified to be involved in regulating plant development and adaption of abiotic stress. Apart from protein-coding genes, more and more evidence indicates that non-coding RNAs (ncRNAs), including small RNAs and long ncRNAs (lncRNAs), can target plant developmental and stress-responsive mRNAs, regulatory genes, DNA regulatory regions, and proteins to regulate the transcription of various genes at the transcriptional, posttranscriptional, and epigenetic level. Currently, the molecular regulatory mechanisms of sRNAs and lncRNAs controlling plant development and abiotic response are being deeply explored. In this review, we summarize the recent research progress of small RNAs and lncRNAs in plants, focusing on the signal factors, expression characters, targets functions, and interplay network of ncRNAs and their targets in plant development and abiotic stress responses. The complex molecular regulatory pathways among small RNAs, lncRNAs, and targets in plants are also discussed. Understanding molecular mechanisms and functional implications of ncRNAs in various abiotic stress responses and development will benefit us in regard to the use of ncRNAs as potential character-determining factors in molecular plant breeding.

Multi-Dimensional Molecular Regulation of Trichome Development in Arabidopsis and Cotton

Plant trichomes are specialized epidermal cells that are widely distributed on plant aerial tissues. The initiation and progression of trichomes are controlled in a coordinated sequence of multiple molecular events. During the past decade, major breakthroughs in the molecular understanding of trichome development were achieved through the characterization of various trichomes defective mutants and trichome-associated genes, which revealed a highly complex molecular regulatory network underlying plant trichome development. This review focuses on the recent millstone in plant trichomes research obtained using genetic and molecular studies, as well as 'omics' analyses in model plant Arabidopsis and fiber crop cotton. In particular, we discuss the latest understanding and insights into the underlying molecular mechanisms of trichomes formation at multiple dimensions, including at the chromatin, transcriptional, post-transcriptional, and post-translational levels. We summarize that the integration of multi-dimensional trichome-associated genes will enable us to systematically understand the molecular regulation network that landscapes the development of the plant trichomes. These advances will enable us to address the unresolved questions regarding the molecular crosstalk that coordinate concurrent and ordered the changes in cotton fiber initiation and progression, together with their possible implications for genetic improvement of cotton fiber.

Identification of key proteins related to high-quality fiber in Upland cotton via proteomics analysis

The dynamics of cotton fiber elongation and microfibirl deposition orientation were delineated; advanced ethylene synthesis and redox reaction homeostasis may be crucial for high-quality fiber formation. Fiber length, strength, and fineness determine the use and commercial value of cotton fiber, but their underlying molecular mechanisms remain obscure. We compared the dynamic change trajectories of length, diameter and microfibril orientation angle of the fibers produced by an introgression line SY6167 which generates high-quality fibers even better than Sea island cotton with those of the common-quality fibers from TM-1 across 5 to 30 days post anthesis (DPA). The proteomes were profiled and compared at six representative time points using 2-DE and MS/MS. 14 proteins differentially expressed inside each of cotton line temporally and significantly different tween the two lines were identified. The dynamic change trajectories of fiber length and microfibril angle are close to "s" and reverse "s" growth curves, respectively. SY6167 and TM-1 fibers entered the logarithmic elongation phase simultaneously at 10 DPA, and SY6167 kept elongating logarithmically for 2 more days than TM-1. In parallel to logarithmic elongation, microfibril orientation angles dived sharply, and SY6167 declined faster for a shorter duration than TM-1. 53% of the identified proteins are related to redox homeostasis, and most of them are expressed at higher levels in SY6167 during logarithmic elongation. 1-Aminocyclopropane-1-Carboxylic Acid Oxidase (ACO) started to accumulate at 16 DPA in SY6167, and its encoding genes were highly expressed at this stage, with a much higher level than TM-1. These findings suggest high-quality fibers are associated with high expression of the proteins related to stress and redox homeostasis, the continuously elevated expression of ethylene synthesis ACO gene may play an essential role.

Inheritance, QTLs, and Candidate Genes of Lint Percentage in Upland Cotton

Cotton (Gossypium spp.) is an important natural fiber plant. Lint percentage (LP) is one of the most important determinants of cotton yield and is a typical quantitative trait with high variation and heritability. Many cotton LP genetic linkages and association maps have been reported. This work summarizes the inheritance, quantitative trait loci (QTLs), and candidate genes of LP to facilitate LP genetic study and molecular breeding. More than 1439 QTLs controlling LP have been reported. Excluding replicate QTLs, 417 unique QTLs have been identified on 26 chromosomes, including 243 QTLs identified at LOD >3. More than 60 are stable, major effective QTLs that can be used in marker-assisted selection (MAS). More than 90 candidate genes for LP have been reported. These genes encode MYB, HOX, NET, and other proteins, and most are preferentially expressed during fiber initiation and elongation. A putative molecular regulatory model of LP was constructed and provides the foundation for the genetic study and molecular breeding of LP.

Comprehensive identification and expression analysis of CRY gene family in Gossypium

BACKGROUND

The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance.

RESULTS

In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments.

CONCLUSION

We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

Blocking IbmiR319a Impacts Plant Architecture and Reduces Drought Tolerance in Sweet Potato

MicroRNA319 (miR319) plays a key role in plant growth, development, and multiple resistance by repressing the expression of targeted TEOSINTE BRANCHED/CYCLOIDEA/PCF (TCP) genes. Two members, IbmiR319a and IbmiR319c, were discovered in the miR319 gene family in sweet potato (Ipomoea batatas [L.] Lam). Here, we focused on the biological function and potential molecular mechanism of the response of IbmiR319a to drought stress in sweet potato. Blocking IbmiR319a in transgenic sweet potato (MIM319) resulted in a slim and tender phenotype and greater sensitivity to drought stress. Microscopic observations revealed that blocking IbmiR319a decreased the cell width and increased the stomatal distribution in the adaxial leaf epidermis, and also increased the intercellular space in the leaf and petiole. We also found that the lignin content was reduced, which led to increased brittleness in MIM319. Quantitative real-time PCR showed that the expression levels of key genes in the lignin biosynthesis pathway were much lower in the MIM319 lines than in the wild type. Ectopic expression of IbmiR319a-targeted genes IbTCP11 and IbTCP17 in Arabidopsis resulted in similar phenotypes to MIM319. We also showed that the expression of IbTCP11 and IbTCP17 was largely induced by drought stress. Transcriptome analysis indicated that cell growth-related pathways, such as plant hormonal signaling, were significantly downregulated with the blocking of IbmiR319a. Taken together, our findings suggest that IbmiR319a affects plant architecture by targeting IbTCP11/17 to control the response to drought stress in sweet potato.

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.