Literature DB >> 35647425

Analysis of Derivatized Wall Teichoic Acids Confirms that a Mutation in Phage-Resistant Listeria monocytogenes Impacts Rhamnose Decoration.

Danielle M Trudelle1, Daniel W Bryan1, Shahla Ray1, John P Munafo1, Thomas G Denes1.   

Abstract

Listeria monocytogenes is a Gram-positive foodborne pathogen that causes listeriosis, an illness that may result in serious health consequences or death. Wall teichoic acids (WTAs) are external cell wall glycopolymers that play many biological roles. Here, the WTA composition was determined for several phage-resistant mutant strains of L. monocytogenes. The strains included wild-type (WT) L. monocytogenes 10403S, and three phage-resistant mutant strains derived from 10403S, consisting of two well-characterized strains and one with unknown impact on cell physiology. Several WTA monomers were prepared from WT 10403S, as analytical standards. The WTA monomer fraction was then isolated from the mutant strains and the corresponding per-trimethylsilylated derivatives were analyzed by gas chromatography-flame ionization detection. WTA monomer, GlcNAc-Rha-Rbo, was detected in 10403S, and not detected in the phage-resistant strains known to lack rhamnose and N-acetylglucosamine; although the expected monomers GlcNAc-Rbo and Rha-Rbo were detected, respectively. GlcNAc-Rha-Rbo was also detected in strain UTK P1-0001, which is known to impact phage adsorption through an undetermined mechanism, albeit at a lower intensity than the WT 10403S, which is consistent with partial loss of function through truncation in RmlC protein. WTA monomers were also detected in an unpurified cell pellet, demonstrating that the method employed in this study can be used to rapidly screen L. monocytogenes without laborious WTA purification. This study lays the groundwork for future studies on WTA compositional analysis to support genomic data, and serves as a foundation for the development of new rapid methods for WTA compositional analysis.
© 2022 The Authors. Published by American Chemical Society.

Entities:  

Year:  2022        PMID: 35647425      PMCID: PMC9134232          DOI: 10.1021/acsomega.1c07403

Source DB:  PubMed          Journal:  ACS Omega        ISSN: 2470-1343


Introduction

Listeria monocytogenes is a foodborne pathogen that is known for its potential to cause listeriosis, a disease that may lead to serious illness or death in the young, old, immunocompromised, or pregnant.[1] Major outbreaks caused by listeriosis have been associated with various ready-to-eat (RTE) foods, including dairy products, produce, and deli meats.[2−4] Listeriosis currently ranks as one of the most deadly and costly foodborne illnesses in the U.S.[5−7]L. monocytogenes is a Gram-positive bacterium, indicated by the thick peptidoglycan (PG) layer comprising approximately 30–40% of its cell wall.[8] This substantial layer of PG supports the presence of glycopolymers (CWGs), which are attached either to the PG itself or to the cell membrane. In Listeria, CWG attached to the cell membrane are lipoteichoic acids (LTAs), while CWGs attached to the PG are wall teichoic acids (WTAs).[9] LTA and WTA share some important functions in the cell, including supporting cell division and morphology, biofilm formation, ion regulation, and virulence.[10,11] Although Listeria can survive without teichoic acids, or even a cell wall as L-form bacteria, such cells require specific growth conditions to do so, and lose functions associated with the cell wall and CWG.[12] Unlike LTA, Listeria WTA display considerable variation within their glycosylation units, acting as the O-antigens for the cell, and are major determinants of the different serotypes (STs).[13−15]L. monocytogenes consists of at least 13 different ST, including 1/2a, 1/2b, 1/2c, 3a, 3b, 3c, 4a, 4b, 4ab, 4c, 4d, 4e and 7.[16] In this study, we focus on the model serotype 1/2a strain 10403S and its mutant derivatives. Over a ten year span, serotype 1/2a strains were responsible for approximately 40% of L. monocytogenes outbreaks in the U.S.[17] The ST 1/2, 3, and 7 strains display a type I WTA structure, consisting of a ribitol (Rbo) backbone with either N-acetylglucosamine (GlcNAc), rhamnose (Rha), or a hydroxyl group (OH) bound to carbons 2 and/or 4. Specifically, the 1/2 ST has a GlcNAc substituent at C2 and a Rha substituent at C4.[15] ST 3 WTA units are affected by a mutation in at least one of the genes necessary for WTA rhamnosylation, and therefore are lacking Rha in this structure.[15] ST 7 strains lack GlcNAc and Rha due to a mutation in at least one of the genes required for WTA rhamnosylation, as well as in at least one of the genes required for WTA N-acetylglucosaminylation.[18] Individual type I WTA units are bound at C1 and C5 of the Rbo by a phosphodiester linkage to form polymer chains approximately 21 units long.[15,19] A pivotal role of WTA is that they serve as the receptors for bacteriophages (or “phages”), viruses that exclusively infect bacteria.[9,20] To replicate, phages inject genetic material into the bacterial host, where the cellular machinery of the host is utilized to produce progeny phage. At the end of this process, the host cell lyses from within, releasing the progeny phage into the environment.[20] A common form of bacterial resistance to Listeria phage infection occurs when the bacteriophages are unable to adsorb to their host. This process has been shown to occur through accumulation of mutations in genes impacting WTA glycosylation; in L. monocytogenes ST 1/2a strains, mutations that affect WTA rhamnosylation or N-acetylglucosaminylation significantly reduce the ability of bacteriophages to adsorb and infect.[18,21,22] Further, mutations causing loss of Rha in WTAs resist infection from almost all Listeria phages tested against it,[23] with only one characterized wild-type phage exception;[24] although phages were able to gain the ability to infect this type of resistant mutant through the process of in vitro evolution.[25] As phages are used as biocontrols in the food industry,[26] improving our knowledge of how L. monocytogenes can resist phage infection is critical for the long-term success of phage-based food safety applications.[27] The structures of Listeria WTA have been analyzed using an array of different analytical methodologies. For the analysis of WTA structural components, methods including gel filtration chromatography,[28,29] permethylation combined with gas chromatography[29,30] with flame ionization detection[31] or mass spectrometry (MS)[19] have been previously reported. For analysis of molecular connectivity within the WTA structure, Smith degradation as well as general oxidation and reduction reactions have been used,[28,29,32] and nuclear magnetic resonance (NMR) spectra have been produced to determine anomeric configurations within the structure of the WTA.[15,19,28,32] More recently, major advances in WTA have been reported, including the employment of ultraperformance liquid chromatography (UPLC) with electrospray ionization (ESI) and tandem MS (MS/MS),[15] and ESI-MS/MS alone.[33] Although the methods used to determine WTA structures have progressively advanced over recent decades, unfortunately, there has not been a significant reduction in the time to obtain pure isolated WTA monomers used in the analysis. The current methods require extraction, purification, and hydrolysis of the WTA polymer before the sample can be analyzed. Although the analysis of the pure WTA is fast employing UPLC-ESI-MS/MS, the instrumentation is not as common in commercial laboratories as gas chromatography-flame ionization detection (GC-FID), nor are WTA monomer reference materials commercially available. The published methods to isolate and purify WTA monomers are time-consuming and laborious. Accordingly, it would be helpful to develop methods with rapid sample preparation that employ common analytical instrumentation; however, pure WTA reference materials, which are currently lacking, are needed for using this approach. The primary goal of this present investigation was to screen several phage-resistant mutant strains of L. monocytogenes for the presence of specific WTA monomers (Table ; Figure ) using per-trimethylsilyl derivatization of the samples and employing gas chromatography - flame ionization detection (GC-FID). Two of the phage-resistant mutants, FSL D4-0014 and FSL D4-0119, have been well-characterized with known impact on WTA’s and were included in the study as controls.[22] The third phage-resistant mutant, UTK P1-0001, has been genetically characterized. This genetic analysis suggested that the strain would express a truncated RmlC protein, which is known to impact rhamnosylation of WTA’s; however, the strain shows a unique resistance pattern to phages through a mechanism of adsorption inhibition, and the impact of the mutation on WTA composition is unknown.[23] Accordingly, the three specific objectives of the present study were to (1) isolate WTA monomers as an analytical standards; (2) analyze four different strains of L. monocytogenes by GC-FID, including (i) wild-type 10403S, (ii) phage-resistant mutant strain FSL D4-0119 (lacking Rha), (iii) phage-resistant mutant strain FSL D4-0014 (lacking GlcNAc), and (iv) phage-resistant mutant strain (UTK P1-0001), which possesses a truncated RmlC protein causing an undetermined impact on cell physiology; and (3) probe the feasibility of using the GC-FID method for WTA analysis of a single colony of L. monocytogenes.
Table 1

