Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation

Genome-Annotated Bacterial Collection of the Barley Rhizosphere Microbiota

A culture collection of 41 bacteria isolated from the rhizosphere of cultivated barley (Hordeum vulgare subsp. vulgare) is available at the Division of Plant Sciences, University of Dundee (UK). The data include information on genes putatively implicated in nitrogen fixation, HCN channels, phosphate solubilization, and linked whole-genome sequences.

Comparative genome analyses of five Vibrio penaeicida strains provide insights into their virulence-related factors

Vibrio penaeicida (family Vibrionaceae) is an important bacterial pathogen that affects Japanese shrimp aquaculture. Only two whole-genome sequences of V. penaeicida are publicly available, which has hampered our understanding of the pathogenesis of shrimp vibriosis caused by this bacterium. To gain insight into the genetic features, evolution and pathogenicity of V. penaeicida, we sequenced five V. penaeicida strains (IFO 15640^T, IFO 15641, IFO 15642, TUMSAT-OK1 and TUMSAT-OK2) and performed comparative genomic analyses. Virulence factors and mobile genetic elements were detected. Furthermore, average nucleotide identities (ANIs), clusters of orthologous groups and phylogenetic relationships were evaluated. The V. penaeicida genome consists of two circular chromosomes. Chromosome I sizes ranged from 4.1 to 4.3 Mb, the GC content ranged from 43.9 to 44.1 %, and the number of predicted protein-coding sequences (CDSs) ranged from 3620 to 3782. Chromosome II sizes ranged from 2.2 to 2.4 Mb, the GC content ranged from 43.5 to 43.8 %, and the number of predicted CDSs ranged from 1992 to 2273. All strains except IFO 15641 harboured one plasmid, having sizes that ranged from 150 to 285 kb. All five genomes had typical virulence factors, including adherence, anti-phagocytosis, flagella-related proteins and toxins (repeats-in-toxin and thermolabile haemolysin). The genomes also contained factors responsible for iron uptake and the type II, IV and VI secretion systems. The genome of strain TUMSAT-OK2 tended to encode more prophage regions than the other strains, whereas the genome of strain IFO 15640^T had the highest number of regions encoding genomic islands. For comparative genome analysis, we used V. penaeicida (strain CAIM 285^T) as a reference strain. ANIs between strain CAIM 285^T and the five V. penaeicida strains were >95 %, which indicated that these strains belong to the same species. Orthology cluster analysis showed that strains TUMSAT-OK1 and TUMSAT-OK2 had the greatest number of shared gene clusters, followed by strains CAIM 285^T and IFO 15640^T. These strains were also the most closely related to each other in a phylogenetic analysis. This study presents the first comparative genome analysis of V. penaeicida and these results will be useful for understanding the pathogenesis of this bacterium.

Genetic Characterization of Multidrug-Resistant E. coli Isolates from Bloodstream Infections in Lithuania

Extraintestinal pathogenic Escherichia coli (ExPEC) isolates are a main cause of bloodstream infections. The aim of this study was to characterize 256 β-lactam–resistant, bacteremia-causing E. coli isolates collected from 12 healthcare institutions in Lithuania in 2014 and 2018. All isolates were identified as E. coli via MALDI-TOF MS and VITEK ®2. In addition, the isolates were analyzed for the presence of 29 resistance genes and 13 virulence genes, divided into phylogenetic groups (A, B1, B2, C, D, E, and F), and characterized using rep-PCR genotyping methods (BOX-PCR and (GTG)5-PCR). Analyzing the results of this study showed tetA-strB-sul2-TEM-NDM-strA-fosA-AIM-sul3-aadA-CTX-M-9 to be the most common resistance gene combination (67.2% of all isolates). Additionally, the most common virulence genes established were fimH (98.4% of all isolates), fyuA (91.8%), and traT (81.3%) and the most common gene combination was fuyA-fimH-iroN (58.6% of all isolates). Next, the isolates were separated into four phylogenetic groups: A, B1, B2, and F, where group A isolates were detected at a significantly higher frequency (79.3% of all isolates). Finally, a total of 235 genotyping profiles were established using rep-PCR methods, and all profiles were separated into fourteen genotypic clusters, with each cluster containing profiles with a variety of virulence and resistance genes not restricted to any specific cluster. The results of this study elucidate E. coli antimicrobial resistance patterns by highlighting the variability and diversity of resistance and virulence genes and providing phylogenetic classification, genetic profiling, and clustering data. These results may improve clinical control of multidrug-resistant infections in healthcare institutions and contribute to the prevention of potential outbreaks.

Integrative Genomics Sheds Light on Evolutionary Forces Shaping the Acidithiobacillia Class Acidophilic Lifestyle

Extreme acidophiles thrive in environments rich in protons (pH values <3) and often high levels of dissolved heavy metals. They are distributed across the three domains of the Tree of Life including members of the Proteobacteria. The Acidithiobacillia class is formed by the neutrophilic genus Thermithiobacillus along with the extremely acidophilic genera Fervidacidithiobacillus, Igneacidithiobacillus, Ambacidithiobacillus, and Acidithiobacillus. Phylogenomic reconstruction revealed a division in the Acidithiobacillia class correlating with the different pH optima that suggested that the acidophilic genera evolved from an ancestral neutrophile within the Acidithiobacillia. Genes and mechanisms denominated as "first line of defense" were key to explaining the Acidithiobacillia acidophilic lifestyle including preventing proton influx that allows the cell to maintain a near-neutral cytoplasmic pH and differ from the neutrophilic Acidithiobacillia ancestors that lacked these systems. Additional differences between the neutrophilic and acidophilic Acidithiobacillia included the higher number of gene copies in the acidophilic genera coding for "second line of defense" systems that neutralize and/or expel protons from cell. Gain of genes such as hopanoid biosynthesis involved in membrane stabilization at low pH and the functional redundancy for generating an internal positive membrane potential revealed the transition from neutrophilic properties to a new acidophilic lifestyle by shaping the Acidithiobacillaceae genomic structure. The presence of a pool of accessory genes with functional redundancy provides the opportunity to "hedge bet" in rapidly changing acidic environments. Although a core of mechanisms for acid resistance was inherited vertically from an inferred neutrophilic ancestor, the majority of mechanisms, especially those potentially involved in resistance to extremely low pH, were obtained from other extreme acidophiles by horizontal gene transfer (HGT) events.

The Phylosymbiosis Pattern Between the Fig Wasps of the Same Genus and Their Associated Microbiota

Microbial communities can be critical for many metazoans, which can lead to the observation of phylosymbiosis with phylogenetically related species sharing similar microbial communities. Most of the previous studies on phylosymbiosis were conducted across the host families or genera. However, it is unclear whether the phylosymbiosis signal is still prevalent at lower taxonomic levels. In this study, 54 individuals from six species of the fig wasp genus Ceratosolen (Hymenoptera: Agaonidae) collected from nine natural populations and their associated microbiota were investigated. The fig wasp species were morphologically identified and further determined by mitochondrial CO1 gene fragments and nuclear ITS2 sequences, and the V4 region of 16S rRNA gene was sequenced to analyze the bacterial communities. The results suggest a significant positive correlation between host genetic characteristics and microbial diversity characteristics, indicating the phylosymbiosis signal between the phylogeny of insect hosts and the associated microbiota in the lower classification level within a genus. Moreover, we found that the endosymbiotic Wolbachia carried by fig wasps led to a decrease in bacterial diversity of host-associated microbial communities. This study contributes to our understanding of the role of host phylogeny, as well as the role of endosymbionts in shaping the host-associated microbial community.

Discovery and Characterization of a New Crustin Antimicrobial Peptide from Amphibalanus amphitrite

Crustins are an antimicrobial peptide (AMP) family that plays an important role in innate immunity in crustaceans. It is important to discover new AMPs from natural sources to expand the current database. Here, we identified and characterized a new crustin family member, named AaCrus1, from Amphibalanus amphitrite. AaCrus1 shares high identity (48.10%) with PvCrus, a Type I crustin of Penaeus vannamei that possesses a whey acidic protein (WAP) domain. AaCrus1 contains 237 amino acids and eight cysteine residues forming conserved 'four-disulfide core' structure. Our recombinant AaCrus1 (rAaCrus 1) could inhibit the growth of two Gram-positive bacteria (Staphylococcus aureus, Bacillus sp. T2) and four Gram-negative bacteria (Vibrio parahaemolyticus, Vibrio harveyi, Vibrio anguillarum, Vibrio alginolyticus) with a minimum inhibitory concentration of 3.5-28 μM. It can further induce agglutination of both Gram-positive and Gram-negative bacteria. rAaCrus1 can bind to bacteria and damage bacterial cell membranes. Furthermore, rAaCrus1 disrupted biofilm development of S. aureus and V. parahaemolyticus. Our discovery and characterization of this new crustin can be further optimized as a good alternative to antibiotics.

