The growth and development of plants are sensitive to their surroundings. Although numerous studies have analyzed plant transcriptomic variation, few have quantified the effect of combinations of factors or identified factor-specific effects. In this study, we performed RNA sequencing (RNA-seq) analysis on tobacco leaves derived from 10 treatment combinations of three groups of ecological factors, i.e., climate factors (CFs), soil factors (SFs), and tillage factors (TFs). We detected 4980, 2916, and 1605 differentially expressed genes (DEGs) that were affected by CFs, SFs, and TFs, which included 2703, 768, and 507 specific and 703 common DEGs (simultaneously regulated by CFs, SFs, and TFs), respectively. GO and KEGG enrichment analyses showed that genes involved in abiotic stress responses and secondary metabolic pathways were overrepresented in the common and CF-specific DEGs. In addition, we noted enrichment in CF-specific DEGs related to the circadian rhythm, SF-specific DEGs involved in mineral nutrient absorption and transport, and SF- and TF-specific DEGs associated with photosynthesis. Based on these results, we propose a model that explains how plants adapt to various ecological factors at the transcriptomic level. Additionally, the identified DEGs lay the foundation for future investigations of stress resistance, circadian rhythm and photosynthesis in tobacco.

Epimedium sagittatum (Sieb. et Zucc) Maxim is a member of the Berberidaceae family of basal eudicot plants, widely distributed and used as a traditional medicinal plant in China for therapeutic effects on many diseases with a long history. Recent data shows that E. sagittatum has a relatively large genome, with a haploid genome size of ~4496 Mbp, divided into a small number of only 12 diploid chromosomes (2n = 2x = 12). However, little is known about Epimedium genome structure and composition. Here we present the analysis of 691 kb of high-quality genomic sequence derived from 672 randomly selected plasmid clones of E. sagittatum genomic DNA, representing ~0.0154% of the genome. The sampled sequences comprised at least 78.41% repetitive DNA elements and 2.51% confirmed annotated gene sequences, with a total GC% content of 39%. Retrotransposons represented the major class of transposable element (TE) repeats identified (65.37% of all TE repeats), particularly LTR (Long Terminal Repeat) retrotransposons (52.27% of all TE repeats). Chromosome analysis and Fluorescence in situ Hybridization of Gypsy-Ty3 retrotransposons were performed to survey the E. sagittatum genome at the cytological level. Our data provide the first insights into the composition and structure of the E. sagittatum genome, and will facilitate the functional genomic analysis of this valuable medicinal plant.

This study aims to provide a theoretical basis on the regulatory mechanism of the α-l,2-fucosyltransferase (FUT1) gene in pigs by analyzing the transcriptional activity of its promoter region. On the basis of the previously obtained promoter sequence, primers upstream and downstream of the gene were designed using the restriction endonucleases KpnI and HindIII respectively, and the recombinant plasmids of the pGL3-promoter were constructed by inserting promoter sequences with partially missing regions. The resultant mutants were observed by transient transfection assay into HEK293 cells, and the transcriptional activity of the promoter region was determined by luciferase activity. The 5'-flanking region of the FUT1 gene (-1150 to +50 bp) exhibited promoter activity. The -1150-bp to -849-bp region showed negative regulation of the gene. The recombinant plasmid pGL3-898 showed the strongest luciferase activity, and the activity showed a decreasing trend when the deleted region was increased. Recombinant plasmids were successfully constructed, verified, and the positive and negative regulation areas and core promoter region were detected, providing a deeper insight into the transcriptional regulatory mechanism of the FUT1 gene.

Toll-like receptor 4 (TLR4) signal activation of macrophages can lead to endotoxin-induced uveitis (EIU). Previously, our research group has demonstrated a higher expression of TLR4 in vivo during EIU than normal. In this study, we analyzed levels of peritoneal macrophage cytokines from C3H/HeN mice with LPS stimulation in vitro to elucidate the effect of TLR4 on cytokines during EIU.

