Cell-Impermeable Buffering Fluorogenic Probes for Live-Cell Super-Resolution Imaging of Plasma Membrane Morphology Dynamics.

Super-resolution fluorescence imaging has emerged as a potent tool for investigating the nanoscale structure and function of the plasma membrane (PM). Nevertheless, the challenge persists in achieving super-resolution imaging of PM dynamics due to limitations in probe photostability and issues with cell internalization staining. Herein, we report assembly-mediated buffering fluorogenic probes BMP-14 and BMP-16 exhibiting fast PM labeling and extended retention time (over 2 h) on PM. The incorporation of alkyl chains proves effective in promoting the aggregation of BMP-14 and BMP-16 into nonfluorescent nanoparticles to realize fluorogenicity and regulate the buffering capacity to rapidly replace photobleached probes ensuring stable long-term super-resolution imaging of PM. Utilizing these PM-buffering probes, we observed dynamic movements of PM filopodia and continuous shrinkage, leading to the formation of extracellular vesicles (EVs) using structured illumination microscopy (SIM). Furthermore, we discovered two distinct modes of EV fusion: one involving fusion through adjacent lipids and the other through filamentous lipid traction. The entire process of EV fusion outside the PM was dynamically tracked. Additionally, BMP-16 exhibited a unique capability of inducing single-molecule fluorescence blinking when used for cell membrane staining. This property makes BMP-16 suitable for the PAINT imaging of cell membranes.

Dynamics and timing of diversification events of ciliated eukaryotes from a large phylogenomic perspective.

Ciliophora, an exceptionally diverse lineage of unicellular eukaryotes, exhibits a remarkable range of species richness across classes in the ciliate Tree of Life. In this study, we have acquired transcriptome and genome data from 40 representative species in seven ciliate classes. Utilizing 247 genes and 105 taxa, we devised a comprehensive phylogenomic tree for Ciliophora, encompassing over 60 % of orders and constituting the most extensive dataset of ciliate species to date. We established a robust phylogenetic framework that encompasses ambiguous taxa and the major classes within the phylum. Our findings support the monophyly of each of two subphyla (Postciliodesmatophora and Intramacronucleata), along with three subclades (Protocruzia, CONTHREEP, and SAPML) nested within Intramacronucleata, and elucidate evolutionary positions among the major classes within the phylum. Drawing on the robust ciliate Tree of Life and three constraints, we estimated the radiation of Ciliophora around 1175 Ma during the middle of the Proterozoic Eon, and most of the ciliate classes diverged from their sister lineage during the latter half of this period. Additionally, based on the time-calibrated tree and species richness pattern, we investigated net diversification rates of Ciliophora and its classes. The global net diversification rate for Ciliophora was estimated at 0.004979 species/Ma. Heterogeneity in net diversification rates was evident at the class level, with faster rates observed in Oligohymenophorea and Spirotrichea than other classes within the subclades CONTHREEP and SAPML, respectively. Notably, our analysis suggests that variations in net diversification rates, rather than clade ages, appear to contribute to the differences in species richness in Ciliophora at the class level.

Genome assembly of a symbiotic balantidia (Balantidium ctenopharyngodoni) in fish hindgut.

Balantidium ctenopharyngodoni is identified as the sole ciliate species that exclusively resides within the hindgut of grass carp with high prevalence and intensity. In this study, the successful cultivation of B. ctenopharyngodoni enabled us to collect enough cells for genome sequencing. Consequently, we acquired a high-quality genome assembly spanning 68.66 Mb, encompassing a total of 22,334 nanochromosomes. Furthermore, we predicted 29,348 protein-coding genes, and 95.5% of them was supported by the RNA-seq data. The trend of GC content in the subtelomeric regions of single-gene chromosomes was similar to other ciliates containing nanochromosomes. A large number of genes encoding carbohydrate-binding modules with affinities for starch and peptidoglycans was identified. The identification of mitochondrion-related organelles (MROs) within genome indicates its well-suited adaptation to the anaerobic conditions in the hindgut environment. In summary, our results will offer resources for understanding the genetic basis and molecular adaptations of balantidia to hindgut of herbivorous fish.

