Chromosome-level genome assembly and demographic history of Euryodendron excelsum in monotypic genus endemic to China

Euryodendron excelsum is in a monotypic genus Euryodendron, endemic to China. It has intermediate morphisms in the Pentaphylacaceae or Theaceae families, which make it distinct. Due to anthropogenic disturbance, E. excelsum is currently found in very restricted and fragmented areas with extremely small populations. Although much research and effort has been applied towards its conservation, its long-term survival mechanisms and evolutionary history remain elusive, especially from a genomic aspect. Therefore, using a combination of long/short whole genome sequencing, RNA sequencing reads, and Hi-C data, we assembled and annotated a high-quality genome for E. excelsum. The genome assembly of E. excelsum comprised 1,059,895,887 bp with 99.66% anchored into 23 pseudo-chromosomes and a 99.0% BUSCO completeness. Comparative genomic analysis revealed the expansion of terpenoid and flavonoid secondary metabolite genes, and displayed a tandem and/or proximal duplication framework of these genes. E. excelsum also displayed genes associated with growth, development, and defence adaptation from whole genome duplication. Demographic analysis indicated that its fluctuations in population size and its recent population decline were related to cold climate changes. The E. excelsum genome assembly provides a highly valuable resource for evolutionary and ecological research in the future, aiding its conservation, management, and restoration.

Genome-based assessment of safety characteristics of Lacticaseibacillus paracasei NY1301 and genomic differences in closely related strains marketed as probiotics

The probiotic attributes of Lacticaseibacillus paracasei NY1301 were comprehensively characterized, and a comparison between the closely related LcA (Actimel) and LcY (Yakult) probiotic strains was conducted using genomic tools. All strains exhibited high genetic similarity and likely shared a common ancestor; differences were primarily expressed as minor chromosomal re-arrangements, substitutions, insertions, and deletions. Compared with LcY, NY1301 exhibited 125 single-nucleotide polymorphisms. NY1301 lacked virulence factors, antibiotic resistance genes, and mutations associated with antibiotic resistance and had a 46-kbp prophage. This prophage is spontaneously induced at low levels and remains in a non-lytic state under standard culture conditions. The observed causal adaptive mutations were likely related to niche adaptation within the respective laboratory or manufacturing processes that occurred during the maintenance of the strains. However, the phenotypic effects of these genomic differences remain unclear. To validate the safety of NY1301, we conducted an open-label trial with healthy participants who consumed excessive amounts of NY1301 (3.0 × 1011 cfu) daily for 28 days. The results of this trial and those of other in vivo studies, coupled with the long history of human consumption without established risks to humans, provide strong evidence confirming the safety of NY1301.

Mechanism and evolutionary analysis of Yarrowia lipolytica CA20 capable of producing erythritol with a high yield based on comparative genomics

Combined mutagenesis is widely applied for the breeding of robust Yarrowia lipolytica used in the production of erythritol. However, the changes of genome after mutagenesis remains unclear. This study aimed to unravel the mechanism involved in the improved erythritol synthesis of CA20 and the evolutionary relationship between different Y. lipolytica by comparative genomics analysis. The results showed that the genome size of Y. lipolytica CA20 was 20,420,510 bp, with a GC content of 48.97%. There were 6330 CDS and 649 ncRNA (non-coding RNA) in CA20 genome. Average nucleotide identity (ANI) analysis showed that CA20 genome possessed high similarity (ANI > 99.50%) with other Y. lipolytica strains, while phylogenetic analysis displayed that CA20 was classified together with Y. lipolytica IBT 446 and Y. lipolytica H222. CA20 shared 5342 core orthologous genes with the 8 strains while harbored 65 specific genes that mainly participated in the substrate and protein transport processes. CA20 contained 166 genes coding for carbohydrate-active enzymes (CAZymes), which was more than that found in other strains (108－137). Notably, 4, 2, and 13 different enzymes belonging to glycoside hydrolases (GHs), glycosyltransferases (GTs), and carbohydrate esterases (CEs), respectively, were only found in CA20. The enzymes involved in the metabolic pathway of erythritol were highly conserved in Y. lipolytica, except for transaldolase (TAL1). In addition, the titer and productivity of erythritol by CA20 were 190.97 g/L and 1.33 g/L/h, respectively, which were significantly higher than that of WT5 wherein 128.61 g/L and 0.92 g/L/h were obtained (P< 0.001). Five frameshift mutation genes and 15 genes harboring nonsynonymous mutation were found in CA20 compared with that of WT5. Most of these genes were involved in the cell division, cell wall synthesis, protein synthesis, and protein homeostasis maintenance. These findings suggested that the genome of Y. lipolytica is conserved during evolution, and the variance of living environment is one important factor leading to genome divergence. The varied number of CAZymes existed in Y. lipolytica is one factor that contributes to the performance difference. The increased synthesis of erythritol by Y. lipolytica CA20 is correlated with the improvement of the stability of cell structure and internal environment. The results of this study provide a basis for the directional breeding of robust strains used in erythritol production.

Comparative Genomic Analysis and Metabolic Potential Profiling of a Novel Culinary-Medicinal Mushroom, Hericium rajendrae (Basidiomycota)

Hericium rajendrae is an emerging species in the genus Hericium with few members. Despite being highly regarded due to its rarity, knowledge about H. rajendrae remains limited. In this study, we sequenced, de novo assembled, and annotated the complete genome of H. rajendrae NPCB A08, isolated from the Qinling Mountains in Shaanxi, China, using the Illumina NovaSeq and Nanopore PromethION technologies. Comparative genomic analysis revealed similarities and differences among the genomes of H. rajendrae, H. erinaceus, and H. coralloides. Phylogenomic analysis revealed the divergence time of the Hericium genus, while transposon analysis revealed evolutionary characteristics of the genus. Gene family variation reflected the expansion and contraction of orthologous genes among Hericium species. Based on genomic bioinformation, we identified the candidate genes associated with the mating system, carbohydrate-active enzymes, and secondary metabolite biosynthesis. Furthermore, metabolite profiling and comparative gene clusters analysis provided strong evidence for the biosynthetic pathway of erinacines in H. rajendrae. This work provides the genome of H. rajendrae for the first time, and enriches the genomic content of the genus Hericium. These findings also facilitate the application of H. rajendrae in complementary drug research and functional food manufacturing, advancing the field of pharmaceutical and functional food production involving H. rajendrae.

A Comprehensive Analysis of Chloroplast Genome Provides New Insights into the Evolution of the Genus Chrysosplenium

Chrysosplenium, a perennial herb in the family Saxifragaceae, prefers to grow in low light and moist environments and is divided into two sections of Alternifolia and Oppositifolia based on phyllotaxy. Although there has been some progress in the phylogeny of Chrysosplenium over the years, the phylogenetic position of some species is still controversial. In this study, we assembled chloroplast genomes (cp genomes) of 34 Chrysosplenium species and performed comparative genomic and phylogenetic analyses in combination with other cp genomes of previously known Chrysosplenium species, for a total of 44 Chrysosplenium species. The comparative analyses revealed that cp genomes of Chrysosplenium species were more conserved in terms of genome structure, gene content and arrangement, SSRs, and codon preference, but differ in genome size and SC/IR boundaries. Phylogenetic analysis showed that cp genomes effectively improved the phylogenetic support and resolution of Chrysosplenium species and strongly supported Chrysosplenium species as a monophyletic taxon and divided into three branches. The results also showed that the sections of Alternifolia and Oppositifolia were not monophyletic with each other, and that C. microspermum was not clustered with other Chrysosplenium species with alternate leaves, but with C. sedakowii into separate branches. In addition, we identified 10 mutational hotspot regions that could serve as potential DNA barcodes for Chrysosplenium species identification. In contrast to Peltoboykinia, the clpP and ycf2 genes of Chrysosplenium were subjected to positive selection and had multiple significant positive selection sites. We further detected a significant positive selection site on the petG gene between the two sections of Chrysosplenium. These evolutionary characteristics may be related to the growth environment of Chrysosplenium species. This study enriches the cp genomes of Chrysosplenium species and provides a reference for future studies on its evolution and origin.

A Chromosome-level Genome Assembly and Evolution Analysis of Andrena camellia (Hymenoptera: Andrenidae)

Andrena camellia, an effective pollinator of the economicallysignificant crop Camellia oleifera, can withstand the toxic pollen of C. oleifera, making A. camellia a crucial for resource conservation and cultivation of C. oleifera. In this study, the whole genome of A. camellia was sequenced on the Oxford Nanopore platform. The assembled genome size was 340.73 Mb including 50 scaffolds (N50=47.435 Mb) and 131 contigs (N50=17.2 Mb). A total of 11, 258 protein-coding genes were annotated, in addition, 1,104 non-coding RNAs were identified. Further analysis that some chromosomes of A. camellia have a high level of synteny with those of Apis mellifera, Osmia bicornis and Andrena minutula. Thus, our reported genome of A. camellia serves as a valuable resource for studying species evolution, behavioral biology, and adaption to toxic pollen of C. oleifera.

2022 updates to the Rat Genome Database: a Findable, Accessible, Interoperable, and Reusable (FAIR) resource

