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

GhLPF1 Associated Network Is Involved with Cotton Lint Percentage Regulation Revealed by the Integrative Analysis of Spatial Transcriptome

Cotton fibers, derived from the epidermis of the ovule, provide a sustainable natural fiber source for the textile industry. Traits related to fiber yield are predominantly determined by molecular regulations in the epidermis of the outer integument (OI) region of the cotton ovule. Here, we identify an R2R3 MYB transcription factor coding gene GhLPF1 within the QTL-LP-ChrA06 locus for lint percentage (LP, percentage of lint to seed cotton) through constructing the 1-Day Post Anthesis Cotton Ovule Spatial Transcriptome Atlas. GhLPF1 is subjected as a downstream target of miR828 during fiber development. The direct downstream genes (DDGs) of GhLPF1 are biased to increased expression in GhLPF1-CR, and are preferentially expressed in OI, so that GhLPF1 is primarily a transcriptional repressor to its DDGs. Population-wide transcriptome analysis confirms that expression variation of GhLPF1-DDGs is significantly biased to negative correlation with LP, among which a type I homeobox protein-coding gene GhHB6 is further validated to be the directly downstream gene of GhLPF1. Given these data, it is demonstrated that GhLPF1 mediates a regulation network in LP as a transcriptional repressor, which makes it a valuable functional marker for fiber-trait improvement application from QTL-LP-ChrA06.

The speed breeding technology of five generations per year in cotton

KEY MESSAGE

Developed a speed breeding technique for cotton that enables up to five generations per year using optimized spectral conditions and immature embryo culture, and created new materials with iaaM gene. Shortening the breeding cycle is an effective way to accelerate crop genetic improvement. Previously we developed an integrated breeding technology for cotton that enabled three to four breeding cycles per year. Here, to further shorten the breeding time, we optimized the light spectrum conditions for cotton development and culture conditions for immature embryo developing into seedling. Under optimized spectrum conditions, JSh929 and ND601 plants exhibited the visible flower buds at 19 and 21 days after emergence (DAE), and the first flower bloomed at around 45 and 46 DAE. Using the optimized immature embryo culture technique, immature embryos of 25-30 days after pollination could develop into fertile plants with cotyledon unfolding at 6 days after culture in vitro. The improved speed breeding technique shortened cotton breeding cycle from about 130 days to a range from 71 to 85 days, an average of 79.5 days, achieving up to around five generations per year. Using this optimized system, we transferred iaaM gene into the high-yield and disease-resistant cultivar JND24, and BC4F3 progenies were obtained within 1.5 years. In addition, the JND24-i3 line was selected with increased lint percentage and improved Micronaire value. These results demonstrate that the optimized speed breeding system offers a rapid and effective way to improve traits of cotton.

NAC transcription factor GbNTL9 modifies the accumulation and organization of cellulose microfibrils to enhance cotton fiber strength

INTRODUCTION

Fiber strength is a critical determinant of fiber quality, with stronger fibers being highly preferred in the cotton textile industry. However, the genetic basis and the specific regulatory mechanism underlying the formation of cotton fiber strength remain largely unknown.

OBJECTIVES

To explore fiber strength-related genes, QTL mapping, map-based cloning, and gene function verification were conducted in a backcross inbred line BS41 derived from interspecific hybridization between upland cotton and sea-island cotton.

METHODS

Upland cotton Emian22 (E22) and an interspecific backcross inbred line (BIL) BS41 were used as parents to construct secondary segregation populations for BSA and QTL mapping of fiber strength. The candidate gene GbNTL9 was identified through map-based cloning and expression analysis. The function of NTL9 was determined through transgenic experiments and cytological observations. The regulatory mechanisms of NTL9 were explored using RNA-seq, RT-qPCR, yeast two-hybrid, bimolecular fluorescence complementation, and yeast one-hybrid.

RESULTS

A major QTL for fiber strength, qFS-A11-1, was mapped to a 14.6-kb genomic region using segregating populations from E22 × BS41. GbNTL9, which encodes a NAC transcription factor, was identified as the candidate gene. Overexpression of both upland cotton genotype NTL9E22 and sea-island genotype NTL9BS41 in upland cotton enhanced fiber strength by facilitating the dense accumulation and orderly organization of cellulose microfibrils within the cell wall. Transcriptomic analysis revealed that NTL9 inhibited the expression of genes involved in secondary wall synthesis, such as CESA4, CESA7, and CESA8, thereby delaying cell wall cellulose deposition and altering the microfibril deposition pattern. NTL9 interacted with MYB6 and functioned as a downstream gene in the ethylene signaling pathway. Additionally, an effective gene marker NTL9-24 was developed to distinguish haplotypes from G. barbadense and G. hirsutum for fiber quality breeding program.

CONCLUSION

Our findings demonstrate that GbNTL9 positively regulates fiber strength through altering the microfibril deposition pattern, and provide a new insight into the molecular mechanism underlying fiber strength.

The Disruptions of Sphingolipid and Sterol Metabolism in the Short Fiber of Ligon-Lintless-1 Mutant Revealed Obesity Impeded Cotton Fiber Elongation and Secondary Cell Wall Deposition

Boosting evidence indicated lipids play important roles in plants. To explore lipid function in cotton fiber development, the lipid composition and content were detected by untargeted and targeted lipidomics. Compared with rapid elongation fibers, the lipid intensity of 16 sub-classes and 56 molecular species decreased, while only 7 sub-classes and 26 molecular species increased in the fibers at the stage of secondary cell wall deposition. Unexpectedly, at the rapid elongation stage, 20 sub-classes and 60 molecular species increased significantly, while only 5 sub-classes and 8 molecular species decreased in the ligon lintless-1 (li-1) mutant compared with its wild-type Texas Maker-1 (TM-1). Particularly, campesteryl, sitosteryl, and total steryl ester increased by 21.8-, 48.7-, and 45.5-fold in the li-1 fibers, respectively. All the molecular species of sphingosine-1-P, phytoceramide-OHFA, and glucosylceramide increased while all sphingosine, phytosphingosine, and glycosyl inositol phospho ceramides decreased in the li-1 fibers. Similarly, the different expression genes between the mutant and wild type were enriched in many pathways involved in the lipid metabolism. Furthermore, the number of lipid droplets also increased in the li-1 leaf and fiber cells when compared with the wild type. These results illuminated that fiber cell elongation being blocked in the li-1 mutant was not due to a lack of lipids, but rather lipid over-accumulation (obesity), which may result from the disruption of sphingolipid and sterol metabolism. This study provides a new perspective for further studying the regulatory mechanisms of fiber development.

Regulatory networks of coresident subgenomes during rapid fiber cell elongation in upland cotton

