Novel Strategy for Human Deep Vein Thrombosis Diagnosis Based on Metabolomics and Stacking Machine Learning.

Deep vein thrombosis (DVT) is a serious health issue that often leads to considerable morbidity and mortality. Diagnosis of DVT in a clinical setting, however, presents considerable challenges. The fusion of metabolomics techniques and machine learning methods has led to high diagnostic and prognostic accuracy for various pathological conditions. This study explored the synergistic potential of dual-platform metabolomics (specifically, gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS)) to expand the detection of metabolites and improve the precision of DVT diagnosis. Sixty-one differential metabolites were identified in serum from DVT patients: 22 from GC-MS and 39 from LC-MS. Among these, five key metabolites were highlighted by SHapley Additive exPlanations (SHAP)-guided feature engineering and then used to develop a stacking diagnostic model. Additionally, a user-friendly interface application system was developed to streamline and automate the application of the diagnostic model, enhancing its practicality and accessibility for clinical use. This work showed that the integration of dual-platform metabolomics with a stacking machine learning model enables faster and more accurate diagnosis of DVT in clinical environments.

Novel ratio-expressions of genes enables estimation of wound age in contused skeletal muscle.

Given that combination with multiple biomarkers may well raise the predictive value of wound age, it appears critically essential to identify new features under the limited cost. For this purpose, the present study explored whether the gene expression ratios provide unique time information as an additional indicator for wound age estimation not requiring the detection of new biomarkers and allowing full use of the available data. The expression levels of four wound-healing genes (Arid5a, Ier3, Stom, and Lcp1) were detected by real-time polymerase chain reaction, and a total of six expression ratios were calculated among these four genes. The results showed that the expression levels of four genes and six ratios of expression changed time-dependent during wound repair. The six expression ratios provided additional temporal information, distinct from the four genes analyzed separately by principal component analysis. The overall performance metrics for cross-validation and external validation of four typical prediction models were improved when six ratios of expression were added as additional input variables. Overall, expression ratios among genes provide temporal information and have excellent potential as predictive markers for wound age estimation. Combining the expression levels of genes with ratio-expression of genes may allow for more accurate estimates of the time of injury.

Development of a screening system of gene sets for estimating the time of early skeletal muscle injury based on second-generation sequencing technology.

The present study is aimed to address the challenge of wound age estimation in forensic science by identifying reliable genetic markers using low-cost and high-precision second-generation sequencing technology. A total of 54 Sprague-Dawley rats were randomly assigned to a control group or injury groups, with injury groups being further divided into time points (4 h, 8 h, 12 h, 16 h, 20 h, 24 h, 28 h, and 32 h after injury, n = 6) to establish rat skeletal muscle contusion models. Gene expression data were obtained using second-generation sequencing technology, and differential gene expression analysis, weighted gene co-expression network analysis (WGCNA) and time-dependent expression trend analysis were performed. A total of six sets of biomarkers were obtained: differentially expressed genes at adjacent time points (127 genes), co-expressed genes most associated with wound age (213 genes), hub genes exhibiting time-dependent expression (264 genes), and sets of transcription factors (TF) corresponding to the above sets of genes (74, 87, and 99 genes, respectively). Then, random forest (RF), support vector machine (SVM) and multilayer perceptron (MLP), were constructed for wound age estimation from the above gene sets. The results estimated by transcription factors were all superior to the corresponding hub genes, with the transcription factor group of WGCNA performed the best, with average accuracy rates of 96% for three models' internal testing, and 91.7% for the highest external validation. This study demonstrates the advantages of the indicator screening system based on second-generation sequencing technology and transcription factor level for wound age estimation.

Pathological Characteristics and Classification of Unstable Coronary Atherosclerotic Plaques.

Important forensic diagnostic indicators of sudden death in coronary atherosclerotic heart disease, such as acute or chronic myocardial ischemic changes, sometimes make it difficult to locate the ischemic site due to the short death process, the lack of tissue reaction time. In some cases, the deceased died of sudden death on the first-episode, resulting in difficulty for medical examiners to make an accurate diagnosis. However, clinical studies on coronary instability plaque revealed the key role of coronary spasm and thrombosis caused by their lesions in sudden coronary death process. This paper mainly summarizes the pathological characteristics of unstable coronary plaque based on clinical medical research, including plaque rupture, plaque erosion and calcified nodules, as well as the influencing factors leading to plaque instability, and briefly describes the research progress and technique of the atherosclerotic plaques, in order to improve the study on the mechanism of sudden coronary death and improve the accuracy of the forensic diagnosis of sudden coronary death by diagnosing different pathologic states of coronary atherosclerotic plaques.

