Extensive intraspecific gene order and gene structural variations in upland cotton cultivars

Pangenome analysis reveals yield- and fiber-related diversity and interspecific gene flow in Gossypium barbadense L

Gossypium barbadense is renowned for its superior fiber quality, particularly its extra-long fibers, although its fiber yield is lower compared to G. hirsutum. Here, to further reveal fiber-related genomic variants of G. barbadense, we de novo assemble 12 genomes of G. barbadense that span the wild-to-domesticated continuum, and construct a graph-based pangenome by integrating these assemblies and 17 publicly available tetraploid cotton genome assemblies. We uncover the divergent evolutionary trajectories and subsequent exchanges between G. barbadense and G. hirsutum through investigation of structural variants (SVs). We perform the SV-based GWAS analysis in G. barbadense and identify four, three, and seven candidate SVs for fiber length, fiber strength, and lint percentage, respectively. Furthermore, we detect the underlying candidate genes and uncover the origin and distribution of favorable alleles, and reveal the tradeoff between lint percentage and fiber quality. These pangenome and trait-associated SVs provide insights into and resources for improving cotton fiber.

Two duplicated GhMML3 genes coordinately control development of lint and fuzz fibers in cotton

Cotton produces two types of fibers: fuzz and lint. Cotton yield is determined by the number of epidermal cells that develop into lint fibers. Despite numerous studies, the genetic and molecular mechanisms that control lint and fuzz fiber development remain unclear. Here, using the recessive naked-seed or fuzzless-linted mutant (n2NSM) in combination with gene editing and complementation, we found that the recessive fuzzless gene n2 encodes the MYBMIXTA-like (MML) transcription factor GhMML3_D12. Overexpression of GhMML3_D12 in n2NSM restored fuzz fiber development, whereas CRISPR-Cas9 knockout of GhMML3_D12 in wild-type cotton (J668) resulted in a fuzzless-linted phenotype. Interestingly, simultaneous edits to GhMML3_D12 and its duplicate GhMML3_A12 resulted in plants with a fiberless (fuzzless-lintless) phenotype. Detailed investigation of the seed fiber phenotypes of segregating progeny derived from a cross between J668 and a fiberless gene-edited mutant of GhMML3 (#mml3s) not only identified progeny that mimicked natural fuzzless and fiberless mutants but also revealed that the duplicated GhMML3_A12 and GhMML3_D12 regulate the development of fuzz and lint fibers in a dose-dependent manner. Comparative transcriptome analysis and single-cell RNA sequencing identified GhMML3 as the central hub of the gene network that regulates fiber initiation and early-stage elongation. The gene regulatory network revealed potential candidate genes and key regulators that may contribute to fiber initiation and development, and a model for the control of lint and fuzz fiber development by GhMML3 was proposed. We also found that the GhMML3_D12 protein can bind directly to the promoters of GhHD-1 and GhMYB25, two key genes involved in fiber initiation, thereby activating their expression. This study provides new insights into the fundamental mechanisms that underlie cotton fiber development.

Candida albicans isolates contain frequent heterozygous structural variants and transposable elements within genes and centromeres

The human fungal pathogen Candida albicans poses a significant burden on global health, causing high rates of mortality and antifungal drug resistance. C. albicans is a heterozygous diploid organism that reproduces asexually. Structural variants (SVs) are an important source of genomic rearrangement, particularly in species that lack sexual recombination. To comprehensively investigate SVs across clinical isolates of C. albicans, we conducted long-read sequencing and genome-wide SV analysis in three distantly related clinical isolates. Our work includes a new, comprehensive analysis of transposable element (TE) composition, location, and diversity. SVs and TEs are frequently close to coding sequences and many SVs are heterozygous, suggesting that SVs might impact gene and allele-specific expression. Most SVs are uniquely present in only one clinical isolate, indicating that SVs represent a significant source of intraspecies genetic variation. We identify multiple, distinct SVs at the centromeres of Chromosome 4 and Chromosome 5, including inversions and transposon polymorphisms. These two chromosomes are often aneuploid in drug-resistant clinical isolates and can form isochromosome structures with breakpoints near the centromere. Further screening of 100 clinical isolates confirms the widespread presence of centromeric SVs in C. albicans, often appearing in a heterozygous state, indicating that SVs are contributing to centromere evolution in C. albicans Together, these findings highlight that SVs and TEs are common across diverse clinical isolates of C. albicans and that the centromeres of this organism are important sites of genome rearrangement.

Chromosome-level genome assembly of the crofton weed (Ageratina adenophora)

Crofton weed (Ageratina adenophora), a significant invasive species, extensively disrupts ecosystem stability, leading to considerable economic losses. However, genetic insights into its invasive mechanisms have been limited by a lack of genomic data. In this study, we present the successful de novo assembly of the triploid genome of A. adenophora, leveraging long-read PacBio Sequel, optical mapping, and Hi-C sequencing. Our assembly resolved into a haplotype-resolved genome comprising 51 chromosomes, with a total size of ~3.82 Gb and a scaffold N50 of 70.8 Mb. BUSCO analysis confirmed the completeness of 97.71% of genes. Genome annotation revealed 3.16 Gb (76.44%) of repetitive sequences and predicted 123,134 protein-coding genes, with 99.03% functionally annotated. The high-quality reference genome will provide valuable genomic resources for future studies on the evolutionary dynamics and invasive adaptations of A. adenophora.

A telomere-to-telomere genome assembly of cotton provides insights into centromere evolution and short-season adaptation

Cotton (Gossypium hirsutum L.) is a key allopolyploid crop with global economic importance. Here we present a telomere-to-telomere assembly of the elite variety Zhongmian 113. Leveraging technologies including PacBio HiFi, Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C, our assembly surpasses previous genomes in contiguity and completeness, resolving 26 centromeric and 52 telomeric regions, 5S rDNA clusters and nucleolar organizer regions. A phylogenetically recent centromere repositioning on chromosome D08 was discovered specific to G. hirsutum, involving deactivation of an ancestral centromere and the formation of a unique, satellite repeat-based centromere. Genomic analyses evaluated favorable allele aggregation for key agronomic traits and uncovered an early-maturing haplotype derived from an 11 Mb pericentric inversion that evolved early during G. hirsutum domestication. Our study sheds light on the genomic origins of short-season adaptation, potentially involving introgression of an inversion from primitively domesticated forms, followed by subsequent haplotype differentiation in modern breeding programs.

Post-polyploidization centromere evolution in cotton

Upland cotton (Gossypium hirsutum) accounts for more than 90% of the world's cotton production and, as an allotetraploid, is a model plant for polyploid crop domestication. In the present study, we reported a complete telomere-to-telomere (T2T) genome assembly of Upland cotton accession Texas Marker-1 (T2T-TM-1), which has a total size of 2,299.6 Mb, and annotated 79,642 genes. Based on T2T-TM-1, interspecific centromere divergence was detected between the A- and D-subgenomes and their corresponding diploid progenitors. Centromere-associated repetitive sequences (CRCs) were found to be enriched for Gypsy-like retroelements. Centromere size expansion, repositioning and structure variations occurred post-polyploidization. It is interesting that CRC homologs were transferred from the diploid D-genome progenitor to the D-subgenome, invaded the A-subgenome and then underwent post-tetraploidization proliferation. This suggests an evolutionary advantage for the CRCs of the D-genome progenitor, presents a D-genome-adopted inheritance of centromere repeats after polyploidization and shapes the dynamic centromeric landscape during polyploidization in polyploid species.

Tandem LTR-retrotransposon structures are common and highly polymorphic in plant genomes

BACKGROUND

LTR-retrotransposons (LTR-RT) are a major component of plant genomes and important drivers of genome evolution. Most LTR-RT copies in plant genomes are defective elements found as truncated copies, nested insertions or as part of more complex structures. The recent availability of highly contiguous plant genome assemblies based on long-read sequences now allows to perform detailed characterization of these complex structures and to evaluate their importance for plant genome evolution.

RESULTS

The detailed analysis of two rice loci containing complex LTR-RT structures showed that they consist of tandem arrays of LTR copies sharing internal LTRs. Our analyses suggests that these LTR-RT tandems are the result of a single insertion and not of the recombination of two independent LTR-RT elements. Our results also suggest that gypsy elements may be more prone to form these structures. We show that these structures are highly polymorphic in rice and therefore have the potential to generate genetic variability. We have developed a computational pipeline (IDENTAM) that scans genome sequences and identifies tandem LTR-RT candidates. Using this tool, we have detected 266 tandems in a pangenome built from the genomes of 76 accessions of cultivated and wild rice, showing that tandem LTR-RT structures are frequent and highly polymorphic in rice. Running IDENTAM in the Arabidopsis, almond and cotton genomes showed that LTR-RT tandems are frequent in plant genomes of different size, complexity and ploidy level. The complexity of differentiating intra-element variations at the nucleotide level among haplotypes is very high, and we found that graph-based pangenomic methodologies are appropriate to resolve these structures.

CONCLUSIONS

Our results show that LTR-RT elements can form tandem arrays. These structures are relatively abundant and highly polymorphic in rice and are widespread in the plant kingdom. Future studies will contribute to understanding how these structures originate and whether the variability that they generate has a functional impact.

Cotton metabolism regulatory network: Unraveling key genes and pathways in fiber development and growth regulation

Cotton (Gossypium hirsutum L.) is one of the world's most important commercial crops. However, the dynamics of metabolite abundance and potential regulatory networks throughout its life cycle remain poorly understood. In this study, we developed a cotton metabolism regulatory network (CMRN) that spans various developmental stages and encompasses 2138 metabolites and 90 309 expressed genesin upland cotton. By integrating high-resolution spatiotemporal metabolome and transcriptome data, we identified 1958 differentially accumulated metabolites and 13 597 co-expressed differentially expressed genes between the dwarf mutant pagoda1 and its wild-type counterpart Zhongmiansuo 24. These metabolites and genes were categorized into seven clusters based on tissue-specific accumulation patterns and gene expression profiles across different developmental stages. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed significant differential enrichment in the fatty acid elongation pathway, particularly in fibers. The differential involvement of genes and metabolites in very-long-chain fatty acid (VLCFA) synthesis led to the identification of GhKCS1b_Dt as a key gene. Overexpression of GhKCS1b_Dt significantly promoted fiber elongation, while its silencing markedly inhibited cotton fiber growth, affirming its positive regulatory role in fiber elongation. This dataset provides a valuable resource for further research into metabolic pathways and gene regulatory networks, offering novel insights for advancing cotton breeding strategies.

Identification of the cysteine-rich transmembrane module CYSTM family in upland cotton and functional analysis of GhCYSTM5_A in cold and drought stresses

Abiotic stress poses adverse impacts on cotton production, raising demands for a better understanding of stress-response mechanisms and developing strategies to improve plant performance to cope with stress. CYSTM (Cysteine-rich transmembrane module) is a widely distributed and conserved family in eukaryotes that performs potential functions in stress tolerance. However, CYSTM genes and their role in stress response is uncharacterized in cotton. Herein, we identified a total of 23 CYSTM genes from upland cotton. They underwent mainly segmental duplications and experienced purifying selection during evolution. Expression profiles revealed GhCYSTMs were closely related to abiotic stress response. Furthermore, GhCYSTM5_A overexpression enhanced the cold and drought tolerance of cotton, while RNAi-mediated knockdown of GhCYSTM5_A decreased stress tolerance. Transcriptome analysis revealed GhCYSTM5_A may contribute to cold and drought tolerance by regulating the expression of oxidative stress-related genes through MAPK signaling. GhCYSTM5_A, localized in the nucleus and cytoplasm interacted with a secreted cysteine-rich peptide GhGASA14. Moreover, GhGASA14 silencing rendered cotton plants vulnerable to cold and drought. These results suggested the potential functions of GhCYSTM genes in abiotic stress and a positive role of GhCYSTM5_A in cold and drought tolerance. This study sheds light on comprehensive characteristics of GhCYSTM, and provides candidate genes for genetic breeding.

Introgression and Mapping of Cotton Leaf Curl Disease Resistance from Wild Gossypium armourianum Kearney into Upland Cotton (G. hirsutum)

Cotton leaf curl disease (CLCuD), caused by the whitefly transmitted geminivirus complex (cotton leaf curl virus [CLCuV] and their satellite molecules), is a serious threat to successful upland cotton production in northwest India and Pakistan. The disease causes significant losses in fiber yield and the quality of cotton. Owing to the regular emergence of resistance-breaking strains of CLCuV, all the previously available CLCuD-resistant germplasms of upland cotton have become compromised, and none of the extant upland cotton cultivars is resistant to this disease. Therefore, alternate sources of CLCuD resistance need to be explored, as genetic resistance is the only pragmatic and tenable management strategy to combat this malady. Here, we report for the first time the introgression and mapping of CLCuD resistance from a related nonprogenitor wild diploid D-genome cotton species, Gossypium armourianum, into upland cotton. A backcross population (G. hirsutum/G. armourianum/G. hirsutum) was developed for this purpose. A single major QTL was found to be associated with resistance to CLCuD and was located on chromosome D01 through the genotyping-by-sequencing technique.

Comparative analysis using a chromosome-scale genome assembly for Funaria hygrometrica suggests greater collinearity in mosses than in seed plants

Mosses, the largest lineage of seed-free plants, have smaller and less variable genome sizes than flowering plants. Nevertheless, whether this difference results from divergent genome dynamics is poorly known. Here, we use newly generated chromosome-scale genome assemblies for Funaria hygrometrica and comparative analysis with other moss and seed plant genomes to investigate moss genome dynamics. Although some aspects of moss genome dynamics are seed plant-like, such as the mechanism of genome size change and de novo gain/loss of genes, moss genomes retain higher synteny, and collinearity over evolutionary time than seed plant genomes. Furthermore, transposable elements and genes are more evenly distributed along chromosomes in mosses than in seed plants, a feature shared with other sequenced seed-free plant genomes. Overall, our findings support the hypothesis that large-scale genome structure and dynamics of mosses and seed plants differ. In particular, our data suggest a lower rate of gene order reshuffling along chromosomes in mosses compared to seed plants. We speculate that such lower rate of structural genomic variation and unique chromosome structure in mosses may contribute to their relatively smaller and less variable genome sizes.

Unraveling key genes and pathways involved in Verticillium wilt resistance by integrative GWAS and transcriptomic approaches in Upland cotton

Verticillium dahliae Kleb, the cause of Verticillium wilt, is a particularly destructive soil-borne vascular disease that affects cotton, resulting in serious decline in fiber quality and causing significant losses in cotton production worldwide. However, the progress in identification of wilt-resistance loci or genes in cotton has been limited, most probably due to the highly complex genetic nature of the trait. Nevertheless, the molecular mechanism behind the Verticillium wilt resistance remains poorly understood. In the present study, we investigated the phenotypic variations in Verticillium tolerance and conducted a genome wide association study (GWAS) among a natural population containing 383 accessions of upland cotton germplasm and performed transcriptomic analysis of cotton genotypes with differential responses to Verticillium wilt. GWAS detected 70 significant SNPs and 116 genes associated with resistance loci in two peak signals on D02 and D11 in E1. The transcriptome analysis identified a total of 2689 and 13289 differentially expressed genes (DEGs) among the Verticillium wilt-tolerant (J46) and wilt-susceptible (J11) genotypes, respectively. The DEGs were predominantly enriched in metabolism, plant hormone signal transduction, phenylpropanoid pathway, MAPK cascade pathway and plant-pathogen interaction pathway in GO and KEGG analyses. The identified DEGs were found to comprise several transcription factor (TF) gene families, primarily including AP2/ERF, ZF, WRKY, NAC and MYB, in addition to pentatricopeptide repeat (PPR) proteins and Resistance (R) genes. Finally, by integrating the two results, 34 candidate genes were found to overlap between GWAS and RNA-seq analyses, associated with Verticillium-wilt resistance, including WRKY, MYB, CYP and RGA. This work contributes to our knowledge of the molecular processes underlying cotton responses to Verticillium wilt, offering crucial insights for additional research into the genes and pathways implicated in these responses and paving the way for developing Verticillium wilt-resistant cotton varieties through accelerated breeding by providing a plethora of candidate genes.

Cotton2035: From genomics research to optimized breeding

Cotton is the world's most important natural fiber crop and serves as an ideal model for studying plant genome evolution, cell differentiation, elongation, and cell wall biosynthesis. The first draft genome assembly for Gossypium raimondii, completed in 2012, marked the beginning of global efforts in studying cotton genomics. Over the past decade, the cotton research community has continued to assemble and refine the genomes for both wild and cultivated Gossypium species. With the accumulation of de novo genome assemblies and resequencing data across virous cotton populations, significant progress has been made in uncovering the genetic basis of key agronomic traits. Achieving the goal of cotton genomics-to-breeding (G2B) will require a deeper understanding of the spatiotemporal regulatory mechanisms involved in genome information storage and expression. We advocate for a cotton ENCODE project to systematically decode the functional elements and regulatory networks within the cotton genome. Technological advances, particularly on single-cell sequencing and high-resolution spatiotemporal omics, will be essential for elucidating these regulatory mechanisms. By integrating multi-omics data, genome editing tools, and artificial intelligence, these efforts will empower the genomics-driven strategies needed for future cotton G2B breeding.

Analysis of the genetic basis of fiber-related traits and flowering time in upland cotton using machine learning

Cotton is an important crop for fiber production, but the genetic basis underlying key agronomic traits, such as fiber quality and flowering days, remains complex. While machine learning (ML) has shown great potential in uncovering the genetic architecture of complex traits in other crops, its application in cotton has been limited. Here, we applied five machine learning models-AdaBoost, Gradient Boosting Regressor, LightGBM, Random Forest, and XGBoost-to identify loci associated with fiber quality and flowering days in cotton. We compared two SNP dataset down-sampling methods for model training and found that selecting SNPs with an Fscale value greater than 0 outperformed randomly selected SNPs in terms of model accuracy. We further performed machine learning quantitative trait loci (mlQTLs) analysis for 13 traits related to fiber quality and flowering days. These mlQTLs were then compared to those identified through genome-wide association studies (GWAS), revealing that the machine learning approach not only confirmed known loci but also identified novel QTLs. Additionally, we evaluated the effect of population size on model accuracy and found that larger population sizes resulted in better predictive performance. Finally, we proposed candidate genes for the identified mlQTLs, including two argonaute 5 proteins, Gh_A09G104100 and Gh_A09G104400, for the FL3/FS2 locus, as well as GhFLA17 and Syntaxin-121 (Gh_D09G143700) for the FSD09_2/FED09_2 locus. Our findings demonstrate the efficacy of machine learning in enhancing the identification of genetic loci in cotton, providing valuable insights for improving cotton breeding strategies.

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.

Natural variations in the Cis-elements of GhRPRS1 contributing to petal colour diversity in cotton

The cotton genus comprises both diploid and allotetraploid species, and the diversity in petal colour within this genus offers valuable targets for studying orthologous gene function differentiation and evolution. However, the genetic basis for this diversity in petal colour remains largely unknown. The red petal colour primarily comes from C, G, K, and D genome species, and it is likely that the common ancestor of cotton had red petals. Here, by employing a clone mapping strategy, we mapped the red petal trait to a specific region on chromosome A07 in upland cotton. Genomic comparisons and phylogenetic analyses revealed that the red petal phenotype introgressed from G. bickii. Transcriptome analysis indicated that GhRPRS1, which encodes a glutathione S-transferase, was the causative gene for the red petal colour. Knocking out GhRPRS1 resulted in white petals and the absence of red spots, while overexpression of both genotypes of GhRPRS1 led to red petals. Further analysis suggested that GhRPRS1 played a role in transporting pelargonidin-3-O-glucoside and cyanidin-3-O-glucoside. Promoter activity analysis indicated that variations in the promoter, but not in the gene body of GhRPRS1, have led to different petal colours within the genus. Our findings provide new insights into orthologous gene evolution as well as new strategies for modifying promoters in cotton breeding.

CRISPR-Cas9-mediated construction of a cotton CDPK mutant library for identification of insect-resistance genes

Calcium-dependent protein kinases (CDPKs) act as key signal transduction enzymes in plants, especially in response to diverse stresses, including herbivory. In this study, a comprehensive analysis of the CDPK gene family in upland cotton revealed that GhCPKs are widely expressed in multiple cotton tissues and respond positively to various biotic and abiotic stresses. We developed a strategy for screening insect-resistance genes from a CRISPR-Cas9 mutant library of GhCPKs. The library was created using 246 single-guide RNAs targeting the GhCPK gene family to generate 518 independent T0 plants. The average target-gene coverage was 86.18%, the genome editing rate was 89.49%, and the editing heritability was 82%. An insect bioassay in the field led to identification of 14 GhCPK mutants that are resistant or susceptible to insects. The mutant that showed the clearest insect resistance, cpk33/74 (in which the homologous genes GhCPK33 and GhCPK74 were knocked out), was selected for further study. Oral secretions from Spodoptera litura induced a rapid influx of Ca2+ in cpk33/74 leaves, resulting in a significant increase in jasmonic acid content. S-adenosylmethionine synthase is an important protein involved in plant stress response, and protein interaction experiments provided evidence for interactions of GhCPK33 and GhCPK74 with GhSAMS1 and GhSAM2. In addition, virus-induced gene silencing of GhSAMS1 and GhSAM2 in cotton impaired defense against S. litura. This study demonstrates an effective strategy for constructing a mutant library of a gene family in a polyploid plant species and offers valuable insights into the role of CDPKs in the interaction between plants and herbivorous insects.

Convergence and divergence of diploid and tetraploid cotton genomes

Polyploidy is an important driving force in speciation and evolution; however, the genomic basis for parallel selection of a particular trait between polyploids and ancestral diploids remains unexplored. Here we construct graph-based pan-genomes for diploid (A2) and allotetraploid (AD1) cotton species, enabled by an assembly of 50 genomes of genetically diverse accessions. We delineate a mosaic genome map of tetraploid cultivars that illustrates genomic contributions from semi-wild forms into modern cultivars. Pan-genome comparisons identify syntenic and hyper-divergent regions of continued variation between diploid and tetraploid cottons, and suggest an ongoing process of sequence evolution potentially linked to the contrasting genome size change in two subgenomes. We highlight 43% of genetic regulatory relationships for gene expression in diploid encompassing sequence divergence after polyploidy, and specifically characterize six underexplored convergent genetic loci contributing to parallel selection of fiber quality. This study offers a framework for pan-genomic dissection of genetic regulatory components underlying parallel selection of desirable traits in organisms.

Genetic Regulatory Perturbation of Gene Expression Impacted by Genomic Introgression in Fiber Development of Allotetraploid Cotton

Interspecific genomic introgression is an important evolutionary process with respect to the generation of novel phenotypic diversity and adaptation. A key question is how gene flow perturbs gene expression networks and regulatory interactions. Here, an introgression population of two species of allopolyploid cotton (Gossypium) to delineate the regulatory perturbations of gene expression regarding fiber development accompanying fiber quality change is utilized. De novo assembly of the recipient parent (G. hirsutum Emian22) genome allowed the identification of genomic variation and introgression segments (ISs) in 323 introgression lines (ILs) from the donor parent (G. barbadense 3-79). It documented gene expression dynamics by sequencing 1,284 transcriptomes of developing fibers and characterized genetic regulatory perturbations mediated by genomic introgression using a multi-locus model. Introgression of individual homoeologous genes exhibiting extreme low or high expression bias can lead to a parallel expression bias in their non-introgressed duplicates, implying a shared yet divergent regulatory fate of duplicated genes following allopolyploidy. Additionally, the IL N182 with improved fiber quality is characterized, and the candidate gene GhFLAP1 related to fiber length is validated. This study outlines a framework for understanding introgression-mediated regulatory perturbations in polyploids, and provides insights for targeted breeding of superior upland cotton fiber.

The genetic basis of leaf hair development in upland cotton (Gossypium hirsutum)

Trichomes, which originate from the epidermal cell of aerial organs, provide plants with defense and secretion functions. Although numerous genes have been implicated in trichome development, the molecular mechanisms underlying trichome cell formation in plants remain incompletely understood. Here, we using genome-wide association study (GWAS) across 1037 diverse accessions in upland cotton (Gossypium hirsutum) to identify three loci associated with leaf pubescence (hair) amount, located on chromosome A06 (LPA1), A08 (LPA2) and A11 (LPA3), respectively. GhHD1, a previously characterized candidate gene, was identified on LPA1 and encodes an HD-Zip transcription factor. For LPA2 and LPA3, we identified two candidate genes, GhGIR1 and GhGIR2, both encoding proteins with WD40 and RING domains that act as inhibitors of leaf hair formation. Expression analysis revealed that GhHD1 was predominantly expressed in hairy accessions, whereas GhGIR1 and GhGIR2 were expressed in hairless accessions. Silencing GhHD1 or overexpressing GhGIR1 in hairy accessions induced in a hairless phenotype, whereas silencing GhGIR2 in hairless accessions resulted in a hairy phenotype. We also demonstrated that GhHD1 interact with both GhGIR1 and GhGIR2, and GhGIR1 can interact with GhGIR2. Further investigation indicated that GhHD1 functions as a transcriptional activator, binding to the promoters of the GhGIR1 and GhGIR2 to active their expression, whereas GhGIR1 and GhGIR2 can suppress the transcriptional activation of GhHD1. Our findings shed light on the intricate regulatory network involving GhHD1, GhGIR1 and GhGIR2 in the initiation and development of plant epidermal hairs in cotton.

Transcriptome Analysis and Identification of Genes Associated with Cotton Seed Size

Cotton seeds, as the main by-product of cotton, are not only an important raw material for edible oil and feed but also a source of biofuel. The quality of cotton seeds directly affects cotton planting and is closely related to the yield and fiber quality. However, the molecular mechanism governing cotton seed size remains largely unexplored. This study investigates the regulatory mechanisms of cotton seed size by focusing on two cotton genotypes, N10 and N12, which exhibit notable phenotypic variations across multiple environments. Developing seeds were sampled at various stages (5, 20, 30, and 35 DPA) and subjected to RNA-seq. Temporal pattern clustering and WGCNA on differentially expressed genes identified 413 candidate genes, including these related to sugar metabolism that were significantly enriched in transcriptional regulation. A genetic transformation experiment indicated that the overexpression of the GhUXS5 gene encoding UDP-glucuronate decarboxylase 5 significantly increased seed size, suggesting an important role of GhUXS5 in regulating cotton seed size. This discovery provides crucial insights into the molecular mechanisms controlling cotton seed size, helping to unravel the complex regulatory network and offering new strategies and targets for cotton breeding to enhance the economic value of cotton seeds and overall cotton yield.

Nematode-resistance loci in upland cotton genomes are associated with structural differences

Reniform and root-knot nematode are two of the most destructive pests of conventional upland cotton, Gossypium hirsutum L., and continue to be a major threat to cotton fiber production in semiarid regions of the Southern United States and Central America. Fortunately, naturally occurring tolerance to these nematodes has been identified in the Pima cotton species (Gossypium barbadense) and several upland cotton varieties (G. hirsutum), which has led to a robust breeding program that has successfully introgressed and stacked these independent resistant traits into several upland cotton lineages with superior agronomic traits, e.g. BAR 32-30 and BARBREN-713. This work identifies the genomic variations of these nematode-tolerant accessions by comparing their respective genomes to the susceptible, high-quality fiber-producing parental line of this lineage: Phytogen 355 (PSC355). We discover several large genomic differences within marker regions that harbor putative resistance genes as well as expression mechanisms shared by the two resistant lines, with respect to the susceptible PSC355 parental line. This work emphasizes the utility of whole-genome comparisons as a means of elucidating large and small nuclear differences by lineage and phenotype.

Genome-wide association study of fiber quality traits in US upland cotton (Gossypium hirsutum L.)

KEY MESSAGE

A GWAS in an elite diversity panel, evaluated across 10 environments, identified genomic regions regulating six fiber quality traits, facilitating genomics-assisted breeding and gene discovery in upland cotton. In this study, an elite diversity panel of 348 upland cotton accessions was evaluated in 10 environments across the US Cotton Belt and genotyped with the cottonSNP63K array, for a genome-wide association study of six fiber quality traits. All fiber quality traits, upper half mean length (UHML: mm), fiber strength (FS: g tex-1), fiber uniformity (FU: %), fiber elongation (FE: %), micronaire (MIC) and short fiber content (SFC: %), showed high broad-sense heritability (> 60%). All traits except FE showed high genomic heritability. UHML, FS and FU were all positively correlated with each other and negatively correlated with FE, MIC and SFC. GWAS of these six traits identified 380 significant marker-trait associations (MTAs) including 143 MTAs on 30 genomic regions. These 30 genomic regions included MTAs identified in at least three environments, and 23 of them were novel associations. Phenotypic variation explained for the MTAs in these 30 genomic regions ranged from 6.68 to 11.42%. Most of the fiber quality-associated genomic regions were mapped in the D-subgenome. Further, this study confirmed the pleiotropic region on chromosome D11 (UHML, FS and FU) and identified novel co-localized regions on D04 (FU, SFC), D05 (UHML, FU, and D06 UHML, FU). Marker haplotype analysis identified superior combinations of fiber quality-associated genomic regions with high trait values (UHML = 32.34 mm; FS = 32.73 g tex-1; FE = 6.75%). Genomic analyses of traits, haplotype combinations and candidate gene information described in the current study could help leverage genetic diversity for targeted genetic improvement and gene discovery for fiber quality traits in cotton.

Terpene synthases GhTPS6 and GhTPS47 participate in resistance to Verticillium dahliae in upland cotton

Terpene synthases (TPSs) are enzymes responsible for catalyzing the production of diverse terpenes, the largest class of secondary metabolites in plants. Here, we identified 107 TPS gene loci encompassing 92 full-length TPS genes in upland cotton (Gossypium hirsutum L.). Phylogenetic analysis showed they were divided into six subfamilies. Segmental duplication and tandem duplication events contributed greatly to the expansion of TPS gene family, particularly the TPS-a and TPS-b subfamilies. Expression profile analysis screened out that GhTPSs may mediate the interaction between cotton and Verticillium dahliae. Three-dimensional structures and subcellular localizations of the two selected GhTPSs, GhTPS6 and GhTPS47, which belong to the TPS-a subfamily, demonstrated similarity in protein structures and nucleus and cytoplasm localization. Virus-induced gene silencing (VIGS) of the two GhTPSs yielded plants characterized by increased wilting and chlorosis, more severe vascular browning, and higher disease index than control plants. Additionally, knockdown of GhTPS6 and GhTPS47 led to the down-regulation of cotton terpene synthesis following V. dahliae infection, indicating that these two genes may positively regulate resistance to V. dahliae through the modulation of disease-resistant terpene biosynthesis. Overall, our study represents a comprehensive analysis of the G. hirsutum TPS gene family, revealing their potential roles in defense responses against Verticillium wilt.

GhRCD1 promotes cotton tolerance to cadmium by regulating the GhbHLH12-GhMYB44-GhHMA1 transcriptional cascade

Heavy metal pollution poses a significant risk to human health and wreaks havoc on agricultural productivity. Phytoremediation, a plant-based, environmentally benign, and cost-effective method, is employed to remove heavy metals from contaminated soil, particularly in agricultural or heavy metal-sensitive lands. However, the phytoremediation capacity of various plant species and germplasm resources display significant genetic diversity, and the mechanisms underlying these differences remain hitherto obscure. Given its potential benefits, genetic improvement of plants is essential for enhancing their uptake of heavy metals, tolerance to harmful levels, as well as overall growth and development in contaminated soil. In this study, we uncover a molecular cascade that regulates cadmium (Cd2+) tolerance in cotton, involving GhRCD1, GhbHLH12, GhMYB44, and GhHMA1. We identified a Cd2+-sensitive cotton T-DNA insertion mutant with disrupted GhRCD1 expression. Genetic knockout of GhRCD1 by CRISPR/Cas9 technology resulted in reduced Cd2+ tolerance in cotton seedlings, while GhRCD1 overexpression enhanced Cd2+ tolerance. Through molecular interaction studies, we demonstrated that, in response to Cd2+ presence, GhRCD1 directly interacts with GhbHLH12. This interaction activates GhMYB44, which subsequently activates a heavy metal transporter, GhHMA1, by directly binding to a G-box cis-element in its promoter. These findings provide critical insights into a novel GhRCD1-GhbHLH12-GhMYB44-GhHMA1 regulatory module responsible for Cd2+ tolerance in cotton. Furthermore, our study paves the way for the development of elite Cd2+-tolerant cultivars by elucidating the molecular mechanisms governing the genetic control of Cd2+ tolerance in cotton.

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.

Genome-Wide Identification of the WRKY Gene Family in Four Cotton Varieties and the Positive Role of GhWRKY31 in Response to Salt and Drought Stress

The WRKY gene family is ubiquitously distributed in plants, serving crucial functions in stress responses. Nevertheless, the structural organization and evolutionary dynamics of WRKY genes in cotton have not been fully elucidated. In this study, a total of 112, 119, 217, and 222 WRKY genes were identified in Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense, respectively. These 670 WRKY genes were categorized into seven distinct subgroups and unequally distributed across chromosomes. Examination of conserved motifs, domains, cis-acting elements, and gene architecture collectively highlighted the evolutionary conservation and divergence within the WRKY gene family in cotton. Analysis of synteny and collinearity further confirmed instances of expansion, duplication, and loss events among WRKY genes during cotton evolution. Furthermore, GhWRKY31 transgenic Arabidopsis exhibited heightened germination rates and longer root lengths under drought and salt stress. Silencing GhWRKY31 in cotton led to reduced levels of ABA, proline, POD, and SOD, along with downregulated expression of stress-responsive genes. Yeast one-hybrid and molecular docking assays confirmed the binding capacity of GhWRKY31 to the W box of GhABF1, GhDREB2, and GhRD29. The findings collectively offer a systematic and comprehensive insight into the evolutionary patterns of cotton WRKYs, proposing a suitable regulatory framework for developing cotton cultivars with enhanced resilience to drought and salinity stress.

Genomic insights into local adaptation of upland cotton in China and Pakistan

KEY MESSAGE

Different kinship and resistance to cotton leaf curl disease (CLCuD) and heat were found between upland cotton cultivars from China and Pakistan. 175 SNPs and 82 InDels loci related to yield, fiber quality, CLCuD, and heat resistance were identified. Elite alleles found in Pakistani accessions aided local adaptation to climatic condition of two countries. Adaptation of upland cotton (Gossypium hirsutum) beyond its center of origin is expected to be driven by tailoring of the genome and genes to enhance yield and quality in new ecological niches. Here, resequencing of 456 upland cotton accessions revealed two distinct kinships according to the associated country. Fiber quality and lint percentage were consistent across kinships, but resistance to cotton leaf curl disease (CLCuD) and heat was distinctly exhibited by accessions from Pakistan, illustrating highly local adaption. A total of 175 SNP and 82 InDel loci related to yield, fiber quality, CLCuD and heat resistance were identified; among them, only two overlapped between Pakistani and Chinese accessions underscoring the divergent domestication and improvement targets in each country. Loci associated with resistance alleles to leaf curl disease and high temperature were largely found in Pakistani accessions to counter these stresses prevalent in Pakistan. These results revealed that breeding activities led to the accumulation of unique alleles and helped upland cotton become adapted to the respective climatic conditions, which will contribute to elucidating the genetic mechanisms that underlie resilience traits and help develop climate-resilient cotton cultivars for use worldwide.

Genomic loci associated with leaf abscission contribute to machine picking and environmental adaptability in upland cotton (Gossypium hirsutum L.)

INTRODUCTION

Defoliation by applying defoliants before machine picking is an important agricultural practice that enhances harvesting efficiency and leads to increased raw cotton purity. However, the fundamental characteristics of leaf abscission and the underlying genetic basis in cotton are not clearly understood.

OBJECTIVES

In this study, we aimed to (1) reveal the phenotypic variations in cotton leaf abscission, (2) discover the whole-genome differentiation sweeps and genetic loci related to defoliation, (3) identify and verify the functions of key candidate genes associated with defoliation, and (4) explore the relationship between haplotype frequency of loci and environmental adaptability.

METHODS

Four defoliation-related traits of 383 re-sequenced Gossypium hirsutum accessions were investigated in four environments. The genome-wide association study (GWAS), linkage disequilibrium (LD) interval genotyping and functional identification were conducted. Finally, the haplotype variation related to environmental adaptability and defoliation traits was revealed.

RESULTS

Our findings revealed the fundamental phenotypic variations of defoliation traits in cotton. We showed that defoliant significantly increased the defoliation rate without incurring yield and fiber quality penalties. The strong correlations between defoliation traits and growth period traits were observed. A genome-wide association study of defoliation traits identified 174 significant SNPs. Two loci (RDR7 on A02 and RDR13 on A13) that significantly associated with the relative defoliation rate were described, and key candidate genes GhLRR and GhCYCD3;1, encoding a leucine-rich repeat (LRR) family protein and D3-type cell cyclin 1 protein respectively, were functional verified by expression pattern analysis and gene silencing. We found that combining of two favorable haplotypes (HapRDR7 and HapRDR13) improved sensitivity to defoliant. The favorable haplotype frequency generally increased in high latitudes in China, enabling adaptation to the local environment.

CONCLUSION

Our findings lay an important foundation for the potentially broad application of leveraging key genetic loci in breeding machine-pickable cotton.

Identification of chromosomes by fluorescence in situ hybridization in Gossypium hirsutum via developing oligonucleotide probes

Discrimination of chromosome is essential for chromosome manipulation or visual chromosome characterization. Oligonucleotide probes can be employed to simplify the procedures of chromosome identification in molecular cytogenetics due to its simplicity, fastness, cost-effectiveness, and high efficiency. So far, however, visual identification of cotton chromosomes remains unsolved. Here, we developed 16 oligonucleotide probes for rapid and accurate identification of chromosomes in Gossypium hirsutum: 9 probes, of which each is able to distinguish individually one pair of chromosomes, and seven probes, of which each distinguishes multiple pairs of chromosomes. Besides the identification of Chrs. A09 and D09, we first find Chr. D08, which carries both 45S and 5S rDNA sequences. Interestingly, we also find Chr. A07 has a small 45S rDNA size, suggesting that the size of this site on Chr. A07 may have reduced during evolution. By the combination of 45S and 5S rDNA sequences and oligonucleotide probes developed, 10 chromosomes (Chrs. 3-7, and 9-13) in A subgenome and 7 (Chrs. 1-2, 4-5, and 7-9) in D subgenome of cotton are able to be recognized. This study establishes cotton oligonucleotide fluorescence in situ hybridization technology for discrimination of chromosomes, which supports and guides for sequence assembling, particularly, for tandem repeat sequences in cotton.

GhWRKY4 binds to the histone deacetylase GhHDA8 promoter to regulate drought and salt tolerance in Gossypium hirsutum

Soil drought and salinization, caused by water deficiency, have become the greatest concerns limiting crop production. Up to now, the WRKY transcription factor and histone deacetylase have been shown to be involved in drought and salt responses. However, the molecular mechanism underlying their interaction remains unclear in cotton. Herein, we identified GhWRKY4, a member of WRKY gene family, which is induced by drought and salt stress and is located in the nucleus. The ectopic expression of GhWRKY4 in Arabidopsis enhanced drought and salt tolerance, and suppressing GhWRKY4 in cotton increased susceptibility to drought and salinity. Subsequently, DAP-seq analysis revealed that the W box element in the promoter of stress-induced genes could potentially be the binding target for GhWRKY4 protein. GhWRKY4 binds to the promoters of GhHDA8 and GhNHX7 via W box element, and the expression level of GhHDA8 was increased in GhWRKY4-silenced plants. In addition, GhHDA8-overexpressed Arabidopsis were found to be hypersensitive to drought and salt stress, while silencing of GhHDA8 enhanced drought and salt tolerance in cotton. The stress-related genes, such as GhDREB2A, GhRD22, GhP5CS, and GhNHX7, were induced in GhHDA8-silenced plants. Our findings indicate that the GhWRKY4-GhHDA8 module regulates drought and salt tolerance in cotton. Collectively, the results provide new insights into the coordination of transcription factors and histone deacetylases in regulating drought and salt stress responses in plants.

High-quality Gossypium hirsutum and Gossypium barbadense genome assemblies reveal the landscape and evolution of centromeres

Centromere positioning and organization are crucial for genome evolution; however, research on centromere biology is largely influenced by the quality of available genome assemblies. Here, we combined Oxford Nanopore and Pacific Biosciences technologies to de novo assemble two high-quality reference genomes for Gossypium hirsutum (TM-1) and Gossypium barbadense (3-79). Compared with previously published reference genomes, our assemblies show substantial improvements, with the contig N50 improved by 4.6-fold and 5.6-fold, respectively, and thus represent the most complete cotton genomes to date. These high-quality reference genomes enable us to characterize 14 and 5 complete centromeric regions for G. hirsutum and G. barbadense, respectively. Our data revealed that the centromeres of allotetraploid cotton are occupied by members of the centromeric repeat for maize (CRM) and Tekay long terminal repeat families, and the CRM family reshapes the centromere structure of the At subgenome after polyploidization. These two intertwined families have driven the convergent evolution of centromeres between the two subgenomes, ensuring centromere function and genome stability. In addition, the repositioning and high sequence divergence of centromeres between G. hirsutum and G. barbadense have contributed to speciation and centromere diversity. This study sheds light on centromere evolution in a significant crop and provides an alternative approach for exploring the evolution of polyploid plants.

Dissecting the major genetic components underlying cotton lint development

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

GhTPPA_2 enhancement of tobacco sugar accumulation and drought tolerance

Trehalose metabolism plays an important role in plant growth and response to abiotic stress. Trehalose-6-phosphate (Tre6P) can help regulate sugar homeostasis and act as an indication signal for intracellular sugar levels. Crop productivity can be greatly increased by altering the metabolic level of endogenous trehalose in plants, which can optimize the source-sink connection. In this study, the upland cotton GhTPP protein family was first homologously compared and 24 GhTPP genes were found. Transcriptome analysis revealed that GhTPP members had obvious tissue expression specificity. Among them, GhTPPA_2 (Gh_A12G223300.1) was predominantly expressed in leaves and bolls. The results of subcellular localization showed that GhTPPA_2 is localized in the chloroplast. Via PlantCare, we analyzed the promoters and found that the expression of GhTPPA_2 may be induced by light, abiotic stress, and hormones such as abscisic acid, ethylene, salicylic acid and jasmonic acid. In addition, GhTPPA_2 was overexpressed (TPPAoe) in tobacco, and we found that the TPPase activity of TPPAoe tobacco increased by 66 %. Soluble sugar content increased by 39 % and starch content increased by 27 %. Whereas, the transgenic tobacco had obvious growth advantages under 100 mM mannitol stress. Transcriptome sequencing results showed that the differential genes between TPPAoe and control were considerably enriched in functions related to photosynthesis, phosphate group metabolism, and carbohydrate metabolism. This study shows that GhTPPA_2 is involved in regulating sugar metabolism, improving soluble sugar accumulation and drought stress tolerance of tobacco, which provides theoretical basis for research on high yield and drought tolerance of crops.

Genomic and epigenomic insights into the mechanism of cold response in upland cotton (Gossypium hirsutum)

Functional genome research, including gene transcriptional and posttranslational modifications of histones, can benefit greatly from a high-quality genome assembly. Histone modification plays a significant role in modulating the responses to abiotic stress in plants. However, there are limited reports on the involvement of dynamic changes in histone modification in cold stress response in upland cotton. In this study, the genome of an elite accession, YM11, with considerable cold stress tolerance was de novo assembled, which yielded a genome of 2343.06 Mb with a contig N50 of 88.96 Mb, and a total of 73,821 protein-coding gene models were annotated. Comparisons among YM11 and five Gossypium allopolyploid cotton assemblies highlighted a large amount of structural variations and presence/absence variations. We analyzed transcriptome and metabolome changes in YM11 seedlings subjected to cold stress. Using the CUT&Tag method, genome-wide H3K4me3 and H3K9ac modification patterns and effect of histone changes on gene expression were profiled during cold stress. Significant and consistently changing histone modifications and the gene expressions were screened, of which transcription factors (TFs) were highlighted. Our results suggest a positive correlation between the changes in H3K4me3, H3K9ac modifications and cold stress-responsive gene activation. This genome assembly and comprehensive analysis of genome-wide histone modifications and gene expression provide insights into the genomic variation and epigenetic responses to cold stress in upland cotton.

Cotton pedigree genome reveals restriction of cultivar-driven strategy in cotton breeding

BACKGROUND

Many elite genes have been identified from the available cotton genomic data, providing various genetic resources for gene-driven breeding. However, backbone cultivar-driven breeding is the most widely applied strategy. Revealing the genetic basis of cultivar-driven strategy's restriction is crucial for transition of cotton breeding strategy.

RESULT

CRI12 is a backbone cultivar in cultivar-driven breeding. Here we sequence the pedigree of CRI12 using Nanopore long-read sequencing. We construct a graphical pedigree genome using the high-quality CRI12 genome and 13,138 structural variations within 20 different pedigree members. We find that low hereditary stability of elite segments in backbone cultivars is a drawback of cultivar-driven strategy. We also identify 623 functional segments in CRI12 for multiple agronomic traits in presence and absence variation-based genome-wide association study on three cohorts. We demonstrate that 25 deleterious segments are responsible for the geographical divergence of cotton in pathogen resistance. We also characterize an elite pathogen-resistant gene (GhKHCP) utilized in modern cotton breeding. In addition, we identify 386 pedigree fingerprint segments by comparing the segments of the CRI12 pedigree with those of a large cotton population.

CONCLUSION

We characterize the genetic patterns of functional segments in the pedigree of CRI12 using graphical genome method, revealing restrictions of cultivar-driven strategies in cotton breeding. These findings provide theoretical support for transitioning from cultivar-driven to gene-driven strategy in cotton breeding.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1

BACKGROUND

Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear.

RESULTS

In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR).