Cotton, an intriguing plant species shaped by polyploidization, evolution, and domestication, holds particular interest due to the complex mechanisms governing fiber traits across its two subgenomes. However, the regulatory elements or transcriptional networks between subgenomes during fiber elongation remain to be fully clarified. Here, we analyzed 1462 cotton fiber samples to reconstruct the gene-expression regulatory networks that influence fiber cell elongation. Inter-subgenome expression quantitative trait loci (eQTLs) largely dictate gene transcription, with a notable tendency for the D subgenome to regulate A-subgenome eGenes. This regulation reveals synchronized homoeologous gene expression driven by co-localized eQTLs and divergent patterns that diminish genetic correlations, thus leading to preferential expression in the A and D subgenomes. Hotspot456 emerged as a key regulator of fiber initiation and elongation, and artificial selection of trans-eQTLs in hotspot456 that positively regulate KCS1 has facilitated cell elongation. Experiments designed to clarify the roles of trans-eQTLs in improved fiber breeding confirmed the inhibition of GhTOL9 by a specific trans-eQTL via GhWRKY28, which negatively affects fiber elongation. We propose a model in which the GhWRKY28-GhTOL9 module regulates this process through the ESCRT (endosomal sorting complex required for transport) pathway. This research significantly advances our understanding of cotton's evolutionary and domestication processes and the intricate regulatory mechanisms that underlie significant plant traits.

Brassinosteroids in cotton: orchestrating fiber development

Cotton cultivation spans over 30 million hectares across 85 countries and regions, with more than half participating in the global cotton textile trade. The elongated cotton fiber cell is an ideal model for studying cell elongation and understanding plant growth and development. Brassinosteroids (BRs), recognized for their role in cell elongation, offer the potential for improving cotton fiber quality and yield. Despite extensive research highlighting BR's positive impact on fiber development, a comprehensive review on this topic has been lacking. This review addresses this gap, providing a detailed analysis of the latest advancements in BR signaling and its effects on cotton fiber development. We explore the complex network of BR biosynthesis components, signaling molecules, and regulators, including crosstalk with other pathways and transcriptional control mechanisms. Additionally, we propose molecular strategies and highlight key genetic elements for optimizing BR-related genes to enhance fiber quality and yield. The review emphasizes the importance of BR homeostasis and the hormonal landscape during cotton fiber development, offering insights into targeted manipulation opportunities and challenges. This consolidation offers a comprehensive understanding of BR's multifaceted roles in fiber development, outlining a strategic approach for BR optimization in cotton fiber quality and yield.

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development

INTRODUCTION

Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations.

OBJECTIVES

The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield.

METHODS

A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies.

RESULTS

An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage.

CONCLUSION

Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

Comparative Lipidomics Analysis Provides New Insights into the Metabolic Basis of Color Formation in Green Cotton Fiber

Green fiber (GF) is a naturally colored fiber. A limited understanding of its color formation mechanism restricts the improvement of colored cotton quality. This experiment used upland cotton green fiber germplasm 1-4560 and genetic inbred line TM-1; the lipid profiles of green fibers at 30 (white stage) and 35 days post-anthesis (DPA) (early greening stage), as well as those of TM-1 at the same stages, were revealed. Among the 109 differential types of lipids (DTLs) unique to GF, the content of phosphatidylserine PS (16:0_18:3) was significantly different at 30 and 35 DPA. It is speculated that this lipid is crucial for the pigment accumulation and color formation process of green fibers. The 197 DTLs unique to TM-1 may be involved in white fiber (WF) development. Among the shared DTLs in GF35 vs. GF30 and WF35 vs. WF30, sulfoquinovosyldiacyl-glycerol SQDG (18:1_18:1) displays a significant difference in the content change between green fibers and white fibers, potentially affecting color formation through changes in content. The enriched metabolic pathways in both comparison groups are relatively conserved. In the most significantly enriched glycerophospholipid metabolic pathway, 1-acyl-sn-glycero-3-phosphocholine (C04230) only appears in white cotton. This indicates differences in the metabolic pathways between white and green fibers, potentially related to different mechanisms of color formation and fiber development. These findings provide a new theoretical basis for studying cotton fiber development and offer important insights into the specific mechanism of green fiber color formation.

Downregulation of the GhROD1 Gene Improves Cotton Fiber Fineness by Decreasing Acyl Pool Saturation, Stimulating Small Heat Shock Proteins (sHSPs), and Reducing H2O2 Production

Cotton fiber is one of the most important natural fiber sources in the world, and lipid metabolism plays a critical role in its development. However, the specific role of lipid molecules in fiber development and the impact of fatty acid alterations on fiber quality remain largely unknown. In this study, we demonstrate that the downregulation of GhROD1, a gene encoding phosphatidylcholine diacylglycerol cholinephosphotransferase (PDCT), results in an improvement of fiber fineness. We found that GhROD1 downregulation significantly increases the proportion of linoleic acid (18:2) in cotton fibers, which subsequently upregulates genes encoding small heat shock proteins (sHSPs). This, in turn, reduces H2O2 production, thus delaying secondary wall deposition and leading to finer fibers. Our findings reveal how alterations in linoleic acid influence cellulose synthesis and suggest a potential strategy to improve cotton fiber quality by regulating lipid metabolism pathways.

GhATL68b regulates cotton fiber cell development by ubiquitinating the enzyme required for β-oxidation of polyunsaturated fatty acids

E3 ligases are key enzymes required for protein degradation. Here, we identified a C3H2C3 RING domain-containing E3 ubiquitin ligase gene named GhATL68b. It is preferentially and highly expressed in developing cotton fiber cells and shows greater conservation in plants than in animals or archaea. The four orthologous copies of this gene in various diploid cottons and eight in the allotetraploid G. hirsutum were found to have originated from a single common ancestor that can be traced back to Chlamydomonas reinhardtii at about 992 million years ago. Structural variations in the GhATL68b promoter regions of G. hirsutum, G. herbaceum, G. arboreum, and G. raimondii are correlated with significantly different methylation patterns. Homozygous CRISPR-Cas9 knockout cotton lines exhibit significant reductions in fiber quality traits, including upper-half mean length, elongation at break, uniformity, and mature fiber weight. In vitro ubiquitination and cell-free protein degradation assays revealed that GhATL68b modulates the homeostasis of 2,4-dienoyl-CoA reductase, a rate-limiting enzyme for the β-oxidation of polyunsaturated fatty acids (PUFAs), via the ubiquitin proteasome pathway. Fiber cells harvested from these knockout mutants contain significantly lower levels of PUFAs important for production of glycerophospholipids and regulation of plasma membrane fluidity. The fiber growth defects of the mutant can be fully rescued by the addition of linolenic acid (C18:3), the most abundant type of PUFA, to the ovule culture medium. This experimentally characterized C3H2C3 type E3 ubiquitin ligase involved in regulating fiber cell elongation may provide us with a new genetic target for improved cotton lint production.