L. monocytogenes Strains and Features

L. monocytogenes StrainFeaturesReference
Wild-type Laboratory Strain
10403SLineage II; 1/2a serotype (GlcNAc and Rha in WTA)(2)
Mutant Strains
FSL D4-001410403S mutant; nonsense mutation in LMRG_00541; deficiency of GlcNAc in WTA; susceptible to phage LP-048 and resistant to LP-125a(22)
FSL D4-011910403S mutant; nonsense mutation in LMRG_00542; deficiency of Rha in WTA; resistant to phages LP-048 and LP-125a(22)
UTK P1-000110403S mutant; frameshift mutation caused by a deletion in LMRG_00544; truncated RmlC protein; susceptible to phage LP-125 and resistant to LP-048a(23)

Resistance shown to occur through mechanism of adsorption inhibition.

Figure 1

Chemical structures of (1) GlcNAc-Rha-Rbo, (2) Rha-Rbo, and (3) GlcNAc-Rbo.

Chemical structures of (1) GlcNAc-Rha-Rbo, (2) Rha-Rbo, and (3) GlcNAc-Rbo. Resistance shown to occur through mechanism of adsorption inhibition.

Results and Discussion

Preparation of WTA Analytical Standards

The objective of this study was to screen several phage-resistant mutant strains of L. monocytogenes for their WTA composition using GC-FID. To accomplish this goal, first, the WTA monomer, GlcNAc-Rha-Rbo, was purified (Figure ). Then the sample was further hydrolyzed, yielding a mixture of GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo as a mixed analytical standard (Figure ). GC-FID analysis of Rbo-5TMS in partial hydrolysates was used to determine a comparable degree of hydrolysis in the samples (Figure ).
Figure 2

EI mass spectrum of trimethylsilyl derivatives of (A) Rbo and (B) Rha, depicting the characteristic fragmentation pattern.

EI mass spectrum of trimethylsilyl derivatives of (A) Rbo and (B) Rha, depicting the characteristic fragmentation pattern. A mixture of GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo was then per-trimethylsilylated, yielding GlcNAc-Rha-Rbo-9TMS, GlcNAc-Rbo-7TMS, and Rha-Rbo-7TMS, respectively (Figure ), for GC–MS and GC-FID analysis. For comparison of WTA monomers in each sample, the areas of each WTA monomer were normalized based on the area of Rbo-5TMS in the 10403S sample (Figure ).
Figure 3

EI mass spectrum of trimethylsilyl derivatives of (A) GlcNAc-Rha-Rbo, (B) Rha-Rbo, and (C) GlcNAc-Rbo, depicting the characteristic fragmentation pattern.

Figure 4

GC-FID analysis of Rbo-5TMS in partial hydrolysate of 10403S, FSL D4-0119, FSL D4-0014, and UTK P1-0001 WTA fractions. GC-FID analysis of Rbo-5TMS in partial hydrolysates was used to determine a comparable degree of hydrolysis in the samples.