Astronauts traveling to Mars will be exposed to high levels of ionizing radiation upon leaving low-Earth orbit. During prolonged space travel, astronauts are exposed to galactic cosmic rays (GCRs) composed of protons; oxygen molecules; and high energy, high mass charged particles. Notably, oxygen molecules can travel through the shielding of spacecraft, potentially impacting 25% of the hippocampus. The aim of the current study was to assess whether 16O-particle radiation induced a behavioral deficit and histological changes in mice. Mice were sent to the National Aeronautics and Space Administration (NASA) Space Radiation Laboratory at Brookhaven National Laboratory and exposed to particulate 16O radiation at doses of 0 and 0.05 Gy. Nine months after irradiation, the mice were tested for novel object recognition and in the Y-maze, after which the animals were sacrificed. The brains were then dissected along the midsagittal plane for Golgi staining. Exposure to 0.05 Gy significantly impaired novel object recognition. However, short term memory and exploratory activity in the Y-maze were not affected. Micromorphometric analysis revealed significant decreases in mushroom spine density in the dentate gyrus and cornu Ammonis-1 and -3 of the hippocampus. Sholl analysis revealed a significant decrease in dendritic complexity in the dentate gyrus. The present data provide evidence that space radiation has deleterious effects on mature neurons associated with hippocampal learning and memory.

Porcine mammary fatty tissues represent an abundant source of natural biomaterial for generation of breast-specific extracellular matrix (ECM). Here we report the extraction of total ECM proteins from pig breast fatty tissues, the fabrication of hydrogel and porous scaffolds from the extracted ECM proteins, the structural properties of the scaffolds (tissue matrix scaffold, TMS), and the applications of the hydrogel in human mammary epithelial cell spatial cultures for cell surface receptor expression, metabolomics characterization, acini formation, proliferation, migration between different scaffolding compartments, and in vivo tumor formation. This model system provides an additional option for studying human breast diseases such as breast cancer.

Petunia, which has been prevalently cultivated in landscaping, is a dicotyledonous herbaceous flower of high ornamental value. Annually, there is a massive worldwide market demand for petunia seeds. The normal development of anther is the necessary prerequisite for the plants to generate seeds. However, the knowledge of petunia anther development processes is still limited. To better understand the mechanisms of petunia anther development, the transcriptomes and metabolomes of petunia anthers at three typical development stages were constructed and then used to detect the gene expression patterns and primary metabolite profiles during the anther development processes. Results suggested that there were many differentially-expressed genes (DEGs) that mainly participated in photosynthesis and starch and sucrose metabolism when DEGs were compared between the different development stages of anthers. In this study, fructose and glucose, which were involved in starch and sucrose metabolism, were taken as the most important metabolites by partial least-squares discriminate analysis (PLS-DA). Additionally, the qRT-PCR analysis of the photosynthetic-related genes all showed decreased expression trends along with the anther development. These pieces of evidence indicated that the activities of energy and carbohydrate metabolic pathways were gradually reduced during all the development stages of anther, which affects the sink strength. Overall, this work provides a novel and comprehensive understanding of the metabolic processes in petunia anthers.

Metformin has been reported to have body weight lowering effects while treating type 2 diabetes. However, limited studies examined the effects of metformin on adipogenesis in vitro, and available data are inconclusive and contradictory. In this study, we examined the effects of a variety of concentrations of metformin on adipocyte differentiation of 3T3-L1 preadipocytes and found metformin exhibits a dual effect on adipogenesis. Metformin at lower concentrations (1.25⁻2.5 mM) significantly induced adipogenesis while at higher concentrations (5⁻10 mM) metformin significantly inhibited adipogenesis in 3T3-L1 cells. The biphasic effect of different doses of metformin on adipogenesis was accompanied by increasing or decreasing the expression of adipogenic and lipogenic genes including peroxisome proliferator-activated receptor (PPARγ), CCAAT/enhancer binding protein α (C/EBPα), and fatty acid synthase (FASN) at both messenger RNA (mRNA) and protein levels. Furthermore, only the higher concentrations of metformin induced the phosphorylation of adenosine 5'-monophosphate (AMP)-activated protein kinase (AMPK), p38, and c-Jun N-terminal kinase (JNK) and reduced the phosphorylation of extracellular regulated protein kinases (ERK) and Akt. Pretreatment with compound C, a specific AMPK inhibitor, significantly countered high concentration of metformin-induced inhibition of adipogenesis. Taken together, these findings demonstrate that the effect of metformin on adipocyte differentiation is biphasic and dose-dependent. Lower concentrations of metformin induce adipogenesis, which could be mediated in an AMPK-independent manner, while higher concentrations of metformin inhibit adipogenesis via AMPK activation.

