Salmonella Enteritidis Isolate Harboring Multiple Efflux Pumps and Pathogenicity Factors, Shows Absence of O Antigen Polymerase Gene

Background

Salmonella enterica is one of the most important causes of gastrointestinal infection in humans, being the great majority of infections related to the consumption of poultry meat and eggs (Foley and Lynne, 2008; EFSA/ECDC, 2015).

In animals, infections caused by serotype Enteritidis are rarely responsible for severe disease with animals frequently becoming asymptomatic carriers, except in the case of young chicks and poults, where outbreaks exhibiting clinical disease are often accompanied by high mortality rates (Foley et al., 2008, 2013). Indeed, S. enterica subsp. enterica serovar Enteritidis (S. Enteritidis) has been responsible for severe disease in industrial poultry farming facilities worldwide, posing a potential hazard for public health (Lutful Kabir, 2010).

In order to be infectious, Salmonella needs to adapt to different niches and conditions, where virulence and heavy-metal-tolerance factors play an important role, through co-selection events and the formation of pathogenicity islands, respectively (Hensel, 2004; Medardus et al., 2014). Furthermore, antibiotic resistance determinants can also facilitate their survival, with ubiquitous chromosomally encoded efflux mechanisms, playing an important role in both intrinsic, and acquired multidrug resistance. Other resistance mechanisms, such as changes in the membrane permeability, enzymatic modification, and target alterations may increase the levels of bacterial resistance, contributing to the success of the infection (Poole, 2004; Delmar et al., 2014; Li et al., 2015).

Both antibiotic susceptibility determination and serotyping constitute very useful tools for the epidemiologic classification of S. enterica isolates. Indeed, in S. enterica, the resistance rates fluctuate according to the serotype and with the antibiotic (Clemente et al., 2015). Classically, serotyping is based on the antigenic reactivity of lipolysaccharide (O antigen) and flagellar proteins (H antigen), followed by a designation using names or formulas (Grimont and Weill, 2007). In this study, we aimed to analyze the genome of a S. Enteritidis isolate responsible for omphalitis in chicks, exploring the molecular features associated with antibiotic resistance and pathogenicity, as well as the ability to spread the respective determinants.

Methods

Bacterial isolate, antibiotic susceptibility testing, and serotyping