EI mass spectrum of trimethylsilyl derivatives of (A) GlcNAc-Rha-Rbo, (B) Rha-Rbo, and (C) GlcNAc-Rbo, depicting the characteristic fragmentation pattern. GC-FID analysis of Rbo-5TMS in partial hydrolysate of 10403S, FSL D4-0119, FSL D4-0014, and UTK P1-0001 WTA fractions. GC-FID analysis of Rbo-5TMS in partial hydrolysates was used to determine a comparable degree of hydrolysis in the samples. The preparation of the WTA reference material from 10403S proceeded as follows. The yields denoted after each purification step were calculated from the prior step. L. monocytogenes, wild-type 10403S, was grown on brain heart infusion (BHI) media (5 L) to a density of 0.8–0.9 OD600. The cells were harvested and autoclaved yielding the autoclaved cell pellet (ACP) (17.5 g). The ACP was then mechanically lysed, yielding isolated cell wall material (CWM) (8.5 g, 47.7% yield). The CWM was then subjected to enzyme treatment and washed, yielding the cell wall carbohydrate fraction (CWCF) (160 mg, 1.9% yield). The CWCF was then hydrolyzed and dialyzed, yielding the crude WTA polymer (cWTA polymer) (39 mg, 24.4% yield). The cWTA polymer was then further purified by anion-exchange chromatography, yielding purified WTA polymer (pWTA polymer) (3 mg, 7.69% yield). The pWTA polymer was then subjected to hydrolysis, using hydrogen fluoride (HF), and lyophilized, yielding purified WTA monomer, GlcNAc-Rha-Rbo (1 mg, 33.3% yield). The purity of GlcNAc-Rha-Rbo was >98%, as confirmed by liquid chromatography-mass spectrometry (LC-MS), high resolution electrospray ionization–time-of-flight mass spectrometry (HRESI-TOFMS), GC-FID, and GC–MS. After the initial isolation of GlcNAc-Rha-Rbo, the material was further partially hydrolyzed. The acid-catalyzed hydrolysis was monitored by GC–MS, and finally yielded a mixture of GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo.

Analytical Confirmation of WTA Analytical Standards

The WTA monomers were subjected to LC-MS, HRESI-TOFMS, GC-FID, and GC–MS for structural confirmation. For GlcNAc-Rha-Rbo, a single peak was observed by LC-MS analysis, and upon TMS derivatization and GC-FID and GC-MS analysis, a single chromatographic peak was also observed. Upon LC-MS analysis, a base ion at m/z 500.2 (100, [M – H]−) and an ion at m/z 1001.5 (18, [2M – H]−) was observed, both consistent with a molecular weight of 501.2 Da. Upon HRESI-TOFMS analysis, a sodium adduct at [M + Na]+m/z 524.1949 was observed (calculated for NaC19H35NO14, m/z 524.1950). These results are consistent with the analytical data for the purified GlcNAc-Rha-Rbo, previously reported as the WTA monomer unit of a ST 1/2a strain of L. monocytogenes.[33] The EI-MS spectrum of the per-trimethylsilylated derivative of GlcNAc-Rha-Rbo, GlcNAc-Rha-Rbo-9TMS, was also in agreement with GlcNAc-Rha-Rbo (Figure ). The retention index (RI) for GlcNAc-Rha-Rbo-9TMS on a DB-5 column was calculated as RI = 3515. After initial isolation, the material was further hydrolyzed yielding a mixture of GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo. Accordingly, the WTA mixture was employed as a reference material to generate RIs of the per-trimethylsilylated derivatives of the WTAs for GC-FID method development. Upon LC-MS analysis of Rha-Rbo, a base ion at m/z 321.1 [100, M + Na]+ and an ion at m/z 298.1 [100, M – H]− was observed, both consistent with a molecular weight of 298.3 Da. Upon HRESI-TOFMS analysis of Rha-Rbo, a sodium adduct at [M + Na]+m/z 321.1150 was observed (calculated for NaC11H22O9, m/z 321.1156). The EI-MS spectrum of the per-trimethylsilylated derivative of Rha-Rbo, Rha-Rbo-7TMS, was also in agreement with Rha-Rbo (Figure ). The RI for Rha-Rbo-7TMS on a DB-5 column was calculated as RI = 2487. Upon LC-MS analysis of GlcNAc-Rbo, a base ion at m/z 378.1 [100, M + Na]+ and an ion at m/z 354.1 [100, M – H]− was observed, both consistent with a molecular weight of 355.3 Da. Upon HRESI-TOFMS analysis of GlcNAc-Rbo, a sodium adduct at [M + Na]+m/z 378.1376 was observed (calculated for NaC13H25NO10, m/z 378.1370. The EI-MS spectrum of the per-trimethylsilylated derivative of GlcNAc-Rbo, GlcNAc-Rbo-7TMS, was also in agreement with GlcNAc-Rbo (Figure ). The RI for GlcNAc-Rbo-7TMS on a DB-5 column was calculated as RI = 2845. The results from these experiments allowed for the development of a GC-FID method to separate all three of the derivatized WTAs with good chromatographic resolution.

Analysis of WTA Monomers in L. monocytogenes Serotype 1/2a Mutants

GC-FID was employed to confirm the presence or absence of wall teichoic acid monomers, including the complete GlcNAc-Rha-Rbo or partial GlcNAc-Rbo and Rha-Rbo variants, in the following strains: (1) L. monocytogenes wild-type strain 10403S (Rha+, GlcNAc+), (2) mutant strain FSL D4-0014 (Rha+, GlcNAc–), (3) mutant strain FSL D4-0119 (Rha–, GlcNAc+), and (4) mutant strain UTK P1-0001 (truncated RmlC protein). As anticipated, GlcNAc-Rha-Rbo was detected in wild-type strain 10403S as the GlcNAc-Rha-Rbo-9TMS derivative. Also, GlcNAc-Rha-Rbo was detected in UTK P1-0001 (truncated RmlC protein), albeit at levels ∼11% those of 10403S (Figure ). In addition, Rha-4TMS and Rha-Rbo-7TMS were detected in both strains; however, Rha and Rha-Rbo were also present at lower levels in the UTK P1-0001 strain (∼13% and ∼11% those of 10403S, respectively). UTK P1-0001 was previously found to have a mutation in a gene coding for RmlC,[23] an epimerase that is involved in the biosynthesis of Rha.[34] This mutation was shown to cause a premature stop codon in the rmlC gene that would only result in the loss of the last three amino acids of RmlC, therefore its function is unlikely to be completely inhibited,[23] which is consistent with our observations of less observed Rha in the strain. Additionally, this strain has previously demonstrated binding to the Listeria phage LP-125, which requires both GlcNAc and Rha for adsorption.[22,23] The results from this study suggest that RmlC in UTK P1-0001 has reduced activity, resulting in less rhamnosylation of the WTA; this would explain the unusual phage resistance phenotype where LP-125 (requiring both GlcNAc and Rha for binding) is able to adsorb to the strain, whereas LP-048 (requiring only Rha) does not appear to efficiently adsorb.[23] Thus, LP-048 binding correlates with the observed reduction in WTA Rha composition in UTK P1-0001. Future studies should explore the potential for this type of mutation to impact biofilm formation, virulence, and resistance to antimicrobials, as rhamnose and other WTA structures have been implicated in these processes for L. monocytogenes.[13,35−37]
Figure 5

Comparison of Rha, Rha-Rbo, GlcNAc-Rbo, and GlcNAc-Rha-Rbo in the partial hydrolysate of 10403S, FSL D4-0119, FSL D4-0014, and UTK P1-0001 WTA fractions. All chromatographic peaks are per-trimethylsilylated derivatives analyzed by GC-FID. The factor denoted in the upper left corner of each trace was calculated by dividing the area of the analyte in the mutant sample by the area of the analyte in the 10403S sample. For comparison of WTA monomers in each mutant strain, the areas of each WTA monomer were normalized based on the area of Rbo-5TMS in the 10403S strain. ND; not detected.

Comparison of Rha, Rha-Rbo, GlcNAc-Rbo, and GlcNAc-Rha-Rbo in the partial hydrolysate of 10403S, FSL D4-0119, FSL D4-0014, and UTK P1-0001 WTA fractions. All chromatographic peaks are per-trimethylsilylated derivatives analyzed by GC-FID. The factor denoted in the upper left corner of each trace was calculated by dividing the area of the analyte in the mutant sample by the area of the analyte in the 10403S sample. For comparison of WTA monomers in each mutant strain, the areas of each WTA monomer were normalized based on the area of Rbo-5TMS in the 10403S strain. ND; not detected. In the FSL D4-0014 (Rha+, GlcNAc–) mutant strain, GlcNAc-Rha-Rbo was not detected, which is consistent with previous results showing that FSL D4-0014 is deficient in GlcNAc, and therefore would be lacking GlcNAc-Rha-Rbo as a WTA monomer. In contrast, Rha-4TMS and Rha-Rbo-7TMS were detected in FSL D4-0014. FSL D4-0014 is a 10403S derived mutant, with a loss of function mutation in a glycosyltransferace gene, lmo1079, responsible for addition of GlcNAc in WTA.[18,22] Phage binding and wheat germ agglutination assays previously confirmed the expected result. In the FSL D4-0119 (Rha–, GlcNAc+) mutant strain, neither GlcNAc-Rha-Rbo nor Rha was detected; although GlcNAc-Rbo was detected at levels ∼25% higher than in 10403S (Figure ). FSL D4-0119 is a 10403S derived mutant with a loss of function mutation (nonsense mutation) in another glycosyltransferase, lmo1080, which is known to be required for rhamnosylation of WTA.[18,22] Consistent with the results, sequencing data and phage spot tests/adsorption assays demonstrated that FSL D4-0119 is deficient in Rha.[22] On the basis of the genomic data, each one of the phage-resistant mutant strains have a mutation in a glycosyltransferase associated with WTA decoration. It has been shown that these genes are consistent hotspots for mutations conferring phage resistance.[18,22,25] In this present study, the WTA fraction of each of the mutant strains tested here was purified and analyzed by GC-FID, confirming the presence (UTK P1-0001; lower levels than 10403S) or absence (FSL D4-0014, FSL D4-0119; not detected) of GlcNAc-Rha-Rbo. This compositional data supports the genomic data and the approach used in this study may be employed as a tool to characterize different strains of L. monocytogenes. To develop high-throughput WTA compositional analysis, more rapid methods are needed to keep pace with the advancements in genomic data sets.

Analysis of WTA Monomers in Less Processed Samples

In the present study, GC-FID confirmation of the presence or absence of GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo from the purified WTA monomer fraction of each L. monocytogenes strain was successful (Figure ); however, the sample preparation for each strain was extensive. The sample preparation required the culture of large quantities of bacteria (5 L; each) and the purification was laborious, including mechanical, chemical, and enzymatic, treatments coupled with chromatographic purification. To assess the feasibility of using the GC-FID method to identify GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo using smaller quantities of bacteria with less processing, samples were collected at four steps throughout the WTA purification process using the wild-type 10403S strain. The four samples included, in decreasing purity, (1) purified GlcNAc-Rha-Rbo, (2) crude WTA polymer (cWTA polymer), (3) CWCF (cell wall carbohydrate fraction), and (4) ACP (autoclaved cell pellet). The three in-process samples, cWTA polymer, CWCF, and ACP were subjected to HF hydrolysis prior to derivatization. All samples, were derivatized with N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA) with 1% trimethylchlorosilane (TMCS) prior to GC-FID analysis. Interestingly, GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo were detected in all the in-process samples analyzed, including the ACP sample that were only steam-killed cells (Figure ).
Figure 6

