Metabolome profiling and transcriptome analysis filling the early crucial missing steps of piperine biosynthesis in Piper nigrum L.

Black pepper (Piper nigrum L.), the world renown as the King of Spices, is not only a flavorsome spice but also a traditional herb. Piperine, a species-specific piper amide, is responsible for the major bioactivity and pungent flavor of black pepper. However, several key steps for the biosynthesis of piperoyl-CoA (acyl-donor) and piperidine (acyl-acceptor), two direct precursors for piperine, remain unknown. In this study, we used guilt-by-association analysis of the combined metabolome and transcriptome, to identify two feruloyldiketide-CoA synthases responsible for the production of the C5 side chain scaffold feruloyldiketide-CoA intermediate, which is considered the first and important step to branch metabolic fluxes from phenylpropanoid pathway to piperine biosynthesis. In addition, we also identified the first two key enzymes for piperidine biosynthesis derived from lysine in P. nigrum, namely a lysine decarboxylase and a copper amine oxidase. These enzymes catalyze the production of cadaverine and 1-piperideine, the precursors of piperidine. In vivo and in vitro experiments verified the catalytic capability of them. In conclusion, our findings revealed enigmatic key steps of piperine biosynthetic pathway and thus provide a powerful reference for dissecting the biosynthetic logic of other piper amides.

ExamPle: explainable deep learning framework for the prediction of plant small secreted peptides.

MOTIVATION: Plant Small Secreted Peptides (SSPs) play an important role in plant growth, development, and plant-microbe interactions. Therefore, the identification of SSPs is essential for revealing the functional mechanisms. Over the last few decades, machine learning-based methods have been developed, accelerating the discovery of SSPs to some extent. However, existing methods highly depend on handcrafted feature engineering, which easily ignores the latent feature representations and impacts the predictive performance. RESULTS: Here, we propose ExamPle, a novel deep learning model using Siamese network and multi-view representation for the explainable prediction of the plant SSPs. Benchmarking comparison results show that our ExamPle performs significantly better than existing methods in the prediction of plant SSPs. Also, our model shows excellent feature extraction ability. Importantly, by utilizing in silicomutagenesis experiment, ExamPle can discover sequential characteristics and identify the contribution of each amino acid for the predictions. The key novel principle learned by our model is that the head region of the peptide and some specific sequential patterns are strongly associated with the SSPs' functions. Thus, ExamPle is expected to be a useful tool for predicting plant SSPs and designing effective plant SSPs. AVAILABILITY AND IMPLEMENTATION: Our codes and datasets are available at https://github.com/Johnsunnn/ExamPle.

Sodium Ion Pre-Intercalation of δ-MnO2 Nanosheets for High Energy Density Aqueous Zinc-Ion Batteries.

With the merits of low cost, environmental friendliness and rich resources, manganese dioxide is considered to be a promising cathode material for aqueous zinc-ion batteries (AZIBs). However, its low ion diffusion and structural instability greatly limit its practical application. Hence, we developed an ion pre-intercalation strategy based on a simple water bath method to grow in situ δ-MnO2 nanosheets on flexible carbon cloth substrate (MnO2), while pre-intercalated Na+ in the interlayer of δ-MnO2 nanosheets (Na-MnO2), which effectively enlarges the layer spacing and enhances the conductivity of Na-MnO2. The prepared Na-MnO2//Zn battery obtained a fairly high capacity of 251 mAh g-1 at a current density of 2 A g-1, a satisfactory cycle life (62.5% of its initial capacity after 500 cycles) and favorable rate capability (96 mAh g-1 at 8 A g-1). Furthermore, this study revealed that the pre-intercalation engineering of alkaline cations is an effective method to boost the properties of δ-MnO2 zinc storage and provides new insights into the construction of high energy density flexible electrodes.

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

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

Comparative proteomic analysis identified proteins and the phenylpropanoid biosynthesis pathway involved in the response to ABA treatment in cotton fiber development.

