Single-Cell Transcriptome Atlas and Regulatory Dynamics in Developing Cotton Anthers.

Plant anthers are composed of different specialized cell types with distinct roles in plant reproduction. High temperature (HT) stress causes male sterility, resulting in crop yield reduction. However, the spatial expression atlas and regulatory dynamics during anther development and in response to HT remain largely unknown. Here, the first single-cell transcriptome atlas and chromatin accessibility survey in cotton anther are established, depicting the specific expression and epigenetic landscape of each type of cell in anthers. The reconstruction of meiotic cells, tapetal cells, and middle layer cell developmental trajectories not only identifies novel expressed genes, but also elucidates the precise degradation period of middle layer and reveals a rapid function transition of tapetal cells during the tetrad stage. By applying HT, heterogeneity in HT response is shown among cells of anthers, with tapetal cells responsible for pollen wall synthesis are most sensitive to HT. Specifically, HT shuts down the chromatin accessibility of genes specifically expressed in the tapetal cells responsible for pollen wall synthesis, such as QUARTET 3 (QRT3) and CYTOCHROME P450 703A2 (CYP703A2), resulting in a silent expression of these genes, ultimately leading to abnormal pollen wall and male sterility. Collectively, this study provides substantial information on anthers and provides clues for heat-tolerant crop creation.

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

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

Histone H3 lysine 27 trimethylation suppresses jasmonate biosynthesis and signaling to affect male fertility under high temperature in cotton.

High-temperature (HT) stress causes male sterility in crops, thus decreasing yields. To explore the possible contribution of histone modifications to male fertility under HT conditions, we defined the histone methylation landscape for the marks histone H3 lysine 27 trimethylation (H3K27me3) and histone H3 lysine 4 trimethylation (H3K4me3) by chromatin immunoprecipitation sequencing (ChIP-seq) in two differing upland cotton (Gossypium hirsutum) varieties. We observed a global disruption in H3K4me3 and H3K27me3 modifications, especially H3K27me3, in cotton anthers subjected to HT. HT affected the bivalent H3K4me3-H3K27me3 modification more than either monovalent modification. We determined that removal of H3K27me3 at the promoters of jasmonate-related genes increased their expression, maintaining male fertility under HT in the HT-tolerant variety at the anther dehiscence stage. Modulating jasmonate homeostasis or signaling resulted in an anther indehiscence phenotype under HT. Chemical suppression of H3K27me3 deposition increased jasmonic acid contents and maintained male fertility under HT. In summary, our study provides new insights into the regulation of male fertility by histone modifications under HT and suggests a potential strategy for improving cotton HT tolerance.

N6-methyladenosine RNA modification regulates cotton drought response in a Ca2+ and ABA-dependent manner.

N6 -methyladenosine (m6 A) is the most prevalent internal modification present in mRNAs, and is considered to participate in a range of developmental and biological processes. Drought response is highly regulated at the genomic, transcriptional and post-transcriptional levels. However, the biological function and regulatory mechanism of m6 A modification in the drought stress response is still poorly understood. We generated a transcriptome-wide m6 A map using drought-resistant and drought-sensitive varieties of cotton under different water deficient conditions to uncover patterns of m6 A methylation in cotton response to drought stress. The results reveal that m6 A represents a common modification and exhibit dramatic changes in distribution during drought stress. More 5'UTR m6 A was deposited in the drought-resistant variety and was associated with a positive effect on drought resistance by regulating mRNA abundance. Interestingly, we observed that increased m6 A abundance was associated with increased mRNA abundance under drought, contributing to drought resistance, and vice versa. The demethylase GhALKBH10B was found to decrease m6 A levels, facilitating the mRNA decay of ABA signal-related genes (GhZEP, GhNCED4 and GhPP2CA) and Ca2+ signal-related genes (GhECA1, GhCNGC4, GhANN1 and GhCML13), and mutation of GhALKBH10B enhanced drought resistance at seedling stage in cotton. Virus-induced gene silencing (VIGS) of two Ca2+ -related genes, GhECA1 and GhCNGC4, reduced drought resistance with the decreased m6 A enrichment on silenced genes in cotton. Collectively, we reveal a novel mechanism of post-transcriptional modification involved in affecting drought response in cotton, by mediating m6 A methylation on targeted transcripts in the ABA and Ca2+ signalling transduction pathways.

High-temperature stress in crops: male sterility, yield loss and potential remedy approaches.