The underlying molecular mechanisms of hormonal regulation of fruit color in fruit-bearing plants

Fruit color is a key feature of fruit quality, primarily influenced by anthocyanin or carotenoid accumulation or chlorophyll degradation. Adapting the pigment content is crucial to improve the fruit's nutritional and commercial value. Genetic factors along with other environmental components (i.e., light, temperature, nutrition, etc.) regulate fruit coloration. The fruit coloration process is influenced by plant hormones, which also play a vital role in various physiological and biochemical metabolic processes. Additionally, phytohormones play a role in the regulation of a highly conserved transcription factor complex, called MBW (MYB-bHLH-WD40). The MBW complex, which consists of myeloblastosis (MYB), basic helix-loop-helix (bHLH), and WD40 repeat (WDR) proteins, coordinates the expression of downstream structural genes associated with anthocyanin formation. In fruit production, the application of plant hormones may be important for promoting coloration. However, concerns such as improper concentration or application time must be addressed. This article explores the molecular processes underlying pigment formation and how they are influenced by various plant hormones. The ABA, jasmonate, and brassinosteroid increase anthocyanin and carotenoid formation, but ethylene, auxin, cytokinin, and gibberellin have positive as well as negative effects on anthocyanin formation. This article establishes the necessary groundwork for future studies into the molecular mechanisms of plant hormones regulating fruit color, ultimately aiding in their effective and scientific application towards fruit coloration.

A telomere-to-telomere cotton genome assembly reveals centromere evolution and a Mutator transposon-linked module regulating embryo development

Assembly of complete genomes can reveal functional genetic elements missing from draft sequences. Here we present the near-complete telomere-to-telomere and contiguous genome of the cotton species Gossypium raimondii. Our assembly identified gaps and misoriented or misassembled regions in previous assemblies and produced 13 centromeres, with 25 chromosomal ends having telomeres. In contrast to satellite-rich Arabidopsis and rice centromeres, cotton centromeres lack phased CENH3 nucleosome positioning patterns and probably evolved by invasion from long terminal repeat retrotransposons. In-depth expression profiling of transposable elements revealed a previously unannotated DNA transposon (MuTC01) that interacts with miR2947 to produce trans-acting small interfering RNAs (siRNAs), one of which targets the newly evolved LEC2 (LEC2b) to produce phased siRNAs. Systematic genome editing experiments revealed that this tripartite module, miR2947-MuTC01-LEC2b, controls the morphogenesis of complex folded embryos characteristic of Gossypium and its close relatives in the cotton tribe. Our study reveals a trans-acting siRNA-based tripartite regulatory pathway for embryo development in higher plants.

The LOB domain protein, a novel transcription factor with multiple functions: A review

The LATERAL ORGAN BOUNDARIES DOMAIN (LBD) protein, named for its LATERAL ORGAN BOUNDARIES (LOB) domain, is a member of a class of specific transcription factors commonly found in plants and is absent from all other groups of organisms. LBD TFs have been systematically identified in about 35 plant species and are involved in regulating various aspects of plant growth and development. However, research on the signaling network and regulatory functions of LBD TFs is insufficient, and only a few members have been studied. Moreover, a comprehensive review of these existing studies is lacking. In this review, the structure, regulatory mechanism and function of LBD TFs in recent years were reviewed in order to better understand the role of LBD TFs in plant growth and development, and to provide a new perspective for the follow-up study of LBD TFs.

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.

Cotton BLH1 and KNOX6 antagonistically modulate fiber elongation via regulation of linolenic acid biosynthesis

BEL1-LIKE HOMEODOMAIN (BLH) proteins are known to function in various plant developmental processes. However, the role of BLHs in regulating plant cell elongation is still unknown. Here, we identify a BLH gene, GhBLH1, that positively regulates fiber cell elongation. Combined transcriptomic and biochemical analyses reveal that GhBLH1 enhances linolenic acid accumulation to promote cotton fiber cell elongation by activating the transcription of GhFAD7A-1 via binding of the POX domain of GhBLH1 to the TGGA cis-element in the GhFAD7A-1 promoter. Knockout of GhFAD7A-1 in cotton significantly reduces fiber length, whereas overexpression of GhFAD7A-1 results in longer fibers. The K2 domain of GhKNOX6 directly interacts with the POX domain of GhBLH1 to form a functional heterodimer, which interferes with the transcriptional activation of GhFAD7A-1 via the POX domain of GhBLH1. Overexpression of GhKNOX6 leads to a significant reduction in cotton fiber length, whereas knockout of GhKNOX6 results in longer cotton fibers. An examination of the hybrid progeny of GhBLH1 and GhKNOX6 transgenic cotton lines provides evidence that GhKNOX6 negatively regulates GhBLH1-mediated cotton fiber elongation. Our results show that the interplay between GhBLH1 and GhKNOX6 modulates regulation of linolenic acid synthesis and thus contributes to plant cell elongation.

Expression profile analysis of cotton fiber secondary cell wall thickening stage

To determine the genes associated with the fiber strength trait in cotton, three different cotton cultivars were selected: Sea Island cotton (Xinhai 32, with hyper-long fibers labeled as HL), and upland cotton (17-24, with long fibers labeled as L, and 62-33, with short fibers labeled as S). These cultivars were chosen to assess fiber samples with varying qualities. RNA-seq technology was used to analyze the expression profiles of cotton fibers at the secondary cell wall (SCW) thickening stage (20, 25, and 30 days post-anthesis (DPA)). The results showed that a large number of differentially expressed genes (DEGs) were obtained from the three assessed cotton cultivars at different stages of SCW development. For instance, at 20 DPA, Sea Island cotton (HL) had 6,215 and 5,364 DEGs compared to upland cotton 17-24 (L) and 62-33 (S), respectively. Meanwhile, there were 1,236 DEGs between two upland cotton cultivars, 17-24 (L) and 62-33 (S). Gene Ontology (GO) term enrichment identified 42 functions, including 20 biological processes, 11 cellular components, and 11 molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified several pathways involved in SCW synthesis and thickening, such as glycolysis/gluconeogenesis, galactose metabolism, propanoate metabolism, biosynthesis of unsaturated fatty acids pathway, valine, leucine and isoleucine degradation, fatty acid elongation pathways, and plant hormone signal transduction. Through the identification of shared DEGs, 46 DEGs were found to exhibit considerable expressional differences at different fiber stages from the three cotton cultivars. These shared DEGs have functions including REDOX enzymes, binding proteins, hydrolases (such as GDSL thioesterase), transferases, metalloproteins (cytochromatin-like genes), kinases, carbohydrates, and transcription factors (MYB and WRKY). Therefore, RT-qPCR was performed to verify the expression levels of nine of the 46 identified DEGs, an approach which demonstrated the reliability of RNA-seq data. Our results provided valuable molecular resources for clarifying the cell biology of SCW biosynthesis during fiber development in cotton.