Macrococcus armenti sp. nov., a novel bacterium isolated from the skin and nasal cavities of healthy pigs and calves

Gram-positive coccoid bacteria were isolated from the nasal cavities of pigs and calves as well as from axillar and inguinal skin regions of pigs. Phylogenetic analysis of seven strains based on complete genome, 16S rRNA, hsp60, dnaJ, rpoB and sodA gene sequences and MALDI-TOF MS profiles revealed that they belonged to the genus Macrococcus with the closest relatedness to Macrococcus canis , Macrococcus caseolyticus subsp. caseolyticus and Macrococcus caseolyticus subsp. hominis . DNA relatedness of the type strain JEK37T with the type strains of M. canis , M. caseolyticus subsp. caseolyticus and M. caseolyticus subsp. hominis was 23.4, 23.1 and 23.0 % by digital DNA–DNA hybridization and 80.39, 80.45 and 80.87 % by average nucleotide identity (ANI) calculations, confirming that they do not belong to the same species. The DNA G+C content of JEK37T was 35.65 mol%. The novel strains can be differentiated from M. canis KM 45013T by the ability to fermentate d-ribose and by the absence of DNAase production and haemolysis, from M. caseolyticus subsp. caseolyticus CCUG 15606T by the ability to fermentate sucrose and from both species by the inability to grow in 9 and 12% NaCl. They differ from M. caseolyiticus subsp. hominis by the presence of α-glucosidase. The most common fatty acids of JEK37T were C14 : 0, C18 : 1 ω9c and C18 : 0. Known polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol, aminolipid, aminoglycolipid, aminophospholipid, glycolipid and phospholipid. Cell-wall peptidoglycan of JEK37T was of type A3α l-Lys–Gly2–L-Ser–Gly (similar to A11.3) and the respiratory quinolone was menaquinone 6. Based on their genotypic and chemotaxonomic characteristics, these strains represent a novel species of the genus Macrococcus , for which we propose the name Macrococcus armenti sp. nov. The type strain is JEK37T (=DSM 112712T=CCOS 1982T).

A holistic genome dataset of bacteria, archaea and viruses of the Pearl River estuary

Estuaries are one of the most important coastal ecosystems. While microbiomes and viromes have been separately investigated in some estuaries, few studies holistically deciphered the genomes and connections of viruses and their microbial hosts along an estuarine salinity gradient. Here we applied deep metagenomic sequencing on microbial and viral communities in surface waters of the Pearl River estuary, one of China’s largest estuaries with strong anthropogenic impacts. Overall, 1,205 non-redundant prokaryotic genomes with ≥50% completeness and ≤10% contamination, and 78,502 non-redundant viral-like genomes were generated from samples of three size fractions and five salinity levels. Phylogenomic analysis and taxonomy classification show that majority of these estuarine prokaryotic and viral genomes are novel at species level according to public databases. Potential connections between the microbial and viral populations were further investigated by host-virus matching. These combined microbial and viral genomes provide an important complement of global marine genome datasets and should greatly facilitate our understanding of microbe-virus interactions, evolution and their implications in estuarine ecosystems.

Measurement(s)	bacteria • Archaea • viruses
Technology Type(s)	Shotgun Sequencing
Sample Characteristic - Organism	estuary metagenome
Sample Characteristic - Environment	subtropical estuarine
Sample Characteristic - Location	Pearl river estuary

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.17139998

Pseudomonas bijieensis Strain XL17 within the P. corrugata Subgroup Producing 2,4-Diacetylphloroglucinol and Lipopeptides Controls Bacterial Canker and Gray Mold Pathogens of Kiwifruit

Kiwifruit worldwide suffers from the devastating diseases of bacterial canker caused by Pseudomonas syringae pv. actinidiae (Psa) and gray mold caused by Botrytis cinerea. Here, an endophytic bacterium XL17 isolated from a rape crown gall was screened out for its potent antagonistic activities against Psa and B. cinerea. Strain XL17 and its cell-free culture filtrate (CF) inhibited the growth of Psa and B. cinerea, Psa-associated leaf necrosis, and B. cinerea-associated kiwifruit necrosis. Electron microscopy showed that XL17 CF could damage the cell structures of Psa and B. cinerea. Genome-based taxonomy revealed that strain XL17 belongs to Pseudomonas bijieensis within the P. corrugata subgroup of the P. fluorescens species complex. Among the P. corrugata subgroup containing 31 genomospecies, the presence of the phl operon responsible for the biosynthesis of the phenolic polyketide 2,4-diacetylphloroglucinol (DAPG) and the absence of the lipopeptide/quorum sensing island can serve as the genetic marker for the determination of a plant-protection life style. HPLC detected DAPG in extracts from XL17 CF. MALDI-TOF-MS analysis revealed that strain XL17 produced cyclic lipopeptides of the viscosin family and orfamide family. Together, phenotypic, genomic, and metabolic analyses identified that P. bijieensis XL17 producing DAPG and cyclic lipopeptides can be used to control bacterial canker and gray mold pathogens of kiwifruit.

Comprehensive Analysis of Imipenemase (IMP)-Type Metallo-β-Lactamase: A Global Distribution Threatening Asia

Antibiotic resistance, particularly beta-lactam resistance, is a major problem worldwide. Imipenemase or IMP-type metallo-β-lactamase (MBL) has become a more prominent enzyme, especially in Asia, since it was discovered in the 1990s in Japan. There are currently 88 variants of IMP-type enzymes. The most commonly identified variant of IMP-type enzymes is IMP−1 variant. IMP-type MBLs have been detected in more than ten species in Enterobacterales. Pseudomonas aeruginosa is the most frequent carrier of IMP-type enzymes worldwide. In Asia, IMP-type MBLs have been distributed in many countries. This work investigated a variety of currently available IMP-type MBLs at both a global level and a regional level. Out of 88 variants of IMP-type MBLs reported worldwide, only 32 variants were found to have susceptibility profiles. Most of the bacterial isolates carrying IMP-type MBLs were resistant to Carbapenems, especially Imipenem and Meropenem, followed by the 3rd-generation cephalosporins, and interestingly, monobactams. Our results comprehensively indicated the distribution of IMP-type MBLs in Asia and raised the awareness of the situation of antimicrobial resistance in the region.

Curtobacterium spp. and Curtobacterium flaccumfaciens: Phylogeny, Genomics-Based Taxonomy, Pathogenicity, and Diagnostics

The genus of Curtobacterium, belonging to the Microbacteriaceae family of the Actinomycetales order, includes economically significant pathogenic bacteria of soybeans and other agricultural crops. Thorough phylogenetic and full-genome analysis using the latest genomic data has demonstrated a complex and contradictory taxonomic picture within the group of organisms classified as the Curtobacterium species. Based on these data, it is possible to delineate about 50 new species and to reclassify a substantial part of the Curtobacterium strains. It is suggested that 53 strains, including most of the Curtobacterium flaccumfaciens pathovars, can compose a monophyletic group classified as C. flaccumfaciens. A genomic analysis using the most recent inventory of bacterial chromosomal and plasmid genomes deposited to GenBank confirmed the possible role of Microbacteriaceae plasmids in pathogenicity and demonstrated the existence of a group of related plasmids carrying virulence factors and possessing a gene distantly related to DNA polymerase found in bacteriophages and archaeal and eukaryotic viruses. A PCR diagnostic assay specific to the genus Curtobacterium was developed and tested. The presented results assist in the understanding of the evolutionary relations within the genus and can lay the foundation for further taxonomic updates.

Genome-Wide Analysis of U-box E3 Ubiquitin Ligase Family in Response to ABA Treatment in Salvia miltiorrhiza

Plant U-box (PUB) proteins are ubiquitin ligases (E3) involved in multiple biological processes and in response to plant stress. However, the various aspects of the genome and the differences in functions between the U-box E3 (UBE3) ubiquitin ligases remain quite obscure in Salvia miltiorrhiza. The 60 UBE3 genes in the S. miltiorrhiza genome were recognized in the present study. The phylogenetic analysis, gene structure, motifs, promoters, and physical and chemical properties of the genes were also examined. Based on the phylogenetic relationship, the 60 UBE3 genes were categorized under six different groups. The U-box domain was highly conserved across the family of UBE3 genes. Analysis of the cis-acting element revealed that the UBE3 genes might play an important role in a variety of biological processes, including a reaction to the abscisic acid (ABA) treatment. To investigate this hypothesis, an ABA treatment was developed for the hairy roots of S. miltiorrhiza. Thirteen out of the UBE3 genes significantly increased after the ABA treatment. The co-expression network revealed that nine UBE3 genes might be associated with phenolic acids or tanshinone biosynthesis. The findings of the present study brought fresh and new understanding to the participation of the UBE3 gene family in plants, specifically in their biological responses mediated by the ABA. In S. miltiorrhiza, this gene family may be crucial during the ABA treatment. Significantly, the results of this study contribute novel information to the understanding of the ubiquitin ligase gene and its role in plant growth.

