Genome assembly of M. spongiola and comparative genomics of the genus Morchella provide initial insights into taxonomy and adaptive evolution.

Morchella spongiola is a highly prized mushroom for its delicious flavor and medical value and is one of the most flourishing, representative, and dominant macrofungi in the Qilian Mountains of the Qinghai-Tibet Plateau subkingdoms (QTPs). However, the understanding of M. spongiola remains largely unknown, and its taxonomy is ambiguous. In this study, we redescribed a unique species of M. spongiola, i.e., micromorphology, molecular data, genomics, and comparative genomics, and the historical biogeography of M. spongiola were estimated for 182 single-copy homologous genes. A high-quality chromosome-level reference genome of M. spongiola M12-10 was obtained by combining PacBio HiFi data and Illumina sequencing technologies; it was approximately 57.1 Mb (contig N50 of 18.14 Mb) and contained 9775 protein-coding genes. Comparative genome analysis revealed considerable conservation and unique characteristics between M. spongiola M12-10 and 32 other Morchella species. Molecular phylogenetic analysis indicated that M. spongiola M12-10 is similar to the M. prava/Mes-7 present in sandy soil near rivers, differentiating from black morels ~ 43.06 Mya (million years ago), and diverged from M. parva/Mes-7 at approximately 12.85 Mya (in the Miocene epoch), which is closely related to the geological activities in the QTPs (in the Neogene). Therefore, M. spongiola is a unique species rather than a synonym of M. vulgaris/Mes-5, which has a distinctive grey-brown sponge-like ascomata. This genome of M. spongiola M12-10 is the first published genome sequence of the species in the genus Morchella from the QTPs, which could aid future studies on functional gene identification, germplasm resource management, and molecular breeding efforts, as well as evolutionary studies on the Morchella taxon in the QTPs.

Isolation of a Virulent Clostridium perfringens Strain from Elaphurus davidianus and Characterization by Whole-Genome Sequence Analysis.

Clostridium perfringens (C. perfringens) is an important veterinary pathogen and a noteworthy threat to human and animal health. Recently, there has been a significant rise in the number of moose fatalities caused by this rare, endemic species in China. Currently, there is an increasing trend in conducting whole-genome analysis of C. perfringens strains originating from pigs and chickens, whereas fewer studies have been undertaken on Elaphurus davidianus-originating strains at the whole-genome level. Our laboratory has identified and isolated five C. perfringens type A from affected Elaphurus davidianus. The current study identified the most potent strain of C. perfringens, which originated from Elaphurus davidianus, and sequenced its genome to reveal virulence genes and pathogenicity. Our findings show that strain CX1-4 exhibits the highest levels of phospholipase activity, hemolytic activity, and mouse toxicity compared to the other four isolated C. perfringens type A strains. The chromosome sequence length of the CX1-4 strain was found to be 3,355,389 bp by complete genome sequencing. The current study unveils the genomic characteristics of C. perfringens type A originating from Elaphurus davidianus. It provides a core foundation for further investigation regarding the prevention and treatment of such infectious diseases in Elaphurus davidianus.

Phylogenetic Analysis and Comparative Genomics of Brucella abortus and Brucella melitensis Strains in Egypt.

Brucellosis is a notifiable disease induced by a facultative intracellular Brucella pathogen. In this study, eight Brucella abortus and eighteen Brucella melitensis strains from Egypt were annotated and compared with RB51 and REV1 vaccines respectively. RAST toolkit in the BV-BRC server was used for annotation, revealing genome length of 3,250,377 bp and 3,285,803 bp, 3289 and 3323 CDS, 48 and 49 tRNA genes, the same number of rRNA (3) genes, 583 and 586 hypothetical proteins, 2697 and 2726 functional proteins for B. abortus and B. melitensis respectively. B. abortus strains exhibit a similar number of candidate genes, while B. melitensis strains showed some differences, especially in the SRR19520422 Faiyum strain. Also, B. melitensis clarified differences in antimicrobial resistance genes (KatG, FabL, MtrA, MtrB, OxyR, and VanO-type) in SRR19520319 Faiyum and (Erm C and Tet K) in SRR19520422 Faiyum strain. Additionally, the whole genome phylogeny analysis proved that all B. abortus strains were related to vaccinated animals and all B. melitensis strains of Menoufia clustered together and closely related to Gharbia, Dameitta, and Kafr Elshiek. The Bowtie2 tool identified 338 (eight B. abortus) and 4271 (eighteen B. melitensis) single nucleotide polymorphisms (SNPs) along the genomes. These variants had been annotated according to type and impact. Moreover, thirty candidate genes were predicted and submitted at GenBank (24 in B. abortus) and (6 in B. melitensis). This study contributes significant insights into genetic variation, virulence factors, and vaccine-related associations of Brucella pathogens, enhancing our knowledge of brucellosis epidemiology and evolution in Egypt.

Whole genome sequencing and analysis of Bacillus sp. TTMP2, a tetramethylpyrazine-producing bacterium.

BACKGROUND: Tetramethylpyrazine has been extensively studied as an anticancer substance and a flavor substance in the fields of medicine and food industry. A strain with high tetramethylpyrazine production was screened from the fermented grains of Danquan winery. Genome sequencing can reveal the potential roles of bacteria by thoroughly examining the connection between genes and phenotypes from a genomic perspective. METHODS AND RESULTS: In this study, whole genome of this strain was sequenced and analyzed. This paper summarized the genomic characteristics of strain TTMP2 and analyzed genes related to the synthesis of tetramethylpyrazine. Bacillus sp. TTMP2 has a complete metabolic pathway for acetoin and tetramethylpyrazine metabolism. Gene function was analyzed by COG annotation, GO annotation, KEGG annotation and functional annotations for lipoproteins, carbohydrate-active enzymes, and pathogen-host interactions. Phylogenetic analysis indicated that Bacillus velezensis had the high homology with Bacillus sp. TTMP2. Genomes of 16 Bacillus species cover all genes of Bacillus, suggesting that genus Bacillus has an open pan-genome and can survive in diverse environments. CONCLUSION: The analysis of genome sequencing data from Bacillus sp. TTMP2 showed that its metabolic characteristics could be deeply understood, indicating that this bacterium had a particular role in tetramethylpyrazine synthesis.

Beyond rational-biosensor-guided isolation of 100 independently evolved bacterial strain variants and comparative analysis of their genomes.

BACKGROUND: In contrast to modern rational metabolic engineering, classical strain development strongly relies on random mutagenesis and screening for the desired production phenotype. Nowadays, with the availability of biosensor-based FACS screening strategies, these random approaches are coming back into fashion. In this study, we employ this technology in combination with comparative genome analyses to identify novel mutations contributing to product formation in the genome of a Corynebacterium glutamicum L-histidine producer. Since all known genetic targets contributing to L-histidine production have been already rationally engineered in this strain, identification of novel beneficial mutations can be regarded as challenging, as they might not be intuitively linkable to L-histidine biosynthesis. RESULTS: In order to identify 100 improved strain variants that had each arisen independently, we performed > 600 chemical mutagenesis experiments, > 200 biosensor-based FACS screenings, isolated > 50,000 variants with increased fluorescence, and characterized > 4500 variants with regard to biomass formation and L-histidine production. Based on comparative genome analyses of these 100 variants accumulating 10-80% more L-histidine, we discovered several beneficial mutations. Combination of selected genetic modifications allowed for the construction of a strain variant characterized by a doubled L-histidine titer (29 mM) and product yield (0.13 C-mol C-mol-1) in comparison to the starting variant. CONCLUSIONS: This study may serve as a blueprint for the identification of novel beneficial mutations in microbial producers in a more systematic manner. This way, also previously unexplored genes or genes with previously unknown contribution to the respective production phenotype can be identified. We believe that this technology has a great potential to push industrial production strains towards maximum performance.

Genome Mining of Novel Targets and Construction of Ladder-shaped melting temperature isothermal amplification Assays for the Identification of Cronobacter sakazakii and Cronobacter malonaticus.

Cronobacter species are potential pathogens that can contaminate powdered infant formula. C. sakazakii and C. malonaticus are the most common species of Cronobacter associated with infections. This study mined new molecular targets for the detection of C. sakazakii and C. malonaticus by using comparative genome approaches. Specific target genes mngB and ompR were obtained and used to detect C. sakazakii and C. malonaticus, respectively. A novel detection method, termed ladder-shape melting temperature isothermal amplification (LMTIA), was developed and evaluated. The detection limit for pure C. sakazakii DNA was 1 pg per reaction and 1 fg per reaction for C. malonaticus. The C. sakazakii, C. malonaticus, and the reference stains were all correctly identified. The amplicons can be successfully visualized and identified by naked eyes when hydroxy naphthol blue dye (HNB dye) was used in the reaction. Therefore, the LMTIA assays developed in this study showed potential application for microorganism identification and detection.

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

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

Genomic characterization and identification of virulence-related genes in Vibrio nigripulchritudo isolated from white leg shrimp Penaeus vannamei.

Vibrio nigripulchritudo causes vibriosis in penaeid shrimps. Here, we used Illumina and Nanopore sequencing technologies to sequence the genomes of three of its strains (TUMSAT-V. nig1, TUMSAT-V. nig2, and TUMSAT-V. nig3) to explore opportunities for disease management. Putative virulence factors and mobile genetic elements were detected while evaluating the phylogenetic relationship of each isolated strain. The genomes consisted of two circular chromosomes (I and II) plus one or two plasmids. The size of chromosome I ranged from 4.02 to 4.07 Mb with an average GC content of 46%, while the number of predicted CDSs ranged from 3563 to 3644. The size of chromosome II ranged from 2.16 to 2.18 Mb, with an average GC content of 45.5%, and the number of predicted CDSs ranged from 1970 to 1987. Numerous virulence genes were identified related to adherence, antiphagocytosis, chemotaxis, motility, iron uptake, quorum sensing, secretion systems, and toxins in all three genomes. Higher numbers of prophages and genomic islands found in TUMSAT-V. nig1 suggest that the strain has experienced numerous horizontal gene transfer events. The presence of antimicrobial resistance genes suggests that the strains have multidrug resistance. Comparative genomic analysis showed that all three strains belonged to the same clade.

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.

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.

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.

The Complete Chloroplast Genome Sequence of Cicer bijugum, Genome Organization, and Comparison with Related Species.

Background: Chickpea is one of Turkey's most significant legumes, and because of its high nutritional value, it is frequently preferred in human nourishment.Chloroplasts, which have their own genetic material, are organelles responsible for photosynthesis in plant cells and their genome contains non-trivial information about the molecular features and evolutionary process of plants. Objective: Current study aimed at revealing complete chloroplast genome sequence of one of the wild type Cicer species, Cicer bijugum, and comparing its genome with cultivated Cicer species, Cicer arietinum, by using bioinformatics analysis tools. Except for Cicer arietinum, there has been no study on the chloroplast genome sequence of Cicer species.Therefore, we targeted to reveal the complete chloroplast genome sequence of wild type Cicer species, Cicer bijugum, and compare the chloroplast genome of Cicer bijugum with the cultivated one Cicer arietinum. Methods: In this study, we sequenced the whole chloroplast genome of Cicer bijugum, one of the wild types of chickpea species, with the help Next Generation Sequencing platform and compared it with the chloroplast genome of the cultivated chickpea species, Cicer arietinum, by using online bioinformatics analysis tools. Results: We determined the size of the chloroplast genome of C. bijugum as 124,804 bp and found that C. bijugum did not contain an inverted repeat region in its chloroplast genome. Comparative analysis of the C. bijugum chloroplast genome uncovered thirteen hotspot regions (psbA, matK, rpoB, rpoC1, rpoC2, psbI, psbK, accD, rps19, ycf2, ycf1, rps15, and ndhF) and seven of them (matK, accD, rps19, ycf1, ycf2, rps15 and ndhF) could potentially be used as strong molecular markers for species identification. It has been determined that C. bijugum was phylogenetically closer to cultivated chickpea as compared to the other species. Conclusion: It is aimed that the data obtained from this study, which is the first study in which whole chloroplast genomes of wild chickpea species were sequenced, will guide researchers in future molecular, evolutionary, and genetic engineering studies with chickpea species.

Genome analysis of clinical isolate of Campylobacter fetus subspecies fetus MMM01 from India reveals genetic determinants of pathogenesis and adaptation.

In this study we report the whole genome sequencing (WGS) based analysis of blood-borne Campylobacter fetus subsp. fetus MMM01 isolated from a diabetic patient to obtain deeper insights in to the virulence and host adaptability. The sequenced genome of C. fetus subsp. fetus MMM01 along with reference genomes retrieved from NCBI was subjected to various in-silico analysis including JSpecies, MLST server, PATRIC server, VFanalyzer, CARD, PHASTER to understand their phylogenetic relation, virulence and antimicrobial resistance profile. The genome had a size of 1,788,790 bp, with a GC content of 33.09%, nearly identical to the reference strain C. fetus subsp. fetus 82-40. The MLST based phylogenetic tree constructed revealed the polyphyletic branching and MMM01 (ST25) was found to be closely related to ST11, both belong to the sap-A serotype which are more common in human infections. VFanalyzer identified 88 protein-coding genes coding for several virulence factors including Campylobacter adhesion to fibronectin, flagellar apparatus, cytolethal distending toxin operons and Campylobacter invasion antigen proteins which enhance the virulence of bacteria along with resistance genes against antibiotics including fluoroquinolone, chloramphenicol, tetracycline, and aminoglycoside in MMM01, which points to enhanced survival and pathogenicity of this zoonotic pathogen. It was interesting to find that MMM01 lacked FGI-II island found in most of the clinical isolates, which encoded CRISPR Cas and prophage II regions. More details about the complexity and evolution of this zoonotic pathogen could be learned from future studies that concentrate on comparative genome analysis using larger genome datasets.

Assembly and Annotation of Red Spruce (Picea rubens) Chloroplast Genome, Identification of Simple Sequence Repeats, and Phylogenetic Analysis in Picea.

We have sequenced the chloroplast genome of red spruce (Picea rubens) for the first time using the single-end, short-reads (44 bp) Illumina sequences, assembled and functionally annotated it, and identified simple sequence repeats (SSRs). The contigs were assembled using SOAPdenovo2 following the retrieval of chloroplast genome sequences using the black spruce (Picea mariana) chloroplast genome as the reference. The assembled genome length was 122,115 bp (gaps included). Comparatively, the P. rubens chloroplast genome reported here may be considered a near-complete draft. Global genome alignment and phylogenetic analysis based on the whole chloroplast genome sequences of Picea rubens and 10 other Picea species revealed high sequence synteny and conservation among 11 Picea species and phylogenetic relationships consistent with their known classical interrelationships and published molecular phylogeny. The P. rubens chloroplast genome sequence showed the highest similarity with that of P. mariana and the lowest with that of P. sitchensis. We have annotated 107 genes including 69 protein-coding genes, 28 tRNAs, 4 rRNAs, few pseudogenes, identified 42 SSRs, and successfully designed primers for 26 SSRs. Mononucleotide A/T repeats were the most common followed by dinucleotide AT repeats. A similar pattern of microsatellite repeats occurrence was found in the chloroplast genomes of 11 Picea species.

Roles of Species-Specific Legumains in Pathogenicity of the Pinewood Nematode Bursaphelenchus xylophilus.

Peptidases are very important to parasites, which have central roles in parasite biology and pathogenesis. In this study, by comparative genome analysis, genome-wide peptidase diversities among plant-parasitic nematodes are estimated. We find that genes encoding cysteine peptidases in family C13 (legumain) are significantly abundant in pine wood nematodes Bursaphelenchus genomes, compared to those in other plant-parasitic nematodes. By phylogenetic analysis, a clade of B. xylophilus-specific legumain is identified. RT-qPCR detection shows that these genes are highly expressed at early stage during the nematode infection process. Utilizing transgene technology, cDNAs of three species-specific legumain were introduced into the Arabidopsis γvpe mutant. Functional complementation assay shows that these B. xylophilus legumains can fully complement the activity of Arabidopsis γVPE to mediate plant cell death triggered by the fungal toxin FB1. Secretory activities of these legumains are experimentally validated. By comparative transcriptome analysis, genes involved in plant cell death mediated by legumains are identified, which enrich in GO terms related to ubiquitin protein transferase activity in category molecular function, and response to stimuli in category biological process. Our results suggest that B. xylophilu-specific legumains have potential as effectors to be involved in nematode-plant interaction and can be related to host cell death.

Telomere-to-telomere genome assembly of asparaginase-producing Trichoderma simmonsii.

BACKGROUND: Trichoderma is a genus of fungi in the family Hypocreaceae and includes species known to produce enzymes with commercial use. They are largely found in soil and terrestrial plants. Recently, Trichoderma simmonsii isolated from decaying bark and decorticated wood was newly identified in the Harzianum clade of Trichoderma. Due to a wide range of applications in agriculture and other industries, genomes of at least 12 Trichoderma spp. have been studied. Moreover, antifungal and enzymatic activities have been extensively characterized in Trichoderma spp. However, the genomic information and bioactivities of T. simmonsii from a particular marine-derived isolate remain largely unknown. While we screened for asparaginase-producing fungi, we observed that T. simmonsii GH-Sj1 strain isolated from edible kelp produced asparaginase. In this study, we report a draft genome of T. simmonsii GH-Sj1 using Illumina and Oxford Nanopore technologies. Furthermore, to facilitate biotechnological applications of this species, RNA-sequencing was performed to elucidate the transcriptional profile of T. simmonsii GH-Sj1 in response to asparaginase-rich conditions. RESULTS: We generated ~ 14 Gb of sequencing data assembled in a ~ 40 Mb genome. The T. simmonsii GH-Sj1 genome consisted of seven telomere-to-telomere scaffolds with no sequencing gaps, where the N50 length was 6.4 Mb. The total number of protein-coding genes was 13,120, constituting ~ 99% of the genome. The genome harbored 176 tRNAs, which encode a full set of 20 amino acids. In addition, it had an rRNA repeat region consisting of seven repeats of the 18S-ITS1-5.8S-ITS2-26S cluster. The T. simmonsii genome also harbored 7 putative asparaginase-encoding genes with potential medical applications. Using RNA-sequencing analysis, we found that 3 genes among the 7 putative genes were significantly upregulated under asparaginase-rich conditions. CONCLUSIONS: The genome and transcriptome of T. simmonsii GH-Sj1 established in the current work represent valuable resources for future comparative studies on fungal genomes and asparaginase production.

Genome sequencing and comparative genome analysis of 6 hypervirulent Klebsiella pneumoniae strains isolated in China.

Hypervirulent Klebsiella pneumoniae (hvKP) has been increasingly reported over the past three decades and causes severe infections. To increase our understanding of hvKP at the genome level, genome sequencing and comparative genome analysis were performed on 6 hvKPs. The whole genome DNA from 6 hvKPs with different capsular serotypes isolated in China was extracted. The genome sequencing and assembly results showed the genome size of the six hvKPs and GC content. Comparative analyses of the genomes revealed the gene homology and genome rearrangement in the 6 hvKPs compared with Klebsiella pneumonia NTUH-K2044. The phylogenetic tree based on full-genome SNPs of the 7 hvKPs showed that NTUH-K2044 formed a single clade, showing distant evolutionary distances with the other six strains, and the non-K1 hvKP strains had a relatively closer phylogenetic relationship. BLAST comparison analysis found that some selected virulence genes had different degrees of deletion in the non-K1 hvKPs. SNP-based virulence gene mutation analysis showed that some virulence genes had different degrees of SNP mutations. The whole-genome sequencing and comparative genome analysis of six hvKP strains with NTUH-K2044 provide us with a basic understanding of the genome composition, genetic polymorphism, evolution and virulence genes of hvKP and a basis for further research on these genes and the pathogenesis of hvKP.

Comparative analysis using the draft genome sequence of California poppy (Eschscholzia californica) for exploring the candidate genes involved in benzylisoquinoline alkaloid biosynthesis.

Genome characterization of California poppy (Eschscholzia californica cv. "Hitoezaki"), which produces pharmaceutically important benzylisoquinoline alkaloids (BIAs), was carried out using the draft genome sequence. The numbers of tRNA and rRNA genes were close to those of the other plant species tested, whereas the frequency of repetitive sequences was distinct from those species. Comparison of the predicted genes with those of Amborella trichopoda, Nelumbo nucifera, Solanum lycopersicum, and Arabidopsis thaliana, and analyses of gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway indicated that the enzyme genes involved in BIA biosynthesis were highly enriched in the California poppy genome. Further comparative analysis using the genome information of Papaver somniferum and Aquilegia coerulea, both BIA-producing plants, revealed that many genes encoding BIA biosynthetic enzymes, transcription factors, transporters, and candidate proteins, possibly related to BIA biosynthesis, were specifically distributed in these plant species.

Biodegradation of phenol by a highly tolerant strain Rhodococcus ruber C1: Biochemical characterization and comparative genome analysis.

A novel phenol-degrading strain was isolated and identified as Rhodococcus ruber C1. The degradation analysis shows that 1806 mg/L of phenol can be completely degraded by strain C1 within 38 h, and the maximum specific growth rate (µmax=1.527 h-1) and maximum specific phenol degradation rate (qmax=3.674 h-1) indicate its excellent phenol metabolism capability. More importantly, phenol can be degraded by strain C1 in the temperature range of 20-45 °C within 72 h, and with longer degradation time, phenol can be completely degraded even at 10, 15 and 50 °C. The whole genome of strain C1 was sequenced, and a comparative genome analysis of strain C1 with 36 other genomes of Rhodococcus was performed. A remarkable gene family expansion occurred during the evolution of Rhodococcus, and a comprehensive evolutionary picture of Rhodococcus at genomic level was presented. Moreover, the copy number of genes involved in phenol metabolism was compared among genus Rhodococcus, and the results demonstrate high phenol degradation capability of strain C1 at genomic level. These findings suggest that Rhodococcus ruber C1 is a bacterium capable of degrading phenol efficiently in the temperature range of 10-50 °C.

Pathogenic and antimicrobial resistance genes in Streptococcus oralis strains revealed by comparative genome analysis.

Streptococcus oralis is an early colonizer bacterium in dental plaques and is considered a potential pathogen of infective endocarditis (IE) disease. In this study, we built a complete genome map of Streptococcus oralis strain SOT, Streptococcus oralis strain SOD and Streptococcus infantis strain SO and performed comparative genomic analysis among these three strains. The results showed that there are five genomic islands (GIs) in strain SOT and one CRISPR in strain SOD. Each genome harbors various pathogenic genes related to diseases and drug resistance, while the antibiotic resistance genes in strains SOT and SOD were quite similar but different from those in strain SO. In addition, we identified 17 main virulence factors and capsule-related genes in three strains. These results suggest the pathogenic potential of Streptococcus strains, which lay a foundation for the prevention and treatment of a Streptococcus oralis infection.