Unraveling the lncRNA-miRNA-mRNA Regulatory Network Involved in Poplar Coma Development through High-Throughput Sequencing

Poplar coma, the fluff-like appendages of seeds originating from the differentiated surface cells of the placenta and funicle, aids in the long-distance dispersal of seeds in the spring. However, it also poses hazards to human safety and causes pollution in the surrounding environment. Unraveling the regulatory mechanisms governing the initiation and development of coma is essential for addressing this issue comprehensively. In this study, strand-specific RNA-seq was conducted at three distinct stages of coma development, revealing 1888 lncRNAs and 52,810 mRNAs. The expression profiles of lncRNAs and mRNAs during coma development were analyzed. Subsequently, potential target genes of lncRNAs were predicted through co-localization and co-expression analyses. Integrating various types of sequencing data, lncRNA-miRNA-TF regulatory networks related to the initiation of coma were constructed. Utilizing identified differentially expressed genes encoding kinesin and actin, lncRNA-miRNA-mRNA regulatory networks associated with the construction and arrangement of the coma cytoskeleton were established. Additionally, relying on differentially expressed genes encoding cellulose synthase, sucrose synthase, and expansin, lncRNA-miRNA-mRNA regulatory networks related to coma cell wall synthesis and remodeling were developed. This study not only enhances the comprehension of lncRNA but also provides novel insights into the molecular mechanisms governing the initiation and development of poplar coma.

Combined genome and transcriptome analysis of elite fiber quality in Gossypium barbadense

Gossypium barbadense, which is one of several species of cotton, is well known for its superior fiber quality. However, the genetic basis of its high-quality fiber remains largely unexplored. Here, we resequenced 269 G. barbadense accessions. Phylogenetic structure analysis showed that the set of accessions was clustered into 3 groups: G1 and G2 mainly included modern cultivars from Xinjiang, China, and G3 was related to widely introduced accessions in different regions worldwide. A genome-wide association study of 5 fiber quality traits across multiple field environments identified a total of 512 qtls (main-effect QTLs) and 94 qtlEs (QTL-by-environment interactions) related to fiber quality, of which 292 qtls and 57 qtlEs colocated with previous studies. We extracted the genes located in these loci and performed expression comparison, local association analysis, and introgression segment identification. The results showed that high expression of hormone-related genes during fiber development, introgressions from Gossypium hirsutum, and the recombination of domesticated elite allelic variation were 3 major contributors to improve the fiber quality of G. barbadense. In total, 839 candidate genes with encoding region variations associated with elite fiber quality were mined. We confirmed that haplotype GB_D03G0092H traced to G. hirsutum introgression, with a 1-bp deletion leading to a frameshift mutation compared with GB_D03G0092B, significantly improved fiber quality. GB_D03G0092H is localized in the plasma membrane, while GB_D03G0092B is in both the nucleus and plasma membrane. Overexpression of GB_D03G0092H in Arabidopsis (Arabidopsis thaliana) significantly improved the elongation of longitudinal cells. Our study systematically reveals the genetic basis of the superior fiber quality of G. barbadense and provides elite segments and gene resources for breeding high-quality cotton cultivars.

GhWER controls fiber initiation and early elongation by regulating ethylene signaling pathway in cotton (Gossypium hirsutum)

Cotton fibers are specialized single-cell trichomes derived from epidermal cells, similar to root hairs and trichomes in Arabidopsis. While the MYB-bHLH-WD40 (MBW) complex has been shown to regulate initiation of both root hairs and trichomes in Arabidopsis, the role of their homologous gene in cotton fiber initiation remains unknown. In this study, we identified a R2R3 MYB transcription factor (TF), GhWER, which exhibited a significant increase in expression within the outer integument of ovule at -1.5 DPA (days post anthesis). Its expression peaked at -1 DPA and then gradually decreased. Knockout of GhWER using CRISPR technology inhibited the initiation and early elongation of fiber initials, resulting in the shorter mature fiber length. Additionally, GhWER interacted with two bHLH TF, GhDEL65 and GhbHLH121, suggesting a potential regulatory complex for fiber development. RNA-seq analysis of the outer integument of the ovule at -1.5 DPA revealed that the signal transduction pathways of ethylene, auxin and gibberellin were affected in the GhWER knockout lines. Further examination demonstrated that GhWER directly activated ethylene signaling genes, including ACS1 and ETR2. These findings highlighted the biological function of GhWER in regulating cotton fiber initiation and early elongation, which has practical significance for improving fiber quality and yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01477-6.

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.

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.

The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum

KEY MESSAGE

GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

In vitro susceptibility of nontuberculous mycobacteria in China

OBJECTIVES

This study aimed to measure the prevalence of resistance to antimicrobial agents, and explore the risk factors associated with drug resistance by using nontuberculous Mycobacteria (NTM) isolates from China.

METHODS

A total of 335 NTM isolates were included in our analysis. Broth dilution method was used to determine in vitro drug susceptibility of NTM isolates.

RESULTS

Clarithromycin (CLA) was the most potent drug for Mycobacterium intracellulare (MI). The resistance rate of 244 MI isolates to CLA was 21%, yielding a minimum inhibitory concentrations (MIC)50 and MIC90 of 8 and 64 mg/L, respectively. 51% of 244 MI isolates exhibited resistance to amikacin (AMK). For 91 Mycobacterium abscessus complex (MABC) isolates, 6 (7%) and 49 (54%) isolates were categorized as resistant to CLA at day 3 and 14, respectively. The resistance rate to CLA for Mycobacterium abscessus subspecies abscessus (MAA) was dramatically higher than that for Mycobacterium abscessus subspecies massiliense (MAM). Additionally, the percentage of patients presenting fever in the CLA-susceptible group was significantly higher than that in the CLA-resistant group.

CONCLUSIONS

Our data demonstrate that approximate one fifth of MI isolates are resistant to CLA. We have identified a higher proportion of CLA-resistant MAA isolates than MAM. The patients caused by CLA-resistant MI are at low risk for presenting with fever relative to CLA-susceptible group.

Genome-Wide Identification and Expression Analysis of 1-Aminocyclopropane-1-Carboxylate Synthase (ACS) Gene Family in Chenopodium quinoa