Specificity of mouse and human Fcgamma receptors and their polymorphic variants for IgG subclasses of different species

Immunoglobulin G (IgG) is the predominant antibody class generated during infections and used for the generation of therapeutic antibodies. Antibodies are mainly characterized in or generated from animal models that support particular infections, respond to particular antigens or allow the generation of hybridomas. Due to the availability of numerous transgenic mouse models and the ease of performing bioassays with human blood cells in vitro, most antibodies from species other than mice and humans are tested in vitro using human cells and/or in vivo using mice. In this process, it is expected, but not yet systematically documented, that IgG from these species interact with human or mouse IgG receptors (FcγRs). In this study, we undertook a systematic assessment of binding specificities of IgG from various species to the families of mouse and human FcγRs, including their polymorphic variants. Our results document the specific binding patterns for each of these IgG (sub)classes, reveal possible caveats of antibody-based immunoassays, and will be a useful reference for the transition from one animal model to preclinical mouse models or human cell-based bioassays. This article is protected by copyright. All rights reserved.

Identification of Olfactory Genes in Monochamus saltuarius and Effects of Bursaphelenchus xylophilus Infestation on Their Expression

The pine wood nematode (PWN) Bursaphelenchus xylophilus has caused disastrous losses of pine forests in many countries, and the success of PWN depends strongly on interactions with its insect vectors. Monochamus saltuarius is a newly recorded vector in Northeast China. Feeding (i.e., immature) and egg-laying (i.e., mature) Monochamus spp. target different host plants, and olfactory cues play important roles regarding host choice. Whether infestation with PWN affects olfactory mechanisms in M. saltuarius related to feeding and oviposition is of interest as this may affect the spread of nematodes to new healthy hosts. However, little is known about molecular mechanisms of the olfactory system of M. saltuarius. We identified chemosensory-related genes in adult M. saltuarius and examined the influence of B. xylophilus on the respective expression patterns. Fifty-three odorant-binding proteins (OBPs), 15 chemosensory proteins, 15 olfactory receptors (ORs), 10 gustatory receptors, 22 ionotropic receptors (IRs), and two sensory neuron membrane proteins were identified, and sex bias among non-infested beetles was mainly found with respect to expression of OBPs. Interestingly, OBPs and ORs were markedly down-regulated in male M. saltuarius infested with B. xylophilus, which may reduce olfactory sensitivity of male M. saltuarius and affect the spreading of B. xylophilus to new hosts. Our results will help understand the interactions between B. xylophilus and M. saltuarius, which may lead to the identification of new control targets in the olfactory system of M. saltuarius.

Pseudonocardia humida sp. nov., an Actinomycete Isolated from Mangrove Soil Showing Distinct Distribution Pattern of Biosynthetic Gene Clusters

A novel actinomycete strain, designated S2-4^T, was isolated from a mangrove soil sample, and a polyphasic approach was employed to determine its taxonomic position. Phylogenetic analysis based on 16S rRNA gene indicated that strain S2-4^T formed a unique clade next to that harboring Pseudonocardia dioxanivorans CB1190^T, which shared the highest sequence similarity (98.20%) with the new isolate. Phylogenetic analysis based on core genes of genomic sequences displayed a different scenario, exhibiting closer phylogenetic relationship of strain S2-4^T with several species with 16S rRNA gene sequence similarities ranging from 96.95 to 98.06%, which was confirmed by the phylogenetic tree reconstructed based on genomic sequences. Further, substantial differences between the genotypic properties of strain S2-4^T and its closest neighbors, including digital DNA-DNA hybridization, average nucleotide identity, and distribution patterns of biosynthetic gene clusters (BGC), indicated the taxonomic position of strain S2-4^T as a novel species of the genus Pseudonocardia. Accordingly, strain S2-4^T was observed to show a different distribution pattern of a predicted BGC encoding ectoine by comparative genomic analysis, which could be strongly linked to its unique habitat distinct from where its close neighbors were isolated. The major cellular fatty acids were iso-C_15:0, C_21:0, and iso-C_16:0. The predominant menaquinone was MK-8(H₄). The polar lipids were composed of diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, phosphatidyl-N-monomethylethanolamine, phosphatidylcholine, phosphatidylinositol, phosphatidylinositol mannosides, and two unidentified glycolipids. Here, we propose a novel species of the genus Pseudonocardia: Pseudonocardia humida sp. nov. with the type strain S2-4^T (= JCM 34291^T = CGMCC 4.7706^T).

Polymerase chain reaction and loop-mediated isothermal amplification targeting lic13162, lic20239, and lipL32 genes for leptospirosis diagnosis

Leptospirosis is a zoonotic disease caused by pathogenic species of Leptospira. Due to the similarity with clinical signs of other febrile diseases, early diagnosis remains challenging. Real-time PCR has been used for direct detection of Leptospira, but it requires thermocyclers and highly trained personnel. Loop-mediated isothermal amplification (LAMP) is a simple and rapid DNA-based assay. Therefore, here we have developed PCR and LAMP assays targeting two novel genes, lic13162 and lic20239, and also lipL32 gene to detect pathogenic Leptospira. Analytical and diagnostic performances were compared with bacterial isolates (including different Leptospira species and serovars) and clinical samples. The results demonstrated that PCR assays targeting lic13162 and lic20239 were successful to amplify Leptospira, but LAMP not. However, both PCR and LAMP targeting lipL32 could detect pathogenic Leptospira. LAMP lipL32 could be performed in 30 min with a detection limit of 156 cells/mL. Diagnostic performance of lipL32-LAMP presented 84.2% sensitivity and 93.2% specificity. In conclusion, lipL32 PCR and LAMP are effective methods to detect pathogenic Leptospira directly from clinical samples.

Genome-wide analysis and expression profiling of Cation/H+ exchanger (CAX) family genes reveal likely functions in cadmium stress responses in poplar

Cadmium, a toxic heavy metal, seriously affects human health and ecological security. The cation/H⁺ exchanger (CAX) family is a unique metal transporter that plays a crucial role in Cd acquisition, transfer, and remission in plants. Although there are many studies related to the genome-wide analysis of Populus trichocarpa, little research has been done on the CAX family genes, especially concerning Cd stress. In this study, genome-wide analysis of the Populus CAX family identified seven stress-related CAX genes. The evolutionary tree indicated that the CaCA family genes were grouped into four clusters. Moreover, seven pairs of genes were derived by segmental duplication in poplars. Cis-acting element analysis identified numerous stress-related elements in the promoters of diverse PtrCAXs. Furthermore, some PtrCAXs were up-regulated by drought, beetle, and mechanical damage, indicating their possible function in regulating stress response. Under cadmium stress, all CAX genes in the roots were up-regulated. Our findings suggest that plants may regulate their response to Cd stress through the TF-CAXs module. Comprehensively investigating the CAX family provides a scientific basis for the phytoremediation of heavy metal pollution by Populus.

Comparative Genome Analysis of Three Komagataeibacter Strains Used for Practical Production of Nata-de-Coco

We determined the whole genome sequences of three bacterial strains, designated as FNDCR1, FNDCF1, and FNDCR2, isolated from a practical nata-de-coco producing bacterial culture. Only FNDCR1 and FNDCR2 strains had the ability to produce cellulose. The 16S rDNA sequence and phylogenetic analysis revealed that all strains belonged to the Komagataeibacter genus but belonged to a different clade within the genus. Comparative genomic analysis revealed cross-strain distribution of duplicated sequences in Komagataeibacter genomes. It is particularly interesting that FNDCR1 has many duplicated sequences within the genome independently of the phylogenetic clade, suggesting that these duplications might have been obtained specifically for this strain. Analysis of the cellulose biosynthesis operon of the three determined strain genomes indicated that several cellulose synthesis-related genes, which are present in FNDCR1 and FNDCR2, were lost in the FNDCF1 strain. These findings reveal important genetic insights into practical nata de coco-producing bacteria that can be used in food development. Furthermore, our results also shed light on the variation in their cellulose-producing abilities and illustrate why genetic traits are unstable for Komagataeibacter and Komagataeibacter-related acetic acid bacteria.

Integrated Analysis of Transcriptome and Proteome to Reveal Pupal Color Switch in Papilio xuthus Butterflies