β-fructofuranosidase (β-FFase) belongs to the glycosyl-hydrolase family 32 (GH32), which can catalyze both the release of β-fructose from β-d-fructofuranoside substrates to hydrolyze sucrose and the synthesis of short-chain fructooligosaccharide (FOS). BmSuc1 has been cloned and identified from the silkworm Bombyx mori as a first animal type of β-FFase encoding gene. It was hypothesized that BmSUC1 plays an important role in the silkworm-mulberry adaptation system. However, there is little information about the enzymatic core sites of BmSUC1. In this study, we mutated three amino acid residues (D63, D181, and E234) that represent important conserved motifs for β-FFase activity in GH32 to alanine respectively by using site-directed mutagenesis. Recombinant proteins of three mutants and wild type BmSUC1 were obtained by using a Bac-to-Bac/BmNPV expression system and BmN cells. Enzymatic activity, kinetic properties, and substrate specificity of the four proteins were analyzed. High Performance Liquid Chromatography (HPLC) was used to compare the hydrolyzing and transfructosylating activities between D181A and wtBmSUC1. Our results revealed that the D63A and E234A mutations lost activity, suggesting that D63 and E234 are key amino acid residues for BmSUC1 to function as an enzyme. The D181A mutation significantly enhanced both hydrolyzing and transfructosylating activities of BmSUC1, indicating that D181 may not be directly involved in catalyzation. The results provide insight into the chemical catalyzation mechanism of BmSUC1 in B. mori. Up-regulated transfructosylating activity of BmSUC1 could provide new ideas for using B. mori β-FFase to produce functional FOS.

Reduced NME1 expression in melanoma cell lines, mouse models of melanoma, and melanoma specimens in human patients is associated with increased metastatic activity. Herein, we investigate the role of NME1 in repair of double-stranded breaks (DSBs) and choice of double-strand break repair (DSBR) pathways in melanoma cells. Using chromatin immunoprecipitation, NME1 was shown to be recruited rapidly and directly to DSBs generated by the homing endonuclease I-PpoI. NME1 was recruited to DSBs within 30 min, in concert with recruitment of ataxia-telangiectasia mutated (ATM) protein, an early step in DSBR complex formation, as well as loss of histone 2B. NME1 was detected up to 5 kb from the break site after DSB induction, suggesting a role in extending chromatin reorganization away from the repair site. shRNA-mediated silencing of NME1 expression led to increases in the homologous recombination (HR) and non-homologous end-joining (NHEJ) pathways of double-strand break repair (DSBR), and reduction in the low fidelity, alternative-NHEJ (A-NHEJ) pathway. These findings suggest low expression of NME1 drives DSBR towards higher fidelity pathways, conferring enhanced genomic stability necessary for rapid and error-free proliferation in invasive and metastatic cells. The novel mechanism highlighted in the current study appears likely to impact metastatic potential and therapy-resistance in advanced melanoma and other cancers.

Epithelial ovarian cancer (EOC) is one of the most lethal gynecologic malignancies. To date, the etiology of this deadly disease remains elusive. FHL2, a member of the four and a half LIM domain family, has been shown to serve either as an oncoprotein or as a tumor suppressor in various cancers. Our previous study showed that FHL2 plays a critical role in the initiation and progression of ovarian granulosa cell tumor via regulating AKT1 transcription. However, direct and systematic evidence of FHL2 in the initiation and progression of EOC remains unclear. In the present study, immunohistochemical analysis from EOC patient tissues showed that positivity and intensity of FHL2 immunosignal were up-regulated in the EOC tissues compared with normal ovary tissues. Knockdown of FHL2 in SKOV-3 cell line reduced cell growth and cell viability, blocked cell cycle progression, and inhibited cell migration. Ectopic expression of FHL2 in IGROV-1 cells which have low endogenous FHL2, promoted cell growth, improved cell viability and enhanced cell migration. Additionally, knock down of FHL2 in the SKOV-3 cell line significantly inhibited anchorage-independent growth indicated by the soft agar assay. In comparison, overexpression of FHL2 in IGROV-1 cell improved the colonies growth in soft agar. Western blot data showed that knockdown of FHL2 downregulated AKT expression level, and upregulated apoptosis related proteins such as cleaved PARP, and cleaved-lamin A. Finally, by employing stable SKOV-3/FHL2 stable knock down cell line, our data clearly showed that knockdown of FHL2 inhibited EOC xenograft initiation in vivo. Taken together, our results showed that FHL2, via regulating cell proliferation, cell cycle, and adhesion, has a critical role in regulating EOC initiation and progression. These results indicate that FHL2 could be a potential target for the therapeutic drugs against EOC.