Ethylene plays an important role in plant development and stress resistance. The rate-limiting enzyme in ethylene biosynthesis is 1-aminocyclopropane-1-carboxylic acid synthase (ACS). C. quinoa (Chenopodium quinoa) is an important food crop known for its strong tolerance to abiotic stresses. However, knowledge regarding the ACS gene family in C. quinoa remains restricted. In this study, we successfully identified 12 ACS genes (CqACSs) from the C. quinoa genome. Through thorough analysis of their sequences and phylogenetic relationships, it was verified that 8 out of these 12 CqACS isozymes exhibited substantial resemblance to ACS isozymes possessing ACS activity. Furthermore, these eight isozymes could be categorized into three distinct groups. The four remaining CqACS genes grouped under category IV displayed notable similarities with AtACS10 and AtACS12, known as amido transferases lacking ACS activity. The CqACS proteins bore resemblance to the AtACS proteins and had the characteristic structural features typically observed in plant ACS enzymes. Twelve CqACS genes were distributed across 8 out of the 18 chromosomes of C. quinoa. The CqACS genes were expanded from segment duplication. Many cis-regulatory elements related with various abiotic stresses, phytohormones, and light were found. The expression patterns of ACS genes varied across different tissues of C. quinoa. Furthermore, the analysis of gene expression patterns under abiotic stress showed that CqACS genes can be responsive to various stresses, implying their potential functions in adapting to various abiotic stresses. The findings from this research serve as a foundation for delving deeper into the functional roles of CqACS genes.

GhXB38D represses cotton fibre elongation through ubiquitination of ethylene biosynthesis enzymes GhACS4 and GhACO1

Ethylene plays an essential role in the development of cotton fibres. Ethylene biosynthesis in plants is elaborately regulated by the activities of key enzymes, 1-aminocyclopropane-1-carboxylate oxidase (ACO) and 1-aminocyclopropane-1-carboxylate synthase (ACS); however, the potential mechanism of post-translational modification of ACO and ACS to control ethylene synthesis in cotton fibres remains unclear. Here, we identify an E3 ubiquitin ligase, GhXB38D, that regulates ethylene biosynthesis during fibre elongation in cotton. GhXB38D gene is highly expressed in cotton fibres during the rapid elongation stage. Suppressing GhXB38D expression in cotton significantly enhanced fibre elongation and length, accompanied by the up-regulation of genes associated with ethylene signalling and fibre elongation. We demonstrated that GhXB38D interacts with the ethylene biosynthesis enzymes GhACS4 and GhACO1 in elongating fibres and specifically mediates their ubiquitination and degradation. The inhibition of GhXB38D gene expression increased the stability of GhACS4 and GhACO1 proteins in cotton fibres and ovules, resulting in an elevated concentration of ethylene. Our findings highlight the role of GhXB38D as a regulator of ethylene synthesis by ubiquitinating ACS4 and ACO1 proteins and modulating their stability. GhXB38D acts as a negative regulator of fibre elongation and serves as a potential target for enhancing cotton fibre yield and quality through gene editing strategy.

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.

Genome-wide association study of fiber yield-related traits uncovers the novel genomic regions and candidate genes in Indian upland cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is a major fiber crop that is cultivated worldwide and has significant economic importance. India harbors the largest area for cotton cultivation, but its fiber yield is still compromised and ranks 22nd in terms of productivity. Genetic improvement of cotton fiber yield traits is one of the major goals of cotton breeding, but the understanding of the genetic architecture underlying cotton fiber yield traits remains limited and unclear. To better decipher the genetic variation associated with fiber yield traits, we conducted a comprehensive genome-wide association mapping study using 117 Indian cotton germplasm for six yield-related traits. To accomplish this, we generated 2,41,086 high-quality single nucleotide polymorphism (SNP) markers using genotyping-by-sequencing (GBS) methods. Population structure, PCA, kinship, and phylogenetic analyses divided the germplasm into two sub-populations, showing weak relatedness among the germplasms. Through association analysis, 205 SNPs and 134 QTLs were identified to be significantly associated with the six fiber yield traits. In total, 39 novel QTLs were identified in the current study, whereas 95 QTLs overlapped with existing public domain data in a comparative analysis. Eight QTLs, qGhBN_SCY_D6-1, qGhBN_SCY_D6-2, qGhBN_SCY_D6-3, qGhSI_LI_A5, qGhLI_SI_A13, qGhLI_SI_D9, qGhBW_SCY_A10, and qGhLP_BN_A8 were identified. Gene annotation of these fiber yield QTLs revealed 2,509 unique genes. These genes were predominantly enriched for different biological processes, such as plant cell wall synthesis, nutrient metabolism, and vegetative growth development in the gene ontology (GO) enrichment study. Furthermore, gene expression analysis using RNAseq data from 12 diverse cotton tissues identified 40 candidate genes (23 stable and 17 novel genes) to be transcriptionally active in different stages of fiber, ovule, and seed development. These findings have revealed a rich tapestry of genetic elements, including SNPs, QTLs, and candidate genes, and may have a high potential for improving fiber yield in future breeding programs for Indian cotton.

Activation of Gossypium hirsutum ACS6 Facilitates Fiber Development by Improving Sucrose Metabolism and Transport

Cotton fiber yield depends on the density of fiber cell initials that form on the ovule epidermis. Fiber initiation is triggered by MYB-MIXTA-like transcription factors (GhMMLs) and requires a sucrose supply. Ethylene or its precursor ACC (1-aminocyclopropane-1-carboxylic acid) is suggested to affect fiber yield. The Gossypium hirsutum (L.) genome contains 35 ACS genes (GhACS) encoding ACC synthases. Here, we explored the role of a GhACS family member in the regulation of fiber initiation. Expression analyses showed that the GhACS6.3 gene pair was specifically expressed in the ovules during fiber initiation (3 days before anthesis to 5 days post anthesis, -3 to 5 DPA), especially at -3 DPA, whereas other GhACS genes were expressed at very low or undetectable levels. The expression profile of GhACS6.3 during fiber initial development was confirmed by qRT-PCR analysis. Transgenic lines overexpressing GhACS6.3 (GhACS6.3-OE) showed increased ACC accumulation in ovules, which promoted the formation of fiber initials and fiber yield components. This was accompanied by increased transcript levels of GhMML3 and increased transcript levels of genes encoding sucrose transporters and sucrose synthase. These findings imply that GhACS6.3 activation is required for fiber initial development. Our results lay the foundation for further research on increasing cotton fiber production.

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.

Single-cell RNA landscape of the special fiber initiation process in Bombax ceiba