Global food security is one of the utmost essential challenges in the 21st century in providing enough food for the growing population while coping with the already stressed environment. High temperature (HT) is one of the main factors affecting plant growth, development and reproduction and causes male sterility in plants. In male reproductive tissues, metabolic changes induced by HT involve carbohydrates, lipids, hormones, epigenetics and reactive oxygen species, leading to male sterility and ultimately reducing yield. Understanding the mechanism and genes involved in these pathways during the HT stress response will provide a new path to improve crops by using molecular breeding and biotechnological approaches. Moreover, this review provides insight into male sterility and integrates this with suggested strategies to enhance crop tolerance under HT stress conditions at the reproductive stage.

Genomic innovation and regulatory rewiring during evolution of the cotton genus Gossypium.

Phenotypic diversity and evolutionary innovation ultimately trace to variation in genomic sequence and rewiring of regulatory networks. Here, we constructed a pan-genome of the Gossypium genus using ten representative diploid genomes. We document the genomic evolutionary history and the impact of lineage-specific transposon amplification on differential genome composition. The pan-3D genome reveals evolutionary connections between transposon-driven genome size variation and both higher-order chromatin structure reorganization and the rewiring of chromatin interactome. We linked changes in chromatin structures to phenotypic differences in cotton fiber and identified regulatory variations that decode the genetic basis of fiber length, the latter enabled by sequencing 1,005 transcriptomes during fiber development. We showcase how pan-genomic, pan-3D genomic and genetic regulatory data serve as a resource for delineating the evolutionary basis of spinnable cotton fiber. Our work provides insights into the evolution of genome organization and regulation and will inform cotton improvement by enabling regulome-based approaches.

Cytochrome P450 mono-oxygenase CYP703A2 plays a central role in sporopollenin formation and ms5ms6 fertility in cotton.

The double-recessive genic male-sterile (ms) line ms5 ms6 has been used to develop cotton (Gossypium hirsutum) hybrids for many years, but its molecular-genetic basis has remained unclear. Here, we identified the Ms5 and Ms6 loci through map-based cloning and confirmed their function in male sterility through CRISPR/Cas9 gene editing. Ms5 and Ms6 are highly expressed in stages 7-9 anthers and encode the cytochrome P450 mono-oxygenases CYP703A2-A and CYP703A2-D. The ms5 mutant carries a single-nucleotide C-to-T nonsense mutation leading to premature chain termination at amino acid 312 (GhCYP703A2-A312aa ), and ms6 carries three nonsynonymous substitutions (D98E, E168K, and G198R) and a synonymous mutation (L11L). Enzyme assays showed that GhCYP703A2 proteins hydroxylate fatty acids, and the ms5 (GhCYP703A2-A312aa ) and ms6 (GhCYP703A2-DD98E,E168K,G198R ) mutant proteins have decreased enzyme activities. Biochemical and lipidomic analyses showed that in ms5 ms6 plants, C12-C18 free fatty acid and phospholipid levels are significantly elevated in stages 7-9 anthers, while stages 8-10 anthers lack sporopollenin fluorescence around the pollen, causing microspore degradation and male sterility. Overall, our characterization uncovered functions of GhCYP703A2 in sporopollenin formation and fertility, providing guidance for creating male-sterile lines to facilitate hybrid cotton production and therefore exploit heterosis for improvement of cotton.

High temperature induces male sterility via MYB66-MYB4-Casein kinase I signaling in cotton.

High temperature (HT) causes male sterility and decreases crop yields. Our previous works have demonstrated that sugar and auxin signaling pathways, Gossypium hirsutum Casein kinase I (GhCKI), and DNA methylation are all involved in HT-induced male sterility in cotton. However, the signaling mechanisms leading to distinct GhCKI expression patterns induced by HT between HT-tolerant and HT-sensitive cotton anthers remain largely unknown. Here, we identified a GhCKI promoter (ProGhCKI) region that functions in response to HT in anthers and found the transcription factor GhMYB4 binds to this region to act as an upstream positive regulator of GhCKI. In the tapetum of early-stage cotton anthers, upregulated expression of GhMYB4 under HT and overexpressed GhMYB4 under normal temperature both led to severe male sterility phenotypes, coupled with enhanced expression of GhCKI. We also found that GhMYB4 interacts with GhMYB66 to form a heterodimer to enhance its binding to ProGhCKI. However, GhMYB66 showed an expression pattern similar to GhMYB4 under HT but did not directly bind to ProGhCKI. Furthermore, HT reduced siRNA-mediated CHH DNA methylations in the GhMYB4 promoter, which enhanced the expression of GhMYB4 in tetrad stage anthers and promoted the formation of the GhMYB4/GhMYB66 heterodimer, which in turn elevated the transcription of GhCKI in the tapetum, leading to male sterility. Overall, we shed light on the GhMYB66-GhMYB4-GhCKI regulatory pathway in response to HT in cotton anthers.

