Efficient production of retinol in Yarrowia lipolytica by increasing stability using antioxidant and detergent extraction.

The demand for bio-based retinol (vitamin A) is currently increasing, however its instability represents a major bottleneck in microbial production. Here, we developed an efficient method to selectively produce retinol in Yarrowia lipolytica. The ß-carotene 15,15'-dioxygenase (BCO) cleaves ß-carotene into retinal, which is reduced to retinol by retinol dehydrogenase (RDH). Therefore, to produce retinol, we first generated ß-carotene-producing strain based on a high-lipid-producer via overexpressing genes including heterologous ß-carotene biosynthetic genes, GGS1F43I mutant of endogenous geranylgeranyl pyrophosphate synthase isolated by directed evolution, and FAD1 encoding flavin adenine dinucleotide synthetase, while deleting several genes previously known to be beneficial for carotenoid production. To produce retinol, 11 copies of BCO gene from marine bacterium 66A03 (Mb.Blh) were integrated into the rDNA sites of the ß-carotene overproducer. The resulting strain produced more retinol than retinal, suggesting strong endogenous promiscuous RDH activity in Y. lipolytica. The introduction of Mb.Blh led to a considerable reduction in ß-carotene level, but less than 5% of the consumed ß-carotene could be detected in the form of retinal or retinol, implying severe degradation of the produced retinoids. However, addition of the antioxidant butylated hydroxytoluene (BHT) led to a >20-fold increase in retinol production, suggesting oxidative damage is the main cause of intracellular retinol degradation. Overexpression of GSH2 encoding glutathione synthetase further improved retinol production. Raman imaging revealed co-localization of retinol with lipid droplets, and extraction of retinol using Tween 80 was effective in improving retinol production. By combining BHT treatment and extraction using Tween 80, the final strain CJ2104 produced 4.86 g/L retinol and 0.26 g/L retinal in fed-batch fermentation in a 5-L bioreactor, which is the highest retinol production titer ever reported. This study demonstrates that Y. lipolytica is a suitable host for the industrial production of bio-based retinol.

Epidemiologic Linkage of COVID-19 Outbreaks at Two University-affiliated Hospitals in the Seoul Metropolitan Area in March 2020.

BACKGROUND: Coronavirus disease 2019 (COVID-19) outbreaks emerged at two university-affiliated hospitals in Seoul (hospital A) and Uijeongbu City (hospital S) in the metropolitan Seoul area in March 2020. The aim of this study was to investigate epidemiological links between the outbreaks using whole genome sequencing (WGS) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). METHODS: Fifteen patients were enrolled in the study, including four non-outbreak (A1-A4) and three outbreak cases (A5-A7) in hospital A and eight cases (S1-S8) in hospital S. Patients' hospital stays, COVID-19 symptoms, and transfer history were reviewed. RNA samples were submitted for WGS and genome-wide single nucleotide variants and phylogenetic relationships were analyzed. RESULTS: The index patient (A5) in hospital A was transferred from hospital S on 26 March. Patients A6 and A7 were the family caregiver and sister, respectively, of the patient who shared a room with A5 for 4 days. Prior to transfer, A5 was at the next bed to S8 in the emergency room on 25 March. Patient S6, a professional caregiver, took care of the patient in the room next to S8's room for 5 days until 22 March and then S5 for another 3 days. WGS revealed that SARS-CoV-2 in A2, A3, and A4 belong to clades V/B.2, S/A, and G/B.1, respectively, whereas that of A5-A7 and S1-S5 are of the V/B.2.1 clade and closely clustered. In particular, SARS-CoV-2 in patients A5 and S5 showed perfect identity. CONCLUSION: WGS is a useful tool to understand epidemiology of SARS-CoV-2. It is the first study to elucidate the role of patient transfer and caregivers as links of nosocomial outbreaks of COVID-19 in multiple hospitals.

ContEst16S: an algorithm that identifies contaminated prokaryotic genomes using 16S RNA gene sequences.