Streamlining sample preparation for rapid analysis of trimethylsilyl derivatives of WTA monmers by GC-FID. (A) Purified WTA monomer (pWTA monomer), (B) crude WTA polymer (cWTA polymer), (C) CWCF, and (D) ACP. The strain used was 10403S.

Streamlining sample preparation for rapid analysis of trimethylsilyl derivatives of WTA monmers by GC-FID. (A) Purified WTA monomer (pWTA monomer), (B) crude WTA polymer (cWTA polymer), (C) CWCF, and (D) ACP. The strain used was 10403S. The calculated signal-to-noise ratio (SNR) for GlcNAc-Rbo-7TMS, Rha-Rbo-7TMS, and GlcNAc-Rha-Rbo-9TMS in the ACP sample was 923:1, 33:1, and 409:1, respectively (see the Supporting Information). The SNR for the per-trimethylsilylated WTA components using the current unoptimized method ranged from 11 to 309 times greater than the standard 3:1 value for the limit of detection (LOD) and 3 to 41 times greater than the standard 10:1 value for the limit of quantification (LOQ). The results of this experiment suggest that the analytical method may be further optimized to serve as a faster and less expensive alternative to the currently published analytical methods for WTA detection and quantitation. Derivatization and GC-FID analysis of ACP samples require significantly less sample, time, reagents, and equipment than purifying WTA samples and analysis by UPLC-ESI-MS/MS.[14] We show that this method can reduce the volume of broth-grown culture of L. monocytogenes needed for analysis from several liters to several milliliters. In addition, the method eliminates time-consuming steps, costly reagents, and specialized equipment. Processing steps that are eliminated include lysing cells with a French Press, enzymatic treatment of the lysates, anion-exchange chromatography, and phosphate testing of eluent fractions, while still detecting WTA monomers above the standard LOQ. However, the limitation to the current method is that more analytical standards would need to be commercially available to develop the method for additional WTA monomers with variable composition. This approach would offer an analytical alternative to rapid protein-binding methods, such as Sumrall et al’s use of phage-derived recombinant affinity proteins to glycotype Listeria.[38] Whereas protein-binding methods may be particularly useful for rapid indirect determination of exposed cell-surface physiology, the analytical methods described here would provide direct quantitative measures of specific WTA components. GC-FID analysis of per-trimethylsilylated derivatives of WTAs with different configurations may also be chromatographically resolved, including WTAs containing diastereomers such as glucose and galactose, and enantiomers, if chiral capillary columns such as a 2,3-diacetoxy-6-O-tert-butyl dimethylsilyl γ-cyclodextrin (Rt-γDEXsa) or a 2,3-di-O-methyl-6-O-tert-butyl dimethylsilyl β-cyclodextrin (Rt-βDEXsm) are employed; however, additional evaluations are needed to test these hypotheses. Furthermore, the results from the GC-FID analysis of the CWCF sample indicates there is little benefit in purifying the sample beyond enzymatic treatment of the cell lysate; by stopping here significant time and resources are saved, and the resulting peak intensity is still about 50% the intensity of the fully purified WTA sample. To develop a robust analytical method that can discriminate a collection of WTA acid monomers with differing carbohydrate moieties and structural configurations, more WTA reference materials would need to be available as analytical standards. In conclusion, in this study, we used GC-FID, GC–MS, and LC-MS to confirm the composition of wall teichoic acids in phage-resistant mutants of L. monocytogenes. Interestingly, we show that the unique phage resistance phenotype observed in 10403S derived mutant UTK P1-0001 is likely due to a significantly diminished relative abundance of Rha decoration in the wall teichoic acids, indicating the relevance of partial loss of function mutations in phage-host interactions. Further, we show that the process of observing wall teichoic acid monomers can be significantly streamlined by analyzing unpurified samples by GC-FID. Taken together, the work presented here provides a roadmap for the rapid characterization of Listeria WTA composition; this can provide valuable information on key aspects of Listeria biology such as determining phage susceptibility patterns, virulence potential, subtyping information, and basic cell physiology.

