Full text
51,698 characters
· extracted from
preprint-html
· click to expand
The Emergence of Antimicrobial and Phage Resistance in Piscine Lactococcus Isolates and Implications for Future Disease Management in Aquaculture. | Authorea try { document.documentElement.classList.add('js'); } catch (e) { } var _gaq = _gaq || []; _gaq.push(['_setAccount', 'G-8VDV14Y67G']); _gaq.push(['_trackPageview']); (function() { var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true; ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js'; var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s); })(); Skip to main content Preprints Collections Wiley Open Research IET Open Research Ecological Society of Japan All Collections About About Authorea FAQs Contact Us Quick Search anywhere Search for preprint articles, keywords, etc. Search Search ADVANCED SEARCH SCROLL This is a preprint and has not been peer reviewed. Data may be preliminary. 3 June 2025 V1 Latest version Share on The Emergence of Antimicrobial and Phage Resistance in Piscine Lactococcus Isolates and Implications for Future Disease Management in Aquaculture. Authors : Adam M. Blanchard 0000-0001-6991-7210 [email protected] , Bailey Secker , Samantha J. Windle , Robert Atterbury , Ha Thanh Dong , Le Thanh Dien , David Huchzermeyer , Bernard Hang'ombe , and Saengchan Senapin Authors Info & Affiliations https://doi.org/10.22541/au.174894333.36446684/v1 281 views 155 downloads Contents Abstract Information & Authors Metrics & Citations View Options References Figures Tables Media Share Abstract Lactococcosis is a major bacterial disease impacting rainbow trout production in South Africa and Southeast Asia, particularly during the summer. This study analysed 15 bacterial isolates from affected aquaculture facilities, revealing that Lactococcus petauri (n=12) was the predominant species, rather than the traditionally recognised L. garvieae (n=3). This finding suggests a possible shift in the aetiology of lactococcosis, with potential implications for disease management. Genotypic analysis revealed that only L. garvieae isolates carried the adhesin gene (adh), which is critical for bacterial adhesion and colonisation. Most isolates possessed sortase-anchored proteins linked to iron uptake, adhesion, and stress resistance, with the LPxTG-6 subgroup unique to L. garvieae . Variations in the capsule operon, including transposase insertions, suggest ongoing horizontal gene transfer, possibly influencing immune evasion. Antimicrobial resistance (AMR) analysis identified efflux pumps ( mdtA, lsaD ) conferring resistance to macrolides and lincosamides, with additional genotypic resistance to erythromycin ( mefA, msrD ) and tetracyclines ( tetS, tetL ). The presence of plasmid-borne tetL raises concerns about potential gene transfer and the persistence of resistance in aquaculture systems. Prophage elements were identified in several isolates, potentially contributing to virulence, and immune modulation. Although phage therapy has shown promise in experimental settings, the presence of viral defence mechanisms may raise a unique challenge. These findings highlight the importance of monitoring pathogen evolution in aquaculture systems and suggest that ongoing genomic surveillance and treatment strategies may need to be adapted to account for the emerging role of L. petauri in lactococcosis outbreaks. Introduction Fish aquaculture is a rapidly growing sector of global food production, contributing significantly to food security, economic development, and livelihoods worldwide 1 . Among the diverse species farmed, economically important fish such as tilapia ( Oreochromis spp.), rainbow trout ( Oncorhynchus mykiss ), and Asian seabass ( Lates calcarifer ) have gained prominence due to their high market demand, nutritional value, and adaptability to various aquaculture systems 2–4 . Tilapia is a staple in both developing and developed countries, prized for its affordability and rapid growth rate. Rainbow trout is highly valued for its premium quality, appealing taste, and contribution to high-value markets. Similarly, Asian seabass is a highly valued species in the aquaculture sector, known for its versatility in farming environments, excellent growth performance, and economic value in both domestic and export markets. The increasing reliance on fish aquaculture to meet the protein needs of a growing human population highlights its importance. However, the sustainability and productivity of the industry are increasingly threatened by bacterial infections 5–7 , which result in significant economic losses and pose challenges to maintaining fish health and welfare. There are several common bacterial diseases affecting farmed fish, including streptococcosis, aeromonasis, columnaris, edwardsiellosis and lactococcosis. Among these, lactococcosis - caused by lactococcosis-causing bacteria (LCB) - is one of the emerging concerns in Africa and Southeast Asia. This disease affects a range of farmed fish, including tilapia and rainbow trout, resulting in economic losses for aquaculture operations 8 . The most common pathogen responsible for lactococcosis is Lactococcus garvieae , but other species such as Lactococcus petauri and Lactococcus formosensis have also been recently reported to contribute to the disease. Lactococcosis is typically characterized by clinical signs such as exophthalmos, acute haemorrhagic septicaemia, with mortality rates between 50-80% 9 . Lactococcosis is classed as an emerging global zoonotic pathogen with cases in marine and terrestrial mammals 10,11 , reptiles 12 , birds 13 , and also rarely in humans 14 . Outbreaks have occurred globally in aquaculture systems 15,16 which have a significant economic and welfare impact to the industry 17 . Outbreaks of disease are normally attributed to L. garvieae . However, classical microbiological identification techniques 18,19 often fail to differentiate this species from related aetiological agents such as L. petauri and L. formosensis. Additionally, pathogenesis is poorly defined and there is a lack of knowledge around virulence determinants in the three species, as most research is focused on clinical isolates from aquaculture 20 . Moreover, recent comparative genomics studies 21,22 have identified other homologous streptococcal virulence genes associated with adhesion 23 , iron transport 24 and capsular polysaccharide 25 . Encapsulation has long been considered a virulence factor in many streptococcal species as it masks the cell surface antigens from the immune system 26,27 . The capsule is immunogenic, and a target for antimicrobial drugs in many bacterial species. Shedding of the capsule on exposure to these compounds can be used as an evasion technique to evade autophagy-mediated killing in macrophages 28 . By comparing pathogenic and non-pathogenic strains of piscine L. garvieae the presence of a capsule has been shown to be crucial for virulence in fish 22 , however, it is sometimes not present, meaning encapsulation is not the sole requirement for virulence 29 . Adhesins are cell-surface components that facilitate the binding of bacteria to host cell receptors. L. garvieae has been found to produce three proteins (PsaA, PavA, and enolase) with high similarity to adhesion virulence factors in other streptococcal species, as well as adhesin clusters 1 and 2 ( adhCI and adhCII ) and adhesin ( adh ) 23,29 . Sortase- anchored proteins contribute to the establishment and persistence of disease, and as such are likely to be virulence factors 30 . Antibiotics have long been deployed to control Lactococcus spp and Streptococcus spp. generally in aquaculture. Most commonly these include erythromycin, oxytetracycline, amoxicillin and low-level doxycycline 31 . However, resistance to erythromycin, oxytetracycline and lincomycin is growing due to the presence of ermB and tetS 32 . Outbreaks are becoming more frequent and control through medicated feed often fails due to a combination of inappetence and antimicrobial resistance (AMR) 15 . This is the first study to date which has analysed Lactococcal genomes isolated from outbreaks of lactococcosis at aquaculture facilities in Zambia, South Africa, Vietnam and Thailand. We have provided a deeper understanding of antimicrobial resistance, virulence factors, abundance of phage resistance mechanisms and the implications for the use of phage therapy in future disease management. Bacterial isolates Lactococcus spp. isolates from Thailand (n=3), Vietnam (n=5), Zambia (n=5), and South Africa (n=2), previously obtained from diseased fish species including tilapia, Asian seabass, and rainbow trout were employed for whole genome sequencing in this study (Table 1). These isolates were previously identified as L. garvieae through 16S rDNA sequencing. Virulence test results from experimental infections using some of these isolates are also presented in Table 1. not-yet-known not-yet-known not-yet-known unknown DNA Extraction DNA extraction of Lactococcus spp. was performed following previously optimised methods33 with some modifications as follows. A single colony was cultured overnight at 30°C in 1.5 mL of tryptic soy broth (Difco). Cells were harvested by centrifugation at 9,000 × g for 5 minutes, followed by washing with 1 mL of distilled water and subsequent collection by centrifugation. The cell pellet was resuspended in 300 µL of EDTA-saline solution (0.01 M EDTA, 0.15 M NaCl, pH 8.0) and thoroughly vortexed. Then, 50 µL of 110 mg/mL lysozyme (SERVA) and 10 µL of 20 mg/mL RNase A (Thermo Fisher Scientific) were added, mixed by vortexing, and incubated at 37°C for 2 hours, with vortexing every 15 minutes. After incubation, 160 µL of 20% SDS and 20 µL of 5 mg/ml proteinase K (Merck) were added, followed by vortexing and incubation at 65°C for 30 minutes. Next, 0.5 volume of 5 M NaCl was added and briefly mixed. An equal volume of phenol (1:1, v/v) was added, vortexed, and centrifuged at 12,000 × g for 15 minutes. The upper aqueous phase was transferred to a new tube, and this extraction step was repeated once more. Then, 2 µL of 20 mg/mL RNase A was added and incubated for 15 minutes, followed by a final phenol extraction. The upper phase was transferred to a new tube. To precipitate the DNA, 0.1 volume of 3 M sodium acetate (pH 5.2) was added and mixed gently, followed by the addition of 2 volumes of cold absolute ethanol. The mixture was incubated at -20°C for 2 hours and then centrifuged for 10 minutes to collect the DNA pellet. The DNA pellet was washed with cold 70% ethanol, air-dried, and resuspended in 20 µl of RNase- and DNase-free water. DNA quality and quantity were measured using a Nanodrop spectrophotometer. Table 1: Provenance of the strains used in the study Thailand Oct 22 Kanchanaburi L. garvieae Red Tilapia 3749 No 34 Thailand May 22 Chachoengsao L. petauri Asian sea bass 3750 No 34 Thailand May 22 Chachoengsao L. petauri Asian sea bass 3752 N/A Zambia Mar 17 Lake Kariba L. petauri Nile tilapia 3816 N/A Zambia Oct 17 Lake Kariba L. petauri Nile tilapia 3817 N/A Zambia Feb 18 Lake Kariba L. petauri Nile tilapia 3818 N/A Zambia Oct 18 Lake Kariba L. petauri Nile tilapia 3819 N/A Zambia Nov 18 Lake Kariba L. petauri Nile tilapia 3820 N/A South Africa Nov 15 Unknown L. petauri Rainbow trout 3821 N/A South Africa May 16 Unknown L. petauri Rainbow trout 3822 N/A Vietnam Aug 22 Ben Tre L. petauri Red Tilapia 3830 Low Vietnam Aug 22 Ben Tre L. garvieae Red Tilapia 3831 Medium Vietnam Aug 22 Ben Tre L. garvieae Red Tilapia 3832 High Vietnam Jan 22 Tien Giang L. petauri Red Tilapia 3833 Medium Vietnam Jan 22 Tien Giang L. petauri Red Tilapia 3834 Low *Species level taxonomic assignment of the 15 Lactococcus genomes is based on digital DNA:DNA hybridisation determined in this study. Sequencing Bacterial DNA samples (1.0 µg each) were sent for whole genome sequencing to BGI Hong Kong Tech Solution NGS Lab via Bangkok Genomics Innovation Co., Ltd. (Thailand). DNA quantity and quality were verified using microplate reading and gel electrophoresis before library preparation. The DNA libraries were then prepared using the short-insert library method and sequenced with DNBseq technology at a PE150 read length generating 4 million reads per sample. Data filtering, including the removal of adaptor sequences and low-quality reads, was performed by the company. Data Analysis Assemblies were created using NF-core Bacass 35 v2.2.0 under short read mode. Briefly, the reads are assessed for quality using FastQC v0.12.0 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and trimmed using FastP 36 v0.23.4. Clean and high quality reads are assembled using Unicycler 37 v0.5.0 and the resulting contigs were checked and filtered for taxonomic assignment using Kraken2 38 v2.1.3. Genome assemblies have been deposited in National Library of Medicine (NCBI) GenBank under the BioProject PRJNA1264975. Taxonomic classification of genomes was calculated using digital DNA:DNA hybridisation with GGDC 39 using reference genomes L. garvieae GCA_016026695, L. petauri GCA_023499275 and L. formosensis GCA_037892375. Relatedness of bacterial genomes was tested using Sourmash 40 in average nucleotide identity (ANI) mode and heatmaps were created in R v4.3.2 41 using Pheatmap (https://github.com/raivokolde/pheatmap). The genomes were parsed for genotypic antimicrobial resistance using the AMR++ and the MegaRes3 database 42 . Virulence gene presence was determined using abricate with a custom database of L. garvieae specific virulence genes as described by Lin et al 21 . Presence of prophage regions was determined using three tools; geNomad v1.8.0 43 , PhageBoost v0.1.7 44 and VIBRANT v1.2.1 45 . The predictions were then consolidated to identify the longest possible region using Biopython 46 . The resulting sequences were processed using CheckV v1.0.3 47 , retaining only sequences identified as ‘medium quality’ or higher. Additionally, where host contamination was identified by CheckV, the trimmed sequence was kept. Sequences from the same host were dereplicated using CheckV supporting code. Finally, the sequences were annotated using Pharokka 48 v1.7.5 followed by phold v0.2.0 (https://github.com/gbouras13/phold). Taxonomic assignment was completed using taxmyPHAGE 49 and proteomic trees were generated using Viral Proteomic Tree (ViPTree). Genomic level comparisons of the prophage were performed using CheckV supporting ANI and clustering code, default clustering cutoffs of 95% ANI over 85% alignment fraction were used as described by Roux et al . 50 . Genomic organisation and protein clustering was performed using LoVis4u v0.1.1 51 , loci were reorientated, and functional annotation was performed for defence, anti-defence, virulence and antimicrobial resistance genes. Prophage were considered ‘intact’ if they contained an integrase, structural proteins and lysis modules in the conserved gene order as assessed by manual curation 52,53 . Finally, the termini of prophage regions of interest were investigated using PhageTermVirome v4.1 54 . Results Genome Assembly The quality control data achieved a mean Phred quality of 39 and de novo assembled following the Bacass 35 short read pipeline to approximately 2.0 Mbp (min 2.0 Mbp max 2.3 Mbp) with a mean N50 of 756 Kbp. Based on Kraken2 38 the final assemblies where classified as Lactococcus species. Further taxonomic assignment was generated through dDDH which yielded a species level classification with probability scores ≥93% to reference species L. garvieae GCA_016026695 and L. petauri GCA_023499275, corresponding to 3 and 12 isolates, respectively (Table 1). Interestingly, all of these isolates had previously been identified as L. garvieae through 16S rDNA sequencing. Genome Comparison Based on average nucleotide identity scores (Figure 1), the isolates seem to form two clusters based on species classification. Cluster 1 is the L. garvieae isolates, where clusters 2, 3 and 4 are all L. petauri where they split into country of origin. Cluster 2 contains isolates from Zambia and South Africa whereas cluster 3 is Vietnamese isolates and cluster 4 is Thai isolates. Figure 1. Average Nucleotide Identify comparison of all isolates. The nucleotide identity score is between 0.94 and 1, representing 94% and 100% similarity between all isolates. The Y axis is coloured by country of origin, and clustered using a Euclidean approach, and the X axis shows host species. Antimicrobial resistance genotype All isolates showed multi-drug genotypic resistance with the presence of mdtA (drug, biocide and metal RND efflux pumps) and lsaD (multi-drug ABC efflux pumps). Two isolates (3833 and 3834) had mefA (MLS erythromycin resistance MFS efflux pumps) and msrD (MLS erythromycin resistance ABC efflux pumps). Only two isolates (3821 and 3822) had plasmid borne tetracycline genotypic resistance with the presence of tetL (accession NF012176.0) and three isolates (3750, 3821 and 3822) had tetS. Also, two isolates (3833 and 3834) had lnuC (lincosamide) genotypic resistance (Supp Table 1). Overall, the isolates could be ascribed to 4 resistance profiles, the most common being MLS efflux pumps (Table 2). Table 2. Antimicrobial Resistance Profiles for the isolates 3749 • • L. garvieae Thailand 3750 • • • L. petauri Thailand 3752 • • L. petauri Thailand 3816 • • L. petauri Zambia 3817 • • L. petauri Zambia 3818 • • L. petauri Zambia 3819 • • L. petauri Zambia 3820 • • L. petauri Zambia 3821 • • • • L. petauri S Africa 3822 • • • • L. petauri S Africa 3830 • • L. petauri Vietnam 3831 • • L. garvieae Vietnam 3832 • • L. garvieae Vietnam 3833 • • • • • L. petauri Vietnam 3834 • • • • • L. petauri Vietnam not-yet-known not-yet-known not-yet-known unknown Gene names and linked function: mdtA Drug and biocide and metal RND efflux pumps, lsaD Multi-drug ABC efflux pumps, lnuC Lincosamide nucleotidyltransferases, mefA MLS erythromycin resistance MFS efflux pumps, msrD MLS erythromycin resistance MFS efflux pumps, tetL Tetracycline resistance MFS efflux pumps, and tetS Tetracycline resistance ribosomal protection proteins. Full data in Supplementary Table 1. Virulence genotype Based on a cutoff of ≥95 identity, all isolates had the three haemolysin genes, the nine iron uptake genes and adhesion genes (enolase, pavA and psaA ) but only three isolates (3749, 3831 and 3832) had adhesin. All isolates had GAPDH, NADH oxidase, phosphoglucomutase, superoxide dismutase and collagenase genes. Only 2 isolates had a complete capsule operon of all 12 loci (3830, 3832), while 2 isolates 3833 and 3834 had 11 genes from the operon, each missing the 1p transposase. All isolates had srtA and LPxTG-1, and nine isolates (3861, 3817, 3818, 3819, 3820, 3821, 3822, 3831 and 3832 had LPxTG-3 and LPxTG-6 while 3749 only had LPxTG-6. Only 2 isolates (3821 and 3832) had LPxTG-2 and LPxTG-4 (Supplementary Table 2). Prophage A total of 144 putative prophage regions were identified in the isolates, of these, 22 were retained after filtering with CheckV. Only two (2/15, 13%) of the Lactococcus genomes did not contain any prophage regions, and between 1-3 prophage regions were identified in the remaining genomes (13/15, 87%). The prophage regions were between 13-69 kbp in length and contained between 29-129 coding sequences. taxmyPHAGE assigned all the prophage to novel genera, however, 3830 prophage 5, 3833 prophage 1, and 3834 prophage 1 were taxonomically similar to Siphovirus -like phage genus Piorkowskivirus . Comparison of the regions using a proteomic tree did not identify clustering based on host species or country of origin (Figure 2). Additionally, when viewed on a proteomic tree with 3,399 sequences from the ViPTree viral database, the prophage were not closely related. However, they did all fall within the area of the tree with viruses of Bacillota (previously Firmicutes ; Supplementary Figure 1). The majority (59%;13/22) of the prophage were similar to other Lactococcus phage, the remaining regions clustered with phage that infect Thermoanaerobacterium , Streptococcus , Enterococcus , and Lactobacillus . Figure 2. Prophage proteomic phylogenetic tree. VipTree was used to generate a proteomic tree of the prophage regions identified as part of the present study and visualised using iToL. The scale represents 1 base difference per 10 bases. The colour bars represent the country of origin and host species. As the prophage regions could not be taxonomically classified, the genomic diversity of the prophage was investigated using pairwise ANI, followed by clustering using thresholds of 95% ANI over 85% alignment fraction, this identified 11 putative virus operational taxonomic units (vOTUs) (Supplementary Figure 2). Similarly, when the annotated prophage were visualised with LoVis4u, protein clustering supported 11 vOTUs (Supplementary Figure 3). This revealed that the proteins within prophage regions were highly conserved within vOTUs. This is particularly interesting for prophage identified in the genomes of L. petauri 3816-20 as these samples were collected from infected Nile tilapia found in Lake Kariba across a period of almost two years suggesting the persistence of L. petauri within the tilapia population over this period. When the prophage regions were interrogated for the presence of antimicrobial resistance and virulence determinants, a number of genes were identified. However, despite the comprehensive methods used to identify prophage regions, only 45% (10/22) were found to contain protein sequences indicative of an intact prophage region. No virulence or AMR genes were identified within an intact prophage region (Figure 3). Despite this, several defence genes were identified associated with intact prophage regions, this would be an important consideration were bacteriophage treatment investigated as an alternative to antibiotics. Specifically, a coding sequence identified as AbiD was found in two prophage regions. Abi systems have been previously reported on lactococcal plasmids 55 . To confirm that AbiD was within the prophage region, PhageTermVirome was used to identify the termini of prophage 3822 prophage 3 and 3821 prophage 1. The termini of 3822 prophage 3 could not be determined however the termini of 3821 prophage 1 was identified and included the coding sequence of AbiD. Restriction-modification (R-M) systems are ubiquitous among bacterial chromosomes and plasmids, these function by degrading DNA while protecting self-DNA with modification 56 . Coding sequences identified as type II DNA methyltransferases and endonucleases were present in multiple prophage regions. Although involved in other phage processes, these are thought to protect the prophage from host R-M systems during the lytic phase 57 . Figure 3. Genomic organisation of the identified intact prophage regions. The genomes are reorientated and homology between coding sequences is shown. Coding sequences are coloured based on presumed function. Discussion This study characterised the genomes of 15 Lactococcus strains isolated from piscine disease outbreaks in Asia and Africa. This revealed the presence of up to five antimicrobial resistance genes and three novel prophage regions per strain. Virulence gene characterisation found that only L. garvieae isolates carried the adhesin gene adh which is critical for bacterial colonisation. Although the prophage regions did not harbour any known virulence or antimicrobial resistance genes, they did contain several genes encoding phage defence systems such as (AbiD and R-M systems). This information is important when considering the potential development of phage therapy for piscine lactococcosis. Interestingly, the isolates from cases of lactococcosis in the present study were not all L. garvieae (n=3), the most commonly recognised aetiological cause of the disease, but a majority were L. petauri (n=12). The predominant detection of L. petauri in this study, supported by recent reports 18,58 , suggests that this species may have been either previously misidentified as L. garvieae , due to their high phenotypic and 16S rRNA gene sequence similarity, or previously overlooked in past surveillance efforts. These findings underscore the importance of implementing both active and retrospective surveillance strategies for Lactococcus spp. in fish populations to improve species-level identification and enhance disease management strategies. Only three isolates (3749, 3831 and 3832) had the adhesin gene ( adh ), which seems specific to the L. garvieae isolates. Adhesins are necessary for the ability of the bacteria to bind to the host cell surface and are required for colonisation and pathogenesis. Despite the essentiality of these loci, it has been documented that some strains of L. garvieae do not possess all three loci adh, pav and psaA 29 and these variations may hint to host specificity 59 . Sortase Proteins linked to Increased Virulence Most of the isolates investigated possessed multiple sortase-anchored proteins, however, the low virulence isolate, 3749 only had LPxTG-6 and isolates 3821 and the high virulence isolate 3832 had LPxTG-2 and LPxTG-4. These sortase anchored proteins in particular are highly associated with the iron uptake operon, adhesion and haemolysis 21 . This may aid in colonisation and stress resistance, enabling adaptation and survival in nutrient-poor or hostile environments. The LPxTG-6 subgroup is unique to isolates 3749, 3821 and 3822 which came from two different locations. There is a dearth of information on the specific functional implications of different sortase anchored protein variants, but due to the high co-occurrence of LPxTG-6 to the iron uptake operon, function could be suggested to play an important role in nutrient uptake, whereas LPxTG-4 seems to be more associated with psaA, which could hint at a roll in initial colonisation 21 . Only isolates 3833 and 3834 had an incomplete complement of 12 genes from the capsule operon, each missing the 1p transposase. The feature of genes encoding transposases or phage-related proteins within a capsule operon may suggest evidence of transposition or horizontal gene transfer. It is likely they alter production of the capsule related proteins rather than complete disruption, serving as a method of immune evasion. This has been documented in other Gram positive bacteria, where the transposase loci in Streptococcus pneumoniae differentiates between serotypes in serogroup 12 60 . Development of Antimicrobial Resistance to Commonly used Antibiotics All isolates processed the efflux pumps mdtA which are associated with macrolide resistance , and lsaD associated with lincosamide resistance. These resistances are commonly found generally in Gram positive bacteria and would not be of cause for concern in first line treatment of infections. However, two isolates showed genotypic resistance to erythromycin ( mefA and msrD ). While not uncommon in Gram positive bacteria, erythromycin has been used in prophylactic feed based medication in aquaculture, which has driven the selective pressure for resistant strains of Lactococcus 61 . Finally, three isolates showed genotypic resistance to tetracyclines with the presence of tetS and two additionally showed plasmid-borne tetracycline resistance from tetL . This resistance is also prevalent in isolates from aquacultural systems, and has been detected more frequently since the increased use of tetracycline for treatment of lactococcosis from 2014 62 . There might be a greater concern from the plasmid-borne nature of tetS and the potential of its ability to transfer if stock is imported from different brood stock providers. Identification of Phage Defence Mechanisms The rapid development of antimicrobial resistance has focused attention on alternative treatments such as bacteriophage therapy, bacterial viruses that infect and replicate within bacterial cells, as an alternative to antibiotics. Only one study has attempted phage therapy against L. garvieae infections in yellowtail fish, which reported significant reductions in mortality in fish treated with phage either injected peritoneally, or through feed 63 . However, timing of phage delivery was critical, with mortality appreciably higher when treatment was delayed to 24 h after bacterial challenge. Despite promising results, there are some limitations of bacteriophage therapy, namely the development of resistance which has been discussed in detail elsewhere (https://doi.org/10.1038/nrmicro2315). One mechanism by which bacteria can be resistant to bacteriophage infection is related to prophage, temperate bacterial viruses which are incorporated in the bacterial genome which can confer phenotypic advantages in environmental niches. These advantages may include increased resistance to acid or antimicrobials, enhanced biofilm production, higher growth rates, and modulation of eukaryotic host immune response 64 . Prophage-encoded proteins of Streptococcus mitis 65 and Enterococcus faecalis 66 facilitate adhesion to human platelets which is crucial for the development of endocarditis, and may be relevant to L. garvieae pathophysiology. Temperate phage are also associated with horizontal gene transfer, via transduction, and conferring superinfection immunity to closely related phage. In the present study, AbiD was identified within the one of the prophage regions. Abortive infection (Abi) systems are diverse in Lactococcus lactis and generally function by causing death of the bacterial host 67 , consequently, the presence of Abi systems may have unexpected effects if bacteriophage were used to treat these infections. Typically, for these reasons, temperate phage are considered unsuitable for use as direct therapeutic agents. However, in cases where strictly-lytic phage are rare or impossible to find, prophage can be used as the basis for developing effective treatments 68 . Most commonly, this includes identifying sequences encoding hydrolytic enzymes (usually endolysins) phage use to escape from infected cells after replication, then expressing and purifying the proteins. Endolysins are particularly effective against Gram positive bacterial pathogens which lack an outer membrane to defend their cell wall from digestion 69 . Lytic and temperate phage which infect L. garvieae have been isolated 42,43 . However, only minimal data is available on the lytic spectrum of these phage – which indicates a very limited host range - so it is difficult to determine if they would be useful against the broad diversity of L. garvieae strains causing disease. Although tilapia is native to Africa, it has been domesticated and cultivated in over 140 countries worldwide. Seed production and aquaculture technologies have advanced significantly in Asia. African countries have a long history of importing tilapia seed from Asian nations. Movement of live fish, which may carry pathogens through subclinical infections for aquaculture, could contribute to the translocation of these pathogens, similar to the dynamics observed with viruses such as TiLV and ISKNV 73,74 . The genomic analysis in this study suggest that the African isolates and Asian isolates form different unique clusters, suggesting their diversification in genetic evolution. Nonetheless, this study provides some fundamental knowledge of genomic characterization of the piscine Lactococcus spp., which is useful for further study on distribution and tracking genomic epidemiology as reported earlier with other pathogens 75 . In summary, the predominance of L. petauri among isolates from lactococcosis outbreaks in fish farms in South Africa and Southeast Asia indicates a potential shift in the etiological landscape of this disease. Comparative genomic analyses revealed distinct differences between L. petauri and L. garvieae in terms of virulence-associated factors, antimicrobial resistance genes, and prophage content. These findings underscore the need for enhanced taxonomic resolution in diagnostic protocols and reinforce the importance of ongoing genomic surveillance. Adaptation of current disease management and treatment strategies is warranted to address the emerging role of L. petauri in aquaculture, particularly in light of its potential for antimicrobial resistance and horizontal gene transfer. Acknowledgements This research was supported by the Budget Bureau and the National Science and Technology Development Agency. We are grateful to Dr Evelien Adriaenssens and Dr Ryan Cook for their helpful discussions regarding analysis of prophage in this study. The authors would like to thank Janchai Wongkaew for her technical assistance. References 1. The State of World Fisheries and Aquaculture 2024 . (FAO, 2024). doi:10.4060/cd0683en.2. Biology and Culture of Asian Seabass Lates Calcarifer . (CRC Press, 2013). doi:10.1201/b15974.3. D’Agaro, E., Gibertoni, P. & Esposito, S. Recent Trends and Economic Aspects in the Rainbow Trout (Oncorhynchus mykiss) Sector. Applied Sciences 12 , 8773 (2022).4. Pullin, R. S. V. An introduction to the tilapias. in Tilapia (eds. Thompson, K., Shoemaker, C. A. & Little, D. C.) 1–25 (CABI, GB, 2025). doi:10.1079/9781800629455.0001.5. Lan, N. G. T. et al. Review of current perspectives and future outlook on bacterial disease prevention through vaccination in Asian seabass ( Lates calcarifer ). Journal of Fish Diseases 47 , e13964 (2024).6. Duman, M., Altun, S., Saticioglu, I. B. & Romalde, J. L. A review of bacterial disease outbreaks in rainbow trout ( Oncorhynchus mykiss ) reported from 2010 to 2022. Journal of Fish Diseases jfd.13886 (2023) doi:10.1111/jfd.13886.7. Soto, E. et al. Tilapia bacterial diseases. in Tilapia (eds. Thompson, K., Shoemaker, C. A. & Little, D. C.) 378–452 (CABI, GB, 2025). doi:10.1079/9781800629455.0009.8. Heckman, T. I. et al. Redefining piscine lactococcosis. Appl Environ Microbiol 90 , e02349-23 (2024).9. Abu‐Elala, N. M., Abd‐Elsalam, R. M. & Younis, N. A. Streptococcosis, Lactococcosis and Enterococcosis are potential threats facing cultured Nile tilapia ( Oreochomis niloticus ) production. Aquac Res 51 , 4183–4195 (2020).10. Thiry, D. et al. Genomic relatedness of a canine Lactococcus garvieae to human, animal and environmental isolates. Research in Veterinary Science 137 , 170–173 (2021).11. Evans, J. J., Pasnik, D. J., Klesius, P. H. & Al-Ablani, S. First Report of Streptococcus Agalactiae And Lactococcus Garvieae from a Wild Bottlenose Dolphin (Tursiops Truncatus). Journal of Wildlife Diseases 42 , 561–569 (2006).12. Capobianco, C. M., Bosch, S. N., Stacy, N. I. & Wellehan, J. F. X. Lactococcus garvieae –associated septicemia in a central bearded dragon. J VET Diagn Invest 36 , 477–480 (2024).13. Barakat, R. K., Griffiths, M. W. & Harris, L. J. Isolation and characterization of Carnobacterium, Lactococcus, and Enterococcus spp. from cooked, modified atmosphere packaged, refrigerated, poultry meat. International Journal of Food Microbiology 62 , 83–94 (2000).14. Colagrossi, L. et al. Evidence of pediatric sepsis caused by a drug resistant Lactococcus garvieae contaminated platelet concentrate. Emerging Microbes & Infections 11 , 1325–1334 (2022).15. Savvidis, G., Anatoliotis, C., Kanaki, Z. & Vafeas, G. Epizootic outbreaks of Lactococcosis disease in rainbow trout, Oncorhynchus mykiss (Walbaum), culture in Greece. Bulletin-European Association of Fish Pathologists 27 , 223 (2007).16. Meyburgh, C., Bragg, R. & Boucher, C. Lactococcus garvieae: an emerging bacterial pathogen of fish. Dis. Aquat. Org. 123 , 67–79 (2017).17. Vendrell, D. et al. Lactococcus garvieae in fish: A review. Comparative Immunology, Microbiology and Infectious Diseases 29 , 177–198 (2006).18. Egger, R. C. et al. Emerging fish pathogens Lactococcus petauri and L. garvieae in Nile tilapia (Oreochromis niloticus) farmed in Brazil. Aquaculture 565 , 739093 (2023).19. Stoppani, N. et al. 16S-23S rRNA Internal Transcribed Spacer Region (ITS) Sequencing: A Potential Molecular Diagnostic Tool for Differentiating Lactococcus garvieae and Lactococcus petauri. Microorganisms 11 , 1320 (2023).20. Gibello, A. et al. The zoonotic potential of Lactococcus garvieae: An overview on microbiology, epidemiology, virulence factors and relationship with its presence in foods. Research in Veterinary Science 109 , 59–70 (2016).21. Lin, Y. et al. Comparative Genomic Analyses of Lactococcus garvieae Isolated from Bovine Mastitis in China. Microbiol Spectr 11 , e02995-22 (2023).22. Morita, H. et al. Complete Genome Sequence and Comparative Analysis of the Fish Pathogen Lactococcus garvieae. PLoS ONE 6 , e23184 (2011).23. Miyauchi, E., Toh, H., Nakano, A., Tanabe, S. & Morita, H. Comparative Genomic Analysis of Lactococcus garvieae Strains Isolated from Different Sources Reveals Candidate Virulence Genes. International Journal of Microbiology 2012 , 1–7 (2012).24. Schmidtke, L. M. & Carson, J. Antigen recognition by rainbow trout (Oncorhynchus mykiss) of whole cell proteins expressed by Lactococcus garvieae when obtained directly from fish and under iron limited culture conditions. Veterinary Microbiology 93 , 63–71 (2003).25. Barnes, A. C. et al. Resistance to serum killing may contribute to differences in the abilities of capsulate and non-capsulated isolates ofLactococcus garvieae to cause disease in rainbow trout (Oncorhynchus mykiss L.). Fish & Shellfish Immunology 12 , 155–168 (2002).26. Gendrin, C. et al. Diminished Capsule Exacerbates Virulence, Blood–Brain Barrier Penetration, Intracellular Persistence, and Antibiotic Evasion of Hyperhemolytic Group B Streptococci. The Journal of Infectious Diseases 217 , 1128–1138 (2018).27. Shi, Y.-A., Lu, S.-L., Noda, T., Chiu, C.-H. & Chiang-Ni, C. Capsule-deficient group A Streptococcus evades autophagy-mediated killing in macrophages. mBio 15 , e00771-24 (2024).28. Kietzman, C. C., Gao, G., Mann, B., Myers, L. & Tuomanen, E. I. Dynamic capsule restructuring by the main pneumococcal autolysin LytA in response to the epithelium. Nat Commun 7 , 10859 (2016).29. Ture, M. & Altinok, I. Detection of putative virulence genes of Lactococcus garvieae. Dis. Aquat. Org. 119 , 59–66 (2016).30. Egan, S. A., Kurian, D., Ward, P. N., Hunt, L. & Leigh, J. A. Identification of Sortase A (SrtA) Substrates in Streptococcus uberis : Evidence for an Additional Hexapeptide (LPXXXD) Sorting Motif. J. Proteome Res. 9 , 1088–1095 (2010).31. Diler, O., Altun, S., Adiloglu, A., Kubilay, A. & Istklt, B. First occurrence of Streptococcosis affecting farmed rainbow trout (Oncorhynchus mykiss) in Turkey. Bulletin-European Association of Fish Pathologists 22 , 21–26 (2002).32. Kawanishi, M. et al. Drug resistance and pulsed-field gel electrophoresis patterns of Lactococcus garvieae isolates from cultured Seriola (yellowtail, amberjack and kingfish) in Japan. Lett Appl Microbiol 40 , 322–328 (2005).33. Salvà Serra, F. et al. A protocol for extraction and purification of high-quality and quantity bacterial DNA applicable for genome sequencing: a modified version of the Marmur procedure. Protocol Exchange (2018) doi:10.1038/protex.2018.084.34. Wongkaew, J. et al. Identification and Pathogenicity of Lactococcus Species in Nile Tilapia () and Asian Sea Bass (). Journal of Fish Diseases n/a , e14113.35. Alexander Peltzer et al. nf-core/bacass: v2.1.0 nf-core/bacass: ‘Navy Steel Swordfish’. Zenodo https://doi.org/10.5281/ZENODO.2669428 (2023).36. Chen, S., Zhou, Y., Chen, Y. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34 , i884–i890 (2018).37. Wick, R. R., Judd, L. M., Gorrie, C. L. & Holt, K. E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol 13 , e1005595 (2017).38. Wood, D. E., Lu, J. & Langmead, B. Improved metagenomic analysis with Kraken 2. Genome Biol 20 , 257 (2019).39. Meier-Kolthoff, J. P., Carbasse, J. S., Peinado-Olarte, R. L. & Göker, M. TYGS and LPSN: a database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic Acids Research 50 , D801–D807 (2022).40. Titus Brown, C. & Irber, L. sourmash: a library for MinHash sketching of DNA. JOSS 1 , 27 (2016).41. R Core Team. R: A Language and Environment for Statistical Computing . (R Foundation for Statistical Computing, Vienna, Austria, 2023).42. Bonin, N. et al. MEGARes and AMR++, v3.0: an updated comprehensive database of antimicrobial resistance determinants and an improved software pipeline for classification using high-throughput sequencing. Nucleic Acids Research 51 , D744–D752 (2023).43. Camargo, A. P. et al. Identification of mobile genetic elements with geNomad. Nat Biotechnol (2023) doi:10.1038/s41587-023-01953-y.44. Sirén, K. et al. Rapid discovery of novel prophages using biological feature engineering and machine learning. NAR Genomics and Bioinformatics 3 , lqaa109 (2021).45. Kieft, K., Zhou, Z. & Anantharaman, K. VIBRANT: automated recovery, annotation and curation of microbial viruses, and evaluation of viral community function from genomic sequences. Microbiome 8 , 90 (2020).46. Cock, P. J. A. et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics 25 , 1422–1423 (2009).47. Nayfach, S. et al. CheckV assesses the quality and completeness of metagenome-assembled viral genomes. Nat Biotechnol 39 , 578–585 (2021).48. Bouras, G. et al. Pharokka: a fast scalable bacteriophage annotation tool. Bioinformatics 39 , btac776 (2023).49. Millard, A. et al. taxmyPHAGE: Automated taxonomy of dsDNA phage genomes at the genus and species level. Preprint at https://doi.org/10.1101/2024.08.09.606593 (2024).50. Roux, S. et al. Minimum Information about an Uncultivated Virus Genome (MIUViG). Nat Biotechnol 37 , 29–37 (2019).51. Egorov, A. A. & Atkinson, G. C. LoVis4u: a locus visualization tool for comparative genomics and coverage profiles. NAR Genomics and Bioinformatics 7 , lqaf009 (2025).52. Ridgway, R., Lu, H., Blower, T. R., Evans, N. J. & Ainsworth, S. Genomic and taxonomic evaluation of 38 Treponema prophage sequences. BMC Genomics 25 , 549 (2024).53. Canchaya, C., Proux, C., Fournous, G., Bruttin, A. & Brüssow, H. Prophage Genomics. Microbiol Mol Biol Rev 67 , 238–276 (2003).54. Garneau, J. R. et al. High-throughput identification of viral termini and packaging mechanisms in virome datasets using PhageTermVirome. Sci Rep 11 , 18319 (2021).55. Mills, S., McAuliffe, O. E., Coffey, A., Fitzgerald, G. F. & Ross, R. P. Plasmids of lactococci – genetic accessories or genetic necessities? FEMS Microbiology Reviews 30 , 243–273 (2006).56. Bernheim, A. & Sorek, R. The pan-immune system of bacteria: antiviral defence as a community resource. Nat Rev Microbiol 18 , 113–119 (2020).57. Murphy, J., Mahony, J., Ainsworth, S., Nauta, A. & van Sinderen, D. Bacteriophage Orphan DNA Methyltransferases: Insights from Their Bacterial Origin, Function, and Occurrence. Applied and Environmental Microbiology 79 , 7547–7555 (2013).58. Vela, A. I. et al. The association of Lactococcus petauri with lactococcosis is older than expected. Aquaculture 578 , 740057 (2024).59. Klemm, P. & Schembri, M. A. Bacterial adhesins: function and structure. International Journal of Medical Microbiology 290 , 27–35 (2000).60. Bentley, S. D. et al. Genetic Analysis of the Capsular Biosynthetic Locus from All 90 Pneumococcal Serotypes. PLoS Genet 2 , e31 (2006).61. Soltani, M., Baldisserotto, B., Hosseini Shekarabi, S. P., Shafiei, S. & Bashiri, M. Lactococcosis a Re-Emerging Disease in Aquaculture: Disease Significant and Phytotherapy. Veterinary Sciences 8 , 181 (2021).62. Torres-Corral, Y. & Santos, Y. Predicting antimicrobial resistance of Lactococcus garvieae: PCR detection of resistance genes versus MALDI-TOF protein profiling. Aquaculture 553 , 738098 (2022).63. Nakai, T. et al. Protective effects of bacteriophage on experimental Lactococcus garvieae infection in yellowtail. Dis. Aquat. Org. 37 , 33–41 (1999).64. Cieślik, M. et al. Temperate Bacteriophages-The Powerful Indirect Modulators of Eukaryotic Cells and Immune Functions. Viruses 13 , 1013 (2021).65. Bensing, B. A., Siboo, I. R. & Sullam, P. M. Proteins PblA and PblB of Streptococcus mitis, which promote binding to human platelets, are encoded within a lysogenic bacteriophage. Infect Immun 69 , 6186–6192 (2001).66. Matos, R. C. et al. Enterococcus faecalis prophage dynamics and contributions to pathogenic traits. PLoS Genet 9 , e1003539 (2013).67. Chopin, M.-C., Chopin, A. & Bidnenko, E. Phage abortive infection in lactococci: variations on a theme. Current Opinion in Microbiology 8 , 473–479 (2005).68. Mondal, S. I., Draper, L. A., Ross, R. P. & Hill, C. Bacteriophage endolysins as a potential weapon to combat Clostridioides difficile infection. Gut Microbes 12 , 1813533 (2020).69. Fischetti, V. A. Bacteriophage endolysins: a novel anti-infective to control Gram-positive pathogens. Int J Med Microbiol 300 , 357–362 (2010).70. Eraclio, G. et al. A Virulent Phage Infecting Lactococcus garvieae, with Homology to Lactococcus lactis Phages. Appl Environ Microbiol 81 , 8358–8365 (2015).71. Ghasemi, S. M., Bouzari, M., Yoon, B. H. & Chang, H.-I. Comparative genomic analysis of Lactococcus garvieae phage WP-2, a new member of Picovirinae subfamily of Podoviridae. Gene 551 , 222–229 (2014).72. Hoai, T. D., Nishiki, I. & Yoshida, T. Properties and genomic analysis of Lactococcus garvieae lysogenic bacteriophage PLgT-1, a new member of Siphoviridae, with homology to Lactococcus lactis phages. Virus Research 222 , 13–23 (2016).73. Ramírez-Paredes, J. G. et al. First detection of infectious spleen and kidney necrosis virus (ISKNV) associated with massive mortalities in farmed tilapia in Africa. Transboundary and Emerging Diseases 68 , 1550–1563 (2021).74. Dong, H. T., Ataguba, G. A., Khunrae, P., Rattanarojpong, T. & Senapin, S. Evidence of TiLV infection in tilapia hatcheries from 2012 to 2017 reveals probable global spread of the disease. Aquaculture 479 , 579–583 (2017).75. Thawornwattana, Y. et al. Tilapia lake virus (TiLV): Genomic epidemiology and its early origin. Transbound Emerg Dis 68 , 435–444 (2021). Supplementary Figures Supplementary Figure 1 Proteomic phylogenetic tree generated using VipTree. VipTree was used to generate a proteomic tree using 22 prophage regions identified as part of the present study in addition to 3399 prokaryotic dsDNA viruses. Supplementary Figure 2 Genomic level comparisons of the prophage using average nucleotide identity (ANI). Clustering shown on the left was performed using cutoffs of 95% ANI over 85% alignment fraction. ANI and clustering was performed using CheckV supporting code. Supplementary Figure 3 Genomic organisation of the 22 prophage regions identified as part of the present study. The genomes are reorientated and homology between coding sequences is shown. Coding sequences are coloured based on presumed function. Supplementary Tables Supplementary Table 1 Antimicrobial resistance genes Supplementary Table 2 Virulence Associated Genes Information & Authors Information Version history V1 Version 1 03 June 2025 Copyright This work is licensed under a Non Exclusive No Reuse License. Keywords antibiotic resistance bacteriophages food-borne pathogens genomics Authors Affiliations Adam M. Blanchard 0000-0001-6991-7210 [email protected] University of Nottingham School of Veterinary Medicine and Science View all articles by this author Bailey Secker University of Nottingham School of Veterinary Medicine and Science View all articles by this author Samantha J. Windle University of Nottingham School of Veterinary Medicine and Science View all articles by this author Robert Atterbury University of Nottingham School of Veterinary Medicine and Science View all articles by this author Ha Thanh Dong Asian Institute of Technology School of Environment Resources and Development View all articles by this author Le Thanh Dien Van Lang University View all articles by this author David Huchzermeyer North-West University View all articles by this author Bernard Hang'ombe University of Zambia School of Veterinary Medicine View all articles by this author Saengchan Senapin National Center for Genetic Engineering and Biotechnology View all articles by this author Metrics & Citations Metrics Article Usage 281 views 155 downloads .FvxKWukQNSOunydq8rnd { width: 100px; } Citations Download citation Adam M. Blanchard, Bailey Secker, Samantha J. Windle, et al. The Emergence of Antimicrobial and Phage Resistance in Piscine Lactococcus Isolates and Implications for Future Disease Management in Aquaculture.. Authorea . 03 June 2025. DOI: https://doi.org/10.22541/au.174894333.36446684/v1 If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download. For more information or tips please see 'Downloading to a citation manager' in the Help menu . Format Please select one from the list RIS (ProCite, Reference Manager) EndNote BibTex Medlars RefWorks Direct import Tips for downloading citations document.getElementById('citMgrHelpLink').addEventListener('click', function() { popupHelp(this.href); return false; }); $(".js__slcInclude").on("change", function(e){ if ($(this).val() == 'refworks') $('#direct').prop("checked", false); $('#direct').prop("disabled", ($(this).val() == 'refworks')); }); View Options View options PDF View PDF Figures Tables Media Share Share Share article link Copy Link Copied! Copying failed. Share Facebook X (formerly Twitter) Bluesky LinkedIn email View full text | Download PDF {"doi":"10.22541/au.174894333.36446684/v1","type":"Article"} Now Reading: Share Figures Tables Close figure viewer Back to article Figure title goes here Change zoom level Go to figure location within the article Download figure Toggle share panel Toggle share panel Share Toggle information panel Toggle information panel Go to previous graphic Go to next graphic Go to previous table Go to next table All figures All tables View all material View all material xrefBack.goTo xrefBack.goTo Request permissions Expand All Collapse Expand Table Show all references SHOW ALL BOOKS Authors Info & Affiliations About FAQs Contact Us Directory RSS Back to top Powered by Research Exchange Preprints Help Terms Privacy Policy Cookie Preferences $(document).ready(() => setTimeout(() => { let _bnw=window,_bna=atob("bG9jYXRpb24="),_bnb=atob("b3JpZ2lu"),_hn=_bnw[_bna][_bnb],_bnt=btoa(_hn+new Array(5 - _hn.length % 4).join(" ")); $.get("/resource/lodash?t="+_bnt); },4000)); (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'a004811a9aea58d3',t:'MTc3OTU0Mzk2MA=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.