Engineering and characterization of carbohydrate-binding modules for imaging cellulose fibrils biosynthesis in plant protoplasts.

Carbohydrate binding modules (CBMs) are noncatalytic domains that assist tethered catalytic domains in substrate targeting. CBMs have therefore been used to visualize distinct polysaccharides present in the cell wall of plant cells and tissues. However, most previous studies provide a qualitative analysis of CBM-polysaccharide interactions, with limited characterization of engineered tandem CBM designs for recognizing polysaccharides like cellulose and limited application of CBM-based probes to visualize cellulose fibrils synthesis in model plant protoplasts with regenerating cell walls. Here, we examine the dynamic interactions of engineered type-A CBMs from families 3a and 64 with crystalline cellulose-I and phosphoric acid swollen cellulose. We generated tandem CBM designs to determine various characteristic properties including binding reversibility toward cellulose-I using equilibrium binding assays. To compute the adsorption (nkon ) and desorption (koff ) rate constants of single versus tandem CBM designs toward nanocrystalline cellulose, we employed dynamic kinetic binding assays using quartz crystal microbalance with dissipation. Our results indicate that tandem CBM3a exhibited the highest adsorption rate to cellulose and displayed reversible binding to both crystalline/amorphous cellulose, unlike other CBM designs, making tandem CBM3a better suited for live plant cell wall biosynthesis imaging applications. We used several engineered CBMs to visualize Arabidopsis thaliana protoplasts with regenerated cell walls using confocal laser scanning microscopy and wide-field fluorescence microscopy. Lastly, we also demonstrated how CBMs as probe reagents can enable in situ visualization of cellulose fibrils during cell wall regeneration in Arabidopsis protoplasts.

Combined metabolomics and machine learning algorithms to explore metabolic biomarkers for diagnosis of acute myocardial ischemia.

Acute myocardial ischemia (AMI) remains the leading cause of death worldwide, and the post-mortem diagnosis of AMI represents a current challenge for both clinical and forensic pathologists. In the present study, the untargeted metabolomics based on ultra-performance liquid chromatography combined with high-resolution mass spectrometry was applied to analyze serum metabolic signatures from AMI in a rat model (n = 10 per group). A total of 28 endogenous metabolites in serum were significantly altered in AMI group relative to control and sham groups. A set of machine learning algorithms, namely gradient tree boosting (GTB), support vector machine (SVM), random forest (RF), logistic regression (LR), and multilayer perceptron (MLP) models, was used to screen the more valuable metabolites from 28 metabolites to optimize the biomarker panel. The results showed that classification accuracy and performance of MLP model were better than other algorithms when the metabolites consisting of L-threonic acid, N-acetyl-L-cysteine, CMPF, glycocholic acid, L-tyrosine, cholic acid, and glycoursodeoxycholic acid. Finally, 17 blood samples from autopsy cases were applied to validate the classification model's value in human samples. The MLP model constructed based on rat dataset achieved accuracy of 88.23%, and ROC of 0.89 for predicting AMI type II in autopsy cases of sudden cardiac death. The results demonstrated that MLP model based on 7 molecular biomarkers had a good diagnostic performance for both AMI rats and autopsy-based blood samples. Thus, the combination of metabolomics and machine learning algorithms provides a novel strategy for AMI diagnosis.

A novel method for determining postmortem interval based on the metabolomics of multiple organs combined with ensemble learning techniques.

Determining postmortem interval (PMI) is one of the most challenging and essential endeavors in forensic science. Developments in PMI estimation can take advantage of machine learning techniques. Currently, applying an algorithm to obtain information on multiple organs and conducting joint analysis to accurately estimate PMI are still in the early stages. This study aimed to establish a multi-organ stacking model that estimates PMI by analyzing differential compounds of four organs in rats. In a total of 140 rats, skeletal muscle, liver, lung, and kidney tissue samples were collected at each time point after death. Ultra-performance liquid chromatography coupled with high-resolution mass spectrometry was used to determine the compound profiles of the samples. The original data were preprocessed using multivariate statistical analysis to determine discriminant compounds. In addition, three interrelated and increasingly complex patterns (single organ optimal model, single organ stacking model, multi-organ stacking model) were established to estimate PMI. The accuracy and generalized area under the receiver operating characteristic curve of the multi-organ stacking model were the highest at 93% and 0.96, respectively. Only 1 of the 14 external validation samples was misclassified by the multi-organ stacking model. The results demonstrate that the application of the multi-organ combination to the stacking algorithm is a potential forensic tool for the accurate estimation of PMI.