Solvatochromic Buffering Fluorescent Probe Resolves the Lipid Transport and Morphological Changes during Lipid Droplet Fusion by Super-Resolution Imaging.

The varied functions of lipid droplets, which encompass the regulation of lipid and energy homeostasis, as well as their association with the occurrence of various metabolic diseases, are intricately linked to their dynamic properties. Super-resolution imaging techniques have emerged to decipher physiological processes and molecular mechanisms on the nanoscale. However, achieving long-term dynamic super-resolution imaging faces challenges due to the need for fluorescent probes with high photostability. This paper introduces LD-CF, a "buffering probe" for imaging lipid droplet dynamics using structured illumination microscopy (SIM). The polarity-sensitive LD-CF eliminates background fluorescence with a "cyan filter" strategy, enabling wash-free imaging of lipid droplets. In the fluorescent "off" state outside droplets, the probes act as a "buffering pool", replacing photobleached probes inside droplets and enabling photostable long-term SIM imaging. With this probe, three modes of lipid droplet fusion were observed, including the discovery of fusion from large to small lipid droplets. Fluorescence intensity tracking also revealed the direction of lipid transport during the lipid droplet fusion.

Multimodal Covariance Network Reflects Individual Cognitive Flexibility.

Cognitive flexibility refers to the capacity to shift between patterns of mental function and relies on functional activity supported by anatomical structures. However, how the brain's structural-functional covarying is preconfigured in the resting state to facilitate cognitive flexibility under tasks remains unrevealed. Herein, we investigated the potential relationship between individual cognitive flexibility performance during the trail-making test (TMT) and structural-functional covariation of the large-scale multimodal covariance network (MCN) using magnetic resonance imaging (MRI) and electroencephalograph (EEG) datasets of 182 healthy participants. Results show that cognitive flexibility correlated significantly with the intra-subnetwork covariation of the visual network (VN) and somatomotor network (SMN) of MCN. Meanwhile, inter-subnetwork interactions across SMN and VN/default mode network/frontoparietal network (FPN), as well as across VN and ventral attention network (VAN)/dorsal attention network (DAN) were also found to be closely related to individual cognitive flexibility. After using resting-state MCN connectivity as representative features to train a multi-layer perceptron prediction model, we achieved a reliable prediction of individual cognitive flexibility performance. Collectively, this work offers new perspectives on the structural-functional coordination of cognitive flexibility and also provides neurobiological markers to predict individual cognitive flexibility.

The P10K database: a data portal for the protist 10 000 genomes project.

Protists, a highly diverse group of microscopic eukaryotic organisms distinct from fungi, animals and plants, exert crucial roles within the earth's biosphere. However, the genomes of only a small fraction of known protist species have been published and made publicly accessible. To address this constraint, the Protist 10 000 Genomes Project (P10K) was initiated, implementing a specialized pipeline for single-cell genome/transcriptome assembly, decontamination and annotation of protists. The resultant P10K database (https://ngdc.cncb.ac.cn/p10k/) serves as a comprehensive platform, collating and disseminating genome sequences and annotations from diverse protist groups. Currently, the P10K database has incorporated 2959 genomes and transcriptomes, including 1101 newly sequenced datasets by P10K and 1858 publicly available datasets. Notably, it covers 45% of the protist orders, with a significant representation (53% coverage) of ciliates, featuring nearly a thousand genomes/transcriptomes. Intriguingly, analysis of the unique codon table usage among ciliates has revealed differences compared to the NCBI taxonomy system, suggesting a need to revise the codon tables used for these species. Collectively, the P10K database serves as a valuable repository of genetic resources for protist research and aims to expand its collection by incorporating more sequenced data and advanced analysis tools to benefit protist studies worldwide.