Pupal color polyphenism in Papilio butterflies, including green, intermediate, or brown, is an excellent study system for understanding phenotypic plasticity. Previous studies suggested that development of brown pupae may be controlled by a hormone called pupal-cuticle-melanizing-hormone (PCMH) which is synthesized and secreted from brain-suboesophageal ganglion and prothoracic ganglion complexes (Br-SG-TG1) during the pre-pupa stage. However, detailed molecular mechanisms of neuroendocrine regulation in pupal color development remain unknown. In this study, we integrated the expression profiles of transcriptome and proteome at pre-pupa stages [2 h after gut purge (T1) and 3 h after forming the garter around the body (T2)] and pigmentation stages [10 h after ecdysis (T3) and 24 h after ecdysis (T4)] to identify important genes and pathways underlying the development of green and brown pupa in the swallowtail butterfly Papilio xuthus. Combined comparisons of each developmental stage and each tissue under green and brown conditions, a total of 1042 differentially expressed genes (DEGs) and 430 different abundance proteins (DAPs) were identified. Weighted gene co-expression network analysis (WGCNA) and enrichment analysis indicate that these DEGs were mainly related to oxidation-reduction, structural constituent of cuticle, and pigment binding. Soft clustering by Mfuzz and enrichment analysis indicate that these DAPs are mainly involved in tyrosine metabolism, insect hormone biosynthesis, and melanogenesis. By homologous alignment, we further identified those genes encoding neuropeptides (51), GPCRs (116), G-proteins (8), cuticular proteins (226), chitinases (16), and chitin deacetylases (8) in the whole genome of P. xuthus and analyzed their expression profiles. Although we identified no gene satisfying with hypothesized expression profile of PCMH, we found some genes in the neuropeptide cascade showed differentially expressed under two pupal color conditions. We also found that Toll signaling pathway genes, juvenile hormone (JH) related genes, and multiple cuticular proteins play important roles in the formation of selective pupal colors during the prepupal-pupal transition. Our data also suggest that both green and brown pupa include complex pigment system that is regulated by genes involved in black, blue, and yellow pigments. Our results provide important insights into the evolution of pupal protective colors among swallowtail butterflies.

Evolutionary and mechanistic diversity of Type I-F CRISPR-associated transposons

Canonical CRISPR-Cas systems utilize RNA-guided nucleases for targeted cleavage of foreign nucleic acids, whereas some nuclease-deficient CRISPR-Cas complexes have been repurposed to direct the insertion of Tn7-like transposons. Here, we established a bioinformatic and experimental pipeline to comprehensively explore the diversity of Type I-F CRISPR-associated transposons. We report DNA integration for 20 systems and identify a highly active subset that exhibits complete orthogonality in transposon DNA mobilization. We reveal the modular nature of CRISPR-associated transposons by exploring the horizontal acquisition of targeting modules and by characterizing a system that encodes both a programmable, RNA-dependent pathway, and a fixed, RNA-independent pathway. Finally, we analyzed transposon-encoded cargo genes and found the striking presence of anti-phage defense systems, suggesting a role in transmitting innate immunity between bacteria. Collectively, this study substantially advances our biological understanding of CRISPR-associated transposon function and expands the suite of RNA-guided transposases for programmable, large-scale genome engineering.

Discrimination between Carbapenem-Resistant and Carbapenem-Sensitive Klebsiella pneumoniae Strains through Computational Analysis of Surface-Enhanced Raman Spectra: a Pilot Study

In clinical settings, rapid and accurate diagnosis of antibiotic resistance is essential for the efficient treatment of bacterial infections. Conventional methods for antibiotic resistance testing are time consuming, while molecular methods such as PCR-based testing might not accurately reflect phenotypic resistance. Thus, fast and accurate methods for the analysis of bacterial antibiotic resistance are in high demand for clinical applications. In this pilot study, we isolated 7 carbapenem-sensitive Klebsiella pneumoniae (CSKP) strains and 8 carbapenem-resistant Klebsiella pneumoniae (CRKP) strains from clinical samples. Surface-enhanced Raman spectroscopy (SERS) as a label-free and noninvasive method was employed for discriminating CSKP strains from CRKP strains through computational analysis. Eight supervised machine learning algorithms were applied for sample analysis. According to the results, all supervised machine learning methods could successfully predict carbapenem sensitivity and resistance in K. pneumoniae, with a convolutional neural network (CNN) algorithm on top of all other methods. Taken together, this pilot study confirmed the application potentials of surface-enhanced Raman spectroscopy in fast and accurate discrimination of Klebsiella pneumoniae strains with different antibiotic resistance profiles. IMPORTANCE With the low-cost, label-free, and nondestructive features, Raman spectroscopy is becoming an attractive technique with great potential to discriminate bacterial infections. In this pilot study, we analyzed surfaced-enhanced Raman spectroscopy (SERS) spectra via supervised machine learning algorithms, through which we confirmed the application potentials of the SERS technique in rapid and accurate discrimination of Klebsiella pneumoniae strains with different antibiotic resistance profiles.

Genomic analysis of 600 vancomycin-resistant Enterococcus faecium reveals a high prevalence of ST80 and spread of similar vanA regions via IS1216E and plasmid transfer in diverse genetic lineages in Ireland

BACKGROUND

Vancomycin-resistant Enterococcus faecium (VREfm) cause a wide range of hospital infections. Ireland has had one of the highest invasive VREfm infection rates in Europe over the last decade, yet little is known about Irish VREfm.

OBJECTIVES

To investigate the population structure of Irish VREfm, explore diversity by analysing the vanA transposon region and compare Irish, Danish and global isolates.

METHODS

E. faecium (n = 648) from five Irish hospitals were investigated, including VREfm [547 rectal screening and 53 bloodstream infection (BSI)] isolates and 48 vancomycin-susceptible (VSEfm) BSI isolates recovered between June 2017 and December 2019. WGS and core-genome MLST (cgMLST) were used to assess population structure. Genetic environments surrounding vanA were resolved by hybrid assembly of short-read (Illumina) and long-read (Oxford Nanopore Technologies) sequences.

RESULTS

All isolates belonged to hospital-adapted clade A1 and the majority (435/648) belonged to MLST ST80. The population structure was highly polyclonal; cgMLST segregated 603/648 isolates into 51 clusters containing mixtures of screening and BSI isolates, isolates from different hospitals, and VREfm and VSEfm. Isolates within clusters were closely related (mean average ≤16 allelic differences). The majority (96.5%) of VREfm harboured highly similar vanA regions located on circular or linear plasmids with multiple IS1216E insertions, variable organization of vanA operon genes and 78.6% harboured a truncated tnpA transposase. Comparison of 648 Irish isolates with 846 global E. faecium from 30 countries using cgMLST revealed little overlap.

CONCLUSIONS

Irish VREfm are polyclonal, yet harbour a characteristic plasmid-located vanA region with multiple IS1216E insertions that may facilitate spread.

Large-Scale Phylogenetic Analysis Reveals a New Genetic Clade among Escherichia coli O26 Strains

Shiga toxin-producing Escherichia coli (STEC) O26 is the predominant non-O157 serogroup causing hemolytic uremic syndrome worldwide. Moreover, the serogroup is highly dynamic and harbors several pathogenic clones. Here, we investigated the phylogenetic relationship of STEC O26 at a global level based on 1,367 strains from 20 countries deposited in NCBI and Enterobase databases. The whole-genome-based analysis identified a new genetic clade, called ST29C4. The new clade was unique in terms of multilocus sequence type (ST29), CRISPR (group Ia), and dominant plasmid gene profile (ehxA+/katP-/espP-/etpD-). Moreover, the combination of multiple typing methods (core genome single nucleotide polymorphism [SNP] typing, CRISPR typing, and virulence genes analysis) demonstrated that this new lineage ST29C4 was in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains. Besides, we observed that ST29C4 harbored extraintestinal pathogenic E. coli (ExPEC)-related virulence gene (VG), tsh, and STEC-associated VG, stx2a, suggesting the emergence of a hybrid pathogen. The ST29C4 strains also exhibited high similarity in stx2a-prophage and integrase with the O104:H4 strain, further demonstrating its potential risk to human health. Collectively, the large-scale phylogenetic analysis extends the understanding of the clonal structure of O26 strains and provides new insights for O26 strain microevolution.

IMPORTANCE Shiga toxin-producing Escherichia coli (STEC) O26 is the second prevalent STEC serogroup only to O157, which can cause a series of diseases ranging from mild diarrhea to life-threatening hemolytic uremic syndrome (HUS). The serogroup is highly diverse and multiple clones are characterized, including ST29C1-C3 and ST21C1-C2. However, the phylogenetic relationship of these clones remains fully unclear. In this study, we revealed a new genetic clade among O26 strains, ST29C4, which was unique in terms of CRISPR, multilocus sequence type (MLST), and plasmid gene profile (PGP). Moreover, the combination of multiple typing methods demonstrated that this new clone was located in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains (i.e., ST29C1-C2 and ST21C1-C2). Overall, the large-scale phylogenetic analysis extends our current understanding of O26 microevolution.

Sentinel Surveillance Reveals Emerging Daptomycin-Resistant ST736 Enterococcus faecium and Multiple Mechanisms of Linezolid Resistance in Enterococci in the United States

