De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer

Probiogenomics of Leuconostoc Mesenteroides Strains F-21 and F-22 Isolated from Human Breast Milk Reveal Beneficial Properties

Bacteria of the Leuconostoc genus are Gram-positive bacteria that are commonly found in raw milk and persist in fermented dairy products and plant food. Studies have already explored the probiotic potential of L. mesenteroides, but not from a probiogenomic perspective, which aims to explore the molecular features responsible for their phenotypes. In the present work, probiogenomic approaches were applied in strains F-21 and F-22 of L. mesenteroides isolated from human milk to assess their biosafety at the molecular level and to correlate molecular features with their potential probiotic characteristics. The complete genome of strain F-22 is 1.99 Mb and presents one plasmid, while the draft genome of strain F-21 is 1.89 Mb and presents four plasmids. A high percentage of average nucleotide identity among other genomes of L. mesenteroides (≥ 96%) corroborated the previous taxonomic classification of these isolates. Genomic regions that influence the probiotic properties were identified and annotated. Both strains exhibited wide genome plasticity, cell adhesion ability, proteolytic activity, proinflammatory and immunomodulation capacity through interaction with TLR-NF-κB and TLR-MAPK pathway components, and no antimicrobial resistance, denoting their potential to be candidate probiotics. Further, the strains showed bacteriocin production potential and the presence of acid, thermal, osmotic, and bile salt resistance genes, indicating their ability to survive under gastrointestinal stress. Taken together, our results suggest that L. mesenteroides F-21 and F-22 are promising candidates for probiotics in the food and pharmaceutical industries.

Gut bacterium Intestinimonas butyriciproducens improves host metabolic health: evidence from cohort and animal intervention studies

BACKGROUND

The human gut microbiome strongly influences host metabolism by fermenting dietary components into metabolites that signal to the host. Our previous work has shown that Intestinimonas butyriciproducens is a prevalent commensal bacterium with the unique ability to convert dietary fructoselysine to butyrate, a well-known signaling molecule with proven health benefits. Dietary fructoselysine is an abundant Amadori product formed in foods during thermal treatment and is part of foods rich in dietary advanced glycation end products which have been associated with cardiometabolic disease. It is therefore of interest to investigate the causal role of this bacterium and fructoselysine metabolism in metabolic disorders.

RESULTS

We assessed associations of I. butyriciproducens with metabolic risk biomarkers at both strain and functional levels using a human cohort characterized by fecal metagenomic analysis. We observed that the level of the bacterial strain as well as fructoselysine fermentation genes were negatively associated with BMI, triglycerides, HbA1c, and fasting insulin levels. We also investigated the fructoselysine degradation capacity within the Intestinimonas genus using a culture-dependent approach and found that I. butyriciproducens is a key player in the butyrogenic fructoselysine metabolism in the gut. To investigate the function of I. butyriciproducens in host metabolism, we used the diet-induced obesity mouse model to mimic the human metabolic syndrome. Oral supplementation with I. butyriciproducens counteracted body weight gain, hyperglycemia, and adiposity. In addition, within the inguinal white adipose tissue, bacterial administration reduced inflammation and promoted pathways involved in browning and insulin signaling. The observed effects may be partly attributable to the formation of the short-chain fatty acids butyrate from dietary fructoselysine, as butyrate plasma and cecal levels were significantly increased by the bacterial strain, thereby contributing to the systemic effects of the bacterial treatment.

CONCLUSIONS

I. butyriciproducens ameliorates host metabolism in the context of obesity and may therefore be a good candidate for new microbiota-therapeutic approaches to prevent or treat metabolic diseases. Video Abstract.

Genomic Analyses of Methicillin-Susceptible and Methicillin-Resistant Staphylococcus pseudintermedius Strains Involved in Canine Infections: A Comprehensive Genotypic Characterization

Staphylococcus pseudintermedius is frequently associated with several bacterial infections in dogs, highlighting a One Health concern due to the zoonotic potential. Given the clinical significance of this pathogen, we performed comprehensive genomic analyses of 28 S. pseudintermedius strains isolated from canine infections throughout whole-genome sequencing using Illumina HiSeq, and compared the genetic features between S. pseudintermedius methicillin-resistant (MRSP) and methicillin-susceptible (MSSP) strains. Our analyses determined that MRSP genomes are larger than MSSP strains, with significant changes in antimicrobial resistance genes and virulent markers, suggesting differences in the pathogenicity of MRSP and MSSP strains. In addition, the pangenome analysis of S. pseudintermedius from canine and human origins identified core and accessory genomes with 1847 and 3037 genes, respectively, which indicates that most of the S. pseudintermedius genome is highly variable. Furthermore, phylogenomic analysis clearly separated MRSP from MSSP strains, despite their infection sites, showing phylogenetic differences according to methicillin susceptibility. Altogether our findings underscore the importance of studying the evolutionary dynamics of S. pseudintermedius, which is crucial for the development of effective prevention and control strategies of resistant S. pseudintermedius infections.

Comparative genetic and pathogenic approaches of Escherichia coli isolated simultaneously from pyometra and urine of bitches

Escherichia coli (E. coli) are widely related to pyometra and cystitis in dogs, and these infections can occur simultaneously. The goal of this study was to determine genetic and pathogenic insights of 14 E. coli isolated simultaneously from pyometra content and bladder urine of seven bitches. To achieve this, in silico and in vitro comparative analyses were conducted. Whole-genome comparisons demonstrated that E. coli isolated from pyometra and urine of the same animal were predominantly genetic extraintestinal E. coli clones belonging to the same Sequence Type and phylogroup. The E. coli clones identified in this study included ST372, ST457, ST12, ST127, ST646, and ST961. Five isolates (35.7%) belonged to the ST12 complex. Except for two E. coli, all other isolates belonged to the B2 Clermont phylogroup. Interestingly, some genomes of E. coli from urine carried more virulence genes than those E. coli from pyometra. Both pyometra and urine E. coli isolates demonstrated a strong affinity for adhering to HeLa and T24 cells, with a low affinity for invading them. However, certain isolates from urine exhibited a greater tendency to adhere to T24 cells in qualitative and quantitative assays compared to isolates from pyometra. In conclusion, this study revealed the high genomic similarity between pyometra and urine E. coli isolates, as well as the virulent capacity of both to colonize endometrial and urothelial cells. The findings of this study underscore the importance of concurrently managing both infections clinically and could potentially contribute to future resources for the prevention of cystitis and pyometra.

Environmental Biofilms from an Urban Community in Salvador, Brazil, Shelter Previously Uncharacterized Saprophytic Leptospira

Biofilms are complex microecosystems with valuable ecological roles that can shelter a variety of microorganisms. Spirochetes from the genus Leptospira have been observed to form biofilms in vitro, in rural environments, and in the kidneys of reservoir rats. The genus Leptospira is composed of pathogenic and non-pathogenic species, and the description of new species is ongoing due to the advent of whole genome sequencing. Leptospires have increasingly been isolated from water and soil samples. To investigate the presence of Leptospira in environmental biofilms, we collected three distinct samples of biofilms formed in an urban setting with poor sanitation: Pau da Lima, in Salvador, Bahia, Brazil. All biofilm samples were negative for the presence of pathogenic leptospires via conventional PCR, but cultures containing saprophytic Leptospira were identified. Whole genomes were generated and analyzed for twenty isolates obtained from these biofilms. For species identification, we used digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) analysis. The obtained isolates were classified into seven presumptive species from the saprophytic S1 clade. ANI and dDDH analysis suggest that three of those seven species were new. Classical phenotypic tests confirmed the novel isolated bacteria as saprophytic Leptospira. The isolates presented typical morphology and ultrastructure according to scanning electron microscopy and formed biofilms under in vitro conditions. Our data indicate that a diversity of saprophytic Leptospira species survive in the Brazilian poorly sanitized urban environment, in a biofilm lifestyle. We believe our results contribute to a better understanding of Leptospira biology and ecology, considering biofilms as natural environmental reservoirs for leptospires.

A Comprehensive Bioinformatics Resource Guide for Genome-Based Antimicrobial Resistance Studies