Feedback-related brain activity in individual decision: evidence from a gambling EEG study.

In this study, based on scalp electroencephalogram (EEG), we conducted cortical source localization and functional network analyses to investigate the underlying mechanism explaining the decision processes when individuals anticipate maximizing gambling benefits, particularly in situations where the decision outcomes are inconsistent with the profit goals. The findings shed light on the feedback monitoring process, wherein incongruity between outcomes and gambling goals triggers a more pronounced medial frontal negativity and activates the frontal lobe. Moreover, long-range theta connectivity is implicated in processing surprise and uncertainty caused by inconsistent feedback conditions, while middle-range delta coupling reflects a more intricate evaluation of feedback outcomes, which subsequently modifies individual decision-making for optimizing future rewards. Collectively, these findings deepen our comprehension of decision-making under circumstances where the profit goals are compromised by decision outcomes and provide electrophysiological evidence supporting adaptive adjustments in individual decision strategies to achieve maximum benefit.

The genome-wide meiotic recombination landscape in ciliates and its implications for crossover regulation and genome evolution.

Meiotic recombination is essential for sexual reproduction and its regulation has been extensively studied in many taxa. However, genome-wide recombination landscape has not been reported in ciliates and it remains unknown how it is affected by the unique features of ciliates: the synaptonemal complex (SC)-independent meiosis and the nuclear dimorphism. Here, we show the recombination landscape in the model ciliate Tetrahymena thermophila by analyzing single-nucleotide polymorphism datasets from 38 hybrid progeny. We detect 1021 crossover (CO) events (35.8 per meiosis), corresponding to an overall CO rate of 9.9 cM/Mb. However, gene conversion by non-crossover is rare (1.03 per meiosis) and not biased towards G or C alleles. Consistent with the reported roles of SC in CO interference, we find no obvious sign of CO interference. CO tends to occur within germ-soma common genomic regions and many of the 44 identified CO hotspots localize at the centromeric or subtelomeric regions. Gene ontology analyses show that CO hotspots are strongly associated with genes responding to environmental changes. We discuss these results with respect to how nuclear dimorphism has potentially driven the formation of the observed recombination landscape to facilitate environmental adaptation and the sharing of machinery among meiotic and somatic recombination.

First complete mitogenome of Massarineae and its contribution to phylogenetic implications in Pleosporales.

Endophytic fungi play an important role in the growth and development of traditional Chinese medicine plants. We isolated a strain of Acrocalymma vagum from the endophytic fungi of the traditional Chinese plants Paris. To accurately identify this endophytic fungal species of interest, we sequenced the mitochondrial genome of A. vagum, which is the first discovered mitochondrial genome in Massarineae. The A. vagum mitochondrial genome consists of a 35,079-bp closed circular DNA molecule containing 36 genes. Then, we compared the general sequence characteristics of A. vagum with those of Pleosporales, and the second structure of the 22 tRNAs was predicted. The phylogenetic relationship of A. vagum was constructed using two different data sets (protein-coding genes and amino acids). The phylogenetic tree shows that A. vagum is located at the root of Pleosporales. The analysis of introns shows that the number of introns increases with the increase in branch length. The results showed that monophyly was confirmed for all families in Pleosporales except for Pleosporaceae. A. vagum is an ancient species in the Pleosporales, and Pleosporaceae may require further revision. In Pleosporales, the number of introns is positively correlated with branch length, providing data for further study on the origin of introns.

FLAIR-hyperintense lesions in anti-MOG-associated encephalitis with seizures overlaying anti-N-methyl-D-aspartate receptor encephalitis: A case report.

