Development and application of a whole genome amplicon sequencing method for infectious salmon anemia virus (ISAV).

Infectious salmon anemia (ISA) is an infectious disease primarily affecting farmed Atlantic salmon, Salmo salar, which is caused by the ISA virus (ISAV). ISAV belongs to the Orthomyxoviridae family. The disease is a serious condition resulting in reduced fish welfare and high mortality. In this study, we designed an amplicon-based sequencing protocol for whole genome sequencing of ISAV. The method consists of 80 ISAV-specific primers that cover 92% of the virus genome and was designed to be used on an Illumina MiSeq platform. The sequencing accuracy was investigated by comparing sequences with previously published Sanger sequences. The sequences obtained were nearly identical to those obtained by Sanger sequencing, thus demonstrating that sequences produced by this amplicon sequencing protocol had an acceptable accuracy. The amplicon-based sequencing method was used to obtain the whole genome sequence of 12 different ISAV isolates from a small local epidemic in the northern part of Norway. Analysis of the whole genome sequences revealed that segment reassortment took place between some of the isolates and could identify which segments that had been reassorted.

Draft genome sequences of 13 Tenacibaculum maritimum isolates from farmed Norwegian Cyclopterus lumpus (lumpfish) and Scophthalmus maximus (turbot).

Thirteen bacterial isolates of Tenacibaculum maritimum were sequenced and assembled. The strains were isolated from four disease outbreaks in farmed marine fish in Norway. Eight isolates were from Cyclopterus lumpus (lumpfish), and five were from Scophthalmus maximus (turbot). Overall, sequence similarity did not correlate with host species or geographic location.

Polyether ionophore resistance in a one health perspective.

Antimicrobial resistance is a major threat to human health and must be approached from a One Health perspective. Use of antimicrobials in animal husbandry can lead to dissemination and persistence of resistance in human pathogens. Polyether ionophores (PIs) have antimicrobial activities and are among the most extensively used feed additives for major production animals. Recent discoveries of genetically encoded PI resistance mechanisms and co-localization of resistance mechanisms against PIs and antimicrobials used in human medicine on transferrable plasmids, have raised concerns that use of PIs as feed additives bear potential risks for human health. This review summarizes the current knowledge on PI resistance and discusses the potential consequences of PI-usage as feed additives in a One Health perspective.

Phylogeography and host specificity of Pasteurellaceae pathogenic to sea-farmed fish in the north-east Atlantic.

The present study was undertaken to address the recent spate of pasteurellosis outbreaks among sea-farmed Atlantic salmon (Salmo salar) in Norway and Scotland, coinciding with sporadic disease episodes in lumpfish (Cyclopterus lumpus) used for delousing purposes in salmon farms. Genome assemblies from 86 bacterial isolates cultured from diseased salmon or lumpfish confirmed them all as bona fide members of the Pasteurellaceae family, with phylogenetic reconstruction dividing them into two distinct branches sharing <88% average nucleotide identity. These branches therefore constitute two separate species, namely Pasteurella skyensis and the as-yet invalidly named "Pasteurella atlantica". Both species further stratify into multiple discrete genomovars (gv.) and/or lineages, each being nearly or fully exclusive to a particular host, geographic region, and/or time period. Pasteurellosis in lumpfish is, irrespective of spatiotemporal origin, linked almost exclusively to the highly conserved "P. atlantica gv. cyclopteri" (Pac). In contrast, pasteurellosis in Norwegian sea-farmed salmon, dominated since the late-1980s by "P. atlantica gv. salmonicida" (Pas), first saw three specific lineages (Pas-1, -2, and -3) causing separate, geographically restricted, and short-lived outbreaks, before a fourth (Pas-4) emerged recently and became more widely disseminated. A similar situation involving P. skyensis (Ps) has apparently been unfolding in Scottish salmon farming since the mid-1990s, where two historic (Ps-1 and -2) and one contemporary (Ps-3) lineages have been recorded. While the epidemiology underlying all these outbreaks/epizootics remains unclear, repeated detection of 16S rRNA gene amplicons very closely related to P. skyensis and "P. atlantica" from at least five cetacean species worldwide raises the question as to whether marine mammals may play a part, possibly as reservoirs. In fact, the close relationship between the studied isolates and Phocoenobacter uteri associated with harbor porpoise (Phocoena phocoena), and their relatively distant relationship with other members of the genus Pasteurella, suggests that both P. skyensis and "P. atlantica" should be moved to the genus Phocoenobacter.

