Non-targeted screening workflows for gas chromatography-high-resolution mass spectrometry analysis and identification of biomagnifying contaminants in biota samples

Analysis of chemical contaminants in fish using high resolution mass spectrometry - A review

High resolution mass spectrometry (HRMS) has become an important tool in environmental and food safety analysis. This review highlights how HRMS has been used to analyze chemical contaminants in fish. Measuring and documenting chemical contaminants in fish serves not only as an indicator of environmental conditions but can also monitor the health of these animals and help protect an important source of human food. The incidence and significance of contaminants including veterinary drugs, human drugs and personal care products, pesticides, persistent organic pollutants, per- and poly fluorinated substances, and marine toxins will be reviewed. The advantage of HRMS over traditional MS is its ability to expand the number of compounds that can be detected and identified. This is true whether HRMS is used for targeted analytes, or more broadly for suspect screening and nontargeted analyses. The classes of compounds, types of fish or seafood, options for data acquisition and analysis, and reports of unexpected findings from recent HMRS methods for chemical contaminants in fish are summarized.

Screening Disinfection Byproducts in Arid-Coastal Wastewater: A Workflow Using GC×GC-TOFMS, Passive Sampling, and NMF Deconvolution Algorithm

Disinfection during tertiary municipal wastewater treatment is a necessary step to control the spread of pathogens; unfortunately, it also gives rise to numerous disinfection byproducts (DBPs), only a few of which are regulated because of the analytical challenges associated with the vast number of potential DBPs. This study utilized polydimethylsiloxane (PDMS) passive samplers, comprehensive two-dimensional gas chromatography (GC×GC) coupled with time-of-flight mass spectrometry (TOFMS), and non-negative matrix factorization (NMF) spectral deconvolution for suspect screening of DBPs in treated wastewater. PDMS samplers were deployed upstream and downstream of the chlorination unit in a municipal wastewater treatment plant located in Abu Dhabi, and their extracts were analyzed using GC×GC-TOFMS. A workflow incorporating a multi-tiered, eight-filter screening process was developed, which successfully enabled the reliable isolation of 22 candidate DBPs from thousands of peaks. The NMF spectral deconvolution improved the match factor score of unknown mass spectra to the reference mass spectra available in the NIST library by 17% and facilitated the identification of seven additional DBPs. The close match of the first-dimension retention index data and the GC×GC elution patterns of DBPs, both predicted using the Abraham solvation model, with their respective experimental counterparts-with the measured data available in the NIST WebBook and the GC×GC elution patterns being those observed for the candidate peaks-significantly enhanced the accuracy of peak assignment. Isotopic pattern analysis revealed a close correspondence for 11 DBPs with clearly visible isotopologues in reference spectra, thereby further strengthening the confidence in the peak assignment of these DBPs. Brominated analogues were prevalent among the detected DBPs, possibly due to seawater intrusion. The fate, behavior, persistence, and toxicity of tentatively identified DBPs were assessed using EPI Suite™ and the CompTox Chemicals Dashboard. This revealed their significant toxicity to aquatic organisms, including developmental, mutagenic, and endocrine-disrupting effects in certain DBPs. Some DBPs also showed activity in various CompTox bioassays, implicating them in adverse molecular pathways. Additionally, 11 DBPs demonstrated high environmental persistence and resistance to biodegradation. This combined approach offers a powerful tool for future research and environmental monitoring, enabling accurate identification and assessment of DBPs and their potential risks.

Evaluation of GC-EI&CI-TOFMS for Nontarget Analysis of Industrial Wastewater Using Hydrophilic-Lipophilic-Balanced SPME

The evaluation of nontarget analysis (NTA) techniques for the monitoring of wastewater is important as wastewater is an anthropogenic pollution source for aquatic ecosystems and a threat to human and environmental health. This study presents the proof-of-concept NTA of industrial wastewater samples. A prototype hydrophilic-lipophilic-balanced (HLB) SPME and gas chromatography interfaced with time-of-flight high-resolution mass spectrometry (GC-TOFMS) with electron ionization (EI) and chemical ionization (CI) in parallel are employed. The HLB-SPME consists of a poly(divinylbenzene-co-N-vinylpyrrolidone) structure, allowing the extraction of hydrophilic as well as lipophilic substances. As the combination of parallel CI and EI data provides a comprehensive data set as a unique feature, this study is strongly focused on the compound identification procedure and confidence reporting of exemplary substances. Furthermore, the use of three different CI reagent ions, including [N2H]+/[N4H]+, [H3O]+, and [NH4]+, enables a broad range of analytes to be ionized in terms of selectivity and softness. The complementary information provided by EI and CI data allows a level 3 identification or higher in 69% of cases. The polarity coverage based on the physicochemical properties of the analytes (such as volatility, water solubility, hydrophilicity, and lipophilicity) was visualized by using Henry's law and octanol-water partitioning constants. In conclusion, the presented approach is shown to be valuable for water analysis and allows enhanced and accelerated compound identification compared to utilizing only one type of ionization.