The use of high-throughput sequencing technologies and bioinformatic tools has greatly transformed microbial genome research. With the help of sophisticated computational tools, it has become easier to perform whole genome assembly, identify and compare different species based on their genomes, and predict the presence of genes responsible for proteins, antimicrobial resistance, and toxins. These bioinformatics resources are likely to continuously improve in quality, become more user-friendly to analyze the multiple genomic data, efficient in generating information and translating it into meaningful knowledge, and enhance our understanding of the genetic mechanism of AMR. In this manuscript, we provide an essential guide for selecting the popular resources for microbial research, such as genome assembly and annotation, antibiotic resistance gene profiling, identification of virulence factors, and drug interaction studies. In addition, we discuss the best practices in computer-oriented microbial genome research, emerging trends in microbial genomic data analysis, integration of multi-omics data, the appropriate use of machine-learning algorithms, and open-source bioinformatics resources for genome data analytics.

Genome analysis of a coral-associated bacterial consortium highlights complementary hydrocarbon degradation ability and other beneficial mechanisms for the host

Here we report the oil degradation genetic potential of six oil-degrading bacteria (ODB), previously used as a bioremediation consortium, isolated from the hydrocoral Millepora alcicornis and seawater. The strains were identified as Halomonas sp. (LC_1), Cobetia sp. (LC_6), Pseudoalteromonas shioyasakiensis (LC_2), Halopseudomonas aestusnigri (LC_3), Shewanella algae (LC_4), and Brucella intermedia (LC_5). The taxonomic identification differed from that of the original paper when we used whole genome gene markers instead of just 16S rRNA gene. Genes responsible for the degradation of aromatic hydrocarbons and n-alkanes were found in all genomes, although different (and complementary) steps of the metabolic pathways were unique to each strain. Genes for naphthalene and toluene degradation were found in various strains. We annotated quinate degradation genes in LC_6, while LC_3 and LC_5 presented genes for biosurfactant and rhamnolipid biosynthesis. We also annotated genes related to beneficial mechanisms for corals, such as genes involved in nitrogen and DMSP metabolism, cobalamin biosynthesis and antimicrobial compounds production. Our findings reinforce the importance of using bacterial consortia for bioremediation approaches instead of single strains, due to their complementary genomic arsenals. We also propose a genome-based framework to select complementary ODB that can provide additional benefits to coral health.

Genomic Characterization of Aeromonas veronii Provides Insights into Taxonomic Assignment and Reveals Widespread Virulence and Resistance Genes throughout the World

Aeromonas veronii is a Gram-negative bacterial species that causes disease in fish and is nowadays increasingly recurrent in enteric infections of humans. This study was performed to characterize newly sequenced isolates by comparing them with complete genomes deposited at the NCBI (National Center for Biotechnology Information). Nine isolates from fish, environments, and humans from the São Francisco Valley (Petrolina, Pernambuco, Brazil) were sequenced and compared with complete genomes available in public databases to gain insight into taxonomic assignment and to better understand virulence and resistance profiles of this species within the One Health context. One local genome and four NCBI genomes were misidentified as A. veronii. A total of 239 virulence genes were identified in the local genomes, with most encoding adhesion, motility, and secretion systems. In total, 60 genes involved with resistance to 22 classes of antibiotics were identified in the genomes, including mcr-7 and cphA. The results suggest that the use of methods such as ANI is essential to avoid misclassification of the genomes. The virulence content of A. veronii from local isolates is similar to those complete genomes deposited at the NCBI. Genes encoding colistin resistance are widespread in the species, requiring greater attention for surveillance systems.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Reference-free phylogeny from sequencing data

Motivation

Clustering of genetic sequences is one of the key parts of bioinformatics analyses. Resulting phylogenetic trees are beneficial for solving many research questions, including tracing the history of species, studying migration in the past, or tracing a source of a virus outbreak. At the same time, biologists provide more data in the raw form of reads or only on contig-level assembly. Therefore, tools that are able to process those data without supervision need to be developed.

Results

In this paper, we present a tool for reference-free phylogeny capable of handling data where no mature-level assembly is available. The tool allows distance calculation for raw reads, contigs, and the combination of the latter. The tool provides an estimation of the Levenshtein distance between the sequences, which in turn estimates the number of mutations between the organisms. Compared to the previous research, the novelty of the method lies in a newly proposed combination of the read and contig measures, a new method for read-contig mapping, and an efficient embedding of contigs.

Comparative genomic analysis of the Dietzia genus: an insight into genomic diversity, and adaptation

Dietzia strains are widely distributed in the environment, presenting an opportunistic role, and some species have undetermined taxonomic characteristics. Here, we propose the existence of errors in the classification of species in this genus using comparative genomics. We performed ANI, dDDH, pangenome and genomic plasticity analyses better to elucidate the phylogenomic relationships between Dietzia strains. For this, we used 55 genomes of Dietzia downloaded from public databases that were combined with a newly sequenced. Sequence analysis of a phylogenetic tree based on genome similarity comparisons and dDDH, ANI analyses supported grouping different Dietzia species into four distinct groups. The pangenome analysis corroborated the classification of these groups, supporting the idea that some species of Dietzia could be reassigned in a possible classification into three distinct species, each containing less variability than that found within the global pangenome of all strains. Additionally, analysis of genomic plasticity based on groups containing Dietzia strains found differences in the presence and absence of symbiotic Islands and pathogenic islands related to their isolation site. We propose that the comparison of pangenome subsets together with phylogenomic approaches can be used as an alternative for the classification and differentiation of new species of the genus Dietzia.

Analysis of Corynebacterium silvaticum genomes from Portugal reveals a single cluster and a clade suggested to produce diphtheria toxin

Background

Corynebacterium silvaticum is a pathogenic, gram-positive bacterial species that causes caseous lymphadenitis in wild boars, domestic pigs and roe deer in Western Europe. It can affect animal production and cause zoonosis. Genome analysis has suggested that one strain from Portugal and one from Austria could probably produce the diphtheria toxin (DT), which inhibits protein synthesis and can cause death.

Methods

To further investigate the species genetic diversity and probable production of DT by Portuguese strains, eight isolates from this country were sequenced and compared to 38 public ones.

Results

Strains from Portugal are monophyletic, nearly identical, form a unique cluster and have 27 out of 36 known Corynebacterium virulence or niche factors. All of them lack a frameshift in the tox gene and were suggested to produce DT. A phylogenetic analysis shows that the species has diverged into two clades. Clade 1 is composed of strains that were suggested to have the ability to produce DT, represented by the monophyletic strains from Portugal and strain 05-13 from Austria. Clade 2 is composed of strains unable to produce DT due to a frameshifted tox gene. The second clade is represented by strains from Austria, Germany and Switzerland. Ten genome clusters were detected, in which strains from Germany are the most diverse. Strains from Portugal belong to an exclusive cluster. The pangenome has 2,961 proteins and is nearly closed (α = 0.968). Exclusive genes shared by clusters 1 and 2, and Portuguese strains are probably not related to disease manifestation as they share the same host but could play a role in their extra-host environmental adaptation. These results show the potential of the species to cause zoonosis, possibly diphtheria. The identified clusters, exclusively shaded genes, and exclusive STs identified in Portugal could be applied in the identification and epidemiology of the species.

Hybrid Genomic Analysis of Salmonella enterica Serovar Enteritidis SE3 Isolated from Polluted Soil in Brazil

In Brazil, Salmonella enterica serovar Enteritidis is a significant health threat. Salmonella enterica serovar Enteritidis SE3 was isolated from soil at the Subaé River in Santo Amaro, Brazil, a region contaminated with heavy metals and organic waste. Illumina HiSeq and Oxford Nanopore Technologies MinION sequencing were used for de novo hybrid assembly of the Salmonella SE3 genome. This approach yielded 10 contigs with 99.98% identity with S. enterica serovar Enteritidis OLF-SE2-98984-6. Twelve Salmonella pathogenic islands, multiple virulence genes, multiple antimicrobial gene resistance genes, seven phage defense systems, seven prophages and a heavy metal resistance gene were encoded in the genome. Pangenome analysis of the S. enterica clade, including Salmonella SE3, revealed an open pangenome, with a core genome of 2137 genes. Our study showed the effectiveness of a hybrid sequence assembly approach for environmental Salmonella genome analysis using HiSeq and MinION data. This approach enabled the identification of key resistance and virulence genes, and these data are important to inform the control of Salmonella and heavy metal pollution in the Santo Amaro region of Brazil.