Materials and Methods

Bacterial Growth Conditions

Working stocks of L. monocytogenes strains (Table ) were stored at −80 °C in BHI (Becton Dickinson, Sparks, MD) broth with 15% (wt/vol) glycerol. L. monocytogenes was streaked onto 1.5% (wt/vol) BHI agar and then incubated at 30 °C for approximately 24 h. Erlenmeyer flasks (125 mL) containing 30 mL of BHI broth were used for overnight (ON) cultures and inoculated with three similar sized colonies from a streak plate stored at 4 °C for less than 2 weeks prior. The inoculated BHI broth was then incubated for 16 h at 30 °C shaking at 160 rpm. ON culture was added to Erlenmeyer flasks 1/4 full with sterile BHI broth in a 1:100 ratio (ON culture to BHI) and incubated at 30 °C shaking at 160 rpm. Cells were grown to an OD600 between approximately 0.8 and 1.0 (GENESYS 30 Visible Light Spectrophotometer; Thermo Scientific, Waltham, MA), then autoclaved for 30 min at 121.1 °C. Flasks were cooled to approximately 4 °C in an ice water bath, then the culture was aliquoted into centrifuge bottles and centrifuged (Beckman J2-HS; Beckman Coulter Life Sciences, Indianapolis, IN) at 7000g for 10 min at 4 °C to pellet cells. The ACPs were collected and frozen for storage at −20 °C.

Chemicals and Reagents

Anhydrous pyridine (reagent grade), BSTFA with 1% TMCS, formic acid, glycerol, glycine, HCl, methanol (HPLC grade), NaCl, NaOH, proteinase K, sodium dodecyl sulfate, Tris base, and hydrofluoric acid (48–51%) were purchased from Thermo Fisher Scientific (Waltham, MA). A mixture of n-alkanes C9–C18 was purchased from Phenomenex (Torrance, CA) and n-alkanes C19, C20, and a mixture of C21–C40 were purchased from Millipore Sigma (Burlington, MA). BHI was purchased from Becton Dickinson (Sparks, MD). MgSO4 was purchased from Acros Organics (Geel, Belgium). DNase and RNase were purchased from Alfa Aesar (Tewksbury, MA). A phosphate standard solution was purchased from Merck (Darmstadt, Germany).

Purification of WTA Monomer Reference Materials

Cell Lysis

Analytical standards of the WTAs (GlcNAc-Rha-Rbo, GlcNAc-Rha, and Rha-Rbo) were isolated following a previously reported procedure with some modifications.[14] The WTA monomers, GlcNAc-Rha and Rha-Rbo, were generated through acid catalyzed hydrolysis of GlcNAc-Rha-Rbo. Frozen cell pellets of L. monocytogenes, wild-type 10403S, previously reported to contain WTA monomer GlcNAc-Rha-Rbo were thawed to room temperature (RT) and resuspended in a saline magnesium (SM) buffer to a density of approximately 0.75 g cells/mL. Cells were lysed by at least two passages through a French Press Pressure Cell (French Press Cell Disrupter; Thermo Electron Corporation, Milford, MA) at 270 MPa. Lysed cells were centrifuged (Eppendorf 5804 R; Eppendorf, Hamburg, Germany) at 1400g for 5 min to remove unbroken cells. The supernatant was collected and centrifuged (Beckman J2-HS; Beckman Coulter Life Sciences, Indianapolis, IN) at 20,000g for 30 min at 4 °C to recover cell walls. Pellets were collected from the supernatant until no more visible solid material remained (9 collections total per strain) which were then washed twice with sterile ultrapure water (20,000g for 30 min at 4 °C). The CWM was pooled into Nalgene Oak Ridge tubes (Thermo Scientific, Waltham, MA), combined and frozen for storage at −20 °C.

Cell Wall Purification