Abscisic acid (ABA) is a plant hormone that plays an important role in cotton fiber development. In this study, the physiological changes and proteomic profiles of cotton (Gossypium hirsutum) ovules were analyzed after 20 days of ABA or ABA inhibitor (ABAI) treatment. The results showed that compared to the control (CK), the fiber length was significantly decreased under ABA treatment and increased under ABAI treatment. Using a tandem mass tags-based quantitative technique, the proteomes of cotton ovules were comprehensively analyzed. A total of 7321 proteins were identified, of which 365 and 69 differentially accumulated proteins (DAPs) were identified in ABA versus CK and ABAI versus CK, respectively. Specifically, 345 and 20 DAPs were up- and down-regulated in the ABA group, and 65 and 4 DAPs were up- and down-regulated in the ABAI group, respectively. The DAPs in the ABA group were mainly enriched in the biosynthesis of secondary metabolites, phenylpropanoid biosynthesis and flavonoid secondary metabolism, whereas the DAPs in the ABAI group were mainly enriched in the indole alkaloid biosynthesis and phenylpropanoid biosynthesis pathways. Moreover, 9 proteins involved in phenylpropanoid biosynthesis were upregulated after ABA treatment, suggesting that this pathway might play important roles in the response to ABA, and 3 auxin-related proteins were upregulated, indicating that auxin might participate in the regulation of fiber development under ABAI treatment.

Identification of Key Aromatic Compounds in Basil (Ocimum L.) Using Sensory Evaluation, Metabolomics and Volatilomics Analysis.

Basil (Ocimum L.) is widely used as a flavor ingredient, however research on basil flavor is limited. In the current study, nine basil species were selected, including Ocimum basilicum L.var. pilosum (Willd.) Benth., Ocimum sanctum, Ocimum basilicum cinnamon, Ocimum gratissimum var. suave, Ocimum tashiroi, Ocimum basilicum, Ocimum americanum, Ocimum basilicum ct linalool, and Ocimum basilicum var. basilicum, and their fragrance and flavor characteristics were assessed by sensory evaluation. The results indicated that Ocimum basilicum var. basilicum and Ocimum gratissimum var. suave have a strong clove smell and exhibited a piquant taste. Metabolomics and volatilomics analyses measured 100 nonvolatile metabolites and 134 volatiles. Differential analysis showed that eugenol, γ-terpinene, germacrene D and malic acid were among the most varied metabolites in basil species. Combined with sensory evaluation results, correlation analysis revealed that ß-pinene and γ-cadinene contributed to the piquant smell, while eugenol and germacrene D contributed to the clove smell, and malic acid and L-(−)-arabitol contributed to the sweet flavor in basil. This study provided comprehensive flavor chemistry profiles of basil species and could be used as a guide for basil flavor improvement. The better understanding of objective sensory attributes and chemical composition of fresh basil could introduce the improved cultivars with preponderant traits, which is also in accordance with the various demands of breeders and growers, food producers, and consumers.

Facile hydrothermal synthesis of cobaltosic sulfide nanorods for high performance supercapacitors.

With high reactivity, electrical conductivity, theoretical specific capacitance and well redox reversibility, transition metal sulfides are considered as a promising anode material for supercapacitors. Hence, we designed a simple two-step hydrothermal process to grow Co4S3 nanorod arrays in situ on flexible carbon cloth substrates. Benefited from the larger specific surface area of nanoarrays, the binder-free Co4S3 electrode demonstrates a higher specific capacity of 1.97 F cm-2 at a current density of 2 mA cm-2, while the Co3O4 electrode has a capacity of only 0.07 F cm-2 at the same current density. Surprisingly, at a high scan rate of 200 mV s-1, the synthesized Co4S3 electrode still maintains almost 100% of its initial capacitance after 5000 cycles. Moreover, when using the prepared Co4S3 and MnO2 electrode as the anode and cathode, the fabricated flexible supercapacitor obtains a high volumetric energy density of 0.87 mW h cm-3 (power density of 0.78 W cm-3) and a peak power density of 0.89 W cm-3 (energy density of 0.50 mW h cm-3). The excellent electrochemical properties imply that there is a large market for the prepared materials in flexible energy storage devices.