Detection, Detrimental Effects, and Transmission Pathways of the Pathogenic Bacterium Acaricomes phytoseiuli in Commercial Predatory Mites

Arthropod pathogens and other microorganisms have been documented from mass production systems aimed at producing natural enemies for pest control. If losses due to pathogens are encountered in such systems, they could lead to uneconomical production of natural enemies and/or a loss of predator quality for effective field control of target pests. Here, we identify and describe the laboratory transmission of a bacterial pathogen, Acaricomes phytoseiuli, in a Chinese strain of the local predatory mite Neoseiulus californicus following contact with Phytoseiulus persimilis, a predatory mite imported for the control of small sap-sucking pests in greenhouses. Analysis of the 16S rRNA gene and whole-genome sequences of A. phytoseiuli isolated from the Chinese strain of N. californicus showed 99.6 and 99.78% similarity, respectively, to the pathogen isolated from a European population (DSM14247 strain). This is the first report of P. persimilis infected with A. phytoseiuli outside Europe and transmitting to a local predatory mite species. A. phytoseiuli severely damaged local N. californicus, leading to a dorso-ventrally flattened body and reduced prey consumption and reproduction as well as early death. Through fluorescence in situ hybridization, A. phytoseiuli was shown to accumulate in the digestive tract of mites and in the oviductal bulb of adult females. Infected males had no obvious symptoms, but they still were able to pass on the infection to healthy females through contact and mating. The pathogen was transmitted vertically to offspring by either infected parent through adherence to eggshells. A. phytoseiuli could also persist in other herbivorous arthropods from the same habitat, suggesting wider potential risks. Our study highlights pathogen risk to predators in local biocontrol industries through pathogen spread from imported material. IMPORTANCE Predatory mites are important natural enemies for biological control of pests, but mass rearing of the mites can be affected by pathogens. In particular, the mite pathogen Acaricomes phytoseiuli may pose a threat to predatory mite production, and we have now detected this pathogen in China. We explored the pathogen's transmission within species, its ability to transfer to a locally important predatory mite species, and symptoms following transfer. The detection of A. phytoseiuli and its ability to transfer to a local predator where it reduces performance highlight the importance of ongoing monitoring and hygiene in the predatory mite industry.

Independent evolution of multi-dominant viral genome species observed in a hepatitis C virus carrier

The viral genome quasispecies composition of hepatitis C virus (HCV) could have important implications to viral pathogenesis and resistance to anti-viral treatment. The purpose of the present study was to profile the HCV RNA quasispecies. We developed a strategy to determine the full-length HCV genome sequences co-existing within a single patient serum by using next-generation sequencing technologies. The isolated viral clones were divided into the groups that can be distinguished by core amino acid 70 substitution. Subsequently, we determined HCV full-length genome sequences of three independent dominant species co-existing in the sequential serum with a 7-year interval. From phylogenetic analysis, these dominant species evolved independently. Our study demonstrated that multiple dominant species co-existed in patient sera and evolved independently.

•

HCV RNA genome forms quasispecies which may contribute viral pathogenesis.

•

A strategy was established to determine the full-length HCV genome sequences co-existing within the sera of a single patient.

•

Multiple dominant viral species co-existed in patient sera and evolved independently.

Isolation and characterization of two novel serotypes of Tibet orbivirus from Culicoides and sentinel cattle in Yunnan Province of China

Tibet orbivirus (TIBOV), a new candidate of Orbivirus genus, was initially isolated from mosquitoes in Tibet in 2009 and subsequently from both Culicoides and mosquitoes in several provinces of China and Japan. Little is known about the origin, genetic diversity, dissemination and pathogenicity of TIBOV, although its potential threat to animal health has been acknowledged. In this study, two viruses, V290/YNSZ and V298/YNJH, were isolated from the Culicoides and sentinel cattle in Yunnan Province. Their genome sequences, cell tropism in mammalian and insect cell lines along with pathogenicity in suckling mice were determined. Genome phylogenetic analyses confirmed their classification as TIBOV species; however, OC1 proteins of the V290/YNSZ and V298/YNJH shared maximum sequence identities of 31.5% and 33.9% with other recognized TIBOV serotypes (TIBOV-1 to TIBOV-4) and formed two monophyletic branches in phylogenetic tree, indicating they represented two novel TIBOV serotypes which were tentatively designated as TIBOV-5 and TIBOV-6. The viruses replicated robustly in BHK, Vero and C6/36 cells and triggered overt clinical symptoms in suckling mice after intracerebral inoculation, causing mortality of 100% and 25%. Cross-sectional epidemiology analysis revealed silent circulation of TIBOV in Yunnan Province with overall prevalence of 16.4% (18/110) in cattle, 10.8% (13/120) in goats and 5.5% (6/110) in swine. The prevalence patterns of four investigated TIBOV serotypes (TIBOV-1, -2, -5 and 6) differed from each one another, with their positive rates ranging from 8.2% (9/110) for TIBOV-2 in cattle to 0.9% (1/110) for TIBOV-1 and TIBOV-5 in cattle and swine. Our findings provided new insights for diversity, pathogenicity and epidemiology of TIBOV and formed a basis for future studies addressing the geographical distribution and the zoonotic potential of TIBOV.

Conventional and Omics Approaches for Understanding the Abiotic Stress Response in Cereal Crops-An Updated Overview

Cereals have evolved various tolerance mechanisms to cope with abiotic stress. Understanding the abiotic stress response mechanism of cereal crops at the molecular level offers a path to high-yielding and stress-tolerant cultivars to sustain food and nutritional security. In this regard, enormous progress has been made in the omics field in the areas of genomics, transcriptomics, and proteomics. Omics approaches generate a massive amount of data, and adequate advancements in computational tools have been achieved for effective analysis. The combination of integrated omics and bioinformatics approaches has been recognized as vital to generating insights into genome-wide stress-regulation mechanisms. In this review, we have described the self-driven drought, heat, and salt stress-responsive mechanisms that are highlighted by the integration of stress-manipulating components, including transcription factors, co-expressed genes, proteins, etc. This review also provides a comprehensive catalog of available online omics resources for cereal crops and their effective utilization. Thus, the details provided in the review will enable us to choose the appropriate tools and techniques to reduce the negative impacts and limit the failures in the intensive crop improvement study.

Unraveling the Genomic Potential of the Thermophilic Bacterium Anoxybacillus flavithermus from an Antarctic Geothermal Environment

Antarctica is a mosaic of extremes. It harbors active polar volcanoes, such as Deception Island, a marine stratovolcano having notable temperature gradients over very short distances, with the temperature reaching up to 100 °C near the fumaroles and subzero temperatures being noted in the glaciers. From the sediments of Deception Island, we isolated representatives of the genus Anoxybacillus, a widely spread genus that is mainly encountered in thermophilic environments. However, the phylogeny of this genus and its adaptive mechanisms in the geothermal sites of cold environments remain unknown. To the best of our knowledge, this is the first study to unravel the genomic features and provide insights into the phylogenomics and metabolic potential of members of the genus Anoxybacillus inhabiting the Antarctic thermophilic ecosystem. Here, we report the genome sequencing data of seven A. flavithermus strains isolated from two geothermal sites on Deception Island, Antarctic Peninsula. Their genomes were approximately 3.0 Mb in size, had a G + C ratio of 42%, and were predicted to encode 3500 proteins on average. We observed that the strains were phylogenomically closest to each other (Average Nucleotide Identity (ANI) > 98%) and to A. flavithermus (ANI 95%). In silico genomic analysis revealed 15 resistance and metabolic islands, as well as genes related to genome stabilization, DNA repair systems against UV radiation threats, temperature adaptation, heat- and cold-shock proteins (Csps), and resistance to alkaline conditions. Remarkably, glycosyl hydrolase enzyme-encoding genes, secondary metabolites, and prophage sequences were predicted, revealing metabolic and cellular capabilities for potential biotechnological applications.