Forensic identification of sudden cardiac death: a new approach combining metabolomics and machine learning.

The determination of sudden cardiac death (SCD) is one of the difficult tasks in the forensic practice, especially in the absence of specific morphological changes in the autopsies and histological investigations. In this study, we combined the metabolic characteristics from corpse specimens of cardiac blood and cardiac muscle to predict SCD. Firstly, ultra-high performance liquid chromatography coupled with high-resolution mass spectrometry (UPLC-HRMS)-based untargeted metabolomics was applied to obtain the metabolomic profiles of the specimens, and 18 and 16 differential metabolites were identified in the cardiac blood and cardiac muscle from the corpses of those who died of SCD, respectively. Several possible metabolic pathways were proposed to explain these metabolic alterations, including the metabolism of energy, amino acids, and lipids. Then, we validated the capability of these combinations of differential metabolites to distinguish between SCD and non-SCD through multiple machine learning algorithms. The results showed that stacking model integrated differential metabolites featured from the specimens showed the best performance with 92.31% accuracy, 93.08% precision, 92.31% recall, 91.96% F1 score, and 0.92 AUC. Our results revealed that the SCD metabolic signature identified by metabolomics and ensemble learning in cardiac blood and cardiac muscle has potential in SCD post-mortem diagnosis and metabolic mechanism investigations.

Postmortem Interval Estimation Using Protein Chip Technology Combined with Multivariate Analysis Methods.

OBJECTIVES: To estimate postmortem interval (PMI) by analyzing the protein changes in skeletal muscle tissues with the protein chip technology combined with multivariate analysis methods. METHODS: Rats were sacrificed for cervical dislocation and placed at 16 ℃. Water-soluble proteins in skeletal muscles were extracted at 10 time points (0 d, 1 d, 2 d, 3 d, 4 d, 5 d, 6 d, 7 d, 8 d and 9 d) after death. Protein expression profile data with relative molecular mass of 14 000-230 000 were obtained. Principal component analysis (PCA) and orthogonal partial least squares (OPLS) were used for data analysis. Fisher discriminant model and back propagation (BP) neural network model were constructed to classify and preliminarily estimate the PMI. In addition, the protein expression profiles data of human skeletal muscles at different time points after death were collected, and the relationship between them and PMI was analyzed by heat map and cluster analysis. RESULTS: The protein peak of rat skeletal muscle changed with PMI. The result of PCA combined with OPLS discriminant analysis showed statistical significance in groups with different time points (P<0.05) except 6 d, 7 d and 8 d after death. By Fisher discriminant analysis, the accuracy of internal cross-validation was 71.4% and the accuracy of external validation was 66.7%. The BP neural network model classification and preliminary estimation results showed the accuracy of internal cross-validation was 98.2%, and the accuracy of external validation was 95.8%. There was a significant difference in protein expression between 4 d and 25 h after death by the cluster analysis of the human skeletal muscle samples. CONCLUSIONS: The protein chip technology can quickly, accurately and repeatedly obtain water-soluble protein expression profiles in rats' and human skeletal muscles with the relative molecular mass of 14 000-230 000 at different time points postmortem. The establishment of multiple PMI estimation models based on multivariate analysis can provide a new idea and method for PMI estimation.

Exploration and Practice of the "One Combination, Two Highlights, Three Combinations, Four in One" Innovative Talents Training Mode in Forensic Medicine.

Talent is one of the basic and strategic supports for building a modern socialist country in all aspects. Since the 1980s, the establishment of forensic medicine major and the cultivation of innovative talents in forensic medicine have become hot topics in higher education in forensic medicine. Over the past 43 years, the forensic medicine team of Shanxi Medical University has adhered to the joint education of public security and colleges, and made collaborative innovation, forming a training mode of "One Combination, Two Highlights, Three Combinations, Four in One" for innovative talents in forensic medicine. It has carried out "5+3/X" integrated reform, and formed a relatively complete talent training innovation mode and management system in teaching, scientific research, identification, major, discipline, team, platform and cultural construction. It has made a historic contribution to China's higher forensic education, accumulated valuable experience for the construction of first-class major and first-class discipline of forensic medicine, and provided strong support for the construction of the national new forensic talent training system. The popularization of this training mode is conducive to the rapid and sustainable development of forensic science, and provides more excellent forensic talents for national building, regional social development and the discipline construction of forensic science.