CONCLUSIONS

Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

Genetic diversity and population structure analyses and genome-wide association studies of photoperiod sensitivity in cotton (Gossypium hirsutum L.)

KEY MESSAGE

Genetic diversity and population structure analyses showed progressively narrowed diversity in US Upland cotton compared to land races. GWAS identified genomic regions and candidate genes for photoperiod sensitivity in cotton. Six hundred fifty-seven accessions that included elite cotton germplasm (DIV panel), lines of a public cotton breeding program (FB panel), and tropical landrace accessions (TLA panel) of Gossypium hirsutum L. were genotyped with cottonSNP63K array and phenotyped for photoperiod sensitivity under long day-length conditions. The genetic diversity analysis using 26,952 polymorphic SNPs indicated a progressively narrowed diversity from the landraces (0.230) to the DIV panel accessions (0.195) and FB panel (0.116). Structure analysis in the US germplasm identified seven subpopulations representing all four major regions of the US cotton belt. Three subpopulations were identified within the landrace accessions. The highest fixation index (FST) of 0.65 was found between landrace accessions of Guatemala and the Plains-type cultivars from Southwest cotton region while the lowest FST values were between the germplasms of Mid-South and Southeastern regions. Genome wide association studies (GWAS) of photoperiod response using 600 phenotyped accessions identified 14 marker trait associations spread across eight Upland cotton chromosomes. Six of these marker trait associations, on four chromosomes (A10, D04, D05, and D06), showed significant epistatic interactions. Targeted genomic analysis identified regions with 19 candidate genes including Transcription factor Vascular Plant One-Zinc Finger 1 (VOZ1) and Protein Photoperiod-Independent Early Flowering 1 (PIE1) genes. Genetic diversity and genome wide analyses of photoperiod sensitivity in diverse cotton germplasms will enable the use of genomic tools to systematically utilize the tropical germplasm and its beneficial alleles for broadening the genetic base in Upland cotton.

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.