Fast anther dehiscence status recognition system established by deep learning to screen heat tolerant cotton.

BACKGROUND: From an economic perspective, cotton is one of the most important crops in the world. The fertility of male reproductive organs is a key determinant of cotton yield. Anther dehiscence or indehiscence directly determines the probability of fertilization in cotton. Thus, rapid and accurate identification of cotton anther dehiscence status is important for judging anther growth status and promoting genetic breeding research. The development of computer vision technology and the advent of big data have prompted the application of deep learning techniques to agricultural phenotype research. Therefore, two deep learning models (Faster R-CNN and YOLOv5) were proposed to detect the number and dehiscence status of anthers. RESULT: The single-stage model based on YOLOv5 has higher recognition speed and the ability to deploy to the mobile end. Breeding researchers can apply this model to terminals to achieve a more intuitive understanding of cotton anther dehiscence status. Moreover, three improvement strategies are proposed for the Faster R-CNN model, where the improved model has higher detection accuracy than the YOLOv5 model. We have made three improvements to the Faster R-CNN model and after the ensemble of the three models and original Faster R-CNN model, R2 of "open" reaches to 0.8765, R2 of "close" reaches to 0.8539, R2 of "all" reaches to 0.8481, higher than the prediction results of either model alone, which are completely able to replace the manual counting results. We can use this model to quickly extract the dehiscence rate of cotton anthers under high temperature (HT) conditions. In addition, the percentage of dehiscent anthers of 30 randomly selected cotton varieties were observed from the cotton population under normal conditions and HT conditions through the ensemble of the Faster R-CNN model and manual counting. The results show that HT decreased the percentage of dehiscent anthers in different cotton lines, consistent with the manual method. CONCLUSIONS: Deep learning technology have been applied to cotton anther dehiscence status recognition instead of manual methods for the first time to quickly screen HT-tolerant cotton varieties. Deep learning can help to explore the key genetic improvement genes in the future, promoting cotton breeding and improvement.

Rapid Identification of Pollen- and Anther-Specific Genes in Response to High-Temperature Stress Based on Transcriptome Profiling Analysis in Cotton.

Anther indehiscence and pollen sterility caused by high temperature (HT) stress have become a major problem that decreases the yield of cotton. Pollen- and anther-specific genes play a critical role in the process of male reproduction and the response to HT stress. In order to identify pollen-specific genes that respond to HT stress, a comparative transcriptome profiling analysis was performed in the pollen and anthers of Gossypium hirsutum HT-sensitive Line H05 against other tissue types under normal temperature (NT) conditions, and the analysis of a differentially expressed gene was conducted in the pollen of H05 under NT and HT conditions. In total, we identified 1111 pollen-specific genes (PSGs), 1066 anther-specific genes (ASGs), and 833 pollen differentially expressed genes (DEGs). Moreover, we found that the late stage of anther included more anther- and pollen-specific genes (APSGs). Stress-related cis-regulatory elements (CREs) and hormone-responsive CREs are enriched in the promoters of APSGs, suggesting that APSGs may respond to HT stress. However, 833 pollen DEGs had only 10 common genes with 1111 PSGs, indicating that PSGs are mainly involved in the processes of pollen development and do not respond to HT stress. Promoters of these 10 common genes are enriched for stress-related CREs and MeJA-responsive CREs, suggesting that these 10 common genes are involved in the process of pollen development while responding to HT stress. This study provides a pathway for rapidly identifying cotton pollen-specific genes that respond to HT stress.

Degradation of de-esterified pctin/homogalacturonan by the polygalacturonase GhNSP is necessary for pollen exine formation and male fertility in cotton.