As a new source of natural fibers, the Bombax ceiba tree can provide thin, light, extremely soft and warm fiber material for the textile industry. Natural fibers are an ideal model system for studying cell growth and differentiation, but the molecular mechanisms that regulate fiber initiation are not fully understood. In B. ceiba, we found that fiber cells differentiate from the epidermis of the inner ovary wall. Each initiated cell then divides into a cluster of fiber cells that eventually develop into mature fibers, a process very different from the classical fiber initiation process of cotton. We used high-throughput single-cell RNA sequencing (scRNA-seq) to examine the special characteristics of fiber initiation in B. ceiba. A total of 15 567 high-quality cells were identified from the inner wall of the B. ceiba ovary, and 347 potential marker genes for fiber initiation cell types were identified. Two major cell types, initiated fiber cells and epidermal cells, were identified and verified by RNA in situ hybridization. A developmental trajectory analysis was used to reconstruct the process of fiber cell differentiation in B. ceiba. Comparative analysis of scRNA-seq data from B. ceiba and cotton (Gossypium hirsutum) confirmed that the additional cell division process in B. ceiba is a novel species-specific mechanism for fiber cell development. Candidate genes and key regulators that may contribute to fiber cell differentiation and division in B. ceiba were identified. This work reveals gene expression signatures during B. ceiba fiber initiation at a single-cell resolution, providing a new strategy and viewpoint for investigation of natural fiber cell differentiation and development.

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.

Cell cycle-dependent kinase inhibitor GhKRP6, a direct target of GhBES1.4, participates in BR regulation of cell expansion in cotton

The steroidal hormone brassinosteroid (BR) has been shown to positively regulate cell expansion in plants. However, the specific mechanism by which BR controls this process has not been fully understood. In this study, RNA-seq and DAP-seq analysis of GhBES1.4 (a core transcription factor in BR signaling) were used to identify a cotton cell cycle-dependent kinase inhibitor called GhKRP6. The study found that GhKRP6 was significantly induced by the BR hormone and that GhBES1.4 directly promoted the expression of GhKRP6 by binding to the CACGTG motif in its promoter region. GhKRP6-silenced cotton plants had smaller leaves with more cells and reduced cell size. Furthermore, endoreduplication was inhibited, which affected cell expansion and ultimately decreased fiber length and seed size in GhKRP6-silenced plants compared with the control. The KEGG enrichment results of control and VIGS-GhKRP6 plants revealed differential expression of genes related to cell wall biosynthesis, MAPK, and plant hormone transduction pathways - all of which are related to cell expansion. Additionally, some cyclin-dependent kinase (CDK) genes were upregulated in the plants with silenced GhKRP6. Our study also found that GhKRP6 could interact directly with a cell cycle-dependent kinase called GhCDKG. Taken together, these results suggest that BR signaling influences cell expansion by directly modulating the expression of cell cycle-dependent kinase inhibitor GhKRP6 via GhBES1.4.

Analysis of Pulmonary Function in Thymoma Subjects: A 20-Year Retrospective Cohort Study

BACKGROUND

Thymoma is the most common tumor of the anterior mediastinum. However, the correlation between thymoma stage and pulmonary function was not assessed. Our objective in this study was to describe the pulmonary function in thymoma subjects stratified with different staging systems.

METHODS

A total of 143 subjects with a diagnosis of thymoma who underwent extended thymectomy for thymoma between January 2001 and December 2019 were reviewed retrospectively. All the subjects experienced pulmonary function tests (PFTs) using Master Screen PFT system and total respiratory resistance measurement.

RESULTS

We evaluated 143 subjects with a diagnosis of thymoma; the significant differences were observed in mean values of vital capacity, inspiratory volume (IC), total lung capacity (TLC), ratio of residual volume to total lung capacity (RV/TLC), forced vital capacity, forced expiratory volume in 1 second, ratio of forced expiratory volume in 1 second to forced vital capacity, peak expiratory flow, peak inspiratory flow, maximum ventilation volume, total airway resistance, and diffusing capacity for carbon monoxide (DLCO) across upper airway obstruction classification. PFTs of subjects with varying Masaoka stages are different. RV and RV/TLC of subjects in stages III and IV were higher than those of normal level, while DLCO of subjects in stage IV was lower than the normal level, and the mean level of IC showed significant difference between stage II and stage III.

DISCUSSION

The pulmonary function patterns of thymoma subjects significantly correlate with tumor location and size rather than clinical Masaoka stage.

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.

Comparison of the Diagnostic Performance of MeltPro and Next-Generation Sequencing in Determining Fluoroquinolone Resistance in Multidrug-Resistant Tuberculosis Isolates

This study systematically investigated the performance of MeltPro and next-generation sequencing in the diagnosis of fluoroquinolone (FQ) resistance among multidrug-resistant tuberculosis patients and explored the relationship between nucleotide alteration and the level of phenotypic susceptibility to FQs. From March 2019 to June 2020, a feasibility and validation study with both MeltPro and next-generation sequencing was performed in 126 patients with multidrug-resistant tuberculosis. Using phenotypic drug susceptibility testing as the gold standard, 95.3% (82 of 86) of ofloxacin-resistant isolates were identified correctly by MeltPro. In addition, whole-genome sequencing was able to detect 83 phenotypically ofloxacin-resistant isolates. The isolates with an individual gyrB mutation outside the quinolone resistance-determining region (QRDR) had minimum inhibitory concentrations (MICs) of ≤2 μg/mL. Despite showing low MICs close to the breakpoint for isolates carrying only gyrA_Ala90Val, the combined mutation gyrB_Asp461Asn caused the ofloxacin MIC to be eight higher than that obtained in Mycobacterium tuberculosis (MTB) isolates with the Ala90Val mutation alone (median, 32 μg/mL; P = 0.038). Heteroresistance was observed in 12 of 88 isolates harboring mutations in the QRDRs. In conclusion, our data show that MeltPro and the whole-genome sequencing assay correctly can identify FQ resistance caused by mutations in the gyrA QRDR. The combined gyrB_Asp461Asn mutation may significantly decrease in vitro FQ susceptibility of MTB isolates with low-level-resistance-associated gyrA mutations.

Brassinosteroids regulate cotton fiber elongation by modulating very-long-chain fatty acid biosynthesis

Brassinosteroid (BR), a growth-promoting phytohormone, regulates many plant growth processes including cell development. However, the mechanism by which BR regulates fiber growth is poorly understood. Cotton (Gossypium hirsutum) fibers are an ideal single-cell model in which to study cell elongation due to their length. Here we report that BR controls cotton fiber elongation by modulating very-long-chain fatty acid (VLCFA) biosynthesis. BR deficiency reduces the expression of 3-ketoacyl-CoA synthases (GhKCSs), the rate-limiting enzymes involved in VLCFA biosynthesis, leading to lower saturated VLCFA contents in pagoda1 (pag1) mutant fibers. In vitro ovule culture experiments show that BR acts upstream of VLCFAs. Silencing of BRI1-EMS-SUPPRESOR 1.4 (GhBES1.4), encoding a master transcription factor of the BR signaling pathway, significantly reduces fiber length, whereas GhBES1.4 overexpression produces longer fibers. GhBES1.4 regulates endogenous VLCFA contents and directly binds to BR RESPONSE ELEMENTS (BRREs) in the GhKCS10_At promoter region, which in turn regulates GhKCS10_At expression to increase endogenous VLCFA contents. GhKCS10_At overexpression promotes cotton fiber elongation, whereas GhKCS10_At silencing inhibits cotton fiber growth, supporting a positive regulatory role for GhKCS10_At in fiber elongation. Overall, these results uncover a mechanism of fiber elongation through crosstalk between BR and VLCFAs at the single-cell level.