The Rat Genome Database (RGD, https://rgd.mcw.edu) has evolved from simply a resource for rat genetic markers, maps, and genes, by adding multiple genomic data types and extensive disease and phenotype annotations and developing tools to effectively mine, analyze, and visualize the available data, to empower investigators in their hypothesis-driven research. Leveraging its robust and flexible infrastructure, RGD has added data for human and eight other model organisms (mouse, 13-lined ground squirrel, chinchilla, naked mole-rat, dog, pig, African green monkey/vervet, and bonobo) besides rat to enhance its translational aspect. This article presents an overview of the database with the most recent additions to RGD's genome, variant, and quantitative phenotype data. We also briefly introduce Virtual Comparative Map (VCMap), an updated tool that explores synteny between species as an improvement to RGD's suite of tools, followed by a discussion regarding the refinements to the existing PhenoMiner tool that assists researchers in finding and comparing quantitative data across rat strains. Collectively, RGD focuses on providing a continuously improving, consistent, and high-quality data resource for researchers while advancing data reproducibility and fulfilling Findable, Accessible, Interoperable, and Reusable (FAIR) data principles.

The First Whole Genome Sequencing of Agaricus bitorquis and Its Metabolite Profiling

Agaricus bitorquis, an emerging wild mushroom with remarkable biological activities and a distinctive oversized mushroom shape, has gained increasing attention in recent years. Despite its status as an important resource of wild edible fungi, knowledge about this mushroom is still limited. In this study, we used the Illumina NovaSeq and Nanopore PromethION platforms to sequence, de novo assemble, and annotate the whole genome and mitochondrial genome (mitogenome) of the A. bitorquis strain BH01 isolated from Bosten Lake, Xinjiang Province, China. Using the genome-based biological information, we identified candidate genes associated with mating type and carbohydrate-active enzymes in A. bitorquis. Cluster analysis based on P450 of basidiomycetes revealed the types of P450 members of A. bitorquis. Comparative genomic, mitogenomic, and phylogenetic analyses were also performed, revealing interspecific differences and evolutionary features of A. bitorquis and A. bisporus. In addition, the molecular network of metabolites was investigated, highlighting differences in the chemical composition and content of the fruiting bodies of A. bitorquis and A. bisporus. The genome sequencing provides a comprehensive understanding and knowledge of A. bitorquis and the genus Agaricus mushrooms. This work provides valuable insights into the potential for artificial cultivation and molecular breeding of A. bitorquis, which will facilitate the development of A. bitorquis in the field of edible mushrooms and functional food manufacture.

A de novo assembled genome of the Tibetan Partridge (Perdix hodgsoniae) and its high-altitude adaptation

The Tibetan Partridge (Perdix hodgsoniae) is an endemic species distributed in high-altitude areas of 3600-5600 m on the Qinghai-Tibet Plateau. To explore how the species is adapted to the high elevation environment, we assembled a draft genome based on both the Illumina and PacBio sequencing platforms with its population genetics and genomics analysis. In total, 134.74 Gb short reads and 30.81 Gb long reads raw data were generated. The 1.05-Gb assembled genome had a contig N50 of 4.56 Mb, with 91.94% complete BUSCOs. The 17 457 genes were annotated, and 11.35% of the genome was composed of repeat sequences. The phylogenetic tree showed that P. hodgsoniae was located at the basal position of the clade, including Golden Pheasant (Chrysolophus pictus), Common Pheasant (Phasianus colchicus), and Mikado Pheasant (Syrmaticus mikado). We found that 1014, 2595, and 2732 of the 6641 one-to-one orthologous genes were under positive selection in P. hodgsoniae, detected using PAML, BUSTED, and aBSREL programs, respectively, of which 965 genes were common under positive selection with 3 different programs. Several positively selected genes and immunity pathways relevant to high-altitude adaptation were detected. Gene family evolution showed that 99 gene families experienced significant expansion events, while 6 gene families were under contraction. The total number of olfactory receptor genes was relatively low in P. hodgsoniae. Genomic data provide an important resource for a further study on the evolutionary history of P. hodgsoniae, which provides a new insight into its high-altitude adaptation mechanisms.

Phenotypic, Genomic, and Transcriptomic Comparison of Industrial Aspergillus oryzae Used in Chinese and Japanese Soy Sauce: Analysis of Key Proteolytic Enzymes Produced by Koji Molds

Aspergillus oryzae, which generates numerous enzymes for the breakdown of raw materials, is an essential koji mold in soy sauce production. For better soy sauce productivity and flavor quality, China and Japan have developed their own industrial A. oryzae strains at distinct evolutionary branches for use in soy sauce production for decades. However, systematic comparison between the two national industrial strains has been poorly conducted, and thus we have not been able to generate adequate knowledge, especially regarding what are the key hydrolytic enzymes produced by A. oryzae during koji production. This study sequenced and assembled three high-quality genome sequences of industrial A. oryzae originating from China and Japan. Based on the genome sequences, a phylogenetic tree analysis was performed and revealed the evolutional distances between the two national industrial koji molds. Meanwhile, a comparative phenotypic analysis revealed that the two national industrial strains differed in growth and catalytic characteristics, particularly in proteolytic enzyme activities. To investigate the molecular mechanism underlying the phenotypic difference, we conducted systematic comparative genome and transcriptome investigations. We found minor differences in the quantity and diversity of proteolytic enzyme genes between Chinese and Japanese koji molds, while the protease secretion ratio and transcriptional level were dissimilar. We identified 58 potential important enzymes associated with high protein breakdown efficiency during industrial koji fermentation by combining comparative phenotypic and transcriptome data. More research is required to confirm the function of these putative key hydrolytic enzymes. IMPORTANCE Aspergillus oryzae is widely used as an industrial koji mold for soy sauce brewing due to its powerful raw material decomposition capability. Although various proteases in A. oryzae have been identified, it remains a challenge to find essential enzymes involved in soy sauce production. Generally, the industrial A. oryzae used in soy sauce brewing has excellent proteolytic activity. Based on this, we analyzed key proteolytic enzymes according to a comparison of the genome and transcriptome between three industrial strains. This study found little difference in gene numbers and mutations of proteolytic enzymes between three industrial A. oryzae strains. However, variations in protease secretion ratio and transcriptome were discovered between industrial strains. Based on that, we generated 58 candidate key proteolytic enzymes. This work comprehensively analyzed three industrial koji molds, revealing genome development under separate artificial domestication and helping in the study of key proteolytic enzymes during soy sauce production.

Genomic Reclassification and Phenotypic Characterization of Pseudomonas putida Strains Deposited in Japanese Culture Collections

Pseudomonas putida is a major species belonging to the genus Pseudomonas. Although several hundred strains of P. putida have been deposited in culture collections, they potentially differ from the genetically defined "true Pseudomonas putida" because many were classified as P. putida based on their phenotypic and metabolic characteristics. A phylogenetic ana-lysis based on the concatenated sequences of the 16S rRNA and rpoD genes revealed that 46 strains of P. putida deposited in Japanese culture collections were classified into nine operational taxonomic units (OTUs) and eleven singletons. The OTU7 strain produces N-acylhomoserine lactone as a quorum-sensing signal. One of the OTU7 strains, JCM 20066, exhibited a ppuI-rsaL-ppuR quorum-sensing system that controls biofilm formation and motility. The P. putida type strain JCM 13063^T and six other strains were classified as OTU4. Classification based on the calculation of whole-genome similarity revealed that three OTU4 strains, JCM 20005, 21368, and 13061, were regarded as the same species as JCM 13063^T and defined as true P. putida. When orthologous genes in the whole-genome sequences of true P. putida strains were screened, PP4_28660 from P. putida NBRC 14164^T (=JCM 13063^T) was present in all true P. putida genome sequences. The internal region of PP4_28660 was successfully amplified from all true P. putida strains using the specific primers designed in this study.

Chromosome-level genome assembly of goose provides insight into the adaptation and growth of local goose breeds

BACKGROUND

Anatidae contains numerous waterfowl species with great economic value, but the genetic diversity basis remains insufficiently investigated. Here, we report a chromosome-level genome assembly of Lion-head goose (Anser cygnoides), a native breed in South China, through the combination of PacBio, Bionano, and Hi-C technologies.

FINDINGS

The assembly had a total genome size of 1.19 Gb, consisting of 1,859 contigs with an N50 length of 20.59 Mb, generating 40 pseudochromosomes, representing 97.27% of the assembled genome, and identifying 21,208 protein-coding genes. Comparative genomic analysis revealed that geese and ducks diverged approximately 28.42 million years ago, and geese have undergone massive gene family expansion and contraction. To identify genetic markers associated with body weight in different geese breeds, including Wuzong goose, Huangzong goose, Magang goose, and Lion-head goose, a genome-wide association study was performed, yielding an average of 1,520.6 Mb of raw data that detected 44,858 single-mucleotide polymorphisms (SNPs). Genome-wide association study showed that 6 SNPs were significantly associated with body weight and 25 were potentially associated. The significantly associated SNPs were annotated as LDLRAD4, GPR180, and OR, enriching in growth factor receptor regulation pathways.

CONCLUSIONS

We present the first chromosome-level assembly of the Lion-head goose genome, which will expand the genomic resources of the Anatidae family, providing a basis for adaptation and evolution. Candidate genes significantly associated with different goose breeds may serve to understand the underlying mechanisms of weight differences.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

Comparative Genomics Reveal the Utilization Ability of Variable Carbohydrates as Key Genetic Features of Listeria Pathogens in Their Pathogenic Lifestyles

BACKGROUND

L. monocytogenes and L. ivanovii, the only two pathogens of Listeria, can survive in various environments, having different pathogenic characteristics. However, the genetic basis of their excellent adaptability and differences in pathogenicity has still not been completely elucidated.

METHODS

We performed a comparative genomic analysis based on 275 L. monocytogenes, 10 L. ivanovii, and 22 non-pathogenic Listeria strains.

RESULTS

Core/pan-genome analysis revealed that 975 gene families were conserved in all the studied strains. Additionally, 204, 242, and 756 gene families existed uniquely in L. monocytogenes, L. ivanovii, and both, respectively. Functional annotation partially verified that these unique gene families were closely related to their adaptability and pathogenicity. Moreover, the protein-protein interaction (PPI) network analysis of these unique gene sets showed that plenty of carbohydrate transport systems and energy metabolism enzymes were clustered in the networks. Interestingly, ethanolamine-metabolic-process-related proteins were significantly enriched in the PPI network of the unique genes of the Listeria pathogens, which can be understood as a determining factor of their pathogenicity.

CONCLUSIONS

The utilization capacity of multiple carbon sources of Listeria pathogens, especially ethanolamine, is the key genetic basis for their ability to adapt to various environments and pathogenic lifestyles.

Organelle Comparative Genome Analysis Reveals Novel Alloplasmic Male Sterility with orf112 in Brassica oleracea L

B. oleracea Ogura CMS is an alloplasmic male-sterile line introduced from radish by interspecific hybridization and protoplast fusion. The introduction of alien cytoplasm resulted in many undesirable traits, which affected the yield of hybrids. Therefore, it is necessary to identify the composition and reduce the content of alien cytoplasm in B. oleracea Ogura CMS. In the present study, we sequenced, assembled, and compared the organelle genomes of Ogura CMS cabbage and its maintainer line. The chloroplast genome of Ogura-type cabbage was completely derived from normal-type cabbage, whereas the mitochondrial genome was recombined from normal-type cabbage and Ogura-type radish. Nine unique regions derived from radish were identified in the mitochondrial genome of Ogura-type cabbage, and the total length of these nine regions was 35,618 bp, accounting for 13.84% of the mitochondrial genome. Using 32 alloplasmic markers designed according to the sequences of these nine regions, one novel sterile source with less alien cytoplasm was discovered among 305 materials and named Bel CMS. The size of the alien cytoplasm in Bel CMS was 21,587 bp, accounting for 8.93% of its mtDNA, which was much less than that in Ogura CMS. Most importantly, the sterility gene orf138 was replaced by orf112, which had a 78-bp deletion, in Bel CMS. Interestingly, Bel CMS cabbage also maintained 100% sterility, although orf112 had 26 fewer amino acids than orf138. Field phenotypic observation showed that Bel CMS was an excellent sterile source with stable 100% sterility and no withered buds at the early flowering stage, which could replace Ogura CMS in cabbage heterosis utilization.

Comparative genomic and transmission analysis of Clostridioides difficile between environmental, animal, and clinical sources in China

Clostridioides difficile is the most common pathogen causing antibiotic-associated diarrhea. Previous studies showed that diverse sources, aside from C. difficile infection (CDI) patients, played a major role in C. difficile hospital transmission. This study aimed to investigate relationships and transmission potential of C. difficile strains from different sources. A prospective study was conducted both in the intensive care unit (ICU) and six livestock farms in China in 2018–2019. Ninety-eight strains from CDI patients (10 isolates), asymptomatic hospitalized carriers (55), the ICU environment (12), animals (14), soil (4), and farmers (3) were collected. Sequence type (ST) 3/ribotype (RT) 001, ST35/RT046, and ST48/RT596 were dominant types, distributed widely in multiple sources. Core-genome single-nucleotide polymorphism (cgSNP) analysis showed that hospital and farm strains shared several common clonal groups (CGs, strains separated by ≤ 2 cgSNPs) (CG4/ST3/RT001, CG7/ST35/RT046, CG11/ST48/RT596). CDI patients, asymptomatic carriers, and the ICU environment strains also shared several common CGs. The number of virulence genes was not statistically different between strains from different sources. Multi-source strains in the same CG carried identical virulence gene sequences, including pathogenicity genes at the pathogenicity locus and adhesion-related genes at S-layer cassette. Resistance genes (ermB, tetM, etc.) were widespread in multiple sources, and multi-source strains in the same CG had similar resistance phenotypes and carried consistent transposons and plasmid types. The study indicated that interspecies and cross-regional transmission of C. difficile occurs between animals, the environment, and humans. Community-associated strains from both farms and asymptomatic hospitalized carriers were important reservoirs of CDI in hospitals.

Genome assembly of primitive cultivated potato Solanum stenotomum provides insights into potato evolution

Genetic diversity is the raw material for germplasm enhancement. Landraces and wild species relatives of potato, which contain a rich gene pool of valuable agronomic traits, can provide insights into the genetic diversity behind the adaptability of the common potato. The diploid plant, Solanum stenotomum (Sst), is believed to have an ancestral relationship with modern potato cultivars and be a potential source of resistance against disease. Sequencing of the Sst genome generated an assembly of 852.85 Mb (N50 scaffold size, 3.7 Mb). Pseudomolecule construction anchored 788.75 Mb of the assembly onto 12 pseudochromosomes, with an anchor rate of 92.4%. Genome annotation yielded 41,914 high-confidence protein-coding gene models and comparative analyses with closely related Solanaceae species identified 358 Sst-specific gene families, 885 gene families with expansion along the Sst lineage, and 149 genes experiencing accelerated rates of protein sequence evolution in Sst, the functions of which were mainly associated with defense responses, particularly against bacterial and fungal infection. Insights into the Sst genome and the genomic variation of cultivated potato taxa are valuable in elaborating the impact of potato evolution in early landrace diploid and facilitate modern potato breeding.

Molecular Characteristics of Novel Phage vB_ShiP-A7 Infecting Multidrug-Resistant Shigella flexneri and Escherichia coli, and Its Bactericidal Effect in vitro and in vivo

In recent years, increasing evidence has shown that bacteriophages (phages) can inhibit infection caused by multidrug-resistant (MDR) bacteria. Here, we isolated a new phage, named vB_ShiP-A7, using MDR Shigella flexneri as the host. vB_ShiP-A7 is a novel member of Podoviridae, with a latency period of approximately 35 min and a burst size of approximately 100 phage particles/cell. The adsorption rate constant of phage vB_ShiP-A7 to its host S. flexneri was 1.405 × 10^–8 mL/min. The vB_ShiP-A7 genome is a linear double-stranded DNA composed of 40,058 bp with 177 bp terminal repeats, encoding 43 putative open reading frames. Comparative genomic analysis demonstrated that the genome sequence of vB_ShiP-A7 is closely related to 15 different phages, which can infect different strains. Mass spectrometry analysis revealed that 12 known proteins and 6 hypothetical proteins exist in the particles of phage vB_ShiP-A7. Our results confirmed that the genome of vB_ShiP-A7 is free of lysogen-related genes, bacterial virulence genes, and antibiotic resistance genes. vB_ShiP-A7 can significantly disrupt the growth of some MDR clinical strains of S. flexneri and Escherichia coli in liquid culture and biofilms in vitro. In addition, vB_ShiP-A7 can reduce the load of S. flexneri by approximately 3–10 folds in an infection model of mice. Therefore, vB_ShiP-A7 is a stable novel phage with the potential to treat infections caused by MDR strains of S. flexneri and E. coli.

Genomic Signatures for Species-Specific Adaptation in Lake Victoria Cichlids Derived from Large-Scale Standing Genetic Variation

The cichlids of Lake Victoria are a textbook example of adaptive radiation, as >500 endemic species arose in just 14,600 years. The degree of genetic differentiation among species is very low due to the short period of time after the radiation, which allows us to ascertain highly differentiated genes that are strong candidates for driving speciation and adaptation. Previous studies have revealed the critical contribution of vision to speciation by showing the existence of highly differentiated alleles in the visual opsin gene among species with different habitat depths. In contrast, the processes of species-specific adaptation to different ecological backgrounds remain to be investigated. Here, we used genome-wide comparative analyses of three species of Lake Victoria cichlids that inhabit different environments-Haplochromis chilotes, H. sauvagei, and Lithochromis rufus-to elucidate the processes of adaptation by estimating population history and by searching for candidate genes that contribute to adaptation. The patterns of changes in population size were quite distinct among the species according to their habitats. We identified many novel adaptive candidate genes, some of which had surprisingly long divergent haplotypes between species, thus showing the footprint of selective sweep events. Molecular phylogenetic analyses revealed that a large fraction of the allelic diversity among Lake Victoria cichlids was derived from standing genetic variation that originated before the adaptive radiation. Our analyses uncovered the processes of species-specific adaptation of Lake Victoria cichlids and the complexity of the genomic substrate that facilitated this adaptation.

Selection of Lactococcus lactis HY7803 for Glutamic Acid Production Based on Comparative Genomic Analysis

Comparative genomic analysis was performed on eight species of lactic acid bacteria (LAB)-Lactococcus (L.) lactis, Lactobacillus (Lb.) plantarum, Lb. casei, Lb. brevis, Leuconostoc (Leu.) mesenteroides, Lb. fermentum, Lb. buchneri, and Lb. curvatus-to assess their glutamic acid production pathways. Glutamic acid is important for umami taste in foods. The only genes for glutamic acid production identified in the eight LAB were for conversion from glutamine in L. lactis and Leu. mesenteroides, and from glucose via citrate in L. lactis. Thus, L. lactis was considered to be potentially the best of the species for glutamic acid production. By biochemical analyses, L. lactis HY7803 was selected for glutamic acid production from among 17 L. lactis strains. Strain HY7803 produced 83.16 pmol/μl glutamic acid from glucose, and exogenous supplementation of citrate increased this to 108.42 pmol/μl. Including glutamic acid, strain HY7803 produced more of 10 free amino acids than L. lactis reference strains IL1403 and ATCC 7962 in the presence of exogenous citrate. The differences in the amino acid profiles of the strains were illuminated by principal component analysis. Our results indicate that L. lactis HY7803 may be a good starter strain for glutamic acid production.

Analysis of the Antifungal Properties of Bacillus velezensis B-4 Through a Bioassay and Complete-Genome Sequencing

The strain B-4, isolated from a field in Changsha (China), presents strong antifungal activities, as identified by the Kirby–Bauer test, especially for pathogens that harm crops. Here, we obtained the complete genome sequence of the strain B-4 by Pacific Biosciences single-molecule real-time sequencing, making it well analyzed for understanding mechanisms and creating biological agents. Its 3,919-kb circular chromosome genome has 3,725 protein-coding genes [coding sequences (CDSs)] and 46.7% guanine–cytosine content. A comparative genome analysis of B-4 with other published strains (including Bacillus velezensis, Bacillus amyloliquefaciens, and Bacillus subtilis) revealed that the strain B-4 is a B. velezensis strain. These different strains have 2,889 CDSs in common, whereas 179 CDSs were found to be unique in the strain B-4, which is a far greater number than that in other strains. Regarding the antifungal activities of B-4, we were specifically concerned with the genes involved in the biosynthesis of secondary metabolites. In total, more than 19.56% of the genome was annotated to 12 gene clusters relating to synthesis of antimicrobial metabolites, which contained various enzyme-encoding operons for non-ribosomal peptide synthetases, polyketide synthases, and lantipeptide synthesis proteins. They were all considered to be related to the production of bacteriostatic substances or stimulation of induced systemic resistance by bacterial metabolites. These situations also present an advantage over those of other strains for biocontrol potential. We provide evidence that the biological control effect of the strain B-4, as demonstrated in antibacterial activity experiments and predicted from the complete genome sequence analysis, provides the basis for research promoting agricultural research on sustainable development, especially the contribution of biotechnology to agriculture.

Comparative genome analysis and mining of secondary metabolites of Paenibacillus polymyxa

Paenibacillus polymyxa is a well-known Gram-positive biocontrol bacterium. It has been reported that many P. polymyxa strains can inhibit bacteria, fungi and other plant pathogens. Paenibacillus polymyxa employs a variety of mechanisms to promote plant growth, so it is necessary to understand the biocontrol ability of bacteria at the genome level. In the present study, thanks to the widespread availability of Paenibacillus genome data and the development of bioinformatics tools, we were able to analyze and mine the genomes of 43 P. polymyxa strains. The strain NCTC4744 was determined not to be P. polymyxa according to digital DNA-DNA hybridization and average nucleotide identity. By analysis of the pan-genome and the core genome, we found that the pan-genome of P. polymyxa was open and that there were 3,192 core genes. In a gene cluster analysis of secondary metabolites, 797 secondary metabolite gene clusters were found, of which 343 are not similar to known clusters and are expected to reveal a large number of new secondary metabolites. We also analyzed the plant growth-promoting genes that were mined and found, surpisingly, that these genes are highly conserved. The results of the present study not only reveal a large number of unknown potential secondary metabolite gene clusters in P. polymyxa, but also suggest that plant growth promotion characteristics are evolutionary adaptations of P. polymyxa to plant-related habitats.

Genomic profile of Brazilian methicillin-resistant Staphylococcus aureus resembles clones dispersed worldwide

PURPOSE

Comparative genomic analysis of strains may help us to better understand the wide diversity of their genetic profiles. The aim of this study was to analyse the genomic features of the resistome and virulome of Brazilian first methicillin-resistant Staphylococcus aureus (MRSA) isolates and their relationship to other Brazilian and international MRSA strains.

METHODOLOGY

The whole genomes of three MRSA strains previously isolated in Vitória da Conquista were sequenced, assembled, annotated and compared with other MRSA genomes. A phylogenetic tree was constructed and the pan-genome and accessory and core genomes were constructed. The resistomes and virulomes of all strains were identified.Results/Key findings. Phylogenetic analysis of all 49 strains indicated different clones showing high similarity. The pan-genome of the analysed strains consisted of 4484 genes, with 31 % comprising the gene portion of the core genome, 47 % comprising the accessory genome and 22 % being singletons. Most strains showed at least one gene related to virulence factors associated with immune system evasion, followed by enterotoxins. The strains showed multiresistance, with the most recurrent genes conferring resistance to beta-lactams, fluoroquinolones, aminoglycosides and macrolides.

CONCLUSIONS

Our comparative genomic analysis showed that there is no pattern of virulence gene distribution among the clones analysed in the different regions. The Brazilian strains showed similarity with clones from several continents.

Characterization a Novel Butyric Acid-Producing Bacterium Collinsella aerofaciens Subsp. Shenzhenensis Subsp. Nov. Microorganisms. 2019;7. [PMID: 30871249 PMCID: PMC6463082 DOI: 10.3390/microorganisms7030078]

Butyrate-producing bacteria can biosynthesize butyrate and alleviate inflammatory diseases. However, few studies have reported that the genus Collinsella has the ability to produce butyric acid. Here, our study depicts a Collinsella strain, which is a rod-shaped obligate anaerobe that is able to produce butyric acid. This microorganism was isolated from a human gut, and the optimal growth conditions were found to be 37 °C on PYG medium with pH 6.5. The 16S rRNA gene sequence demonstrated that this microorganism shared 99.93% similarity with C. aerofaciens ATCC 25986^T, which was higher than the threshold (98.65%) for differentiating two species. Digital DNA⁻DNA hybridization and average nucleotide identity values also supported that this microorganism belonged to the species C. aerofaciens. Distinct phenotypic characteristics between TF06-26 and the type strain of C. aerofaciens, such as the fermentation of D-lactose, D-fructose and D-maltose, positive growth under pH 5 and 0.2% (w/v) cholate, suggested this strain was a novel subspecies. Comparative genome analysis revealed that butyric acid kinase and phosphate butyryltransferase enzymes were coded exclusively by this strain, indicating a specific butyric acid-producing function of this C. aerofaciens subspecies within the genus Collinsella. Thus, Collinsella aerofaciens subsp. shenzhenensis subsp. nov. was proposed, with set strain TF06-26^T (=CGMCC 1.5216^T = DSM 105138^T) as the type strain.

Ancrocorticia populi gen. nov., sp. nov, isolated from the symptomatic bark of Populus × euramericana canker

A Gram-staining positive facultative anaerobic, non-motile strain, sk1b4^T , was isolated from canker of symptomatic bark tissue of a Populus × euramericana. 16S rRNA gene sequence analyses showed that strain sk1b4^T shared the highest similarity with Arcanobacterium phocisimile (94.1%). Within the phylogenetic tree, the novel isolate formed a distinct branch from Actinobaculum, Arcanobacterium, and Trueperella. The percentage of conserved proteins calculated from genomic sequence indicated a low level of relatedness between the novel strain and its phylogenetic neighbors. Growth of the novel strain occurred at temperatures between 10 and 41°C, and within a pH range of 6.0-9.0; optimal growth occurred at 30°C and at pH 6.0-9.0. Growth also occurred within a NaCl concentration of 1%-5% (w/v). The major fatty acids of the strain were C_14:0 , C_16:0 , and C_18:1 ω9c_, and major polar lipids were glycolipid, phosphatidylinositol mannoside, phospholipid, diphosphatidylglycerol, and phosphatidylglycerol. Respiratory quinone was absent. On the basis of phenotypic and genotypic characteristics, we propose that the novel isolate should be classified as a novel species in a new genus: Ancrocorticia populi gen. nov., sp. nov. The type strain is sk1b4^T (=CFCC 14564^T = KCTC 39919^T ).

Virulence and genomic features of a bla CTX-M-3 and bla CTX-M-14 coharboring hypermucoviscous Klebsiella pneumoniae of serotype K2 and ST65

Background

Capsular serotype K2 Klebsiella pneumoniae of sequence type (ST) 65 has been recognized as a hypervirulent clone. Simultaneous presence of different bla_CTX-M genes has never been reported in this clone. In the present study, the genetic characteristics and virulence phenotype of a CTX-M-3 and CTX-M-14 coproducing ST65 K. pneumoniae human isolate, KP_06, that caused an intracranial infection, are evaluated.

Methods

The potential virulence of KP_06 was assayed by in vitro and in vivo methods. The molecular biology and whole-genome sequencing technology were used to analyze the genomic features associated with the virulence of this strain.

Results

The KP_06 exhibited typical features of hypervirulent K. pneumoniae (hvKP), showing hypermucoviscosity phenotype and belonging to K2 and ST65. Apart from virulence genes linked to hvKP, including rmpA, rmpA2, and clb cluster and genes encoding siderophores, it was found to harbor a ~170 kb pLVPK-like virulence plasmid. In contrast to most hvKP, KP_06 was resistant to cephalosporins and the coexistence of bla_CTX-M-3 and bla_CTX-M-14 was detected. Further experiments demonstrated that this strain was classified as a nonbiofilm producer and serum sensitivity (grade 1) and killed only 30% of Galleria mellonella inoculated with 1×10⁶ colony-forming unit of the specimen within 48 hours, suggesting relatively low virulence. Comparative genomic analysis of KP_06 with five K2 hypermucoviscous K. pneumoniae (HMKP) revealed seven unique orthologies with varied function in this strain. Intriguingly, the virulence genes identified in KP-06 were unexpectedly more diverse than those observed in five other K2 HMKP strains.

Conclusion

Our data support the notion that neither virulence-associated genes (clusters) nor the pLVPK-like virulence plasmid is sufficient for the hypervirulence of K. pneumoniae. Future studies aiming to explore the virulence of K. pneumoniae should take genome-based profile together with experimental work. The detailed mechanism involving in the impaired virulence of KP_06 remains to be further explored.

The Infection Process of Yersinia ruckeri: Reviewing the Pieces of the Jigsaw Puzzle

Finding the keys to understanding the infectious process of Yersinia ruckeri was not a priority for many years due to the prompt development of an effective biotype 1 vaccine which was used mainly in Europe and USA. However, the gradual emergence of outbreaks in vaccinated fish, which have been reported since 2003, has awakened interest in the mechanism of virulence in this pathogen. Thus, during the last two decades, a large number of studies have considerably enriched our knowledge of many aspects of the pathogen and its interaction with the host. By means of both conventional and a variety of novel strategies, such as cell GFP tagging, bioluminescence imaging and optical projection tomography, it has been possible to determine three putative Y. ruckeri infection routes, the main point of entry for the bacterium being the gill lamellae. Moreover, a wide range of potential virulence factors have been highlighted by specific gene mutagenesis strategies or genome-wide transposon/plasmid insertion-based screening approaches, such us in vivo expression technology (IVET) and signature tagged mutagenesis (STM). Finally, recent proteomic and whole genomic analyses have allowed many of the genes and systems that are potentially implicated in the organism's pathogenicity and its adaptation to the host environmental conditions to be elucidated. Altogether, these studies contribute to a better understanding of the infectious process of Y. ruckeri in fish, which is crucial for the development of more effective strategies for preventing or treating enteric redmouth disease (ERM).

Characterization and Genomic Study of Phage vB_EcoS-B2 Infecting Multidrug-Resistant Escherichia coli

The potential of bacteriophage as an alternative antibacterial agent has been reconsidered for control of pathogenic bacteria due to the widespread occurrence of multi-drug resistance bacteria. More and more lytic phages have been isolated recently. In the present study, we isolated a lytic phage named vB_EcoS-B2 from waste water. VB_EcoS-B2 has an icosahedral symmetry head and a long tail without a contractile sheath, indicating that it belongs to the family Siphoviridae. The complete genome of vB_EcoS-B2 is composed of a circular double stranded DNA of 44,283 bp in length, with 54.77% GC content. vB_EcoS-B2 is homologous to 14 relative phages (such as Escherichia phage SSL-2009a, Escherichia phage JL1, and Shigella phage EP23), but most of these phages exhibit different gene arrangement. Our results serve to extend our understanding toward phage evolution of family Siphoviridae of coliphages. Sixty-five putative open reading frames were predicted in the complete genome of vB_EcoS-B2. Twenty-one of proteins encoded by vB_EcoS-B2 were determined in phage particles by Mass Spectrometry. Bacteriophage genome and proteome analysis confirmed the lytic nature of vB_EcoS-B2, namely, the absence of toxin-coding genes, islands of pathogenicity, or genes through lysogeny or transduction. Furthermore, vB_EcoS-B2 significantly reduced the growth of E. coli MG1655 and also inhibited the growth of several multi-drug resistant clinical stains of E. coli. Phage vB_EcoS-B2 can kill some of the MRD E. coli entirely, strongly indicating us that it could be one of the components of phage cocktails to treat multi-drug resistant E. coli. This phage could be used to interrupt or reduce the spread of multi-drug resistant E. coli.

Comparative genomic analyses reveal the features for adaptation to nematodes in fungi

Nematophagous (NP) fungi are ecologically important components of the soil microbiome in natural ecosystems. Esteya vermicola (Ev) has been reported as a NP fungus with a poorly understood evolutionary history and mechanism of adaptation to parasitism. Furthermore, NP fungal genomic basis of lifestyle was still unclear. We sequenced and annotated the Ev genome (34.2 Mbp) and integrated genetic makeup and evolution of pathogenic genes to investigate NP fungi. The results revealed that NP fungi had some abundant pathogenic genes corresponding to their niche. A number of gene families involved in pathogenicity were expanded, and some pathogenic orthologous genes underwent positive selection. NP fungi with diverse morphological features exhibit similarities of evolutionary convergence in attacking nematodes, but their genetic makeup and microscopic mechanism are different. Endoparasitic NP fungi showed similarity in large number of transporters and secondary metabolite coding genes. Noteworthy, expanded families of transporters and endo-beta-glucanase implied great genetic potential of Ev in quickly perturbing nematode metabolism and parasitic behavior. These results facilitate our understanding of NP fungal genomic features for adaptation to nematodes and lay a solid theoretical foundation for further research and application.

Comparative Genome Analysis of Rathayibacter tritici NCPPB 1953 with Rathayibacter toxicus Strains Can Facilitate Studies on Mechanisms of Nematode Association and Host Infection

Rathayibacter tritici, which is a Gram positive, plant pathogenic, non-motile, and rod-shaped bacterium, causes spike blight in wheat and barley. For successful pathogenesis, R. tritici is associated with Anguina tritici, a nematode, which produces seed galls (ear cockles) in certain plant varieties and facilitates spread of infection. Despite significant efforts, little research is available on the mechanism of disease or bacteria-nematode association of this bacterium due to lack of genomic information. Here, we report the first complete genome sequence of R. tritici NCPPB 1953 with diverse features of this strain. The whole genome consists of one circular chromosome of 3,354,681 bp with a GC content of 69.48%. A total of 2,979 genes were predicted, comprising 2,866 protein coding genes and 49 RNA genes. The comparative genomic analyses between R. tritici NCPPB 1953 and R. toxicus strains identified 1,052 specific genes in R. tritici NCPPB 1953. Using the BlastKOALA database, we revealed that the flexible genome of R. tritici NCPPB 1953 is highly enriched in 'Environmental Information Processing' system and metabolic processes for diverse substrates. Furthermore, many specific genes of R. tritici NCPPB 1953 are distributed in substrate-binding proteins for extracellular signals including saccharides, lipids, phosphates, amino acids and metallic cations. These data provides clues on rapid and stable colonization of R. tritici for disease mechanism and nematode association.

CsSNP: A Web-Based Tool for the Detecting of Comparative Segments SNPs

SNP (single nucleotide polymorphism) is a popular tool for the study of genetic diversity, evolution, and other areas. Therefore, it is necessary to develop a convenient, utility, robust, rapid, and open source detecting-SNP tool for all researchers. Since the detection of SNPs needs special software and series steps including alignment, detection, analysis and present, the study of SNPs is limited for nonprofessional users. CsSNP (Comparative segments SNP, http://biodb.sdau.edu.cn/cssnp/ ) is a freely available web tool based on the Blat, Blast, and Perl programs to detect comparative segments SNPs and to show the detail information of SNPs. The results are filtered and presented in the statistics figure and a Gbrowse map. This platform contains the reference genomic sequences and coding sequences of 60 plant species, and also provides new opportunities for the users to detect SNPs easily. CsSNP is provided a convenient tool for nonprofessional users to find comparative segments SNPs in their own sequences, and give the users the information and the analysis of SNPs, and display these data in a dynamic map. It provides a new method to detect SNPs and may accelerate related studies.