Escherichia coli multilocus sequence type 38 from humans and broiler production represent distinct monophyletic groups.

Escherichia coli belonging to multilocus sequence type 38 (ST38) is a well-known cause of extra-intestinal infections in humans, and are frequently associated with resistance to extended-spectrum cephalosporins (ESCs). Resistance to carbapenems, mediated by blaOXA-genes has also been reported in this ST. Recently, the European Centre for Disease Prevention and Control (ECDC) released a rapid risk assessment on the increased detection of OXA-244 producing E. coli ST38 in humans, requesting further knowledge to determine the source. ST38 is also one of the most common STs among ESC-resistant E. coli from broiler production. Our aim was to investigate the genetic characteristics and relationship between E. coli ST38 from broiler production and humans, and to investigate if there has been a potential spillover between these sources. A total of 288 E. coli ST38 genomes isolated from humans in Europe (collected 2009-2019) and from Nordic broiler production (collected 2011-2014) were analyzed. The results showed distinct monophyletic clades associated to humans and broiler production. Furthermore, there were differences in the ESC resistance genes present in E. coli ST38 from the two sources. The blaOXA-244 gene was not present in E. coli from broiler production. Our results show that ST38 from humans and broiler production belong to well-separated clades, and suggest that the increased detection of OXA-244-producing E. coli ST38 in humans is not associated with spillover from broiler production.

Draft Genome Sequences of 23 Tenacibaculum Isolates from Farmed Norwegian Lumpfish, Cyclopterus lumpus.

Draft genome sequences of 23 Tenacibaculum sp. strains that were isolated from Cyclopterus lumpus (lumpfish) were investigated to elucidate possible routes of transmission between Salmo salar (Atlantic salmon) and lumpfish.

Closed Genome Sequences of Providencia alcalifaciens Isolates from Dogs.

Eight Providencia alcalifaciens isolates from eight different dogs in Norway with acute hemorrhagic diarrhea were sequenced. Based on Illumina and Oxford Nanopore Technologies sequencing, all of the genomes were complete and closed after hybrid assembly.

Tenacibaculosis in Norwegian Atlantic salmon (Salmo salar) cage-farmed in cold sea water is primarily associated with Tenacibaculum finnmarkense genomovar finnmarkense.

Skin conditions associated with Tenacibaculum spp. constitute a significant threat to the health and welfare of sea-farmed Atlantic salmon (Salmo salar L.) in Norway. Fifteen presumptive tenacibaculosis outbreaks distributed along the Norwegian coast during the late winter and spring of 2018 were investigated. Bacteriological culture confirmed the presence of Tenacibaculum spp. Seventy-six isolates cultured from individual fish were selected and subjected to whole-genome sequencing and MALDI-TOF MS analysis. Average nucleotide identity and MALDI-TOF analyses confirmed the presence of T. finnmarkense and T. dicentrarchi, with further division of T. finnmarkense into genomovars (gv.) finnmarkense and ulcerans. Core genome multilocus sequence typing (cgMLST) and single-nucleotide polymorphism (SNP) analyses identified the presence of a genetically conserved cluster of gv. finnmarkense isolates against a background of relatively genetically diverse gv. finnmarkense and gv. ulcerans isolates in 13 of the 15 studied cases. This clustering strongly suggests a link between T. finnmarkense gv. finnmarkense and development of clinical tenacibaculosis in sea-farmed Norwegian salmon in the late winter and spring. Analysis of 25 Tenacibaculum isolates collected during the spring of 2019 from similar cases identified a similar distribution of genotypes. Low water temperatures were common to all cases, and most incidences involved relatively small fish shortly after sea transfer, suggesting that these fish are particularly predisposed to Tenacibaculum infection.

An Official Outbreak Investigation of Acute Haemorrhagic Diarrhoea in Dogs in Norway Points to Providencia alcalifaciens as a Likely Cause.