Estimating the time of skeletal muscle contusion based on the spatial distribution of neutrophils: a practical approach to forensic problems.

The study aimed to explore the neutrophil's spatial distributions used to estimate the histological age of contused skeletal muscle, and assessed the accuracy of various indicators, such as the proportion of neutrophils, "neutrophil mean distance," and distribution of neutrophils in areas of "contiguous contour lines." Fifty-five Sprague-Dawley rats were divided randomly into a control group and contusion groups at 1, 1.5, 2, 3, 4, and 6 h, as well as 1, 3, 5, and 15 days, post-injury (n = 5 per group). Nuclei and neutrophils were detected by hematoxylin and eosin (HE) staining and immunohistochemical (IHC) staining. At 0-24 h after injury, the distribution of neutrophils at distances of 100, 200, 300, 400, 500, and 600 µm from adjacent blood vessels was determined, and the best samples were screened to estimate wound age. To estimate wound age as accurately as possible, Fisher discriminant analysis (FDA) of the proportion of neutrophils, neutrophil mean distance, and distribution of neutrophils was performed, and 100.0% and 95.0% of the original and cross-validated cases were correctly classified, respectively. The spatial distribution of neutrophils at different distances from adjacent blood vessels showed a strong correlation with the histological age of contusion skeletal muscle, and the combination of the proportion of neutrophils, neutrophil mean distance, and distribution of neutrophils could be used to accurately estimate wound age.

Fabrication of zinc-loaded silicon carbide nanocomposite for in vitro cell viability and in vivo wound dressing care.

AIM: In this investigation, Zinc-silicon carbide (Zn-SiC) materials were fabricated by a simple approach by using Zn nanoparticles (Zn-NPs) loaded on silicon carbide (SiC) with enhanced antibacterial and healing activity. METHODS: Zn-NPs loaded on SiC fabricated by the DIY laser melting technique. The TEM and Zeta-sizer confirmed the morphology and size of the nanoparticles. The characterisation was done using Fourier transforms infra-red spectroscopy, and X-ray diffraction, thermogravimetric analysis. Further, the fabricated nanoparticles were evaluated for their mechanical properties and biocompatibility under storage conditions. In vivo wound healing was measured by observing a percentage reduction in the wound. RESULTS: Zn-SiC NPs have 54.6 ± 5.25 nm mean particle size, -15.9 ± 2.35 mV zeta potential with 0.187 ± 0.05 polydispersity index. The nanoparticles showed good biocompatibility and in vivo wound healing properties. CONCLUSIONS: These results strongly support the possibility of using these Zn particles loaded on SiC NPs as a promising wound healing agent after caesarean section.

Comparison of Medical Dispute Resolution Mechanisms in China and Abroad.

Medical disputes are one of the common problems concerned by the whole world. All countries and regions have established their own medical dispute resolution mechanisms, in accordance with their own national conditions. Medical dispute identification opinions, as one of the important bases for identifying the responsibilities of both doctors and patients, play a pivotal role in the process of dispute settlement. A reasonable medical dispute resolution mechanism and standardized medical dispute identification model can help resolve disputes flexibly and reduce the conflict between doctors and patients. This paper briefly compares the medical dispute resolution mechanism and identification mode of China and several other representative countries (the United States, Britain, France, Germany, Italy, Japan, etc.), and discusses their respective characteristics and shortcomings, to bring some enlightenment to the medical dispute resolution and identification in our country.

Application of Linear Regression Model of Gpnmb Gene in Rat Injury Time Estimation.