Benchmarking and Assessment of Eight De Novo Genome Assemblers on Viral Next-Generation Sequencing Data, Including the SARS-CoV-2

Viral genomics has become crucial in clinical diagnostics and ecology, not to mention to stem the COVID-19 pandemic. Whole-genome sequencing (WGS) is pivotal in gaining an improved understanding of viral evolution, genomic epidemiology, infectious outbreaks, pathobiology, clinical management, and vaccine development. Genome assembly is one of the crucial steps in WGS data analyses. A series of different assemblers has been developed with the advent of high-throughput next-generation sequencing (NGS). Various studies have reported the evaluation of these assembly tools on distinct datasets; however, these lack data from viral origin. In this study, we performed a comparative evaluation and benchmarking of eight de novo assemblers: SOAPdenovo, Velvet, assembly by short sequences (ABySS), iterative De Bruijn graph assembler (IDBA), SPAdes, Edena, iterative virus assembler, and VICUNA on the viral NGS data from distinct Illumina (GAIIx, Hiseq, Miseq, and Nextseq) platforms. WGS data of diverse viruses, that is, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), dengue virus 3, human immunodeficiency virus 1, hepatitis B virus, human herpesvirus 8, human papillomavirus 16, rhinovirus A, and West Nile virus, were utilized to assess these assemblers. Performance metrics such as genome fraction recovery, assembly lengths, NG50, N50, contig length, contig numbers, mismatches, and misassemblies were analyzed. Overall, three assemblers, that is, SPAdes, IDBA, and ABySS, performed consistently well, including for genome assembly of SARS-CoV-2. These assembly methods should be considered and recommended for future studies of viruses. The study also suggests that implementing two or more assembly approaches should be considered in viral NGS studies, especially in clinical settings. Taken together, the benchmarking of eight de novo genome assemblers reported in this study can inform future public health and ecology research concerning the viruses, the COVID-19 pandemic, and viral outbreaks.

Polyphasic Analysis Reveals Potential Petroleum Hydrocarbon Degradation and Biosurfactant Production by Rare Biosphere Thermophilic Bacteria From Deception Island, an Active Antarctic Volcano

Extreme temperature gradients in polar volcanoes are capable of selecting different types of extremophiles. Deception Island is a marine stratovolcano located in maritime Antarctica. The volcano has pronounced temperature gradients over very short distances, from as high as 100°C in the fumaroles to subzero next to the glaciers. These characteristics make Deception a promising source of a variety of bioproducts for use in different biotechnological areas. In this study, we isolated thermophilic bacteria from sediments in fumaroles at two geothermal sites on Deception Island with temperatures between 50 and 100°C, to evaluate the potential capacity of these bacteria to degrade petroleum hydrocarbons and produce biosurfactants under thermophilic conditions. We isolated 126 thermophilic bacterial strains and identified them molecularly as members of genera Geobacillus, Anoxybacillus, and Brevibacillus (all in phylum Firmicutes). Seventy-six strains grew in a culture medium supplemented with crude oil as the only carbon source, and 30 of them showed particularly good results for oil degradation. Of 50 strains tested for biosurfactant production, 13 showed good results, with an emulsification index of 50% or higher of a petroleum hydrocarbon source (crude oil and diesel), emulsification stability at 100°C, and positive results in drop-collapse, oil spreading, and hemolytic activity tests. Four of these isolates showed great capability of degrade crude oil: FB2_38 (Geobacillus), FB3_54 (Geobacillus), FB4_88 (Anoxybacillus), and WB1_122 (Geobacillus). Genomic analysis of the oil-degrading and biosurfactant-producer strain FB4_88 identified it as Anoxybacillus flavithermus, with a high genetic and functional diversity potential for biotechnological applications. These initial culturomic and genomic data suggest that thermophilic bacteria from this Antarctic volcano have potential applications in the petroleum industry, for bioremediation in extreme environments and for microbial enhanced oil recovery (MEOR) in reservoirs. In addition, recovery of small-subunit rRNA from metagenomes of Deception Island showed that Firmicutes is not among the dominant phyla, indicating that these low-abundance microorganisms may be important for hydrocarbon degradation and biosurfactant production in the Deception Island volcanic sediments.

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Uziflies (Family: Tachinidae) are dipteran endoparasites of sericigenous insects which cause major economic loss in the silk industry globally. Here, we are presenting the first full mitogenome of Blepharipa sp. (Acc: KY644698, 15,080 bp, A + T = 78.41%), a dipteran parasitoid of Muga silkworm (Antheraea assamensis) found in the Indian states of Assam and Meghalaya. This study has confirmed that Blepharipa sp. mitogenome gene content and arrangement is similar to other Tachinidae and Sarcophagidae flies of Oestroidea superfamily, typical of ancestral Diptera. Although, Calliphoridae and Oestridae flies have undergone tRNA translocation and insertion, forming unique intergenic spacers (IGS) and overlapping regions (OL) and a few of them (IGS, OL) have been conserved across Oestroidea flies. The Tachinidae mitogenomes exhibit more AT content and AT biased codons in their protein-coding genes (PCGs) than the Oestroidea counterpart. About 92.07% of all (3722) codons in PCGs of this new species have A/T in their 3rd codon position. The high proportion of AT and repeats in the control region (CR) affects sequence coverage, resulting in a short CR (Blepharipa sp.: 168 bp) and a smaller tachinid mitogenome. Our research unveils those genes with a high AT content had a reduced effective number of codons, leading to high codon usage bias. The neutrality test shows that natural selection has a stronger influence on codon usage bias than directed mutational pressure. This study also reveals that longer PCGs (e.g., nad5, cox1) have a higher codon usage bias than shorter PCGs (e.g., atp8, nad4l). The divergence rates increase nonlinearly as AT content at the 3rd codon position increases and higher rate of synonymous divergence than nonsynonymous divergence causes strong purifying selection. The phylogenetic analysis explains that Blepharipa sp. is well suited in the family of insectivorous tachinid maggots. It's possible that biased codon usage in the Tachinidae family reduces the effective number of codons, and purifying selection retains the core functions in their mitogenome, which could help with efficient metabolism in their endo-parasitic life style and survival strategy.

The Space-Exposed Kombucha Microbial Community Member Komagataeibacter oboediens Showed Only Minor Changes in Its Genome After Reactivation on Earth

Komagataeibacter is the dominant taxon and cellulose-producing bacteria in the Kombucha Microbial Community (KMC). This is the first study to isolate the K. oboediens genome from a reactivated space-exposed KMC sample and comprehensively characterize it. The space-exposed genome was compared with the Earth-based reference genome to understand the genome stability of K. oboediens under extraterrestrial conditions during a long time. Our results suggest that the genomes of K. oboediens IMBG180 (ground sample) and K. oboediens IMBG185 (space-exposed) are remarkably similar in topology, genomic islands, transposases, prion-like proteins, and number of plasmids and CRISPR-Cas cassettes. Nonetheless, there was a difference in the length of plasmids and the location of cas genes. A small difference was observed in the number of protein coding genes. Despite these differences, they do not affect any genetic metabolic profile of the cellulose synthesis, nitrogen-fixation, hopanoid lipids biosynthesis, and stress-related pathways. Minor changes are only observed in central carbohydrate and energy metabolism pathways gene numbers or sequence completeness. Altogether, these findings suggest that K. oboediens maintains its genome stability and functionality in KMC exposed to the space environment most probably due to the protective role of the KMC biofilm. Furthermore, due to its unaffected metabolic pathways, this bacterial species may also retain some promising potential for space applications.

Bacteriocin Producing Streptococcus agalactiae Strains Isolated from Bovine Mastitis in Brazil