A brassinosteroid transcriptional regulatory network participates in regulating fiber elongation in cotton

Brassinosteroids (BRs) participate in the regulation of plant growth and development through BRI1-EMS-SUPPRESSOR1 (BES1)/BRASSINAZOLE-RESISTANT1 (BZR1) family transcription factors. Cotton (Gossypium hirsutum) fibers are highly elongated single cells, and BRs play a vital role in the regulation of fiber elongation. However, the mode of action on how BR is involved in the regulation of cotton fiber elongation remains unexplored. Here, we generated GhBES1.4 over expression lines and found that overexpression of GhBES1.4 promoted fiber elongation, whereas silencing of GhBES1.4 reduced fiber length. DNA affinity purification and sequencing (DAP-seq) identified 1,531 target genes of GhBES1.4, and five recognition motifs of GhBES1.4 were identified by enrichment analysis. Combined analysis of DAP-seq and RNA-seq data of GhBES1.4-OE/RNAi provided mechanistic insights into GhBES1.4-mediated regulation of cotton fiber development. Further, with the integrated approach of GWAS, RNA-seq, and DAP-seq, we identified seven genes related to fiber elongation that were directly regulated by GhBES1.4. Of them, we showed Cytochrome P450 84A1 (GhCYP84A1) and 3-hydroxy-3-methylglutaryl-coenzyme A reductase 1 (GhHMG1) promote cotton fiber elongation. Overall, the present study established the role of GhBES1.4-mediated gene regulation and laid the foundation for further understanding the mechanism of BR participation in regulating fiber 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.

Genome scale analysis of 1-aminocyclopropane-1-carboxylate oxidase gene family in G. barbadense and its functions in cotton fiber development

A class of proteins, 1-aminocyclopropane-1-carboxylate oxidase (ACO), is required in the final step of production of ethylene from its immediate precursor 1-aminocyclopropane-1-carboxylic acid (ACC). Despite the crucial and regulatory role of ACO gene family in the fiber development, it has not been thoroughly analyzed and annotated in G. barbadense genome. In the present study, we have identified and characterized all isoforms of ACO gene family from genomes of Gossypium arboreum, G. barbadense, G. hirsutum and G. raimondii. Phylogenetic analysis classified all ACO proteins into six distinct groups on the basis of maximum likelihood. Gene locus analysis and circos plots indicated the distribution and relationship of these genes in cotton genomes. Transcriptional profiling of ACO isoforms in G. arboreum, G. barbadense and G. hirsutum fiber development exhibited the highest expression in G. barbadense during early fiber elongation. Moreover, the accumulation of ACC was found highest in developing fibers of G. barbadense in comparison with other cotton species. ACO expression and ACC accumulation correlated with the fiber length in cotton species. Addition of ACC to the ovule cultures of G. barbadense significantly increased fiber elongation while ethylene inhibitors hindered fiber elongation. These findings will be helpful in dissecting the role of ACOs in cotton fiber development and pave a way towards genetic manipulations for fiber quality improvement.

Recent progression and future perspectives in cotton genomic breeding

Upland cotton is an important global cash crop for its long seed fibers and high edible oil and protein content. Progress in cotton genomics promotes the advancement of cotton genetics, evolutionary studies, functional genetics, and breeding, and has ushered cotton research and breeding into a new era. Here, we summarize high-impact genomics studies for cotton from the last 10 years. The diploid Gossypium arboreum and allotetraploid Gossypium hirsutum are the main focus of most genetic and genomic studies. We next review recent progress in cotton molecular biology and genetics, which builds on cotton genome sequencing efforts, population studies, and functional genomics, to provide insights into the mechanisms shaping abiotic and biotic stress tolerance, plant architecture, seed oil content, and fiber development. We also suggest the application of novel technologies and strategies to facilitate genome-based crop breeding. Explosive growth in the amount of novel genomic data, identified genes, gene modules, and pathways is now enabling researchers to utilize multidisciplinary genomics-enabled breeding strategies to cultivate "super cotton", synergistically improving multiple traits. These strategies must rise to meet urgent demands for a sustainable cotton industry.

Postharvest chilling diminishes melon flavor via effects on volatile acetate ester biosynthesis

In postharvest handling systems, refrigeration can extend fruit shelf life and delay decay via slowing ripening progress; however, it selectively alters the biosynthesis of flavor-associated volatile organic compounds (VOCs), which results in reduced flavor quality. Volatile esters are major contributors to melon fruit flavor. The more esters, the more consumers enjoy the melon fruit. However, the effects of chilling on melon flavor and volatiles associated with consumer liking are yet to be fully understood. In the present study, consumer sensory evaluation showed that chilling changed the perception of melon fruit. Total ester content was lower after chilling, particularly volatile acetate esters (VAEs). Transcriptomic analysis revealed that transcript abundance of multiple flavor-associated genes in fatty acid and amino acid pathways was reduced after chilling. Additionally, expression levels of the transcription factors (TFs), such as NOR, MYB, and AP2/ERF, also were substantially downregulated, which likely altered the transcript levels of ester-associated pathway genes during cold storage. VAE content and expression of some key genes recover after transfer to room temperature. Therefore, chilling-induced changes of VAE profiles were consistent with expression patterns of some pathway genes that encode specific fatty acid- and amino acid-mobilizing enzymes as well as TFs involved in fruit ripening, metabolic regulation, and hormone signaling.

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.

Attenuated Cytokine-Induced Memory-Like Natural Killer Cell Responses to Mycobacterium tuberculosis in Tuberculosis Patients

Purpose

This study aimed to investigate the phenotype, proliferation and functional alterations of cytokine-induced memory-like natural killer (CIML NK) cells from healthy subjects and TB patients, and assessed the efficacy of CIML NK cells in response to H37Rv-infected U937 cells in vitro.

Methods

Fresh peripheral blood mononuclear cells (PBMCs) were isolated from healthy people and tuberculosis patients and activated for 16h using low-dose IL-15, or IL-12, IL-15, IL-18 combination or IL-12, IL-15, IL-18 and MTB H37Rv lysates, respectively, followed by low-dose IL-15 maintenance for another 7 days. Then, the PBMCs were co-cultured with K562 and H37Rv-infected U937, and the purified NK cells were co-cultured with H37Rv infected U937. The phenotype, proliferation and response function of CIML NK cells were assessed using flow cytometry. Finally, colony forming units were enumerated to confirm the survival of intracellular MTB.

