Downregulation of apoptotic repressor AVEN exacerbates cardiac injury after myocardial infarction.

Myocardial infarction (MI) is a leading cause of heart failure (HF), associated with morbidity and mortality worldwide. As an essential part of gene expression regulation, the role of alternative polyadenylation (APA) in post-MI HF remains elusive. Here, we revealed a global, APA-mediated, 3' untranslated region (3' UTR)-lengthening pattern in both human and murine post-MI HF samples. Furthermore, the 3' UTR of apoptotic repressor gene, AVEN, is lengthened after MI, contributing to its downregulation. AVEN knockdown increased cardiomyocyte apoptosis, whereas restoration of AVEN expression substantially improved cardiac function. Mechanistically, AVEN 3' UTR lengthening provides additional binding sites for miR-30b-5p and miR-30c-5p, thus reducing AVEN expression. Additionally, PABPN1 (poly(A)-binding protein 1) was identified as a potential regulator of AVEN 3' UTR lengthening after MI. Altogether, our findings revealed APA as a unique mechanism regulating cardiac injury in response to MI and also indicated that the APA-regulated gene, AVEN, holds great potential as a critical therapeutic target for treating post-MI HF.

Approximate estimation of cell-type resolution transcriptome in bulk tissue through matrix completion.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for uncovering cellular heterogeneity. However, the high costs associated with this technique have rendered it impractical for studying large patient cohorts. We introduce ENIGMA (Deconvolution based on Regularized Matrix Completion), a method that addresses this limitation through accurately deconvoluting bulk tissue RNA-seq data into a readout with cell-type resolution by leveraging information from scRNA-seq data. By employing a matrix completion strategy, ENIGMA minimizes the distance between the mixture transcriptome obtained with bulk sequencing and a weighted combination of cell-type-specific expression. This allows the quantification of cell-type proportions and reconstruction of cell-type-specific transcriptomes. To validate its performance, ENIGMA was tested on both simulated and real datasets, including disease-related tissues, demonstrating its ability in uncovering novel biological insights.

MicroRNAs 24 and 27 Suppress Allergic Inflammation and Target a Network of Regulators of T Helper 2 Cell-Associated Cytokine Production.

MicroRNAs (miRNAs) are important regulators of cell fate decisions in immune responses. They act by coordinate repression of multiple target genes, a property that we exploited to uncover regulatory networks that govern T helper-2 (Th2) cells. A functional screen of individual miRNAs in primary T cells uncovered multiple miRNAs that inhibited Th2 cell differentiation. Among these were miR-24 and miR-27, miRNAs coexpressed from two genomic clusters, which each functioned independently to limit interleukin-4 (IL-4) production. Mice lacking both clusters in T cells displayed increased Th2 cell responses and tissue pathology in a mouse model of asthma. Gene expression and pathway analyses placed miR-27 upstream of genes known to regulate Th2 cells. They also identified targets not previously associated with Th2 cell biology which regulated IL-4 production in unbiased functional testing. Thus, elucidating the biological function and target repertoire of miR-24 and miR-27 reveals regulators of Th2 cell biology.

ipaQTL-atlas: an atlas of intronic polyadenylation quantitative trait loci across human tissues.

Functional interpretation of disease-associated non-coding variants remains a significant challenge in the post-GWAS era. Our recent study has identified 3'UTR alternative polyadenylation (APA) quantitative trait loci (3'aQTLs) and connects APA events with QTLs as a major driver of human traits and diseases. Besides 3'UTR, APA events can also occur in intron regions, and increasing evidence has connected intronic polyadenylation with disease risk. However, systematic investigation of the roles of intronic polyadenylation in human diseases remained challenging due to the lack of a comprehensive database across a variety of human tissues. Here, we developed ipaQTL-atlas (http://bioinfo.szbl.ac.cn/ipaQTL) as the first comprehensive portal for intronic polyadenylation. The ipaQTL-atlas is based on the analysis of 15 170 RNA-seq data from 838 individuals across 49 Genotype-Tissue Expression (GTEx v8) tissues and contains ∼0.98 million SNPs associated with intronic APA events. It provides an interface for ipaQTLs search, genome browser, boxplots, and data download, as well as the visualization of GWAS and ipaQTL colocalization results. ipaQTL-atlas provides a one-stop portal to access intronic polyadenylation information and could significantly advance the discovery of APA-associated disease susceptibility genes.

