Early archaeological evidence of wheat and cotton from medieval Ile-Ife, Nigeria.

This study reports the earliest directly dated occurrence of archaeological wheat and cotton in the humid forests of West Africa. These are the first archaeobotanical results from the medieval urban center of Ile-Ife, southwestern Nigeria, best known for its famous artworks. Both wheat and cotton likely spread through trans-Saharan trade networks that laid the foundation for later European trade systems. Forty-eight (48) grains of free-threshing wheat (Triticum aestivum/durum) represent the largest assemblage of wheat recovered in sub-Saharan West Africa, which is surprising given that wheat cannot be cultivated locally. Larger quantities of cotton (Gossypium sp.) recovered from late 12th- to early 13th-century CE contexts suggest earlier and more widespread use than wheat. Cotton may have been cultivated and manufactured into cloth locally. The quick adoption of these exotic crops illustrates the active negotiation of prestige through culinary and adornment practices, as well as a high degree of agricultural experimentation.

Promotion of apoplastic oxidative burst by artificially selected GhCBSX3A enhances Verticillium dahliae resistance in upland cotton.

Verticillium wilt (VW) is a devasting disease affecting various plants, including upland cotton, a crucial fiber crop. Despite its impact, the genetic basis underlying cotton's susceptibility or defense against VW remains unclear. Here, we conducted a genome-wide association study on VW phenotyping in upland cotton and identified a locus on A13 that is significantly associated with VW resistance. We then identified a cystathionine ß-synthase domain gene at A13 locus, GhCBSX3A, which was induced by Verticillium dahliae. Functional analysis, including expression silencing in cotton and overexpression in Arabidopsis thaliana, confirmed that GhCBSX3A is a causal gene at the A13 locus, enhancing SAR-RBOHs-mediated apoplastic oxidative burst. We found allelic variation on the TATA-box of GhCBSX3A promoter attenuated its expression in upland cotton, thereby weakening VW resistance. Interestingly, we discovered that altered artificial selection of GhCBSX3A_R (an elite allele for VW) under different VW pressures during domestication and other improved processes allows specific human needs to be met. Our findings underscore the importance of GhCBSX3A in response to VW, and we propose a model for defense-associated genes being selected depending on the pathogen's pressure. The identified locus and gene serve as promising targets for VW resistance enhancement in cotton through genetic engineering.

GhMYB30-GhMUR3 affects fiber elongation and secondary wall thickening in cotton.

Xyloglucan, an important hemicellulose, plays a crucial role in maintaining cell wall structure and cell elongation. However, the effects of xyloglucan on cotton fiber development are not well understood. GhMUR3 encodes a xyloglucan galactosyltransferase that is essential for xyloglucan synthesis and is highly expressed during fiber elongation. In this study, we report that GhMUR3 participates in cotton fiber development under the regulation of GhMYB30. Overexpression GhMUR3 affects the fiber elongation and cell wall thickening. Transcriptome showed that the expression of genes involved in secondary cell wall synthesis was prematurely activated in OE-MUR3 lines. In addition, GhMYB30 was identified as a key regulator of GhMUR3 by Y1H, Dual-Luc, and electrophoretic mobility shift assay (EMSA) assays. GhMYB30 directly bound the GhMUR3 promoter and activated GhMUR3 expression. Furthermore, DAP-seq of GhMYB30 was performed to identify its target genes in the whole genome. The results showed that many target genes were associated with fiber development, including cell wall synthesis-related genes, BR-related genes, reactive oxygen species pathway genes, and VLCFA synthesis genes. It was demonstrated that GhMYB30 may regulate fiber development through multiple pathways. Additionally, GhMYB46 was confirmed to be a target gene of GhMYB30 by EMSA, and GhMYB46 was significantly increased in GhMYB30-silenced lines, indicating that GhMYB30 inhibited GhMYB46 expression. Overall, these results revealed that GhMUR3 under the regulation of GhMYB30 and plays an essential role in cotton fiber elongation and secondary wall thickening. Additionally, GhMYB30 plays an important role in the regulation of fiber development and regulates fiber secondary wall synthesis by inhibiting the expression of GhMYB46.

Yield-related quantitative trait loci identification and lint percentage hereditary dissection under salt stress in upland cotton.

