High salinity is one of the main factors limiting cotton growth and productivity. The genes that regulate salt stress in TM-1 upland cotton were monitored using microarray and real-time PCR (RT-PCR) with samples taken from roots. Microarray analysis showed that 1503 probe sets were up-regulated and 1490 probe sets were down-regulated in plants exposed for 3h to 100mM NaCl, and RT-PCR analysis validated 42 relevant/related genes. The distribution of enriched gene ontology terms showed such important processes as the response to water stress and pathways of hormone metabolism and signal transduction were induced by the NaCl treatment. Some key regulatory gene families involved in abiotic and biotic sources of stress such as WRKY, ERF, and JAZ were differentially expressed. Our transcriptome analysis might provide some useful insights into salt-mediated signal transduction pathways in cotton and offer a number of candidate genes as potential markers of tolerance to salt stress.

Polygalacturonase-inhibiting protein (PGIP), belonging to a group of plant defence proteins, specifically inhibits endopolygalacturonases secreted by pathogens. Herein, we showed that purified GhPGIP1 is a functional inhibitor of Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, the two fungal pathogens causing cotton wilt. Transcription of GhPGIP1 was increased in cotton upon infection, wounding, and treatment with defence hormone and H2O2. Resistance by GhPGIP1 was examined by its virus-induced gene silencing in cotton and overexpression in Arabidopsis. GhPGIP1-silenced cotton was highly susceptible to the infections. GhPGIP1 overexpression in transgenic Arabidopsis conferred resistance to the infection, accompanied by enhanced expression of pathogenesis-related proteins (PRs), isochorismate synthase 1 (ICS1), enhanced disease susceptibility 1 (EDS1), and phytoalexin-deficient 4 (PAD4) genes. Transmission electron microscopy revealed cell wall alteration and cell disintegration in plants inoculated with polygalacturonase (PGs), implying its role in damaging the cell wall. Docking studies showed that GhPGIP1 interacted strongly with C-terminal of V. dahliae PG1 (VdPG1) beyond the active site but weakly interacted with C-terminal of F. oxysporum f. sp. vasinfectum (FovPG1). These findings will contribute towards the understanding of the roles of PGIPs and in screening potential combat proteins with novel recognition specificities against evolving pathogenic factors for countering pathogen invasion.

TOR (Target of Rapamycin) kinase is an evolutionarily conserved protein kinase, which integrates stress-related cues with growth and metabolic outputs. Long non-coding RNAs (lncRNAs) play a vital role in the regulation of eukaryotic genes. However, little is known about TOR's function in regulating the expression of lncRNAs in plants. In this study, four putative homologous genes encoding the TOR protein were identified by utilizing the recently completed cotton genome. Pharmacological experiments with TOR inhibitor AZD8055 and on silencing GhTOR genes resulted in obvious cotton growth retardation, indicating the conserved role of TOR in plant growth. The expression pattern analyses in different tissues reveal that TOR may play a role in root development, and the transcript levels of TOR genes were changed under different stress conditions. Importantly, we found TOR may be a key player in regulating the expression of long non-coding RNAs (lncRNAs). A total of 10,315 lncRNAs were discovered in cotton seedlings, 90.7% of which were long intergenic ncRNAs. Moreover, we identified the differentially expressed lncRNAs, of which 296 were significantly upregulated and 105 were downregulated in TOR inactivated plants. GO and KEGG analyses of differentially expressed lncRNA neighboring genes reveal that these differentially expressed lncRNA-targeted genes are involved in many life processes, including stress response, glutathione, and ribosomes in cotton. A series of differentially expressed lncRNAs potentially involved in plant stress response was identified under TOR inhibition. Collectively, these results suggest that cotton TOR proteins may directly modulate the expression of putative stress-related lncRNAs and eventually play a potential role in the cotton stress response.