Antibiotic resistance is one of the biggest health challenges of our time. We are now facing a post-antibiotic era in which microbial infections, currently treatable, could become fatal. In this scenario, antimicrobial peptides such as bacteriocins represent an alternative solution to traditional antibiotics because they are produced by many organisms and can inhibit bacteria, fungi, and/or viruses. Herein, we assessed the antimicrobial activity and biotechnological potential of 54 Streptococcus agalactiae strains isolated from bovine mastitis. Deferred plate antagonism assays revealed an inhibition spectrum focused on species of the genus Streptococcus—namely, S. pyogenes, S. agalactiae, S. porcinus, and S. uberis. Three genomes were successfully sequenced, allowing for their taxonomic confirmation via a multilocus sequence analysis (MLSA). Virulence potential and antibiotic resistance assessments showed that strain LGMAI_St_08 is slightly more pathogenic than the others. Moreover, the mreA gene was identified in the three strains. This gene is associated with resistance against erythromycin, azithromycin, and spiramycin. Assessments for secondary metabolites and antimicrobial peptides detected the bacteriocin zoocin A. Finally, comparative genomics evidenced high similarity among the genomes, with more significant similarity between the LGMAI_St_11 and LGMAI_St_14 strains. Thus, the current study shows promising antimicrobial and biotechnological potential for the Streptococcus agalactiae strains.

Paenibacillus piscarius sp. nov., a novel nitrogen-fixing species isolated from the gut of the armored catfish Parotocinclus maculicauda

A Gram-positive, nitrogen-fixing and endospore-forming strain, designated P121^T, was isolated from the gut of the armored catfish (Parotocinclus maculicauda) and identified as a member of the genus Paenibacillus based on the sequences of the 16S rRNA encoding gene, rpoB, gyrB and nifH genes and phenotypic analyses. The most closely related species to strain P121^T were Paenibacillus rhizoplanae DSM 103993^T, Paenibacillus silagei DSM 101953^T and Paenibacillus borealis DSM 13188^T, with similarity values of 98.9, 98.3 and 97.6%, respectively, based on 16S rRNA gene sequences. Genome sequencing revealed a genome size of 7,513,698 bp, DNA G + C content of 53.9 mol% and the presence of the structural nitrogenase encoding genes (nifK, nifD and nifH) and of other nif genes necessary for nitrogen fixation. Digital DNA-DNA hybridization (dDDH) experiments and average nucleotide identity (ANI) analyses between strain P121^T and the type strains of the closest species demonstrated that the highest values were below the thresholds of 70% dDDH (42.3% with P. borealis) and 95% ANI (84.28% with P. silagei) for bacterial species delineation, indicating that strain P121^T represents a distinct species. Its major cellular fatty acid was anteiso-C_15:0 (42.4%), and the major isoprenoid quinone was MK-7. Based on physiological, genomic, biochemical and chemotaxonomic characteristics, we propose that strain P121^T represents a novel species for which the name Paenibacillus piscarius sp. nov. is proposed (type strain = DSM 25072 = LFB-Fiocruz 1636).

SLAG: A Program for Seeded Local Assembly of Genes in Complex Genomes

Although finished genomes have become more common, there is still a need for assemblies of individual genes or chromosomal regions when only unassembled reads are available. slag (Seeded Local Assembly of Genes) fulfils this need by performing iterative local assembly based on cycles of matching‐read retrieval with blast and assembly with cap3, phrap, spades, canu or unicycler. The target sequence can be nucleotide or protein. Read fragmentation allows slag to use phrap or cap3 to assemble long reads at lower coverage (e.g., 5×) than is possible with canu or unicycler. In simple, nonrepetitive genomes, a slag assembly can cover a whole chromosome, but in complex genomes the growth of target‐matching contigs is limited as additional reads are consumed by consensus contigs consisting of repetitive elements. Apart from genomic complexity, contig length and correctness depend on read length and accuracy. With pyrosequencing or Illumina reads, slag‐assembled contigs are accurate enough to allow design of PCR primers, while contigs assembled from Oxford Nanopore or pre‐HiFi Pacific Biosciences long reads are generally only accurate enough to design baiting sequences for further targeted sequencing. In an application with real reads, slag successfully extended sequences for four wheat genes, which were verified by cloning and Sanger sequencing of overlapping amplicons. slag is a robust alternative to atram2 for local assemblies, especially for read sets with less than 20× coverage. slag is freely available at https://github.com/cfcrane/SLAG.

Genome Sequence of Pseudomonas sp. Strain LAP_36, A Rhizosphere Bacterium Isolated from King George Island, Antarctica

Pseudomonas sp. strain LAP_36 was isolated from rhizosphere soil from Deschampsia antarctica on King George Island, South Shetland Islands, Antarctica. Here, we report on its draft genome sequence, which consists of 8,794,771 bp with 60.0% GC content and 8,011 protein-coding genes.

A bioinformatic platform to integrate target capture and whole genome sequences of various read depths for phylogenomics

The increasing availability of short‐read whole genome sequencing (WGS) provides unprecedented opportunities to study ecological and evolutionary processes. Although loci of interest can be extracted from WGS data and combined with target sequence data, this requires suitable bioinformatic workflows. Here, we test different assembly and locus extraction strategies and implement them into secapr, a pipeline that processes short‐read data into multilocus alignments for phylogenetics and molecular ecology analyses. We integrate the processing of data from low‐coverage WGS (<30×) and target sequence capture into a flexible framework, while optimizing de novo contig assembly and loci extraction. Specifically, we test different assembly strategies by contrasting their ability to recover loci from targeted butterfly protein‐coding genes, using four data sets: a WGS data set across different average coverages (10×, 5× and 2×) and a data set for which these loci were enriched prior to sequencing via target sequence capture. Using the resulting de novo contigs, we account for potential errors within contigs and infer phylogenetic trees to evaluate the ability of each assembly strategy to recover species relationships. We demonstrate that choosing multiple sizes of kmer simultaneously for assembly results in the highest yield of extracted loci from de novo assembled contigs, while data sets derived from sequencing read depths as low as 5× recovers the expected species relationships in phylogenetic trees. By making the tested assembly approaches available in the secapr pipeline, we hope to inspire future studies to incorporate complementary data and make an informed choice on the optimal assembly strategy.

Discriminating Clonotypes of Influenza A Virus Genes by Nanopore Sequencing

Influenza viruses still pose a serious threat to humans, and we have not yet been able to effectively predict future pandemic strains and prepare vaccines in advance. One of the main reasons is the high genetic diversity of influenza viruses. We do not know the individual clonotypes of a virus population because some are the majority and others make up only a small fraction of the population. First-generation (FGS) and next-generation sequencing (NGS) technologies have inherent limitations that are unable to resolve a minority clonotype’s information in the virus population. Third-generation sequencing (TGS) technologies with ultra-long reads have the potential to solve this problem but have a high error rate. Here, we evaluated emerging direct RNA sequencing and cDNA sequencing with the MinION platform and established a novel approach that combines the high accuracy of Illumina sequencing technology and long reads of nanopore sequencing technology to resolve both variants and clonotypes of influenza virus. Furthermore, a new program was written to eliminate the effect of nanopore sequencing errors for the analysis of the results. By using this pipeline, we identified 47 clonotypes in our experiment. We conclude that this approach can quickly discriminate the clonotypes of virus genes, allowing researchers to understand virus adaptation and evolution at the population level.

Complete Genome Sequence of the Polysaccharide-Degrading Rumen Bacterium Pseudobutyrivibrio xylanivorans MA3014 Reveals an Incomplete Glycolytic Pathway

Bacterial species belonging to the genus Pseudobutyrivibrio are important members of the rumen microbiome contributing to the degradation of complex plant polysaccharides. Pseudobutyrivibrio xylanivorans MA3014 was selected for genome sequencing to examine its ability to breakdown and utilize plant polysaccharides. The complete genome sequence of MA3014 is 3.58 Mb, consists of three replicons (a chromosome, chromid, and plasmid), has an overall G + C content of 39.6%, and encodes 3,265 putative protein-coding genes (CDS). Comparative pan-genomic analysis of all cultivated and currently available P. xylanivorans genomes has revealed a strong correlation of orthologous genes within this rumen bacterial species. MA3014 is metabolically versatile and capable of growing on a range of simple mono- or oligosaccharides derived from complex plant polysaccharides such as pectins, mannans, starch, and hemicelluloses, with lactate, butyrate, and formate as the principal fermentation end products. The genes encoding these metabolic pathways have been identified and MA3014 is predicted to encode an extensive range of Carbohydrate-Active enZYmes with 78 glycoside hydrolases, 13 carbohydrate esterases, and 54 glycosyl transferases, suggesting an important role in solubilization of plant matter in the rumen.