Thanks to the recent advancement of DNA sequencing technology, the cost and time of prokaryotic genome sequencing have been dramatically decreased. It has repeatedly been reported that genome sequencing using high-throughput next-generation sequencing is prone to contaminations due to its high depth of sequencing coverage. Although a few bioinformatics tools are available to detect potential contaminations, these have inherited limitations as they only use protein-coding genes. Here we introduce a new algorithm, called ContEst16S, to detect potential contaminations using 16S rRNA genes from genome assemblies. We screened 69 745 prokaryotic genomes from the NCBI Assembly Database using ContEst16S and found that 594 were contaminated by bacteria, human and plants. Of the predicted contaminated genomes, 8 % were not predicted by the existing protein-coding gene-based tool, implying that both methods can be complementary in the detection of contaminations. A web-based service of the algorithm is available at www.ezbiocloud.net/tools/contest16s.

Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies.

The recent advent of DNA sequencing technologies facilitates the use of genome sequencing data that provide means for more informative and precise classification and identification of members of the Bacteria and Archaea. Because the current species definition is based on the comparison of genome sequences between type and other strains in a given species, building a genome database with correct taxonomic information is of paramount need to enhance our efforts in exploring prokaryotic diversity and discovering novel species as well as for routine identifications. Here we introduce an integrated database, called EzBioCloud, that holds the taxonomic hierarchy of the Bacteria and Archaea, which is represented by quality-controlled 16S rRNA gene and genome sequences. Whole-genome assemblies in the NCBI Assembly Database were screened for low quality and subjected to a composite identification bioinformatics pipeline that employs gene-based searches followed by the calculation of average nucleotide identity. As a result, the database is made of 61 700 species/phylotypes, including 13 132 with validly published names, and 62 362 whole-genome assemblies that were identified taxonomically at the genus, species and subspecies levels. Genomic properties, such as genome size and DNA G+C content, and the occurrence in human microbiome data were calculated for each genus or higher taxa. This united database of taxonomy, 16S rRNA gene and genome sequences, with accompanying bioinformatics tools, should accelerate genome-based classification and identification of members of the Bacteria and Archaea. The database and related search tools are available at www.ezbiocloud.net/.

A large-scale evaluation of algorithms to calculate average nucleotide identity.

Average nucleotide identity (ANI) is a category of computational analysis that can be used to define species boundaries of Archaea and Bacteria. Calculating ANI usually involves the fragmentation of genome sequences, followed by nucleotide sequence search, alignment, and identity calculation. The original algorithm to calculate ANI used the BLAST program as its search engine. An improved ANI algorithm, called OrthoANI, was developed to accommodate the concept of orthology. Here, we compared four algorithms to compute ANI, namely ANIb (ANI algorithm using BLAST), ANIm (ANI using MUMmer), OrthoANIb (OrthoANI using BLAST) and OrthoANIu (OrthoANI using USEARCH) using >100,000 pairs of genomes with various genome sizes. By comparing values to the ANIb that is considered a standard, OrthoANIb and OrthoANIu exhibited good correlation in the whole range of ANI values. ANIm showed poor correlation for ANI of <90%. ANIm and OrthoANIu runs faster than ANIb by an order of magnitude. When genomes that are larger than 7 Mbp were analysed, the run-times of ANIm and OrthoANIu were shorter than that of ANIb by 53- and 22-fold, respectively. In conclusion, ANI calculation can be greatly sped up by the OrthoANIu method without losing accuracy. A web-service that can be used to calculate OrthoANIu between a pair of genome sequences is available at http://www.ezbiocloud.net/tools/ani . For large-scale calculation and integration in bioinformatics pipelines, a standalone JAVA program is available for download at http://www.ezbiocloud.net/tools/orthoaniu .

Mucosal and salivary microbiota associated with recurrent aphthous stomatitis.