An outbreak investigation was initiated in September 2019, following a notification to the Norwegian Food Safety Authority (NFSA) of an unusually high number of dogs with acute haemorrhagic diarrhoea (AHD) in Oslo. Diagnostic testing by reporting veterinarians had not detected a cause. The official investigation sought to identify a possible common cause, the extent of the outbreak and prevent spread. Epidemiological data were collected through a survey to veterinarians and interviews with dog owners. Diagnostic investigations included necropsies and microbiological, parasitological and toxicological analysis of faecal samples and food. In total, 511 dogs with acute haemorrhagic diarrhoea were registered between 1 August and 1 October. Results indicated a common point source for affected dogs, but were inconclusive with regard to common exposures. A notable finding was that 134 of 325 faecal samples (41%) cultured positive for Providencia alcalifaciens. Whole genome sequencing (WGS) of 75 P. alcalifaciens isolates from 73 dogs revealed that strains from 51 dogs belonged to the same WGS clone. Findings point to P. alcalifaciens as implicated in the outbreak, but investigations are needed to reveal the pathogenic potential of P. alcalifaciens in dogs and its epidemiology.

One Health Surveillance Codex: promoting the adoption of One Health solutions within and across European countries.

Cross-sector communication, collaboration and knowledge exchange are still significant challenges for practical adoption of the One Health paradigm. To address these needs the "One Health Surveillance Codex" (OHS Codex) was established to provide a framework for the One Health community to continuously share practical solutions (e.g. tools, technical resources, guidance documents and experiences) applicable for national and international stakeholders from different One Health Surveillance sectors. Currently, the OHS Codex provides a number of resources that support the adoption of the OH paradigm in areas linked to the harmonization and interpretation of surveillance data. The OHS Codex framework comprises four high-level "action" principles, which respectively support collaboration, knowledge exchange, data interoperability, and dissemination. These principles match well with priority areas identified in the "Tripartite Guide to Addressing Zoonotic Diseases in Countries" published by WHO, FAO and OIE. Within each of the four principles, the OHS Codex provides a collection of useful resources as well as pointers to success stories for the application of these resources. As the OHS Codex is designed as an open community framework, it will continuously evolve and adapt to the needs of the OH community in the future.

Tenacibaculum piscium sp. nov., isolated from skin ulcers of sea-farmed fish, and description of Tenacibaculum finnmarkense sp. nov. with subdivision into genomovars finnmarkense and ulcerans.

Results of previous multilocus sequence and whole-genome-based analyses have suggested that a homogeneous group of isolates belonging to the genus Tenacibaculum, represented by strain TNO020T and associated with skin ulcer development in sea-farmed fish, represents an as-yet-undescribed species. Comparative whole-genome analysis performed in the present study clustered five isolates, including TNO020T, in a distinct lineage within the genus Tenacibaculum. Phenotypic differences, high intra-cluster average nucleotide identity (ANI) values and low ANI values with other Tenacibaculum species support the proposal of a novel species, for which we propose the name Tenacibaculum piscium sp. nov. with strain TNO020T (=CCUG 73833T=NCIMB 15240T) as the type strain. Further, large-scale genome analyses confirmed the existence of two different phylogenetic lineages within 'T. finnmarkense', a species effectively but not validly published previously. ANI values just above the species delineation threshold of 95-96 % confirmed that both lineages belong to the same species. This result was also supported by DNA-DNA hybridization values. Phenotypically, the two conspecific lineages are distinguishable by differences in growth temperature range and ability to degrade l-proline. For the group of isolates already commonly known as 'T. finnmarkense', we propose the name Tenacibaculum finnmarkense sp. nov., with strain TNO006T (=CCUG 73831T=NCIMB 15238T) as the type strain. We further propose the subdivision of T. finnmarkense sp. nov. into two genomovars, T. finnmarkense genomovar finnmarkense with strain TNO006T (=CCUG 73831T=NCIMB 15238T) as the type strain and T. finnmarkense genomovar ulcerans with strain TNO010T (=CCUG 73832T=NCIMB 15239T) as the type strain.

Comparative Genome Analyses of Wild Type- and Quinolone Resistant Escherichia coli Indicate Dissemination of QREC in the Norwegian Broiler Breeding Pyramid.