Comprehensive characterization and prioritization of halogenated organic compounds in fish and their implications for exposure

Fish are an important pollution indicator for biomonitoring of halogenated organic compounds (HOCs) in aquatic environments, and HOCs in fish may pose health threats to consumers. This study performed nontarget and comprehensive analyses of HOCs in fish from an e-waste recycling zone by gas chromatography-high-resolution mass spectrometry, and further prioritized their human exposure risks. A total of 1652 formulas of HOCs were found in the fish, of which 1222, 117, and 313 were organochlorines, organobromines, and organochlorine-bromines, respectively. The total concentrations of HOCs were 15.4-18.7 μg/g (wet weight), and organobromines were the predominant (14.1-16.8 μg/g). Of the HOCs, 41 % were elucidated with tentative structures and divided into 13 groups. The estimated total daily exposures of HOCs via dietary consumption of the fish for local adult residents were 3082-3744 ng/kg bw/day. The total exposures were dominated by several groups of HOCs with the following contribution order: polyhalogenated biphenyls and their derivatives > polyhalogenated diphenyl ethers > halo- (H-)alkanes/olefines > H-benzenes > H-dioxins > H-polycyclic aromatic hydrocarbons > H-phenols. The comprehensive characterization and prioritization results provide an overview of the species and distributions of HOCs in edible fish, and propose an inventory of crucial HOCs associated with high exposure risks.

Non-targeted analysis of lipidic extracts by high-resolution mass spectrometry to characterise the chemical exposome: Comparison of four clean-up strategies applied to egg

Biota samples are used to monitor chemical stressors and their impact on the ecosystem and to describe dietary chemical exposure. These complex matrices require an extraction step followed by clean-up to avoid damaging sensitive analytical instruments based on chromatography coupled to mass spectrometry. While interest for non-targeted analysis (NTA) is increasing, there is no versatile or generic sample preparation for a wide range of contaminants suitable for a diversity of biotic matrices. Among the contaminants' variety, persistent contaminants are mostly hydrophobic (mid- to non-polar) and bio-magnify through the lipidic fraction. During their extraction, lipids are generally co-extracted, which may cause matrix effect during the analysis such as hindering the acquired signal. The aim of this study was to evaluate the efficacy of four clean-up methods to selectively remove lipids from extracts prior to NTA. We evaluated (i) gel permeation chromatography (GPC), (ii) Captiva EMR-lipid cartridge (EMR), (iii) sulphuric acid degradation (H2SO4) and (iv) polydimethyl siloxane (PDMS) for their efficiency to remove lipids from hen egg extracts. Gas and liquid chromatography coupled with high-resolution mass spectrometry fitted with either electron ionisation or electrospray ionisation sources operating in positive and negative modes were used to determine the performances of the clean-up methods. A set of 102 chemicals with a wide range of physico-chemical properties that covers the chemical space of mid- to non-polar contaminants, was used to assess and compare recoveries and matrix effects. Matrix effects, that could hinder the mass spectrometer signal, were lower for extracts cleaned-up with H2SO4 than for the ones cleaned-up with PDMS, EMR and GPC. The recoveries were satisfactory for both GPC and EMR while those determined for PDMS and H2SO4 were low due to poor partitioning and degradation/dissociation of the compounds, respectively. The choice of the clean-up methods, among those assessed, should be a compromise that takes into account the matrix under consideration, the levels and the physico-chemical properties of the contaminants.

What is in the fish? Collaborative trial in suspect and non-target screening of organic micropollutants using LC- and GC-HRMS

A collaborative trial involving 16 participants from nine European countries was conducted within the NORMAN network in efforts to harmonise suspect and non-target screening of environmental contaminants in whole fish samples of bream (Abramis brama). Participants were provided with freeze-dried, homogenised fish samples from a contaminated and a reference site, extracts (spiked and non-spiked) and reference sample preparation protocols for liquid chromatography (LC) and gas chromatography (GC) coupled to high resolution mass spectrometry (HRMS). Participants extracted fish samples using their in-house sample preparation method and/or the protocol provided. Participants correctly identified 9-69 % of spiked compounds using LC-HRMS and 20-60 % of spiked compounds using GC-HRMS. From the contaminated site, suspect screening with participants' own suspect lists led to putative identification of on average ∼145 and ∼20 unique features per participant using LC-HRMS and GC-HRMS, respectively, while non-target screening identified on average ∼42 and ∼56 unique features per participant using LC-HRMS and GC-HRMS, respectively. Within the same sub-group of sample preparation method, only a few features were identified by at least two participants in suspect screening (16 features using LC-HRMS, 0 features using GC-HRMS) and non-target screening (0 features using LC-HRMS, 2 features using GC-HRMS). The compounds identified had log octanol/water partition coefficient (KOW) values from -9.9 to 16 and mass-to-charge ratios (m/z) of 68 to 761 (LC-HRMS and GC-HRMS). A significant linear trend was found between log KOW and m/z for the GC-HRMS data. Overall, these findings indicate that differences in screening results are mainly due to the data analysis workflows used by different participants. Further work is needed to harmonise the results obtained when applying suspect and non-target screening approaches to environmental biota samples.

Temporal trends of lipophilic organic contaminants in blue mussel (1994-2017) and eelpout (1994-2017) from the southern Baltic Sea

A time-trend study was carried out for two important Baltic Sea species, blue mussel (1994-2017, 11 samples) and eelpout (1994-2017, 11 samples), to track the changes in levels of regulated persistent organic pollutants (POPs) and show potential increases in the levels of the contaminants of emerging concern (CECs). It was carried out utilizing gas chromatography-high-resolution mass spectrometry (GC-HRMS) based non-target screening (NTS). Data were acquired in two modes - electron ionization (EI) and electron capture negative ion chemical ionization (ECNI) - to widen the contaminant coverage, and treated using a fast semi-automated NTS data processing workflow. The study revealed that >250 tentatively identified compounds show statistically significant temporal trends in Baltic blue mussel and eelpout. A large number of regulated substances, including but not limited to PCBs, DDTs and other organochlorine pesticides (OCPs), chlorobenzenes, and many polybrominated diphenyl ethers (PBDEs), showed significant declining trends, as was expected. Their rates of decline were in good agreement with previously reported data. In contrast, increasing trends were observed for many CECs, some polycyclic aromatic compounds (PAHs), and hydrocarbons. The CEC group included, among others, four compounds, namely, one personal care product ingredient, 2-ethylhexyl stearate, one brominated compound 1,2,3,5-tetrabromobenzene and two intermediates 4-isopropoxyaniline and bilobol dimethyl ether, that were reported in marine biota for the first time to the best of our knowledge. Several compounds, including four CECs and two unknown brominated compounds, showed levels considerably higher than the common legacy pollutants (CB-153 and BDE-99), which might be taken into consideration for future monitoring and risk assessment. In addition, this work revealed the presence of a plethora of organoiodinated compounds that exhibited statistically significant temporal trends in the samples under study, which could be of future interest.

Cracked and shucked: GC-APCI-IMS-HRMS facilitates identification of unknown halogenated organic chemicals in French marine bivalves

High resolution mass spectrometry (HRMS)-based non-target analysis coupled with ion mobility spectrometry (IMS) is gaining momentum due to its ability to provide complementary information which can be useful in the identification of unknown organic chemicals in support of efforts in unraveling the complexity of the chemical exposome. The chemical exposome in the marine environment, though not as well studied as its freshwater counterparts, is not foreign to chemical diversity specially when it comes to potentially bioaccumulative and bioactive polyhalogenated organic contaminants and natural products. In this work we present in detail how we utilized IMS-HRMS coupled with gas chromatographic separation and atmospheric pressure chemical ionization (APCI) to annotate polyhalogenated organic chemicals in French bivalves collected from 25 sites along the French coasts. We describe how we used open cheminformatic tools to exploit isotopologue patterns, isotope ratios, Kendrick mass defect (Cl scale), and collisional cross section (CCS), in order to annotate 157 halogenated features (level 1: 54, level 2: 47, level 3: 50, and level 4: 6). Grouping the features into 11 compound classes was facilitated by a KMD vs CCS plot which showed co-clustering of potentially structurally-related compounds. The features were semi-quantified to gain insight into the distribution of these halogenated features along the French coast, ultimately allowing us to differentiate between sites that are more anthropologically impacted versus sites that are potentially biodiverse.

Thin-layer chromatography coupled with HPLC-DAD/UHPLC-HRMS for target and non-target determination of emerging halogenated organic contaminants in animal-derived foods

In this study, a simple, solvent-saving, and sensitive method was established using high-performance liquid chromatography coupled with high-resolution mass spectrometry to quantitively determine 16 emerging halogenated organic contaminants (HOCs), including polyhalogenated carbazoles, halogenated phenols, and tetrabromobisphenol-A analogs, and to qualitatively identify non-target HOCs in animal-derived food samples. The sample extracts were physically frozen to remove most lipids and further purified by thin-layer chromatography according to the targets polarity. The method detection limit and method quantification limit of 16 HOCs were in the range of 0.003-9.168 and 0.010-30.601 ng·g-1 dry weight, respectively. The recoveries were within 69.1-111.0 %, the intra/inter-day precisions were 0.1-6.1 % and 0.1-6.7 %, and the matrix effects were between -12.1 and 10.8 %, all within the acceptable range. Finally, 16 HOCs were detected in nine actual samples in range of not detected-307.22 ng·g-1 dry weight. Moreover, five bromides and two chlorides were identified by using non-target analysis in animal-derived foods.

Comprehensive non-target screening of biomagnifying organic contaminants in the Baltic Sea food web

High-resolution mass spectrometry (HRMS) based non-target screening (NTS) is a powerful approach for the simultaneous determination of multiple environmental contaminant classes in complex biota samples. In this study, trophic biomagnification factor (TMF) directed NTS was performed to find and (tentatively) identify known, emerging, and new chemical contaminants that are persistent and biomagnify in Baltic Sea biota. The investigated food web included seven species: one filter feeder (blue mussel, Mytilus edulis), two fish (eelpout, Zoarces viviparous; herring, Clupea harengus), two marine mammals (harbor porpoise, Phocoena phocoena; grey seal, Halichoerus grypus) and two birds (guillemot, Uria aalge; white-tailed sea eagle, Haliaeetus albicilla). The NTS procedure included extraction with organic solvent mixtures, two-step high-resolution gel permeation chromatography clean-up, Florisil® fractionation, gas chromatography (GC) HRMS analysis in electron ionization (EI) and electron capture negative ion chemical ionization (ECNI) modes, and NTS data processing. The latter was performed differently for the EI and ECNI data: the EI data were treated using a flexible and highly automated TMF-directed NTS workflow, whereas the ECNI data were treated with a simpler and less automated workflow that specifically screened for brominated compounds. The two workflows collectively revealed biomagnification (statistically significant TMF values) of >250 tentatively identified compounds, including legacy persistent organic pollutants (POPs), such as PCBs and PCB-related compounds, DDT and its metabolites, and organochlorine pesticides (OCPs), contaminants of emerging concern (CECs), and halogenated natural products (HNPs). Among the tentatively identified CECs, nine have not previously been reported in environmental biota samples. These included four polymer additives (used as antioxidants, rubber additives or plasticizers) and two cosmetic product additives (ethyl myristate and isopropyl palmitate). The CECs should be prioritized for future structure verification and quantification using reference standards.

Enhancing spectral quality in complex environmental matrices: Supporting suspect and non-target screening in zebra mussels with ion mobility

Identification of bioaccumulating contaminants of emerging concern (CECs) via suspect and non-target screening remains a challenging task. In this study, ion mobility separation with high-resolution mass spectrometry (IM-HRMS) was used to investigate the effects of drift time (DT) alignment on spectrum quality and peak annotation for screening of CECs in complex sample matrices using data independent acquisition (DIA). Data treatment approaches (Binary Sample Comparison) and prioritisation strategies (Halogen Match, co-occurrence of features in biota and the water phase) were explored in a case study on zebra mussel (Dreissena polymorpha) in Lake Mälaren, Sweden's largest drinking water reservoir. DT alignment evidently improved the fragment spectrum quality by increasing the similarity score to reference spectra from on average (±standard deviation) 0.33 ± 0.31 to 0.64 ± 0.30 points, thus positively influencing structure elucidation efforts. Thirty-two features were tentatively identified at confidence level 3 or higher using MetFrag coupled with the new PubChemLite database, which included predicted collision cross-section values from CCSbase. The implementation of predicted mobility data was found to support compound annotation. This study illustrates a quantitative assessment of the benefits of IM-HRMS on spectral quality, which will enhance the performance of future screening studies of CECs in complex environmental matrices.

A time-trend guided non-target screening study of organic contaminants in Baltic Sea harbor porpoise (1988-2019), guillemot (1986-2019), and white-tailed sea eagle (1965-2017) using gas chromatography-high-resolution mass spectrometry

The rate of decline in regulated persistent organic pollutant (POP) concentrations in Baltic Sea biota has leveled off in recent years, with new contaminants frequently being discovered. There is, therefore, a need for comprehensive approaches to study occurrence and temporal trends of a wide range of environmental contaminants, including legacy POPs, contaminants of emerging concern (CECs), and new contaminants. In the current work, non-target screening (NTS) workflows were developed and used for, to the best of our knowledge, the first time-trend directed NTS of biota using gas chromatography-high-resolution mass spectrometry (GC-HRMS). To maximize contaminant coverage, both electron ionization (EI) and electron capture negative ion chemical ionization (ECNI) were used. The EI data were treated using highly automated workflows to find, prioritize, and tentatively identify contaminants with statistically significant temporal trends. The ECNI data were manually processed and reviewed prior to time-trend analysis. Altogether, more than 300 tentatively identified contaminants were found to have significant temporal trends in samples of Baltic guillemot, harbor porpoise, or white-tailed sea eagle. Significant decreases were found for many regulated chemicals, as could be expected, such as PCBs, polychlorinated terphenyls, chlorobenzenes, toxaphenes, DDT, other organochlorine pesticides, and tri- and tetra- bromodiphenyl ethers (BDEs). The rate of decline of legacy POPs agreed well with data reported from targeted analyses. Significant increases were observed for small polycyclic aromatic hydrocarbons, heptaBDEs, CECs, and terpenes and related compounds. The CECs included, among others, one plasticizer tributyl acetylcitrate (ATBC), two antioxidants 2,6-bis(1,1-dimethylethyl)phenol and 2,6-bis(tert-butyl)-4-(4-morpholinyl-methyl)phenol, and two compounds used in polymer production, trimethyl isocyanurate and 2-mercaptobenzothiazole, which had not previously been reported in biota. Their increased concentrations in biota indicate increased use and release. The increase in ATBC may be linked to increased use of it as a substitute for di-2-ethylhexyl phthalate (DEHP), which has been phased out over the last decade.

Biomagnification and Temporal Trends of New and Emerging Dechloranes and Related Transformation Products in Baltic Sea Biota

To enhance knowledge of the environmental distribution and temporal trends of dechloranes and their transformation products (TPs) we performed suspect screening of Baltic Sea biota (eelpout, herring, harbor porpoise, guillemot and white-tailed sea eagle). Evaluation of new and "digitally frozen" gas chromatography/high-resolution mass spectrometry data revealed 31 compounds: five dechloranes (Dechlorane [Mirex], Dechlorane 602, Dechlorane 603, and syn-/anti-Dechlorane Plus [DP]), three isomers, and 23 TPs. Six new Dechlorane 603 TPs and two new DP TPs were detected, including one hydroxy-TP. Some TPs occurred at much higher concentrations than the parent compounds (e.g., Dechlorane 603 TPs were >10-fold more abundant than their parent). Concentrations of contaminants in the most contaminated species (white-tailed sea eagle) changed little over the period 1965-2017. Slow declines were detected for most compounds (median, 2% per year), although concentrations of DP and DP-TPs increased by 1% per year. Ten contaminants biomagnify, and the trophic magnification factors for TPs of Mirex, Dechlorane 602 and Dechlorane 603 (8.2 to 17.8) were similar to the parent compounds (6.6 to 12.4) and higher than that of DP (2.4, nonsignificant). The results are discussed in relation to the current review of DP for potential listing under the Stockholm Convention on POPs.

Application of a novel prioritisation strategy using non-target screening for evaluation of temporal trends (1969-2017) of contaminants of emerging concern (CECs) in archived lynx muscle tissue samples

Most environmental monitoring studies of contaminants of emerging concern (CECs) focus on aquatic species and target specific classes of CECs. Even with wide-scope target screening methods, relevant CECs may be missed. In this study, non-target screening (NTS) was used for tentative identification of potential CECs in muscle tissue of the terrestrial top predator Eurasian lynx (Lynx lynx). Temporal trend analysis was applied as a prioritisation tool for archived samples, using univariate statistical tests (Mann-Kendall and Spearman rank). Pooled lynx muscle tissue collected from 1969 to 2017 was analysed with an eight-point time series using a previously validated screening workflow. Following peak detection, peak alignment, and blank subtraction, 12,941 features were considered for statistical analysis. Prioritisation by time-trend analysis detected 104 and 61 features with statistically significant increasing and decreasing trends, respectively. Following probable molecular formula assignment and elucidation with MetFrag, two compounds with increasing trends, and one with a decreasing trend, were tentatively identified. These results show that, despite low expected concentration levels and high matrix effects in terrestrial species, it is possible to prioritise CECs in archived lynx samples using NTS and univariate statistical approaches.