BACKGROUND: Recurrent aphthous stomatitis (RAS) is a common oral mucosal disorder of unclear etiopathogenesis. Although recent studies of the oral microbiota by high-throughput sequencing of 16S rRNA genes have suggested that imbalances in the oral microbiota may contribute to the etiopathogenesis of RAS, no specific bacterial species associated with RAS have been identified. The present study aimed to characterize the microbiota in the oral mucosa and saliva of RAS patients in comparison with control subjects at the species level. RESULTS: The bacterial communities of the oral mucosa and saliva from RAS patients with active lesions (RAS, n = 18 for mucosa and n = 8 for saliva) and control subjects (n = 18 for mucosa and n = 7 for saliva) were analyzed by pyrosequencing of the 16S rRNA genes. There were no significant differences in the alpha diversity between the controls and the RAS, but the mucosal microbiota of the RAS patients showed increased inter-subject variability. A comparison of the relative abundance of each taxon revealed decreases in the members of healthy core microbiota but increases of rare species in the mucosal and salivary microbiota of RAS patients. Particularly, decreased Streptococcus salivarius and increased Acinetobacter johnsonii in the mucosa were associated with RAS risk. A dysbiosis index, which was developed using the relative abundance of A. johnsonii and S. salivarius and the regression coefficients, correctly predicted 83 % of the total cases for the absence or presence of RAS. Interestingly, A. johnsonii substantially inhibited the proliferation of gingival epithelial cells and showed greater cytotoxicity against the gingival epithelial cells than S. salivarius. CONCLUSION: RAS is associated with dysbiosis of the mucosal and salivary microbiota, and two species associated with RAS have been identified. This knowledge may provide a diagnostic tool and new targets for therapeutics for RAS.

Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species.

Despite recent advances in commercially optimized identification systems, bacterial identification remains a challenging task in many routine microbiological laboratories, especially in situations where taxonomically novel isolates are involved. The 16S rRNA gene has been used extensively for this task when coupled with a well-curated database, such as EzTaxon, containing sequences of type strains of prokaryotic species with validly published names. Although the EzTaxon database has been widely used for routine identification of prokaryotic isolates, sequences from uncultured prokaryotes have not been considered. Here, the next generation database, named EzTaxon-e, is formally introduced. This new database covers not only species within the formal nomenclatural system but also phylotypes that may represent species in nature. In addition to an identification function based on Basic Local Alignment Search Tool (blast) searches and pairwise global sequence alignments, a new objective method of assessing the degree of completeness in sequencing is proposed. All sequences that are held in the EzTaxon-e database have been subjected to phylogenetic analysis and this has resulted in a complete hierarchical classification system. It is concluded that the EzTaxon-e database provides a useful taxonomic backbone for the identification of cultured and uncultured prokaryotes and offers a valuable means of communication among microbiologists who routinely encounter taxonomically novel isolates. The database and its analytical functions can be found at http://eztaxon-e.ezbiocloud.net/.

Impact of enrofloxacin on the human intestinal microbiota revealed by comparative molecular analysis.

The indigenous human intestinal microbiota could be disrupted by residues of antibiotics in foods as well as therapeutically administered antibiotics to humans. These disruptions may lead to adverse health outcomes. To observe the possible impact of residues of antibiotics at concentrations below therapeutic levels on human intestinal microbiota, we performed studies using in vitro cultures of fecal suspensions from three individuals with 10 different concentrations (0, 0.1, 0.5, 1, 5, 10, 15, 25, 50 and 150 µg/ml) of the fluoroquinolone, enrofloxacin. The bacterial communities of the control and enrofloxacin dosed fecal samples were analyzed by denaturing gradient gel electrophoresis (DGGE) and pyrosequencing. In addition, changes of functional gene expression were analyzed by a pyrosequencing-based random whole-community mRNA sequencing method. Although each individual had a unique microbial composition, the communities of all individuals were affected by enrofloxacin. The proportions of two phyla, namely, Bacteroidetes and Proteobacteria, were significantly reduced with increasing concentrations of enrofloxacin exposure, while the proportion of Firmicutes increased. Principal Coordinate Analysis (PCoA) using the Fast UniFrac indicated that the community structures of intestinal microbiota were shifted by enrofloxacin. Most of the mRNA transcripts and the anti-microbial drug resistance genes increased with increasing concentrations of enrofloxacin. 16S rRNA gene pyrosequencing of control and enrofloxacin treated fecal suspensions provided valuable information of affected bacterial taxa down to the species level, and the community transcriptomic analyses using mRNA revealed the functional gene expression responses of the changed bacterial communities by enrofloxacin.