Changes in primary metabolites and volatile organic compounds in cotton seedling leaves exposed to silver ions and silver nanoparticles revealed by metabolomic analysis.

In the area of climate change, nanotechnology provides handy tools for improving crop production and assuring sustainability in global agricultural system. Due to excellent physiological and biochemical properties, silver nanoparticles (AgNPs) have been widely studied for potential use in agriculture. However, there are concerns about the mechanism of the toxic effects of the accumulation of AgNPs on crop growth and development. In this study, the impacts of AgNPs on cotton (Gossypium hirsutum) seedlings were evaluated by integrating physiological and comprehensive metabolomic analyses. Potting-soil-grown, two-week-old cotton seedlings were foliar-exposed to 5 mg/plant AgNP or 0.02 mg/plant Ag+ (equivalent to the free Ag+ released from AgNPs). Primary metabolites and volatile organic compounds (VOCs) were identified by gas chromatography-mass spectrometry (GC-MS) and solid-phase microextraction (SPME) GC-MS, respectively. AgNPs inhibited the photosynthetic capacity of the cotton leaves. The metabolic spectrum analysis identified and quantified 73 primary metabolites and 45 VOCs in cotton leaves. Both treatments significantly changed the metabolite profiles of plant leaves. Among the primary metabolites, AgNPs induced marked changes in amino acids, sugars and sugar alcohols. Among the VOCs, 13 volatiles, mainly aldehydes, alkanes and terpenoids, were specifically altered only in response to AgNPs. In summary, our study showed that the comprehensive influence of AgNPs on primary metabolites and VOCs was not merely attributed to the released Ag+ but was caused by AgNP-specific effects on cotton leaves. These results provide important knowledge about the physiological and chemical changes in cotton leaves upon exposure to AgNPs and offer a new insight for supporting the sustainable use of AgNPs in agriculture.

Small RNA sequencing reveals various microRNAs involved in piperine biosynthesis in black pepper (Piper nigrum L.).

BACKGROUND: Black pepper (Piper nigrum L.), an important and long-cultivated spice crop, is native to South India and grown in the tropics. Piperine is the main pungent and bioactive alkaloid in the berries of black pepper, but the molecular mechanism for piperine biosynthesis has not been determined. MicroRNAs (miRNAs), which are classical endogenous noncoding small RNAs, play important roles in regulating secondary metabolism in many species, but less is known regarding black pepper or piperine biosynthesis. RESULTS: To dissect the functions of miRNAs in secondary metabolism especially in piperine biosynthesis, 110 known miRNAs, 18 novel miRNAs and 1007 individual targets were identified from different tissues of black pepper by small RNA sequencing. qRT-PCR and 5'-RLM-RACE experiments were conducted to validate the reliability of the sequencing data and predicted targets. We found 3 miRNAs along with their targets including miR166-4CL, miR396-PER and miR397-CCR modules that are involved in piperine biosynthesis. CONCLUSION: MiRNA regulation of secondary metabolism is a common phenomenon in plants. Our study revealed new miRNAs that regulate piperine biosynthesis, which are special alkaloids in the piper genus, and they might be useful for future piperine genetic improvement of black pepper.

Transcriptome analysis reveals MYB and WRKY transcription factors involved in banana (Musa paradisiaca AA) magnesium deficiency.

