Next Article in Journal
Presence of Methicillin Resistant Staphylococcus aureus (MRSA) in Feces of the Small Indian Mongooses (Urva auropunctata) on Saint Kitts and Nevis, West Indies
Next Article in Special Issue
Short Tryptamine-Based Peptoids as Potential Therapeutics for Microbial Keratitis: Structure-Function Correlation Studies
Previous Article in Journal
Occurrence of Extended Spectrum Cephalosporin-, Carbapenem- and Colistin-Resistant Gram-Negative Bacteria in Fresh Vegetables, an Increasing Human Health Concern in Algeria
Previous Article in Special Issue
Antimicrobial Random Peptide Mixtures Eradicate Acinetobacter baumannii Biofilms and Inhibit Mouse Models of Infection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomewide Analysis and Biological Characterization of Cathelicidins with Potent Antimicrobial Activity and Low Cytotoxicity from Three Bat Species

by
Munjeong Choi
,
Hye-sun Cho
,
Byeongyong Ahn
,
Somasundaram Prathap
,
Soundrarajan Nagasundarapandian
and
Chankyu Park
*
Department of Stem Cell and Regenerative Biotechnology, Konkuk University, Hwayang-dong, Gwangjin-gu, Seoul 05029, Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Antibiotics 2022, 11(8), 989; https://doi.org/10.3390/antibiotics11080989
Submission received: 1 July 2022 / Revised: 19 July 2022 / Accepted: 20 July 2022 / Published: 22 July 2022
(This article belongs to the Special Issue Peptide-Based Antibiotics: Challenges and Opportunities)

Abstract

:
Cathelicidins are potent antimicrobial peptides with broad spectrum antimicrobial activity in many vertebrates and an important component of the innate immune system. However, our understanding of the genetic variations and biological characteristics of bat cathelicidins is limited. In this study, we performed genome-level analysis of the antimicrobial peptide cathelicidins from seven bat species in the six families, listed 19 cathelicidin-like sequences, and showed that the number of functional cathelicidin genes differed among bat species. Based on the identified biochemical characteristics of bat cathelicidins, three cathelicidins, HA-CATH (from Hipposideros armiger), ML-CATH (from Myotis lucifugus), and PD-CATH (from Phyllostomus discolor), with clear antimicrobial signatures were chemically synthesized and evaluated antimicrobial activity. HA-CATH showed narrow-spectrum antibacterial activity against a panel of 12 reference bacteria, comprising 6 Gram-negative and 6 Gram-positive strains. However, ML-CATH and PD-CATH showed potent antibacterial activity against a broad spectrum of Gram-negative and Gram-positive bacteria with minimum inhibitory concentration (MIC) of 1 and 3 μg/mL, respectively, against Staphylococcus aureus. ML-CATH and PD-CATH also showed antifungal activities against Candida albicans and Cryptococcus cuniculi with MIC of 5 to 40 μg/mL, respectively, and 80% inhibition of the metabolism of Mucor hiemalis hyphae at 80 μg/mL, while displaying minimal cytotoxicity to HaCaT cells. Taken together, although the spectrum and efficacy of bat cathelicidins were species-dependent, the antimicrobial activity of ML-CATH and PD-CATH was comparable to that of other highly active cathelicidins in vertebrates while having negligible cytotoxicity to mammalian cells. ML-CATH and PD-CATH can be exploited as promising candidates for the development of antimicrobial therapeutics.

1. Introduction

Chiroptera (bats) are the second largest order of mammals following rodents [1], and they possess unique biological features, such as laryngeal echolocation [2], vocal learning [3], the ability to fly [4], and exceptional longevity [5], that are the rarest among all mammals. Additionally, bats are known as reservoir hosts of many zoonotic viruses, including the Nipah virus [6], Ebola virus [7] and coronavirus [8]. With no observed pathological signs of the disease, they can better tolerate viral infections than most other mammals because of their adaptation to antiviral immune responses [9]. In particular, the similarity between severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2), the causative agent of current coronavirus disease 2019 outbreaks, and bat-borne coronavirus (Bat CoV RATG13) illuminated bats as special reservoirs for viral pathogens that cause severe diseases in other mammals [10]. Currently, various efforts to develop therapeutics against SARS-CoV-2 are being made [11,12].
Comparative analyses of the immune systems of bats and other mammals have been conducted, which suggest that bats have evolved novel mechanisms to limit proinflammatory responses induced by viral infections, while maintaining type I interferon (IFN) responses to restrict viral propagation [13,14,15]. Many severe viral infections cause excessive inflammation-associated pathogenesis in humans [16]. Studies have also shown the conservation of cytosolic RNA sensors, including pattern recognition receptors such as Toll-like receptors [14]. To better understand the genetic and evolutionary mechanisms underlying the unique adaptations of bats, studies on diverse bat genomes, including the Bat1K initiative, are being conducted [9,17,18]. However, the complex pathogenicity of viral infections in bats and its impact on bat fitness have not yet been clearly elucidated [19].
Antimicrobial peptides (AMPs), also known as host defense peptides, are important components of the innate immune system that are found in various organisms [20]. In addition to their broad range of antimicrobial activities against microbes, including bacteria, parasites, fungi, and virus [21,22,23,24,25], they act as immunomodulatory molecules [26,27,28]. Cathelicidins are potent AMPs consisting of an evolutionarily conserved cathelin-like domain (CLD) and an active mature peptide region at the C-terminus [29]. The release of mature peptides via enzymatic cleavage by specific proteases results in the generation of functional AMPs [30,31]. The expression levels of AMPs, including cathelicidins, are associated with the anti-inflammatory responses and clearance of pathogens [26,32,33].
Recent studies have revealed various functions other than antibacterial activity of vertebrate cathelicidins, particularly for LL-37 and CRAMP of humans and mice, respectively [34,35]. For example, a defective intestinal expression of CRAMP in non-obese diabetic (NOD) mice induced the colonic dysbiosis and could contribute to autoimmunity in the pancreas in turn [36]. Moreover, endogenous CRAMP and LL-37 are suggested as inactivators of Zika virus by inhibiting viral replication and directly inactivating viral virions, thus exhibiting therapeutic potential against the viral infection in vivo [37]. The protective role of LL-37 against Candida auris known for causing severe infectious disease with high mortality attributed to multidrug-resistance against all major classes of antifungal drugs was also reported [38]. In addition, in vivo free-radical scavenging ability of a cathelicidin was reported from frogs [39].
The availability of genome sequences from diverse vertebrates provides a window of opportunity to identify genetic and functional diversity of cathelicidins [40,41,42]. However, understanding the biological effects of newly discovered cathelicidins from the high throughput analysis is slowly progressing. In this study, we conducted genome-wide in silico analysis of cathelicidin genes from seven bat species in the six families, Rousettus aegyptiacus, Hipposideros armiger, Miniopterus natalensis, Myotis lucifugus, Rhinolophus ferrumequinum, Phyllostomus discolor, and Desmodus rotundus, to maximize the genetic diversity and experimentally characterize their biological effects, including antibacterial, antifungal, and mammalian cytotoxic effects, in addition to the recent bioinformatics analysis of bat cathelicidins [43]. Our results contribute to a better understanding of the possible roles of AMPs in the antimicrobial defense of bats and can be exploited for the development of antimicrobial therapeutics.

2. Results

2.1. Identification of 27 AMP Family Genes via the In Silico Analysis of R. ferrumequinum Genome

To understand the genetic characteristics of the AMP subgenome in bats, we conducted in silico identification of classical AMP encoding genes from the reference genome (mRhiFer1_v1.p, GCF_004115265.1) of R. ferrumequinum as a representative species. We identified a total of 27 putative AMP genes from R. ferrumequinum, including three alpha-defensins, 21 beta-defensins, one cathelicidin, one hepcidin, and one LEAP-2, using National Center for Biotechnology Information (NCBI) genome blast with query coverage > 50% and e-value < 0.001 (Table S1). The number of AMP genes identified from our analysis was consistent with the gene numbers reported in R. ferrumequinum for each AMP family in the NCBI protein database without finding additional AMP genes.

2.2. Differences in the Number of Functional Cathelicidin Genes among Seven Bat Species Belonging to Different Families

To determine the number of functional cathelicidin genes in bat genomes, 64 unique cathelicidin protein sequences from 36 vertebrate species were blasted using the tblastn mode against the genomes of seven different bat species, R. aegyptiacus, H. armiger, M. natalensis, M. lucifugus, R. ferrumequinum, P. discolor and D. rotundus. Consequently, 12 cathelicidin-like sequences were identified in the analysis with query coverage > 50% and e-value < 0.001 (Table S2). In addition, BLASTP analysis using the same queries against nonredundant protein sequences of each bat species (taxid:9407, 186990, 291302, 59463, 59479, 89673, and 9430) resulted in the identification of 19 matches of cathelicidin-like sequences (Table S3). Further analysis showed that the difference in the number of matches was due to the discrepancy in the annotation of cathelicidin genes between the current genome assemblies and available sequences of bat cathelicidins in the NCBI protein database (Table S4). When the cathelicidin sequences of bats in public databases were taxonomically classified, the number of cathelicidin genes varied across bat species. A single cathelicidin-like gene has been detected in the genomes of H. armiger (LOC109375616) and M. natalensis (LOC107529421). Two putative cathelicidin isoforms, XP_036075573.1 and KAF6473429.1, have been reported in R. aegyptiacus, but they seem to result from two different annotations of an identical gene. In D. rotundus, the antimicrobial activity-containing region was predicted only from one (XP_024421797.1) of the two cathelicidin-like sequences (XP_024421797.1 and XP_024421798.1). Similarly, two cathelicidin-like sequences, KAF6312594.1 and XP_032988513.1, were encoded from an identical locus in R. ferrumequinum, and the structure of KAF6312595.1 lacks an exon. Therefore, five out of the seven species seem to have a single functional cathelicidin gene in their genomes, and more than one cathelicidin could be present in M. lucifugus and P. discolor considering the presence of multiple sequences with antimicrobial signatures, as described in Table S4.