Improvement of image quality using amplitude-based respiratory gating in PET-computed tomography scanning.

OBJECTIVES: This study sought to provide data supporting the expanded clinical use of respiratory gating by assessing the diagnostic accuracy of breathing motion correction using amplitude-based respiratory gating. METHODS: A respiratory movement tracking device was attached to a PET-computed tomography scanner, and images were obtained in respiratory gating mode using a motion phantom that was capable of sensing vertical motion. Specifically, after setting amplitude changes and intervals according to the movement cycle using a total of nine combinations of three waveforms and three amplitude ranges, respiratory motion-corrected images were reconstructed using the filtered back projection method. After defining areas of interest in the acquired images in the same image planes, statistical analyses were performed to compare differences in standardized uptake value (SUV), lesion volume, full width at half maximum (FWHM), signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). RESULTS: SUVmax increased by 89.9%, and lesion volume decreased by 27.9%. Full width at half maximum decreased by 53.9%, signal-to-noise ratio increased by 11% and contrast-to-noise ratio increased by 16.3%. Optimal results were obtained when using a rest waveform and 35% duty cycle, in which the change in amplitude in the respiratory phase signal was low, and a constant level of long breaths was maintained. CONCLUSIONS: These results demonstrate that respiratory-gated PET-CT imaging can be used to accurately correct for SUV changes and image distortion caused by respiratory motion, thereby providing excellent imaging information and quality.

Complete Multipartite Genome Sequence of the Cupriavidus basilensis Type Strain, a 2,6-Dichlorophenol-Degrading Bacterium.

We report the complete 8.94-Mb genome sequence of the type strain of Cupriavidus basilensis (DSM 11853T = CCUG 49340T = RK1T), formed by two chromosomes and six putative plasmids, which offers insights into its chloroaromatic-biodegrading capabilities.

Emotional well-being and gut microbiome profiles by enterotype.

With increasing attention being paid to improving emotional well-being, recent evidence points to gut microbiota as a key player in regulating mental and physical health via bidirectional communication between the brain and gut. Here, we examine the association between emotional well-being and gut microbiome profiles (i.e., gut microbiome composition, diversity, and the moderating role of the enterotypes) among healthy Korean adults (n = 83, mean age = 48.9, SD = 13.2). The research was performed using high-throughput 16S rRNA gene sequencing to obtain gut microbiome profiles, as well as a self-report survey that included the Positive Affect Negative Affect Schedule (PANAS). The cluster-based analysis identified two enterotypes dominated by the genera Bacteroides (n = 49) and Prevotella (n = 34). Generalized linear regression analysis reveals significant associations between positive emotion and gut microbiome diversity (Shannon Index) among participants in the Prevotella dominant group, whereas no such relationship emerged among participants in the Bacteroides group. Moreover, a novel genus from the family Lachnospiraceae is associated with emotional well-being scores, both positive and negative. Together, the current findings highlight the enterotype-specific links between the gut microbiota community and emotion in healthy adults and suggest the possible roles of the gut microbiome in promoting mental health.

Improved Metagenomic Taxonomic Profiling Using a Curated Core Gene-Based Bacterial Database Reveals Unrecognized Species in the Genus Streptococcus.