MAIN CONCLUSION: The banana development was inhibited under the long-term magnesium deficiency (MD) stress, resulting in the leaf chlorosis. MYB108 and WRKY75 are involved in regulating the growth and development of banana leaves and roots under long-term MD. Magnesium deficiency (MD) causes plant growth inhibition, ageing acceleration, yield reduction and quality decline of banana (Musa paradisiaca AA), but the molecular regulatory mechanisms underlying the changes in response to long-term MD conditions remain unknown. In this study, a long-term MD experiment was performed with banana seedlings at the four-leaf stage. Compared to those in the control group, the growth of leaves and roots of seedlings in the long-term MD treatment experimental groups was inhibited, and the Mg content and chlorophyll contents were decreased. Leaves and roots of seedlings from the control and experimental groups were subsequently collected for RNA sequencing to identify the genes that respond to long-term MD. More than 50 million reads were identified from each sample, resulting in the detection of 3500 and 948 differentially expressed genes (DEGs) in the leaves and roots, respectively. MYB and WRKY transcription factors (TFs) involved in plant stress responses were selected for further analysis, and 102 MYB and 149 WRKY TFs were differentially expressed. Furthermore, two highly differentially expressed candidate genes, MYB108 and WRKY75, were functionally analyzed using Arabidopsis mutants grown under long-term MD conditions. The results showed that the density of root hairs on the wild type (WT) was than that on the myb108 and wrky75 mutants under MD, implying that the mutants were more sensitive to MD than the WT. This research broadens our understanding the underlying molecular mechanism of banana seedlings adapted to the long-term MD condition.

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.

The Chromosome-Level Genome of Miracle Fruit (Synsepalum dulcificum) Provides New Insights Into the Evolution and Function of Miraculin.

Miracle fruit (Synsepalum dulcificum) is a rare valuable tropical plant famous for a miraculous sweetening glycoprotein, miraculin, which can modify sour flavors to sweet flavors tasted by humans. Here, we present a chromosome-level high-quality genome of S. dulcificum with an assembly genome size of ∼550 Mb, contig N50 of ∼14.14 Mb, and 37,911 annotated protein-coding genes. Phylogenetic analysis revealed that S. dulcificum was most closely related to Camellia sinensis and Diospyros oleifera, and that S. dulcificum diverged from the Diospyros genus ∼75.8 million years ago (MYA), and that C. sinensis diverged from Synsepalum ∼63.5 MYA. Ks assessment and collinearity analysis with S. dulcificum and other species suggested that a whole-genome duplication (WGD) event occurred in S. dulcificum and that there was good collinearity between S. dulcificum and Vitis vinifera. On the other hand, transcriptome and metabolism analysis with six tissues containing three developmental stages of fleshes and seeds of miracle fruit revealed that Gene Ontology (GO) terms and metabolic pathways of "cellular response to chitin," "plant-pathogen interaction," and "plant hormone signal transduction" were significantly enriched during fruit development. Interestingly, the expression of miraculin (Chr10G0299340) progressively increased from vegetative organs to reproductive organs and reached an incredible level in mature fruit flesh, with an fragments per kilobase of transcript per million (FPKM) value of ∼113,515, which was the most highly expressed gene among all detected genes. Combining the unique signal peptide and the presence of the histidine-30 residue together composed the main potential factors impacting miraculin's unique properties in S. dulcificum. Furthermore, integrated analysis of weighted gene coexpression network analysis (WGCNA), enrichment and metabolite correlation suggested that miraculin plays potential roles in regulating plant growth, seed germination and maturation, resisting pathogen infection, and environmental pressure. In summary, valuable genomic, transcriptomic, and metabolic resources provided in this study will promote the utilization of S. dulcificum and in-depth research on species in the Sapotaceae family.

Identification and functional analysis of the CorA/MGT/MRS2-type magnesium transporter in banana.

Magnesium (Mg) plays an irreplaceable role in plant growth and development. Mg transporters, especially CorA/MGT/MRS2 family proteins, played a vital role in regulating Mg content in plant cells. Although extensive work has been conducted in model crops, such as Arabidopsis, rice, and maize, the relevant information is scarce in tropical crops. In this study, 10 MaMRS2 genes in banana (Musa acuminata) were isolated from its genome and classified into five distinct clades. The putative physiochemical properties, chromosome location, gene structure, cis-acting elements, and duplication relationships in between these members were analyzed. Complementary experiments revealed that three MaMRS2 gene members (MaMRS2-1, MaMRS2-4, MaMRS2-7), from three distinct phylogenetic branches, were capable of restoring the function of Mg transport in Salmonella typhimurium mutants. Semi-quantitative RT-PCR showed that MaMRS2 genes were differentially expressed in banana cultivar 'Baxijiao' (Musa spp. AAA Cavendish) seedlings. The result was confirmed by real-time PCR analysis, in addition to tissue specific expression, expression differences among MaMRS2 members were also observed under Mg deficiency conditions. These results showed that Mg transporters may play a versatile role in banana growth and development, and our work will shed light on the functional analysis of Mg transporters in banana.

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.