Enterococcus faecalis and faecium with resistance to daptomycin and/or linezolid are emerging globally. We present the genomic characterization of daptomycin- and linezolid-resistant E. faecalis and E. faecium surveillance isolates from the United States, 2013-2016. Daptomycin resistance was low among E. faecalis (2/364, 0.5%) and E. faecium (17/344, 5%). The majority (71%, 12/17) of daptomycin-resistant E. faecium isolates belonged to the emerging ST736 clone and contained mutations in liaFSR and cls previously associated with resistance. However, 1/2 E. faecalis and 3/17 E. faecium did not contain these mutations previously associated with daptomycin resistance. Linezolid resistance was rare among E. faecalis (1/364, 0.3%) and E. faecium (2/344, 0.6%). These two E. faecium isolates, one of which was also resistant to daptomycin and vancomycin, contained the 23S rRNA nucleotide mutation (G2576T) associated with linezolid resistance. Long-read sequencing revealed the linezolid-resistant E. faecalis isolate contained chromosomal- and plasmid-encoded copies of optrA. The chromosomal optrA was located on the recently described Tn6674 multiresistance transposon. The second copy of optrA was encoded on an ∼65 kb mosaic plasmid, with component regions sharing high sequence identity to optrA-encoding multiresistance plasmids of animal origin. The optrA-encoding plasmid contained open reading frames predicted to encode proteins associated with a pheromone-responsive plasmid transfer system, and filter mating experiments confirmed the plasmid was conjugative. Continued surveillance of enterococci is necessary to assess the prevalence and trends of daptomycin and linezolid resistance in the United States, characterize resistance mechanisms and how they transfer, and monitor for emerging sequence types associated with resistance.

Identifying temporal and spatial patterns of variation from multimodal data using MEFISTO

Factor analysis is a widely used method for dimensionality reduction in genome biology, with applications from personalized health to single-cell biology. Existing factor analysis models assume independence of the observed samples, an assumption that fails in spatio-temporal profiling studies. Here we present MEFISTO, a flexible and versatile toolbox for modeling high-dimensional data when spatial or temporal dependencies between the samples are known. MEFISTO maintains the established benefits of factor analysis for multimodal data, but enables the performance of spatio-temporally informed dimensionality reduction, interpolation, and separation of smooth from non-smooth patterns of variation. Moreover, MEFISTO can integrate multiple related datasets by simultaneously identifying and aligning the underlying patterns of variation in a data-driven manner. To illustrate MEFISTO, we apply the model to different datasets with spatial or temporal resolution, including an evolutionary atlas of organ development, a longitudinal microbiome study, a single-cell multi-omics atlas of mouse gastrulation and spatially resolved transcriptomics.

The resistome of commensal Escherichia coli isolated from broiler carcasses “produced without the use of antibiotics”

Several strategies have been in place in food animal production to reduce the unnecessary use of antimicrobial agents. Beyond the monitoring of their use, the evaluation of the effect of these strategies on the occurrence and types of antimicrobial resistance (AMR) associated genes is crucial to untangle the potential emergence and spread of AMR to humans through the food chain. In the present study, the occurrence of these genes was evaluated in commensal Escherichia coli isolated from broiler carcasses “produced without the use of antibiotics” in 3 antibiotic-free (AB-free) farms in Italy in 2019. Sequenced data were analyzed along with publicly available genomes of E. coli collected in Italy from the broiler food chain from previous years (2017–2018). The genetic relationships among all 93 genomes were assessed on de novo assemblies by in silico MLST and SNP calling. Moreover, the resistomes of all genomes were investigated. According to SNP calling, genomes were gathered in three clades. Clade A encompassed, among others, ST117, ST8070, and ST1011 genomes. ST10 belonged to clade B, whereas Clade C included ST58, ST297, ST1101, and ST23 among others. Regarding the occurrence of AMR genes, a statistically significant lower occurrence of these genes in the genomes of this study in comparison to the public genomes was observed considering the whole group of genes as well as genes specifically conferring resistance to aminoglycosides, β-lactams, phenicols, trimethoprim, and lincosamides. Moreover, significant reductions were observed by comparing the whole group of AMR associated mutations, as well as those specifically for fluoroquinolones and fosfomycin resistance. Although the identification of 3° generation cephalosporin resistance associated genes in AB-free E. coli is a concern, this study provides a first indication of the impact of a more prudent use of antimicrobial agents on the occurrence of AMR genes in Italian broiler production chain. More studies are needed in next years on a higher number of genomes to confirm this preliminary observation.

Comparative Study in Zebrafish and Medaka Unravels the Mechanisms of Tissue Regeneration

Tissue regeneration has been in the spotlight of research for its fascinating nature and potential applications in human diseases. The trait of regenerative capacity occurs diversely across species and tissue contexts, while it seems to decline over evolution. Organisms with variable regenerative capacity are usually distinct in phylogeny, anatomy, and physiology. This phenomenon hinders the feasibility of studying tissue regeneration by directly comparing regenerative with non-regenerative animals, such as zebrafish (Danio rerio) and mice (Mus musculus). Medaka (Oryzias latipes) is a fish model with a complete reference genome and shares a common ancestor with zebrafish approximately 110–200 million years ago (compared to 650 million years with mice). Medaka shares similar features with zebrafish, including size, diet, organ system, gross anatomy, and living environment. However, while zebrafish regenerate almost every organ upon experimental injury, medaka shows uneven regenerative capacity. Their common and distinct biological features make them a unique platform for reciprocal analyses to understand the mechanisms of tissue regeneration. Here we summarize current knowledge about tissue regeneration in these fish models in terms of injured tissues, repairing mechanisms, available materials, and established technologies. We further highlight the concept of inter-species and inter-organ comparisons, which may reveal mechanistic insights and hint at therapeutic strategies for human diseases.

Widespread prevalence and molecular epidemiology of tet(X4) and mcr-1 harboring Escherichia coli isolated from chickens in Pakistan

The emergence and spread of plasmid-mediated tigecycline resistance gene tet(X4) and colistin resistance gene mcr-1 in Escherichia coli (E. coli) pose a potential threat to public health, due to the importance of colistin and tigecycline for treating serious clinical infections. However, the characterization of bacteria coharboring both genes was few reported. Here, we described the molecular epidemiology of tet(X4) and mcr-1 harboring E. coli strains of chicken origin in Pakistan, with methods including PCR, antimicrobial susceptibility testing, DNA transfer assays, plasmid replicon typing, whole-genome sequencing and bioinformatics analysis. The tet(X4) gene was identified in 36 isolates exhibiting high levels of tigecycline resistance (MICs, 16-128 mg/L). Worryingly, 24 of the 36 tet(X4)-bearing isolates were confirmed as colistin resistance, positive for plasmid-borne mcr-1. We observed the prevalence of tet(X4)-bearing IncFII plasmid with mcr-1-bearing IncI2 plasmid in 12 E. coli isolates, with a high co-transfer frequency except for one strain PK8233, in which tet(X4)- and mcr-1-bearing plasmids were non-transferable. Coexistence of tet(X4)-bearing IncFII plasmid with mcr-1-carrying multidrug-resistant (MDR) IncHI2 plasmid was also identified in 10 E. coli isolates, and a relatively low co-transfer frequency was obtained except PK8575, in which mcr-1 was non-transferable. The transferability of pPK8275-tetX in PK8275 and pPK8233-tetX in PK8233, that could transfer from E. coli J53 to C600 by conjugation, was interfered by certain factors in PK8275 and PK8233. This may provide new insights to prevent and control the spread of antibiotic resistance genes. Two strains were reported to co-carry tet(X4)-positive IncQ1 plasmid and mcr-1-positive IncI2 plasmid. Convergence of tet(X4) and mcr-1 genes in E. coli by conjugative or mobilizable plasmids may lead to potentially widespread transmission of such resistance genes, which may incur antibiotic-resistance crisis globally.

Eco-phylogenetic analyses reveal divergent evolution of vitamin B12 metabolism in the marine bacterial family 'Psychromonadaceae'

Cobalamin (vitamin B₁₂ ) is an essential micronutrient required by both prokaryotes and eukaryotes. Nevertheless, with high genetic and metabolic cost, de novo cobalamin biosynthesis is exclusive to a subset of prokaryotic taxa. Many Cyanobacterial and Archaeal taxa have been implicated in de novo cobalamin biosynthesis in epi- and mesopelagic ocean respectively. However, the contributions of Gammaproteobacteria particularly the family 'Psychromonadaceae' is largely unknown. Through phylo-pangenomic analyses using concatenated single-copy proteins and homologous gene clusters respectively, the phylogenies within 'Psychromonadaceae' recapitulate both their taxonomic delineations and environmental distributions. Moreover, uneven distribution of cobalamin de novo biosynthetic operon and cobalamin-dependent light-responsive regulon were observed, and of which the linkages to the environmental conditions where cobalamin availability and light regime can be varied respectively were discussed, suggesting the impacts of ecological divergence in shaping their disparate cobalamin-related metabolisms. Functional analysis demonstrated a varying degree of cobalamin dependency for both central metabolic processes and cobalamin-mediated light-responsive regulation, and underlying sequence characteristics of cis- and trans-regulatory elements were revealed. Our findings emphasized the potential roles of cobalamin in shaping the ecological distributions and driving the metabolic evolution in the marine bacterial family 'Psychromonadaceae', and have further implications for an improved understanding of nutritional interdependencies and community metabolism modulated by cobalamin.