Detection of epistasis interaction loci for fiber quality-related trait via 3VmrMLM in upland cotton

Cotton fiber quality-related traits, such as fiber length, fiber strength, and fiber elongation, are affected by complex mechanisms controlled by multiple genes. Determining the QTN-by-QTN interactions (QQIs) associated with fiber quality-related traits is therefore essential for accelerating the genetic enhancement of cotton breeding. In this study, a natural population of 1,245 upland cotton varieties with 1,122,352 SNPs was used for detecting the main-effect QTNs and QQIs using the 3V multi-locus random-SNP-effect mixed linear model (3VmrMLM) method. A total of 171 significant main-effect QTNs and 42 QQIs were detected, of which 22 were both main-effect QTNs and QQIs. Of the detected 42 QQIs, a total of 13 significant loci and 5 candidate genes were reported in previous studies. Among the three interaction types, the AD interaction type has a preference for the trait of FE. Additionally, the QQIs have a substantial impact on the enhancement predictability for fiber quality-related traits. The study of QQIs is crucial for elucidating the genetic mechanism of cotton fiber quality and enhancing breeding efficiency.

Quantitative trait loci and candidate genes for yield-related traits of upland cotton revealed by genome-wide association analysis under drought conditions

BACKGROUND

Due to the influence of extreme weather, the environment in China's main cotton-producing areas is prone to drought stress conditions, which affect the growth and development of cotton and lead to a decrease in cotton yield.