Small RNA sequencing revealed various microRNAs involved in ethylene-triggered flowering process in Aechmea fasciata.

Ethylene-triggered flowering is a common phenomenon in plants of the family Bromeliaceae, but its molecular mechanism remains unclear. As a classical group of small RNAs, microRNAs play an essential role in the regulation of flowering. In this study, we found that various miRNAs participate in the ethylene-triggered flowering process in Aechmea fasciata via small RNA sequencing using juvenile and adult plants treated with ethylene for 24 hours. Finally, 63 known miRNAs, 52 novel miRNAs and 1721 target genes were identified or predicted. Expression changes of specific miRNAs were validated by qRT-PCR and northern blotting. Some predicted targets, including SPL, GAMYB and ARF, were verified in RLM-RACE experiments. Gene Ontology (GO) and KEGG analysis showed that numerous developmental and RNA-related processes were enriched. Integrated analysis of the transcriptomic data with small RNA sequencing revealed that numerous miRNAs and targets involved in ethylene-triggered flowering in A. fasciata. Our study is helpful for illuminating the molecular basis of the ethylene-triggered flowering phenomenon in Bromeliaceae.

A simple and efficient cloning system for CRISPR/Cas9-mediated genome editing in rice.

Rapidly growing genetics and bioinformatics studies provide us with an opportunity to obtain a global view of the genetic basis of traits, but also give a challenge to the function validation of candidate genes. CRISPR/Cas9 is an emerging and efficient tool for genome editing. To construct expression clones for the CRISPR/Cas9, most current methods depend on traditional cloning using Gateway reaction or specific type IIS restriction enzymes and DNA ligation, based on multiple steps of PCR. We developed a system for introducing sgRNA expression cassette(s) directly into plant binary vectors in one step. In this system, one sgRNA expression cassette(s) is generated by an optimized multiplex PCR, in which an overlapping PCR took place. Whilst, two sgRNA expression cassettes were amplified in a single round of PCR. Subsequently, an LR or Golden gate reaction was set up with unpurified PCR product and befitting destination vector. We are able to construct expression clones within 36 h, which greatly improves efficiency and saves cost. Furthermore, the efficiency of this system was verified by an agrobacterium-mediated genetic transformation in rice. The system reported here provides a much more efficient and simpler procedure to construct expression clones for CRISPR/Cas9-mediated genome editing.

Characterization of genes in guar gum biosynthesis based on quantitative RNA-sequencing in guar bean (Cyamopsis tetragonoloba).

Guar gum is an important raw material in the food, textile and oil industries, but the biosynthesis of guar gum remains unclear. To illuminate the genes involved in guar gum biosynthesis, guar beans from 30 and 40 days after flowering (DAF) were used for RNA sequencing in this study. A total of 2,535 and 2,724 preferentially expressed genes were found in 30 and 40 DAF endosperm, and 3,720 and 2,530 preferentially expressed genes were found in 30 and 40 DAF embryos, respectively. Of these, mannan synthase genes, α-galactosyltransferase genes and cellulose synthase genes were preferentially expressed in the endosperm from 30 and 40 DAF. The high expression level of these glycometabolism genes in endosperm is consistent with the expectation that the main component of guar gum is galactomannan. We believe that genes related to guar gum biosynthesis found in this study will be useful for both new variety development via genetic engineering and synthetic biology research on guar gum biosynthesis in the future.

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.