Comprehensive mapping of alternative polyadenylation site usage and its dynamics at single-cell resolution.

Alternative polyadenylation (APA) plays an important role in posttranscriptional gene regulation such as transcript stability and translation efficiency. However, our knowledge about APA dynamics at the single-cell level is largely unexplored. Here, we developed single-cell polyadenylation sequencing, a strand-specific approach for sequencing the 3' end of transcripts, to investigate the landscape of APA at the single-cell level. By analyzing several cell lines, we found many genes using multiple polyA sites in bulk data are prone to use only one polyA site in each single cell. Interestingly, cell cycle genes were significantly enriched in genes with high variation in polyA site usages. Furthermore, the 414 genes showing a polyA site usage switch after cell synchronization enriched cell cycle genes, while the differentially expressed genes after cell synchronization did not enrich cell cycle genes. We further identified 812 genes showing polyA site usage changes between neighboring cell cycles, which were grouped into six clusters, with cell phase-specific functional categories enriched in each cluster. Deletion of one polyA site in MSL1 and SCCPDH results in slower and faster cell cycle progression, respectively, supporting polyA site usage switch played an important role in cell cycle. These results indicate that APA is an important layer for cell cycle regulation.

Systematic evaluation of the effect of polyadenylation signal variants on the expression of disease-associated genes.

Single nucleotide variants (SNVs) within polyadenylation signals (PASs), a specific six-nucleotide sequence required for mRNA maturation, can impair RNA-level gene expression and cause human diseases. However, there is a lack of genome-wide investigation and systematic confirmation tools for identifying PAS variants. Here, we present a computational strategy to integrate the most reliable resources for discovering distinct genomic features of PAS variants and also develop a credible and convenient experimental tool to validate the effect of PAS variants on expression of disease-associated genes. This approach will greatly accelerate the deciphering of PAS variation-related human diseases.

Cancer-associated dynamics and potential regulators of intronic polyadenylation revealed by IPAFinder using standard RNA-seq data.

Intronic polyadenylation (IpA) usually leads to changes in the coding region of an mRNA, and its implication in diseases has been recognized, although at its very beginning status. Conveniently and accurately identifying IpA is of great importance for further evaluating its biological significance. Here, we developed IPAFinder, a bioinformatic method for the de novo identification of intronic poly(A) sites and their dynamic changes from standard RNA-seq data. Applying IPAFinder to 256 pan-cancer tumor/normal pairs across six tumor types, we discovered 490 recurrent dynamically changed IpA events, some of which are novel and derived from cancer-associated genes such as TSC1, SPERD2, and CCND2 Furthermore, IPAFinder revealed that IpA could be regulated by factors related to splicing and m6A modification. In summary, IPAFinder enables the global discovery and characterization of biologically regulated IpA with standard RNA-seq data and should reveal the biological significance of IpA in various processes.

Characterization and optimization of 5´ untranslated region containing poly-adenine tracts in Kluyveromyces marxianus using machine-learning model.