Efficient assembly consensus algorithms for divergent contig sets

Assembly is a fundamental task in genome sequencing, and many assemblers have been made available in the last decade. Because of the wide range of possible choices, it can be hard to determine which tool or parameter to use for a specific genome sequencing project. In this paper, we propose a consensus approach that takes the best parts of several contigs datasets produced by different methods, and combines them into a better assembly. This amounts to orienting and ordering sets of contigs, which can be viewed as an optimization problem where the aim is to find an alignment of two fragmented strings that maximizes an arbitrary scoring function between matched characters. In this work, we investigate the computational complexity of this problem. We first show that it is NP-hard, even in an alphabet with only two symbols and with all scores being either 0 or 1. On the positive side, we propose an efficient, quadratic time algorithm that achieves approximation factor 3.

Draft Genome Sequences of 13 Vibrio cholerae Strains from the Rio Grande Delta

Vibrio cholerae is the etiologic agent of cholera, an acute and often fatal diarrheal disease that affects millions globally. We report the draft genome sequences of 13 non-O1/O139 V. cholerae strains isolated from the Rio Grande Delta in Texas. These genomes will aid future analyses of environmental serovars.

Novel putative bluetongue virus serotype 29 isolated from inapparently infected goat in Xinjiang of China

Bluetongue virus (BTV) is the 'type' species of the genus Orbivirus causing bluetongue (BT) in sheep, bovine and other ruminants. Twenty-four serotypes and several atypical serotypes of BTV were identified worldwide. In present study, a novel strain of BTV (V196/XJ/2014) was isolated from an asymptomatic sentinel goat in Yuli County, Xinjiang of China. Serotype identification of this isolate exhibited uniform negative results by serotype-specific conventional RT-PCR and real-time RT-PCR for BTV-1 to BTV-27, and virus neutralization tests using reference sera of BTV-1 to BTV-24. Genomic analysis showed V196/XJ/2014 grouped with atypical serotypes of BTV-25 to BTV-28, BTV-X/XJ1407, BTV-X/ITL2015 and BTV-Y/TUN2017, while segment 2 and VP2 protein of V196/XJ/2014 shared <63.4%/61.4% nucleic acids and amino acids sequence identities with other recognized BTV serotypes and its segment 2 formed a separate 'nucleotype' in phylogenetic tree. These results indicated V196/XJ/2014 does not belong to any reported serotypes of BTV. Further studies of infectivity and pathogenicity showed that goats infected with V196/XJ/2014 did not exhibit observed clinical symptoms, but high level of virus amplification and homologous neutralization antibodies were detected post-infection. Our studies suggested a novel putative serotype of BTV-29 was isolated in Xinjiang of China, which expands our knowledge about the diversity of BTV.

Probiogenomics of Lactobacillus delbrueckii subsp. lactis CIDCA 133: In Silico, In Vitro, and In Vivo Approaches

Lactobacillus delbrueckii subsp. lactis CIDCA 133 (CIDCA 133) has been reported as a potential probiotic strain, presenting immunomodulatory properties. This study investigated the possible genes and molecular mechanism involved with a probiotic profile of CIDCA 133 through a genomic approach associated with in vitro and in vivo analysis. Genomic analysis corroborates the species identification carried out by the classical microbiological method. Phenotypic assays demonstrated that the CIDCA 133 strain could survive acidic, osmotic, and thermic stresses. In addition, this strain shows antibacterial activity against Salmonella Typhimurium and presents immunostimulatory properties capable of upregulating anti-inflammatory cytokines Il10 and Tgfb1 gene expression through inhibition of Nfkb1 gene expression. These reported effects can be associated with secreted, membrane/exposed to the surface and cytoplasmic proteins, and bacteriocins-encoding genes predicted in silico. Furthermore, our results showed the genes and the possible mechanisms used by CIDCA 133 to produce their beneficial host effects and highlight its use as a probiotic microorganism.

Olea Europaea Geminivirus: A Novel Bipartite Geminivirid Infecting Olive Trees

In 2014, high-throughput sequencing of libraries of total DNA from olive trees allowed the identification of two geminivirus-like contigs. After conventional resequencing of the two genomic DNAs, their analysis revealed they belonged to the same viral entity, for which the provisional name of Olea europaea geminivirus (OEGV) was proposed. Although DNA-A showed a genome organization similar to that of New World begomoviruses, DNA-B had a peculiar ORF arrangement, consisting of a movement protein (MP) in the virion sense and a protein with unknown function on the complementary sense. Phylogenetic analysis performed either on full-length genome or on coat protein, replication associated protein (Rep), and MP sequences did not endorse the inclusion of this virus in any of the established genera in the family Geminiviridae. A survey of 55 plants revealed that the virus is widespread in Apulia (Italy) with 91% of the samples testing positive, although no correlation of OEGV with a disease or specific symptoms was encountered. Southern blot assay suggested that the virus is not integrated in the olive genome. The study of OEGV-derived siRNA obtained from small RNA libraries of leaves and fruits of three different cultivars, showed that the accumulation of the two genomic components is influenced by the plant genotype while virus-derived-siRNA profile is in line with other geminivirids reported in literature. Single-nucleotide polymorphism (SNP) analysis unveiled a low intra-specific variability.

Characterization of the first vaginal Lactobacillus crispatus genomes isolated in Brazil

BACKGROUND

Lactobacillus crispatus is the dominant species in the vaginal microbiota associated with health and considered a homeostasis biomarker. Interestingly, some strains are even used as probiotics. However, the genetic mechanisms of L. crispatus involved in the control of the vaginal microbiome and protection against bacterial vaginosis (BV) are not entirely known. To further investigate these mechanisms, we sequenced and characterized the first four L. crispatus genomes from vaginal samples from Brazilian women and used genome-wide association study (GWAS) and comparative analyses to identify genetic mechanisms involved in healthy or BV conditions and selective pressures acting in the vaginal microbiome.

METHODS

The four genomes were sequenced, assembled using ten different strategies and automatically annotated. The functional characterization was performed by bioinformatics tools comparing with known probiotic strains. Moreover, it was selected one representative strain (L. crispatus CRI4) for in vitro detection of phages by electron microscopy. Evolutionary analysis, including phylogeny, GWAS and positive selection were performed using 46 public genomes strains representing health and BV conditions.

RESULTS

Genes involved in probiotic effects such as lactic acid production, hydrogen peroxide, bacteriocins, and adhesin were identified. Three hemolysins and putrescine production were predicted, although these features are also present in other probiotic strains. The four genomes presented no plasmids, but 14 known families insertion sequences and several prophages were detected. However, none of the mobile genetic elements contained antimicrobial resistance genes. The genomes harbor a CRISPR-Cas subtype II-A system that is probably inactivated due to fragmentation of the genes csn2 and cas9. No genomic feature was associated with a health condition, perhaps due to its multifactorial characteristic. Five genes were identified as under positive selection, but the selective pressure remains to be discovered. In conclusion, the Brazilian strains investigated in this study present potential protective properties, although in vitro and in vivo studies are required to confirm their efficacy and safety to be considered for human use.

Engineered yeast genomes accurately assembled from pure and mixed samples

Yeast whole genome sequencing (WGS) lacks end-to-end workflows that identify genetic engineering. Here we present Prymetime, a tool that assembles yeast plasmids and chromosomes and annotates genetic engineering sequences. It is a hybrid workflow-it uses short and long reads as inputs to perform separate linear and circular assembly steps. This structure is necessary to accurately resolve genetic engineering sequences in plasmids and the genome. We show this by assembling diverse engineered yeasts, in some cases revealing unintended deletions and integrations. Furthermore, the resulting whole genomes are high quality, although the underlying assembly software does not consistently resolve highly repetitive genome features. Finally, we assemble plasmids and genome integrations from metagenomic sequencing, even with 1 engineered cell in 1000. This work is a blueprint for building WGS workflows and establishes WGS-based identification of yeast genetic engineering.