Quinolones are important antimicrobials for both humans and animals, and resistance toward these compounds is a serious threat to public health. In Norway, quinolone resistant E. coli (QREC) have been detected at low levels in a high proportion of broiler flocks, even without the use of quinolones in rearing of broilers. Due to the pyramidal structure of broiler breeding, QREC isolates may be disseminated from grandparent animals down through the pyramid. However, quinolone resistance can also develop in wild type E. coli through specific chromosomal mutations, and by horizontal acquisition of plasmid-mediated quinolone resistance genes. The goal of this study was to determine whether QREC is disseminated through the broiler breeding pyramid or developed locally at some stage in the broiler production chain. For this purpose, we whole genome sequenced wild type- and QREC isolates from broiler and parent flocks that had been isolated in the Norwegian monitoring program for antimicrobial resistance in feed, food and animals (NORM-VET) between 2006 and 2017, from 22 different production sites. The sequencing data was used for typing of the isolates, phylogenetic analysis and identification of relevant resistance mechanisms. Highly similar QREC isolates were identified within major sequence types from multiple production sites, suggesting dissemination of QREC isolates in the broiler production chain. The occurrence of potential resistance development among the WT E. coli was low, indicating that this may be a rare phenomenon in the Norwegian broiler production. The results indicate that the majority of the observed QREC at the bottom of the broiler production pyramid originates from parent or grandparent animals. These results highlight the importance of surveillance at all levels of the broiler production pyramid and of implementation of proper biosecurity measures to control dissemination of QREC.

Dissemination of Quinolone-Resistant Escherichia coli in the Norwegian Broiler and Pig Production Chains and Possible Persistence in the Broiler Production Environment.

In Norway, the use of quinolones in livestock populations is very low, and prophylactic use is prohibited. Despite this, quinolone-resistant Escherichia coli (QREC) isolates are present at low levels in several animal species. The source of these QREC isolates is unknown. The aim of this study was to characterize and compare QREC isolates from different animal species to identify putative factors that may promote the occurrence of QREC. A total of 280 QREC isolates, from broilers, pigs, red foxes, and wild birds, were whole-genome sequenced and analyzed. Well-known chromosomal and plasmid-mediated resistance mechanisms were identified. In addition, mutations in marR, marA, and rpoB causing novel amino acid substitutions in their respective proteins were detected. Phylogenetic analyses were used to determine the relationships between the isolates. Quinolone resistance mechanism patterns appeared to follow sequence type groups. Similar QREC isolates with similar resistance mechanism patterns were detected from the samples, and further phylogenetic analysis indicated close evolutionary relationships between specific isolates from different sources. This suggests the dissemination of highly similar QREC isolates between animal species and also the persistence of QREC strains within the broiler production chain. This highlights the importance of both control measures at the top of the production chain as well as biosecurity measures to avoid the further dissemination and persistence of QREC in these environments.IMPORTANCE Since antimicrobial usage is low in Norwegian animal husbandry, Norway is an ideal country to study antimicrobial resistance in the absence of selective pressure from antimicrobial usage. In particular, the usage of quinolones is very low, which makes it possible to investigate the spread and development of quinolone resistance in natural environments. Comparison of quinolone-resistant E. coli (QREC) isolates from livestock and wild animals in light of this low quinolone usage provides new insights into the development and dissemination of QREC in both natural and production environments. With this information, preventive measures may be taken to prevent further dissemination within Norwegian livestock and between other animals, thus maintaining the favorable situation in Norway.

Multilocus Variable-Number Tandem-Repeat Analysis of Yersinia ruckeri Confirms the Existence of Host Specificity, Geographic Endemism, and Anthropogenic Dissemination of Virulent Clones.

A multilocus variable-number tandem-repeat analysis (MLVA) assay was developed for epizootiological study of the internationally significant fish pathogen Yersinia ruckeri, which causes yersiniosis in salmonids. The assay involves amplification of 10 variable-number tandem-repeat (VNTR) loci in two five-plex PCRs, followed by capillary electrophoresis. A collection of 484 Y. ruckeri isolates, originating from various biological sources and collected from four continents over 7 decades, was analyzed. Minimum-spanning-tree cluster analysis of MLVA profiles separated the studied population into nine major clonal complexes and a number of minor clusters and singletons. The major clonal complexes could be associated with host species, geographic origin, and serotype. A single large clonal complex of serotype O1 isolates dominating the yersiniosis situation in international rainbow trout farming suggests anthropogenic spread of this clone, possibly related to transport of fish. Moreover, subclustering within this clonal complex indicates putative transmission routes and multiple biotype shift events. In contrast to the situation in rainbow trout, Y. ruckeri strains associated with disease in Atlantic salmon appear as more or less geographically isolated clonal complexes. A single complex of serotype O1 exclusive to Norway was found to be responsible for almost all major yersiniosis outbreaks in modern Norwegian salmon farming, and site-specific subclustering further indicates persistent colonization of freshwater farms in Norway. Identification of genetically diverse Y. ruckeri isolates from clinically healthy fish and environmental sources also suggests the widespread existence of less-virulent or avirulent strains.IMPORTANCE This comprehensive population study substantially improves our understanding of the epizootiological history and nature of an internationally important fish-pathogenic bacterium. The MLVA assay developed and presented represents a high-resolution typing tool particularly well suited for Yersinia ruckeri infection tracing, selection of strains for vaccine inclusion, and risk assessment. The ability of the assay to separate isolates into geographically linked and/or possibly host-specific clusters reflects its potential utility for maintenance of national biosecurity. The MLVA is internationally applicable and robust, and it provides clear, unambiguous, and easily interpreted results. Typing is reasonably inexpensive, with a moderate technological requirement, and may be completed from a harvested colony within a single working day. As the resulting MLVA profiles are readily portable, any Y. ruckeri strain may rapidly be placed in a global epizootiological context.