Salinity is frequently mentioned as a major constraint in worldwide agricultural production. Lint percentage (LP) is a crucial yield-component in cotton lint production. While the genetic factors affect cotton yield in saline soils are still unclear. Here, we employed a recombinant inbred line population in upland cotton (Gossypium hirsutum L.) and investigated the effects of salt stress on five yield and yield component traits, including seed cotton yield per plant, lint yield per plant, boll number per plant, boll weight, and LP. Between three datasets of salt stress (E1), normal growth (E2), and the difference values dataset of salt stress and normal conditions (D-value), 87, 82, and 55 quantitative trait loci (QTL) were detectable, respectively. In total, five QTL (qLY-Chr6-2, qBNP-Chr4-1, qBNP-Chr12-1, qBNP-Chr15-5, qLP-Chr19-2) detected in both in E1 and D-value were salt related QTL, and three stable QTL (qLP-Chr5-3, qLP-Chr13-1, qBW-Chr5-5) were detected both in E1 and E2 across 3 years. Silencing of nine genes within a stable QTL (qLP-Chr5-3) highly expressed in fiber developmental stages increased LP and decreased fiber length (FL), indicating that multiple minor-effect genes clustered on Chromosome 5 regulate LP and FL. Additionally, the difference in LP caused by Gh_A05G3226 is mainly in transcription level rather than in the sequence difference. Moreover, silencing of salt related gene (GhDAAT) within qBNP-Chr4-1 decreased salt tolerance in cotton. Our findings shed light on the regulatory mechanisms underlining cotton salt tolerance and fiber initiation.

SR45a plays a key role in enhancing cotton resistance to Verticillium dahliae by alternative splicing of immunity genes.

Alternative splicing (AS) of pre-mRNAs increases the diversity of transcriptome and proteome and plays fundamental roles in plant development and stress responses. However, the prevalent changes in AS events and the regulating mechanisms of plants in response to pathogens remain largely unknown. Here, we show that AS changes are an important mechanism conferring cotton immunity to Verticillium dahliae (Vd). GauSR45a, encoding a serine/arginine-rich RNA binding protein, was upregulated expression and underwent AS in response to Vd infection in Gossypium australe, a wild diploid cotton species highly resistant to Vd. Silencing GauSR45a substantially reduced the splicing ratio of Vd-induced immune-associated genes, including GauBAK1 (BRI1-associated kinase 1) and GauCERK1 (chitin elicitor receptor kinase 1). GauSR45a binds to the GAAGA motif that is commonly found in the pre-mRNA of genes essential for PTI, ETI, and defense. The binding between GauSR45a and the GAAGA motif in the pre-mRNA of BAK1 was enhanced by two splicing factors of GauU2AF35B and GauU1-70 K, thereby facilitating exon splicing; silencing either AtU2AF35B or AtU1-70 K decreased the resistance to Vd in transgenic GauSR45a Arabidopsis. Overexpressing the short splicing variant of BAK1GauBAK1.1 resulted in enhanced Verticillium wilt resistance rather than the long one GauBAK1.2. Vd-induced far more AS events were in G. barbadense (resistant tetraploid cotton) than those in G. hirsutum (susceptible tetraploid cotton) during Vd infection, indicating resistance divergence in immune responses at a genome-wide scale. We provided evidence showing a fundamental mechanism by which GauSR45a enhances cotton resistance to Vd through global regulation of AS of immunity genes.

GhEXL3 participates in brassinosteroids regulation of fiber elongation in Gossypium hirsutum.

Cotton fiber (Gossypium hirsutum) serves as an ideal model for investigating the molecular mechanisms of plant cell elongation at the single-cell level. Brassinosteroids (BRs) play a crucial role in regulating plant growth and development. However, the mechanism by which BR influences cotton fiber elongation remains incompletely understood. In this study, we identified EXORDIUM-like (GhEXL3) through transcriptome analysis of fibers from BR-deficient cotton mutant pagoda 1 (pag1) and BRI1-EMS-SUPPRESSOR 1 (GhBES1.4, encoding a central transcription factor of BR signaling) overexpression cotton lines. Knockout of GhEXL3 using CRISPR/Cas9 was found to impede cotton fiber elongation, while its overexpression promoted fiber elongation, suggesting a positive regulatory function for GhEXL3 in fiber elongation. Furthermore, in vitro ovule culture experiments revealed that the overexpression of GhEXL3 partially counteracted the inhibitory effects of brassinazole (BRZ) on cotton fiber elongation, providing additional evidence of GhEXL3 involvement in BR signaling pathways. Moreover, our findings demonstrate that GhBES1.4 directly binds to the E-box (CACGTG) motif in the GhEXL3 promoter region and enhances its transcription. RNA-seq analysis revealed that overexpression of GhEXL3 upregulated the expression of EXPs, XTHs, and other genes associated with fiber cell elongation. Overall, our study contributes to understanding the mechanism by which BR regulates the elongation of cotton fibers through the direct modulation of GhEXL3 expression by GhBES1.4.