Shotgun metagenomics is of great importance in order to understand the composition of the microbial community associated with a sample and the potential impact it may exert on its host. For clinical metagenomics, one of the initial challenges is the accurate identification of a pathogen of interest and ability to single out that pathogen within a complex community of microorganisms. However, in absence of an accurate identification of those microorganisms, any kind of conclusion or diagnosis based on misidentification may lead to erroneous conclusions, especially when comparing distinct groups of individuals. When comparing a shotgun metagenomic sample against a reference genome sequence database, the classification itself is dependent on the contents of the database. Focusing on the genus Streptococcus, we built four synthetic metagenomic samples and demonstrated that shotgun taxonomic profiling using the bacterial core genes as the reference database performed better in both taxonomic profiling and relative abundance prediction than that based on the marker gene reference database included in MetaPhlAn2. Additionally, by classifying sputum samples of patients suffering from chronic obstructive pulmonary disease, we showed that adding genomes of genomospecies to a reference database offers higher taxonomic resolution for taxonomic profiling. Finally, we show how our genomospecies database is able to identify correctly a clinical stool sample from a patient with a streptococcal infection, proving that genomospecies provide better taxonomic coverage for metagenomic analyses.

Comparative Genomic Analysis of the 2016 Vibrio cholerae Outbreak in South Korea.

In August 2016, South Korea experienced a cholera outbreak that caused acute watery diarrhea in three patients. This outbreak was the first time in 15 years that an outbreak was not linked to an overseas source. To identify the cause and to study the epidemiological implications of this outbreak, we sequenced the whole genome of Vibrio cholerae isolates; three from each patient and one from a seawater sample. Herein we present comparative genomic data which reveals that the genome sequences of these four isolates are very similar. Interestingly, these isolates form a monophyletic clade with V. cholerae strains that caused an outbreak in the Philippines in 2011. The V. cholerae strains responsible for the Korean and Philippines outbreaks have almost identical genomes in which two unique genomic islands are shared, and they both lack SXT elements. Furthermore, we confirm that seawater is the likely source of this outbreak, which suggests the necessity for future routine surveillance of South Korea's seashore.

UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction.

Genome-based phylogeny plays a central role in the future taxonomy and phylogenetics of Bacteria and Archaea by replacing 16S rRNA gene phylogeny. The concatenated core gene alignments are frequently used for such a purpose. The bacterial core genes are defined as single-copy, homologous genes that are present in most of the known bacterial species. There have been several studies describing such a gene set, but the number of species considered was rather small. Here we present the up-to-date bacterial core gene set, named UBCG, and software suites to accommodate necessary steps to generate and evaluate phylogenetic trees. The method was successfully used to infer phylogenomic relationship of Escherichia and related taxa and can be used for the set of genomes at any taxonomic ranks of Bacteria. The UBCG pipeline and file viewer are freely available at https://www.ezbiocloud.net/tools/ubcg and https://www.ezbiocloud.net/tools/ubcg_viewer , respectively.

A single gene of a commensal microbe affects host susceptibility to enteric infection.

Indigenous microbes inside the host intestine maintain a complex self-regulating community. The mechanisms by which gut microbes interact with intestinal pathogens remain largely unknown. Here we identify a commensal Escherichia coli strain whose expansion predisposes mice to infection by Vibrio cholerae, a human pathogen. We refer to this strain as 'atypical' E. coli (atEc) because of its inability to ferment lactose. The atEc strain is resistant to reactive oxygen species (ROS) and proliferates extensively in antibiotic-treated adult mice. V. cholerae infection is more severe in neonatal mice transplanted with atEc compared with those transplanted with a typical E. coli strain. Intestinal ROS levels are decreased in atEc-transplanted mice, favouring proliferation of ROS-sensitive V. cholerae. An atEc mutant defective in ROS degradation fails to facilitate V. cholerae infection when transplanted, suggesting that host infection susceptibility can be regulated by a single gene product of one particular commensal species.

Complete Genome Sequence of Middle East Respiratory Syndrome Coronavirus KOR/KNIH/002_05_2015, Isolated in South Korea.