The isolate (LV60) was recovered from a sample collected from the yolk sac of a chick with omphalitis, under the scope of the “Salmonella National Control Programme in food-producing animals and food of animal origin for bacteriological diagnosis, serotype identification and antibiotic susceptibility testing.” The guidelines of the Commission Decision (CD), 2007/407/EC were followed. LV60 was tested for its antimicrobial resistance through the determination of minimum inhibitory concentrations (MICs) using the agar dilution method, as previously described (Clemente et al., 2013) and according to the European Committee on Antimicrobial Susceptibility Testing (EUCAST) guidelines (http://www.eucast.org/). Briefly, a panel of 11 antibiotic compounds was tested in a 2-fold concentration series over the following ranges: ampicillin and tetracycline (0.5–64 μg/mL), gentamicin and trimethoprim (0.25–32 μg/mL), ciprofloxacin (0.008–8 μg/mL), cefotaxime (0.06–8 μg/mL), nalidixic acid and streptomycin (2–512 μg/mL), chloramphenicol (2–256 μg/mL), florfenicol (1–128 μg/mL) and sulphamethoxazole (8–1024 μg/mL). The epidemiological cut-off values recommended by EUCAST to Salmonella spp. were used for the interpretation of susceptibility testing results. Quality control was performed using the Escherichia coli ATCC 25922 strain. LV60 isolate was then serotyped by the slide agglutination method for its O and H antigens using the method of Kauffman-White scheme (Grimont and Weill, 2007).

Whole genome sequencing (WGS), assembly, and annotation

Genomic DNA was extracted using DNeasy Blood and Tissue Kit (Qiagen), and DNA quantification was performed by Qubit Fluorometric Quantitation (Life Technologies), according to with the manufacturer's instructions. The genome was sequenced using a double strategy of 454 (Roche) and MiSeq (Illumina) sequencing.

Five hundred nanograms of bacterial DNA were fragmented by nebulization, followed by adaptor ligation to create double stranded DNA libraries and sequenced on a 454 GS FLX Titanium according to the standard manufacturer's instructions (Roche-454 Life Sciences). The second genome library was prepared from 1 ng of genomic DNA using the Nextera XT DNA Sample Preparation Kit (Illumina) and sequenced on the Illumina MiSeq sequencer (Illumina) using paired-end 2 × 150 bp reads.

First quality evaluation of raw read sequences and their corresponding quality values were assigned by the FastQC software. Reads were then trimmed and filtered according to quality criteria, and de novo assembled with Ray, version 2.3.1 (Boisvert et al., 2010). Contigs were searched for identity through blastn (http://blast.ncbi.nlm.nih.gov/Blast.cgi) against the nr/nt NCBI database to identify the closest bacterial genome and/or plasmid. Therefore, LV60 genome was mapped against the bacterial genome of S. Enteritidis strain p125109 and its plasmid (NC_011294 and HG970000, respectively) using GS Mapper version 2.9 (Roche). Additionally SNV (single nucleotide variants) and structural variants were also detected with the GS Mapper (Roche, version 2.9).

Structural and functional annotation was performed using PGP (Prokaryotic Genome Prediction) (Egas et al., 2014), an in-house developed pipeline. Taxonomy identification was performed by BLASTP search against the NCBI GenBank non-redundant (nr) database of the 16 s rRNA sequence gene, identified in the previous step and confirmed using RNAmmer v1.2 (Lagesen et al., 2007).

The final data was submitted in the DDBJ/EMBL/GenBank databases, using the Sequin software tool (http://www.ncbi.nlm.nih.gov/Sequin/). This dataset, which includes files in Genbank (LIHI01.1.gbff.gz), Fasta (LIHI01.1.fsa_nt.gz), and ASN.1 (LIHI01.1.bbs.gz) formats, can be accessed and/or reused at http://www.ncbi.nlm.nih.gov/nuccore/LIHI00000000.

In silico analyses

CLC genomics workbench 8.0 (QIAGEN, Aarhus), PathogenFinder 1.1, ResFinder 2.1, PlasmidFinder 1.3, and MLST 1.8 (MultiLocus Sequence Typing) were used to estimate the number of pathogenicity determinants, acquired antibiotic resistance genes, plasmids and the MLST using the S. Enteritidis genome (Larsen et al., 2012; Zankari et al., 2012; Cosentino et al., 2013; Carattoli et al., 2014). SeqSero tool was used for Salmonella serotyping by whole genome sequencing (Zhang et al., 2015).

PHAST search web tool was applied to detect, identify and annotate prophage sequences (Zhou et al., 2011). ISsaga was used for the high throughput identification and semiautomatic annotation of insertion sequences in the genome (Varani et al., 2011). The presence of molecular determinants of antimicrobial resistance was predicted based on homology and SNP models using the Comprehensive Antibiotic Resistance Database (CARD; https://card.mcmaster.ca/analyze/rgi), through Resistance Gene Identifier software (RGI; McArthur et al., 2013).

Results

LV60 isolate was serotyped as S. Enteritidis, using the method of Kauffman-White scheme, and found to be wild-type to all the antibiotics tested, except tetracycline.

The de novo assembly yielded 4.977 Mbp distributed in 83 contigs (largest contig with 970,921 bp) with a N50 of 491,005 bp. Overall, the structural and functional annotation with PGP detected 97 tRNA genes, 7 rRNA genes and identified 4656 mRNA genes.

From mapping against the bacterial genome of S. Enteritidis strain p125109, the main difference between the two genomes was the absence of the O-antigen polymerase gene wzy in the LV60 isolate, which in S. Enteritidis is located outside the O antigen gene cluster (Liu et al., 2014). The coding sequence of wzy gene was searched against the assembled genome using blastn, confirming its absence. The flanking regions of wzy gene, which coded for a disrupted membrane and a hypothetical protein, were also absent. The wzy gene is involved in the Wzx/Wzy-dependent pathway, which constitutes the predominant pathway for O-antigen production in Gram-negative bacteria, specifically in Salmonella (Hong et al., 2015).

However, in this study, the absence of the wzy gene did not compromised the use of a high-throughput genome sequencing serotype determination method (Zhang et al., 2015), which corroborated the result obtained by the gold standard method. Indeed, this method, based on the detection of O and H antigens encoding genes, predicted an antigenic profile 9:g,m:- based on the O-9,46 wbaV gene, which encodes to the O-antigen tyvelosyl transferase. Furthermore, the S. Enteritidis serotype was confirmed by the presence of sdf gene (Salmonella difference fragment virulence gene), a characteristic marker of commonly circulating S. enterica serovar Enteritidis (Agron et al., 2001).

Sixty-one SNVs were detected between LV60 and the S. Enteritidis strain p125109. The SNVs that resulted in amino acid substitutions are represented in Table 1. In silico analysis with ResFinder tool did not reveal the presence of any acquired antibiotic resistance genes (90% identity and 40% minimum length) or plasmids (95% identity). However, the RGI analysis, using the perfect algorithm, showed the presence of a Salmonella-specific MerR-like gold (Au) sensor- GolS—involved in Au resistance (Pontel et al., 2007). This constitutes a matter of concern since antibacterial biocides and metals can contribute to the development and maintenance of antibiotic resistance in bacterial communities through mechanisms of cross- or co-resistance (Baker-Austin et al., 2006; Lemire et al., 2013; Pal et al., 2015).

Furthermore, the RGI strict algorithm, which detects previously unknown variants of known antimicrobial resistance genes, identified 52 genes involved in efflux, transport, and permeability, which might justify the low-level tetracycline resistance identified by phenotypic methods (Table 2). Resistance to additional classes of antibiotics such as fluoroquinolones, aminoglycosides, and chloramphenicol were bioinformatically predicted. Indeed, efflux pumps are often associated with discrete decreases in antibiotic susceptibility that may not necessarily reflect an alteration in interpretation categories (Fernández and Hancock, 2012). Genes responsible for the intrinsic resistance to benzylpenicillin, glycopeptides, macrolides, and rifampicin were also detected.

The total number of pathogenicity determinants present in the genome of S. Enteritidis LV60, matching 1164 pathogenic families, showed a 94.1% certainty of the isolate being a human pathogen. Here we highlight the presence of Salmonella Pathogenicity Island 4, which usually encodes a non-fimbrial adhesion and the cognate type 1 secretion system (Gerlach et al., 2007).

The use of complementary web tools assigned this isolate to ST11, which according with MLST data (http://mlst.warwick.ac.uk/) is commonly found among CTX-M-14 and CTX-M-15-producing S. Enteritidis human isolates (Kim et al., 2011; Bado et al., 2012). In this study, the identification of ST11 in an isolate of animal origin, together with other pathogenicity determinants may suggest its zoonotic potential.

We also identified 6 prophage regions, among which three were incomplete and three were intact. The last included prophage regions reaching the lengths of 64.3, 49.2, and 31.7 Kb, and encoding 42, 78, and 66 DNA coding sequences, respectively.

Overall, 33 different IS were detected within the genome, which were distributed as follows: 27.03% of IS3 family, 18.92% of IS256 family, 13.51% of IS unclassified elements, 10.81% of IS200/IS605 complex, and of ISL3 family, 8.11% of IS481 family, 5.41% of IS630 family, and 2.7% of IS1 and IS110 families. All identified structures (pathogenicity island, prophages, ISs) constitute a multiplicity of pathogenicity factors in LV60 S. Enteritidis isolate and contribute for the fitness of the isolate in different environments; its presence may also suggest the possibility of acquisition of other factors by different mechanisms, including resistance genes e.g., by horizontal gene transfer, contributing to its biological diversity and genetic evolution.

Conclusion

The detection of an avian S. Enteritidis isolate harboring multiple efflux pumps, pathogenicity factors, a variety of mobile genetic elements and heavy-metal-tolerance genes raises concerns regarding the dissemination of infection in birds and potential risk of zoonotic transmission.

This study demonstrated the added value of WGS as a routine tool for surveillance programs directed to food-producing animals, which might complement sanitary measures, essential to prevent the spread of Salmonella infections among animals. It also proved to have an added value as a complementary typing method. Moreover, the simultaneous detection of putative Au resistance, intrinsic antibiotic resistant genes, and mobile genetic elements, underline this method as a helpful resource to follow the spread and evolution of antibiotic resistance in this species by genomic comparison studies.

Data access

This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession LIHI00000000. The version described in this paper is version LIHI01000000.

Funding

DJ has received research funding from Fundação para a Ciência e a Tecnologia (FCT, grant number SFRH/BD/80001/2011). VM was supported by FCT fellowship (grant SFRH/BPD/77486/2011), financed by the European Social Funds (COMPETE-FEDER) and national funds of the Portuguese Ministry of Education and Science (POPH-QREN). We thank the support of FCT grant number PEst-OE/AGR/UI0211/2011-2014 and UID/MULTI/00211/2013.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

• 1

  AgronP. G.WalkerR. L.KindeH.SawyerS. J.HayesD. C.WollardJ.et al. (2001). Identification by subtractive hybridization of sequences specific for Salmonella enterica serovar Enteritidis. Appl. Environ. Microbiol.67, 4984–4991. 10.1128/AEM.67.11.4984-4991.2001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  BadoI.García-FulgueirasV.CordeiroN. F.BetancorL.CaiataL.SeijaV.et al. (2012). First human isolate of Salmonella enterica serotype Enteritidis harboring bla_CTX−M−14 in South America. Antimicrob. Agents Chemother.56, 2132–2134. 10.1128/AAC.05530-11

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  Baker-AustinC.WrightM. S.StepanauskasR.McArthurJ. V. (2006). Co-selection of antibiotic and metal resistance. Trends Microbiol.14, 176–182. 10.1016/j.tim.2006.02.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BoisvertS.LavioletteF.CorbeilJ. (2010). Ray: simultaneous assembly of reads from a mix of high-throughput sequencing technologies. J. Comput. Biol.17, 1519–1533. 10.1089/cmb.2009.0238

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  CarattoliA.ZankariE.García-FernándezA.Voldby LarsenM.LundO.VillaL.et al. (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother.58, 3895–3903. 10.1128/AAC.02412-14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  ClementeL.ManageiroV.FerreiraE.Jones-DiasD.CorreiaI.ThemudoP.et al. (2013). Occurrence of extended-spectrum β-lactamases among isolates of Salmonella enterica subsp. enterica from food-producing animals and food products, in Portugal. Int. J. Food Microbiol.167, 221–228. 10.1016/j.ijfoodmicro.2013.08.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  ClementeL.ManageiroV.Jones-DiasD.CorreiaI.ThemudoP.AlbuquerqueT.et al. (2015). Antimicrobial susceptibility and oxymino-β-lactam resistance mechanisms in Salmonella enterica and Escherichia coli isolates from different animal sources. Res. Microbiol.166, 574–583. 10.1016/j.resmic.2015.05.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  CosentinoS.Voldby LarsenM.Møller AarestrupF.LundO. (2013). PathogenFinder-distinguishing friend from foe using bacterial whole genome sequence data. PLoS ONE8:e77302. 10.1371/journal.pone.0077302

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  DelmarJ. A.SuC. C.YuE. W. (2014). Bacterial multidrug efflux transporters. Annu. Rev. Biophys.43, 93–117. 10.1146/annurev-biophys-051013-022855

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  EgasC.BarrosoC.FroufeH. J.PachecoJ.AlbuquerqueL.Da CostaM. S. (2014). Complete genome sequence of the radiation-resistant bacterium Rubrobacter radiotolerans RSPS-4. Stand. Genomic Sci.9, 1062–1075. 10.4056/sigs.5661021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  European Food Safety Authority/European Center Disease Control (EFSA/ECDC) (2015). EU Summary Report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2013. EFSA J. 13:4036. 10.2903/j.efsa.2015.4036

  • CrossRef
  • Google Scholar

• 12

  FernándezL.HancockR. E. W. (2012). Adaptive and mutational resistance: role of porins and efflux pumps in drug resistance. Clin. Microbiol. Rev.25, 661–681. 10.1128/CMR.00043-12

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  FoleyS. L.JohnsonT. J.RickeS. C.NayakR.DanzeisenJ. (2013). Salmonella pathogenicity and host adaptation in chicken-associated serovars. Microbiol. Mol. Biol. Rev.77, 582–607. 10.1128/MMBR.00015-13

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  FoleyS. L.LynneA. M. (2008). Food animal-associated Salmonella challenges: pathogenicity and antimicrobial resistance. J. Anim. Sci.86, E173–E187. 10.2527/jas.2007-0447

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  FoleyS. L.LynneA. M.NayakR. (2008). Salmonella challenges: prevalence in swine and poultry and potential pathogenicity of such isolates. J. Anim. Sci.86, E149–E162. 10.2527/jas.2007-0464

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  GerlachR. G.JäckelD.StecherB.WagnerC.LupasA.HardtW. D.et al. (2007). Salmonella Pathogenicity Island 4 encodes a giant non-fimbrial adhesin and the cognate type 1 secretion system. Cell Microbiol.9, 1834–1850. 10.1111/j.1462-5822.2007.00919.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  GrimontP. A.WeillF. X. (2007). Antigenic Formulae of the Salmonella Serovars, 9th Edn. Paris: Institute Pasteur; WHO Collaborating Centre for Reference and Research on Salmonella.

  • Google Scholar

• 18

  HenselM. (2004). Evolution of pathogenicity islands of Salmonella enterica. Int. J. Med. Microbiol.294, 95–102. 10.1016/j.ijmm.2004.06.025

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  HongY.MorcillaV. A.LiuM. A.RussellE. L.ReevesP. R. (2015). Three Wzy polymerases are specific for particular forms of an internal linkage in otherwise identical O units. Microbiology161, 1639–1647. 10.1099/mic.0.000113

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  KimY.BaeI. K.JeongS. H.LeeC. H.LeeH. K.AhnJ.et al. (2011). Occurrence of IncFII plasmids carrying the bla_CTX−M−15 gene in Salmonella enterica serovar Enteritidis sequence type 11 in Korea. Diagn. Microbiol. Infect. Dis.71, 171–173. 10.1016/j.diagmicrobio.2011.05.004

  • CrossRef
  • Google Scholar

• 21

  LagesenK.HallinP. F.RødlandE.StærfeldtH. H.RognesT.UsseryD. W. (2007). RNammer: consistent annotation of rRNA genes in genomic sequences. Nucleic Acids Res.35, 3100–3108. 10.1093/nar/gkm160

  • CrossRef
  • Google Scholar

• 22

  LarsenM. V.CosentinoS.RasmussenS.FriisC.HasmanH.MarvigR. L.et al. (2012). Multilocus sequence typing of total-genome-sequenced bacteria. J. Clin. Microbiol.50, 1355–1361. 10.1128/JCM.06094-11

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  LemireJ. A.HarrisonJ. J.TurnerR. J. (2013). Antimicrobial activity of metals: mechanisms, molecular targets and applications. Nat. Rev. Microbiol.11, 371–384. 10.1038/nrmicro3028

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  LiX.-Z.PlésiatP.NikaidoH. (2015). The challenge of efflux-mediated antibiotic resistance in gram-negative bacteria. Clin. Microbiol. Rev.28, 337–418. 10.1128/CMR.00117-14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  LiuB.KnirelY. A.FengL.PerepelovA. V.SenchenkovaS. N.ReevesP. R.et al. (2014). Structural diversity in Salmonella O antigens and its genetic basis. FEMS Microbiol. Rev.38, 56–89. 10.1111/1574-6976.12034

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  Lutful KabirS. M. (2010). Avian colibacillosis and salmonellosis: a closer look at epidemiology, pathogenesis, diagnosis, control and public health concerns. Int. J. Environ. Res. Public Health7, 89–114. 10.3390/ijerph7010089

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  McArthurA. G.WaglechnerN.NizamF.YanA.AzadM. A.BaylayA. J.et al. (2013). The comprehensive antibiotic resistance database. Antimicrob. Agents Chemother.57, 3348–3357. 10.1128/AAC.00419-13

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  MedardusJ. J.MollaB. Z.NicolM.MorrowW. M.Rajala-SchultzP. J.KazwalaR.et al. (2014). In-feed use of heavy metal micronutrients in U.S. Swine production systems and its role in persistence of multidrug-resistant salmonellae. Appl. Environ. Microbiol.80, 2317–2325. 10.1128/AEM.04283-13

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  PalC.Bengtsson-PalmeJ.KristianssonE.LarssonD. G. (2015). Co-occurrence of resistance genes to antibiotics, biocides and metals reveals novel insights into their co-selection potential. BMC Genomics16:964. 10.1186/s12864-015-2153-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  PontelL. B.AuderoM. E. P.EsparizM.ChecaS. K.SonciniF. C. (2007). GolS controls the response to gold by the hierarchical induction of Salmonella-specific genes that include a CBA efflux-coding operon. Mol. Microbiol.66, 814–825. 10.1111/j.1365-2958.2007.05963.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  PooleK. (2004). Efflux-mediated multiresistance in Gram-negative bacteria. Clin. Microbiol. Infect.10, 12–26. 10.1111/j.1469-0691.2004.00763.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  VaraniA. M.SiguierP.GourbeyreE.CharneauV.ChandlerM. (2011). ISsaga is an ensemble of web-based methods for high throughput identification and semi-automatic annotation of insertion sequences in prokaryotic genomes. Genome Biol.12, R30. 10.1186/gb-2011-12-3-r30

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  ZankariE.HasmanH.CosentinoS.VestergaardM.RasmussenS.LundO.et al. (2012). Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother.67, 2640–2644. 10.1093/jac/dks261

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  ZhangS.YinY.JonesM. B.ZhangZ.Deatherage KaiserB. L.DinsmoreB. A.et al. (2015). Salmonella serotype determination utilizing high-throughput genome sequencing data. J. Clin. Microbiol.53, 1685–1692. 10.1128/JCM.00323-15

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  ZhouY.LiangY.LynchK. H.DennisJ. J.WishartD. S. (2011). PHAST: a fast phage search tool. Nucleic Acids Res.39, W347–W352. 10.1093/nar/gkr485

  • Pubmed Abstract
  • CrossRef
  • Google Scholar