Suppression of plant immunity by Verticillium dahliae effector Vd6317 through AtNAC53 association.

Verticillium dahliae, a soil-borne fungal pathogen, compromises host innate immunity by secreting a plethora of effectors, thereby facilitating host colonization and causing substantial yield and quality losses. The mechanisms underlying the modulation of cotton immunity by V. dahliae effectors are predominantly unexplored. In this study, we identified that the V. dahliae effector Vd6317 inhibits plant cell death triggered by Vd424Y and enhances PVX viral infection in Nicotiana benthamiana. Attenuation of Vd6317 significantly decreased the virulence of V. dahliae, whereas ectopic expression of Vd6317 in Arabidopsis and cotton enhanced susceptibility to V. dahliae infection, underscoring Vd6317's critical role in pathogenicity. We observed that Vd6317 targeted the Arabidopsis immune regulator AtNAC53, thereby impeding its transcriptional activity on the defense-associated gene AtUGT74E2. Arabidopsis nac53 and ugt74e2 mutants exhibited heightened sensitivity to V. dahliae compared to wild-type plants. A mutation at the conserved residue 193L of Vd6317 abrogated its interaction with AtNAC53 and reduced the virulence of V. dahliae, which was partially attributable to a reduction in Vd6317 protein stability. Our findings unveil a hitherto unrecognized regulatory mechanism by which the V. dahliae effector Vd6317 directly inhibits the plant transcription factor AtNAC53 activity to suppress the expression of AtUGT74E2 and plant defense.

Synergistic interplay of redox homeostasis and polysaccharide synthesis promotes cotton fiber elongation.

Cell polarity is the foundation of cell development and tissue morphogenesis. The investigation of polarized growth provides opportunities to gain profound insights into morphogenesis and tissue functionality in organisms. Currently, there are still many mysteries surrounding the mechanisms that regulate polarized cell growth. Cotton fiber cells serve as an excellent model for studying polarized growth, and provide important clues for unraveling the molecular mechanisms, signaling pathways, and regulatory networks of polarized growth. In this study, we characterized two functional genes, GhMDHAR1AT/DT and GhDHAR2AT/DT with predominant expression during fiber elongation. Loss of function of both genes contributed to a significant increase in fiber length. Transcriptomic data revealed up-regulated expression of antioxidant genes in CRISPR mutant lines, along with delayed expression of secondary wall-related genes and temporally prolonged expression of primary wall-related genes. Experimental evidence demonstrated that the increase in GSH content and glutathione peroxidase (GPX) enzyme activity led to enhanced total antioxidant capacity (T-AOC), resulting in reduced H2O2 levels, which contributed to the extension of fiber elongation stage in CRISPR mutant lines. Moreover, the increased polysaccharide synthesis in CRISPR mutant lines was found to provide an abundant supply of raw materials for fiber cell wall elongation, suggesting that synergistic interplay between redox homeostasis and polysaccharide synthesis in fiber cells may facilitate cell wall remodeling and fiber elongation. This study provides valuable insights for deciphering the mechanisms of cell polarized growth and improving cotton fiber quality.

Phosphatidic acid interacts with an HD-ZIP transcription factor GhHOX4 to influence its function in fiber elongation of cotton (Gossypium hirsutum).