RESULTS

In this study, 188 upland cotton germplasm resources were phenotyped for data of 8 traits (including 3 major yield traits) under drought conditions in three environments for two consecutive years. Correlation analysis revealed significant positive correlations between the three yield traits. Genetic analysis showed that the estimated heritability of the seed cotton index (SC) under drought conditions was the highest (80.81%), followed by that of boll weight (BW) (80.64%) and the lint cotton index (LC) (70.49%) With genome-wide association study (GWAS) analysis, a total of 75 quantitative trait loci (QTLs) were identified, including two highly credible new QTL hotspots. Three candidate genes (Gh_D09G064400, Gh_D10G261000 and Gh_D10G254000) located in the two new QTL hotspots, QTL51 and QTL55, were highly expressed in the early stage of fiber development and showed significant correlations with SC, LC and BW. The expression of three candidate genes in two extreme materials after drought stress was analyzed by qRT-PCR, and the expression of these two materials in fibers at 15, 20 and 25 DPA. The expression of these three candidate genes was significantly upregulated after drought stress and was significantly higher in drought-tolerant materials than in drought-sensitive materials. In addition, the expression levels of the three candidate genes were higher in the early stage of fiber development (15 DPA), and the expression levels in drought-tolerant germplasm were higher than those in drought-sensitive germplasm. These three candidate genes may play an important role in determining cotton yield under drought conditions.