The full genome sequence of a Middle East respiratory syndrome coronavirus (MERS-CoV) was identified from cultured and isolated in Vero cells. The viral genome sequence has high similarity to 53 human MERS-CoVs, ranging from 99.5% to 99.8% at the nucleotide level.

Comparative genomics of Neisseria weaveri clarifies the taxonomy of this species and identifies genetic determinants that may be associated with virulence.

A group of bacterial strains formerly known as CDC group M-5 are opportunistic pathogens to humans. In 1993, a name, Neisseria weaveri, was proposed by two independent studies to harbor CDC group M-5 strains, namely N. weaveri Holmes et al. 1993 and N. weaveri Andersen et al. 1993, with two different 'type' strains. However, no study has been conducted on to the relatedness of the two 'type' strains, although the close relationship of the two taxa has long been accepted unofficially. Formally, the status of the name N. weaveri Andersen et al. 1993 is illegitimate because it is a later homonym of N. weaveri Holmes et al., 1993; but the name of the strain is still validly published. In this study, we attempt to resolve the confusion caused by the apparent duplication of the species N. weaveri (with different type strains) using whole genome shotgun sequencing. We also sought to gain insight into the genetic characteristics of N. weaveri by conducting comparative genomics. On the basis of genomic similarities revealed through a comparative genomic study, we propose that N. weaveri Andersen et al. 1993 should be re-classified as a later heterotypic synonym of N. weaveri Holmes et al., 1993.

Fungal community associated with genetically modified poplar during metal phytoremediation.

Due to the increasing demand for phytoremediation, many transgenic poplars have been developed to enhance the bioremediation of heavy metals. However, structural changes to indigenous fungal communities by genetically modified organisms (GMO) presents a major ecological issue, due to the important role of fungi for plant growth in natural environments. To evaluate the effect of GM plant use on environmental fungal soil communities, extensive sequencing-based community analysis was conducted, while controlling the influence of plant clonality, plant age, soil condition, and harvesting season. The rhizosphere soils of GM and wild type (WT) poplars at a range of growth stages were sampled together with unplanted, contaminated soil, and the fungal community structures were investigated by pyrosequencing the D1/D2 region of the 28S rRNA gene. The results show that the overall structure of the rhizosphere fungal community was not significantly influenced by GM poplars. However, the presence of GM specific taxa, and faster rate of community change during poplar growth, appeared to be characteristic of the GM plant-induced effects on soil-born fungal communities. The results of this study provide additional information about the potential effects of GM poplar trees aged 1.5-3 years, on the soil fungal community.

Comparative approach to capture bacterial diversity of coastal waters.

Despite the revolutionary advancements in DNA sequencing technology and cultivation techniques, few studies have been done to directly compare these methods. In this study, a 16S rRNA gene-based, integrative approach combining culture-independent techniques with culture-dependent methods was taken to investigate the bacterial community structure of coastal seawater collected from the Yellow Sea, Korea. For culture-independent studies, we used the latest model pyrosequencer, Roche/454 Genome Sequencer FLX Titanium. Pyrosequencing captured a total of 52 phyla including 27 candidate divisions from the water column, whereas the traditional cloning approach captured only 15 phyla including 2 candidate divisions. In addition, of 878 genera retrieved, 92.1 % of the sequences were unique to pyrosequencing. For culture-dependent analysis, plate culturing, plate washing, enrichment, and high-throughput culturing (HTC) methods were applied. Phylogenetic analysis showed that the plate-washing clones formed a cluster devoid of any previously cultured representatives within the family Rhodobacteraceae. One HTC isolate (SF293) fell into the OM182 clade, which was not recovered by other culturing methods described here. By directly comparing the sequences obtained from cultures with those from culture-independent work, we found that only 33% of the culture sequences were identical to those from clone libraries and pyrosequences. This study presents a detailed comparison of common molecular and cultivation techniques available in microbial ecology. As different methods yielded different coverage, we suggest choosing the approach after carefully examining the scientific questions being asked.