Male-sterile lines are very important for selective breeding, and anther dehiscence defect is an effective way to generate male-sterile lines. Although several bHLH-family proteins in Arabidopsis have been characterized, little is known about the role of bHLH-family proteins in cotton. Here, we isolated a novel bHLH protein from cotton (Gossypium hirsutum), named GhBEE1-Like. Protein domain analysis showed that GhBEE1-Like contained a basic domain and an HLH domain. Subcellular localization analysis revealed that GhBEE1-Like was a nuclear-localized protein. Expression pattern analysis showed GhBEE1-Like was highly expressed in floral organs, and its expression was induced by the active brassinosteroid (BR) substance 24-epi-BL. GhBEE1-Like overexpression in Arabidopsis resulted in two types of transgenic lines, one with normal anther dehiscence and the other with defective anther dehiscence. Semi-qRT-PCR and qRT-PCR analyses revealed that GhBEE1-Like transcript levels acted as a check-point determining how anther dehiscence proceeds in these transgenic lines; regulated transcript levels result in normal anther dehiscence, whereas uncontrolled transcript levels lead to anther indehiscence. These results suggest that GhBEE1-Like plays an important role via its accumulation in regulating anther dehiscence. Therefore, controlling the level of GhBEE1-Like expression in cotton could be a convenient tool for generating male-sterile lines to use in selective breeding.

Metabolic pathways are interconnected and yet relatively independent. Genes involved in metabolic modules are required for the modules to run. Study of the relationships between genes and metabolic modules improves the understanding of metabolic pathways in plants. The WIN transcription factor activates the cuticle biosynthesis pathway and promotes cuticle biosynthesis. The relationship between the WIN transcription factor and other metabolic pathways is unknown. Our aim was to determine the relationships between the main genes involved in cuticle biosynthesis and those involved in other metabolic pathways. We did this by cloning a cotton WIN gene, GhWIN2, and studying its influence on other pathways.

More than 200 plants have been suffering from Verticillium wilt caused by Verticillium dahliae (V. dahliae) across the world. The target of rapamycin (TOR) is a lethal gene and controls cell growth and development in various eukaryotes, but little is known about TOR signaling in V. dahliae. Here, we found that V. dahliae strain is hypersensitive to rapamycin in the presence of rapamycin binding protein VdFKBP12 while the deletion mutant aaavdfkbp12 is insensitive to rapamycin. Heterologous expressing VdFKBP12 in Arabidopsis conferred rapamycin sensitivity, indicating that VdFKBP12 can bridge the interaction between rapamycin and TOR across species. The key across species of TOR complex 1 (TORC1) and TORC2 have been identified in V. dahliae, suggesting that TOR signaling pathway is evolutionarily conserved in eukaryotic species. Furthermore, the RNA-seq analysis showed that ribosomal biogenesis, RNA polymerase II transcription factors and many metabolic processes were significantly suppressed in rapamycin treated cells of V. dahliae. Importantly, transcript levels of genes associated with cell wall degrading enzymes (CWEDs) were dramatically down-regulated in TOR-inhibited cells. Further infection assay showed that the pathogenicity of V. dahliae and occurrence of Verticillium wilt can be blocked in the presence of rapamycin. These observations suggested that VdTOR is a key target of V. dahliae for controlling and preventing Verticillium wilt in plants.

Polyploidy is a common evolutionary occurrence in plants. Recently, published genomes of allotetraploid G. hirsutum and its donors G. arboreum and G. raimondii make cotton an accessible polyploid model. This study used chromatin immunoprecipitation with high-throughput sequencing (ChIP-Seq) to investigate the genome-wide distribution of H3K4me3 in G. arboreum and G. hirsutum, and explore the conservation and variation of genome structures between diploid and allotetraploid cotton. Our results showed that H3K4me3 modifications were associated with active transcription in both cottons. The H3K4me3 histone markers appeared mainly in genic regions and were enriched around the transcription start sites (TSSs) of genes. We integrated the ChIP-seq data of H3K4me3 with RNA-seq and ESTs data to refine the genic structure annotation. There were 6,773 and 12,773 new transcripts discovered in G. arboreum and G. hirsutum, respectively. Furthermore, co-expression networks were linked with histone modification and modularized in an attempt to explain differential H3K4me3 enrichment correlated with changes in gene transcription during cotton development and evolution. Taken together, we have combined epigenomic and transcriptomic datasets to systematically discover functional genes and compare them between G. arboreum and G. hirsutum, which may be beneficial for studying diploid and allotetraploid plants with large genomes and complicated evolution.