CONCLUSIONS

This study is helpful for understanding the regulatory genes affecting cotton yield under drought conditions and provides germplasm and candidate gene resources for breeding high-yield cotton varieties under these conditions.

Single-cell resolution analysis reveals the preparation for reprogramming the fate of stem cell niche in cotton lateral meristem

BACKGROUND

Somatic embryogenesis is a major process for plant regeneration. However, cell communication and the gene regulatory network responsible for cell reprogramming during somatic embryogenesis are still largely unclear. Recent advances in single-cell technologies enable us to explore the mechanism of plant regeneration at single-cell resolution.

RESULTS

We generate a high-resolution single-cell transcriptomic landscape of hypocotyl tissue from the highly regenerable cotton genotype Jin668 and the recalcitrant TM-1. We identify nine putative cell clusters and 23 cluster-specific marker genes for both cultivars. We find that the primary vascular cell is the major cell type that undergoes cell fate transition in response to external stimulation. Further developmental trajectory and gene regulatory network analysis of these cell clusters reveals that a total of 41 hormone response-related genes, including LAX2, LAX1, and LOX3, exhibit different expression patterns in the primary xylem and cambium region of Jin668 and TM-1. We also identify novel genes, including CSEF, PIS1, AFB2, ATHB2, PLC2, and PLT3, that are involved in regeneration. We demonstrate that LAX2, LAX1 and LOX3 play important roles in callus proliferation and plant regeneration by CRISPR/Cas9 editing and overexpression assay.