Normal pairing and exchanging is an important basis to evaluate the genetic relationship between homologous chromosomes in a wheat background. The pairing behavior between 6V#2 and 6V#4, two chromosomes from different Dasypyrum villosum accessions, is still not clear. In this study, two wheat alien substitution lines, 6V#2 (6A) and 6V#4 (6D), were crossed to obtain the F1 hybrids and F2 segregating populations, and the testcross populations were obtained by using the F1 as a parent crossed with wheat variety Wan7107. The chromosomal behavior at meiosis in pollen mother cells (PMCs) of the F1 hybrids was observed using a genomic in situ hybridization (GISH) technique. Exchange events of two alien chromosomes were investigated in the F2 populations using nine polymerase chain reaction (PCR) markers located on the 6V short arm. The results showed that the two alien chromosomes could pair with each other to form ring- or rod-shaped bivalent chromosomes in 79.76% of the total PMCs, and most were pulled to two poles evenly at anaphase I. Investigation of the F2 populations showed that the segregation ratios of seven markers were consistent with the theoretical values 3:1 or 1:2:1, and recombinants among markers were detected. A genetic linkage map of nine PCR markers for 6VS was accordingly constructed based on the exchange frequencies and compared with the physical maps of wheat and barley based on homologous sequences of the markers, which showed that conservation of sequence order compared to 6V was 6H and 6B > 6A > 6D. In the testcross populations with 482 plants, seven showed susceptibility to powdery mildew (PM) and lacked amplification of alien chromosomal bands. Six other plants had amplification of specific bands of both the alien chromosomes at multiple sites, which suggested that the alien chromosomes had abnormal separation behavior in about 1.5% of the PMCs in F1, which resulted in some gametes containing two alien chromosomes. In addition, three new types of chromosome substitution were developed. This study lays a foundation for alien allelism tests and further assessment of the genetic relationship among 6V#2, 6V#4, and their wheat homoeologous chromosomes.

Fat deposition, which influences pork production, meat quality and growth efficiency, is an economically important trait in pigs. Numerous studies have demonstrated that stearoyl-CoA desaturase (SCD), a key enzyme that catalyzes the conversion of saturated fatty acids into monounsaturated fatty acids, is associated with fatty acid composition in pigs. As SCD was observed to be significantly induced in 3T3-L1 preadipocytes differentiation, we hypothesized that it plays a role in porcine adipocyte differentiation and fat deposition. In this study, we revealed that SCD is highly expressed in adipose tissues from seven-day-old piglets, compared to its expression in tissues from four-month-old adult pigs. Moreover, we found that SCD and lipogenesis-related genes were induced significantly in differentiated porcine adipocytes. Using CRISPR/Cas9 technology, we generated SCD-/- porcine embryonic fibroblasts (PEFs) and found that the loss of SCD led to dramatically decreased transdifferentiation efficiency, as evidenced by the decreased expression of known lipid synthesis-related genes, lower levels of oil red O staining and significantly lower levels of triglyceride content. Our study demonstrates the critical role of SCD expression in porcine adipocyte differentiation and paves the way for identifying it as the promising candidate gene for less fat deposition in pigs.

Four cyclodextrins (CD) including β-cyclodextrin (β-CD), γ-cyclodextrin (γ-CD), heptakis-O-(2-hydroxypropyl)-β-cyclodextrin (HP-β-CD), and heptakis-O-(2, 6-di-O-methyl)-β-cyclodextrin (DM-β-CD) were used as solubilizer to study the solubility enhancement of myricetin. The results of the phase solubility study showed that the presence of CDs could enhance the solubility of myricetin by forming 1:1 complexes. Among all CDs, HP-β-CD had the highest solubilization effect to myricetin. The concentration of myricetin could be 1.60 × 10-4 moL/L when the presence of HP-β-CD reached 1.00 × 10-2 moL/L, which was 31.45 times higher than myricetin's aqueous solubility. Subsequently, the HP-β-CD:myricetin complex was characterized by Fourier transform infrared spectroscopy (FT-IR), X-ray diffraction (XRD), and thermogravimetric analysis (TGA). In order to get an insight of the plausible structure of the complex, molecular docking was used to study the complexation process of HP-β-CD and myricetin. In the complex, the A ring and C ring of myricetin were complexed into the hydrophobic cavity of HP-β-CD, while the ring B was located at the wide rim of HP-β-CD. Four hydrogen bonding interactions were found between HP-β-CD and -OH groups of the guest in the HP-β-CD: myricetin complex. The complexation energy (△E) for the host-guest interactions was calculated with a negative sign, indicating the formation of the complex was an exergonic process. A 30-ns molecular dynamics simulation was conducted to the HP-β-CD: myricetin complex. Calculation results showed that no large structural deformation or position change were observed during the whole simulation time span. The average root-mean-square deviation (RMSD) changes of the host and guest were 2.444 and 1.145 Å, respectively, indicating the complex had excellent stability.