The pollen wall exine provides a protective layer for the male gametophyte and is largely composed of sporopollenin, which comprises fatty acid derivatives and phenolics. However, the biochemical nature of the external exine is poorly understood. Here, we show that the male sterile line 1355A of cotton mutated in NO SPINE POLLEN (GhNSP) leads to defective exine formation. The GhNSP locus was identified through map-based cloning and confirmed by genetic analysis (co-segregation test and allele prediction using the CRISPR/Cas9 system). In situ hybridization showed that GhNSP is highly expressed in tapetum. GhNSP encodes a polygalacturonase protein homologous to AtQRT3, which suggests a function for polygalacturonase in pollen exine formation. These results indicate that GhNSP is functionally different from AtQRT3, the latter has the function of microspore separation. Biochemical analysis showed that the percentage of de-esterified pectin was significantly increased in the 1355A anthers at developmental stage 8. Furthermore, immunofluorescence studies using antibodies to the de-esterified and esterified homogalacturonan (JIM5 and JIM7) showed that the Ghnsp mutant exhibits abundant of de-esterified homogalacturonan in the tapetum and exine, coupled with defective exine formation. The characterization of GhNSP provides new understanding of the role of polygalacturonase and de-esterified homogalacturonan in pollen exine formation.

A male-sterile mutant with necrosis-like dark spots on anthers was generated in cotton.

Although conventional hybrid breeding has paved the way for improving cotton production and other properties, it is undoubtedly time and labor consuming, while the cultivation of male sterile line can fix the problem. Here, we induced male sterile mutants by simultaneously editing three cotton EXCESS MICROSPOROCYTES1 (GhEMS1) genes by CRISPR/Cas9. Notably, the GhEMS1 genes are homologous to AtEMS1 genes, which inhibit the production of middle layer and tapetum cells as well, leading to male sterility in cotton. Interestingly, there are necrosis-like dark spots on the surface of the anthers of GhEMS1s mutants, which is different from AtEMS1 mutant whose anther surface is clean and smooth, suggesting that the function of EMS1 gene has not been uncovered yet. Moreover, we have detected mutations in GhEMS1 genes from T0 to T3 mutant plants, which had necrosis-like dark spots as well, indicating that the mutation of the three GhEMS1 genes could be stably inherited. Dynamic transcriptomes showed plant hormone pathway and anther development genetic network were differential expression in mutant and wild-type anthers. And the lower level of IAA content in the mutant anthers than that in the wild type at four anther developmental stages may be the reason for the male sterility. This study not only facilitates the exploration of the basic research of cotton male sterile lines, but also provides germplasms for accelerating the hybrid breeding in cotton.

Nanopore-Based Comparative Transcriptome Analysis Reveals the Potential Mechanism of High-Temperature Tolerance in Cotton (Gossypium hirsutum L.).

Extreme high temperatures are threatening cotton production around the world due to the intensification of global warming. To cope with high-temperature stress, heat-tolerant cotton cultivars have been bred, but the heat-tolerant mechanism remains unclear. This study selected heat-tolerant ('Xinluzao36') and heat-sensitive ('Che61-72') cultivars of cotton treated with high-temperature stress as plant materials and performed comparative nanopore sequencing transcriptome analysis to reveal the potential heat-tolerant mechanism of cotton. Results showed that 120,605 nonredundant sequences were generated from the raw reads, and 78,601 genes were annotated. Differentially expressed gene (DEG) analysis showed that a total of 19,600 DEGs were screened; the DEGs involved in the ribosome, heat shock proteins, auxin and ethylene signaling transduction, and photosynthesis pathways may be attributed to the heat tolerance of the heat-tolerant cotton cultivar. This study also predicted a total of 5118 long non-coding RNAs (lncRNAs)and 24,462 corresponding target genes. Analysis of the target genes revealed that the expression of some ribosomal, heat shock, auxin and ethylene signaling transduction-related and photosynthetic proteins may be regulated by lncRNAs and further participate in the heat tolerance of cotton. This study deepens our understandings of the heat tolerance of cotton.

Cotton pan-genome retrieves the lost sequences and genes during domestication and selection.

BACKGROUND: Millennia of directional human selection has reshaped the genomic architecture of cultivated cotton relative to wild counterparts, but we have limited understanding of the selective retention and fractionation of genomic components. RESULTS: We construct a comprehensive genomic variome based on 1961 cottons and identify 456 Mb and 357 Mb of sequence with domestication and improvement selection signals and 162 loci, 84 of which are novel, including 47 loci associated with 16 agronomic traits. Using pan-genome analyses, we identify 32,569 and 8851 non-reference genes lost from Gossypium hirsutum and Gossypium barbadense reference genomes respectively, of which 38.2% (39,278) and 14.2% (11,359) of genes exhibit presence/absence variation (PAV). We document the landscape of PAV selection accompanied by asymmetric gene gain and loss and identify 124 PAVs linked to favorable fiber quality and yield loci. CONCLUSIONS: This variation repertoire points to genomic divergence during cotton domestication and improvement, which informs the characterization of favorable gene alleles for improved breeding practice using a pan-genome-based approach.