Distribution of Beta-Lactamase Producing Gram-Negative Bacterial Isolates in Isabela River of Santo Domingo, Dominican Republic

Bacteria carrying antibiotic resistance genes (ARGs) are naturally prevalent in lotic ecosystems such as rivers. Their ability to spread in anthropogenic waters could lead to the emergence of multidrug-resistant bacteria of clinical importance. For this study, three regions of the Isabela river, an important urban river in the city of Santo Domingo, were evaluated for the presence of ARGs. The Isabela river is surrounded by communities that do not have access to proper sewage systems; furthermore, water from this river is consumed daily for many activities, including recreation and sanitation. To assess the state of antibiotic resistance dissemination in the Isabela river, nine samples were collected from these three bluedistinct sites in June 2019 and isolates obtained from these sites were selected based on resistance to beta-lactams. Physico-chemical and microbiological parameters were in accordance with the Dominican legislation. Matrix-assisted laser desorption ionization-time of flight mass spectrometry analyses of ribosomal protein composition revealed a total of 8 different genera. Most common genera were as follows: Acinetobacter (44.6%) and Escherichia (18%). Twenty clinically important bacterial isolates were identified from urban regions of the river; these belonged to genera Escherichia (n = 9), Acinetobacter (n = 8), Enterobacter (n = 2), and Klebsiella (n = 1). Clinically important multi-resistant isolates were not obtained from rural areas. Fifteen isolates were selected for genome sequencing and analysis. Most isolates were resistant to at least three different families of antibiotics. Among beta-lactamase genes encountered, we found the presence of blaTEM, blaOXA, blaSHV, and blaKPC through both deep sequencing and PCR amplification. Bacteria found from genus Klebsiella and Enterobacter demonstrated ample repertoire of antibiotic resistance genes, including resistance from a family of last resort antibiotics reserved for dire infections: carbapenems. Some of the alleles found were KPC-3, OXA-1, OXA-72, OXA-132, CTX-M-55, CTX-M-15, and TEM-1.

Fundamentally different repetitive element composition of sex chromosomes in Rumex acetosa

BACKGROUND AND AIMS

Dioecious species with well-established sex chromosomes are rare in the plant kingdom. Most sex chromosomes increase in size but no comprehensive analysis of the kind of sequences that drive this expansion has been presented. Here we analyse sex chromosome structure in common sorrel (Rumex acetosa), a dioecious plant with XY1Y2 sex determination, and we provide the first chromosome-specific repeatome analysis for a plant species possessing sex chromosomes.

METHODS

We flow-sorted and separately sequenced sex chromosomes and autosomes in R. acetosa using the two-dimensional fluorescence in situ hybridization in suspension (FISHIS) method and Illumina sequencing. We identified and quantified individual repeats using RepeatExplorer, Tandem Repeat Finder and the Tandem Repeats Analysis Program. We employed fluorescence in situ hybridization (FISH) to analyse the chromosomal localization of satellites and transposons.

KEY RESULTS

We identified a number of novel satellites, which have, in a fashion similar to previously known satellites, significantly expanded on the Y chromosome but not as much on the X or on autosomes. Additionally, the size increase of Y chromosomes is caused by non-long terminal repeat (LTR) and LTR retrotransposons, while only the latter contribute to the enlargement of the X chromosome. However, the X chromosome is populated by different LTR retrotransposon lineages than those on Y chromosomes.

CONCLUSIONS

The X and Y chromosomes have significantly diverged in terms of repeat composition. The lack of recombination probably contributed to the expansion of diverse satellites and microsatellites and faster fixation of newly inserted transposable elements (TEs) on the Y chromosomes. In addition, the X and Y chromosomes, despite similar total counts of TEs, differ significantly in the representation of individual TE lineages, which indicates that transposons proliferate preferentially in either the paternal or the maternal lineage.

Methods for Proteogenomics Data Analysis, Challenges, and Scalability Bottlenecks: A Survey

Big Data Proteogenomics lies at the intersection of high-throughput Mass Spectrometry (MS) based proteomics and Next Generation Sequencing based genomics. The combined and integrated analysis of these two high-throughput technologies can help discover novel proteins using genomic, and transcriptomic data. Due to the biological significance of integrated analysis, the recent past has seen an influx of proteogenomic tools that perform various tasks, including mapping proteins to the genomic data, searching experimental MS spectra against a six-frame translation genome database, and automating the process of annotating genome sequences. To date, most of such tools have not focused on scalability issues that are inherent in proteogenomic data analysis where the size of the database is much larger than a typical protein database. These state-of-the-art tools can take more than half a month to process a small-scale dataset of one million spectra against a genome of 3 GB. In this article, we provide an up-to-date review of tools that can analyze proteogenomic datasets, providing a critical analysis of the techniques' relative merits and potential pitfalls. We also point out potential bottlenecks and recommendations that can be incorporated in the future design of these workflows to ensure scalability with the increasing size of proteogenomic data. Lastly, we make a case of how high-performance computing (HPC) solutions may be the best bet to ensure the scalability of future big data proteogenomic data analysis.

Features of sRNA biogenesis in rice revealed by genetic dissection of sRNA expression level

We previously conducted a QTL analysis of small RNA (sRNA) abundance in flag leaves of an immortalized rice F₂ (IMF2) population by aligning sRNA reads to the reference genome to quantify the expression levels of sRNAs. However, this approach missed about half of the sRNAs as only 50% of all sRNA reads could be uniquely aligned to the reference genome. Here, we quantified the expression levels of sRNAs and sRNA clusters without the use of a reference genome. QTL analysis of the expression levels of sRNAs and sRNA clusters confirmed the feasibility of this approach. sRNAs and sRNA clusters with identified QTLs were then aligned to the high-quality parental genomes of the IMF2 population to resolve the identified QTLs into local vs. distant regulation mode. We were able to detect new QTL hotspots by considering sRNAs aligned to multiple positions of the parental genomes and sRNAs unaligned to the parental genomes. We found that several local-QTL hotspots were caused by sequence variations in long inverted repeats, which probably function as precursors of sRNAs, between the two parental genomes. The expression levels of these sRNAs were significantly associated with the presence/absence of the long inverted repeats in the IMF2 population. Moreover, we found that the variations in whole-genome sRNA species composition among different IMF2s were attributed to sRNA biogenesis genes including OsDCL2b and OsRDR2. Our results highlight that genetic dissection of sRNA expression is a promising approach to disclose new components functioning in sRNA biogenesis and new mechanisms of sRNA biogenesis.

Complete Genome Sequence of the Acetic Acid Bacterium Acetobacter aceti NBRC 14818

We report here the complete genome sequence of the acetic acid bacterium Acetobacter aceti type strain NBRC 14818. The genome comprises a chromosome of 3,596,270 bp with 57.1% GC content and four plasmids/phages of 63,279 bp, 25,755 bp, 4,858 bp, and 2,964 bp, with an average GC content of 57.0%.

Draft Genome Sequence of Simplicillium aogashimaense 72-15.1, a Putative Endophyte of Brachiaria brizantha

Here, we report a 29-Mb draft genome sequence of strain 72-15.1 of Simplicillium aogashimaense (Cordycipitaceae, Hypocreales). Strain 72-15.1 was a unique fungal isolate obtained from asymptomatic tillers of the tropical grass Brachiaria brizantha.

Here, we report a 29-Mb draft genome sequence of strain 72-15.1 of Simplicillium aogashimaense (Cordycipitaceae, Hypocreales). Strain 72-15.1 was a unique fungal isolate obtained from asymptomatic tillers of the tropical grass Brachiaria brizantha.

Improvement of gibberellin production by a newly isolated Fusarium fujikuroi mutant

AIMS

To obtain and investigate the potential mechanism for GA₃ production in Fusarium fujikuroi GA-251, a high GA₃ producer.

METHODS AND RESULTS