BACKGROUND: The 5´ untranslated region (5´ UTR) plays a key role in regulating translation efficiency and mRNA stability, making it a favored target in genetic engineering and synthetic biology. A common feature found in the 5´ UTR is the poly-adenine (poly(A)) tract. However, the effect of 5´ UTR poly(A) on protein production remains controversial. Machine-learning models are powerful tools for explaining the complex contributions of features, but models incorporating features of 5´ UTR poly(A) are currently lacking. Thus, our goal is to construct such a model, using natural 5´ UTRs from Kluyveromyces marxianus, a promising cell factory for producing heterologous proteins. RESULTS: We constructed a mini-library consisting of 207 5´ UTRs harboring poly(A) and 34 5´ UTRs without poly(A) from K. marxianus. The effects of each 5´ UTR on the production of a GFP reporter were evaluated individually in vivo, and the resulting protein abundance spanned an approximately 450-fold range throughout. The data were used to train a multi-layer perceptron neural network (MLP-NN) model that incorporated the length and position of poly(A) as features. The model exhibited good performance in predicting protein abundance (average R2 = 0.7290). The model suggests that the length of poly(A) is negatively correlated with protein production, whereas poly(A) located between 10 and 30 nt upstream of the start codon (AUG) exhibits a weak positive effect on protein abundance. Using the model as guidance, the deletion or reduction of poly(A) upstream of 30 nt preceding AUG tended to improve the production of GFP and a feruloyl esterase. Deletions of poly(A) showed inconsistent effects on mRNA levels, suggesting that poly(A) represses protein production either with or without reducing mRNA levels. CONCLUSION: The effects of poly(A) on protein production depend on its length and position. Integrating poly(A) features into machine-learning models improves simulation accuracy. Deleting or reducing poly(A) upstream of 30 nt preceding AUG tends to enhance protein production. This optimization strategy can be applied to enhance the yield of K. marxianus and other microbial cell factories.

Preliminary study on AI-assisted diagnosis of bone remodeling in chronic maxillary sinusitis.

OBJECTIVE: To construct the deep learning convolution neural network (CNN) model and machine learning support vector machine (SVM) model of bone remodeling of chronic maxillary sinusitis (CMS) based on CT image data to improve the accuracy of image diagnosis. METHODS: Maxillary sinus CT data of 1000 samples in 500 patients from January 2018 to December 2021 in our hospital was collected. The first part is the establishment and testing of chronic maxillary sinusitis detection model by 461 images. The second part is the establishment and testing of the detection model of chronic maxillary sinusitis with bone remodeling by 802 images. The sensitivity, specificity and accuracy and area under the curve (AUC) value of the test set were recorded, respectively. RESULTS: Preliminary application results of CT based AI in the diagnosis of chronic maxillary sinusitis and bone remodeling. The sensitivity, specificity and accuracy of the test set of 93 samples of CMS, were 0.9796, 0.8636 and 0.9247, respectively. Simultaneously, the value of AUC was 0.94. And the sensitivity, specificity and accuracy of the test set of 161 samples of CMS with bone remodeling were 0.7353, 0.9685 and 0.9193, respectively. Simultaneously, the value of AUC was 0.89. CONCLUSION: It is feasible to use artificial intelligence research methods such as deep learning and machine learning to automatically identify CMS and bone remodeling in MSCT images of paranasal sinuses, which is helpful to standardize imaging diagnosis and meet the needs of clinical application.

Comprehensive analysis reveals distinct mutational signature and its mechanistic insights of alcohol consumption in human cancers.

Alcohol consumption is a critical risk factor for multiple types of cancer. A genome can be attacked and acquire numerous somatic mutations in the environment of alcohol exposure. Mutational signature has the capacity illustrating the complex somatic mutation patterns in cancer genome. Recent studies have discovered distinct mutational signatures associating with alcohol consumption in liver and esophageal cancers. However, their prevalence among diverse cancers, impact of genetic background and origin of alcohol-induced mutational signatures remain unclear. By a comprehensive bioinformatics analysis on somatic mutations from patients of four cancer types with drinking information, we identified nine mutational signatures (signatures B-J), among which signature J (similar to COSMIC signature 16) was distinctive to alcohol drinking. Signature J was associated with HNSC, ESCA and LIHC but not PAAD. Interestingly, patients with mutated allele rs1229984 in ADH1B had lower level of signature J while mutated allele rs671 in ALDH2 exhibited higher signature J abundance, suggesting acetaldehyde is one cause of signature J. Intriguingly, somatic mutations of three potential cancer driver genes (TP53, CUL3 and NSD1) were found the critical contributors for increased mutational load of signature J in alcohol consumption patients. Furthermore, signature J was enriched with early accumulated clonal mutations compared to mutations derived from late tumor growth. This study systematically characterized alcohol-related mutational signature and indicated mechanistic insights into the prevalence, origin and gene-environment interaction regarding the risk oncogenic mutations associated with alcohol intake.

Strategies adopted by Salmonella to survive in host: a review.