Upland cotton, the mainly cultivated cotton species in the world, provides over 90% of natural raw materials (fibers) for the textile industry. The development of cotton fibers that are unicellular and highly elongated trichomes on seeds is a delicate and complex process. However, the regulatory mechanism of fiber development is still largely unclear in detail. In this study, we report that a homeodomain-leucine zipper (HD-ZIP) IV transcription factor, GhHOX4, plays an important role in fiber elongation. Overexpression of GhHOX4 in cotton resulted in longer fibers, while GhHOX4-silenced transgenic cotton displayed a "shorter fiber" phenotype compared with wild type. GhHOX4 directly activates two target genes, GhEXLB1D and GhXTH2D, for promoting fiber elongation. On the other hand, phosphatidic acid (PA), which is associated with cell signaling and metabolism, interacts with GhHOX4 to hinder fiber elongation. The basic amino acids KR-R-R in START domain of GhHOX4 protein are essential for its binding to PA that could alter the nuclear localization of GhHOX4 protein, thereby suppressing the transcriptional regulation of GhHOX4 to downstream genes in the transition from fiber elongation to secondary cell wall (SCW) thickening during fiber development. Thus, our data revealed that GhHOX4 positively regulates fiber elongation, while PA may function in the phase transition from fiber elongation to SCW formation by negatively modulating GhHOX4 in cotton.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds.

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton.

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Molecular traits of MAPK kinases and the regulatory mechanism of GhMAPKK5 alleviating drought/salt stress in cotton.

Mitogen-activated protein kinase kinases (MAPKKs) play a critical role in the mitogen-activated protein kinase (MAPK) signaling pathway, transducing external stimuli into intracellular responses and enabling plant adaptation to environmental challenges. Most research has focused on the model plant Arabidopsis (Arabidopsis thaliana). The systematic analysis and characterization of MAPKK genes across different plant species, particularly in cotton (Gossypium hirsutum), are somewhat limited. Here, we identified MAPKK family members from 66 different species, which clustered into 5 different sub-groups, and MAPKKs from four cotton species clustered together. Through further bioinformatic and expression analysis, GhMAPKK5 was identified as the most responsive MAPKK member to salt and drought stress among the 23 MAPKKs identified in Gossypium hirsutum. Silencing GhMAPKK5 in cotton through virus-induced gene silencing (VIGS) led to quicker wilting under salt and drought conditions, while overexpressing GhMAPKK5 in Arabidopsis enhanced root growth and seed germination under these stresses, demonstrating GhMAPKK5's positive role in stress tolerance. Transcriptomics and Yeast-Two-Hybrid assays revealed a MAPK cascade signal module comprising GhMEKK (Mitogen-activated protein kinase kinase kinases)3/8/31-GhMAPKK5-GhMAPK11/23. This signaling cascade may play a role in managing drought and salt stress by regulating transcription factor genes, such as WRKYs, which are involved in the biosynthesis and transport pathways of ABA, proline, and RALF. This study is highly important for further understanding the regulatory mechanism of MAPKK in cotton, contributing to its stress tolerance and offering potential in targets for genetic enhancement.

Using a cellulose-complementary oligosaccharide as a tool to probe exposed cellulosic surfaces in cotton fibres and growing plant cell walls.

Cellulosic microfibrils in plant cell walls are largely ensheathed and probably tethered by hydrogen-bonded hemicelluloses. Ensheathing may vary developmentally as hemicelluloses are peeled to enable cell expansion. We characterised a simple method to quantify ensheathed versus naked cellulosic surfaces based on the ability to adsorb a radiolabelled 'cellulose-complementary oligosaccharide', [3H]cellopentaitol. Filter-paper (cellulose) adsorbed 40% and >80% of aqueous 5nM [3H]cellopentaitol within ~1h and ~20h respectively. When [3H]cellopentaitol was rapidly dried onto filter-paper, ~50% of it was desorbable by water, whereas after ~1d annealing in aqueous medium the adsorption became too strong to be reversible in water. 'Strongly' adsorbed [3H]cellopentaitol was, however, ~98% desorbed by 6M NaOH, ~50% by 0.2M cellobiose, and ~30% by 8M urea, indicating a role for hydrogen-bonding reinforced by complementarity of shape. Gradual adsorption was promoted by kosmotropes (1.4M Na2SO4 or 30% methanol), and inhibited by chaotropes (8M urea), supporting a role for hydrogen-bonding. [3H]Cellopentaitol adsorption was strongly competed by non-radioactive cello-oligosaccharides (Cell2-6), the IC50 (half-inhibitory concentration) being highly size-dependent: Cell2, ~70 mM; Cell3, ~7 mM; and Cell4-6, ~0.05 mM. Malto-oligosaccharides (400mM) had no effect, confirming the role of complementarity. The quantity of adsorbed [3H]cellopentaitol was proportional to mass of cellulose. Of seven cottons tested, wild-type Gossypium arboreum fibres were least capable of adsorbing [3H]cellopentaitol, indicating ensheathment of their microfibrillar surfaces, confirmed by their resistance to cellulase digestion, and potentially attributable to a high glucuronoarabinoxylan content. In conclusion, [3H]cellopentaitol adsorption is a simple, sensitive and quantitative way of titrating 'naked' cellulose surfaces.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton.

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a "genome shock", leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton (Gossypium hirsutum and Gossypium barbadense, AADD) and its extant AA (Gossypium arboreum) and DD (Gossypium raimondii) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii. In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense, but absent from a resynthesized hybrid of G. arboreum and G. raimondii, suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis-regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Pangenome analysis reveals transposon-driven genome evolution in cotton.