Melatonin functions as a plant hormone/regulator in the regulation of growth and development. However, the underlying mechanisms are still unclear. In this study, we found that a high dose of melatonin inhibited hypocotyl elongation in a dose-dependent manner in Arabidopsis. An expression profile analysis showed that hypocotyl growth inhibition by melatonin was involved in reprograming the expression of cell elongation genes and brassinosteroid (BRs) biosynthetic genes. Furthermore, similar to BR biosynthetic inhibitor brassinazole (BRZ), a high concentration of melatonin upregulated BR-biosynthetic genes and downregulated BR-induced genes involved in cell elongation, while melatonin was inefficient in brassinazole-resistant mutants like the bzr1-1D and bes1-D in hypocotyl inhibition. The comparative expression profile analysis showed an opposite expression mode in the co-regulated genes between melatonin and BZR1 or melatonin and brassinolide (BL). Additionally, exogenous BL rescued the repressive phenotype of BR biosynthesis-deficient mutant like det2-1 even in the presence of high-dose melatonin, but not BR receptor mutant bri1-5 or signal transduction mutant bin2-1. A biochemical analysis further confirmed that melatonin reduced endogenous BR levels in a dose-dependent manner in Arabidopsis. Taken together, these results indicate that melatonin inhibits BR biosynthesis but does not block BR signaling in the inhibition of hypocotyl elongation and extends insights on the role of melatonin in cross-talking with plant hormone signaling.

Extrafloral nectaries are a defence trait that plays important roles in plant-animal interactions. Gossypium species are characterized by cellular grooves in leaf midribs that secret large amounts of nectar. Here, with a panel of 215 G. arboreum accessions, we compared extrafloral nectaries to nectariless accessions to identify a region of Chr12 that showed strong differentiation and overlapped with signals from GWAS of nectaries. Fine mapping of an F2 population identified GaNEC1, encoding a PB1 domain-containing protein, as a positive regulator of nectary formation. An InDel, encoding a five amino acid deletion, together with a nonsynonymous substitution, was predicted to cause 3D structural changes in GaNEC1 protein that could confer the nectariless phenotype. mRNA-Seq analysis showed that JA-related genes are up-regulated and cell wall-related genes are down-regulated in the nectary. Silencing of GaNEC1 led to a smaller size of foliar nectary phenotype. Metabolomics analysis identified more than 400 metabolites in nectar, including expected saccharides and amino acids. The identification of GaNEC1 helps establish the network regulating nectary formation and nectar secretion, and has implications for understanding the production of secondary metabolites in nectar. Our results will deepen our understanding of plant-mutualism co-evolution and interactions, and will enable utilization of a plant defence trait in cotton breeding efforts.

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

Cystathionine β-synthase (CBS) domains containing proteins (CDCPs) form a large family and play roles in development via regulation of the thioredoxin system as well as abiotic and biotic stress responses of plant. However, the comprehensive study of CBS genes remained elusive in cotton. Here, we identified 237 CBS genes in 11 plant species and the phylogenetic analysis categorized CBS genes into four groups. Whole-genome or segmental with dispersed duplication events contributed to GhCBS gene family expansion. Moreover, orthologous/paralogous genes among three cotton species (G. hirsutum, G. arboreum, and G. raimondii) were detected from the syntenic map among eight plant species. Strong purifying selection for dicotyledonous and monocotyledonous CBS genes, and cis-elements related to plant growth and development, abiotic and hormonal response were observed. Transcriptomic data and qRT-PCR validation of 12 GhCBS genes indicated their critical role in ovule development as most of the genes showed high enrichment. Further, some of GhCBS (GhCBS5, GhCBS16, GhCBS17, GhCBS24, GhCBS25, GhCBS26, and GhCBS52) genes were regulated under various abiotic and hormonal treatments for different time points and involve in ovule and fiber development which provided key genes for future cotton breeding programs. In addition, transgenic tobacco plants overexpressing GhCBS4 transiently exhibited higher water and chlorophyll content indicating improved tolerance toward drought stress. Overall, this study provides the characterization of GhCBS genes for plant growth, abiotic and hormonal stresses, thereby, intimating their significance in cotton molecular breeding for resistant cultivars.

Immune responses that result in asthma exacerbation are associated with allergen or viral exposure. Identification of common immune factors will be beneficial for the development of uniformed targeted therapy. We employed a House Dust Mite (HDM) mouse model of asthma and challenged allergic HDM mice with allergens (HDM, cockroach extract (CRE)) or respiratory syncytial virus (RSV). Purified lung immune cells underwent high-dimensional single-cell RNA deep sequencing (scRNA-seq) to generate an RNA transcriptome. Gene silencing with siRNA was employed to confirm the efficacy of scRNA-seq analysis. scRNA-seq UMAP analysis portrayed an array of cell markers within individual immune clusters. SCENIC R analysis showed an increase in regulon number and activity in CD11b- DC cells. Analysis of conserved regulon factors further identified Creb5 as a shared regulon between the exacerbation groups. Creb5 siRNAs attenuated HDM, CRE or RSV-induced asthma exacerbation. scRNA-seq multidimensional analysis of immune clusters identified gene pathways that were conserved between the exacerbation groups. We propose that these analyses provide a strong framework that could be used to identify specific therapeutic targets in multifaceted pathologies.