A phosphorus-limitation induced, functionally conserved DUF506 protein is a repressor of root hair elongation in plants

Root hairs (RHs) function in nutrient and water acquisition, root metabolite exudation, soil anchorage and plant-microbe interactions. Longer or more abundant RHs are potential breeding traits for developing crops that are more resource-use efficient and can improve soil health. While many genes are known to promote RH elongation, relatively little is known about genes and mechanisms that constrain RH growth. Here we demonstrate that a DOMAIN OF UNKNOWN FUNCTION 506 (DUF506) protein, AT3G25240, negatively regulates Arabidopsis thaliana RH growth. The AT3G25240 gene is strongly and specifically induced during phosphorus (P)-limitation. Mutants of this gene, which we call REPRESSOR OF EXCESSIVE ROOT HAIR ELONGATION 1 (RXR1), have much longer RHs, higher phosphate content and seedling biomass, while overexpression of the gene exhibits opposite phenotypes. Co-immunoprecipitation, pull-down and bimolecular fluorescence complementation (BiFC) analyses reveal that RXR1 physically interacts with a RabD2c GTPase in nucleus, and a rabd2c mutant phenocopies the rxr1 mutant. Furthermore, N-terminal variable region of RXR1 is crucial for inhibiting RH growth. Overexpression of a Brachypodium distachyon RXR1 homolog results in repression of RH elongation in Brachypodium. Taken together, our results reveal a novel DUF506-GTPase module with a prominent role in repression of plant RH elongation especially under P stress.

Exploring the role of framework mutations in enabling breadth of a cross-reactive antibody (CR3022) against the SARS-CoV-2 RBD and its variants of concern

Cross-reactive and broadly neutralizing antibodies against surface proteins of diverse strains of rapidly evolving viral pathogens like SARS-CoV-2 can prevent infection and therefore are crucial for the development of effective universal vaccines. While antibodies typically incorporate mutations in their complementarity determining regions during affinity maturation, mutations in the framework regions have been reported as players in determining properties of broadly neutralizing antibodies against HIV and the Influenza virus. We propose an increase in the cross-reactive potential of CR3022 against the emerging SARS- CoV-2 variants of concern through enhanced conformational flexibility. In this study, we use molecular dynamics simulations, in silico mutagenesis, structural modeling, and docking to explore the role of light chain FWR mutations in CR3022, a SARS-CoV anti-spike (S)-protein antibody cross-reactive to the S-protein receptor binding domain of SARS-CoV-2. Our study shows that single substitutions in the light chain framework region of CR3022 with conserved epitopes across SARS-CoV strains allow targeting of diverse antibody epitope footprints that align with the epitopes of recently-categorized neutralizing antibody classes while enabling binding to more than one strain of SARS-CoV-2. Our study has implications for rapid and evolution-based engineering of broadly neutralizing antibodies and reaffirms the role of framework mutations in effective change of antibody orientation and conformation via improved flexibility.Communicated by Ramaswamy H. Sarma.

Negative Modulation of the Angiogenic Cascade Induced by Allosteric Kinesin Eg5 Inhibitors in a Gastric Adenocarcinoma In Vitro Model

Eg5 is a kinesin essential in bipolar spindle formation, overexpressed in tumours, thus representing a new target in cancer therapy. We aimed at evaluating the anti-cancer activity of Eg5 thiadiazoline inhibitors 2 and 41 on gastric adenocarcinoma cells (AGS), focusing on the modulation of angiogenic signalling. Docking studies confirmed a similar interaction with Eg5 to that of the parent compound K858. Thiadiazolines were also tested in combination with Hesperidin (HSD). Cell cycle analysis reveals a reduction of G1 and S phase percentages when 41 is administered as well as HSD in combination with K858. Western blot reveals Eg5 inhibitors capability to reduce PI3K, p-AKT/Akt and p-Erk/Erk expressions; p-Akt/Akt ratio is even more decreased in HSD+2 sample than the p-Erk/Erk ratio in HSD+41 or K858. VEGF expression is reduced when HSD+2 and HSD+41 are administered with respect to compounds alone, after 72 h. ANGPT2 gene expression increases in cells treated with 41 and HSD+2 compared to K858. The wound-healing assay highlights a reduction in the cut in HSD+2 sample compared to 2 and HSD. Thus, Eg5 inhibitors appear to modulate angiogenic signalling by controlling VEGF activity even better if combined with HSD. Overall, Eg5 inhibitors can represent a promising starting point to develop innovative anti-cancer strategies.

Genomic Changes within a Subset of IncI2 Plasmids Associated with Dissemination of mcr-1 Genes and Other Important Antimicrobial Resistance Determinants

IncI2 plasmids appear to have only recently become associated with resistance genes; however, their tendency to carry resistance to the antibiotics of last resort and their widespread distribution increase their relative importance. In this study, we describe lineages within this plasmid family that have an increased likelihood of acquisition of antimicrobial resistance genes. Globally distributed mcr-1-carrying IncI2 plasmids were found to cluster with other IncI2 plasmids carrying extended-spectrum beta-lactamase genes, and separately from the non-resistant IncI2 plasmids. In addition, insertion sequence (IS) elements with no direct association with the acquired resistance genes also clustered with the resistance plasmids in the phylogenetic tree. In recognition of the biased sequencing of resistant plasmids globally, the analysis was also performed on resistant and non-resistant IncI2 plasmids sequenced in the USA through government surveillance efforts that do not rely on antibiotic selection. This analysis confirmed a distinct clustering associated with both resistance and mobile elements and identified possible genomic changes in core genes that correlate with increased acquisition of foreign DNA. This work highlights a potential genetic mechanism for increased uptake of foreign DNA within this prevalent family of plasmids.

Epidemiology and Genetic Analysis of SARS-CoV-2 in Myanmar during the Community Outbreaks in 2020

We aimed to analyze the situation of the first two epidemic waves in Myanmar using the publicly available daily situation of COVID-19 and whole-genome sequencing data of SARS-CoV-2. From March 23 to December 31, 2020, there were 33,917 confirmed cases and 741 deaths in Myanmar (case fatality rate of 2.18%). The first wave in Myanmar from March to July was linked to overseas travel, and then a second wave started from Rakhine State, a western border state, leading to the second wave spreading countrywide in Myanmar from August to December 2020. The estimated effective reproductive number (R_t) nationwide reached 6–8 at the beginning of each wave and gradually decreased as the epidemic spread to the community. The whole-genome analysis of 10 Myanmar SARS-CoV-2 strains together with 31 previously registered strains showed that the first wave was caused by GISAID clade O or PANGOLIN lineage B.6 and the second wave was changed to clade GH or lineage B.1.36.16 with a close genetic relationship with other South Asian strains. Constant monitoring of epidemiological situations combined with SARS-CoV-2 genome analysis is important for adjusting public health measures to mitigate the community transmissions of COVID-19.

Population structure analysis and laboratory monitoring of Shigella by core-genome multilocus sequence typing

The laboratory surveillance of bacillary dysentery is based on a standardised Shigella typing scheme that classifies Shigella strains into four serogroups and more than 50 serotypes on the basis of biochemical tests and lipopolysaccharide O-antigen serotyping. Real-time genomic surveillance of Shigella infections has been implemented in several countries, but without the use of a standardised typing scheme. Here, we study over 4000 reference strains and clinical isolates of Shigella, covering all serotypes, with both the current serotyping scheme and the standardised EnteroBase core-genome multilocus sequence typing scheme (cgMLST). The Shigella genomes are grouped into eight phylogenetically distinct clusters, within the E. coli species. The cgMLST hierarchical clustering (HC) analysis at different levels of resolution (HC2000 to HC400) recognises the natural population structure of Shigella. By contrast, the serotyping scheme is affected by horizontal gene transfer, leading to a conflation of genetically unrelated Shigella strains and a separation of genetically related strains. The use of this cgMLST scheme will facilitate the transition from traditional phenotypic typing to routine whole-genome sequencing for the laboratory surveillance of Shigella infections.

Lab-based surveillance of Shigella has traditionally been based on serotyping but increasing availability of whole genome sequencing could enable higher resolution typing. Here, the authors apply a core genome multilocus sequence typing scheme to Shigella sequence data and describe its population structure.

Microbial Consortia Involved in Traditional Sicilian Sourdough: Characterization of Lactic Acid Bacteria and Yeast Populations