BACKGROUND: Transposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution. RESULTS: Here, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes. CONCLUSIONS: Our pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton.

BACKGROUND: Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors. RESULTS: In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms. CONCLUSIONS: These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility.

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Cotton peroxisome-localized lysophospholipase counteracts the toxic effects of Verticillium dahliae NLP1 and confers wilt resistance.

Plasma membrane represents a critical battleground between plants and attacking microbes. Necrosis-and-ethylene-inducing peptide 1 (Nep1)-like proteins (NLPs), cytolytic toxins produced by some bacterial, fungal and oomycete species, are able to target on lipid membranes by binding eudicot plant-specific sphingolipids (glycosylinositol phosphorylceramide) and form transient small pores, causing membrane leakage and subsequent cell death. NLP-producing phytopathogens are a big threat to agriculture worldwide. However, whether there are R proteins/enzymes that counteract the toxicity of NLPs in plants remains largely unknown. Here we show that cotton produces a peroxisome-localized enzyme lysophospholipase, GhLPL2. Upon Verticillium dahliae attack, GhLPL2 accumulates on the membrane and binds to V. dahliae secreted NLP, VdNLP1, to block its contribution to virulence. A higher level of lysophospholipase in cells is required to neutralize VdNLP1 toxicity and induce immunity-related genes expression, meanwhile maintaining normal growth of cotton plants, revealing the role of GhLPL2 protein in balancing resistance to V. dahliae and growth. Intriguingly, GhLPL2 silencing cotton plants also display high resistance to V. dahliae, but show severe dwarfing phenotype and developmental defects, suggesting GhLPL2 is an essential gene in cotton. GhLPL2 silencing results in lysophosphatidylinositol over-accumulation and decreased glycometabolism, leading to a lack of carbon sources required for plants and pathogens to survive. Furthermore, lysophospholipases from several other crops also interact with VdNLP1, implying that blocking NLP virulence by lysophospholipase may be a common strategy in plants. Our work demonstrates that overexpressing lysophospholipase encoding genes have great potential for breeding crops with high resistance against NLP-producing microbial pathogens.

High-resolution sequencing of nine elite upland cotton cultivars uncovers genic variations and breeding improvement targets.

Structural variations (SVs) are critical factors affecting genome evolution and important traits. However, identification results and functional analyses of SVs in upland cotton are rare. Here, based on the genetic relationships, breeding history and cumulative planting area of upland cotton in China, nine predominant cultivars from the past 60 years (1950s-2010s) were selected for long read sequencing to uncover genic variations and breeding improvement targets for this crop. Based on the ZM24 reference genome, 0.88-1.47 × 104 SVs per cultivar were identified, and an SV set was constructed. SVs affected the expression of a large number of genes during fiber elongation, and a transposable element insertion resulted in the glandless phenotype in upland cotton. Six widespread inversions were identified based on nine draft genomes and high-throughput chromosome conformation capture data. Multiple haplotype blocks that were always associated with aggregated SVs were demonstrated to play a pivotal role in the agronomic traits of upland cotton and drove its adaptation to the northern planting region. Exotic introgression was the source of these haplotype blocks and increased the genetic diversity of upland cotton. Our results enrich the genome resources of upland cotton, and the identified SVs will promote genetic and breeding research in cotton.