Drought and high salinity are key limiting factors for cotton production. Therefore, research is increasingly focused on the underlying stress response mechanisms of cotton. We first identified and cloned a novel gene encoding the 525 amino acids in cotton, namely GhWRKY6. qRT-PCR analysis indicated that GhWRKY6 was induced by NaCl, PEG 6000 and ABA. Analyses of germination rate and root length indicated that overexpression of GhWRKY6 in Arabidopsis resulted in hypersensitivity to ABA, NaCl, and PEG 6000. In contrast, the loss-of-function mutant wrky6 was insensitive and had slightly longer roots than the wild-type did under these treatment conditions. Furthermore, GhWRKY6 overexpression in Arabidopsis modulated salt- and drought-sensitive phenotypes and stomatal aperture by regulating ABA signaling pathways, and reduced plant tolerance to abiotic stress through reactive oxygen species (ROS) enrichment, reduced proline content, and increased electrolytes and malondialdehyde (MDA). The expression levels of a series of ABA-, salt- and drought-related marker genes were altered in overexpression seedlings. Virus-induced gene silencing (VIGS) technology revealed that down-regulation of GhWRKY6 increased salt tolerance in cotton. These results demonstrate that GhWRKY6 is a negative regulator of plant responses to abiotic stress via the ABA signaling pathway.

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

Plants undergo a phase transition from vegetative to reproductive development that triggers floral induction. Genes containing an AAI (α-amylase inhibitor) domain form a large gene family, but there have been no comprehensive analyses of this gene family in any plant species. Here, we identified 336 AAI genes from nine plant species including122 AAI genes in cotton (Gossypium hirsutum). The AAI gene family has evolutionarily conserved amino acid residues throughout the plant kingdom. Phylogenetic analysis classified AAI genes into five major clades with significant polyploidization and showing effects of genome duplication. Our study identified 42 paralogous and 216 orthologous gene pairs resulting from segmental and whole-genome duplication, respectively, demonstrating significant contributions of gene duplication to expansion of the cotton AAI gene family. Further, GhAAI66 was preferentially expressed in flower tissue and as responses to phytohormone treatments. Ectopic expression of GhAAI66 in Arabidopsis and silencing in cotton revealed that GhAAI66 triggers a phase transition to induce early flowering. Further, GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis of RNA sequencing data and qRT-PCR (quantitative reverse transcription-PCR) analysis indicated that GhAAI66 integrates multiple flower signaling pathways including gibberellin, jasmonic acid, and floral integrators to trigger an early flowering cascade in Arabidopsis. Therefore, characterization of the AAI family provides invaluable insights for improving cotton breeding.

The CONSTANS (CO)-like gene family has been well studied for its role in the regulation of plant flowering time. However, their role remains poorly understood in cotton. To better understand the possible roles of CO-like in cotton, we performed a comprehensive genome-wide analysis of CO-like genes in cotton. Phylogenetic tree analysis showed that CO-like genes naturally clustered into three groups. Segmental duplication and whole genome duplication (WGD), which occurred before polyploidy, were important contributors to its expansion within the At ("t" indicates tetraploid) and Dt subgenomes, particularly in Group III. Long-terminal repeat retroelements were identified as the main transposable elements accompanying 18 genes. The genotype of GhCOL12_Dt displayed low diversity; it was a candidate involved in domestication. Selection pressure analyses indicated that relaxed purifying selection might have provided the main impetus during the evolution of CO-like genes in upland cotton. In addition, the high expression in the torus and calycle indicated that CO-like genes might affect flowering. The genes from Group II, and those from Group III involved in segmental duplication or WGD, might play important roles in response to drought and salt stress. Overall, this comprehensive genome-wide study of the CO-like gene family would facilitate further detailed studies in cotton.