2.3. In Silico Determination of the Antimicrobial Activity Core Regions of Three Bat Cathelicidins

We conducted an in silico prediction of the antimicrobial activity core region for the 19 bat cathelicidin-like sequences to define the signal peptide, conserved CLD domain, and mature peptide region using two different prediction tools, antimicrobial sequence scanning system (AMPA) and database of antimicrobial activity and structure of peptides (DBAASP) (Table S4). Consequently, five putatively functional sequences, KAF6473429.1, XP_019486615.1, XP_006108360.2, XP_035886276.1, and XP_028374415.1, with the complete pre-pro protein structure of cathelicidins and predicted presence of core antimicrobial activity regions were determined (Table S5). Among them, three predicted antimicrobial activity core regions corresponding to the cathelicidins of H. armiger (XP_019486615.1), M. lucifugus (XP_006108360.2) and P. discolor (XP_028374415.1) which were named as ΔHA-CATH, ΔML-CATH, and ΔPD-CATH where “Δ” indicate a deletion of several amino acids from the nascent mature peptide region, namely HA-CATH, ML-CATH, and PD-CATH, respectively, were selected and subjected to the analysis of biochemical features using Antimicrobial Peptide Database (APD3). The three antimicrobial activity core regions or short form peptide regions (ΔHA-CATH, ΔML-CATH, and ΔPD-CATH) consist of 18 to 25 amino acids with 2.08 to 2.93 kDa in molecular weight, whereas the entire mature or long form peptides, HA-CATH, ML-CATH, and PD-CATH, were composed of 35 to 40 residues with 3.98 to 4.64 kDa in molecular weight (Table 1). The ratios of hydrophobic residues and net charges for the short- and long-form peptides were 31 to 50% and +3 to +6, and 30 to 40% and +4 to +8, respectively. Their amino acid sequence similarity to known cathelicidins was highest with cathelicidins from sheep (SMAP-29), dogs (K9CATH), frogs (Temporin-CPb, Palustrin-2CG1), and fish (TP4, HKPLP) and the level ranged from 38.64 to 51.22%.

2.4. Confirmation of the Antibacterial Activity of Three Bat Cathelicidins

Six bat-cathelicidin peptides (ΔHA-CATH, ΔML-CATH, ΔPD-CATH, HA-CATH, ML-CATH, and PD-CATH) were synthesized, and antibacterial activities were evaluated against a panel of 12 reference bacteria, comprised of 6 Gram-negative strains, Escherichia coli, Pseudomonas aeruginosa, Salmonella enterica serovar Typhimurium, Acinetobacter baumannii, Klebsiella pneumoniae subsp. pneumoniae, and Enterobacter cloacae subsp. cloacae, and 6 Gram-positive strains, Staphylococcus aureus, Bacillus cereus, Enterococcus faecalis, Streptococcus agalactiae, Streptococcus dysgalactiae, and Streptococcus equi subsp. zooepidemicus, which are associated with human opportunistic and/or nosocomial pathogens (Table 2). ML-CATH and PD-CATH showed highly potent antibacterial activity against most strains in our bacterial panel with minimum inhibitory concentration (MIC) of 1–35 μg/mL (0.3–8.8 μM) and 3–30 μg/mL (0.7–7.3 μM) except P. aeruginosa (>40 μg/mL or 9.8 μM). In contrast, HA-CATH possessed a moderate antibacterial activity with MIC of 5–34 μg/mL (1.1–7.3 μM) against E. coli, A. baumannii, K. pneumoniae among Gram-negative strains and S. aureus and B. cereus among Gram-positive strains, indicating that the antimicrobial potency and spectrum of bat cathelicidins varied among species. However, ΔHA-CATH, ΔML-CATH, and ΔPD-CATH, the short-form peptides corresponding to the predicted antibacterial core region, did not show any observable antibacterial activity against all strains of the reference bacterial panel, even at a high concentration (>40 μg/mL or 13.7 μM), indicating the prediction of active core region for the bat cathelicidins was not accurate.

2.5. Confirmation of the Antifungal Activity of Three Bat Cathelicidins

Six bat-cathelicidin peptides (ΔHA-CATH, ΔML-CATH, ΔPD-CATH, HA-CATH, ML-CATH, and PD-CATH) were evaluated for the presence of antifungal activity against three different fungal species, Candida albicans (C. albicans), Cryptococcus cuniculi (C. cuniculi), and Mucor hiemalis (M. hiemalis) (Table 3; Figure 1). The antifungal effect of the yeast-type fungi, C. albicans and C. cuniculi, was estimated in a manner similar to the antibacterial activity method. However, the antifungal activity of M. hiemalis, a mold-type fungus, was estimated based on the inhibition of the metabolic rate against hyphae. The treatment of ML-CATH and PD-CATH showed effective antifungal activities against the two yeast forms fungi with MIC of 5 to 40 μg/mL (1.2 to 10.1 μM) and M. hiemalis with a similar rate of metabolic inhibition (~80%) to fluvastatin sodium. In addition, intermediate levels of antifungal activity were observed against C. albicans and C. cuniculi from the treatment of ΔPD-CATH which did not show antibacterial activity against our reference bacterial panel, indicating the inherently strong antimicrobial potency of PD-CATH compared to the other bat cathelicidins tested in this study and differences in the effect of ΔPD-CATH against bacterial and fungal cell walls. In addition, HA-CATH only showed slight antifungal activity (MIC = 45 μg/mL or 9.7 μM) against C. cuniculi, indicating limited activity against fungi which is consistent with the antibacterial activity of HA-CATH (Table 2). The antifungal activity of the tested peptides seems to be higher for yeast than mold-type fungi. Consistent with the results of antibacterial analysis, ΔHA-CATH, ΔML-CATH, and ΔPD-CATH did not show antifungal effects against the reference strains, except ΔPD-CATH against C. cuniculi.

2.6. Negligible Cytotoxicity of Bat Cathelicidins against Human Keratinocytes

The level of cytotoxicity is important for the therapeutic use of AMPs as antimicrobial agents. To estimate the cytotoxicity of HA-CATH, ML-CATH, and PD-CATH in mammalian cells, the viability of HaCaT cells was evaluated for each peptide (Table 4). The results showed that no cytotoxicity was observed at 64 μg/mL concentration, and only less than 3% of cells were affected at 160 μg/mL for both HA-CATH and ML-CATH. For PD-CATH, the cytotoxicity was not observed at 64 μg/mL, but the cell viability was decreased to 10.4% at 160 μg/mL concentration, showing significant cytotoxicity at the high concentration. Therefore, the cytotoxicity on mammalian cells showed a positive correlation with the potency of bat cathelicidins. In addition, melittin, an AMP known to have a potent cell lytic activity, showed severe cell lysis, with 7.2% cell viability at 64 μg/mL concentration. These results indicate that the cytotoxic effect of HA-CATH and ML-CATH is at a minimal level, and given the effective bactericidal concentration of PD-CATH (3–30 μg/mL), the cytotoxicity is within the acceptable range for therapeutic use.

3. Discussion