RATIONALE: FLAIR-hyperintense lesions in anti-myelin oligodendrocyte glycoprotein (MOG)-associated encephalitis with seizures (FLAMES) is a rare clinical phenotype of anti-MOG; immunoglobulin G-associated disease is often misdiagnosed as viral encephalitis in the early stages. Anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis is an autoimmune encephalitis caused by antibodies targeting the GluN1 subunit of the NMDAR. The coexistence of anti-NMDAR encephalitis and FLAMES is very rare. PATIENT CONCERNS: A 20-year-old female patient initially presented with seizures accompanied by daytime sleepiness. DIAGNOSES: Magnetic resonance imaging revealed FLAIR-hyperintense lesions in unilateral cerebral cortex. NMDAR antibodies was positive in the cerebrospinal fluid and MOG antibodies in the serum. INTERVENTIONS: Steroid therapy was administrated. OUTCOMES: The symptoms completely relieved. At 6-month follow-up, the patient's condition remained stable. Magnetic resonance imaging showed no abnormalities in the unilateral cerebral cortex. CONCLUSION: When a patient with anti-NMDAR encephalitis or FLAMES is encountered in clinical practice, the coexistence of these diseases with double-positive anti-NMDAR and MOG antibodies should be considered and adopt appropriate evaluation and treatment.

Risk Factors for Residual Back Pain After PVP Treatment for osteoporotic Thoracolumbar Compression Fractures: A Retrospective Cohort Study.

OBJECTIVES: To explore the risk factors of residual back pain after percutaneous vertebroplasty (PVP) in patients with osteoporotic vertebral compression fracture (OVCF). METHODS: We retrospectively reviewed the records of 675 patients with OVCF treated with PVP from January 2015 to January 2020. Postoperative back pain intensity was assessed by the VAS score. Residual back pain was defined as the presence of postoperative moderate-severe pain (average VAS score≥4), and the variables included patient characteristics, baseline symptoms, imaging data and operation-related factors. Risk factors were identified with univariate and multivariate logistic regression analysis. RESULTS: Residual back pain occurred in 46 of the 675 patients included in the study, with an incidence rate of 6.8%. Multivariate logistic regression analysis showed that low Pre-BMD (OR = 3.576, P = 0.041), multiple vertebral fractures (OR = 2.795, P = 0.026), posterior fascia injury (OR = 4.083, P = 0.032), cement diffusion volume rate <0.2 (OR = 3.507, P = 0.013), facet joint violation (OR = 11.204, P < 0.001), and depression (OR = 3.562, P = 0.035) were positively correlated with residual back pain after PVP. CONCLUSIONS: Low pre-BMD (pre-bone mineral density), multiple vertebral fractures, posterior fascia injury, cement diffusion volume rate <0.2, facet joint violation and depression were the independent risk factors of residual back pain after PVP.

Disease-specific resting-state EEG network variations in schizophrenia revealed by the contrastive machine learning.

Given a multitude of genetic and environmental factors, when investigating the variability in schizophrenia (SCZ) and the first-degree relatives (R-SCZ), latent disease-specific variation is usually hidden. To reliably investigate the mechanism underlying the brain deficits from the aspect of functional networks, we newly iterated a framework of contrastive variational autoencoders (cVAEs) applied in the contrasts among three groups, to disentangle the latent resting-state network patterns specified for the SCZ and R-SCZ. We demonstrated that the comparison in reconstructed resting-state networks among SCZ, R-SCZ, and healthy controls (HC) revealed network distortions of the inner-frontal hypoconnectivity and frontal-occipital hyperconnectivity, while the original ones illustrated no differences. And only the classification by adopting the reconstructed network metrics achieved satisfying performances, as the highest accuracy of 96.80% ± 2.87%, along with the precision of 95.05% ± 4.28%, recall of 98.18% ± 3.83%, and F1-score of 96.51% ± 2.83%, was obtained. These findings consistently verified the validity of the newly proposed framework for the contrasts among the three groups and provided related resting-state network evidence for illustrating the pathological mechanism underlying the brain deficits in SCZ, as well as facilitating the diagnosis of SCZ.

Physicochemical and Functional Similarity Assessment Between Proposed Bevacizumab Biosimilar BAT1706 and Reference Bevacizumab.

BACKGROUND: BAT1706 is a proposed biosimilar of bevacizumab, a vascular endothelial growth factor A (VEGF-A)-targeting biologic used to treat several different cancers, including metastatic colorectal cancer. A comprehensive physicochemical and functional similarity assessment is a key component of demonstrating biosimilarity between a reference biologic and a proposed biosimilar. Here we report the physicochemical and functional similarity of BAT1706 and reference bevacizumab sourced from both the United States (US-bevacizumab) and the European Union (EU-bevacizumab). METHOD: A large range of product attributes, including primary and higher order structure, post-translational modifications, purity, stability, and potency, were characterized for BAT1706 and EU/US-bevacizumab using sensitive state-of-the-art analytical techniques. Up to 18 lots of US- and 29 lots of EU-bevacizumab, and 10 unique drug substance lots of BAT1706, were assessed. RESULT: BAT1706 was shown to have an identical amino acid sequence and an indistinguishable higher-order structure compared with EU/US-bevacizumab. BAT1706 and EU/US-bevacizumab also exhibited similar post-translational modifications, glycan profiles, and charge variants. Potency, assessed using a wide range of bioassays, was also shown to be comparable between BAT1706 and EU/US-bevacizumab, with statistical equivalence demonstrated for VEGF-A binding and neutralizing activity. CONCLUSION: Overall, this extensive comparability exercise demonstrated BAT1706 to match EU/US-bevacizumab in terms of all physicochemical and functional attributes assessed.

Edge-centric functional network predicts risk propensity in economic decision-making: evidence from a resting-state fMRI study.

Despite node-centric studies revealing an association between resting-state functional connectivity and individual risk propensity, the prediction of future risk decisions remains undetermined. Herein, we applied a recently emerging edge-centric method, the edge community similarity network (ECSN), to alternatively describe the community structure of resting-state brain activity and to probe its contribution to predicting risk propensity during gambling. Results demonstrated that inter-individual variability of risk decisions correlates with the inter-subnetwork couplings spanning the visual network (VN) and default mode network (DMN), cingulo-opercular task control network, and sensory/somatomotor hand network (SSHN). Particularly, participants who have higher community similarity of these subnetworks during the resting state tend to choose riskier and higher yielding bets. And in contrast to low-risk propensity participants, those who behave high-risky show stronger couplings spanning the VN and SSHN/DMN. Eventually, based on the resting-state ECSN properties, the risk rate during the gambling task is effectively predicted by the multivariable linear regression model at the individual level. These findings provide new insights into the neural substrates of the inter-individual variability in risk propensity and new neuroimaging metrics to predict individual risk decisions in advance.

Comparative Genomic and Transcriptomic Profiling Revealed the Molecular Basis of Starch Promoting the Growth and Proliferation of Balantioides coli.

Carbohydrates are the main source of nutrition for B. coli, supplying energy for cell growth and development. The research aimed at investigating the mechanism of starch on the growth and replication of B. coli. Single-cell separation was used to isolate single trophozoites of B. coli under a stereomicroscope, transcriptomic profiling was conducted based on the SMART-seq2 single-cell RNA-seq method. Comparative genomic analysis was performed on B. coli and eight other ciliates to obtain specific and expanded gene families of B. coli. GO and KEGG enrichment analysis were used to analyze the key genes of B. coli under the action of starch in the present study. The results of single-cell RNA-seq depicts starch affected the growth and replication of B. coli in two ways: (1) the cell cycle was positively promoted by the activation of the cAMP/PKA signaling pathway via glycolysis; (2) the cell autophagy was suppressed through the PI3K/AKT/mTOR pathway. Genes involved in endocytosis, carbohydrate utilization, and the cAMP/PKA signaling pathway were highly enriched in both specific and expanded gene families of B. coli. Starch can be ingested and hydrolyzed into glucose, in turn affecting various biological processes of B. coli. The molecular mechanism of the effect of starch on the growth and proliferation of B. coli by promoting cell cycle and inhibiting the autophagy of trophozoites has been elucidated in our study.

iGDP: An integrated genome decontamination pipeline for wild ciliated microeukaryotes.

Ciliates are a large group of ubiquitous and highly diverse single-celled eukaryotes that play an essential role in the functioning of microbial food webs. However, their genomic diversity is far from clear due to the need to develop cultivation methods for most species, so most research is based on wild organisms that almost invariably contain contaminants. Here we establish an integrated Genome Decontamination Pipeline (iGDP) that combines homology search, telomere reads-assisted and clustering approaches to filter contaminated ciliate genome assemblies from wild specimens. We benchmarked the performance of iGDP using genomic data from a contaminated ciliate culture and the results showed that iGDP could recall 91.9% of the target sequences with 96.9% precision. We also used a synthetic dataset to offer guidelines for the application of iGDP in the removal of various groups of contaminants. Compared with several popular metagenome binning tools, iGDP could show better performance. To further validate the effectiveness of iGDP on real-world data, we applied it to decontaminate genome assemblies of three wild ciliate specimens and obtained their genomes with high quality comparable to that of previously well-studied model ciliate genomes. It is anticipated that the newly generated genomes and the established iGDP method will be valuable community resources for detailed studies on ciliate biodiversity, phylogeny, ecology and evolution. The pipeline (https://github.com/GWang2022/iGDP) can be implemented automatically to reduce manual filtering and classification and may be further developed to apply to other microeukaryotes.

Giant proteins in a giant cell: Molecular basis of ultrafast Ca2+-dependent cell contraction.

The giant single-celled eukaryote, Spirostomum, exhibits one of the fastest movements in the biological world. This ultrafast contraction is dependent on Ca2+ rather than ATP and therefore differs to the actin-myosin system in muscle. We obtained the high-quality genome of Spirostomum minus from which we identified the key molecular components of its contractile apparatus, including two major Ca2+ binding proteins (Spasmin 1 and 2) and two giant proteins (GSBP1 and GSBP2), which act as the backbone and allow for the binding of hundreds of spasmins. The evidence suggests that the GSBP-spasmin protein complex is the functional unit of the mesh-like contractile fibrillar system, which, coupled with various other subcellular structures, provides the mechanism for repetitive ultrafast cell contraction and extension. These findings improve our understanding of the Ca2+-dependent ultrafast movement and provide a blueprint for future biomimicry, design, and construction of this kind of micromachine.

Protective effects of budesonide on LPS-induced podocyte injury by modulating macrophage M1/M2 polarization: Evidence from in vitro and in silico studies.

Budesonide (Bud), one of the most widely used lung medicines, is currently used as a repurposing medicine for immunoglobulin A nephropathy (IgAN) treatment. The progression of IgAN is related to inflammation involving macrophages and podocytes. The present study aimed to explore the effects of Bud on classically activated (M1)/alternatively activated (M2) macrophage polarization and podocyte injury under lipopolysaccharide (LPS)-induced inflammatory stress in vitro. Anti-inflammatory bioinformation of Bud was identified based on the Gene Expression Omnibus database. RAW264.7 cells were treated with normal medium, LPS, curcumin (Cur, positive control), or Bud 5, 10, or 20 µM. The expression levels of inducible nitric oxide synthase (iNOS), TNF-α, mannose receptor (CD206) and arginase (Arg)-1 were quantified by western blotting. The collected supernatants from macrophages were termed (Nor)MS, (LPS)MS, (Cur)MS and (Bud)MS. The TNF-α, IL-1ß and nitric oxide (NO) levels in the supernatants were evaluated by ELISA and Griess assay. The podocytes were cultured in different supernatants and their survival rates were assessed by bromodeoxyuridine assay. TNF signaling is an important pathway by which Bud exerts anti-inflammatory activities. Compared with the LPS group, 5, 10 and 20 µM Bud significantly increased Arg-1 and decreased iNOS expression (Six: P<0.05) and 20 µM Bud significantly increased Arg-1 and CD206 and decreased iNOS and TNF-α expression (Four: P<0.05). Cur significantly decreased iNOS and TNF-α expression (Two: P<0.05). Compared with LPS, 5, 10 and 20 µM Bud and Cur significantly decreased TNF-α, IL-1ß and NO levels (All: P<0.05). The podocyte survival rates of (Bud)MS and (Cur)MS were significantly higher than those of (LPS)MS (Four: P<0.05). The protective effect of Bud on podocyte injury is related to its modulation of M1/M2 polarization.

Comparative genomic analysis of Babesia duncani responsible for human babesiosis.

BACKGROUND: Human babesiosis, caused by parasites of the genus Babesia, is an emerging and re-emerging tick-borne disease that is mainly transmitted by tick bites and infected blood transfusion. Babesia duncani has caused majority of human babesiosis in Canada; however, limited data are available to correlate its genomic information and biological features. RESULTS: We generated a B. duncani reference genome using Oxford Nanopore Technology (ONT) and Illumina sequencing technology and uncovered its biological features and phylogenetic relationship with other Apicomplexa parasites. Phylogenetic analyses revealed that B. duncani form a clade distinct from B. microti, Babesia spp. infective to bovine and ovine species, and Theileria spp. infective to bovines. We identified the largest species-specific gene family that could be applied as diagnostic markers for this pathogen. In addition, two gene families show signals of significant expansion and several genes that present signatures of positive selection in B. duncani, suggesting their possible roles in the capability of this parasite to infect humans or tick vectors. CONCLUSIONS: Using ONT sequencing and Illumina sequencing technologies, we provide the first B. duncani reference genome and confirm that B. duncani forms a phylogenetically distinct clade from other Piroplasm parasites. Comparative genomic analyses show that two gene families are significantly expanded in B. duncani and may play important roles in host cell invasion and virulence of B. duncani. Our study provides basic information for further exploring B. duncani features, such as host-parasite and tick-parasite interactions.

Identification of human mental workload levels in a language comprehension task with imbalance neurophysiological data.

BACKGROUND AND OBJECTIVE: Operator's capability for accurately comprehending verbal commands is critically important to maintain the performance of human-machine interaction. It can be evaluated by human mental workload measured with electroencephalography (EEG). However, the time duration of different workload conditions within a task session is unequal due to varied psychophysiological processes across individuals. It leads to data imbalance of the EEG for training workload classifiers. METHODS: In this study, we propose an EEG feature oversampling technique, Gaussian-SMOTE based feature ensemble (GSMOTE-FE), for workload recognition with imbalanced classes. First, artificial EEG instances are drawn from a Gaussian distribution in the margin between the minority and majority workload classes. Tomek links are detected as clues to remove redundant feature vectors. Then, we embed a feature selection module based on the GINI importance while an ensemble classifier committee with bootstrap aggregating is used to further enhance classification performance. RESULTS: We validate the GSMOTE-FE framework based on an experiment that simulates operators to understand the correct meaning of the instructions in the Chinese language. Participants' EEG signals and reaction time data were both recorded to validate the proposed workload classifier. Workload classification accuracy and Macro-F1 values are 0.6553 and 0.5862, respectively. Corresponding G-mean and AUC achieve at 0.5757 and 0.5958, respectively. CONCLUSIONS: The performance of the GSMOTE-FE is demonstrated to be comparable with the advanced oversampling techniques. The workload classifier has the capability to indicate low and high levels of the task demand of the Chinese language understanding task.