Occurrence of quinolone resistant E. coli originating from different animal species in Norway.

The aim of this study was to describe and compare the occurrence of quinolone resistant Escherichia coli (QREC) in various animal species in relation to human population density. Data from the Norwegian monitoring programme for antimicrobial resistance in feed, food and animals from 2006 to 2016 was compiled and analysed. In total, 4568 E. coli isolates were included in this study. The isolates originated from broilers, layers, cattle, turkeys, dogs, wild birds, red foxes, reindeer, sheep, horses and pigs. Data regarding the geographical location of sampling was obtained for 4050 of these isolates and used to categorize the isolates depending on the human population density of the area. In total, 1.4% of the isolates were categorized as quinolone resistant. Compared to most European countries, there was an overall low occurrence of QREC in various animal species in Norway, though with an interspecies variation with the highest occurrence in broilers and wild birds (p < 0.05). Human population density was not associated with the occurrence of QREC. Since fluoroquinolones are not used prophylactically and in almost negligent amounts in various species in Norway, the interspecies variation in the occurrence of QREC suggests that other factors than fluoroquinolone use may be important in the development of QREC.

NucDiff: in-depth characterization and annotation of differences between two sets of DNA sequences.

BACKGROUND: Comparing sets of sequences is a situation frequently encountered in bioinformatics, examples being comparing an assembly to a reference genome, or two genomes to each other. The purpose of the comparison is usually to find where the two sets differ, e.g. to find where a subsequence is repeated or deleted, or where insertions have been introduced. Such comparisons can be done using whole-genome alignments. Several tools for making such alignments exist, but none of them 1) provides detailed information about the types and locations of all differences between the two sets of sequences, 2) enables visualisation of alignment results at different levels of detail, and 3) carefully takes genomic repeats into consideration. RESULTS: We here present NucDiff, a tool aimed at locating and categorizing differences between two sets of closely related DNA sequences. NucDiff is able to deal with very fragmented genomes, repeated sequences, and various local differences and structural rearrangements. NucDiff determines differences by a rigorous analysis of alignment results obtained by the NUCmer, delta-filter and show-snps programs in the MUMmer sequence alignment package. All differences found are categorized according to a carefully defined classification scheme covering all possible differences between two sequences. Information about the differences is made available as GFF3 files, thus enabling visualisation using genome browsers as well as usage of the results as a component in an analysis pipeline. NucDiff was tested with varying parameters for the alignment step and compared with existing alternatives, called QUAST and dnadiff. CONCLUSIONS: We have developed a whole genome alignment difference classification scheme together with the program NucDiff for finding such differences. The proposed classification scheme is comprehensive and can be used by other tools. NucDiff performs comparably to QUAST and dnadiff but gives much more detailed results that can easily be visualized. NucDiff is freely available on https://github.com/uio-cels/NucDiff under the MPL license.

Use of high throughput sequencing and light microscopy show contrasting results in a study of phytoplankton occurrence in a freshwater environment.

Assessing phytoplankton diversity is of primary importance for both basic and applied ecological studies. Following the advances in molecular methods, phytoplankton studies are switching from using classical microscopy to high throughput sequencing approaches. However, methodological comparisons of these approaches have rarely been reported. In this study, we compared the two methods, using a unique dataset of multiple water samples taken from a natural freshwater environment. Environmental DNA was extracted from 300 water samples collected weekly during 20 years, followed by high throughput sequencing of amplicons from the 16S and 18S rRNA hypervariable regions. For each water sample, phytoplankton diversity was also estimated using light microscopy. Our study indicates that species compositions detected by light microscopy and 454 high throughput sequencing do not always match. High throughput sequencing detected more rare species and picoplankton than light microscopy, and thus gave a better assessment of phytoplankton diversity. However, when compared to light microscopy, high throughput sequencing of 16S and 18S rRNA amplicons did not adequately identify phytoplankton at the species level. In summary, our study recommends a combined strategy using both morphological and molecular techniques.

Pathogens in Urine from a Female Patient with Overactive Bladder Syndrome Detected by Culture-independent High Throughput Sequencing: A Case Report.

INTRODUCTION: Overactive bladder syndrome (OAB) is described as urgency, with or without urgency incontinence. A range of medical conditions shares the symptoms of OAB, however the diagnosis is contingent on the exclusion of urinary tract infection (UTI). Knowing that urine dipstick and routine culture of bacteria can miss UTI diagnosis caused by low-count bacteriuria or "difficult-to-culture" pathogens, we examined a case of OAB with a culture-independent approach. CASE PRESENTATION: A 61-year-old Norwegian female with a long history of urinary symptoms and a diagnosis of OAB was selected as a suitable subject for a culture-independent 16S rDNA analysis on the patient´s urine. The patient's medical records showed no history of recurrent UTI, however, when the urine specimen was sent to routine culture at the time of study it showed a significant bacteriuria caused by a single bacterium, and the patient was prescribed antibiotics. The 16S rDNA analysis revealed not one, but many different bacteria, including a considerable amount of fastidious bacteria, indicating a polymicrobial state. One year later, the subject was still experiencing severe symptoms, and a follow-up analysis was performed. This time the urine-culture was negative, however, the 16S rDNA profile was quite similar to that of the first sample, again displaying a complex bacterial profile. CONCLUSION: The use of 16S rDNA pyrosequencing and sequence analysis to uncover "difficult-to-culture" bacteria should be considered when examining patients with chronic urinary symptoms. These methods may contribute to further elucidation of the etiology of overactive bladder syndrome and other urinary syndromes.

Tiling array study of MNNG treated Escherichia coli reveals a widespread transcriptional response.

The alkylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG) is known to trigger the adaptive response by inducing the ada-regulon - consisting of three DNA repair enzymes Ada, AlkB, AlkA and the enigmatic AidB. We have applied custom designed tiling arrays to study transcriptional changes in Escherichia coli following a MNNG challenge. Along with the expected upregulation of the adaptive response genes (ada, alkA and alkB), we identified a number of differentially expressed transcripts, both novel and annotated. This indicates a wider regulatory response than previously documented. There were 250 differentially-expressed and 2275 similarly-expressed unannotated transcripts. We found novel upregulation of several stress-induced transcripts, including the SOS inducible genes recN and tisAB, indicating a novel role for these genes in alkylation repair. Furthermore, the ada-regulon A and B boxes were found to be insufficient to explain the regulation of the adaptive response genes after MNNG exposure, suggesting that additional regulatory elements must be involved.

From green to red: horizontal gene transfer of the phycoerythrin gene cluster between Planktothrix strains.

Horizontal gene transfer is common in cyanobacteria, and transfer of large gene clusters may lead to acquisition of new functions and conceivably niche adaption. In the present study, we demonstrate that horizontal gene transfer between closely related Planktothrix strains can explain the production of the same oligopeptide isoforms by strains of different colors. Comparison of the genomes of eight Planktothrix strains revealed that strains producing the same oligopeptide isoforms are closely related, regardless of color. We have investigated genes involved in the synthesis of the photosynthetic pigments phycocyanin and phycoerythrin, which are responsible for green and red appearance, respectively. Sequence comparisons suggest the transfer of a functional phycoerythrin gene cluster generating a red phenotype in a strain that is otherwise more closely related to green strains. Our data show that the insertion of a DNA fragment containing the 19.7-kb phycoerythrin gene cluster has been facilitated by homologous recombination, also replacing a region of the phycocyanin operon. These findings demonstrate that large DNA fragments spanning entire functional gene clusters can be effectively transferred between closely related cyanobacterial strains and result in a changed phenotype. Further, the results shed new light on the discussion of the role of horizontal gene transfer in the sporadic distribution of large gene clusters in cyanobacteria, as well as the appearance of red and green strains.