Frozen cell wall materials were thawed to RT, weighed, and resuspended in 10 mM Tris–Cl (pH 7.6) for a combined volume of 36 867 μL. The DNase (Alfa Aesar, Tewksbury, MA) working solution was prepared by mixing DNase powder (lyophilized by manufacturer in 2.5 mM calcium acetate and 2.5 mM magnesium sulfate) with sterile ultrapure water to a concentration of 10 mg/mL. The RNase (Alfa Aesar, Tewksbury, MA) working solution was prepared by mixing lyophilized RNase powder with 100 mM Tris–Cl (pH 7.6) to a concentration of 10 mg/mL. The proteinase K (Fisher Scientific, Fair Lawn, NJ) working solution was prepared by mixing lyophilized proteinase K powder with 10 mM Tris–Cl (pH 7.6) to a concentration of 10 mg/mL. Cell wall materials were mixed with DNase and RNase working solutions (376 μL each) together with each enzyme at a final concentration of 100 μg/mL, then inverted 20 times and incubated at 25 °C for 3.5 h with two inversions. Following this, 380 μL of the proteinase K working solution was added for a final concentration of 100 μg/mL, inverted 20 times, and incubated at 25 °C for 2 h with inversions every 30 min. After enzyme treatments, cell walls were pelleted by centrifugation at 20,000g for 30 min at 4 °C (Beckman J2-HS; Beckman Coulter Life Sciences, Indianapolis, IN). The supernatant was discarded, and pellets were stored overnight at 4 °C. Following this, pellets were resuspended with 30 mL of a 4% (w/v) sodium dodecyl sulfate (SDS) solution and aliquoted into glass tubes (∼10–12 mL per tube). Tubes were incubated in water for 30 min at 100 °C. After cooling to RT, sample aliquots were recombined into Nalgene Oak Ridge tubes and SDS-insoluble material (cell pellet) was collected by centrifugation at 20,000g for 30 min at 20 °C (Beckman J2-HS; Beckman Coulter Life Sciences, Indianapolis, IN). Detergent was removed after washing the pellet five times with sterile ultrapure water at 20,000g for 30 min at 20 °C. The resulting CWCF was then resuspended in 5 mL of sterile ultrapure water and transferred into 50 mL centrifuge tubes, frozen on an angle at −20 °C, then lyophilized (VirTis Advantage Plus EL-85; SP Scientific, Gardiner, NY) and stored at −20 °C with desiccant.

Wall Teichoic Acid Extraction

The lyophilized carbohydrate fraction of treated cells was mixed with 25 mM glycine/HCl buffer (pH 2.5) in Reacti-Vials (Thermo Fisher Scientific Inc., Waltham, MA) and hydrolyzed for 10 min at 100 °C. After cooling to RT, samples were centrifuged (Avanti J-26 XP; Beckman Coulter Life Sciences, Indianapolis, IN) at 30,000g for 30 min at 4 °C to pellet insoluble materials. The supernatant was collected, and the pellet was resuspended in the same buffer. Hydrolysis and centrifugation were repeated twice; all collected supernatant was pooled and dialyzed (20 mL D-Tube Dialyzer Mega, MWCO 3.5 kDa; MilliporeSigma, Burlington, MA) at 4 °C against 2 L of ultrapure water for approximately 24 h (with one change of water at 12 h) to remove the buffer. The cWTA polymer solution was then frozen on a slant at −20 °C, lyophilized, and stored at −20 °C with a desiccant.

Wall Teichoic Acid Polymer Purification

The cWTA polymer (10 mg) was dissolved in a starting buffer (750 μL of 10 mM Tris–HCl, pH 7.5) and manually loaded onto the ÄKTA pure (GE Healthcare Bio-Sciences AB, Uppsala, Sweden) LC system. WTA purification was performed with anion-exchange chromatography using a HiTrap DEAE FF Column (5 mL; GE Healthcare Bio-Sciences AB, Uppsala, Sweden). The column was first equilibrated with two column volumes of the starting buffer at a flow rate of 5 mL/min, then 20 column volumes of fractions were collected at a flow rate of 1 mL/min into tubes by elution using a linear gradient of 0–1 M NaCl solution. Starting and elution buffers were filtered through a 0.45 μm Nylon membrane (Whatman-GE Healthcare, Buckinghamshire, United Kingdom) before use in the ÄKTA system. All glassware used for buffer preparation and storage was acid washed with a 10% HCl solution and rinsed with deionized (DI) and ultrapure water before use to prevent phosphate contamination.

Phosphate Standard Curve Preparation

A working stock solution of phosphate standard with a concentration of 10 mg/L PO4 was prepared from a phosphate standard solution (Merck, Darmstadt, Germany) with a concentration of 1000 mg/L PO4. From the working stock solution, 5 mL phosphate standards were prepared in concentrations ranging from 0 to 5 mg/L PO4. A phosphate test kit (Spectroquant; Merck, Darmstadt, Germany) was applied to each standard as well as an ultrapure water blank as per manufacturers’ instructions. Absorbencies of each standard were read in 10 mm cuvettes with a spectrophotometer (GENESYS 30 Visible Spectrophotometer; Thermo Fisher Scientific, Madison, WI) at a wavelength of 690 nm. A standard curve was developed using Microsoft Excel (Version 1811) to establish a linear regression formula (y = 0.1613x + 0.0094) for determining unknown phosphate values based on sample absorbencies at 690 nm.

Determination of Wall Teichoic Acid Containing Fractions

Fractions obtained after WTA purification were tested for UV activity using 1 mL samples aliquoted into acid-washed quartz cuvettes. Absorbance was read at a wavelength of 205 nm using the NanoDrop One Spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA) against a blank containing a 1:1 mixture of 10 mM Tris–HCl, pH 7.5 and 1 M NaCl. A 30 μL subsample was taken from fractions showing absorbency at 205 nm, diluted to a total volume of 5 mL with ultrapure water, and treated with decomposition reagent (NANOCOLOR NanOx Metal; Macherey-Nagel, Düren, Germany) as per manufacturer’s instructions. A blank of ultrapure water was treated the same way. A phosphate test kit (Spectroquant; Merck, Darmstadt, Germany) was then applied to the treated subsamples as well as the blank as per manufacturer’s instructions. Absorbencies were read by spectrophotometry (GENESYS 30 Visible Light Spectrophotometer; Thermo Scientific, Waltham, MA) at a wavelength of 690 nm to calculate phosphate concentration using the formula obtained from the phosphate standard curve (y = 0.1613x + 0.0094). Fractions calculated as having over 25 mg/L PO4 were dialyzed (20 mL D-Tube Dialyzer Mega, MWCO 3.5 kDa; MilliporeSigma, Burlington, MA) at 4 °C against 2 L ultrapure water for 24 h (with one change of water at 12 h) to remove the buffer. The pWTA polymer solution was then frozen on a slant at −20 °C, lyophilized, and stored at −20 °C with desiccant.