OBJECTIVES: To investigate the effects of injury time, postmortem interval (PMI) and postmortem storage temperature on mRNA expression of glycoprotein non-metastatic melanoma protein B (Gpnmb), and to establish a linear regression model between Gpnmb mRNA expression and injury time, to provide aimed at providing potential indexes for injury time estimation. METHODS: Test group SD rats were anesthetized and subjected to blunt contusion and randomly divided into 0 h, 4 h, 8 h, 12 h, 16 h, 20 h and 24 h groups after injury, with 18 rats in each group. After cervical dislocation, 6 rats in each group were collected and stored at 0 ℃, 16 ℃ and 26 ℃, respectively. The muscle tissue samples of quadriceps femoris injury were collected at 0 h, 12 h and 24 h postmortem at the same temperature. The grouping method and treatment method of the rats in the validation group were the same as above. The expression of Gpnmb mRNA in rat skeletal muscle was detected by RT-qPCR. The Pearson correlation coefficient was used to evaluate the correlation between Gpnmb mRNA expression and injury time, PMI, and postmortem storage temperature. SPSS 25.0 software was used to construct a linear regression model, and the validation group data was used for the back-substitution test. RESULTS: The expression of Gpnmb mRNA continued to increase with the prolongation of injury time, and the expression level was highly correlated with injury time (P<0.05), but had little correlation with PMI and postmortem storage temperature (P>0.05). The linear regression equation between injury time (y) and Gpnmb mRNA relative expression (x) was y=0.611 x+4.489. The back-substitution test proved that the prediction of the model was accurate. CONCLUSIONS: The expression of Gpnmb mRNA is almost not affected by the PMI and postmortem storage temperature, but is mainly related to the time of injury. Therefore, a linear regression model can be established to infer the time of injury.

Research Progress of Automatic Diatom Test by Artificial Intelligence.

Diatom test is the main laboratory test method in the diagnosis of drowning in forensic medicine. It plays an important role in differentiating the antemortem drowning from the postmortem drowning and inferring drowning site. Artificial intelligence (AI) automatic diatom test is a technological innovation in forensic drowning diagnosis which is based on morphological characteristics of diatom, the application of AI algorithm to automatic identification and classification of diatom in tissues and organs. This paper discusses the morphological diatom test methods and reviews the research progress of automatic diatom recognition and classification involving AI algorithms. AI deep learning algorithm can assist diatom testing to obtain objective, accurate, and efficient qualitative and quantitative analysis results, which is expected to become a new direction of diatom testing research in the drowning of forensic medicine in the future.

Comparison among Four Deep Learning Image Classification Algorithms in AI-based Diatom Test.

OBJECTIVES: To select four algorithms with relatively balanced complexity and accuracy among deep learning image classification algorithms for automatic diatom recognition, and to explore the most suitable classification algorithm for diatom recognition to provide data reference for automatic diatom testing research in forensic medicine. METHODS: The "diatom" and "background" small sample size data set (20 000 images) of digestive fluid smear of corpse lung tissue in water were built to train, validate and test four convolutional neural network (CNN) models, including VGG16, ResNet50, InceptionV3 and Inception-ResNet-V2. The receiver operating characteristic curve (ROC) of subjects and confusion matrixes were drawn, recall rate, precision rate, specificity, accuracy rate and F1 score were calculated, and the performance of each model was systematically evaluated. RESULTS: The InceptionV3 model achieved much better results than the other three models with a balanced recall rate of 89.80%, a precision rate of 92.58%. The VGG16 and Inception-ResNet-V2 had similar diatom recognition performance. Although the performance of diatom recall and precision detection could not be balanced, the recognition ability was acceptable. ResNet50 had the lowest diatom recognition performance, with a recall rate of 55.35%. In terms of feature extraction, the four models all extracted the features of diatom and background and mainly focused on diatom region as the main identification basis. CONCLUSIONS: Including the Inception-dependent model, which has stronger directivity and targeting in feature extraction of diatom. The InceptionV3 achieved the best performance on diatom identification and feature extraction compared to the other three models. The InceptionV3 is more suitable for daily forensic diatom examination.

Estimating Postmortem Interval Using Intestinal Microbiota Diversity Based on 16S rRNA High-throughput Sequencing Technology. / 基于16S rRNAé«˜é€šé‡æµ‹åºæŠ€æœ¯çš„è‚ é“èŒç¾¤å¤šæ ·æ€§æŽ¨æ–­æ­»äº¡æ—¶é—´.

OBJECTIVES: To explore the correlation between intestinal microbiota and postmortem interval(PMI) in rats by using 16S rRNA high-throughput sequencing technology. METHODS: Rats were killed by anesthesia and placed at 16 ℃, and DNA was extracted in caecum at 14 time points of 0, 1, 2, 3, 5, 7, 9, 12, 15, 18, 21, 24, 27 and 30 d after death. The 16S rRNA high-throughput sequencing technology was used to detect intestinal microbiota in rat cecal contents, and the results were used to analyze the rat intestinal microbiota diversity and differences. RESULTS: The total number of intestinal microbial communities did not change significantly within 30 days after death, but the diversity showed an upward trend. A total of 119 bacterial communities were significantly changed at 13 time points after death. The models for PMI estimation were established by using partial least squares (PLS) regression at all time points, before 9 days and after 12 days, reaching an R2 of 0.795, 0.767 and 0.445, respectively; and the root mean square errors (RMSEs) were 6.57, 1.96 and 5.37 d, respectively. CONCLUSIONS: Using 16S rRNA high-throughput sequencing technology, the composition and structure of intestinal microbiota changed significantly within 30 d after death. In addition, the established PLS regression model suggested that the PMI was highly correlated with intestinal microbiota composition, showing a certain time series change.