Results

CIML NK phenotypes from TB patients were similar to healthy controls. CIML NK cells undergo higher rates of proliferation after IL-12/15/18 pre-activation. Moreover, the poor expansion potential of CIML NK cells co-stimulated with MTB lysates. CIML NK cells from healthy individuals showed enhanced IFN-γ functional to H37Rv infected U937 cells, along with significantly enhanced killing of H37Rv. However, the CIML NK cells from TB patients show attenuated IFN-γ production and now enhanced the ability of killing intracellular MTB compared to those from healthy donors after co-cultured with H37Rv infected U937.

Conclusion

CIML NK cells from healthy individuals exist the increased ability of IFN-γ secretion and boosted anti-MTB activity in vitro, which from TB patients show impaired IFN-γ production and no enhanced anti-MTB activity compared to those from healthy donors. Additionally, we observe the poor expansion potential of CIML NK cells co-stimulated with antigens from MTB. These results open up new possibilities for NK cell-based anti-tuberculosis immunotherapeutic strategies.

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.

The physiological and transcriptomic study of secondary growth in Neolamarckia cadamba stimulated by the ethylene precursor ACC

Though many biological roles of ethylene have been investigated intensively, the molecular mechanism of ethylene's action in woody plants remains unclear. In this study, we investigated the effects of exogenous 1-aminocyclopropane-1-carboxylic acid (ACC), the precursor of ethylene, on the growth of Neolamarckia cadamba seedlings, a fast-growing tropical tree. After 14 days of ACC treatment, the plants showed a reduced physiological morphology while stem diameter increased; however, this did not occur after the addition of 1-MCP. Meanwhile, the lignin content of N. cadamba also increased. Transcriptome analysis revealed that the expression of the ethylene biosynthesis and signaling genes ACC oxidase (ACO) and ethylene insensitive 3 (EIN3) were up-regulated mainly at the 6th hour and the 3rd day of the ACC treatment, respectively. The transcription levels of transcription factors, mainly in the basic helix-loop-helix (bHLH), ethylene response factor (ERF), WRKY and v-myb avian myeloblastosis viral oncogene homolog (MYB) families, involved in the ethylene signaling and secondary growth also increased significantly. Furthermore, in accordance to the increased lignification of the stem, the transcriptional level of key enzymes in the phenylalanine pathway were elevated after the ACC treatment. Our results revealed the physiological and molecular mechanisms underlying the secondary growth stimulated by exogenous ACC treatment on N. cadamba seedlings.

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.

A very-long-chain fatty acid synthesis gene, SD38, influences plant height by activating ethylene biosynthesis in rice

As an important trait in crop breeding, plant height is associated with lodging resistance and yield. With the identification and cloning of several semi-dwarfing genes, increasing numbers of semi-dwarf cultivars have emerged, which has led to a 'green revolution' in rice (Oryza sativa) production. In this study, we identified a rice semi-dwarf mutant, semi-dwarf 38 (sd38), which showed significantly reduced cell length. SD38 encodes a fatty acid elongase, β-ketoacyl-CoA synthase, which is involved in the synthesis of very-long-chain fatty acids (VLCFAs). Expression analysis showed that SD38 was localized on the membrane of the endoplasmic reticulum, and was expressed in all analyzed tissues with differential abundance. The mutation of SD38 affected lipid metabolism in the sd38 mutant. A functional complementarity test in Saccharomyces cerevisiae indicated that SD38 was capable of complementing the deficiency of ELO3p activity in BY4741-elo3 knockout yeast cells by participating in the synthesis of C24:0 VLCFA. Significant changes were observed in the expression of genes involved in ethylene synthesis, which resulted in reduced content of the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) in the sd38 mutant. Exogenously supplied VLCFA (C24:0) increased the expression levels of OsACS3, OsACS4, and OsACO7 and the plant height of sd38 mutant seedlings, similar to the effect of exogenous application of ACC and ethephon. These results reveal a relationship among VLCFAs, ethylene biosynthesis, and plant height and improve our understanding of plant height development in crops.

Status and prospects of genome-wide association studies in cotton

Over the last two decades, the use of high-density SNP arrays and DNA sequencing have allowed scientists to uncover the majority of the genotypic space for various crops, including cotton. Genome-wide association study (GWAS) links the dots between a phenotype and its underlying genetics across the genomes of populations. It was first developed and applied in the field of human disease genetics. Many areas of crop research have incorporated GWAS in plants and considerable literature has been published in the recent decade. Here we will provide a comprehensive review of GWAS studies in cotton crop, which includes case studies on biotic resistance, abiotic tolerance, fiber yield and quality traits, current status, prospects, bottlenecks of GWAS and finally, thought-provoking question. This review will serve as a catalog of GWAS in cotton and suggest new frontiers of the cotton crop to be studied with this important tool.

GhBZR3 suppresses cotton fiber elongation by inhibiting very-long-chain fatty acid biosynthesis

The BRASSINAZOLE-RESISTANT (BZR) transcription factor is a core component of brassinosteroid (BR) signaling and is involved in the development of many plant species. BR is essential for the initiation and elongation of cotton fibers. However, the mechanism of BR-regulating fiber development and the function of BZR is poorly understood in Gossypium hirsutum L. (cotton). Here, we identified a BZR family transcription factor protein referred to as GhBZR3 in cotton. Overexpression of GhBZR3 in Arabidopsis caused shorter root hair length, hypocotyl length, and hypocotyl cell length, indicating that GhBZR3 negatively regulates cell elongation. Pathway enrichment analysis from VIGS-GhBZR3 cotton plants found that fatty acid metabolism and degradation might be the regulatory pathway that is primarily controlled by GhBZR3. Silencing GhBZR3 expression in cotton resulted in taller plant height as well as longer fibers. The very-long-chain fatty acid (VLCFA) content was also significantly increased in silenced GhBZR3 plants compared with the wild type. The GhKCS13 promoter, a key gene for VLCFA biosynthesis, contains two GhBZR3 binding sites. The results of yeast one-hybrid, electrophoretic mobility shift, and luciferase assays revealed that GhBZR3 directly interacted with the GhKCS13 promoter to suppress gene expression. Taken together, these results indicate that GhBZR3 negatively regulates cotton fiber development by reducing VLCFA biosynthesis. This study not only deepens our understanding of GhBZR3 function in cotton fiber development, but also highlights the potential of improving cotton fiber length and plant growth using GhBZR3 and its related genes in future cotton breeding programs.