An enhanced photosynthesis and carbohydrate metabolic capability contributes to heterosis of the cotton (Gossypium hirsutum) hybrid 'Huaza Mian H318', as revealed by genome-wide gene expression analysis.

BACKGROUND: Heterosis has been exploited for decades in different crops due to resulting in dramatic increases in yield, but relatively little molecular evidence on this topic was reported in cotton. RESULTS: The elite cotton hybrid variety 'Huaza Mian H318' (H318) and its parental lines were used to explore the source of its yield heterosis. A four-year investigation of yield-related traits showed that the boll number of H318 showed higher stability than that of its two parents, both in suitable and unsuitable climate years. In addition, the hybrid H318 grew faster and showed higher fresh and dry weights than its parental lines at the seedling stage. Transcriptome analysis of seedlings identified 17,308 differentially expressed genes (DEGs) between H318 and its parental lines, and 3490 extremely changed DEGs were screened out for later analysis. Most DEGs (3472/3490) were gathered between H318 and its paternal line (4-5), and only 64 DEGs were found between H318 and its maternal line (B0011), which implied that H318 displays more similar transcriptional patterns to its maternal parent at the seedling stage. GO and KEGG analyses showed that these DEGs were highly enriched in photosynthesis, lipid metabolic, carbohydrate metabolic and oxidation-reduction processes, and the expression level of these DEGs was significantly higher in H318 relative to its parental lines, which implied that photosynthesis, metabolism and stress resistances were enhanced in H318. CONCLUSION: The enhanced photosynthesis, lipid and carbohydrate metabolic capabilities contribute to the heterosis of H318 at the seedling stage, and establishes a material foundation for subsequent higher boll-setting rates in complex field environments.

A combination of genome-wide and transcriptome-wide association studies reveals genetic elements leading to male sterility during high temperature stress in cotton.

Global warming has reduced the productivity of many field-grown crops, as the effects of high temperatures can lead to male sterility in such plants. Genetic regulation of the high temperature (HT) response in the major crop cotton is poorly understood. We determined the functionality and transcriptomes of the anthers of 218 cotton accessions grown under HT stress. By analyzing transcriptome divergence and implementing a genome-wide association study (GWAS), we identified three thermal tolerance associated loci which contained 75 protein coding genes and 27 long noncoding RNAs, and provided expression quantitative trait loci (eQTLs) for 13 132 transcripts. A transcriptome-wide association study (TWAS) confirmed six causal elements for the HT response (three genes overlapped with the GWAS results) which are involved in protein kinase activity. The most susceptible gene, GhHRK1, was confirmed to be a previously uncharacterized negative regulator of the HT response in both cotton and Arabidopsis. These functional variants provide a new understanding of the genetic basis for HT tolerance in male reproductive organs.

Efficient CRISPR/Cas9 mediated Pooled-sgRNAs assembly accelerates targeting multiple genes related to male sterility in cotton.

BACKGROUND: Upland cotton (Gossypium hirsutum), harboring a complex allotetraploid genome, consists of A and D sub-genomes. Every gene has multiple copies with high sequence similarity that makes genetic, genomic and functional analyses extremely challenging. The recent accessibility of CRISPR/Cas9 tool provides the ability to modify targeted locus efficiently in various complicated plant genomes. However, current cotton transformation method targeting one gene requires a complicated, long and laborious regeneration process. Hence, optimizing strategy that targeting multiple genes is of great value in cotton functional genomics and genetic engineering. RESULTS: To target multiple genes in a single experiment, 112 plant development-related genes were knocked out via optimized CRISPR/Cas9 system. We optimized the key steps of pooled sgRNAs assembly method by which 116 sgRNAs pooled together into 4 groups (each group consisted of 29 sgRNAs). Each group of sgRNAs was compiled in one PCR reaction which subsequently went through one round of vector construction, transformation, sgRNAs identification and also one round of genetic transformation. Through the genetic transformation mediated Agrobacterium, we successfully generated more than 800 plants. For mutants identification, Next Generation Sequencing technology has been used and results showed that all generated plants were positive and all targeted genes were covered. Interestingly, among all the transgenic plants, 85% harbored a single sgRNA insertion, 9% two insertions, 3% three different sgRNAs insertions, 2.5% mutated sgRNAs. These plants with different targeted sgRNAs exhibited numerous combinations of phenotypes in plant flowering tissues. CONCLUSION: All targeted genes were successfully edited with high specificity. Our pooled sgRNAs assembly offers a simple, fast and efficient method/strategy to target multiple genes in one time and surely accelerated the study of genes function in cotton.