Sourdough is one of the oldest starters traditionally used for making baked goods, offering several advantages to the sensory, rheology, and shelf life of final products. The present study investigated, for the first time, the microbiota of spontaneously fermented Maiorca dough samples collected from bakeries located in Sicily (Italy). Four sourdough samples (M1, M2, M3, and M4), were produced using Triticum vulgare Host. var. albidum Koern (Maiorca grain) were subjected to LAB and yeasts isolation and identification at the species level. The in-depth characterization of the lactobacilli population revealed that Lactiplantibacillus plantarum and Levilactobacillus brevis unquestionably dominated the Maiorca sourdough ecosystem. Concerning the yeasts community, high species diversity was found. Saccharomyces cerevisiae and Wickerhamomyces anomalus were the most frequently isolated species. In addition, Torulaspora delbrueckii, Pichia kluyveri, Candida boidinii, and Candida diddensiae were also detected. Investigations on both pro-technological and functional traits of the isolated strains could lead to the selection of starters for the production of baked goods.

Genomic Characterization of Carbapenem-Non-susceptible Pseudomonas aeruginosa Clinical Isolates From Saudi Arabia Revealed a Global Dissemination of GES-5-Producing ST235 and VIM-2-Producing ST233 Sub-Lineages

Carbapenem-resistant P. aeruginosa has become a major clinical problem due to limited treatment options. However, studies assessing the trends in the molecular epidemiology and mechanisms of antibiotic resistance in this pathogen are lacking in Saudi Arabia. Here, we reported the genome characterization in a global context of carbapenem non-susceptible clinical isolates from a nationally representative survey. The antibiotic resistance profiles of the isolates (n = 635) collected over 14 months between March 2018 and April 2019 from different geographical regions of Saudi Arabia showed resistance rates to relevant β-lactams, aminoglycosides and quinolones ranging between 6.93 and 27.56%. Overall, 22.52% (143/635) of the isolates exhibited resistance to both imipenem and meropenem that were mainly explained by porin loss and efflux overexpression. However, 18.18% of resistant isolates harbored genes encoding GES (69.23%), VIM (23.07%), NDM (3.85%) or OXA-48-like (3.85%) carbapenemases. Most common GES-positive isolates produced GESs −5, −15 or −1 and all belonged to ST235 whereas the VIM-positive isolates produced mainly VIM-2 and belonged to ST233 or ST257. GES and VIM producers were detected at different sampling periods and in different surveyed regions. Interestingly, a genome-wide comparison revealed that the GES-positive ST235 and VIM-2-positive ST233 genomes sequenced in this study and those available through public databases from various locations worldwide, constituted each a phylogenetically closely related sub-lineage. Profiles of virulence determinants, antimicrobial resistance genes and associated mobile elements confirmed relatedness within each of these two different sub-lineages. Sequence analysis located the bla_GES gene in nearly all studied genomes (95.4%) in the same integrative conjugative element that also harbored the acc(6′)-Ib, aph(3′)-XV, aadA6, sul1, tet(G), and catB resistance genes while bla_VIM–2 in most (98.89%) ST233-positive genomes was co-located with aac(6′)-I1, dfrB-5, and aac(3′)-Id in the same class I integron. The study findings revealed the global spread of GES-5 ST235 and VIM-2 ST233 sub-lineages and highlighted the importance of routine detection of rare β-lactamases.

Detecting Flavobacterial Fish Pathogens in the Environment via High-Throughput Community Analysis

Diseases caused by the fish pathogens Flavobacterium columnare and Flavobacterium psychrophilum are major contributors of preventable losses in the aquaculture industry. The persistent and difficult-to-control infections caused by these bacteria make timely intervention and prophylactic elimination of pathogen reservoirs important measures to combat these disease-causing agents. In this study, we present two independent assays for detecting these pathogens in a range of environmental samples. Natural water samples were inoculated with F. columnare and F. psychrophilum over 5 orders of magnitude, and pathogen levels were detected using Illumina MiSeq sequencing and droplet digital PCR. Both detection methods accurately identified pathogen-positive samples and showed good agreement in quantifying each pathogen. Additionally, the real-world application of these approaches was demonstrated using environmental samples collected at a rainbow trout (Oncorhynchus mykiss) aquaculture facility. These results show that both methods can serve as useful tools for surveillance efforts in aquaculture facilities, where the early detection of these flavobacterial pathogens may direct preventative measures to reduce disease occurrence. IMPORTANCE Early detection of a deadly disease outbreak in a population can be the difference between mass mortality or mitigated effects. In the present study, we evaluated and compared two molecular techniques for detecting economically impactful aquaculture pathogens. We demonstrate that one of these techniques, 16S rRNA gene sequencing using Illumina MiSeq technology, provides the ability to accurately detect two freshwater fish pathogens, F. columnare and F. psychrophilum, while simultaneously profiling the native microbial community. The second technique, droplet digital PCR, is commonly used for pathogen detection, and the results obtained using the assays we designed with this method served to validate those obtained using the MiSeq method. These two methods offer distinct advantages. The MiSeq method pairs pathogen detection and microbial community profiling to answer immediate and long-term fish health concerns, while the droplet digital PCR method provides fast and highly sensitive detection that is useful for surveillance and rapid clinical responses.

Synthesis of Nuclear and Chloroplast Data Combined With Network Analyses Supports the Polyploid Origin of the Apple Tribe and the Hybrid Origin of the Maleae-Gillenieae Clade

Plant biologists have debated the evolutionary origin of the apple tribe (Maleae; Rosaceae) for over a century. The "wide-hybridization hypothesis" posits that the pome-bearing members of Maleae (base chromosome number x = 17) resulted from a hybridization and/or allopolyploid event between progenitors of other tribes in the subfamily Amygdaloideae with x = 8 and x = 9, respectively. An alternative "spiraeoid hypothesis" proposed that the x = 17 of Maleae arose via the genome doubling of x = 9 ancestors to x = 18, and subsequent aneuploidy resulting in x = 17. We use publicly available genomic data-448 nuclear genes and complete plastomes-from 27 species representing all major tribes within the Amygdaloideae to investigate evolutionary relationships within the subfamily containing the apple tribe. Specifically, we use network analyses and multi-labeled trees to test the competing wide-hybridization and spiraeoid hypotheses. Hybridization occurred between an ancestor of the tribe Spiraeeae (x = 9) and an ancestor of the clade Sorbarieae (x = 9) + Exochordeae (x = 8) + Kerrieae (x = 9), giving rise to the clade Gillenieae (x = 9) + Maleae (x = 17). The ancestor of the Maleae + Gillenieae arose via hybridization between distantly related tribes in the Amygdaloideae (i.e., supporting the wide hybridization hypothesis). However, some evidence supports an aspect of the spiraeoid hypothesis-the ancestors involved in the hybridization event were likely both x = 9, so genome doubling was followed by aneuploidy to result in x = 17 observed in Maleae. By synthesizing existing genomic data with novel analyses, we resolve the nearly century-old mystery regarding the origin of the apple tribe. Our results also indicate that nuclear gene tree-species tree conflict and/or cytonuclear conflict are pervasive at several other nodes in subfamily Amygdaloideae of Rosaceae.

Comparative Transcriptome Analysis in Taraxacum koksaghyz to Identify Genes that Determine Root Volume and Root Length

The Russian dandelion (Taraxacum koksaghyz, family Asteraceae) produces large amounts of natural rubber in the laticifers of its roots. This species has been proposed as an alternative source of natural rubber to augment or partly replace the rubber tree (Hevea brasiliensis) but domestication would require genetic improvement to increase rubber yields and agronomic optimization to facilitate harvesting and processing. Optimization has focused thus far on the size and shape of the roots, the primary storage organ for natural rubber and inulin. However, the corresponding genetic factors are poorly understood. Here we describe the comparative transcriptomic analysis of root tissues from T. koksaghyz plant sets featuring different root sizes and shapes, aiming to identify differentially expressed genes correlating with root length or root diameter in the upper root and root tip. The resulting datasets revealed multiple candidate genes for each trait and root part, including a glucan endo-1,3-β-d-glucosidase, an allene oxide synthase 3, and a TIFY10A/JAZ1 homolog. These three genes were tested by qRT-PCR in outdoor-grown plants with diverse root morphology, and the expression of two genes correlated with the appropriate root morphotype, confirming the effectiveness of our method. We evaluated the candidate genes to gain insight into their potential functions in root development. Such candidate genes could be suitable for marker-assisted breeding programs in the future.

ITS DNA Barcoding Reveals That Halophila stipulacea Still Remains the Only Non-Indigenous Seagrass of the Mediterranean Sea