Salmonella, a Gram-negative bacterium that infects humans and animals, causes diseases ranging from gastroenteritis to severe systemic infections. Here, we discuss various strategies used by Salmonella against host cell defenses. Epithelial cell invasion largely depends on a Salmonella pathogenicity island (SPI)-1-encoded type 3 secretion system, a molecular syringe for injecting effector proteins directly into host cells. The internalization of Salmonella into macrophages is primarily driven by phagocytosis. After entering the host cell cytoplasm, Salmonella releases many effectors to achieve intracellular survival and replication using several secretion systems, primarily an SPI-2-encoded type 3 secretion system. Salmonella-containing vacuoles protect Salmonella from contacting bactericidal substances in epithelial cells and macrophages. Salmonella modulates the immunity, metabolism, cell cycle, and viability of host cells to expand its survival in the host, and the intracellular environment of Salmonella-infected cells promotes its virulence. This review provides insights into how Salmonella subverts host cell defenses for survival.

Comprehensive characterization of somatic variants associated with intronic polyadenylation in human cancers.

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.

Independent component analysis based gene co-expression network inference (ICAnet) to decipher functional modules for better single-cell clustering and batch integration.

With the tremendous increase of publicly available single-cell RNA-sequencing (scRNA-seq) datasets, bioinformatics methods based on gene co-expression network are becoming efficient tools for analyzing scRNA-seq data, improving cell type prediction accuracy and in turn facilitating biological discovery. However, the current methods are mainly based on overall co-expression correlation and overlook co-expression that exists in only a subset of cells, thus fail to discover certain rare cell types and sensitive to batch effect. Here, we developed independent component analysis-based gene co-expression network inference (ICAnet) that decomposed scRNA-seq data into a series of independent gene expression components and inferred co-expression modules, which improved cell clustering and rare cell-type discovery. ICAnet showed efficient performance for cell clustering and batch integration using scRNA-seq datasets spanning multiple cells/tissues/donors/library types. It works stably on datasets produced by different library construction strategies and with different sequencing depths and cell numbers. We demonstrated the capability of ICAnet to discover rare cell types in multiple independent scRNA-seq datasets from different sources. Importantly, the identified modules activated in acute myeloid leukemia scRNA-seq datasets have the potential to serve as new diagnostic markers. Thus, ICAnet is a competitive tool for cell clustering and biological interpretations of single-cell RNA-seq data analysis.

The Role of Tumor Metabolic Reprogramming in Tumor Immunity.

The occurrence and development of tumors require the metabolic reprogramming of cancer cells, namely the alteration of flux in an autonomous manner via various metabolic pathways to meet increased bioenergetic and biosynthetic demands. Tumor cells consume large quantities of nutrients and produce related metabolites via their metabolism; this leads to the remodeling of the tumor microenvironment (TME) to better support tumor growth. During TME remodeling, the immune cell metabolism and antitumor immune activity are affected. This further leads to the escape of tumor cells from immune surveillance and therefore to abnormal proliferation. This review summarizes the regulatory functions associated with the abnormal biosynthesis and activity of metabolic signaling molecules during the process of tumor metabolic reprogramming. In addition, we provide a comprehensive description of the competition between immune cells and tumor cells for nutrients in the TME, as well as the metabolites required for tumor metabolism, the metabolic signaling pathways involved, and the functionality of the immune cells. Finally, we summarize current research targeted at the development of tumor immunotherapy. We aim to provide new concepts for future investigations of the mechanisms underlying the metabolic reprogramming of tumors and explore the association of these mechanisms with tumor immunity.

[Contamination status and exposure risk of mycotoxins in Coicis Semen].