Fifty genome sequences of diverse bat species with varying degrees of coverage are currently available in the public database (https://www.ncbi.nlm.nih.gov/genome, accessed on 22 June 2022). Considering the extreme phylogenetic diversity of bats, comparative analysis of their AMP genes may reveal interesting outcomes regarding the genetic and functional diversity of these molecules. The number of classical AMP genes of R. ferrumequinum was similar to that of naked mole rats when compared to other mammals (Table S6) [25,44,45,46,47,48]. However, only limited information is available on the genetic and biological characteristics of bat AMPs, including cathelicidins. In addition, the results of our analysis to identify cathelicidin orthologs from the seven selected bat species showed few differences in annotations compared to the currently predicted protein sequences of cathelicidins in NCBI (Table S4), suggesting the need of more cDNA information, especially for isoform-like sequences. Cathelicidins are probably the most potent AMPs in mammals, with variations in their activities in cases of multiple cathelicidins within a species [49]. Although few studies reported bioinformatic prediction of bat cathelicidins based on available sequencing information and experimental characterization of their antimicrobial activity [43,50], the number of tested bacterial strains and peptides was highly limited, and other characteristics such as fungicidal and metabolically inhibitory activities against fungi and cytotoxicity have not been validated experimentally. To the best of our knowledge, this is the first report of the antimicrobial activity of bat cathelicidins against diverse bacteria and fungi associated with human health and food pathogens.
A previous study briefly reported the bactericidal effect of a bat cathelicidin, Ml-LN-35, which is identical to ML-CATH, against E. coli and B. cereus [50]. The results were consistent with the results in this study with MICs ranging from 0.5 to 0.8 μM (Table 2). In contrast, antibacterial activity against S. aureus ATCC 6538 was not detected in the previous study, while the peptide showed highly potent anti-staphylococcal activity (1 μg/mL or 0.3 μM) as shown in Table 2. This inconsistency might be attributed to different experimental condition including media composition and peptide synthesis with some modification such as acetylation and amidation in the N- and C-terminal end, respectively, which can affect the nature and structure of AMPs and influence biological activity and stability of the molecules in turn [51,52,53].
In our functional analysis of bat cathelicidins, all three mature peptide region-derived peptides showed antibacterial activity against diverse bacterial strains. However, the truncated peptides corresponding to the predicted antimicrobial activity core region as a minimal unit for biological activity did not show substantial antibacterial activity, although the two different forms of peptides shared common characteristics in hydrophobicity and positive net charge (Table 1). This is somewhat different from the results of our previous studies to characterize cathelicidins of non-bat species [25,47,54].
The alteration of primary sequences and/or physicochemical properties of α-helical cationic AMPs, especially cathelicidins, could affect the formation of amphipathic α-helical structures, hydrophobicity, the density and distribution of charge locally or overall, and minimal length, influencing biological activities, including antimicrobial and mammalian cytotoxicity [55,56,57,58]. For example, a previous study showed that IG-13, corresponding to the minimal active domain of LL-37, a well-known human cathelicidin, was approximately 30-fold less active than LL-37 [58]. The smallest fragment, LL-19, did not form a clear α-helix in the buffer, although the peptide region originated from the α-helical region [55]. In addition, the hydrophobic tail length of the peptides could affect their ability to insert and penetrate the bacterial membrane [57]. Therefore, the loss of antimicrobial activity in the predicted antimicrobial activity core region peptides 18 to 25 amino acids in size of bat cathelicidins could be influenced by peptide length. Further studies are necessary to determine the amino acids critical for maintaining the antimicrobial activity of nascent mature peptides.
Among the mature peptide regions of the three cathelicidins, the specificity and strength of the antibacterial activities of ML-CATH and PD-CATH were highly similar (Table 2). This could be attributed to the similarities in peptide length, hydrophobicity, and net charge between the two peptides despite their sequence differences (Table 1; Figure S1).
In addition to the extensive functional studies on LL-37 and CRAMP, recent interests in AMPs as potent antibiotic agents lead to the identification and characterization of novel cathelicidins from diverse species including Python bivittatus (Burmese python), Heterocephalus glaber (naked mole rat), Monodelphis domestica (gray short-tailed opossum), Lates calcarifer (Asian sea bass) and Sarcophilus harrisii (Tasmanian devil) [25,47,54,59,60]. These studies suggested that cathelicidins are promising candidates to develop antimicrobial therapeutics considering potent broad-spectrum bactericidal activity, minimal cytotoxicity to mammalian cells, and inducibility of lower antibiotic resistance compared to classical antibiotics. Interestingly, bat cathelicidins, ML-CATH and PD-CATH, showed much broader and stronger antibacterial activities than all the peptides described above (Table 2). Furthermore, ML-CATH and PD-CATH showed milder toxicity compared to ModoCath1 which is one of the cathelicidins showing broad spectrum with high potency. The viability of human keratinocyte was affected at a five-fold higher concentration (160 µg/mL) in PD-CATH comparing to ModoCath1 (32 µg/mL) (Table 4). These characteristics could make ML-CATH and PD-CATH more suitable for pharmaceutical applications than other cathelicidins. However, further studies are still required to overcome intrinsic disadvantages of AMPs including high production cost, low stability in vivo, and deficiency of proper delivery system. An effective production system of recombinant cathelicidins using an engineered green fluorescent protein (GFP) as a tag protein was attempted for large-scale production of AMPs [61]. The conjugation of AMPs with nanoparticles has been suggested to overcome poor stability of AMPs in biological fluids and proteolytic degradation and to control pharmacokinetics [62,63,64].
In mammals, the copy number of cathelicidin genes varies among species, and a single functional cathelicidin gene is present in dogs, mice, and humans, whereas more than one functional gene is found in pigs, goats, sheep, cattle, and whales [65]. Intraspecies copy number variation has also been reported in pigs [48]. It is interesting to note that the antimicrobial potency of ML-CATH and PD-CATH, obtained from multiple cathelicidin gene-carrying bat species, was higher and more broadly effective against diverse bacterial strains than HA-CATH, obtained from species carrying only a single cathelicidin gene. The expansion of cathelicidin genes in certain bat species may be due to the functional importance of the molecule for their survival. It may also enhance the antimicrobial potency of paralogous genes in parallel.
The variations in the number of cathelicidin genes among different bat species can be attributed to the large phylogenetic diversity of bats in the order Chiroptera [66]. Moreover, in our analysis of the seven bat species, all annotated cathelicidin genes were mapped to the CDC25A and NME6 interval, the evolutionarily conserved cathelicidin residing region in mammals [48], except M. lucifugus with less complete genome information, suggesting the evolutionary conservation of synteny in bats, similar to other mammals. Therefore, the copy number of cathelicidins in bats could vary from a single copy to multiple copies depending on the species, which is consistent with a previous analysis [43]. However, further refinement of the current bat genome assemblies and annotations is needed to determine the copy number of cathelicidin family genes more accurately.
Endogenous AMPs also showed antiviral activity against diverse viruses, including the human influenza virus, human immunodeficiency virus, rabies virus, West Nile virus (WNV), and herpes simplex virus (HSV) [25,67,68,69,70]. Among them, gramicidin S and melittin from Bacillus brevis and bee venom, respectively, showed noteworthy anti-SARS-CoV-2 effects with viral particle reduction compared to remdesivir, a repurposed drug used to treat the infection [71]. Antiviral activities against HSV and WNV were also reported for other cathelicidins, such as LL-37, Pb-CATH4, Hg-CATH, and ModoCath5 [25,70,72]. Because of the therapeutic potential of these peptides, it will be interesting to evaluate the antiviral activity of bat cathelicidins, which also showed negligible cytotoxicity to mammalian cells in this study. The antibacterial and antifungal activities of HA-CATH, ML-CATH, and PD-CATH observed in this study can serve as references for the potency and antimicrobial spectrum analyses of bat cathelicidins.

4. Materials and Methods

4.1. In Silico Identification of Cathelicidin-like Sequences in the Genome of Seven Bat Species

Four hundred sequences belonging to five major AMP families of vertebrates, alpha-defensins, beta-defensins, cathelicidin, hepcidin, and LEAP-2, were obtained similarly to the previous studies [25,47]. Briefly, the sequences of 2477 nonredundant AMPs were downloaded from UniProtKB/Swiss-Prot (http://www.uniprot.org/uniprot/, accessed on 16 March 2022; [73]) using the query “antimicrobial peptide AND reviewed: yes”. After removing sequences originated from nonvertebrates, the sequences corresponding to five major AMP families described above were collected (Table S7) [25]. The sequences were blasted using the NCBI blastp mode (https://www.ncbi.nlm.nih.gov/, accessed on 16 March 2022) against two different reference genomes of R. ferrumequinum (mRhiFer1_v1.p, GCF_004115265.1 and mRhiFer1.p, GCA_014108255.1) with the genome coverage of 52.8 x and 73.3 x, respectively, in NCBI. In addition, 64 cathelicidin-like sequences from 400 vertebrate AMPs with both the CLD and mature domains were subjected to NCBI blastp and tblastn analyses against the reference genomes of additional six bat species with the coverage of 7 x to 218.6 x in NCBI genome database (https://www.ncbi.nlm.nih.gov/genome, accessed on 18 March 2022), including R. aegyptiacus (mRouAeg1.p, GCF_014176215.1), H. armiger (ASM189008v1, GCF_001890085.1), M. natalensis (Mnat.v1, GCF_001595765.1), M. lucifugus (Myoluc2.0, GCF_014108235.1), P. discolor (mPhyDis1.pri.v3, GCF_004126475.2, and mPhyDis1_v1.p, GCF_004126475.1), and D. rotundus (ASM294091v2, GCF_002940915.1). The prediction of genetic structures and positional confirmation of putative cathelicidin sequences in bat genomes were conducted using the NCBI (https://www.ncbi.nlm.nih.gov/genome/gdv, accessed on 28 March 2022) and the Bat1K genome browser (https://genome-public.pks.mpg.de, accessed on 1 April 2022). The last accession of genome databases was 1 May 2022.

4.2. In Silico Prediction of Putative Cathelicidins for Their Biological Activities

The in silico determination of the signal peptide and CLD regions was conducted using the SignalP 4.1 server (https://services.healthtech.dtu.dk/service.php?SignalP-4.1, accessed on 16 April 2022; [74]) and HMMER (https://www.ebi.ac.uk/Tools/hmmer/, accessed on 17 April 2022; [75]), respectively. AMPA (https://tcoffee.crg.eu/apps/ampa/do, accessed on 16 April 2022; [76]) and DBAASP (https://dbaasp.org/, accessed on 16 April 2022; [77]) were used to predict the potential antimicrobial activity domains using the default window size and threshold value. Protein secondary structures were predicted using the PSIPRED workbench (http://bioinf.cs.ucl.ac.uk/psipred/, accessed on 16 April 2022; [78]). The biochemical characteristics of predicted mature peptide sequences, including hydrophobicity, net charge, molecular weight, and similarity to known AMPs, were estimated using APD3 (https://aps.unmc.edu/, accessed on 1 May 2022; [79]).

4.3. Peptide Synthesis

Peptides corresponding to the predicted active core and mature peptide regions of cathelicidins were generated via solid-phase peptide synthesis, purified via high-performance liquid chromatography (HPLC), and confirmed by mass spectroscopy using a commercial service (GenScript, Piscataway Township, NJ, USA). Peptide names were given according to the acronym of binomial nomenclature of the species followed by “CATH” indicating cathelicidin. The peptide names and sequences for each species are ΔHA-CATH (N-LLRRGGRKIGQGLERIGQRIQGF-C) and HA-CATH (N-ILGRLRDLLRRGGRKIGQGLERIGQRIQGFFSNREPMEES-C) for H. armiger, ΔML-CATH (N-GIFILKHRRPIGRGIEIT-C) and ML-CATH (N-LNPLIKAGIFILKHRRPIGRGIEITGRGIKKFFSK-C) for M. lucifugus, and ΔPD-CATH (N-IAGRIAGKLIGDAINRHRERNRQRR-C) and PD-CATH (N-ILGPALRIGGRIAGRIAGKLIGDAINRHRERNRQRRG-C) for P. discolor. The synthesized peptides were dissolved in distilled water at a concentration of 4 mg/mL.

4.4. Evaluation of Antibacterial Activity

The antibacterial activities of the chemically synthesized peptides were evaluated against a panel of bacteria consisting of 6 Gram-negative strains, including E. coli ATCC 25922 (American Type Culture Collection, Manassas, VA, USA), P. aeruginosa ATCC 27853, S. Typhimurium ATCC 14028, A. baumannii KCTC (Korean Collection for Type Cultures, Jeongeup, Korea) 23254, K. pneumoniae KCTC 1726, and E. cloacae ATCC 13047, and 6 Gram-positive strains, including S. aureus ATCC 6538, B. cereus ATCC 10876, E. faecalis ATCC 29212, S. agalactiae ATCC 27956, S. dysgalactiae ATCC 27957, and S. equi ATCC 43079. Ampicillin (Sigma-Aldrich, St. Louis, MO, USA), chloramphenicol (Sigma-Aldrich), and gentamicin sulfate (Sigma-Aldrich) were used as controls for antibacterial activity. The MIC was determined by a colorimetric method using the Microbial Viability Assay Kit-WST (Dojindo, Kumamoto, Japan) according to the manufacturer’s manual and the Clinical and Laboratory Standards Institute (CLSI) guidelines (2018). Four colonies of each bacterium were inoculated into 5 mL of Luria-Bertani (LB) broth medium (BD Bioscience, Franklin Lakes, NJ, USA) or brain heart infusion broth (BHIB; BD Bioscience) at 37 °C for 6 h with shaking at 220 rpm. E. faecalis and streptococci were cultured in BHIB because of their slow growth, whereas all other bacteria were cultured in LB medium. The cells were washed by sterile saline (0.9% NaCl) twice and adjusted in a single well of a 96-well plate at a cell density of 1.5 × 105 CFU/well in accordance with 0.5 McFarland standard. Subsequently, 180 µL/well of fresh Mueller-Hinton broth (MHB; BD Bioscience), except for E. faecalis and streptococci using BHIB was added to the plate using BHIB. Final concentrations of 1–40 µg/mL (0.3–19.3 μM) of each peptide and reference antibiotic were added to each well. The plates were then incubated at 37 °C for 6 h. Subsequently, 10 µL of the coloring reagent was added, and the mixture was incubated at 37 °C for 2 h. The absorbance of each well was measured at 450 nm using a microplate spectrophotometer (Bio-Rad, xMark spectrophotometer, Hercules, CA, USA). MIC values were determined when the difference in the absorbance values between treatments and blanks (media and coloring reagent only) was less than 0.05. Experiments were performed in triplicate.

4.5. Evaluation of Antifungal Activity

The antifungal activities of the peptides were evaluated against two yeast fungi, C. albicans KCTC 7270 and C. cuniculi KCTC 17232, and a mold fungus, M. hiemalis KCTC 26779. C. albicans and C. cuniculi were cultured in yeast mold broth (YMB; BD Science) agar at 35 °C and 25 °C, respectively. Cells were prepared in YMB for C. albicans and RPMI 1640 medium (with l-glutamine, without phenol red and sodium bicarbonate, Sigma Aldrich, Saint Louis, MO, USA) and adjusted to pH 7.0 with 0.165 M 3-(N-morpholino) propanesulfonic acid (MOPS) (iNtRON Biotechnology, Seoul, Korea), as previously described [80] for C. cuniculi. According to the 0.5 McFarland standard, cells at 1.5 × 105 CFU/well were seeded into a 96-well plate containing 180 µL of fresh medium. Subsequently, peptides at final concentrations of 5–45 ug/mL (1.2–19.3 μM) peptides were added to each well. Ciclopirox (≥98% (HPLC), Sigma Aldrich) was used as a reference for the antifungal activity. The antifungal assay was conducted in the same manner as the antibacterial assay described above, according to the manufacturer’s protocol. Briefly, C. albicans and C. cuniculi were incubated at 35 °C for 6 h and 25 °C for 47 h, respectively, and 10 µL of the coloring reagent was added and further incubated for 18 h and 25 h. The MIC values were determined in the same way as for the antibacterial assay.
To evaluate the antifungal activity of peptides against M. hiemalis, the metabolic activity of hyphae was measured using the same colorimetric assay as described above, with slight modifications to a previous study [81]. Briefly, M. hiemalis was streaked onto malt extract broth (BD Biosciences) agar in a Petri dish and cultured at 25 °C for seven days. Then, conidia were collected from the dish by adding phosphate buffered saline (PBS) containing 0.05% Tween-80 (Sigma Aldrich) and filtered through sterilized 40 um cell strainer (SPL Life Science, Pocheon, Korea) to remove hyphae. After washing twice with PBS and centrifugation, 5 × 104 conidia were counted using a hemocytometer, suspended in 50 µL of RPMI 1640 as described above, and seeded in wells of a 96-well plate. After 16 h of incubation at 25 °C for hyphae formation, the plate was centrifuged at 3000× g and the media were changed. The peptides were then added at final concentrations of 20, 40, and 80 µg/mL. Fluvastatin sodium (70 µM; Sigma-Aldrich) was used as a reference for antifungal activity. For the coloring reaction, 10 µL of the coloring solution (Microbial Viability Assay Kit-WST, Dojindo) was added to each well and incubated for 2 h. Subsequently, the plate was centrifuged, and the supernatant was transferred to a new 96-well plate. The absorbance at 450 nm (for the treated group and control) and 650 nm (for the background) was measured for each well using a microplate reader (xMarkTM spectrophotometer; Bio-Rad). The metabolic activity of the fungal cells was calculated using the following equation. All experiments were performed in triplicate.
Fungal metabolic activity (%) = 100 × (OD (treated group) − OD (background))/(OD (nevative control) − OD (Background))

4.6. In Vitro Cytotoxicity Assay

HaCaT cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM; HycloneTM, Logan, UT, USA) with 10% fetal bovine serum (FBS; HycloneTM) and 1% penicillin/streptomycin (Pen-Strep; HycloneTM) up to 80% confluence at 37 ℃ and 5% CO2. Cells were detached from the plate by adding 0.25% trypsin-EDTA solution (GibcoTM; Carlsbad, CA, USA). In a 96-well plate, 2 × 104 cells were seeded per well and cultured for 24 h. Subsequently, the medium was replaced with 100 μL DMEM (HycloneTM) containing 10% heat-inactivated FBS. Peptides were added to each well at concentrations of 64 μg/mL and 160 μg/mL. As references for complete cell lysis, 64 μg/mL melittin (Sigma-Aldrich) and Triton X-100 (Sigma-Aldrich) were used. Untreated cells were used as negative controls. After 24 h incubation at 37 ℃ and 5 % CO2, 10 μL of coloring solution (Cell Proliferation Reagent WST-1TM; Sigma Aldrich) was added to each well according to the manufacturer’s instructions. The absorbance at 450 nm (treated group and control) and 650 nm (background) was measured in each well using a microplate reader (xMarkTM spectrophotometer; Bio-Rad). The cell viability was calculated using the following equation. All experiments were performed in triplicate.
Cell viability (%) = 100 × (OD (treated group) − OD(Background))/(OD (negative control) − OD (Background))

4.7. Statistical Analysis

Statistical analysis was performed using Student’s paired t-test with a two-tailed distribution using RStudio (https://www.R-project.org/, accessed on 13 July 2022).

5. Conclusions

In summary, species-specific differences in genetic and functional characteristics of cathelicidins from seven bat species were studied. We listed 19 cathelicidin-like sequences from seven bat species and tested the biological activity of three cathelicidins (HA-CATH, ML-CATH, and PD-CATH). Among them, ML-CATH and PD-CATH showed potent antibacterial activity against a broad spectrum of Gram-negative and Gram-positive bacteria. The potency of antibacterial activity of ML-CATH and PD-CATH was comparable to that of other highly active cathelicidins reported in vertebrates. ML-CATH and PD-CATH also showed antifungal activities against C. albicans, C. cuniculi, and M. hiemalis. The cytotoxicity on HaCaT cells showed a positive correlation with the potency of bat cathelicidins, but given the effective bactericidal concentration, the cytotoxicity of bat cathelicidins is within the acceptable range for therapeutic exploitation. ML-CATH and PD-CATH are promising candidates for the development of peptide antibiotics, but further studies including peptide production and in vivo activity are in need.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/antibiotics11080989/s1, Figure S1: Comparison of the amino acid sequences and properties of the three bat cathelicidins evaluated in this study, Table S1: List of the putative antimicrobial peptide (AMP) genes of Rhinolophus ferrumequinum, Table S2: Results of the National Center for Biotechnology Information (NCBI) tblastn analysis of the genomes of seven bat species using 64 vertebrate cathelicidin sequences as queries, Table S3: Results of the NCBI blastp analysis of the genomes of seven bat species using 64 vertebrate cathelicidin sequences as queries, Table S4: Results of the in silico analysis of annotated putative cathelicidin-like sequences of seven bat species in NCBI, Table S5: In silico prediction of the antimicrobial activity-containing regions of five functional bat cathelicidins, Table S6: Comparison of the AMP family gene numbers of six mammalian species, Table S7: A list of antimicrobial peptides queries for the blastp analysis of antimicrobial peptides.

Author Contributions

Conceptualization, M.C., H.-s.C. and C.P.; Writing—Original Draft Preparation, M.C., H.-s.C. and C.P.; Investigation, M.C., H.-s.C. and B.A.; Formal Analysis, M.C., H.-s.C. and B.A.; Resources, S.P. and S.N.; C.P. provided advice with data interpretation. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture and Forestry (IPET) through the High Value-Added Food Technology Development Program, funded by the Ministry of Agriculture, Food and Rural Affairs (MAFRA) [322017-3], and the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) [2021R1A2C3010505].

Data Availability Statement

The dataset analyzed during this study can be found in the NCBI database under the following accession numbers: XM_036219680.1, XM_019631070.1, XM_016202809.1, XM_006108300.3, XM_006108299.3, XM_006108298.3, XM_006108738.3, XM_033132622.1, XM_036030383.1, XM_028518614.2, XM_024566029.2, XM_024566030.2, XP_032960802.1, XP_032960796.1, XP_032960740.1, XP_032961053.1, XP_032991305.1, XP_032960777.1, XP_032960718.1, XP_032960717.1, XP_032959658.1, XP_032959919.1, XP_032959522.1, XP_032959179.1, XP_032952314.1, XP_032950360.1, XP_032951348.1, XP_032951380.1, XP_032951378.1, XP_032951381.1, XP_032951382.1, XP_032951383.1, XP_032951385.1, XP_032950352.1, XP_032959847.1, XP_032991304.1, XP_032988513.1, XP_032984445.1, XP_032952959.1, XP_036075573.1, KAF6473429.1, XP_019486615.1, XP_016058295.1, XP_006108362.1, XP_006108361.1, XP_006108360.2, XP_006108800.3, KAF6312595.1, KAF6312594.1, XP_035886276.1, XP_028374415.1, KAF6098810.1, KAF6098809.1, KAF6098808.1, KAF6100376.1, XP_024421797.1, and XP_024421798.1. Additional data supporting the conclusions of this study are available in the Supplementary Materials File.

Conflicts of Interest

The authors declare that they have competing financial interests. The peptide sequences, ML-CATH and PD-CATH were subject to domestic and foreign patent applications by Konkuk University.

References

  1. Jones, K.E. Chiroptera (Bats). In eLS, 1st ed.; Wiley: Hoboken, NJ, USA, 2006. [Google Scholar]
  2. Teeling, E.C.; Jones, G.; Rossiter, S.J. Phylogeny, Genes, and Hearing: Implications for the Evolution of Echolocation in Bats. In Bat Bioacoustics; Fenton, M.B., Grinnell, A.D., Popper, A.N., Fay, R.R., Eds.; Springer: New York, NY, USA, 2016; pp. 25–54. [Google Scholar]
  3. Prat, Y.; Azoulay, L.; Dor, R.; Yovel, Y. Crowd vocal learning induces vocal dialects in bats: Playback of conspecifics shapes fundamental frequency usage by pups. PLoS Biol. 2017, 15, e2002556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Simmons, N.B. An Eocene Big Bang for Bats. Science 2005, 307, 527–528. [Google Scholar] [CrossRef] [PubMed]
  5. Huang, Z.; Jebb, D.; Teeling, E.C. Blood miRNomes and transcriptomes reveal novel longevity mechanisms in the long-lived bat, Myotis myotis. BMC Genom. 2016, 17, 906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Drexler, J.F.; Corman, V.M.; Müller, M.A.; Maganga, G.D.; Vallo, P.; Binger, T.; Gloza-Rausch, F.; Cottontail, V.M.; Rasche, A.; Yordanov, S.; et al. Bats host major mammalian paramyxoviruses. Nat. Commun. 2012, 3, 796. [Google Scholar] [CrossRef] [Green Version]
  7. Goldstein, T.; Anthony, S.J.; Gbakima, A.; Bird, B.H.; Bangura, J.; Tremeau-Bravard, A.; Belaganahalli, M.N.; Wells, H.L.; Dhanota, J.K.; Liang, E.; et al. The discovery of Bombali virus adds further support for bats as hosts of ebolaviruses. Nat. Microbiol. 2018, 3, 1084–1089. [Google Scholar] [CrossRef]
  8. Cui, J.; Li, F.; Shi, Z.-L. Origin and evolution of pathogenic coronaviruses. Nat. Rev. Microbiol. 2019, 17, 181–192. [Google Scholar] [CrossRef] [Green Version]
  9. Jebb, D.; Huang, Z.; Pippel, M.; Hughes, G.M.; Lavrichenko, K.; Devanna, P.; Winkler, S.; Jermiin, L.S.; Skirmuntt, E.C.; Katzourakis, A.; et al. Six reference-quality genomes reveal evolution of bat adaptations. Nature 2020, 583, 578–584. [Google Scholar] [CrossRef]
  10. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [Green Version]
  11. Matin, M.M.; Uzzaman, M.; Chowdhury, S.A.; Bhuiyan, M.M.H. In vitro antimicrobial, physicochemical, pharmacokinetics and molecular docking studies of benzoyl uridine esters against SARS-CoV-2 main protease. J. Biomol. Struct. Dyn. 2022, 40, 3668–3680. [Google Scholar] [CrossRef]
  12. Yang, H.; Rao, Z. Structural biology of SARS-CoV-2 and implications for therapeutic development. Nat. Rev. Microbiol. 2021, 19, 685–700. [Google Scholar] [CrossRef]
  13. Wang, L.-F.; Gamage, A.M.; Chan, W.O.Y.; Hiller, M.; Teeling, E.C. Decoding bat immunity: The need for a coordinated research approach. Nat. Rev. Immunol. 2021, 21, 269–271. [Google Scholar] [CrossRef]
  14. Banerjee, A.; Baker, M.L.; Kulcsar, K.; Misra, V.; Plowright, R.; Mossman, K. Novel Insights into Immune Systems of Bats. Front. Immunol. 2020, 11, 26. [Google Scholar] [CrossRef]
  15. Hayman, D.T.S. Bat tolerance to viral infections. Nat. Microbiol. 2019, 4, 728–729. [Google Scholar] [CrossRef]
  16. Bhattacharyya, S. Inflammation During Virus Infection: Swings and Roundabouts. In Dynamics of Immune Activation in Viral Diseases; Bramhachari, P.V., Ed.; Springer: Singapore, 2020; pp. 43–59. [Google Scholar]
  17. Teeling, E.; Vernes, S.C.; Davalos, L.M.; Ray, D.A.; Gilbert, M.T.P.; Myers, E.; Bat1K Consortium. Bat biology, genomes, and the Bat1K project: To generate chromosome-level genomes for all living bat species. Annu. Rev. Anim. Biosci. 2018, 6, 23–46. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, G.; Cowled, C.; Shi, Z.; Huang, Z.; Bishop-Lilly, K.A.; Fang, X.; Wynne, J.W.; Xiong, Z.; Baker, M.L.; Zhao, W.; et al. Comparative Analysis of Bat Genomes Provides Insight into the Evolution of Flight and Immunity. Science 2013, 339, 456–460. [Google Scholar] [CrossRef] [Green Version]
  19. Brook, C.E.; Dobson, A.P. Bats as ‘special’ reservoirs for emerging zoonotic pathogens. Trends Microbiol. 2015, 23, 172–180. [Google Scholar] [CrossRef]
  20. Zhan, L.-J.; Gallo, R.L. Antimicrobial peptides. Curr. Biol. 2016, 26, R14–R19. [Google Scholar] [CrossRef]
  21. Nizet, V.; Ohtake, T.; Lauth, X.; Trowbridge, J.; Rudisill, J.; Dorschner, R.A.; Pestonjamasp, V.; Piraino, J.; Huttner, K.; Gallo, R.L. Innate antimicrobial peptide protects the skin from invasive bacterial infection. Nature 2001, 414, 454–457. [Google Scholar] [CrossRef]
  22. Mor, A. Multifunctional host defense peptides: Antiparasitic activities. FEBS J. 2009, 276, 6474–6482. [Google Scholar] [CrossRef]
  23. De Ullivarri, M.F.; Arbulu, S.; Garcia-Gutierrez, E.; Cotter, P.D. Antifungal Peptides as Therapeutic Agents. Front. Cell. Infect. Microbiol. 2020, 10, 105. [Google Scholar] [CrossRef]
  24. Carriel-Gomes, M.C.; Kratz, J.M.; Barracco, M.A.; Bachére, E.; Barardi, C.R.M.; Simões, C.M.O. In vitro antiviral activity of antimicrobial peptides against herpes simplex virus 1, adenovirus, and rotavirus. Mem. Inst. Oswaldo Cruz 2007, 102, 469–472. [Google Scholar] [CrossRef]
  25. Cho, H.S.; Yum, J.; Larivière, A.; Lévêque, N.; Le, Q.V.C.; Ahn, B.; Jeon, H.; Hong, K.; Soundrarajan, N.; Kim, J.H.; et al. Opossum Cathelicidins Exhibit Antimicrobial Activity Against a Broad Spectrum of Pathogens Including West Nile Virus. Front. Immunol. 2020, 11, 347. [Google Scholar] [CrossRef] [Green Version]
  26. Beaumont, P.E.; McHugh, B.; Findlay, E.G.; Mackellar, A.; MacKenzie, K.J.; Gallo, R.L.; Govan, J.R.W.; Simpson, A.J.; Davidson, D.J. Cathelicidin Host Defence Peptide Augments Clearance of Pulmonary Pseudomonas aeruginosa Infection by Its Influence on Neutrophil Function In Vivo. PLoS ONE 2014, 9, e99029. [Google Scholar] [CrossRef]
  27. Choi, K.-Y.G.; Napper, S.; Mookherjee, N. Human cathelicidin LL-37 and its derivative IG-19 regulate interleukin-32-induced inflammation. Immunology 2014, 143, 68–80. [Google Scholar] [CrossRef]
  28. Niyonsaba, F.; Ushio, H.; Nakano, N.; Ng, W.; Sayama, K.; Hashimoto, K.; Nagaoka, I.; Okumura, K.; Ogawa, H. Antimicrobial Peptides Human β-Defensins Stimulate Epidermal Keratinocyte Migration, Proliferation and Production of Proinflammatory Cytokines and Chemokines. J. Investig. Dermatol. 2007, 127, 594–604. [Google Scholar] [CrossRef] [Green Version]
  29. Tomasinsig, L.; Zanetti, M. The Cathelicidins-Structure, Function and Evolution. Curr. Protein Pept. Sci. 2005, 6, 23–34. [Google Scholar] [CrossRef]
  30. Sørensen, O.E.; Follin, P.; Johnsen, A.H.; Calafat, J.; Tjabringa, G.S.; Hiemstra, P.; Borregaard, N. Human cathelicidin, hCAP-18, is processed to the antimicrobial peptide LL-37 by extracellular cleavage with proteinase 3. Blood 2001, 97, 3951–3959. [Google Scholar] [CrossRef] [Green Version]
  31. Gudmundsson, G.H.; Agerberth, B.; Odeberg, J.; Bergman, T.; Olsson, B.; Salcedo, R. The Human Gene FALL39 and Processing of the Cathelin Precursor to the Antibacterial Peptide LL-37 in Granulocytes. Eur. J. Biochem. 1996, 238, 325–332. [Google Scholar] [CrossRef] [PubMed]
  32. Van Harten, R.M.; Van Woudenbergh, E.; Van Dijk, A.; Haagsman, H.P. Cathelicidins: Immunomodulatory Antimicrobials. Vaccines 2018, 6, 63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Nagasundarapandian, S.; Cho, H.-S.; Prathap, S.; Kang, M.; Choi, M.; Lee, Y.; Jeon, H.; Song, H.; Kim, J.-H.; Park, C. Cathelicidin ΔPb-CATH4 derived from Python bivittatus accelerates the healing of Staphylococcus aureus-infected wounds in mice. Amino Acids 2021, 53, 313–317. [Google Scholar] [CrossRef] [PubMed]
  34. Turner, J.; Cho, Y.; Dinh, N.-N.; Waring, A.J.; Lehrer, R.I. Activities of LL-37, a Cathelin-Associated Antimicrobial Peptide of Human Neutrophils. Antimicrob. Agents Chemother. 1998, 42, 2206–2214. [Google Scholar] [CrossRef] [Green Version]
  35. Gallo, R.L.; Kim, K.J.; Bernfield, M.; Kozak, C.A.; Zanetti, M.; Merluzzi, L.; Gennaro, R. Identification of CRAMP, a Cathelin-related Antimicrobial Peptide Expressed in the Embryonic and Adult Mouse. J. Biol. Chem. 1997, 272, 13088–13093. [Google Scholar] [CrossRef] [Green Version]
  36. Liang, W.; Enée, E.; Andre-Vallee, C.; Falcone, M.; Sun, J.; Diana, J. Intestinal Cathelicidin Antimicrobial Peptide Shapes a Protective Neonatal Gut Microbiota Against Pancreatic Autoimmunity. Gastroenterology 2022, 162, 1288–1302.e16. [Google Scholar] [CrossRef]
  37. Liu, Z.; Wu, J.; Qin, Z.; Dong, C.; Yang, H.; Sun, J.; Xu, W.; Wei, L. Endogenous cathelicidin is required for protection against ZIKV-caused testis damage via inactivating virons. Antivir. Res. 2022, 198, 105248. [Google Scholar] [CrossRef]
  38. Rather, I.A.; Sabir, J.S.M.; Asseri, A.H.; Ali, S. Antifungal Activity of Human Cathelicidin LL-37, a Membrane Disrupting Peptide, by Triggering Oxidative Stress and Cell Cycle Arrest in Candida auris. J. Fungi 2022, 8, 204. [Google Scholar] [CrossRef]
  39. Feng, G.; Wei, L.; Che, H.; Shen, Y.; Mi, K.; Bian, H.; Yang, H.; Wu, J.; Mu, L. Cathelicidin-NV from Nanorana ventripunctata effectively protects HaCaT cells, ameliorating ultraviolet B-induced skin photoaging. Peptides 2022, 150, 170712. [Google Scholar] [CrossRef]
  40. Wang, A.; Zhang, F.; Guo, Z.; Chen, Y.; Zhang, M.; Yu, H.; Wang, Y. Characterization of a Cathelicidin from the Colubrinae Snake, Sinonatrix annularis. Zool. Sci. 2019, 36, 68–76. [Google Scholar] [CrossRef]
  41. Kalita, D.J. Characterization of cathelicidin gene from buffalo (Bubalus bubalis). Afr. J. Biotechnol. 2015, 14, 758–763. [Google Scholar]
  42. Dalla Valle, L.; Benato, F.; Paccanaro, M.C.; Alibardi, L. Bioinformatic and molecular characterization of cathelicidin-like peptides isolated from the green lizard Anolis carolinensis (Reptilia: Lepidosauria: Iguanidae). Ital. J. Zool. 2013, 80, 177–186. [Google Scholar] [CrossRef]
  43. De la Lastra, J.M.P.; Asensio-Calavia, P.; González-Acosta, S.; Baca-González, V.; Morales-delaNuez, A. Bioinformatic Analysis of Genome-Predicted Bat Cathelicidins. Molecules 2021, 26, 1811. [Google Scholar] [CrossRef]
  44. Belov, K.; Sanderson, C.E.; Deakin, J.E.; Wong, E.S.; Assange, D.; McColl, K.A.; Gout, A.; de Bono, B.; Barrow, A.D.; Speed, T.P.; et al. Characterization of the opossum immune genome provides insights into the evolution of the mammalian immune system. Genome Res. 2007, 17, 982–991. [Google Scholar] [CrossRef] [PubMed]
  45. Choi, M.K.; Le, M.T.; Nguyen, D.T.; Choi, H.; Kim, W.; Kim, J.H.; Chun, J.; Hyeon, J.; Seo, K.; Park, C. Genome-level identification, gene expression, and comparative analysis of porcine β-defensin genes. BMC Genet. 2012, 13, 98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Kang, M.; Ahn, B.; Yum, J.; Cho, H.; Choi, M.; Hong, K.; Choi, Y.; Kim, J.; Park, C. Influence of habitat change from land to sea on the evolution of antimicrobial peptide gene families, including β-defensin gene clusters, in mammals. J. Zool. Syst. Evol. Res. 2021, 59, 510–521. [Google Scholar] [CrossRef]
  47. Cho, H.S.; Soundrarajan, N.; Le Van Chanh, Q.; Jeon, H.; Cha, S.Y.; Kang, M.; Ahn, B.; Hong, K.; Song, H.; Kim, J.H.; et al. The novel cathelicidin of naked mole rats, Hg-CATH, showed potent antimicrobial activity and low cytotoxicity. Gene 2018, 676, 164–170. [Google Scholar] [CrossRef]
  48. Ahn, B.; Jeon, H.; Cho H-s Nagasundarapandian, S.; Park, C. Sequence polymorphisms of PR39 cathelicidins and extensive copy variations in commercial pig breeds. Gene 2022, 822, 146323. [Google Scholar] [CrossRef]
  49. Zanetti, M. Cathelicidins, multifunctional peptides of the innate immunity. J. Leukoc. Biol. 2004, 75, 39–48. [Google Scholar] [CrossRef]
  50. Pérez de Lastra, J.M.; Garrido-Orduña, C.; Borges, L.; Tejera, D.; Marino, M.; González-Guerra, E.; Borges, A.A. Antimicrobial activity of cathelicidins of mammals from avian, aquatic and terrestrial environments. In Proceedings of the IV International Conference on Antimicrobial Research ICAR, Torremolinos, Málaga, Spain, 29 June–1 July 2016. [Google Scholar]
  51. Taylor, S.W.; Craig, A.G.; Fischer, W.H.; Park, M.; Lehrer, R.I. Styelin D, an Extensively Modified Antimicrobial Peptide from Ascidian Hemocytes. J. Biol. Chem. 2000, 275, 38417–38426. [Google Scholar] [CrossRef] [Green Version]
  52. Dennison, S.R.; Harris, F.; Bhatt, T.; Singh, J.; Phoenix, D.A. The effect of C-terminal amidation on the efficacy and selectivity of antimicrobial and anticancer peptides. Mol. Cell Biochem. 2009, 332, 43. [Google Scholar] [CrossRef] [Green Version]
  53. Simmaco, M.; Kreil, G.; Barra, D. Bombinins, antimicrobial peptides from Bombina species. Biochim. Biophys. Acta Biomembr. 2009, 1788, 1551–1555. [Google Scholar] [CrossRef] [Green Version]
  54. Kim, D.; Soundrarajan, N.; Lee, J.; Cho, H.-S.; Choi, M.; Cha, S.-Y.; Ahn, B.; Jeon, H.; Le, M.T.; Song, H.; et al. Genomewide Analysis of the Antimicrobial Peptides in Python bivittatus and Characterization of Cathelicidins with Potent Antimicrobial Activity and Low Cytotoxicity. Antimicrob. Agents Chemother. 2017, 61, e00530-17. [Google Scholar] [CrossRef] [Green Version]
  55. Nagant, C.; Pitts, B.; Nazmi, K.; Vandenbranden, M.; Bolscher, J.G.; Stewart, P.S.; Dehaye, J.-P. Identification of Peptides Derived from the Human Antimicrobial Peptide LL-37 Active against Biofilms Formed by Pseudomonas aeruginosa Using a Library of Truncated Fragments. Antimicrob. Agents Chemother. 2012, 56, 5698–5708. [Google Scholar] [CrossRef] [Green Version]
  56. Travis, S.M.; Anderson, N.N.; Forsyth, W.R.; Espiritu, C.; Conway, B.D.; Greenberg, E.P.; McCray, P.; Lehrer, R.I.; Welsh, M.J.; Tack, B.F. Bactericidal Activity of Mammalian Cathelicidin-Derived Peptides. Infect. Immun. 2000, 68, 2748–2755. [Google Scholar] [CrossRef] [Green Version]
  57. Van Dijk, A.; Molhoek, E.M.; Veldhuizen, E.J.A.; Bokhoven, J.L.M.T.v.; Wagendorp, E.; Bikker, F.; Haagsman, H.P. Identification of chicken cathelicidin-2 core elements involved in antibacterial and immunomodulatory activities. Mol. Immunol. 2009, 46, 2465–2473. [Google Scholar] [CrossRef]
  58. Hertog, A.L.D.; Marle, J.V.; Veerman, E.C.I.; Valentijn-Benz, M.; Nazmi, K.; Kalay, H.; Grün, C.H.; Van’t Hof, W.; Bolscher, J.G.; Nieuw Amerongen, A.V. The human cathelicidin peptide LL-37 and truncated variants induce segregation of lipids and proteins in the plasma membrane of Candida albicans. Biol. Chem. 2006, 387, 1495–1502. [Google Scholar] [CrossRef]
  59. Taheri, B.; Mohammadi, M.; Nabipour, I.; Momenzadeh, N.; Roozbehani, M. Identification of novel antimicrobial peptide from Asian sea bass (Lates calcarifer) by in silico and activity characterization. PLoS ONE 2018, 13, e0206578. [Google Scholar] [CrossRef]
  60. Peel, E.; Cheng, Y.; Djordjevic, J.T.; Fox, S.; Sorrell, T.C.; Belov, K. Cathelicidins in the Tasmanian devil (Sarcophilus harrisii). Sci. Rep. 2016, 6, 35019. [Google Scholar] [CrossRef] [Green Version]
  61. Soundrarajan, N.; Cho, H.-S.; Ahn, B.; Choi, M.; Thong, L.M.; Choi, H.; Cha, S.-Y.; Kim, J.-H.; Park, C.-K.; Seo, K.; et al. Green fluorescent protein as a scaffold for high efficiency production of functional bacteriotoxic proteins in Escherichia coli. Sci. Rep. 2016, 6, 20661. [Google Scholar] [CrossRef] [Green Version]
  62. Arqué, X.; Torres, M.D.T.; Patiño, T.; Boaro, A.; Sánchez, S.; de la Fuente-Nunez, C. Autonomous Treatment of Bacterial Infections in Vivo Using Antimicrobial Micro- and Nanomotors. ACS Nano 2022, 16, 7547–7558. [Google Scholar] [CrossRef]
  63. Brandelli, A. Nanostructures as Promising Tools for Delivery of Antimicrobial Peptides. Mini Rev. Med. Chem. 2012, 12, 731–741. [Google Scholar] [CrossRef]
  64. Piras, A.M.; Maisetta, G.; Sandreschi, S.; Gazzarri, M.; Bartoli, C.; Grassi, L.; Esin, S.; Chiellini, F.; Batoni, G. Chitosan nanoparticles loaded with the antimicrobial peptide temporin B exert a long-term antibacterial activity in vitro against clinical isolates of Staphylococcus epidermidis. Front. Microbiol. 2015, 6, 372. [Google Scholar] [CrossRef] [Green Version]
  65. Zanetti, M. The role of cathelicidins in the innate host defenses of mammals. Curr. Issues Mol. Biol. 2005, 7, 179–196. [Google Scholar]
  66. Jones, K.E.; Purvis, A.; MacLarnon, A.N.N.; Bininda-Emonds, O.R.P.; Simmons, N.B. A phylogenetic supertree of the bats (Mammalia: Chiroptera). Biol. Rev. 2002, 77, 223–259. [Google Scholar] [CrossRef]
  67. Barlow, P.G.; Svoboda, P.; Mackellar, A.; Nash, A.A.; York, I.A.; Pohl, J.; Davidson, D.J.; Donis, R.O. Antiviral Activity and Increased Host Defense against Influenza Infection Elicited by the Human Cathelicidin LL-37. PLoS ONE 2011, 6, e25333. [Google Scholar] [CrossRef]
  68. Zhang, H.; Zhao, Q.; Bhattacharya, S.; Waheed, A.A.; Tong, X.; Hong, A.; Heck, S.; Curreli, F.; Goger, M.; Cowburn, D.; et al. A Cell-penetrating Helical Peptide as a Potential HIV-1 Inhibitor. J. Mol. Biol. 2008, 378, 565–580. [Google Scholar] [CrossRef] [Green Version]
  69. Ben Mechlia, M.; Belaid, A.; Castel, G.; Jallet, C.; Mansfield, K.L.; Fooks, A.R.; Hani, K.; Tordo, N. Dermaseptins as potential antirabies compounds. Vaccine 2019, 37, 4694–4700. [Google Scholar] [CrossRef]
  70. Damour, A.; Garcia, M.; Cho, H.-S.; Larivière, A.; Lévêque, N.; Park, C.; Bodet, C. Characterisation of Antiviral Activity of Cathelicidins from Naked Mole Rat and Python bivittatus on Human Herpes Simplex Virus 1. Pharmaceuticals 2021, 14, 715. [Google Scholar] [CrossRef]
  71. Enayathullah, M.G.; Parekh, Y.; Banu, S.; Ram, S.; Nagaraj, R.; Kumar, B.K.; Idris, M.M. Gramicidin S and melittin: Potential anti-viral therapeutic peptides to treat SARS-CoV-2 infection. Sci. Rep. 2022, 12, 3446. [Google Scholar] [CrossRef] [PubMed]
  72. Roy, M.; Lebeau, L.; Chessa, C.; Damour, A.; Ladram, A.; Oury, B.; Boutolleau, D.; Bodet, C.; Lévêque, N. Comparison of Anti-Viral Activity of Frog Skin Anti-Microbial Peptides Temporin-Sha and [K3]SHa to LL-37 and Temporin-Tb against Herpes Simplex Virus Type 1. Viruses 2019, 11, 77. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Boutet, E.; Lieberherr, D.; Tognolli, M.; Schneider, M.; Bairoch, A. UniProtKB/Swiss-Prot. In Plant Bioinformatics: Methods and Protocols; Edwards, D., Ed.; Humana Press: Totowa, NJ, USA, 2007; pp. 89–112. [Google Scholar]
  74. Petersen, T.N.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods 2011, 8, 785–786. [Google Scholar] [CrossRef] [PubMed]
  75. Potter, S.C.; Luciani, A.; Eddy, S.R.; Park, Y.; Lopez, R.; Finn, R.D. HMMER web server: 2018 update. Nucleic. Acids Res. 2018, 46, W200–W204. [Google Scholar] [CrossRef] [Green Version]
  76. Torrent, M.; Nogués, V.M.; Boix, E. A theoretical approach to spot active regions in antimicrobial proteins. BMC Bioinform. 2009, 10, 373. [Google Scholar] [CrossRef] [Green Version]
  77. Pirtskhalava, M.; Amstrong, A.A.; Grigolava, M.; Chubinidze, M.; Alimbarashvili, E.; Vishnepolsky, B.; Gabrielian, A.; Rosenthal, A.; Hurt, D.E.; Tartakovsky, M. DBAASP v3: Database of antimicrobial/cytotoxic activity and structure of peptides as a resource for development of new therapeutics. Nucleic Acids Res. 2021, 49, D288–D297. [Google Scholar] [CrossRef]
  78. McGuffin, L.J.; Bryson, K.; Jones, D.T. The PSIPRED protein structure prediction server. Bioinformatics 2000, 16, 404–405. [Google Scholar] [CrossRef]
  79. Wang, G.; Li, X.; Wang, Z. APD3: The antimicrobial peptide database as a tool for research and education. Nucleic Acids Res. 2016, 44, D1087–D1093. [Google Scholar] [CrossRef] [Green Version]
  80. Cantón, E.; Espinel-Ingroff, A.; Pemán, J. Trends in antifungal susceptibility testing using CLSI reference and commercial methods. Expert Rev. Anti-Infect. Ther. 2009, 7, 107–119. [Google Scholar] [CrossRef]
  81. Ballard, E.; Yucel, R.; Melchers, W.J.G.; Brown, A.J.P.; Verweij, P.E.; Warris, A. Antifungal Activity of Antimicrobial Peptides and Proteins against Aspergillus fumigatus. J. Fungi 2020, 6, 65. [Google Scholar] [CrossRef]
Figure 1. Metabolic activity after the addition of bat cathelicidin-derived peptides against Mucor hiemalis KCTC 26779. ΔHA-CATH, ΔML-CATH, and ΔPD-CATH are the predicted antimicrobial activity core region peptides of Hipposideros armiger (accession number, XP_019486615.1), Myotis lucifugus (accession number, XP_006108360.2), and Phyllostomus discolor (accession number, XP_028374415.1) cathelicidins, respectively. HA-CATH, ML-CATH, and PD-CATH are the predicted nascent mature region peptides. Statistical significance was indicated above bars; a (p < 0.01) and b (p < 0.001).
Figure 1. Metabolic activity after the addition of bat cathelicidin-derived peptides against Mucor hiemalis KCTC 26779. ΔHA-CATH, ΔML-CATH, and ΔPD-CATH are the predicted antimicrobial activity core region peptides of Hipposideros armiger (accession number, XP_019486615.1), Myotis lucifugus (accession number, XP_006108360.2), and Phyllostomus discolor (accession number, XP_028374415.1) cathelicidins, respectively. HA-CATH, ML-CATH, and PD-CATH are the predicted nascent mature region peptides. Statistical significance was indicated above bars; a (p < 0.01) and b (p < 0.001).
Antibiotics 11 00989 g001
Table 1. Biochemical characteristics of bat cathelicidin-derived peptides.
Table 1. Biochemical characteristics of bat cathelicidin-derived peptides.
Bat Cathelicidin-Derived PeptidesSequenceLength<H> azb (+)Molecular Weight (Da)Similarity
AMP cSource(%) d
ΔHA-CATHLLRRGGRKIGQGLERIGQRIQGF233152609.08SMAP-29Ovis aries46.87
HA-CATHILGRLRDLLRRGGRKIGQGLERIGQRIQGFFSNREPMEES403044640.37K9CATHCanis familiaris51.22
ΔML-CATHGIFILKHRRPIGRGIEIT185032076.52Temporin-CPbLithobates capito40
ML-CATHLNPLIKAGIFILKHRRPIGRGIEITGRGIKKFFSK354083975.87Palustrin-2CG1Amolops chunganensis38.64
ΔPD-CATHIAGRIAGKLIGDAINRHRERNRQRR253662927.37TP4Oreochromis niloticus44.83
PD-CATHILGPALRIGGRIAGRIAGKLIGDAINRHRERNRQRRG373584088.79HKPLPHippocampus kuda43.59
a Ratio of hydrophobic residues. b Charge. c Most similar peptides among known antimicrobial peptides (AMPs). d Percentage of identical residues, including gaps, in the globally aligned sequence.
Table 2. Minimum inhibitory concentrations of three bat cathelicidin-derived peptides against various bacterial strains.
Table 2. Minimum inhibitory concentrations of three bat cathelicidin-derived peptides against various bacterial strains.
StrainMinimum Inhibitory Concentration (μg/mL, μM)
HA-CATHML-CATHPD-CATHchloramphenicolAmpicillinGentamicin
Gram-negative bacteriaEscherichia coli ATCC 2592218 (3.9)2 (0.5)7 (1.7)3 (9.3)5 (14.3)1 (2.1)
Pseudomonas aeruginosa ATCC 27853>40 (8.6)>40 (10.1)>40 (9.8)80 (247.6)>640 (1831.7)1 (2.1)
Salmonella enterica serovar Typhimurium ATCC 14028>40 (8.6)21 (5.3)22 (5.4)5 (15.5)>80 (228.8)1 (2.1)
Acinetobacter baumannii KCTC 232545 (1.1)4 (1.0)4 (1.0)38 (117.6)>80 (228.8)10 (21.0)
Klebsiella pneumoniae subsp. pneumoniae KCTC 172634 (7.3)12 (3.0)13 (3.2)> 80 (247.6)>80 (228.8)19 (39.80
Enterobacter cloacae subsp. cloacae ATCC 13047>40 (8.6)35 (8.8)30 (7.3)5 (15.5)>80 (228.8)1 (2.1)
Gram-positive bacteriaStaphylococcus aureus ATCC 653826 (5.6)1 (0.3)3 (0.7)7.5 (23.2)2 (5.7)1 (2.1)
Bacillus cereus ATCC 1087625 (5.4)3 (0.8)6 (1.5)10 (31.0)80 (228.8)1 (2.1)
Enterococcus faecalis ATCC 29212>40 (8.6)5 (1.3)6 (1.5)10 (31.0)10 (28.6)90 (189.0)
Streptococcus agalactiae ATCC 27956>40 (8.6)8 (2.0)8 (2.0)5 (15.5)4 (11.4)75 (157.5)
Streptococcus dysgalactiae ATCC 27957>40 (8.6)14 (3.5)18 (4.4)4 (12.4)2 (5.7)15 (31.5)
Streptococcus equi subsp. zooepidemicus ATCC 43079>40 (8.6)9 (2.3)17 (4.2)5 (15.5)2 (5.7)45 (94.5)
Table 3. Minimum inhibitory concentrations of six bat cathelicidin-derived peptides against two yeast form fungal strains.
Table 3. Minimum inhibitory concentrations of six bat cathelicidin-derived peptides against two yeast form fungal strains.
StrainMIC (μg/mL, μM)
ΔHA-CATHHA-CATHΔML-CATHML-CATHΔPD-CATHPD-CATHCiclopirox a
Candida albicans
KCTC 7270
>40 (15.3)>40 (8.6)>40 (19.3)40 (10.1)25 (8.5)5 (1.2)3.5 (16. 9)
Cryptococcus cuniculi
KCTC 17232
>40 (15.3)45 (9.7)>40 (19.3)5 (1.3)15 (5.1)5 (1.2)0.5 (2.4)
a Reference for antifungal activity.
Table 4. Viability of HaCaT cells after the addition of bat cathelicidin-derived peptides.
Table 4. Viability of HaCaT cells after the addition of bat cathelicidin-derived peptides.
TreatmentConcentration (μg/mL)Cell Viability ± SD (%)
ΔHA-CATH64100.5 ± 0.5
16097.7 ± 4.3
HA-CATH64100 ± 0.3
16097.5 ± 0.6
ΔML-CATH64100.6 ± 1.3
16096.0 ± 0.8
ML-CATH6499.7 ± 0.8
16097 ± 1.8
ΔPD-CATH64100.3 ± 0.3
16098.4 ± 2.3
PD-CATH64100.5 ± 0.2
16010.4 ± 0.4 *
Melittin647.2 ± 0.03 *
Triton X-100-8.9 ± 0.05 *
Negative control-99.5 ± 0.83
Note: Statistical significance; * p < 0.001.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Choi, M.; Cho, H.-s.; Ahn, B.; Prathap, S.; Nagasundarapandian, S.; Park, C. Genomewide Analysis and Biological Characterization of Cathelicidins with Potent Antimicrobial Activity and Low Cytotoxicity from Three Bat Species. Antibiotics 2022, 11, 989. https://doi.org/10.3390/antibiotics11080989

AMA Style

Choi M, Cho H-s, Ahn B, Prathap S, Nagasundarapandian S, Park C. Genomewide Analysis and Biological Characterization of Cathelicidins with Potent Antimicrobial Activity and Low Cytotoxicity from Three Bat Species. Antibiotics. 2022; 11(8):989. https://doi.org/10.3390/antibiotics11080989

Chicago/Turabian Style

Choi, Munjeong, Hye-sun Cho, Byeongyong Ahn, Somasundaram Prathap, Soundrarajan Nagasundarapandian, and Chankyu Park. 2022. "Genomewide Analysis and Biological Characterization of Cathelicidins with Potent Antimicrobial Activity and Low Cytotoxicity from Three Bat Species" Antibiotics 11, no. 8: 989. https://doi.org/10.3390/antibiotics11080989

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop