Genome Phylogeny Research Articles

Closed genomes uncover a saltwater species of Candidatus Electronema and shed new light on the boundary between marine and freshwater cable bacteria

Cable bacteria of the Desulfobulbaceae family are centimeter-long filamentous bacteria, which are capable of conducting long-distance electron transfer. Currently, all cable bacteria are classified into two candidate genera: Candidatus Electronema, typically found in freshwater environments, and Candidatus Electrothrix, typically found in saltwater environments. This taxonomic framework is based on both 16S rRNA gene sequences and metagenome-assembled genome (MAG) phylogenies. However, most of the currently available MAGs are highly fragmented, incomplete, and thus likely miss key genes essential for deciphering the physiology of cable bacteria. Also, a closed, circular genome of cable bacteria has not been published yet. To address this, we performed Nanopore long-read and Illumina short-read shotgun sequencing of selected environmental samples and a single-strain enrichment of Ca. Electronema aureum. We recovered multiple cable bacteria MAGs, including two circular and one single-contig. Phylogenomic analysis, also confirmed by 16S rRNA gene-based phylogeny, classified one circular MAG and the single-contig MAG as novel species of cable bacteria, which we propose to name Ca. Electronema halotolerans and Ca. Electrothrix laxa, respectively. The Ca. Electronema halotolerans, despite belonging to the previously recognized freshwater genus of cable bacteria, was retrieved from brackish-water sediment. Metabolic predictions showed several adaptations to a high salinity environment, similar to the “saltwater” Ca. Electrothrix species, indicating how Ca. Electronema halotolerans may be the evolutionary link between marine and freshwater cable bacteria lineages.

Salmonella Heidelberg and Salmonella Minnesota in Brazilian broilers: Genomic characterization of third-generation cephalosporin and fluoroquinolone-resistant strains.

Salmonella serovars Heidelberg and Minnesota encoding antimicrobial resistance to third-generation cephalosporins and fluoroquinolones are often detected in poultry/poultry meat. We analysed the genomes of 10 Salmonella Heidelberg (SH) and 4 Salmonella Minnesota (SM) from faecal isolates of Brazilian poultry. These featured virulent and multidrug-resistant characteristics, with AmpC beta-lactamase (blaCMY-2 ) predominance (9/14), for all SM (4/4) and some SH (3/10) located on IncC plasmid replicons. IncC carrying blaCTX-M-2 was only detected among SH (3/10). Mutation in the gyrA/parC genes was present in all SH, whereas SM harboured parC mutation plus qnrB19 on ColRNAI plasmids (3/4). In silico resistance overall corroborated with phenotypic results. Core genome phylogenies showed close clustering and high similarities between the Brazilian and poultry meat/food isolates from Europe, and to human isolates from European countries with documented import of Brazilian poultry meat. Conjugation assays with SM successfully transferred blaCMY-2 , and qnrB19 to an Escherichia coli recipient. The findings reinforce the ongoing antimicrobial resistance acquisition of SH and Minnesota and the risks for disseminating resistant strains and/or mobile elements which may increasingly affect importing countries and the need for controlling AMR in major poultry-exporting countries like Brazil.

Genomic and Metabolite Profiling Reveal a Novel Streptomyces Strain, QHH-9511, from the Qinghai-Tibet Plateau.

The prevalence of superbugs, represented by methicillin-resistant Staphylococcus aureus (MRSA), has become a serious clinical and public safety concern with rising incidence in hospitals. Polyketides with diverse chemical structures harbor many antimicrobial activities, including those of rifampin and rapamycin against MRSA. Streptomyces sp. QHH-9511 was isolated from a niche habitat in the Qinghai-Tibet Plateau and used to produce antibacterial metabolites. Herein, an integrated approach combining genome mining and metabolic analysis were employed to decipher the chemical origin of the antibacterial components with pigmented properties in strain QHH-9511, a novel Streptomyces species from a lichen symbiont on the Qinghai-Tibet Plateau. Genomic phylogeny assembled at the chromosome level revealed its unique evolutionary state. Further genome mining uncovered 36 candidate gene clusters, most of which were uncharacterized. Meanwhile, based on liquid chromatography coupled to diode array detection mass spectrometry, a series of granaticins, BSMs, chromones, phaeochromycins, and related molecules were discovered by using the Global Natural Product Social molecular networking platform. Subsequently, several pigment compounds were isolated and identified by high-resolution mass spectrometry and/or nuclear magnetic resonance, among which the structure-activity relationships of seven aromatic polyketides showed that the fused lactone ring of the C-2 carboxyl group could increase antibacterial activity. Genetic experiments indicated that all seven aromatic polyketides are a series of metabolic shunts produced by a single type II polyketide synthase (PKS) cluster. Comparative genomic analysis of granaticin producers showed that the granaticin gene cluster is widely distributed. This study provides an efficient method to combine genome mining and metabolic profiling techniques to uncover bioactive metabolites derived from specific habitats, while deepening our understanding of aromatic polyketide biosynthesis. IMPORTANCE Undescribed microorganisms from special habitats are being screened for anti-superbug drug molecules. In a project to screen actinomycetes for anti-MRSA activity, we isolated a Streptomyces strain from Qinghai Lake lichens. The phylogeny based on the genome assembled at the chromosome level revealed this strain's unique evolutionary state. The chemical origins of the antibacterial components with pigment properties in strain QHH-9511 were determined using an integrated approach combining genome mining and metabolic analysis. Further genome mining uncovered 36 secondary metabolite gene clusters, the majority of which were previously unknown. A series of aromatic compounds were discovered using molecular network analysis, separation, and extraction. Genetic experiments revealed that all seven aromatic polyketides are a series of metabolic shunts produced by a single cluster of type II PKSs. This study describes a method for identifying novel Streptomyces from specific habitats by combining genome mining with metabolic profiling techniques.

Comparative phylogenomics of Trueperella pyogenes reveals host-based distinction of strains.

Trueperella pyogenes, an opportunistic pathogen causes various ailments in different animals. Different strains from different animals have distinct characters phenotypically and genotypically. Hence understanding the strains in a particular geographical location helps in framing the preventive measures. Comparative genomics of all the available T. pyogenes genome in the NCBI was conducted to understand the relatedness among strains. Whole genome phylogeny showed host associated clustering of strains recovered from swine lungs. Core genome phylogeny also showed host associated clustering mimicking whole genome phylogeny results. MLST analysis showed that there was higher diversity among cattle strains. Multidimensional scaling revealed five swine clusters, two cattle and buffalo clusters. Pangenome analysis also showed that T. pyogenes had an open genome with 57.09% accessory genome. Host specific genes were identified by pangenome analysis, and (R)-citramalate synthase was specific for swine strains of Asian origin. Host specifc genes identified by pangenome analysis can be exploited for developing a molecular assay to specifically identify the strains. The study shows that MLST having higher discriminatory power can be used as an epidemiological tool for strain discrimination of T. pyogenes.

Molecular characterization of a novel victorivirus isolated from the phytopathogenic fungus Phaeobotryon rhois.

Phaeobotryon rhois is an important pathogenic fungus that causes dieback and canker disease of woody hosts. A novel mycovirus, tentatively named "Phaeobotryon rhois victorivirus 1" (PrVV1), was identified in P. rhois strain SX8-4. The PrVV1 has a double-stranded RNA (dsRNA) genome that is 5,224 base pairs long and contains two open reading frames (ORF1 and ORF2), which overlap at a AUGA sequence. ORF1 encodes a polypeptide of 786 amino acids (aa) that contains the conserved coat protein (CP) domain of victoriviruses, while ORF2, encodes a large polypeptide of 826 aa that contains the conserved RNA-dependent RNA polymerase (RdRp) domain of victoriviruses. Our analysis of genomic structure, homology, and phylogeny indicated that PrVV1 is a novel member of the genus Victorivirus in the family Totiviridae. This is the first report of the complete genome sequence of a victorivirus that infects P. rhois.

Complete chloroplast genome sequencing and phylogeny of wild Atractylodes lancea from Yuexi, Anhui province

This study investigated the choroplast genome sequence of wild Atractylodes lancea from Yuexi in Anhui province by high-throughput sequencing, followed by characterization of the genome structure, which laid a foundation for the species identification, analysis of genetic diversity, and resource conservation of A. lancea. To be specific, the total genomic DNA was extracted from the leaves of A. lancea with the improved CTAB method. The chloroplast genome of A. lancea was sequenced by the high-throughput sequencing technology, followed by assembling by metaSPAdes and annotation by CPGAVAS2. Bioiformatics methods were employed for the analysis of simple sequence repeats(SSRs), inverted repeat(IR) border, codon bias, and phylogeny. The results showed that the whole chloroplast genome of A. lancea was 153 178 bp, with an 84 226 bp large single copy(LSC) and a 18 658 bp small single copy(SSC) separated by a pair of IRs(25 147 bp). The genome had the GC content of 37.7% and 124 genes: 87 protein-coding genes, 8 rRNA genes, and 29 tRNA genes. It had 26 287 codons and encoded 20 amino acids. Phylogenetic analysis showed that Atractylodes species clustered into one clade and that A. lancea had close genetic relationship with A. koreana. This study established a method for sequencing the chloroplast genome of A. lancea and enriched the genetic resources of Compositae. The findings are expected to lay a foundation for species identification, analysis of genetic diversity, and resource conservation of A. lancea.

Evolutionary relationships of mitogenomes in a recently radiated Old World avian family

Environmentally heterogeneous mountains provide opportunities for rapid diversification and speciation. The family Prunellidae (accentors) is a group of birds comprising primarily mountain specialists that have recently radiated across the Palearctic region. This rapid diversification poses challenges to resolving their phylogeny. Herein we sequenced the complete mitogenomes and estimated the phylogeny using all 12 (including 28 individuals) currently recognized species of Prunellidae. We reconstructed the mitochondrial genome phylogeny using 13 protein-coding genes of 12 species and 2 Eurasian Tree Sparrows (Passer montanus). Phylogenetic relationships were estimated using a suite of analyses: maximum likelihood, maximum parsimony and the coalescent-based SVDquartets. Divergence times were estimated by implementing a Bayesian relaxed clock model in BEAST2. Based on the BEAST time-calibrated tree, we implemented an ancestral area reconstruction using RASP v.4.3. Our phylogenies based on the maximum likelihood, maximum parsimony and SVDquartets approaches support a clade of large-sized accentors (subgenus Laiscopus) to be sister to all other accentors with small size (subgenus Prunella). In addition, the trees also support the sister relationship of P. immaculata and P. rubeculoides ​+ ​P.atrogularis with 100% bootstrap support, but the relationships among the remaining eight species in the Prunella clade are poorly resolved. These species cluster in different positions in the three phylogenetic trees and the nodes are often poorly supported. The five nodes separating the seven species diverged simultaneously within less than half million years (i.e., between 2.71 and 3.15 million years ago), suggesting that the recent radiation is likely responsible for rampant incomplete lineage sorting and gene tree conflicts. Ancestral area reconstruction indicates a central Palearctic region origin for Prunellidae. Our study highlights that whole mitochondrial genome phylogeny can resolve major lineages within Prunellidae but is not sufficient to fully resolve the relationship among the species in the Prunella clade that almost simultaneously diversify during a short time period. Our results emphasize the challenge to reconstruct reliable phylogenetic relationship in a group of recently radiated species.

Mitogenomic Codon Usage Patterns of Superfamily Certhioidea (Aves, Passeriformes): Insights into Asymmetrical Bias and Phylogenetic Implications

The superfamily Certhioidea currently comprises five families. Due to the rapid diversification, the phylogeny of Certhioidea is still controversial. The advent of next generation sequencing provides a unique opportunity for a mitogenome-wide study. Here, we first provided six new complete mitogenomes of Certhioidea (Certhia americana, C. familiaris, Salpornis spilonota, Cantorchilus leucotis, Pheugopedius coraya, and Pheugopedius genibarbis). We further paid attention to the genomic characteristics, codon usages, evolutionary rates, and phylogeny of the Certhioidea mitogenomes. All mitogenomes we analyzed displayed typical ancestral avian gene order with 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and one control region (CR). Our study indicated the strand-biased compositional asymmetry might shape codon usage preferences in mitochondrial genes. In addition, natural selection might be the main factor in shaping the codon usages of genes. Additionally, evolutionary rate analyses indicated all mitochondrial genes were under purifying selection. Moreover, MT-ATP8 and MT-CO1 were the most rapidly evolving gene and conserved genes, respectively. According to our mitophylogenetic analyses, the monophylies of Troglodytidae and Sittidae were strongly supported. Importantly, we suggest that Salpornis should be separated from Certhiidae and put into Salpornithidae to maintain the monophyly of Certhiidae. Our findings are useful for further evolutionary studies within Certhioidea.

Identification and Genomic Characterization of Escherichia albertii in Migratory Birds from Poyang Lake, China

Escherichia albertii is an emerging zoonotic foodborne enteropathogen leading to human gastroenteritis outbreaks. Although E. albertii has been isolated from birds which have been considered as the potential reservoirs of this bacterium, its prevalence in migratory birds has rarely been described. In this study, E. albertii in migratory birds from Poyang Lake was investigated and characterized using whole genome sequencing. Eighty-one fecal samples from nine species of migratory birds were collected and 24/81 (29.6%) tested PCR-positive for E. albertii-specific genes. A total of 47 isolates was recovered from 18 out of 24 PCR-positive samples. All isolates carried eae and cdtB genes. These isolates were classified into eight E. albertii O-genotypes (EAOgs) (including three novel EAOgs) and three E. albertii H-genotypes (EAHgs). Whole genome phylogeny separated migratory bird-derived isolates into different lineages, some isolates in this study were phylogenetically closely grouped with poultry-derived or patient-derived strains. Our findings showed that migratory birds may serve as an important reservoir for heterogeneous E. albertii, thereby acting as potential transmission vehicles of E. albertii to humans.

2197. Genomic Surveillance for Human Adenovirus Infections among Department of Defense Beneficiaries in the 2018-2019 Season

Abstract Background Human adenoviruses (family Adenoviridae, genus Mastadenovirus) (HAdV) species B, C and E cause acute respiratory illnesses (ARI). In particular, respiratory HAdV types 3, 4, 7, 11, 14, 21 and 55 were associated with ARI outbreaks in the US and in other countries. The risk of outbreaks can be effectively controlled by implementation of the live HAdV types 4 and 7 vaccines or mitigated by screening and preventive measures. Methods During the influenza season 2018 – 2019, the DoD Global Respiratory Pathogen Surveillance Program (DoDGRS) received 24,300 specimens from its surveillance network. In total, 322 (1.3%) respiratory specimens were positive for HAdV in molecular assays. HAdV samples that produced positive cytopathic effects (CPE) in subsequent viral cultivation were subjected to next-gen sequencing (NGS) for genome sequence assembly, HAdV genome typing, whole genome phylogeny, and sequence comparative analyses to identify sequence insertions, deletions and amino acid mutations. Results Out of 166 viral culture samples available for NGS, whole genome sequences were obtained for 161 isolates. A variety types of HAdV were identified, including HAdV-1 (N = 15), HAdV-2 (N = 26), HAdV-3 (N = 60), HAdV-4 (N = 13), HAdV-5 (N = 6), HAdV-6 (N = 3), HAdV-7 (N = 12), HAdV-14 (N = 10), HAdV-55 (N = 1), and HAdV-56 (N = 1). HAdV types 4, 7 and 14, which were reported causing ARI outbreaks in the US military and other congregate settings, were found in clustered cases in several US states. Comparative sequence analyses of these isolates revealed the emergence of novel genetic mutations despite the stability of adenovirus genomes. Conclusion Routine clinical respiratory pathogen diagnostics or molecular surveillance of respiratory diseases detect adenovirus without determination of HAdV type. Human adenoviruses are genetically diverse and can infect a number of tissues with severities varied from mild to fatal. Enhanced genomic surveillance of HAdV in the US and worldwide will shed light on prevalence, genetic divergence, and viral evolution of human adenoviruses and inform timely risk assessment and countermeasures. Disclosures All Authors: No reported disclosures.

A comparative genomics approach for identifying genetic factors in Escherichia coli isolates associated with bovine diseases.

Escherichia coli are ubiquitously present bacterial pathogens that cause septicaemia, diarrhoea and other clinical illness in farm animals. Many pathogen factors can be associated with disease conditions. Currently, studies inferring E. coli genetic factors associated with infection in bovines are limited. Hence, the present study envisaged to determine the pathogen genetic factors associated with bovine disease conditions. The comparative genomic analysis involved genome sequence data of 135 diseased and 145 healthy bovine origin E. coli strains. Phylogroups A and C, as well as pathotypes ExPEC and EPEC, were found to have a strong connection with bovine disease strains. STEC strains, including EHEC, seem to play a less important role in bovine disease. Sequence types (STs) predominant among strains from diarrhoeal origin were ST 301 (CC 165) and ST 342. Correlation of core genome phylogeny with accessory gene based clustering, phylogroups and pathotypes indicated lineage specific virulence factors mostly associated with disease conditions. Comparative genomic analysis was applied to infer genetic factors significant in bovine disease origin E. coli strains. Isolates from bovine disease origin were enriched for the phylogroups A and C, and for the pathotypes ExPEC and EPEC. However, there was minimal evidence of STEC involvement. The study also indicated predominant genetic lineages and virulence genes (pap, sfa and afa) associated with disease origin strains. SIGNIFICANCE AND IMPACT OF STUDY;: The study revealed significant pathotypes, phylgroups, serotypes and sequence types associated with bovine disease conditions. These identified genetic factors can be applied for disease diagnosis, implementing vaccine and therapeutic measures. In addition, E. coli isolates from the bovine species revealed a complex pattern of disease epidemiology.

Within-Host Genotypic and Phenotypic Diversity of Contemporaneous Carbapenem-Resistant Klebsiella pneumoniae from Blood Cultures of Patients with Bacteremia.

It is unknown whether bacterial bloodstream infections (BSIs) are commonly caused by single organisms or mixed microbial populations. We hypothesized that contemporaneous carbapenem-resistant Klebsiella pneumoniae (CRKP) strains from blood cultures of individual patients are genetically and phenotypically distinct. We determined short-read whole-genome sequences of 10 sequence type 258 (ST258) CRKP strains from blood cultures in each of 6 patients (Illumina HiSeq). Strains clustered by patient by core genome and pan-genome phylogeny. In 5 patients, there was within-host strain diversity by gene mutations, presence/absence of antibiotic resistance or virulence genes, and/or plasmid content. Accessory gene phylogeny revealed strain diversity in all 6 patients. Strains from 3 patients underwent long-read sequencing for genome completion (Oxford Nanopore) and phenotypic testing. Genetically distinct strains within individuals exhibited significant differences in carbapenem and other antibiotic responses, capsular polysaccharide (CPS) production, mucoviscosity, and/or serum killing. In 2 patients, strains differed significantly in virulence during mouse BSIs. Genetic or phenotypic diversity was not observed among strains recovered from blood culture bottles seeded with index strains from the 3 patients and incubated in vitro at 37°C. In conclusion, we identified genotypic and phenotypic variant ST258 CRKP strains from blood cultures of individual patients with BSIs, which were not detected by the clinical laboratory or in seeded blood cultures. The data suggest a new paradigm of CRKP population diversity during BSIs, at least in some patients. If validated for BSIs caused by other bacteria, within-host microbial diversity may have implications for medical, microbiology, and infection prevention practices and for understanding antibiotic resistance and pathogenesis. IMPORTANCE The long-standing paradigm for pathogenesis of bacteremia is that, in most cases, a single organism passes through a bottleneck and establishes itself in the bloodstream (single-organism hypothesis). In keeping with this paradigm, standard practice in processing positive microbiologic cultures is to test single bacterial strains from morphologically distinct colonies. This study is the first genome-wide analysis of within-host diversity of Klebsiella pneumoniae strains recovered from individual patients with bloodstream infections (BSIs). Our finding that positive blood cultures comprised genetically and phenotypically heterogeneous carbapenem-resistant K. pneumoniae strains challenges the single-organism hypothesis and suggests that at least some BSIs are caused by mixed bacterial populations that are unrecognized by the clinical laboratory. The data support a model of pathogenesis in which pressures in vivo select for strain variants with particular antibiotic resistance or virulence attributes and raise questions about laboratory protocols and treatment decisions directed against single strains.

Revision of Quedius sensu stricto (Coleoptera: Staphylinidae)

We here present the first integrative revision of the subgenus Quedius Stephens, 1829 sensu stricto where taxonomic decisions are based on morphology, genomic phylogeny (published elsewhere) and single locus DNA evidence. The subgenus is restricted to the Holarctic region and includes some of the largest, most robust members of the genus Quedius. For species delimitation, a total of 200 COI barcodes covering nearly all species within the subgenus were evaluated through phylogenetic, cluster and network analysis. Taxonomic, distributional and bionomic data, all hitherto published and new are synthesized and an identification key is constructed for all species. All 20 species of Quedius s. str. are divided into five species groups based on genomic phylogeny, morphology, distributional patterns and practical considerations. Overall, there was good congruence among various lines of evidence. As exceptions, high COI barcode variability was found within the wingless and patchily distributed Q. unicolor Kiesenwetter, 1847 and Q. sundukovi Smetana, 2003 (3 and 5 BOLD BINs, respectively) from the Q. molochinus-group without corresponding morphological or geographic patterns. On the contrary, in the Q. pallipes-group very little divergence in the COI barcode was found between clearly morphologically separated species. From the taxonomically challenging Nearctic Quedius molochinoides-group, Quedius altanai Hansen & Brunke sp. nov. is described as new species from the central and southern Rockies, rendering Q. lanei Hatch, 1957 to be restricted to the eastern foothills of the Cascades, Sierra Nevada and Blue Mountains. In the case of Q. subunicolor and Q. altaicus, continuous morphological and COI variation firmly established the earlier suspected synonymy Quedius subunicolor Korge, 1961 = Quedius altaicus Korge, 1962 syn. nov. Other new synonyms established here are: Quedius pallipes Lucas, 1849 = Quedius simplicifrons Fairmaire, 1861 syn. nov. = Quedius levasseuri Coiffait, 1964 syn. nov.; Quedius hispanicus Bernhauer, 1898 stat. reinstituted = Quedius cobosi Coiffait, 1964 syn. nov.; Quedius fuliginosus (Gravenhorst, 1802) = Quedius latus Hochhuth, 1851 syn. nov. = Quedius viduus Sawada, 1965, syn. nov. Quedius gracilis Stephens, 1832 syn. rev. is moved from synonymy with Q. fuliginosus and placed into synonymy with Q. curtipennis Bernhauer, 1908, the younger being the valid name due to prevailing usage in accordance with the ICZN Article 23.9; Quedius hammianus Sharp, 1911 syn. revid., Q. secundus Last, 1952 syn. revid., and Q. rufulus Blümml, 1898 syn. revid., are moved from synonymy with Q. simplicifrons and placed into synonymy with Q. hispanicus Bernhauer, 1898; Quedius sardous Gridelli, 1924 syn. revid., and Q. leonhardi Bernhauer, 1914 syn. revid., are moved from synonymy with Q. molochinus (Gravenhorst, 1806) to synonymy with Q. pallipes. Human mediated dispersals, most likely through the historical transport of ship ballast, are noted as regular phenomenon for some members of the subgenus. Quedius fuliginosus, Q. curtipennis and Q. molochinus, introduced from the Palearctic to the Nearctic, are the most noticeable examples.

Paradox versus paradigm: A disconnect between understanding and management of freshwater cyanobacterial harmful algae blooms

AbstractFreshwater cyanobacterial harmful algae blooms (cHABs) are a major threat to human and environmental health and are increasing globally in frequency and severity. To manage this threat in a timely manner, science must focus on increasing our ability to predict the growth and toxigenicity of specific taxa of cyanobacteria.Recent molecular research has revealed striking genomic and metabolic diversity among the many morphologically indistinguishable sub‐species and strains of cyanobacteria. Assemblage‐level molecular metabolic capability surveys promise to improve our ability to predict cyanobacterial responses to environmental forcing, but many of these cutting‐edge techniques are not widely available or cost‐effective enough to be employed in routine monitoring programmes to support management decisions. Taxonomic ambiguity, cryptic functional specialisation, incongruence between genomic capability and phylogeny, and genomic flexibility impose severe challenges to our ability to ascribe autecological attributes at a level of taxonomic resolution that is attainable under current management strategies (i.e. Linnaean species). This lack of knowledge prohibits reliable predictions of species' responses to environmental stressors.Cyanobacterial species comprise consortia of metabolically diverse, morphologically indistinct strains that span a range of ecological specialisation. Under current, broadly applied taxonomic concepts, these species functionally embody a generalist ecological strategy—persisting and/or proliferating where other specialised competitors are negatively impacted.We postulate that within current management frameworks, characterising of cyanobacterial species as competing generalists, as well as considering abundance trajectories of well‐characterised, non‐cyanobacterial specialist phytoplankton will generate more scalable, mechanistic, and management‐relevant insight into increasing cHAB frequency and severity in suitable time frames.Here we recommend that cHAB management considers the competitive framework of phytoplankton communities, including cyanobacteria, wherein diverse environmental changes lead to deterministic responses by readily identifiable, documented specialist taxa. Characterising these changes in community structure will quantify the relative importance of altered stressors and resource availability that can be exploited by metabolically flexible cyanobacteria.

Whole mitochondrial genome phylogeny of Drosophilidae

A total of 241 mitochondrial genomes were assembled and annotated from the SRA database to reconstruct a mtDNA genome phylogeny for the genus Drosophila, the family Drosophilidae, and close relatives. The resulting mtDNA genome phylogeny is largely congruent with previous higher-level analyses of Drosophila species with the exception of the relationships between the melanogaster, montium, anannassae, saltans and obscura groups. Although relationships within these species groups are congruent between nuclear and mtDNA studies, the mtDNA genome phylogeny of the groups is different when compared to earlier studies. Monophyly of known species groups within the genus Drosophila are highly supported and, as in previous work, the genera Lordiphosa, Hirtodrosophila, Zaprionus and Scaptomya are all imbedded within the genus Drosophila. Incongruence and partitioned support analyses indicate that DNA sequences are better at resolving the phylogeny than their translated protein sequences. Such analyses also indicate that genes on the minus strand of the circular molecule (Lrrna, Srrna, ND4, ND4L and ND5) provide most of the support for the overall phylogenetic hypothesis.

Phylogenetic analyses of Salmonella detected along the broiler production chain in Trinidad and Tobago

This study was conducted to determine the phylogenies of Salmonella strains isolated from cross-sectional studies conducted at hatcheries, broiler farms, processing plants, and retail outlets (broiler production chain) in Trinidad and Tobago over 4 yr (2016-2019). Whole-genome sequencing (WGS) was used to characterize Salmonella isolates. Core genome phylogenies of 8 serovars of public health significance were analyzed for similarities in origin and relatedness. In addition, Salmonella strains isolated from human salmonellosis cases in Trinidad were analyzed for their relatedness to the isolates detected along the broiler production chain. The common source of these isolates of diverse serovars within farms, within processing plants, between processing plants and retail outlets, and among farm-processing plant-retail outlet continuum was well-supported (bootstrap value >70%) by the core genome phylogenies for the respective serovars. Also, genome analyses revealed clustering of Salmonella serovars of regional (intra-Caribbean) and international (extra-Caribbean) origin. Similarly, strains of S. Enteritidis and S. Infantis isolated from human clinical salmonellosis in 2019 from Trinidad and Tobago clustered with our processing plant isolates recovered in 2018. This study is the first phylogenetic analysis of Salmonella isolates using WGS from the broiler industry in the Caribbean region. The use of WGS confirmed the genetic relatedness and transmission of Salmonella serovars contaminating chickens in broiler processing, and retailing in the country, with zoonotic and food safety implications for humans.

Seventy-eight entire mitochondrial genomes and nuclear rRNA genes provide insight into the phylogeny of the hard ticks, particularly the Haemaphysalis species, Africaniella transversale and Robertsicus elaphensis

Hoogstraal and Kim (1985) proposed from morphology, three groups of Haemaphysalis subgenera: (i) the “structurally advanced”; (ii) the “structurally intermediate”; and (iii) the “structurally primitive” subgenera. Nuclear gene phylogenies, however, did not indicate monophyly of these morphological groups but alas, only two mitochondrial (mt) genomes from the “structurally intermediate” subgenera had been sequenced. The phylogeny of Haemaphysalis has not yet been resolved. We aimed to resolve the phylogeny of the genus Haemaphysalis, with respect to the subgenus Alloceraea. We presented 15 newly sequenced and annotated mt genomes from 15 species of ticks, five species of which have not been sequenced before, and four new 18S rRNA and 28S rRNA nuclear gene sequences. Our datasets were constructed from 10 mt protein-coding genes, cox1, and the 18S and 28S nuclear rRNA genes. We found a 132-bp insertion between tRNA-Glu (E) gene and the nad1 gene in the mt genome of Haemaphysalis (Alloceraea) inermis that resembles insertions in H. (Alloceraea) kitaokai and Rhipicephalus (Boophilus) geigyi. Our mt phylogenies had the three species of Amblyomma (Aponomma) we sequenced embedded in the main clade of Amblyomma: Am. (Aponomma) fimbriatum, Am. (Aponomma) gervaisi and Am. (Aponomma) latum. This is further support for the hypothesis that the evolution of eyes appears to have occurred in the most-recent-common-ancestor of Amblyocephalus (i.e. Amblyomminae plus Rhipicephalinae) and that eyes were subsequently lost in the most-recent-common-ancestor of the subgenus Am. (Aponomma). Either Africaniella transversale or Robertsicus elaphensis, or perhaps Af. transversale plus Ro. elaphensis, appear to be the sister-group to the rest of the metastriate Ixodida. Our cox1 phylogenies did not indicate monophyly of the “structurally primitive”, “structurally intermediate” nor the “structurally advanced” groups of Haemaphysalis subgenera. Indeed, the subgenus Alloceraea may be the only monophyletic subgenus of the genus Haemaphysalis sequenced thus far. All of our mt genome and cox1 phylogenies had the subgenus Alloceraea in a clade that was separate from the rest of the Haemaphysalis ticks. If Alloceraea is indeed the sister to the rest of the Haemaphysalis subgenera this would resonate with the argument of Hoogstraal and Kim (1985), that Alloceraea was a subgenus of “primitive” Haemaphysalis. Alectorobius capensis from Japan had a higher genetic-identity to A. sawaii, which was also from Japan, than to the A. capensis from South Africa. This indicates that A. capensis from Japan may be a cryptic species with respect to the A. capensis from South Africa.

Identification and genomic analysis of a Vibrio cholerae strain isolated from a patient with bloodstream infection

Vibrio cholerae is a bacterium ubiquitous in aquatic environments which can cause widespread infection worldwide. V. cholerae gradually became a rare species of bacteria in clinical microbiology laboratories with the control of the cholera epidemic. In this study, we isolated a V. cholerae strain, named VCHL017, from the blood of an elderly patient without gastrointestinal symptoms. The patient had a history of hookworm infection and multiple myeloma. Furthermore, she was immunocompromised, and received long-term chemotherapy and antimicrobial agents. VCHL017 was inoculated on blood agar and thiosulfate citrate bile salt sucrose plates (TCBS) to observe morphological characteristics. Then this isolate was identified by matrix-assisted laser desorption/ionization time-of-flight spectrometry (MALDI-TOF MS). The minimum inhibitory concentrations (MICs) for cefazolin, ceftazidime, cefepime, meropenem, tetracycline, ciprofloxacin, chloramphenicol, and gentamicin of VCHL017 were determined by the microbroth dilution method. PCR and serum agglutination tests were used to determine whether the serogroups of the isolate belonged to the O1/O139 and cholera toxin encoding genes. Finally, the genomic features and phylogeny of VCHL017 were analyzed by whole genome sequencing (WGS). VCHL017 was a non-O1/O139 V cholerae strain that did not carry the ctxA gene. Antimicrobial susceptibility tests revealed that VCHL017 was susceptive to chloramphenicol and tetracycline. Although it did not carry the genes encoding the cholera toxin, WGS indicated that VCHL017 carried a variety of other virulence factors. By calculating the average nucleotide identity (ANI), we precisely identified the species of VCHL017 as V. cholerae. There are also A171S and A202S missense mutations in gyrA of VCHL017. The phylogenetic analysis indicated that VCHL017 was closely related to V. cholerae strains isolated from aquatic environments. Our results suggest that continuous monitoring is necessary for non-O1/O139 V cholerae strains isolated from outside the digestive tract, which could be pathogenic through multiple virulence factors.

Recombination in sarbecovirus lineage and mutations/insertions in spike protein are linked to the emergence and adaptation of SARS-CoV-2.

The outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Wuhan city, China in December 2019 and thereafter its spillover across the world has created a global pandemic and public health crisis. Right after, there has been intense interest in understanding how the SARS-CoV-2 originated and evolved. This paper also aims to shed light on the origin and evolution of SARS-CoV- 2. A consensus result based on whole genome phylogeny, gene tree analysis, and genetic similarity study revealed that SARS-CoV-2 evolved from Bat-CoV-RaTG13. Furthermore, recombination analysis indicated that probable origin of SARS-CoV-2 is the results of ancestral intra-species recombination events between bat coronaviruses belonging to Sarbecovirus sub-genus. Multiple sequence alignment (MSA) revealed the insertion of four amino acid residues "PRRA" (Proline-Arginine-Arginine-Alanine) to the S1/S2 site in the spike protein of SARS-CoV-2, and structural modeling of spike protein of bat-CoV-RaTG13 also shows a high number of mutations at one of the receptor binding domains (RBD). Acquisition of the furin cleavage sites ("PRRA") along with high number of mutations at one of its RBD is probably responsible for the adaptation of SARS-CoV-2 into human systems. Furthermore, the codon adaptation index (CAI) was used to quantify the magnitude of adaptive efficacy of SARS-CoV-2 in human host in comparison with SARS-CoV. The CAI result showed a relatively less adaptive efficacy of the newly emerged SARS-CoV-2 to the human systems, which might be an indication of its mild clinical severity and progression compared to SARS-CoVs.

Identification, comprehensive characterization, and comparative genomics of the HERV-K(HML8) integrations in the human genome

Around 8% of the human genome is composed by Human Endogenous Retroviruses (HERVs), ancient viral sequences inherited from the primate germ line after their infection by now extinct retroviruses. Given the still underexplored physiological and pathological roles of HERVs, it is fundamental to increase our information about the genomic composition of the different groups, to lay reliable foundation for functional studies. Among HERVs, the most characterized elements belong to the beta-like superfamily HERV-K, comprising 10 groups (HML1-10) with HML2 being the most recent and studied one. Among HMLs, the HML8 group is the only one still lacking a comprehensive genomic description. In the present work, we investigated HML8 sequences’ distribution in the human genome (GRCh38/hg38), identifying 23 novel proviruses and characterizing the overall 78 HML8 proviruses in terms of genome structure, phylogeny, and integration pattern. HML8 elements were significantly enriched in human chromosomes 8 and X (p<0.005) while chromosomes 17 and 20 showed fewer integrations than expected (p<0.025 and p<0.005, respectively). Phylogenetic analyses classified HML8 members into 3 clusters, corresponding to the three LTR types MER11A, MER11B and MER11C. Besides different LTR types, common signatures in the internal structure suggested the potential existence of three different ancestral HML8 variants. Accordingly, time of integration estimation coupled with comparative genomics revealed that these three clusters have a different time of integration in the primates’ genome, with MER11C elements being significantly younger than MER11A- and MER11B associated proviruses (p<0.005 and p<0.05, respectively). Approximately 30% of the HML8 elements were found co-localized within human genes, sometimes in exonic portions and with the same orientation, deserving further studies for their possible effects on gene expression. Overall, we provide the first detailed picture of the HML8 group distribution and variety among the genome, creating the backbone for the specific analysis of their transcriptional activity in healthy and diseased conditions.