By investigating the contamination status and predicting the exposure risk of mycotoxin in Coicis Semen, we aim to provide guidance for the safety supervision of Chinese medicinal materials and the formulation(revision) of mycotoxin limit standards. The content of 14 mycotoxins in the 100 Coicis Semen samples collected from five major markets of Chinese medicinal materials in China was determined by UPLC-MS/MS. The probability evaluation model based on Monte Carlo simulation method was established after Chi-square test and One-way ANOVA of the sample contamination data. Health risk assessment was performed on the basis of margin of exposure(MOE) and margin of safety(MOS). The results showed that zearalenone(ZEN), aflatoxin B_1(AFB_1), deoxynivalenol(DON), sterigmatocystin(ST), and aflatoxin B_2(AFB_2) in the Coicis Semen samples had the detection rates of 84%, 75%, 36%, 19%, and 18%, and the mean contamination levels of 117.42, 4.78, 61.16, 6.61, and 2.13 µg·kg~(-1), respectively. According to the limit standards in the Chinese Pharmacopoeia(2020 edition), AFB_1, AFs and ZEN exceeded the standards to certain extents, with the over-standard rates of 12.0%, 9.0%, and 6.0%, respectively. The exposure risks of Coicis Semen to AFB_1, AFB2, ST, DON, and ZEN were low, while 86% of the samples were contaminated with two or more toxins, which needs more attention. It is suggested that the research on the combined toxicity of different mycotoxins should be strengthened to accelerate the cumulative exposure assessment of mixed contaminations and the formulation(revision) of toxin limit standards.

H3K4 Methyltransferase Smyd3 Mediates Vascular Smooth Muscle Cell Proliferation, Migration, and Neointima Formation.

[Figure: see text].

Follicular CXCR5- expressing CD8(+) T cells curtail chronic viral infection.

During chronic viral infection, virus-specific CD8(+) T cells become exhausted, exhibit poor effector function and lose memory potential. However, exhausted CD8(+) T cells can still contain viral replication in chronic infections, although the mechanism of this containment is largely unknown. Here we show that a subset of exhausted CD8(+) T cells expressing the chemokine receptor CXCR5 has a critical role in the control of viral replication in mice that were chronically infected with lymphocytic choriomeningitis virus (LCMV). These CXCR5(+) CD8(+) T cells were able to migrate into B-cell follicles, expressed lower levels of inhibitory receptors and exhibited more potent cytotoxicity than the CXCR5(-) [corrected] subset. Furthermore, we identified the Id2-E2A signalling axis as an important regulator of the generation of this subset. In patients with HIV, we also identified a virus-specific CXCR5(+) CD8(+) T-cell subset, and its number was inversely correlated with viral load. The CXCR5(+) subset showed greater therapeutic potential than the CXCR5(-) [corrected] subset when adoptively transferred to chronically infected mice, and exhibited synergistic reduction of viral load when combined with anti-PD-L1 treatment. This study defines a unique subset of exhausted CD8(+) T cells that has a pivotal role in the control of viral replication during chronic viral infection.

Median nerve block increases the success rate of radial artery cannulation in women with gestational hypertension undergoing cesarean section.

BACKGROUND: The radial artery cannulation helps to maintain the stability of maternal hemodynamics and reduce complications, however, it is difficult for women with gestational hypertension. Ultrasound-guided median nerve block can cause arterial vasodilation, which may improve the success rate of radial artery cannulation. METHODS: Ninety-two women with gestational hypertension and risks of intra-operative bleeding undergoing cesarean section following failed ultrasound-guided cannulation were identified and randomized into the median nerve block group and control group. Median nerve block was performed under the guidance of ultrasound in the middle forearm and 5 ml of 0.5% lidocaine was injected. Subcutaneous local block was administered in the control group. The ultrasound-guided radial artery cannulation was performed ten minutes after blocking. Baseline measurements (T1) were performed after 10 minutes of rest. All variables were measured again at 10 (T2) and 30 (T3) minutes after median nerve block or local block. The primary outcome was the success rate of radial artery cannulation within 10 minutes after blocking. The puncture time, number of attempts, the overall complications, and ultrasonographic measurements including radial artery diameter and cross-sectional area were recorded before (T1), 10 minutes (T2) after, and 30 minutes (T3) after block. RESULTS: A total of 92 pregnant women were identified and completed the follow-up. As compared to control group, the first-attempt success rate of radial artery cannulation was significantly higher (95.7% vs78.3%, p = 0.027) and procedure time to success was significantly shorter (118 ± 19 s vs 172 ± 66 s, p < 0.001) in median nerve group. Median nerve group also had a significantly less overall number of attempts (p = 0.024). Compared with control group, the diameter and cross-sectional area of radial artery increased significantly at the T2 and T3 points in median nerve group (p < 0.001), as well as percentage change of radial artery diameter and CSA. No difference was observed in the overall complication at chosen radial artery, which including vasospasm (21.7% vs 28.3%; p = 0.470) and hematoma (4.3% vs 8.7%; p = 0.677). CONCLUSIONS: Ultrasound-guided median nerve block can increase the first-attempt success rate of chosen radial artery cannulation in women with gestational hypertension and risks of intra-operative bleeding undergoing cesarean section following failed radial artery cannulation, and especially for those anesthesiologists with less experienced in radial artery cannulation. TRIAL REGISTRATION: ChiCTR2100052862; http://www.chictr.org.cn , Principal investigator: MEN, Date of registration: 06/11/2021.