CONCLUSIONS

This study provides novel insights on the role of the regulatory network in cell fate transition and reprogramming during plant regeneration driven by somatic embryogenesis.

Comprehensive Evaluation and Transcriptome Analysis Reveal the Salt Tolerance Mechanism in Semi-Wild Cotton (Gossypium purpurascens)

Elevated salinity significantly threatens cotton growth, particularly during the germination and seedling stages. The utilization of primitive species of Gossypium hirsutum, specifically Gossypium purpurascens, has the potential to facilitate the restoration of genetic diversity that has been depleted due to selective breeding in modern cultivars. This investigation evaluated 45 G. purpurascens varieties and a salt-tolerant cotton variety based on 34 morphological, physiological, and biochemical indicators and comprehensive salt tolerance index values. This study effectively identified a total of 19 salt-tolerant and two salt-resistant varieties. Furthermore, transcriptome sequencing of a salt-tolerant genotype (Nayanmian-2; NY2) and a salt-sensitive genotype (Sanshagaopao-2; GP2) revealed 2776, 6680, 4660, and 4174 differentially expressed genes (DEGs) under 0.5, 3, 12, and 24 h of salt stress. Gene ontology enrichment analysis indicated that the DEGs exhibited significant enrichment in biological processes like metabolic (GO:0008152) and cellular (GO:0009987) processes. MAPK signaling, plant-pathogen interaction, starch and sucrose metabolism, plant hormone signaling, photosynthesis, and fatty acid metabolism were identified as key KEGG pathways involved in salinity stress. Among the DEGs, including NAC, MYB, WRKY, ERF, bHLH, and bZIP, transcription factors, receptor-like kinases, and carbohydrate-active enzymes were crucial in salinity tolerance. Weighted gene co-expression network analysis (WGCNA) unveiled associations of salt-tolerant genotypes with flavonoid metabolism, carbon metabolism, and MAPK signaling pathways. Identifying nine hub genes (MYB4, MYB105, MYB36, bZIP19, bZIP43, FRS2 SMARCAL1, BBX21, F-box) across various intervals offered insights into the transcriptional regulation mechanism of salt tolerance in G. purpurascens. This study lays the groundwork for understanding the important pathways and gene networks in response to salt stress, thereby providing a foundation for enhancing salt tolerance in upland cotton.