Salt stress seriously affects yield and quality of crops. The fruit of Lycium barbarum (LBF) is extensively used as functional food due to its rich nutrient components. It remains unclear how salt stress influences the quality of LBF. In this study, we identified 71 differentially accumulated metabolites (DAMs) and 1396 differentially expressed genes (DEGs) among ripe LBF with and without 300 mM of NaCl treatment. Pearson correlation analysis indicated that the metabolomic changes caused by salt stress were strongly related to oxidoreductases; hydrolases; and modifying enzymes, in particular, acyltransferases, methyltransferases and glycosyltransferases. Further analysis revealed that salt stress facilitated flavonoid glycosylation and carotenoid esterification by boosting the expression of structural genes in the biosynthetic pathways. These results suggested that salt stress prompts the modification of flavonoids and carotenoids to alleviate ROS damage, which in turn improves the quality of LBF. Our results lay a solid foundation for uncovering the underlying molecular mechanism of salt stress orchestrating LBF quality, and the candidate genes identified will be a valuable gene resource for genetic improvement of L. barbarum.

The abscisic acid (ABA) increase and auxin decline are both indicators of ripening initiation in grape berry, and norisoprenoid accumulation also starts at around the onset of ripening. However, the relationship between ABA, auxin, and norisoprenoids remains largely unknown, especially at the transcriptome level. To investigate the transcriptional and posttranscriptional regulation of the ABA and synthetic auxin 1-naphthaleneacetic acid (NAA) on norisoprenoid production, we performed time-series GC-MS and RNA-seq analyses on Vitis vinifera L. cv. Cabernet Sauvignon grape berries from pre-veraison to ripening. Higher levels of free norisoprenoids were found in ABA-treated mature berries in two consecutive seasons, and both free and total norisoprenoids were significantly increased by NAA in one season. The expression pattern of known norisoprenoid-associated genes in all samples and the up-regulation of specific alternative splicing isoforms of VviDXS and VviCRTISO in NAA-treated berries were predicted to contribute to the norisoprenoid accumulation in ABA and NAA-treated berries. Combined weighted gene co-expression network analysis (WGCNA) and DNA affinity purification sequencing (DAP-seq) analysis suggested that VviGATA26, and the previously identified switch genes of myb RADIALIS (VIT_207s0005g02730) and MAD-box (VIT_213s0158g00100) could be potential regulators of norisoprenoid accumulation. The positive effects of ABA on free norisoprenoids and NAA on total norisoprenoid accumulation were revealed in the commercially ripening berries. Since the endogenous ABA and auxin are sensitive to environmental factors, this finding provides new insights to develop viticultural practices for managing norisoprenoids in vineyards in response to changing climates.

Salt stress seriously restricts crop yield and quality, leading to an urgent need to understand its effects on plants and the mechanism of plant responses. Although phytohormones are crucial for plant responses to salt stress, the role of phytohormone signal transduction in the salt stress responses of stress-resistant species such as Sophora alopecuroides has not been reported. Herein, we combined transcriptome and metabolome analyses to evaluate expression changes of key genes and metabolites associated with plant hormone signal transduction in S. alopecuroides roots under salt stress for 0 h to 72 h. Auxin, cytokinin, brassinosteroid, and gibberellin signals were predominantly involved in regulating S. alopecuroides growth and recovery under salt stress. Ethylene and jasmonic acid signals may negatively regulate the response of S. alopecuroides to salt stress. Abscisic acid and salicylic acid are significantly upregulated under salt stress, and their signals may positively regulate the plant response to salt stress. Additionally, salicylic acid (SA) might regulate the balance between plant growth and resistance by preventing reduction in growth-promoting hormones and maintaining high levels of abscisic acid (ABA). This study provides insight into the mechanism of salt stress response in S. alopecuroides and the corresponding role of plant hormones, which is beneficial for crop resistance breeding.