Measurement of HONO flux using the aerodynamic gradient method over an agricultural field in the Huaihe River Basin, China.

To investigate nitrous acid (HONO) levels and potential HONO sources above crop rotation fields. The HONO fluxes were measured by the aerodynamic gradient (AG) method from 14 December 2019 to 2 January 2020 over an agricultural field in the Huaihe River Basin. The ambient HONO levels were measured at two different heights (0.15 and 1.5 m), showing a typical diurnal cycle with low daytime levels and high nighttime levels. The upward HONO fluxes were mostly observed during the day, whereas deposition dominated at night. The diurnal variation of HONO flux followed solar radiation, with a noontime maximum of 0.2 nmol/(m2∙sec). The average upward HONO flux of 0.06 ± 0.17 nmol/(m2∙sec) indicated that the agricultural field was a net source for atmospheric HONO. The higher HONO/NO2 ratio and NO2-to-HONO conversion rate close to the surface suggested that nocturnal HONO was formed and released near the ground. The unknown HONO source was derived from the daytime HONO budget analysis, with an average strength of 0.31 ppbV/hr at noontime. The surface HONO flux, which was highly correlated with the photolysis frequency J(NO2) (R2 = 0.925) and the product of J(NO2) × NO2 (R2 = 0.840), accounted for ∼23% of unknown daytime HONO source. The significant correlation between HONO fluxes and J(NO2) suggests a light-driven HONO formation mechanism responsible for the surface HONO flux during daytime.

[Methods for synchronous detection of 14 mycotoxins in Pseudostellariae Radix and investigation on its contamination].

This study aimed to establish a method for synchronous detection of 14 mycotoxins in Pseudostellariae Radix and investigate its contamination with mycotoxins, so as to provide technical guidance for monitoring the quality of Chinese medicinal materials and medication safety. The sample was extracted with 80% acetonitrile in an oscillator for 1 h, purified using the modified QuEChERS purifying agent(0.1 g PSA + 0.3 g C_(18) + 0.3 g MgSO_4), and separated on a Waters HSS T3 chromatographic column(2.1 mm×100 mm, 1.8 µm). The gradient elution was carried out with 0.1% formic acid in water and acetonitrile, followed by the scanning in the multi-reaction monitoring(MRM) mode and the analysis of mycotoxin contamination in 26 Pseudostellariae Radix samples. The recovery rates of the established method were within the range of 82.17%-113.6%, with the RSD values less than 7% and the limits of quantification(LOQ) being 0.019-0.976 µg·kg~(-1). The detection rate of 14 mycotoxins in 26 batches of medicinal materials was 53.85%. The detection rate of sterigmatocystin(ST) was the highest, followed by those of zearalenone(ZEN), aflatoxin G_2(AFG_2), fumonisin B_1(FB_1), HT-2 toxin, and nivalenol(NIV). Their respective detection rates were 38.46%, 26.92%, 23.08%, 11.54%, 11.54%, and 7.69%, with the pollution ranges being 1.48-69.65, 0.11-31.05, 0.11-0.66, 0.28-0.83, 20.86-42.56, and 0.46-1.84 µg·kg~(-1), respectively. The established method for the detection of 14 mycotoxins is accurate, fast and reliable. The research results have very important practical significance for guiding the monitoring and prevention and control of exogenous fungal contamination of Chinese medicinal materials.