Novel insights into wound age estimation: combined with "up, no change, or down" system and cosine similarity in python environment.

Wound age estimation is a complex, multifactorial issue. It is considered to have great practical significance that combining multi-biomarkers and multi-methods for injury time estimation. We optimized our earlier "up, no change, or down" model by adding data on the expression levels of mRNAs encoding ABHD2, MAD2L2, and ARID5A, and we converted the relative quantitative expression levels of seven genes into a vector rather than a color model. We used Python to derive the cosine similarity (CS) between a test set and the vector matrix; the highest similarity most accurately reflected the injury time. For the optimized model, the internal and external verifications were approximately 0.71 and 0.66, respectively. The good double-blinded results indicated that the model was stable and reliable. In summary, we used a vector matrix and cosine similarities derived by Python to mine the levels of genes expressed in contused skeletal muscle. We are the first to combine several biomarkers and methods for wound age estimation.

Temporal expression of wound healing-related genes inform wound age estimation in rats after a skeletal muscle contusion: a multivariate statistical model analysis.

Although many time-dependent parameters involved in wound healing have been exhaustively investigated, establishing an objective and reliable means for estimating wound age remains a challenge. In this study, 78 Sprague-Dawley rats were divided randomly into a control group and contusion groups at 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, and 48 h post-injury (n = 6 per group). The expression of 35 wound healing-related genes was explored in contused skeletal muscle by real-time polymerase chain reaction. Differences between the groups were assessed by partial least squares discriminant analysis (PLS-DA). The results show that the samples were classified into three groups by wound age (4-12, 16-24, and 28-48 h). A Fisher discriminant analysis model of 14 selected genes was constructed, and 94.9% cross-validated grouped cases were correctly classified. A PLS regression analysis using 14 genes showed reasonable internal predictive validity, with a root mean squared error of cross-validation of approximately 8 h. To examine whether the prediction models were capable of analyzing new (ungrouped) cases, an external validation was carried out using the expression data from an additional 30 rats. Approximately 76.7% of ungrouped cases were correctly classified, which was a lower proportion than that for cross-validation. Similarly, the prediction results of the PLS model showed lower relatively external predictive validity (root mean squared error of prediction = 11 h) than internal predictive validity. Although the prediction results were less accurate than expected, the gene expression modeling and multivariate analyses showed great potential for estimating injury time. These multivariate methods may be valuable when devising future wound time estimation strategies.

Pharmacokinetics of meperidine (pethidine) in rabbit oral fluid: correlation with plasma concentrations after controlled administration.

Oral fluid assays for quantifying drugs are useful in forensic toxicology and drug monitoring. Compared with blood and urine specimens, oral fluid collection is simple, non-invasive, and more difficult to adulterate. Therefore, we investigated whether meperidine and its metabolites could be detected in oral fluid and whether there was a predictable relationship between oral fluid and plasma concentrations. Male New Zealand white rabbits (n = 10) were administered meperidine hydrochloride (20 mg/kg, intravenous). Then, plasma and oral fluid were collected at various time points up to 10 h after administration. We developed a simple and sensitive gas chromatography-mass spectrometry method for the determination of meperidine and normeperidine in oral fluid and plasma. We estimated the apparent pharmacokinetic parameters for meperidine in oral fluid and plasma and determined the ratio and correlation between oral fluid and plasma concentrations. The results demonstrate that this method has excellent specificity, linearity, precision, and recovery. Meperidine and normeperidine were detected in both body fluids; meperidine was the most abundant analyte in oral fluid. The oral fluid-to-plasma drug concentration ratios did not differ significantly over time (p > 0.05). In addition, oral fluid and plasma levels of meperidine and normeperidine were significantly correlated over time (r = 0.713 and 0.725, respectively; p < 0.05). These results provide context for interpreting meperidine and metabolite concentrations in oral fluid and support the utility of oral fluid as an alternative matrix in clinical and forensic testing.