Application of inert gas-mediated ionization for qualitative screening of chlorinated aromatics in house dust by comprehensive two-dimensional gas chromatography-high-resolution time-of-flight mass spectrometry

The development of highly selective and sensitive analytical methods for the nontarget screening of persistent organic pollutants such as halogenated compounds in environmental samples is a challenging task. Soft ionization mass spectrometry has emerged as a powerful technique for obtaining essential molecular information, and it is expected to reveal compounds that remain hidden with conventional fragmentation techniques such as electron ionization (EI). In this study, a soft ionization method based on electron capture negative ionization using an inert gas was developed for the nontarget screening of chlorinated aromatics in environmental samples and was applied to comprehensive two-dimensional gas chromatography-high-resolution time-of-flight mass spectrometry (GC × GC-HRToFMS). In particular, argon (Ar) and helium (He) were evaluated as inert moderating gases, and were compared against the conventional methane (CH₄). The optimal ionization conditions, including the flow rate and ion source temperature, were investigated based on the molecular ion intensities of highly chlorinated aromatics decachlorobiphenyl and octachlorodibenzofuran. Ar-mediated soft ionization provided the best sensitivity to molecular ions among the three gases at a low flow rate (0.1 mL min^-1) and low ion source temperature, and more selective detection of molecular ions (i.e., less fragmentation) was obtained with the inert gases than with CH₄. This method is also applicable to other chlorinated aromatics such as tetra- to nonachlorobiphenyls, tetra- to heptachlorinated dibenzofurans, pentachlorobenzene, and hexachlorobenzene. To demonstrate the applicability of the proposed method to a wide range of chlorinated aromatics in environmental samples, both Ar-mediated soft ionization and conventional EI were applied to GC × GC-HRToFMS for analysis of a crude extract of house dust. Soft ionization enabled the selective and sensitive detection of molecular ions for minor amounts of chlorinated aromatics, even in complex matrices. Furthermore, the extracted ion chromatograms of halide anions (Cl^- or Br^-) were useful for screening other chlorinated or brominated compounds in the environmental samples. The results suggest that combining information on halide anions obtained by soft ionization and the structural information provided by EI would constitute a powerful approach for the comprehensive identification of chlorinated aromatics.

Non-targeted analysis of vulgarisins by using collisional dissociation mass spectrometry for the discovery of analogues from Prunella vulgaris

Vulgarisins are members of diterpenoids with rare 5/6/4/5 ring skeleton from Prunella vulgaris Linn. (P. vulgaris). Their molecular scaffolds comprise different hydroxylation and degree of esterification. Vulgarisins have attracted many attentions in the fields of food and medicine for their potent bioactivities. Firstly, four reference compounds were analyzed by higher-energy collisional dissociation mass spectrometry (HCD MS/MS) and the fragmentation patterns for molecular scaffold were summarized. And then, a high-performance liquid chromatography/electrospray ionization/high-resolution mass spectrometry (HPLC-ESI-HR-MS) method was adopted to investigate the P. vulgaris extracts. Finally, the proposed analysis results were successfully applied to facilitate the discovery of the vulgarisins analogues from P. vulgaris. For the four reference compounds, the sodium adduct was the predominate ion in full scan. A specific fragmentation pathway of [M+Na]⁺ ions leads to produce diagnostic ions of vulgarisins at m/z 325 under HCD, which was formed through consecutive-side chains lost. Twenty-three diterpenoids, including 18 vulgarisins analogues, were identified or tentatively characterized in the botanical extracts of P. vulgaris based on their elemental constituents and characteristic fragment ion profiles. Two new vulgarisins analogues in the plant were isolated and their structures were illustrated based on extensive spectroscopic analysis using ¹D and ²D nuclear magnetic resonance (NMR) spectroscopy. The HCD MS/MS method, including the profiles of the diagnostic ions induced by characteristic fragmentation, is an effective technique for the discovery of vulgarisins analogues in P. vulgaris. The expected fragmentation pattern knowledge will also facilitate the analysis of other natural products.