Acetylation of GhCaM7 enhances cotton resistance to Verticillium dahliae

Protein lysine acetylation is an important post-translational modification mechanism involved in cellular regulation in eukaryotes. Calmodulin (CaM) is a ubiquitous Ca2+ sensor in eukaryotes and is crucial for plant immunity, but it is so far unclear whether acetylation is involved in CaM-mediated plant immunity. Here, we found that GhCaM7 is acetylated upon Verticillium dahliae (V. dahliae) infection and a positive regulator of V. dahliae resistance. Overexpressing GhCaM7 in cotton and Arabidopsis enhances V. dahliae resistance and knocking-down GhCaM7 makes cotton more susceptible to V. dahliae. Transgenic Arabidopsis plants overexpressing GhCaM7 with mutation at the acetylation site are more susceptible to V. dahliae than transgenics overexpressing the wild-type GhCaM7, implying the importance of the acetylated GhCaM7 in response to V. dahliae infection. Yeast two-hybrid, bimolecular fluorescent complementation, luciferase complementation imaging, and coimmunoprecipitation assays demonstrated interaction between GhCaM7 and an osmotin protein GhOSM34 that was shown to have a positive role in V. dahliae resistance. GhCaM7 and GhOSM34 are co-localized in the cell membrane. Upon V. dahliae infection, the Ca2+ content reduces almost instantly in plants with downregulated GhCaM7 or GhOSM34. Down regulating GhOSM34 enhances accumulation of Na+ and increases cell osmotic pressure. Comparative transcriptomic analyses between cotton plants with an increased or reduced expression level of GhCaM7 and wild-type plants indicate the involvement of jasmonic acid signaling pathways and reactive oxygen species in GhCaM7-enabled disease resistance. Together, these results demonstrate the involvement of CaM protein in the interaction between cotton and V. dahliae, and more importantly, the involvement of the acetylated CaM in the interaction.