Target of Rapamycin (TOR) is an eukaryotic protein kinase and evolutionally conserved from the last eukaryotic common ancestor (LECA) to humans. The growing evidences have shown that TOR signaling acts as a central controller of cell growth and development. The downstream effectors of TOR have been well-identified in yeast and animals by using the immunosuppression agent rapamycin. However, less is known about TOR in plants. This is largely due to the fact that plants are insensitive to rapamycin. In this study, AZD8055 (AZD), the novel ATP-competitive inhibitor of TOR, was employed to decipher the downstream effectors of TOR in Arabidopsis. One AZD insensitive mutant, T O R - i nhibitor i n sensitive- 1 (trin1), was screened from 10,000 EMS-induced mutation seeds. The cotyledons of trin1 can turn green when its seeds were germinated on ½ MS medium supplemented with 2 μM AZD, whereas the cotyledons greening of wild-type (WT) can be completely blocked at this concentration. Through genetic mapping, TRIN1 was mapped onto the long arm of chromosome 2, between markers SGCSNP26 and MI277. Positional cloning revealed that TRIN1 was an allele of ABI4, which encoded an ABA-regulated AP2 domain transcription factor. Plants containing P35S::TRIN1 or P35S::TRIN1-GUS were hypersensitive to AZD treatment and displayed the opposite phenotype observed in trin1. Importantly, GUS signaling was significantly enhanced in P35S::TRIN1-GUS transgenic plants in response to AZD treatment, indicating that suppression of TOR resulted in the accumulation of TRIN1. These observations revealed that TOR controlled seed-to-seedling transition by negatively regulating the stability of TRIN1 in Arabidopsis. For the first time, TRIN1, the downstream effector of TOR signaling, was identified through a chemical genetics approach.

Members of the NF-YB transcription factor gene family play important roles in diverse processes related to plant growth and development, such as seed development, drought tolerance, and flowering time. However, the function of NF-YB genes in cotton remains unclear. A total of 23, 24, and 50 NF-YB genes were identified in Gossypium arboreum (G. arboreum), Gossypium raimondii (G. raimondii), and G. hirsutum, respectively. A systematic phylogenetic analysis was carried out in G. arboretum, G. raimondii, G. hirsutum, Arabidopsis thaliana, cacao, rice and, sorghum, where the 150 NF-YB genes were divided into five groups (α-ε). Of these groups, α is the largest clade, and γ contains the LEC1 type NF-YB proteins. Syntenic analyses revealed that paralogues of NF-YB genes in G. hirsutum exhibited good collinearity. Owing to segmental duplication within the A sub-genome (At) and D sub-genome (Dt), there was an expanded set of NF-YB genes in G. hirsutum. Furthermore, we investigated the structures of exons, introns, and conserved motifs of NF-YB genes in upland cotton. Most of the NF-YB genes had only one exon, and the genes from the same clade exhibited a similar motif pattern. Expression data show that most NF-YB genes were expressed ubiquitously, and only a few genes were highly expressed in specific tissues, as confirmed by quantitative real-time PCR (qRT-PCR) analysis. The overexpression of GhDNF-YB22 gene, predominantly expressed in embryonic tissues, indicates that GhDNF-YB22 may affect embryogenesis in cotton. This study is the first comprehensive characterization of the GhNF-YB gene family in cotton, and showed that NF-YB genes could be divided into five clades. The duplication events that occurred over the course of evolution were the major impetus for NF-YB gene expansion in upland cotton. Collectively, this work provides insight into the evolution of NF-YB in cotton and further our knowledge of this commercially important species.

Asthma exacerbations are associated with heightened asthma symptoms, which can result in hospitalization in severe cases. However, the molecular immunologic processes that determine the course of an exacerbation remain poorly understood, impeding the progression of development of effective therapies.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been spreading globally for over 2 years, causing serious contagious disease and incalculable damage. The introduction of vaccines has slowed the spread of SARS-CoV-2 to some extent, but there remains a need for specific and effective treatment. The high chemical diversity and safety profiles of natural products make them a potential source of effective anti-SARS-CoV-2 drugs. Cotton plant is one of the most important economic and medical crops and is the source of a large number of antiviral phytochemicals. In this work, we used SARS-CoV-2 main protein (Mpro ) as the target to identify potential anti-SARS-CoV-2 natural products in cotton. An in vitro assay showed that of all cotton tissues examined, cotton flower extracts (CFs) exhibited optimal inhibitory effects against Mpro . We proceeded to use the CF metabolite database to screen natural Mpro inhibitors by combining virtual screening and biochemical assays. We identified that several CF natural products, including astragalin, myricitrin, and astilbin, significantly inhibited Mpro with half-maximal inhibitory concentrations (IC50s) of 0.13, 10.73, and 7.92 μm, respectively. These findings may serve as a basis for further studies into the suitability of cotton as a source of potential therapeutics for SARS-CoV-2.