The bulbil is an important vegetative reproductive organ in triploid Lilium lancifolium whose development is promoted by cytokinins. Type-B response regulators (RRs) are critical regulators that mediate primary cytokinin responses and promote cytokinin-induced gene expression. However, the function of cytokinin type-B Arabidopsis RRs (ARRs) in regulating bulbil formation is unclear. In this study, we identified five type-B LlRRs, LlRR1, LlRR2, LlRR10, LlRR11 and LlRR12, in L. lancifolium for the first time. The five LlRRs encode proteins of 715, 675, 573, 582 and 647 amino acids. All of the regulators belong to the B-I subfamily, whose members typically contain a conserved CheY-homologous receiver (REC) domain and an Myb DNA-binding (MYB) domain at the N-terminus. As transcription factors, all five type-B LlRRs localize at the nucleus and are widely expressed in plant tissues, especially during axillary meristem (AM) formation. Functional analysis showed that type-B LlRRs are involved in bulbil formation in a functionally redundant manner and can activate LlRR9 expression. In summary, our study elucidates the process by which cytokinins regulate bulbil initiation in L. lancifolium through type-B LlRRs and lays a foundation for research on the molecular mechanism of bulbil formation in the lily.

Cadmium is a carcinogen that can induce ER stress, DNA damage, oxidative stress and cell death. The yeast mitogen-activated protein kinase (MAPK) signalling pathways paly crucial roles in response to various stresses. Here, we demonstrate that the unfolded protein response (UPR) pathway, the high osmolarity glycerol (HOG) pathway and the cell wall integrity (CWI) pathway are all essential for yeast cells to defend against the cadmium-induced toxicity, including the elevated ROS and cell death levels induced by cadmium. We show that the UPR pathway is required for the cadmium-induced phosphorylation of HOG_MAPK Hog1 but not for CWI_MAPK Slt2, while Slt2 but not Hog1 is required for the activation of the UPR pathway through the transcription factors of Swi6 and Rlm1. Moreover, deletion of HAC1 and IRE1 could promote the nuclear accumulation of Hog1, and increase the cytosolic and bud neck localisation of Slt2, indicating crucial roles of Hog1 and Slt2 in regulating the cellular process in the absence of UPR pathway. Altogether, our findings highlight the significance of these two MAPK pathways of HOG and CWI and their interrelationship with the UPR pathway in responding to cadmium-induced toxicity in budding yeast.

The function and regulation of lipid metabolic genes are essential for plant male reproduction. However, expression regulation of lipid metabolic genic male sterility (GMS) genes by noncoding RNAs is largely unclear. Here, we systematically predicted the microRNA regulators of 34 maize white brown complex members in ATP-binding cassette transporter G subfamily (WBC/ABCG) genes using transcriptome analysis. Results indicate that the ZmABCG26 transcript was predicted to be targeted by zma-miR164h-5p, and their expression levels were negatively correlated in maize B73 and Oh43 genetic backgrounds based on both transcriptome data and qRT-PCR experiments. CRISPR/Cas9-induced gene mutagenesis was performed on ZmABCG26 and another lipid metabolic gene, ZmFAR1. DNA sequencing, phenotypic, and cytological observations demonstrated that both ZmABCG26 and ZmFAR1 are GMS genes in maize. Notably, ZmABCG26 proteins are localized in the endoplasmic reticulum (ER), chloroplast/plastid, and plasma membrane. Furthermore, ZmFAR1 shows catalytic activities to three CoA substrates in vitro with the activity order of C12:0-CoA > C16:0-CoA > C18:0-CoA, and its four key amino acid sites were critical to its catalytic activities. Lipidomics analysis revealed decreased cutin amounts and increased wax contents in anthers of both zmabcg26 and zmfar1 GMS mutants. A more detailed analysis exhibited differential changes in 54 monomer contents between wild type and mutants, as well as between zmabcg26 and zmfar1. These findings will promote a deeper understanding of miRNA-regulated lipid metabolic genes and the functional diversity of lipid metabolic genes, contributing to lipid biosynthesis in maize anthers. Additionally, cosegregating molecular markers for ZmABCG26 and ZmFAR1 were developed to facilitate the breeding of male sterile lines.