High day and night temperatures distinctively disrupt fatty acid and jasmonic acid metabolism, inducing male sterility in cotton.

High temperature stress is an inevitable environmental factor in certain geographical regions. To study the effect of day and night high temperature stress on male reproduction, the heat-sensitive cotton line H05 was subjected to high temperature stress. High day/normal night (HN) and normal day/high night (NH) temperature treatments were compared with normal day/normal night (NN) temperature as a control. At the anther dehiscence stage, significant differences were observed, with a reduction in flower size and filament length, and sterility in pollen, seen in NH more than in HN. A total of 36 806 differentially expressed genes were screened, which were mainly associated with fatty acid and jasmonic acid (JA) metabolic pathways. Fatty acid and JA contents were reduced more in NH than HN. Under NH, ACYL-COA OXIDASE 2 (ACO2), a JA biosynthesis gene, was down-regulated. Interestingly, aco2 CRISPR-Cas9 mutants showed male sterility under the NN condition. The exogenous application of methyl jasmonate to early-stage buds of mutants rescued the sterile pollen and indehiscent anther phenotypes at the late stage. These data show that high temperature at night may affect fatty acid and JA metabolism in anthers by suppressing GhACO2 and generate male sterility more strongly than high day temperature.

Proteomic analysis reveals that sugar and fatty acid metabolisms play a central role in sterility of the male-sterile line 1355A of cotton.

Cotton (Gossypium spp.) is one of the most important economic crops and exhibits yield-improving heterosis in specific hybrid combinations. The genic male-sterility system is the main strategy used for producing heterosis in cotton. To better understand the mechanisms of male sterility in cotton, we carried out two-dimensional electrophoresis (2-DE) and label-free quantitative proteomics analysis in the anthers of two near-isogenic lines, the male-sterile line 1355A and the male-fertile line 1355B. We identified 39 and 124 proteins that were significantly differentially expressed between these two lines in the anthers at the tetrad stage (stage 7) and uninucleate pollen stage (stage 8), respectively. Gene ontology-based analysis revealed that these differentially expressed proteins were mainly associated with pyruvate, carbohydrate, and fatty acid metabolism. Biochemical analysis revealed that in the anthers of line 1355A, glycolysis was activated, which was caused by a reduction in fructose, glucose, and other soluble sugars, and that accumulation of acetyl-CoA was increased along with a significant increase in C14:0 and C18:1 free fatty acids. However, the activities of pyruvate dehydrogenase and fatty acid biosynthesis were inhibited and fatty acid ß-oxidation was activated at the translational level in 1355A. We speculate that in the 1355A anther, high rates of glucose metabolism may promote fatty acid synthesis to enable anther growth. These results provide new insights into the molecular mechanism of genic male sterility in upland cotton.

Reference genome sequences of two cultivated allotetraploid cottons, Gossypium hirsutum and Gossypium barbadense.

Allotetraploid cotton species (Gossypium hirsutum and Gossypium barbadense) have long been cultivated worldwide for natural renewable textile fibers. The draft genome sequences of both species are available but they are highly fragmented and incomplete1-4. Here we report reference-grade genome assemblies and annotations for G. hirsutum accession Texas Marker-1 (TM-1) and G. barbadense accession 3-79 by integrating single-molecule real-time sequencing, BioNano optical mapping and high-throughput chromosome conformation capture techniques. Compared with previous assembled draft genomes1,3, these genome sequences show considerable improvements in contiguity and completeness for regions with high content of repeats such as centromeres. Comparative genomics analyses identify extensive structural variations that probably occurred after polyploidization, highlighted by large paracentric/pericentric inversions in 14 chromosomes. We constructed an introgression line population to introduce favorable chromosome segments from G. barbadense to G. hirsutum, allowing us to identify 13 quantitative trait loci associated with superior fiber quality. These resources will accelerate evolutionary and functional genomic studies in cotton and inform future breeding programs for fiber improvement.