- Research article
- Open Access
Identification of novel antimicrobial resistance genes from microbiota on retail spinach
BMC Microbiologyvolume 13, Article number: 272 (2013)
Drug resistance genes and their mobile genetic elements are frequently identified from environmental saprophytic organisms. It is widely accepted that the use of antibiotics in animal husbandry selects for drug resistant microorganisms, which are then spread from the farm environment to humans through the consumption of contaminated food products. We wished to identify novel drug resistance genes from microbial communities on retail food products. Here, we chose to study the microbial communities on retail spinach because it is commonly eaten raw and has previously been associated with outbreaks of bacterial infections.
We created metagenomic plasmid libraries from microbiota isolated from retail spinach samples. We identified five unique plasmids that increased resistance to antimicrobial drugs in the E. coli host. These plasmids were identified in E. coli that grew on plates that contained ampicillin (pAMP), aztreonam (pAZT), ciprofloxacin (pCIP), trimethoprim (pTRM), and trimethoprim-sulfamethoxazole (pSXT). We identified open reading frames with similarity to known classes of drug resistance genes in the DNA inserts of all 5 plasmids. These drug resistance genes conferred resistance to fluoroquinolones, cephalosporins, and trimethoprim, which are classes of antimicrobial drugs frequently used to treat human Gram negative bacterial infections. These results show that novel drug resistance genes are found in microbiota on retail produce items.
Here we show that microbiota of retail spinach contains DNA sequences previously unidentified as conferring antibiotic resistance. Many of these novel sequences show similarity to genes found in species of bacteria, which have previously been identified as commensal or saprophytic bacteria found on plants. We showed that these resistance genes are capable of conferring clinically relevant levels of resistance to antimicrobial agents. Food saprophytes may serve as an important reservoir for new drug-resistance determinants in human pathogens.
The spread of antimicrobial resistance genes has made previously manageable bacterial infections increasingly more difficult to treat. In addition, there has been a gradual decline in the development of new antimicrobial drugs, especially against Gram negative bacterial pathogens. The identification of genes in Gram negative bacteria that confer resistance to cephalosporins, carbapenems, and fluoroquinolones has created fears that we are returning to the pre-antibiotic era . These multidrug-resistant infections often occur in hospitals and are frequently caused by species belonging to the normal microbiota of the human host [1, 2]. This suggests that the microbiota of the patients themselves is the reservoir for many of the organisms that cause hospital acquired infections. Furthermore, recent work has demonstrated that the intestinal microbiota of humans and food animals are a reservoir of drug resistance genes . Consequently, a better understanding of how drug resistance genes enter the human microbiota is imperative to better prevent drug resistant infections.
Drug resistance genes and their mobile genetic elements are frequently identified from environmental saprophytic organisms. These include samples taken from soil, water, and wild animals [4–6]. Additionally, these genes have been identified in environmental samples from ancient and pristine environments – samples that have never been exposed to human activity [4, 5, 7]. Due to the great diversity of antibiotic resistance genes found in environment, it has been hypothesized that environmental microbes serve as a reservoir of drug resistance genes and that a few then enter human pathogens [4, 5]. These drug resistance genes are spread between bacteria via mobile genetic elements, such as plasmids, transposons, and integrons [4–6, 8]. The detection of mobile genetic elements and drug resistance genes in the environment has led to numerous studies and policies to address the effects of environmental exposure to antimicrobial agents on human pathogens [9, 10].
It is widely accepted that the use of antibiotics in animal husbandry selects for drug resistant microorganisms, which are then spread from the farm environment to humans through the consumption of contaminated food products . Numerous studies of bacterial pathogens in food products, such as Campylobacter and Salmonella, have demonstrated that the use of antimicrobial drugs in agriculture can result in drug resistant infections in humans [9–13]. However, the majority of studies have been limited to species of zoonotic pathogens that cause foodborne disease and these studies are frequently done as part of national surveillance programs for food safety. Species of bacteria that are not usually considered foodborne pathogens, but nonetheless are found in both the human and food product microbiota, are usually not included in studies of drug resistant bacteria in retail food products. Studies of microbiota of animals demonstrated that commensal organisms are a reservoir of antimicrobial drug resistance genes . These studies include the identification of antimicrobial resistance genes from animal feces including chickens and cows [3, 14–16]. Animal manure is frequently used as fertilizer in agriculture and may contribute to the spread of drug resistance genes. The spread of drug resistance genes by commensal bacteria on food products is an area that requires further study .
Produce items, which are frequently eaten raw, are one way consumers are exposed to microbiota on retail food products . Previous work has shown that the normal microbiota of retail produce items harbors clinically relevant drug resistance genes . However, previous studies have relied on PCR based methods to identify known drug resistance genes, which limits the number and types of drug resistance that could potentially be identified. Other studies used functional metagenomic libraries to identify novel antimicrobial resistance genes from environmental samples in a sequence independent manner [19, 20]. We wanted to apply this sequence independent approach to investigate the presence of antimicrobial resistance genes on retail spinach. We chose to study the microbial communities on retail spinach because it is commonly eaten raw and has previously been associated with outbreaks of bacterial infections .
To do this, we made two metagenomic plasmid libraries with DNA isolated from the microbiota of retail spinach. One library was made from a cultured sample of spinach microbiota while the other was made in a culture independent manner. We then screened these libraries for their ability to confer resistance to antibiotics to an E. coli host.
Results and discussion
Isolation of antibiotic resistant clones
The first plasmid library, which was constructed from a cultured sample, contained 160 Mb of inserted DNA. The second library, which was constructed from an uncultured sample, contained 140 Mb of inserted DNA. We first constructed a cultured library because we wished to enrich for microbial DNA to increase the chances of cloning DNA sequences that contained drug-resistance genes. We constructed a library from an uncultured sample because we wanted to identify potential drug-resistance genes from bacterial organisms that cannot be cultivated in artificial medium. From the cultivated library, we isolated four different antimicrobial resistance-conferring clones. From the uncultivated library, we isolated one additional antimicrobial resistance-conferring clone.
The mean size of the plasmid DNA inserts in both libraries was two Kb. We identified five unique plasmids that conferred increased drug resistance (minimum inhibitory concentration or MIC) to the host E. coli. Each plasmid was named after the antimicrobial agent to which it conferred resistance (Table 1). These plasmids were identified in E. coli that grew on plates that contained ampicillin (pAMP), aztreonam (pAZT), ciprofloxacin (pCIP), trimethoprim (pTRM), and trimethoprim-sulfamethoxazole (pSXT). The plasmids pAMP, pAZT, pCIP, and pTRM were isolated from the library made from a cultured sample. The plasmid pSXT was isolated from the library made from an uncultured sample.
pAMP increased the MIC of ampicillin 4 fold (4 ug/ml to 16 ug/ml). pAZT increased the MIC of the host strain 96 fold to aztreonam (.125 ug/ml to 12 ug/ml), 10 fold to cefepime (.096 ug/ml to 1 ug/ml), and 6 fold to piperacillin (2 ug/ml to 12 ug/ml). Additionally, pAZT encoded an ESBL phenotype as measured by the ceftazidime, ceftazidime/clavulanic acid ESBL Etest, (TZ 16 ug/ml and TZL 1 ug/ml). pCIP caused a 62-fold increase in resistance to ciprofloxacin (<.002 ug/ml to .125 ug/ml) as well as a 31-fold increase in resistance to levofloxacin (.012 ug/ml to .38 ug/ml). The MIC of trimethoprim for E. coli carrying pTRM increased >258 fold (.124 ug/ml to >32 ug/ml) and 7 fold to trimethoprim-sulfamethoxazole (.064 ug/ml to .5 ug/ml). pSXT caused a 15-fold increase in resistance to trimethoprim-sulfamethoxazole (.064 ug/ml to 1 ug/ml) as well as an >256-fold increase in resistance to trimethoprim alone (.125 ug/ml to >32 ug/ml).
Identification of antibiotic resistance genes and phylogenetic analysis
In pAZT we identified a sequence with 94% identity at the nucleotide level to a gene that encodes penicillin-binding protein 1A identified in Bacillus subtilis subsp. Spizizenii [GenBank, gb|CP002905.1|]. The expression of altered penicillin binding proteins are known to confer resistance to beta-lactam and cephalosporin antibiotics in various clinically important pathogens . However, the ability of this sequence to confer clinically relevant levels of cephalosporin resistance or an ESBL phenotype has not been previously reported.
In pAMP we identified a beta-lactamase gene with 71% identity to the ERP-1 gene that encodes a class A extended spectrum beta-lactamase found in Erwinia persicin [GenBank, gb|AY077733.1|] . When transformed into an E. coli host, ERP-1 was reported to increase resistance to penicillins and cephalosporins, including piperacillin, cefotaxime, and ceftazidime . Surprisingly, we found that pAMP did not increase the MIC of the host E. coli to piperacillin, ceftazidime, or cefotaxime. Also, pAMP did not increase the MIC of this host strain to any tested cephalosporin, monobactam, or carbapenem (cefepime, aztreonam, or imipenem). However, this isolate tested positive for the presence of a beta-lactamase by the nitrocefin assay. These results suggest that the novel sequence we identified in pAMP is distinct from ERP-1 in terms of the spectrum of drug-resistance phenotype it encodes. However, it is also possible that these observations are artifacts due to poor expression of the gene in a heterologous host.
Surprisingly, when the pCIP DNA sequence was submitted to BlastN, only two other sequences in the NCBI non-redundant nucleotide database were identified. The sequences were part of whole genome sequences of Exiguobacterium antarcticum and Exiguobacterium sibiricum. These species were identified in a frozen Antarctic lake and a core sample of the Siberian tundra [24, 25]. When submitted to BlastP, a fluoroquinolone resistance protein from Oceanobacillus sp. Ndiop was identified. This quinolone resistance protein was a predicted pentapeptide repeat protein (PRP) [GenBank, ref|ZP_10910075.1|]. One known class of plasmid-mediated quinolone resistance conferring sequences is called QNR. QNRs are pentapeptide repeat proteins and have been associated with extended spectrum beta-lactamases . In addition to the PRP, the pCIP DNA insert contained a second open reading frame that showed similarity to a beta-lactamase domain containing protein (Figure 1). However, pCIP did not increase the MIC of the host E. coli to any of the tested beta-lactam antibiotics (ampicillin and piperacillin), cephalosporin antibiotics (ceftazidime, cefotaxime, cefepime), or carbapenem antibiotics (imipenem). Additionally, pCIP did not test positive for beta-lactamase production by the nitrocefin assay. This indicated that the predicted beta-lactamase is either not expressed or does not function as a beta-lactamase.
We subcloned the PRP sequence in pCIP (pPRP:1B) in order to confirm that the predicted PRP was able to confer resistance to ciprofloxacin. As compared to the empty vector, the pPRP:1B increased the MIC of ciprofloxacin of the host E. coli 7 fold (.016 to .125) and the MIC of levofloxacin 5 fold (.047 to .25 ug/ml). The MIC conferred by pPRP:1B to the host E. coli is consistent with previously reported MICs from other QNR sequences found in human pathogens .
We found that the pTRM DNA insert has 84% identity to a region of the Erwinia pyrifoliae DSM 12163 complete genome [GenBank, emb|FN392235.1]. This region contains the folA gene, which encodes a dihydrofolate reductase (DHFR). DHFR is the target of trimethoprim . The expression of a DHFR that is not susceptible to trimethoprim is a well-known mechanism of resistance . However, the acquisition of the folA gene from Erwinia pyrifoliae DSM 12163 has not previously been shown to confer resistance to trimethoprim.
Similarly, the DNA insert from pSXT has 88% identity to a region of the Pseudomonas fluorescens SBW25 complete genome . This region also encodes a predicted dihydrofolate reductase. This sequence from Pseudomonas fluorescens SBW25 has not previously been shown to confer resistance to trimethoprim or trimethoprim sulfamethoxazole.
Here we show that microbiota of retail spinach contains previously unidentified antibiotic resistance-conferring genes and that functional metagenomic libraries can be used to screen retail food products for drug resistance genes in a sequence independent manner. Furthermore, due to the limited amount of DNA that can be cloned into a plasmid library and the requirement that the drug resistance gene be expressed in a heterologous host, it is likely we only identified a fraction of drug resistance conferring genes present in our spinach samples.
Although none was identical in DNA sequence, many of these novel sequences show sequence similarity to genes found in species of bacteria that have previously been identified as commensal or saprophytic bacteria found on plants [23, 28]. This suggests the sequences we identified are not the result of contamination from animals or humans. We showed that these resistance genes are capable of conferring clinically relevant levels of resistance to commonly used classes of antimicrobial agents, including cephalosporins and fluoroquinolones.
The novel antimicrobial resistance genes we identified include beta-lactamases, a pentapeptide repeat protein, a penicillin binding protein, and putative dihydrofolate reductase genes. These types of resistance mechanisms are some of the most common and clinically problematic mechanisms of drug resistance found in pathogens . We do not know at this time if these genes will become clinically important, and one limitation of this study is that we did not analyze these genes for their potential for horizontal transfer to human pathogens. Recent functional genomics analysis of environmental soil samples has not only identified drug-resistance genes with identical nucleotide sequences from human pathogens, but also mobile gene sequences providing evidence for possible horizontal gene transfers . Further studies using sequence independent methods to identify antimicrobial resistance genes from retail food products should be done to better understand the role of saprophytes as a reservoir for new drug-resistance genes.
Metagenome plasmid library construction and screening
Two metagenomic plasmid libraries of spinach microbiota were constructed. One was based on DNA extracted from cultured bacteria and the other was based on DNA extracted from uncultured spinach wash. The metagenomic DNA used to create the first library was obtained by washing twenty five grams of bagged “baby spinach” in PBS. A description of the spinach used to create the library has been previously published . Briefly, the spinach was purchased from a supermarket located in Berkeley, California in 2007. They included organic as well as non-organic spinach. One milliliter of the PBS wash was then used to inoculate 50 ml TSB. This culture was grown at 37°C with shaking overnight. The culture was then centrifuged at 10,000 x g for 10 minutes. DNA was extracted from the resulting pellet by the phenol chloroform method.
The metagenomic DNA used to create the second library was obtained by washing six bags of “baby spinach” in two liters of PBS. Six different brands of spinach were purchased from three retailers located in Berkeley, California in 2011. The spinach was incubated in PBS at room temperature for two hours. The resulting wash was then filtered through sterilized cheesecloth and a sterilized coffee filter to remove spinach debris. The filtered wash was then centrifuged at 10,000 x g for 20 minutes. DNA was extracted directly from the resulting pellet with the Gnome DNA isolation Kit, MP Biomedical.
Two plasmid libraries with metagenomic DNA inserts were constructed in pSMART-LC kan vector in the E. coli host, E. Cloni (Lucigen corp., Middleton WI). The pSMART vector confers resistance to kanamycin and has transcriptional terminators flanking the cloning sites. Consequently, transcription of the cloned sequences requires a native promoter.
E. coli clones containing the two plasmid libraries were then screened for resistance to antimicrobial agents on Mueller Hinton agar plates containing one of the following 16 antimicrobial agents: ampicillin, carbenicillin, ticarcillin, amoxicillin/clavulanic acid, ticarcillin/clavulanic acid, cefotaxime, ceftazidime, aztreonam, meropenem, gentamicin, nalidixic acid, ciprofloxacin, trimethoprim, trimethoprim-sulfamethoxazole, chloramphenicol, or tetracycline. The phenotype of resistance to antimicrobial agents was confirmed by retransforming the recombinant plasmid into E. Cloni. Growth on Muller Hinton agar containing kanamycin was used as a positive control for transformation. The acquisition of drug resistance from the transformation of the plasmid was demonstrated by growth on Mueller Hinton agar containing the corresponding antimicrobial agent. We used the empty vector, pSMART, in the E. coli host as a negative control for antibiotic stability.
Antimicrobial susceptibility testing
The MIC of each E. coli clone was determined by Etest (Biomerieux, France) according to manufactures recommendations. All Etests were repeated in at least two independent experiments. ATCC 29522, ATCC 700603, and ATCC 35218 were used for control as recommended by the manufacturer. In accordance with manufactures recommendations, less than a fourfold difference in MIC was considered to be with in the expected margin of error for this test.
A colony was spotted onto a sterile Petri dish and then covered with nitrocefin, as previously described . E. Cloni containing the empty vector (pSMART) was used as a negative control.
Sequencing and data analysis: The sequence of the DNA insert in the resistance conferring plasmid was determined by primer walking at the University of California, Berkeley sequencing facility. The sequences were assembled with Geneious Version 5.6, (Biomatters, New Zealand). The insert sequences were then submitted to ORFinder and the BLAST suit of programs at NCBI . The nucleotide sequences of the insert from each plasmid have been deposited in Genbank with the following accession numbers: pAMP: KF791056, pAZT: KF791057, pCIP: KF791058, pTRM: KF791059, pSXT: KF791060.
Standard protocols for ligation independent cloning in to vector 1B, QB3 Macrolab, University of California, Berkeley were used.
Boucher HW, Talbot GH, Bradley JS, Edwards JE, Gilbert D, Rice LB, Scheld M, Spellberg B, Bartlett J: Bad bugs, No drugs: No ESKAPE! an update from the infectious diseases society of America. Clin Infect Dis. 2009, 48 (1): 1-12. 10.1086/595011.
Choffnes ER, Relman DA, Mack A: Antibiotic resistance: Implications for global health and novel intervention strategies: Workshop summary. 2010, Washington D.C: National Academies Press
Penders J, Stobberingh EE, Savelkoul PHM, Wolffs P: The human microbiome as a reservoir of antimicrobial resistance. Front Microbiol. 2013, 4: 86-
Wright GD: Antibiotic resistance in the environment: a link to the clinic?. Curr Opin Microbiol. 2010, 13 (5): 589-594. 10.1016/j.mib.2010.08.005.
Allen HK, Donato J, Wang HH, Cloud-Hansen KA, Davies J, Handelsman J: Call of the wild: antibiotic resistance genes in natural environments. Nat Rev Microbiol. 2010, 8 (4): 251-259. 10.1038/nrmicro2312.
Martinez JL: Bottlenecks in the transmission of antibiotic resistance from natural ecosystems to human bacterial pathogens. Front Microbiol. 2011, 2: 265-
D’Costa VM, King CE, Kalan L, Morar M, Sung WWL, Schwarz C, Froese D, Zazula G, Calmels F, Debruyne R, et al: Antibiotic resistance is ancient. Nature. 2011, 477 (7365): 457-461. 10.1038/nature10388.
Levy SB, Marshall B: Antibacterial resistance worldwide: causes, challenges and responses. Nat Med. 2004, 10: S122-S129. 10.1038/nm1145.
Expert Meeting on Animal Feed Impact on Food, Safety: ReportAnimal feed impact on food safety: report of the FAO/ WHO Expert Meeting: 8-12 October 2007. 2008, Rome, Italy: FAO and WHO
Hawser S: Surveillance programmes and antibiotic resistance: worldwide and regional monitoring of antibiotic resistance trends. Antibiotic Resistance Volume 211. Edited by: Coates A. 2012, Berlin Heidelberg: Springer, 31-43.
Hawkey PM: The growing burden of antimicrobial resistance. J Antimicrob Chemother. 2008, 62 (suppl 1): i1-i9.
Angulo FJ, Baker NL, Olsen SJ, Anderson A, Barrett TJ: Antimicrobial use in agriculture: controlling the transfer of antimicrobial resistance to humans. Semin Pediatr Infect Dis. 2004, 15 (2): 78-85. 10.1053/j.spid.2004.01.010.
U.S. Department of Health and Human Services, Food and Drug Administration: Drug Administration. National antimicrobial resistance monitoring system–enteric bacteria (NARMS): 2008 executive report. 2011, Rockville
Zhou W, Wang Y, Lin J: Functional cloning and characterization of antibiotic resistance genes from the chicken Gut microbiome. Appl Environ Microbiol. 2012, 78 (8): 3028-3032. 10.1128/AEM.06920-11.
Durso LM, Harhay GP, Bono JL, Smith TPL: Virulence-associated and antibiotic resistance genes of microbial populations in cattle feces analyzed using a metagenomic approach. J Microbiol Methods. 2011, 84 (2): 278-282. 10.1016/j.mimet.2010.12.008.
Qu A, Brulc JM, Wilson MK, Law BF, Theoret JR, Joens LA, Konkel ME, Angly F, Dinsdale EA, Edwards RA, et al: Comparative metagenomics reveals host specific metavirulomes and horizontal gene transfer elements in the chicken cecum microbiome. PLoS ONE. 2008, 3 (8): e2945-10.1371/journal.pone.0002945.
Lynch MF, Tauxe RV, Hedberg CW: The growing burden of foodborne outbreaks due to contaminated fresh produce: risks and opportunities. Epidemiology & Infection. 2009, 137 (Special Issue 03): 307-315.
Raphael E, Wong LK, Riley LW: Extended-spectrum beta-lactamase gene sequences in gram-negative saprophytes on retail organic and nonorganic spinach. Appl Environ Microbiol. 2011, 77 (5): 1601-1607. 10.1128/AEM.02506-10.
Allen HK, Moe LA, Rodbumrer J, Gaarder A, Handelsman J: Functional metagenomics reveals diverse [beta]-lactamases in a remote Alaskan soil. ISME J. 2008, 3 (2): 243-251.
Donato JJ, Moe LA, Converse BJ, Smart KD, Berklein FC, McManus PS, Handelsman J: Metagenomic analysis of apple orchard soil reveals antibiotic resistance genes encoding predicted bifunctional proteins. Appl Environ Microbiol. 2010, 76 (13): 4396-4401. 10.1128/AEM.01763-09.
Wendel AM, Johnson DH, Sharapov U, Grant J, Archer JR, Monson T, Koschmann C, Davis JP: Multistate outbreak of Escherichia coli O157:H7 infection associated with consumption of packaged spinach, August–September 2006: The Wisconsin investigation. Clin Infect Dis. 2009, 48 (8): 1079-1086. 10.1086/597399.
Zhao S, Duncan M, Tomberg J, Davies C, Unemo M, Nicholas RA: Genetics of chromosomally mediated intermediate resistance to ceftriaxone and cefixime in Neisseria gonorrhoeae. Antimicrob Agents Chemother. 2009, 53 (9): 3744-3751. 10.1128/AAC.00304-09.
Vimont S, Poirel L, Naas T, Nordmann P: Identification of a chromosome-borne expanded-spectrum class A β-lactamase from Erwinia persicina. Antimicrob Agents Chemother. 2002, 46 (11): 3401-3405. 10.1128/AAC.46.11.3401-3405.2002.
Vishnivetskaya T, Kathariou S, McGrath J, Gilichinsky D, Tiedje JM: Low-temperature recovery strategies for the isolation of bacteria from ancient permafrost sediments. Extremophiles. 2000, 4 (3): 165-173. 10.1007/s007920070031.
Carneiro AR, Ramos RTJ, Dall’Agnol H, Pinto AC, de Castro SS, Santos AR, Guimarães LC, Almeida SS, Baraúna RA, das Graças DA, et al: Genome sequence of Exiguobacterium antarcticum B7, isolated from a biofilm in Ginger Lake, King George Island, Antarctica. J Bacteriol. 2012, 194 (23): 6689-6690. 10.1128/JB.01791-12.
Nordmann P, Poirel L: Emergence of plasmid-mediated resistance to quinolones in Enterobacteriaceae. J Antimicrob Chemother. 2005, 56 (3): 463-469. 10.1093/jac/dki245.
Antimicrobial Agents: Antibacterials and Antifungals. Edited by: Bryskier A. 2005, Washington: ASM Press
Silby M, Cerdeno-Tarraga A, Vernikos G, Giddens S, Jackson R, Preston G, Zhang X-X, Moon C, Gehrig S, Godfrey S, et al: Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens. Genome Biol. 2009, 10 (5): R51-10.1186/gb-2009-10-5-r51.
Forsberg KJ, Reyes A, Wang B, Selleck EM, Sommer MO, Dantas G: The shared antibiotic resistome of soil bacteria and human pathogens. Science. 2012, 337 (6098): 1107-1111. 10.1126/science.1220761.
Bourgault AM, Rosenblatt JE: Characterization of anaerobic gram-negative bacilli by using rapid slide tests for beta-lactamase production. J Clin Microbiol. 1979, 9 (6): 654-656.
Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden T: NCBI BLAST: a better web interface. Nucleic Acids Res. 2008, 36 (Web Server issue): W5-W9.
We would like to acknowledge the invaluable contributions of Eva Raphael and Olivera Marjanovic to the coordination of this study and the laboratory work required to create the first metagenomic library. This project was supported in part by NIAID/NIH grant R01AI059523 and T32AI007620.
The authors declare that they have no competing interests.
HB carried out the laboratory and bioinformatics studies and drafted the manuscript. LWR conceived of the study, and participated in its design and coordination and helped to draft the manuscript. Both authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.