Hydrofluoric Acid Hydrolysis

Samples (2 mg; pWTA polymer) were subjected to hydrolysis using 200 μL of HF (48–51%) at 0 °C for 24 h prior to evaporation over NaOH pellets in a chamber under vacuum. The reaction was monitored yielding the first intact GlcNAc-Rha-Rbo monomer. Then the reaction was allowed to proceed for a longer time, yielding a mixture of GlcNAc-Rha and Rha-Rbo as partial hydrolysis products. Samples were subjected to a vacuum until HF evaporation was complete; the sample was then mixed with 500 μL of ultrapure water, frozen at −20 °C, lyophilized, and stored at −20 °C with a desiccant. GC-FID analysis of Rbo-5TMS in partial hydrolysates was used to determine a comparable degree of hydrolysis in the samples.

GlcNAc-Rha-Rbo

C19H35NO14; HR-ToF-MS (ESI+) m/z 524.1949 ([M + Na]+, measured); m/z 524.1950 ([M + Na]+, calculated for NaC19H35NO14; ESI+-MS, m/z 524.1 [100, 2 M + Na]+, 204.1 (25), 502.3 [18, M + H]+, 356.2 (18); ESI–-MS, m/z 500.2 [100, M – H]−, 1001.5 [18, 2 M – H]−, 336.2 (12), 276.0 (5), 202.1 (3). GlcNAc-Rha-Rbo-9TMS: MS (EI) m/z (%) 73 (100), 204 (97), 173 (85), 217 (73), 330 (59), 147 (46), 103 (33), 420 (30), 363 (30), 273 (21), 451 (10), 240 (8), 115 (6), 492 (2) (Figure ); RI on DB-5, 3515.

Rha-Rbo

C11H22O9; HR-ToF-MS (ESI+) m/z 321.1150 ([M + Na] +, measured); m/z 321.1156 ([M + Na] +, calculated for NaC11H22O9; ESI+-MS, m/z 321.1 [100, M + Na]+, 299.1 [16, M + H]+, 175 (8); ESI–-MS, m/z 298.1 [100, M – H]−, 595.2 [12, 2 M – H]−. Rha-Rbo-7TMS: MS (EI) m/z (%) 204 (100), 73 (96), 217 (49), 147 (47), 129 (39), 103 (21), 243 (16), 363 (14), 273 (12), 333 (10), 423 (5), 115 (5) (Figure ); RI on DB-5, 2487.

GlcNAc-Rbo

C13H25NO10; HR-ToF-MS (ESI+) m/z 378.1376 ([M + Na]+, measured); m/z 378.1370 ([M + Na] +, calculated for NaC13H25NO10; ESI+-MS, m/z 378.1 [100, M + Na]+, 356.2 [18, M + H]+, 226.1 (10); ESI–-MS, m/z 354.1 [100, M – H]−, 707.1 [18, 2 M – H]−. GlcNAc-Rbo-7TMS: MS (EI) m/z (%) 73 (100), 173 (80), 217 (58), 330 (50), 147 (45), 129 (40), 103 (29), 246 (19), 420 (14), 271 (7) (Figure ); RI on DB-5, 2845.

Liquid Chromatography-Mass Spectrometry (LC-MS)

LC-MS analysis of the purified WTA monomer of 10403S was performed with an Agilent 1260 series HPLC system (Agilent Technologies Inc., Santa Clara, CA). The system was equipped with an autosampler, a BIN Pump SL binary pump, a TCC SL thermostated column compartment, and a DADSL diode array detector, interfaced to a 6410 triple-quadrupole LC-MS mass selective detector equipped with an API-ESI ionization source. Prior to injection, the sample was dissolved in methanol to a concentration of 1 mg/mL. Chromatographic separations for 10 μL injection volumes were performed using a Gemini column (250 × 4.6 mm i.d.; 5.0 μm particle size) (Phenomenex, Torrance, CA). The column temperature was set at 25 °C and operated at a 1.0 mL/min flow rate. DI water with 0.1% formic acid (A) and methanol with 0.1% formic acid (B) were employed in the binary mobile phase with a linear gradient of 5–55% over 50 min; 55–90% over 5 min; elution at 90% for 5 min, followed by re-equilibration over 10 min. Data acquisition and analysis were performed using Mass Hunter Workstation Data Acquisition, Qualitative Analysis, and Quantitative Analysis software. LC-MS analysis was performed in both negative and positive ion mode with ionization parameters set at capillary voltage, 3.5 kV; nebulizer pressure, 35 psi; drying gas flow, 13.0 mL/min; drying gas temperature, 350 °C; and mass scan range, m/z 300–2000. Quantitative analysis of the sample was performed in negative ion mode with the same ionization parameters as described above.

Sample Preparation for GC-FID and GC-MS Analysis

Hydrofluoric acid hydrolysis of samples was conducted as described above. 50 μL of 1,3,5-trihydroxybenzene (400 μg/mL DI H20) was added as an internal standard. Samples (10 mg; each of each strain), including in-process ACP, CWCF, and cWTA polymer (for the streamlining experiment) were subjected to hydrolysis using 200 μL of HF (48–51%) at 0 °C for 20 h prior to evaporation over NaOH pellets in a chamber under vacuum. Samples were subjected to a vacuum until HF evaporation was complete; the sample was then mixed with 500 μL of ultrapure water, frozen at −20 °C, lyophilized, and stored at −20 °C with a desiccant, prior to derivatization and GC-FID and GC-MS analysis.

Sample Derivatization

Prior to GC-FID and GC–MS analysis, standards and samples were derivatized similarly to the methods used by Munafo et al.[39] with modifications. Approximately 1 mg of the purified WTA monomer sample for each strain and approximately 1 mg each of the streamlining samples (i.e., ACP, CWCF, and cWTA polymer) were derivatized in Reacti-Vials (Thermo Fisher Scientific Inc., Waltham, MA) at 70 °C for 1 h using a mixture of 7 parts anhydrous pyridine and 3 parts BSTFA with 1% TMCS for a total volume of 100 μL (Thermo Scientific, Bellefonte, PA). Standards of Rha and Rbo were derivatized similarly, using 1 mg of each with a total volume of 1 mL derivatizing reagents (Figure ). Samples were then analyzed by GC-FID and GC–MS. For comparison of WTA monomers in each mutant strain, the areas of each WTA monomer were normalized based on the area of Rbo-5TMS in the 10403S sample using the following equation: Areanormlized = (Analytesample × Rbo-5TMS10403s)/Rbo-5TMSsample.

Gas Chromatography-Flame Ionization Detection (GC-FID)

GC-FID (6890 Series; Agilent Technologies, Santa Clara, CA) analysis was conducted similarly to the methods used by Munafo et al.[39] with some modifications. GC-FID was performed by manual injection (Hamilton Company, Bonaduz, Switzerland) with 1 μL of the derivatized sample and a split ratio of 1:10. The inlet temperature was 250 °C. Helium was used as the carrier gas at a flow rate of 1.5 mL/min. The oven was set to an initial temperature of 80 °C (held for 1 min) with ramp of 6 °C/min to 300 °C and held for 15 min. The total run time was 53.67 min. The column used was an HP-5 with capillary size 30.0 m × 0.32 mm × 0.25 μm (Agilent). The detector temperature was 250 °C. Data was analyzed using GC ChemStation Rev. A 10.02 [1757] software (Agilent). Linear RIs were determined for each analyte using the retention times of the analyte and n-alkanes (C9–C40) by linear interpolation.

Gas Chromatography–Mass Spectrometry (GC–MS)

GC–MS was performed on an Agilent 6890 series gas chromatograph (Agilent Technologies, Santa Clara, CA, USA) coupled to an Agilent 5973 mass spectrometer detector. The capillary column used for chromatographic separation was a fused silica DB-5 column (30 m × 0.25 mm × 0.25 μm; Agilent). A 1 μL split/spitless injection (1:1 split) was made by an autosampler using a 10 μL syringe. The inlet temperature was 250 °C. Helium was used as a carrier gas with a constant flow of 1.5 mL/min. The oven temperature was initially held at 80 °C for 1 min followed by an increase in temperature at a rate of 6 °C/min until the oven temperature reached 280 °C and held at this temperature for 90 min. The mass spectrometer detector was coupled to the GC via a transfer line heated at 280 °C and operated in electron ionization (EI) mode at 70 eV. The total run time was 124.338 min. The detector scan range was set to m/z 50–800.

High Resolution Electrospray Ionization-Time-of-Flight Mass Spectrometry (HRESI-TOFMS)

High resolution mass spectra were recorded on a Bruker BioTOF II ESI under the following conditions: source temperature, 150 °C; acceleration voltage, 8500; mass resolution, 10 000 fwhm; scan range, m/z 100–1400; drying gas, N2.
  36 in total

1.  Cross-resistance to phage infection in Listeria monocytogenes serotype 1/2a mutants.

Authors:  Danielle M Trudelle; Daniel W Bryan; Lauren K Hudson; Thomas G Denes
Journal:  Food Microbiol       Date:  2019-06-05       Impact factor: 5.516

2.  Listeria monocytogenes L-forms respond to cell wall deficiency by modifying gene expression and the mode of division.

Authors:  Simone Dell'Era; Carmen Buchrieser; Elisabeth Couvé; Barbara Schnell; Yves Briers; Markus Schuppler; Martin J Loessner
Journal:  Mol Microbiol       Date:  2009-06-23       Impact factor: 3.501

3.  RmlC, the third enzyme of dTDP-L-rhamnose pathway, is a new class of epimerase.

Authors:  M F Giraud; G A Leonard; R A Field; C Berlind; J H Naismith
Journal:  Nat Struct Biol       Date:  2000-05

4.  Multistate outbreak of listeriosis associated with cantaloupe.

Authors:  Jeffrey T McCollum; Alicia B Cronquist; Benjamin J Silk; Kelly A Jackson; Katherine A O'Connor; Shaun Cosgrove; Joe P Gossack; Susan S Parachini; Neena S Jain; Paul Ettestad; Mam Ibraheem; Venessa Cantu; Manjiri Joshi; Tracy DuVernoy; Norman W Fogg; James R Gorny; Kathryn M Mogen; Charlotte Spires; Paul Teitell; Lavin A Joseph; Cheryl L Tarr; Maho Imanishi; Karen P Neil; Robert V Tauxe; Barbara E Mahon
Journal:  N Engl J Med       Date:  2013-09-05       Impact factor: 91.245

5.  l-Rhamnosylation of wall teichoic acids promotes efficient surface association of Listeria monocytogenes virulence factors InlB and Ami through interaction with GW domains.

Authors:  Filipe Carvalho; Sandra Sousa; Didier Cabanes
Journal:  Environ Microbiol       Date:  2018-07-31       Impact factor: 5.491

6.  The cell wall binding domain of Listeria bacteriophage endolysin PlyP35 recognizes terminal GlcNAc residues in cell wall teichoic acid.

Authors:  Marcel R Eugster; Martina C Haug; Simona G Huwiler; Martin J Loessner
Journal:  Mol Microbiol       Date:  2011-08-02       Impact factor: 3.501

7.  Epidemic listeriosis associated with Mexican-style cheese.

Authors:  M J Linnan; L Mascola; X D Lou; V Goulet; S May; C Salminen; D W Hird; M L Yonekura; P Hayes; R Weaver
Journal:  N Engl J Med       Date:  1988-09-29       Impact factor: 91.245

8.  Learning from listeria: safer food for all.

Authors:  Peter J Salama; Peter K Ben Embarek; Jayshree Bagaria; Ibrahima S Fall
Journal:  Lancet       Date:  2018-06-09       Impact factor: 79.321

Review 9.  Location, synthesis and function of glycolipids and polyglycerolphosphate lipoteichoic acid in Gram-positive bacteria of the phylum Firmicutes.

Authors:  Nathalie T Reichmann; Angelika Gründling
Journal:  FEMS Microbiol Lett       Date:  2011-03-25       Impact factor: 2.742

10.  Rapid analysis of Listeria monocytogenes cell wall teichoic acid carbohydrates by ESI-MS/MS.

Authors:  Marcel R Eugster; Martin J Loessner
Journal:  PLoS One       Date:  2011-06-30       Impact factor: 3.240

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.