Non-indigenous species (NIS) are one of the major threats to the native marine ecosystems of the Mediterranean Sea. Halophila stipulacea was the only exotic seagrass of the Mediterranean until 2018, when small patches of a species morphologically identified as Halophila decipiens were reported in Salamina Island, Greece. Given the absence of reproductive structures during the identification and the taxonomic ambiguities known to lead to misidentifications on this genus, we reassessed the identity of this new exotic record using DNA barcoding (rbcL, matK and ITS) and the recently published taxonomic key. Despite their morphologic similarity to H. decipiens based on the new taxonomic key, the specimens showed no nucleotide differences with H. stipulacea specimens (Crete) for the three barcodes and clustered together on the ITS phylogenetic tree. Considering the high species resolution of the ITS region and the common morphological variability within the genus, the unequivocal genetic result suggests that the Halophila population found in Salamina Island most likely corresponds to a morphologically variant H. stipulacea. Our results highlight the importance of applying an integrated taxonomic approach (morphological and molecular) to taxonomically complex genera such as Halophila, in order to avoid overlooking or misreporting species range shifts, which is essential for monitoring NIS introductions.

Redefining the coenzyme A transferase superfamily with a large set of manually-annotated proteins

The coenzyme A (CoA) transferases are a superfamily of proteins central to the metabolism of acetyl-CoA and other CoA thioesters. They are diverse group, catalyzing over a hundred biochemical reactions and spanning all three domains of life. A deeply rooted idea, proposed two decades ago, is these enzymes fall into three families (I, II, III). Here we find they fall into different families, which we achieve by analyzing all CoA transferases characterized to date. We manually annotated 94 CoA transferases with functional information (including rates of catalysis for 208 reactions) from 97 publications. This represents all enzymes we could find in the primary literature, and it is double the number annotated in four protein databases (BRENDA, KEGG, MetaCyc, UniProt). We found family I transferases are not closely related to each other in terms of sequence, structure, and reactions catalyzed. This family is not even monophyletic. These problems are solved by regrouping the three families into six, including one family with many non-CoA transferases. The problem (and solution) became apparent only by analyzing our large set of manually-annotated proteins. It would have been missed if we had used the small number of proteins annotated in UniProt and other databases. Our work is important to understanding the biology of CoA transferases. It also warns investigators doing phylogenetic analyses of proteins to go beyond information in databases. This article is protected by copyright. All rights reserved.

Occurrence of Rotavirus A Genotypes and Other Enteric Pathogens in Diarrheic Suckling Piglets from Spanish Swine Farms

Simple Summary

Neonatal diarrhea is a major cause of economic losses in the swine industry worldwide and has significant impact in Spain, which is one of the biggest pork producers globally. Multiple infectious agents can contribute to this condition, with some viruses such as species A rotavirus (RVA) playing a major role. Studies on their occurrence and genetic diversity are essential for development of RVA vaccines. In this study, fecal samples from diarrheic suckling piglets originating from farms distributed throughout Spain were analyzed for RVA and four other common enteric pathogens using molecular methods. The individual prevalence was 89.4%, 64.4%, 44.9%, 33.7% and 4.4% for Clostridiumperfringens, Clostridioides (formerly Clostridium) difficile, species A rotavirus, species C rotavirus and porcine epidemic diarrhea virus, respectively. Most specimens (96.9%) were positive for at least one of the target pathogens and concurrent infections were common. The molecular characterization of RVA positive specimens of specific genes used for genotyping revealed the extensive genetic diversity of RVA strains circulating in swine herds in Spain. Comparison with genotypes contained in the commercial vaccine available in Spain showed differences in the identity of the predominant RVA genotypes from diarrheic piglets in the sampled pig farms. These findings contribute to the surveillance of RVA strains circulating in swine herds in Spain and may help optimize target vaccine design.

Abstract

Species A rotavirus (RVA) is a major viral pathogen causing diarrhea in suckling piglets. Studies on its genetic heterogeneity have implications for vaccine efficacy in the field. In this study, fecal samples (n = 866) from diarrheic piglets younger than 28 days were analyzed over a two-year period (2018–2019). Samples were submitted from 426 farms located in 36 provinces throughout Spain and were tested using real-time PCR (qPCR) and reverse transcription real-time PCR (RT-qPCR) for five enteric pathogens. The individual prevalence was 89.4%, 64.4%, 44.9%, 33.7% and 4.4% for Clostridiumperfringens, Clostridioides (formerly Clostridium) difficile, species A rotavirus, species C rotavirus and porcine epidemic diarrhea virus, respectively. Most specimens (96.9%) were positive for at least one of the target pathogens, and more than 80% of samples harbored mixed infections. Nucleotide sequencing of 70 specimens positive for RVA revealed the presence of the VP7 genotypes G4, G9, G3, G5, G11 and the VP4 genotypes P7, P23, P6 and P13, with the combinations G4P7 and G9P23 being the most prevalent, and especially in the areas with the highest pig population. The study shows the extensive genetic diversity of RVA strains as well as discrepancies with the genotypes contained in the vaccine available in Spain, and multiple amino acid differences in antigenic epitopes of different G- and P- genotypes with the vaccine strains. Further investigations are needed to determine the efficacy of the vaccine to confer clinical protection against heterologous strains.

Comparative Genome Analysis Reveals Accumulation of Single-Nucleotide Repeats in Pathogenic Escherichia Lineages

Homopolymeric tracts (HPTs) can lead to phase variation and DNA replication slippage, driving adaptation to environmental changes and evolution of genes and genomes. However, there is limited information on HPTs in Escherichia; therefore, we conducted a comprehensive cross-strain search for HPTs in Escherichia genomes. We determined the HPT genomic distribution and identified a pattern of high-frequency HPT localization in pathogenic Escherichia lineages. Notably, HPTs localized near transcriptional regulatory genes. Additionally, excessive repeats accumulated in toxin-coding genes. Moreover, the genomic localization of some HPTs might be derived from exogenous DNA, such as that of bacteriophages. Altogether, our findings may prove useful for understanding the role of HPTs in Escherichia genomes.

Whole-Genome Sequence of a Porcine Circovirus Type 2 Strain Detected in Assam, India

Porcine circovirus-associated disease caused by porcine circovirus type 2 (PCV2) is a vital threat to the global pig industry. In this study, we have characterized the complete genome sequence of a PCV2 isolate, namely, Assam-01, belonging to the genotype PCV2d.

Discovery of Two Novel Oxidases Using a High-Throughput Activity Screen

Discovery of novel enzymes is a challenging task, yet a crucial one, due to their increasing relevance as chemical catalysts and biotechnological tools. In our work we present a high-throughput screening approach to discovering novel activities. A screen of 96 putative oxidases with 23 substrates led to the discovery of two new enzymes. The first enzyme, N-acetyl-D-hexosamine oxidase (EC 1.1.3.29) from Ralstonia solanacearum, is a vanillyl alcohol oxidase-like flavoprotein displaying the highest activity with N-acetylglucosamine and N-acetylgalactosamine. Before our discovery of the enzyme, its activity was an orphan one - experimentally characterized but lacking the link to amino acid sequence. The second enzyme, from an uncultured marine euryarchaeota, is a long-chain alcohol oxidase (LCAO, EC 1.1.3.20) active with a range of fatty alcohols, with 1-dodecanol being the preferred substrate. The enzyme displays no sequence similarity to previously characterised LCAOs, and thus is a completely novel representative of a protein with such activity.

Identification of methicillin-resistant Staphylococcus aureus ST8 isolates in China with potential high virulence

Methicillin-resistant Staphylococcus aureus (MRSA) ST8 strains have spread worldwide, causing outbreaks in various regions. However, this clone has only been sporadically reported in China. Consequently, detailed information regarding the phylogeny and potential virulence of S. aureus ST8 strains in China remains unknown. In this study, we characterized six ST8 strains collected from three tertiary hospitals in China, including three MRSA (MR50, MR526, and MR254) and three MSSA (H78, H849 and H863). Whole genome sequencing and phylogenetic analysis showed that the six strains formed two separate clusters, including two (MR50 and MR526) and four (MR254, H78, H849 and H863) isolates, respectively. Among them, MR50 and MR526 harboured spa t008, SCCmec IVa, arginine catabolic mobile element, and were phylogenetically close to the epidemic USA300 strains, while other four strains belonged to spa t9101 and formed a unique branch. MR254 carried a novel hybrid SCCmec element (namely SCCmec254). Same as the USA300 prototype strain LAC, the China S. aureus ST8 strains produced weak biofilms except MR254. Among them, MR254 had significantly stronger haemolysis ability and higher α-toxin levels than others, while MR526 showed comparable haemolysis and α-toxin production levels as USA300-LAC. In mouse skin abscess model, MR254 showed particularly strong invasions, accompanied by necrosis, while MR526 exhibited similar infection levels as USA300-LAC. These data suggested that the China MRSA ST8 isolates (e.g. MR254 and MR526) were highly virulent, displaying higher or similar virulence potential as the epidemic USA300 strain. Active surveillance should be enacted to closely monitor the further spread of these hyper-virulent MRSA strains in China.