High-quality genome assembly of Verticillium dahliae VD991 allows for screening and validation of pathogenic genes

Verticillium dahliae (V. dahliae) is a notorious soil-borne pathogen causing Verticillium wilt in more than 400 dicotyledonous plants, including a wide range of economically important crops, such as cotton, tomato, lettuce, potato, and romaine lettuce, which can result in extensive economic losses. In the last decade, several studies have been conducted on the physiological and molecular mechanisms of plant resistance to V. dahliae. However, the lack of a complete genome sequence with a high-quality assembly and complete genomic annotations for V. dahliae has limited these studies. In this study, we produced a full genomic assembly for V. dahliae VD991 using Nanopore sequencing technology, consisting of 35.77 Mb across eight pseudochromosomes and with a GC content of 53.41%. Analysis of the genome completeness assessment (BUSCO alignment: 98.62%; Illumina reads alignment: 99.17%) indicated that our efforts resulted in a nearly complete and high-quality genomic assembly. We selected 25 species closely related to V. dahliae for evolutionary analysis, confirming the evolutionary relationship between V. dahliae and related species, and the identification of a possible whole genome duplication event in V. dahliae. The interaction between cotton and V. dahliae was investigated by transcriptome sequencing resulting in the identification of many genes and pathways associated with cotton disease resistance and V. dahliae pathogenesis. These results will provide new insights into the pathogenic mechanisms of V. dahliae and contribute to the cultivation of cotton varieties resistant to Verticillium wilt.

Identification and characterization of the LDAP family revealed GhLDAP2_Dt enhances drought tolerance in cotton

Lipid droplet-associated proteins (LDAPs) play essential roles in tissue growth and development and in drought stress responses in plants. Cotton is an important fiber and cash crop; however, the LDAP family has not been characterized in cotton. In this study, a total of 14, six, seven, and seven genes were confirmed as LDAP family members in Gossypium hirsutum, Gossypium raimondii, Gossypium arboreum, and Gossypium stocksii, respectively. Additionally, expansion in the LDAP family occurred with the formation of Gossypium, which is mirrored in the number of LDAPs found in five Malvaceae species (Gossypioides kirkii, Bombax ceiba, Durio zibethinus, Theobroma cacao, and Corchorus capsularis), Arabidopsis thaliana, and Carica papaya. The phylogenetic tree showed that the LDAP genes in cotton can be divided into three groups (I, II, and III). The analysis of gene structure and conserved domains showed that LDAPs derived from group I (LDAP1/2/3) are highly conserved during evolution, while members from groups II and III had large variations in both domains and gene structures. The gene expression pattern analysis of LDAP genes showed that they are expressed not only in the reproductive organs (ovule) but also in vegetative organs (root, stem, and leaves). The expression level of two genes in group III, GhLDAP6_At/Dt, were significantly higher in fiber development than in other tissues, indicating that it may be an important regulator of cotton fiber development. In group III, GhLDAP2_At/Dt, especially GhLDAP2_Dt was strongly induced by various abiotic stresses. Decreasing the expression of GhLDAP2_Dt in cotton via virus-induced gene silencing increased the drought sensitivity, and the over-expression of GhLDAP2_Dt led to increased tolerance to mannitol-simulated osmotic stress at the germination stage. Thus, we conclude that GhLDAP2_Dt plays a positive role in drought tolerance.

Genome-wide association analysis reveals a novel pathway mediated by a dual-TIR domain protein for pathogen resistance in cotton

BACKGROUND

Verticillium wilt is one of the most devasting diseases for many plants, leading to global economic loss. Cotton is known to be vulnerable to its fungal pathogen, Verticillium dahliae, yet the related genetic mechanism remains unknown.

RESULTS

By genome-wide association studies of 419 accessions of the upland cotton, Gossypium hirsutum, we identify ten loci that are associated with resistance against Verticillium wilt. Among these loci, SHZDI1/SHZDP2/AYDP1 from chromosome A10 is located on a fragment introgressed from Gossypium arboreum. We characterize a large cluster of Toll/interleukin 1 (TIR) nucleotide-binding leucine-rich repeat receptors in this fragment. We then identify a dual-TIR domain gene from this cluster, GhRVD1, which triggers an effector-independent cell death and is induced by Verticillium dahliae. We confirm that GhRVD1 is one of the causal gene for SHZDI1. Allelic variation in the TIR domain attenuates GhRVD1-mediated resistance against Verticillium dahliae. Homodimerization between TIR1-TIR2 mediates rapid immune response, while disruption of its αD- and αE-helices interface eliminates the autoactivity and self-association of TIR1-TIR2. We further demonstrate that GhTIRP1 inhibits the autoactivity and self-association of TIR1-TIR2 by competing for binding to them, thereby preventing the resistance to Verticillium dahliae.

CONCLUSIONS

We propose the first working model for TIRP1 involved self-association and autoactivity of dual-TIR domain proteins that confer compromised pathogen resistance of dual-TIR domain proteins in plants. The findings reveal a novel mechanism on Verticillium dahliae resistance and provide genetic basis for breeding in future.

Genome-wide identification of terpenoid synthase family genes in Gossypium hirsutum and functional dissection of its subfamily cadinene synthase A in gossypol synthesis

Plant terpenoid synthase (TPS) family genes participate in metabolite synthesis, hormones, gossypol, etc. Here, we genome-widely identified TPS family genes in 12 land plant species. Four hundred and thirty TPS-related genes were divided into seven subfamilies. The TPS-c in Bryophytes was suggested to be the earliest subfamily, followed by the TPS-e/f and TPS-h presence in ferns. TPS-a, the largest number of genes, was derived from monocotyledonous and dicotyledonous plants. Collinearity analysis showed that 38 out of the 76 TPS genes in G. hirsutum were collinear within G. arboreum and G. raimondii. Twenty-one GhTPS-a genes belong to the cadinene synthase (GhCDN) subfamily and were divided into five groups, A, B, C, D, and E. The special cis-elements in the promoters of 12 GhCDN-A genes suggested that the JA and ethylene signaling pathways may be involved in their expression regulation. When 12 GhCDN-A genes were simultaneously silenced through virus-induced gene silencing, the glandular color of GhCDN-A-silenced plants was lighter than that of the control, supported by a gossypol content decrease based on HPLC testing, suggesting that GhCDN-A subgroup genes participate in gossypol synthesis. According to RNA-seq analysis, gossypol synthesis-related genes and disease-resistant genes in the glandular variety exhibited upregulated expression compared to the glandless variety, whereas hormone signaling-related genes were downregulated. All in all, these results revealed plant TPS gene evolution rules and dissected the TPS subfamily, GhCDN-A, function in gossypol synthesis in cotton.

Structural variation (SV)-based pan-genome and GWAS reveal the impacts of SVs on the speciation and diversification of allotetraploid cottons

Structural variations (SVs) have long been described as being involved in the origin, adaption, and domestication of species. However, the underlying genetic and genomic mechanisms are poorly understood. Here, we report a high-quality genome assembly of Gossypium barbadense acc. Tanguis, a landrace that is closely related to formation of extra-long-staple (ELS) cultivated cotton. An SV-based pan-genome (Pan-SV) was then constructed using a total of 182 593 non-redundant SVs, including 2236 inversions, 97 398 insertions, and 82 959 deletions from 11 assembled genomes of allopolyploid cotton. The utility of this Pan-SV was then demonstrated through population structure analysis and genome-wide association studies (GWASs). Using segregation mapping populations produced through crossing ELS cotton and the landrace along with an SV-based GWAS, certain SVs responsible for speciation, domestication, and improvement in tetraploid cottons were identified. Importantly, some of the SVs presently identified as associated with the yield and fiber quality improvement had not been identified in previous SNP-based GWAS. In particular, a 9-bp insertion or deletion was found to associate with elimination of the interspecific reproductive isolation between Gossypium hirsutum and G. barbadense. Collectively, this study provides new insights into genome-wide, gene-scale SVs linked to important agronomic traits in a major crop species and highlights the importance of SVs during the speciation, domestication, and improvement of cultivated crop species.