Fusarium fujikuroi IMI 58289 was bred with Cobalt-60 (⁶⁰ Co) radiation and lithium chloride treatment. The best mutant strain GA-251 was obtained for the subsequent optimization of fermentation conditions. The yield of GA₃ by GA-251 was 2100 mg l^-1 , while the wild-type strain was 100 mg l^-1 , which is a 21-fold increase in the yield. To elucidate the mechanism of high GA₃ yield of GA-251, the genome was sequenced and compared with wild-type strain IMI 58289. The results showed 2295 single nucleotide polymorphisms, 1242 small indels and 30 structural variants. These mutations were analysed and enriched in the MAPK signalling pathway, the mRNA surveillance pathway and endocytosis. The potential reasons for the improved GA₃ biosynthesis were investigated.

CONCLUSIONS

The potential mechanism of high GA₃ yield was attributed to endocytosis pathway and histone modification proteins family.

SIGNIFICANCE AND IMPACT OF THE STUDY

A mutant strain GA-251 in this work that could potentially be utilized in the industrial yield of GA₃ . The comparative genome analysis would shed light onto the mechanism of yield improvement and be a theoretical guide for further metabolic engineering.

Modeling of the Coral Microbiome: the Influence of Temperature and Microbial Network

Host-associated microbial communities are shaped by extrinsic and intrinsic factors to the holobiont organism. Environmental factors and microbe-microbe interactions act simultaneously on the microbial community structure, making the microbiome dynamics challenging to predict. The coral microbiome is essential to the health of coral reefs and sensitive to environmental changes. Here, we develop a dynamic model to determine the microbial community structure associated with the surface mucus layer (SML) of corals using temperature as an extrinsic factor and microbial network as an intrinsic factor. The model was validated by comparing the predicted relative abundances of microbial taxa to the relative abundances of microbial taxa from the sample data. The SML microbiome from Pseudodiploria strigosa was collected across reef zones in Bermuda, where inner and outer reefs are exposed to distinct thermal profiles. A shotgun metagenomics approach was used to describe the taxonomic composition and the microbial network of the coral SML microbiome. By simulating the annual temperature fluctuations at each reef zone, the model output is statistically identical to the observed data. The model was further applied to six scenarios that combined different profiles of temperature and microbial network to investigate the influence of each of these two factors on the model accuracy. The SML microbiome was best predicted by model scenarios with the temperature profile that was closest to the local thermal environment, regardless of the microbial network profile. Our model shows that the SML microbiome of P. strigosa in Bermuda is primarily structured by seasonal fluctuations in temperature at a reef scale, while the microbial network is a secondary driver.IMPORTANCE Coral microbiome dysbiosis (i.e., shifts in the microbial community structure or complete loss of microbial symbionts) caused by environmental changes is a key player in the decline of coral health worldwide. Multiple factors in the water column and the surrounding biological community influence the dynamics of the coral microbiome. However, by including only temperature as an external factor, our model proved to be successful in describing the microbial community associated with the surface mucus layer (SML) of the coral P. strigosa The dynamic model developed and validated in this study is a potential tool to predict the coral microbiome under different temperature conditions.

Stem rust resistance in wheat is suppressed by a subunit of the mediator complex

Stem rust is an important disease of wheat that can be controlled using resistance genes. The gene SuSr-D1 identified in cultivar ‘Canthatch’ suppresses stem rust resistance. SuSr-D1 mutants are resistant to several races of stem rust that are virulent on wild-type plants. Here we identify SuSr-D1 by sequencing flow-sorted chromosomes, mutagenesis, and map-based cloning. The gene encodes Med15, a subunit of the Mediator Complex, a conserved protein complex in eukaryotes that regulates expression of protein-coding genes. Nonsense mutations in Med15b.D result in expression of stem rust resistance. Time-course RNAseq analysis show a significant reduction or complete loss of differential gene expression at 24 h post inoculation in med15b.D mutants, suggesting that transcriptional reprogramming at this time point is not required for immunity to stem rust. Suppression is a common phenomenon and this study provides novel insight into suppression of rust resistance in wheat.

Stem rust is an important disease of wheat and resistance present in some cultivars can be suppressed by the SuSr-D1 locus. Here the authors show that SuSr-D1 encodes a subunit of the Mediator Complex and that nonsense mutations are sufficient to abolish suppression and confer stem rust resistance.

Plasmid-mediated metronidazole resistance in Clostridioides difficile

Metronidazole was until recently used as a first-line treatment for potentially life-threatening Clostridioides difficile (CD) infection. Although cases of metronidazole resistance have been documented, no clear mechanism for metronidazole resistance or a role for plasmids in antimicrobial resistance has been described for CD. Here, we report genome sequences of seven susceptible and sixteen resistant CD isolates from human and animal sources, including isolates from a patient with recurrent CD infection by a PCR ribotype (RT) 020 strain, which developed resistance to metronidazole over the course of treatment (minimal inhibitory concentration [MIC] = 8 mg L⁻¹). Metronidazole resistance correlates with the presence of a 7-kb plasmid, pCD-METRO. pCD-METRO is present in toxigenic and non-toxigenic resistant (n = 23), but not susceptible (n = 563), isolates from multiple countries. Introduction of a pCD-METRO-derived vector into a susceptible strain increases the MIC 25-fold. Our finding of plasmid-mediated resistance can impact diagnostics and treatment of CD infections.

Cases of C. difficile (CD) resistant to metronidazole have been reported but the mechanism remains enigmatic. Here the authors identify a plasmid, which correlates with metronidazole resistance status in a large international collection of CD isolates, and demonstrate that the plasmid can confer metronidazole resistance.

Insect genomes: progress and challenges

In the wake of constant improvements in sequencing technologies, numerous insect genomes have been sequenced. Currently, 1219 insect genome-sequencing projects have been registered with the National Center for Biotechnology Information, including 401 that have genome assemblies and 155 with an official gene set of annotated protein-coding genes. Comparative genomics analysis showed that the expansion or contraction of gene families was associated with well-studied physiological traits such as immune system, metabolic detoxification, parasitism and polyphagy in insects. Here, we summarize the progress of insect genome sequencing, with an emphasis on how this impacts research on pest control. We begin with a brief introduction to the basic concepts of genome assembly, annotation and metrics for evaluating the quality of draft assemblies. We then provide an overview of genome information for numerous insect species, highlighting examples from prominent model organisms, agricultural pests and disease vectors. We also introduce the major insect genome databases. The increasing availability of insect genomic resources is beneficial for developing alternative pest control methods. However, many opportunities remain for developing data-mining tools that make maximal use of the available insect genome resources. Although rapid progress has been achieved, many challenges remain in the field of insect genomics.

Mapping Ethanol Tolerance in Budding Yeast Reveals High Genetic Variation in a Wild Isolate

Ethanol tolerance, a polygenic trait of the yeast Saccharomyces cerevisiae, is the primary factor determining industrial bioethanol productivity. Until now, genomic elements affecting ethanol tolerance have been mapped only at low resolution, hindering their identification. Here, we explore the genetic architecture of ethanol tolerance, in the F6 generation of an Advanced Intercrossed Line (AIL) mapping population between two phylogenetically distinct, but phenotypically similar, S. cerevisiae strains (a common laboratory strain and a wild strain isolated from nature). Under ethanol stress, 51 quantitative trait loci (QTLs) affecting growth and 96 QTLs affecting survival, most of them novel, were identified, with high resolution, in some cases to single genes, using a High-Resolution Mapping Package of methodologies that provided high power and high resolution. We confirmed our results experimentally by showing the effects of the novel mapped genes: MOG1, MGS1, and YJR154W. The mapped QTLs explained 34% of phenotypic variation for growth and 72% for survival. High statistical power provided by our analysis allowed detection of many loci with small, but mappable effects, uncovering a novel “quasi-infinitesimal” genetic architecture. These results are striking demonstration of tremendous amounts of hidden genetic variation exposed in crosses between phylogenetically separated strains with similar phenotypes; as opposed to the more common design where strains with distinct phenotypes are crossed. Our findings suggest that ethanol tolerance is under natural evolutionary fitness-selection for an optimum phenotype that would tend to eliminate alleles of large effect. The study provides a platform for development of superior ethanol-tolerant strains using genome editing or selection.