A whole-genome assembly of Drosophila

A Newly Isolated Strain Lysobacter brunescens YQ20 and Its Performance on Wool Waste Biodegradation

Wool keratin is difficult to degrade as comparing to feathers because of its tough secondary structure. In order to develop an approach for high-value utilization of wool fiber waste by keratinolytic microorganisms, which is produced from shearing, weaving, and industrial processing of wool, screening of wool-degrading bacterium with high degradation efficiency were performed in this study. To this end, Lysobacter brunescens YQ20 was identified and characterized. The optimized conditions for wool degradation were pH 9.0 and 37°C with 20% liquid volume of Erlenmeyer flask. After fermentation, 15 essential amino acids were detected when wool fiber waste was fermented. The total amino acids produced from 1% wool per hour were 13.7 mg/L. The concentration was 8.6-fold higher than that produced by the strain Stenotrophomonas maltophilia BBE11-1, which had previously been reported to have the highest wool-degrading capacity. Our study reports the first Lysobacter strain that exhibits efficient wool degradation and yields higher concentrations of amino acids than previously reported strains. Whole-genome sequencing indicated that there were 18 keratinase-like genes in the genome of YQ20, which exhibited a long evolutionary distance from those of Bacillus. Therefore, L. brunescens YQ20 may have applications in the environmentally friendly management of wool waste as fertilizer in agriculture.

The complete sequence of a human genome

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

Transfer of Human Microbiome to Drosophila Gut Model

Laboratory animals with human microbiome have increasingly been used to study the role of bacteria and host interaction. Drosophila melanogaster, as a model of microbiota-host interaction with high reproductive efficiency and high availability, has always been lacking studies of interaction with human gut microbiome. In this study, we attempted to use antibiotic therapy and human fecal exposure strategy to transfer the human microbiome to the drosophila. The method includes depleting the original intestinal bacteria using a broad-spectrum antibiotic and then introducing human microorganisms by a diet supplemented with donor’s fecal samples. The sequencing results showed that 80–87.5% of the OTUs (Operational Taxonomic Units) from donor feces were adopted by the recipient drosophila following 30 days of observation. In comparison to females, the male recipient drosophila inherited more microbiota from the donor feces and had significantly increased lifespan as well as improved vertical climbing ability. Furthermore, distinctly differential expression patterns for age and insulin-like signaling-related genes were obtained for the male vs. female recipients. Only the male drosophila offspring acquired the characteristics of the donor fecal microbiota.

Comparative genomics and evolutionary analysis of Lactococcus garvieae isolated from human endocarditis

Lactococcus garvieae is a well-known pathogen of fish, but is rarely involved in infections in humans and other mammals. In humans, the main clinical manifestation of L. garvieae infections is endocarditis usually related to the ingestion of contaminated food, such as undercooked fish and shellfish. This study presents the first complete genomic sequence of a clinical L. garvieae strain isolated from a patient with endocarditis and its comparative analysis with other genomes. This human isolate contains a circular chromosome of 2 099 060 bp and one plasmid of 50 557 bp. In comparison with other fully sequenced L. garvieae strains, the chromosomal DNA of L. garvieae Lg-Granada carries a low proportion of insertion sequence elements and a higher number of putative prophages. Our results show that, in general, L. garvieae is a highly recombinogenic species with an open pangenome in which almost 30 % of its genome has undergone horizontal transfers. Within the genus Lactococcus, L. lactis is the main donor of genetic components to L. garvieae but, taking Lg-Granada as a representative, this bacterium tends to import more genes from Bacilli taxa than from other Lactococcus species.

Voltage-gated proton channels in polyneopteran insects

Voltage‐gated proton channels (H_V1) are expressed in eukaryotes, including basal hexapods and polyneopteran insects. However, currently, there is little known about H_V1 channels in insects. A characteristic aspartate (Asp) that functions as the proton selectivity filter (SF) and the RxWRxxR voltage‐sensor motif are conserved structural elements in H_V1 channels. By analysing Transcriptome Shotgun Assembly (TSA) databases, we found 33 polyneopteran species meeting these structural requirements. Unexpectedly, an unusual natural variation Asp to glutamate (Glu) at SF was found in Phasmatodea and Mantophasmatodea. Additionally, we analysed the expression and function of H_V1 in the phasmatodean stick insect Extatosoma tiaratum (Et). EtH_V1 is strongly expressed in nervous tissue and shows pronounced inward proton conduction. This is the first study of a natural occurring Glu within the SF of a functional H_V1 and might be instrumental in uncovering the physiological function of H_V1 in insects.

The Genome of the "Sea Vomit" Didemnum vexillum

Tunicates are the sister group of vertebrates and thus occupy a key position for investigations into vertebrate innovations as well as into the consequences of the vertebrate-specific genome duplications. Nevertheless, tunicate genomes have not been studied extensively in the past, and comparative studies of tunicate genomes have remained scarce. The carpet sea squirt Didemnum vexillum, commonly known as "sea vomit", is a colonial tunicate considered an invasive species with substantial ecological and economical risk. We report the assembly of the D. vexillum genome using a hybrid approach that combines 28.5 Gb Illumina and 12.35 Gb of PacBio data. The new hybrid scaffolded assembly has a total size of 517.55 Mb that increases contig length about eightfold compared to previous, Illumina-only assembly. As a consequence of an unusually high genetic diversity of the colonies and the moderate length of the PacBio reads, presumably caused by the unusually acidic milieu of the tunic, the assembly is highly fragmented (L50 = 25,284, N50 = 6539). It is sufficient, however, for comprehensive annotations of both protein-coding genes and non-coding RNAs. Despite its shortcomings, the draft assembly of the "sea vomit" genome provides a valuable resource for comparative tunicate genomics and for the study of the specific properties of colonial ascidians.

A review of computational tools for generating metagenome-assembled genomes from metagenomic sequencing data

Metagenomic sequencing provides a culture-independent avenue to investigate the complex microbial communities by constructing metagenome-assembled genomes (MAGs). A MAG represents a microbial genome by a group of sequences from genome assembly with similar characteristics. It enables us to identify novel species and understand their potential functions in a dynamic ecosystem. Many computational tools have been developed to construct and annotate MAGs from metagenomic sequencing, however, there is a prominent gap to comprehensively introduce their background and practical performance. In this paper, we have thoroughly investigated the computational tools designed for both upstream and downstream analyses, including metagenome assembly, metagenome binning, gene prediction, functional annotation, taxonomic classification, and profiling. We have categorized the commonly used tools into unique groups based on their functional background and introduced the underlying core algorithms and associated information to demonstrate a comparative outlook. Furthermore, we have emphasized the computational requisition and offered guidance to the users to select the most efficient tools. Finally, we have indicated current limitations, potential solutions, and future perspectives for further improving the tools of MAG construction and annotation. We believe that our work provides a consolidated resource for the current stage of MAG studies and shed light on the future development of more effective MAG analysis tools on metagenomic sequencing.

Nearly 30 Years of Animal Models to Study Amyotrophic Lateral Sclerosis: A Historical Overview and Future Perspectives

Amyotrophic lateral sclerosis (ALS) is a fatal, multigenic, multifactorial, and non-cell autonomous neurodegenerative disease characterized by upper and lower motor neuron loss. Several genetic mutations lead to ALS development and many emerging gene mutations have been discovered in recent years. Over the decades since 1990, several animal models have been generated to study ALS pathology including both vertebrates and invertebrates such as yeast, worms, flies, zebrafish, mice, rats, guinea pigs, dogs, and non-human primates. Although these models show different peculiarities, they are all useful and complementary to dissect the pathological mechanisms at the basis of motor neuron degeneration and ALS progression, thus contributing to the development of new promising therapeutics. In this review, we describe the up to date and available ALS genetic animal models, classified by the different genetic mutations and divided per species, pointing out their features in modeling, the onset and progression of the pathology, as well as their specific pathological hallmarks. Moreover, we highlight similarities, differences, advantages, and limitations, aimed at helping the researcher to select the most appropriate experimental animal model, when designing a preclinical ALS study.

Sex without crossing over in the yeast Saccharomycodes ludwigii

BACKGROUND

Intermixing of genomes through meiotic reassortment and recombination of homologous chromosomes is a unifying theme of sexual reproduction in eukaryotic organisms and is considered crucial for their adaptive evolution. Previous studies of the budding yeast species Saccharomycodes ludwigii suggested that meiotic crossing over might be absent from its sexual life cycle, which is predominated by fertilization within the meiotic tetrad.

RESULTS

We demonstrate that recombination is extremely suppressed during meiosis in Sd. ludwigii. DNA double-strand break formation by the conserved transesterase Spo11, processing and repair involving interhomolog interactions are required for normal meiosis but do not lead to crossing over. Although the species has retained an intact meiotic gene repertoire, genetic and population analyses suggest the exceptionally rare occurrence of meiotic crossovers in its genome. A strong AT bias of spontaneous mutations and the absence of recombination are likely responsible for its unusually low genomic GC level.

CONCLUSIONS

Sd. ludwigii has followed a unique evolutionary trajectory that possibly derives fitness benefits from the combination of frequent mating between products of the same meiotic event with the extreme suppression of meiotic recombination. This life style ensures preservation of heterozygosity throughout its genome and may enable the species to adapt to its environment and survive with only minimal levels of rare meiotic recombination. We propose Sd. ludwigii as an excellent natural forum for the study of genome evolution and recombination rates.

Reconstruction of evolving gene variants and fitness from short sequencing reads

Directed evolution can generate proteins with tailor-made activities. However, full-length genotypes, their frequencies and fitnesses are difficult to measure for evolving gene-length biomolecules using most high-throughput DNA sequencing methods, as short read lengths can lose mutation linkages in haplotypes. Here we present Evoracle, a machine learning method that accurately reconstructs full-length genotypes (R² = 0.94) and fitness using short-read data from directed evolution experiments, with substantial improvements over related methods. We validate Evoracle on phage-assisted continuous evolution (PACE) and phage-assisted non-continuous evolution (PANCE) of adenine base editors and OrthoRep evolution of drug-resistant enzymes. Evoracle retains strong performance (R² = 0.86) on data with complete linkage loss between neighboring nucleotides and large measurement noise, such as pooled Sanger sequencing data (~US$10 per timepoint), and broadens the accessibility of training machine learning models on gene variant fitnesses. Evoracle can also identify high-fitness variants, including low-frequency 'rising stars', well before they are identifiable from consensus mutations.

Trends in biological data integration for the selection of enzymes and transcription factors related to cellulose and hemicellulose degradation in fungi

Fungi are key players in biotechnological applications. Although several studies focusing on fungal diversity and genetics have been performed, many details of fungal biology remain unknown, including how cellulolytic enzymes are modulated within these organisms to allow changes in main plant cell wall compounds, cellulose and hemicellulose, and subsequent biomass conversion. With the advent and consolidation of DNA/RNA sequencing technology, different types of information can be generated at the genomic, structural and functional levels, including the gene expression profiles and regulatory mechanisms of these organisms, during degradation-induced conditions. This increase in data generation made rapid computational development necessary to deal with the large amounts of data generated. In this context, the origination of bioinformatics, a hybrid science integrating biological data with various techniques for information storage, distribution and analysis, was a fundamental step toward the current state-of-the-art in the postgenomic era. The possibility of integrating biological big data has facilitated exciting discoveries, including identifying novel mechanisms and more efficient enzymes, increasing yields, reducing costs and expanding opportunities in the bioprocess field. In this review, we summarize the current status and trends of the integration of different types of biological data through bioinformatics approaches for biological data analysis and enzyme selection.

Genomic Studies of White-Rot Fungus Cerrena unicolor SP02 Provide Insights into Food Safety Value-Added Utilization of Non-Food Lignocellulosic Biomass

Cerrena unicolor is an ecologically and biotechnologically important wood-degrading basidiomycete with high lignocellulose degrading ability. Biological and genetic investigations are limited in the Cerrena genus and, thus, hinder genetic modification and commercial use. The aim of the present study was to provide a global understanding through genomic and experimental research about lignocellulosic biomass utilization by Cerrena unicolor. In this study, we reported the genome sequence of C. unicolor SP02 by using the Illumina and PacBio 20 platforms to obtain trustworthy assembly and annotation. This is the combinational 2nd and 3rd genome sequencing and assembly of C. unicolor species. The generated genome was 42.79 Mb in size with an N50 contig size of 2.48 Mb, a G + C content of 47.43%, and encoding of 12,277 predicted genes. The genes encoding various lignocellulolytic enzymes including laccase, lignin peroxidase, manganese peroxidase, cytochromes P450, cellulase, xylanase, α-amylase, and pectinase involved in the degradation of lignin, cellulose, xylan, starch, pectin, and chitin that showed the C. unicolor SP02 potentially have a wide range of applications in lignocellulosic biomass conversion. Genome-scale metabolic analysis opened up a valuable resource for a better understanding of carbohydrate-active enzymes (CAZymes) and oxidoreductases that provide insights into the genetic basis and molecular mechanisms for lignocellulosic degradation. The C. unicolor SP02 model can be used for the development of efficient microbial cell factories in lignocellulosic industries. The understanding of the genetic material of C. unicolor SP02 coding for the lignocellulolytic enzymes will significantly benefit us in genetic manipulation, site-directed mutagenesis, and industrial biotechnology.

Cas9-mediated genome editing reveals a significant contribution of calcium signaling pathways to anhydrobiosis in Pv11 cells

Pv11 is an insect cell line established from the midge Polypedilum vanderplanki, whose larval form exhibits an extreme desiccation tolerance known as anhydrobiosis. Pv11 itself is also capable of anhydrobiosis, which is induced by trehalose treatment. Here we report the successful construction of a genome editing system for Pv11 cells and its application to the identification of signaling pathways involved in anhydrobiosis. Using the Cas9-mediated gene knock-in system, we established Pv11 cells that stably expressed GCaMP3 to monitor intracellular Ca2+ mobilization. Intriguingly, trehalose treatment evoked a transient increase in cytosolic Ca2+ concentration, and further experiments revealed that the calmodulin-calcineurin-NFAT pathway contributes to tolerance of trehalose treatment as well as desiccation tolerance, while the calmodulin-calmodulin kinase-CREB pathway conferred only desiccation tolerance on Pv11 cells. Thus, our results show a critical contribution of the trehalose-induced Ca2+ surge to anhydrobiosis and demonstrate temporally different roles for each signaling pathway.

Transporters and tubule crystals in the insect Malpighian tubule

The insect renal (Malpighian) tubules are functionally homologous to the mammalian kidney. Accumulating evidence indicates that renal tubule crystals form in a manner similar to mammalian kidney stones. In Drosophila melanogaster, crystals can be induced by diet, toxic substances, or genetic mutations that reflect circumstances influencing or eliciting kidney stones in mammals. Incredibly, many mammalian proteins have distinct homologs in Drosophila, and the function of most homologs have been demonstrated to recapitulate their mammalian and human counterparts. Here, we discuss the present literature establishing Drosophila as a nephrolithiasis model. This insect model may be used to investigate and understand the etiology of kidney stone diseases, especially with regard to calcium oxalate, calcium phosphate and xanthine or urate crystallization.

Chromatin Organization and Function in Drosophila

Eukaryotic genomes are packaged into high-order chromatin structures organized in discrete territories inside the cell nucleus, which is surrounded by the nuclear envelope acting as a barrier. This chromatin organization is complex and dynamic and, thus, determining the spatial and temporal distribution and folding of chromosomes within the nucleus is critical for understanding the role of chromatin topology in genome function. Primarily focusing on the regulation of gene expression, we review here how the genome of Drosophila melanogaster is organized into the cell nucleus, from small scale histone–DNA interactions to chromosome and lamina interactions in the nuclear space.

Complete Genome Sequence of Lactobacillus helveticus JCM 1004, an Aminopeptidase-Producing Lactic Acid Bacterium

We report the complete genome sequence of Lactobacillus helveticus JCM 1004, an aminopeptidase-producing lactic acid bacterium. The genome consists of a circular chromosome which comprises 2,261,280 bp, with a G+C content of 37.56%. The genome was predicted to harbor 13 rRNA genes, 64 tRNA genes, and 2,462 protein-coding sequences.

New approaches for metagenome assembly with short reads

In recent years, the use of longer range read data combined with advances in assembly algorithms has stimulated big improvements in the contiguity and quality of genome assemblies. However, these advances have not directly transferred to metagenomic data sets, as assumptions made by the single genome assembly algorithms do not apply when assembling multiple genomes at varying levels of abundance. The development of dedicated assemblers for metagenomic data was a relatively late innovation and for many years, researchers had to make do using tools designed for single genomes. This has changed in the last few years and we have seen the emergence of a new type of tool built using different principles. In this review, we describe the challenges inherent in metagenomic assemblies and compare the different approaches taken by these novel assembly tools.

An ancient interleukin-16-like molecule regulates hemocyte proliferation via integrin β1 in invertebrates

Interleukins (ILs) are cytokines with crucial functions in innate and adaptive immunity. IL genes are only found in vertebrates, except for IL-16, which has been cloned in some arthropod species. However, the function of this gene in invertebrates is unknown. In the present study, an IL-16-like gene (EsIL-16) was identified from the Chinese mitten crab Eriocheir sinensis. EsIL-16 was predicted to encode a precursor (proEsIL-16) that shares similarities with pro-IL-16 proteins from insects and vertebrates. We show that caspase-3 processes proEsIL-16 into an approximately 144-kDa N-terminal prodomain with nuclear import activity and an approximately 34-kDa mature peptide that might be secreted into the extracellular region. EsIL-16 mRNA could be detected in all analyzed tissues and was significantly upregulated after immune challenge both in vitro and in vivo. T7 phage display library screening suggested potential binding activity between EsIL-16 and integrin, which was confirmed by coimmunoprecipitation assay. Interestingly, EsIL-16 promoted cell proliferation via integrin β1 in primary cultured crab hemocytes and Drosophila S2 cells. Furthermore, the interaction between EsIL-16 and integrin β1 was necessary to efficiently protect the host from bacterial infection. To our knowledge, this study revealed integrin β1 as a receptor for IL-16 and the function of this interaction in hemocyte proliferation in invertebrates for the first time. These results provide new insights into the regulation of innate immune responses in invertebrates and shed the light on the evolution of ILs within the animal kingdom.

Genomic architecture constrained placental mammal X Chromosome evolution

Susumu Ohno proposed that the gene content of the mammalian X Chromosome should remain highly conserved due to dosage compensation. X Chromosome linkage (gene order) conservation is widespread in placental mammals but does not fall within the scope of Ohno's prediction and may be an indirect result of selection on gene content or selection against rearrangements that might disrupt X-Chromosome inactivation (XCI). Previous comparisons between the human and mouse X Chromosome sequences have suggested that although single-copy X Chromosome genes are conserved between species, most ampliconic genes were independently acquired. To better understand the evolutionary and functional constraints on X-linked gene content and linkage conservation in placental mammals, we aligned a new, high-quality, long-read X Chromosome reference assembly from the domestic cat (incorporating 19.3 Mb of targeted BAC clone sequence) to the pig, human, and mouse assemblies. A comprehensive analysis of annotated X-linked orthologs in public databases demonstrated that the majority of ampliconic gene families were present on the ancestral placental X Chromosome. We generated a domestic cat Hi-C contact map from an F1 domestic cat/Asian leopard cat hybrid and demonstrated the formation of the bipartite structure found in primate and rodent inactivated X Chromosomes. Conservation of gene order and recombination patterns is attributable to strong selective constraints on three-dimensional genomic architecture necessary for superloop formation. Species with rearranged X Chromosomes retain the ancestral order and relative spacing of loci critical for superloop formation during XCI, with compensatory inversions evolving to maintain these long-range physical interactions.

Empirical evaluation of methods for de novo genome assembly

Technologies for next-generation sequencing (NGS) have stimulated an exponential rise in high-throughput sequencing projects and resulted in the development of new read-assembly algorithms. A drastic reduction in the costs of generating short reads on the genomes of new organisms is attributable to recent advances in NGS technologies such as Ion Torrent, Illumina, and PacBio. Genome research has led to the creation of high-quality reference genomes for several organisms, and de novo assembly is a key initiative that has facilitated gene discovery and other studies. More powerful analytical algorithms are needed to work on the increasing amount of sequence data. We make a thorough comparison of the de novo assembly algorithms to allow new users to clearly understand the assembly algorithms: overlap-layout-consensus and de-Bruijn-graph, string-graph based assembly, and hybrid approach. We also address the computational efficacy of each algorithm's performance, challenges faced by the assem- bly tools used, and the impact of repeats. Our results compare the relative performance of the different assemblers and other related assembly differences with and without the reference genome. We hope that this analysis will contribute to further the application of de novo sequences and help the future growth of assembly algorithms.

Applications of CRISPR-Cas9 as an Advanced Genome Editing System in Life Sciences

Targeted nucleases are powerful genomic tools to precisely change the target genome of living cells, controlling functional genes with high exactness. The clustered regularly interspaced short palindromic repeats associated protein 9 (CRISPR-Cas9) genome editing system has been identified as one of the most useful biological tools in genetic engineering that is taken from adaptive immune strategies for bacteria. In recent years, this system has made significant progress and it has been widely used in genome editing to create gene knock-ins, knock-outs, and point mutations. This paper summarizes the application of this system in various biological sciences, including medicine, plant science, and animal breeding.

Analysis of the Taxonomy and Pathogenic Factors of Pectobacterium aroidearum L6 Using Whole-Genome Sequencing and Comparative Genomics

Soft rot pectobacteria are devastating plant pathogens with a global distribution and a broad host range. Pectobacterium aroidearum L6, previously isolated from leaves of Syngonium podophyllum, is a pectolytic bacterial pathogen that causes typical soft rot on S. podophyllum. There is a shortage for genome data of P. aroidearum, which seriously hinders research on classification and pathogenesis of Pectobacterium. We present here the complete genome sequence of P. aroidearum L6. The L6 strain carries a single 4,995,896-bp chromosome with 53.10% G + C content and harbors 4,306 predicted protein-coding genes. We estimated in silico DNA-DNA hybridization and average nucleotide identity values in combination with the whole-genome-based phylogeny from 19 Pectobacterium strains including P. aroidearum L6. The results showed that L6 and PC1 formed a population distinct from other populations of the Pectobacterium genus. Phylogenetic analysis based on 16S rRNA and genome sequences showed a close evolutionary relationship among Pectobacterium species. Overall, evolutionary analysis showed that L6 was in the same branch with PC1. In comparison with 18 Pectobacterium spp. reference pathogens, strain L6 had 2,712 gene families, among which 1,632 gene families were identified as orthologous to those strains, as well as 1 putative unique gene family. We discovered 478 genes, 10.4% of the total of predicted genes, that were potentially related to pathogenesis using the Virulence Factors of Pathogenic Bacteria database. A total of 25 genes were related to toxins, 35 encoded plant cell-wall degrading enzymes, and 122 were involved in secretion systems. This study provides a foundation for a better understanding of the genomic structure of P. aroidearum and particularly offers information for the discovery of potential pathogenic factors and the development of more effective strategies against this pathogen.

An anchored chromosome-scale genome assembly of spinach improves annotation and reveals extensive gene rearrangements in euasterids

Spinach (Spinacia oleracea L.) is a member of the Caryophyllales family, a basal eudicot asterid that consists of sugar beet (Beta vulgaris L. subsp. vulgaris), quinoa (Chenopodium quinoa Willd.), and amaranth (Amaranthus hypochondriacus L.). With the introduction of baby leaf types, spinach has become a staple food in many homes. Production issues focus on yield, nitrogen-use efficiency and resistance to downy mildew (Peronospora effusa). Although genomes are available for the above species, a chromosome-level assembly exists only for quinoa, allowing for proper annotation and structural analyses to enhance crop improvement. We independently assembled and annotated genomes of the cultivar Viroflay using short-read strategy (Illumina) and long-read strategies (Pacific Biosciences) to develop a chromosome-level, genetically anchored assembly for spinach. Scaffold N50 for the Illumina assembly was 389 kb, whereas that for Pacific BioSciences was 4.43 Mb, representing 911 Mb (93% of the genome) in 221 scaffolds, 80% of which are anchored and oriented on a sequence-based genetic map, also described within this work. The two assemblies were 99.5% collinear. Independent annotation of the two assemblies with the same comprehensive transcriptome dataset show that the quality of the assembly directly affects the annotation with significantly more genes predicted (26,862 vs. 34,877) in the long-read assembly. Analysis of resistance genes confirms a bias in resistant gene motifs more typical of monocots. Evolutionary analysis indicates that Spinacia is a paleohexaploid with a whole-genome triplication followed by extensive gene rearrangements identified in this work. Diversity analysis of 75 lines indicate that variation in genes is ample for hypothesis-driven, genomic-assisted breeding enabled by this work.

Conservation genomics reveals pesticide and pathogen exposure in the declining bumble bee Bombus terricola

In recent years, many pollinators have experienced large population declines, which threaten food security and the stability of natural ecosystems. Bumble bees are particularly important because their ability to “buzz” pollinate and tolerate cooler conditions make them critical pollinators for certain plants and regions. Here, we apply a conservation genomics approach to study the vulnerable Bombus terricola. We sequenced RNA from 30 worker abdomens, 18 of which were collected from agricultural sites and 12 of which were collected from nonagricultural sites. We found transcriptional signatures associated with exposure to insecticides, with gene expression patterns suggesting that bumble bees were exposed to neonicotinoids and/or fipronil—two compounds known to negatively impact bees. We also found transcriptional signatures associated with pathogen infections. In addition to the transcriptomic analysis, we carried out a metatranscriptomic analysis and detected five pathogens in the abdomens of workers, three of which are common in managed honey bee and bumble bee colonies. Our conservation genomics study provides functional support for the role of pesticides and pathogen spillover in the decline of B. terricola. We demonstrate that conservation genomics is an invaluable tool which allows researchers to quantify the effects of multiple stressors that impact pollinator populations in the wild.

Genome assembly using quantum and quantum-inspired annealing

Recent advances in DNA sequencing open prospects to make whole-genome analysis rapid and reliable, which is promising for various applications including personalized medicine. However, existing techniques for de novo genome assembly, which is used for the analysis of genomic rearrangements, chromosome phasing, and reconstructing genomes without a reference, require solving tasks of high computational complexity. Here we demonstrate a method for solving genome assembly tasks with the use of quantum and quantum-inspired optimization techniques. Within this method, we present experimental results on genome assembly using quantum annealers both for simulated data and the [Formula: see text]X 174 bacteriophage. Our results pave a way for a significant increase in the efficiency of solving bioinformatics problems with the use of quantum computing technologies and, in particular, quantum annealing might be an effective method. We expect that the new generation of quantum annealing devices would outperform existing techniques for de novo genome assembly. To the best of our knowledge, this is the first experimental study of de novo genome assembly problems both for real and synthetic data on quantum annealing devices and quantum-inspired techniques.

Levenshtein Distance, Sequence Comparison and Biological Database Search

Levenshtein edit distance has played a central role-both past and present-in sequence alignment in particular and biological database similarity search in general. We start our review with a history of dynamic programming algorithms for computing Levenshtein distance and sequence alignments. Following, we describe how those algorithms led to heuristics employed in the most widely used software in bioinformatics, BLAST, a program to search DNA and protein databases for evolutionarily relevant similarities. More recently, the advent of modern genomic sequencing and the volume of data it generates has resulted in a return to the problem of local alignment. We conclude with how the mathematical formulation of Levenshtein distance as a metric made possible additional optimizations to similarity search in biological contexts. These modern optimizations are built around the low metric entropy and fractional dimensionality of biological databases, enabling orders of magnitude acceleration of biological similarity search.

A method for achieving complete microbial genomes and improving bins from metagenomics data

Metagenomics facilitates the study of the genetic information from uncultured microbes and complex microbial communities. Assembling complete genomes from metagenomics data is difficult because most samples have high organismal complexity and strain diversity. Some studies have attempted to extract complete bacterial, archaeal, and viral genomes and often focus on species with circular genomes so they can help confirm completeness with circularity. However, less than 100 circularized bacterial and archaeal genomes have been assembled and published from metagenomics data despite the thousands of datasets that are available. Circularized genomes are important for (1) building a reference collection as scaffolds for future assemblies, (2) providing complete gene content of a genome, (3) confirming little or no contamination of a genome, (4) studying the genomic context and synteny of genes, and (5) linking protein coding genes to ribosomal RNA genes to aid metabolic inference in 16S rRNA gene sequencing studies. We developed a semi-automated method called Jorg to help circularize small bacterial, archaeal, and viral genomes using iterative assembly, binning, and read mapping. In addition, this method exposes potential misassemblies from k-mer based assemblies. We chose species of the Candidate Phyla Radiation (CPR) to focus our initial efforts because they have small genomes and are only known to have one ribosomal RNA operon. In addition to 34 circular CPR genomes, we present one circular Margulisbacteria genome, one circular Chloroflexi genome, and two circular megaphage genomes from 19 public and published datasets. We demonstrate findings that would likely be difficult without circularizing genomes, including that ribosomal genes are likely not operonic in the majority of CPR, and that some CPR harbor diverged forms of RNase P RNA. Code and a tutorial for this method is available at https://github.com/lmlui/Jorg and is available on the DOE Systems Biology KnowledgeBase as a beta app.

Complete Genome Sequence of Anaerostipes caccae Strain L1-92T, a Butyrate-Producing Bacterium Isolated from Human Feces

Anaerostipes caccae strain L1-92^T is a well-known butyrate-producing bacterium that has been isolated from human feces. In this announcement, we present the complete genome sequence of A. caccae strain L1-92^T, which comprises 3,590,719 bp with a G+C content of 44.30%. The genome harbors 3,369 predicted protein-coding genes.

Characterization of Proteus vulgaris Strain P3M, a Foodborne Multidrug-Resistant Bacterium Isolated from Penaeus vannamei in China

Proteus vulgaris is an important foodborne opportunistic pathogen, both environmentally and clinically. The use of appropriate antibiotics has significant therapeutic effects, but has led to the emergence and spread of drug-resistant strains. In this study, a P. vulgaris strain, designated "P3M," was isolated from Penaeus vannamei in Tianjin, China. The whole genome of P3M was sequenced, generating detailed information, including the key genes involved in important metabolic pathways and their physiological functions. A total of 218 antibiotic resistance genes (ARGs) were predicted in the genome. The determination of various minimum inhibitory concentrations indicated that P3M is a multidrug-resistant (MDR) bacterium, with significant resistance to 16 antibiotics in seven categories. Determination of fractional inhibitory concentration index showed that the combination of ciprofloxacin plus tetracycline exhibited synergistic antimicrobial activity. Bioinformatics and phylogenetic analyses detected the presence of two two-component systems that mediate multidrug resistance and several mobile genetic elements involved in the horizontal transfer of ARGs in P3M. P. vulgaris strains represent a serious challenge to clinicians and infection control teams for its ubiquity worldwide and close relevance with human life. To the best of our knowledge, we report the first isolation and characterization of an important foodborne MDR P. vulgaris strain, and this study will provide necessary theoretical basis for the selection and clinical use of the appropriate antibiotics.

Benchmarking of next and third generation sequencing technologies and their associated algorithms for de novo genome assembly

Genome assemblers are computational tools for de novo genome assembly, based on a plenitude of primary sequencing data. The quality of genome assemblies is estimated by their contiguity and the occurrences of misassemblies (duplications, deletions, translocations or inversions). The rapid development of sequencing technologies has enabled the rise of novel de novo genome assembly strategies. The ultimate goal of such strategies is to utilise the features of each sequencing platform in order to address the existing weaknesses of each sequencing type and compose a complete and correct genome map. In the present study, the hybrid strategy, which is based on Illumina short paired‑end reads and Nanopore long reads, was benchmarked using MaSuRCA and Wengan assemblers. Moreover, the long‑read assembly strategy, which is based on Nanopore reads, was benchmarked using Canu or PacBio HiFi reads were benchmarked using Hifiasm and HiCanu. The assemblies were performed on a computational cluster with limited computational resources. Their outputs were evaluated in terms of accuracy and computational performance. PacBio HiFi assembly strategy outperforms the other ones, while Hi‑C scaffolding, which is based on chromatin 3D structure, is required in order to increase continuity, accuracy and completeness when large and complex genomes, such as the human one, are assembled. The use of Hi‑C data is also necessary while using the hybrid assembly strategy. The results revealed that HiFi sequencing enabled the rise of novel algorithms which require less genome coverage than that of the other strategies making the assembly a less computationally demanding task. Taken together, these developments may lead to the democratisation of genome assembly projects which are now approachable by smaller labs with limited technical and financial resources.

Efficient hybrid de novo assembly of human genomes with WENGAN

Generating accurate genome assemblies of large, repeat-rich human genomes has proved difficult using only long, error-prone reads, and most human genomes assembled from long reads add accurate short reads to polish the consensus sequence. Here we report an algorithm for hybrid assembly, WENGAN, that provides very high quality at low computational cost. We demonstrate de novo assembly of four human genomes using a combination of sequencing data generated on ONT PromethION, PacBio Sequel, Illumina and MGI technology. WENGAN implements efficient algorithms to improve assembly contiguity as well as consensus quality. The resulting genome assemblies have high contiguity (contig NG50: 17.24-80.64 Mb), few assembly errors (contig NGA50: 11.8-59.59 Mb), good consensus quality (QV: 27.84-42.88) and high gene completeness (BUSCO complete: 94.6-95.2%), while consuming low computational resources (CPU hours: 187-1,200). In particular, the WENGAN assembly of the haploid CHM13 sample achieved a contig NG50 of 80.64 Mb (NGA50: 59.59 Mb), which surpasses the contiguity of the current human reference genome (GRCh38 contig NG50: 57.88 Mb).

Complete genome sequence of the deep South China Sea-derived Streptomyces niveus SCSIO 3406, the producer of cytotoxic and antibacterial marfuraquinocins

Streptomyces niveus SCSIO 3406 was isolated from a sediment sample collected from South China Sea at a depth of 3536 m. Four new sesquiterpenoid naphthoquinones, marfuraquinocins A-D, and two new geranylated phenazines, i. e. phenaziterpenes A and B, were isolated from the fermentation broth of the strain. Here, we present its genome sequence, which contains 7,990,492 bp with a G+C content of 70.46% and harbors 7088 protein-encoding genes. The genome sequence analysis revealed the presence of a 28,787 bp gene cluster encoding for 24 open reading frames including 1,3,6,8-tetrahydroxynaphthalene synthase and monooxygenase, seven phenazine biosynthesis proteins, two prenyltransferases and a squalene-hopene cyclase. These genes are known to be necessary for the biosynthesis of both marfuraquinocins and phenaziterpenes. Outside the gene cluster (and scattered around the genome), there are seven genes belonging to the methylerythritol phosphate pathway for the biosynthesis of the essential primary metabolite, isopentenyl diphosphate, as well as six geranyl diphosphate/farnesyl diphosphate synthase genes. The strain S. niveus SCSIO 3406 showed type I PKS, type III PKS and nonribosomal peptide synthetase cluster. The sequence will provide the genetic basis for better understanding of biosynthesis mechanism of the above mentioned six compounds and for the construction of improved strain for the industrial production of antimicrobial agents.

The mitochondrial genome of Telenomus remus (Hymenoptera: Platygastridae)

Telenomus remus Nixon, 1937 is an important parasitoid of lepidopterans. We sequenced the mitochondrial genome of T. remus, 15,500 bp in size, and possessed all 37 typical mitochondrial genes. A few tRNAs show gene arrangements compared with the ancestral gene order, mainly involving in the four tRNA clusters (E-C-Y-Q-I-A, D-K, N-F-S1-R, and M-V). The nucleotide sequences of 13 protein-coding genes of this sequence and another seven species from Platygastridae were used for phylogenetic analysis by MrBayes, with two species from Cynipoidea as an outgroup. The topology demonstrated that T. remus was most closely related to Telenomus sp.

Complete genome sequencing of Bacillus sp. TK-2, analysis of its cold evolution adaptability

To date, a large number of Bacillus species from different sources have been identified. However, there are few investigations on genome information and evolutionary insights of Bacillus species from cold environments. Bacillus sp. TK-2, isolated from the soil of Changbai Mountain, is a gram-positive bacterium with cold adaptation characteristics. In this study, we present the annotated complete genome sequence of Bacillus sp. TK-2. The genome comprised 5,286,177 bp with a GC content of 35.88%, 5293 protein-encoding genes, 32 rRNA, and 77 tRNA. Numerous genes related to cold adaptation were detected in the genome of Bacillus sp. TK-2, mainly involving in energy supply, regulation of cell membrane fluidity, antioxidant, and molecular chaperones. In addition, the strain TK-2 classified in the Bacillus groups was distributed on a terminal branch with Bacillus cereus A1 by Blastn and phylogenetic analysis in NCBI database. Complete genome sequences of the strain TK-2 and Bacillus cereus A1 were compared by the online tool "Average Nucleotide Identity", showing that the average nucleotide identity of these two strains was 98.26%. In parallel, A comparative analysis of the genomes of both Bacillus sp. TK-2 and Bacillus cereus A1 was conducted. Through the analysis of core and specific genes with cd-hit, it was found that the two strains had 5691 pan gene, 4524 core gene, and 1167 specific gene clusters. Among the 624 specific gene clusters of Bacillus sp. TK-2, some cold tolerance genes were detected, which implied the unique adaptability of Bacillus sp. TK-2 in long-term low temperature environments. Importantly, enzyme-encoding genes related to the degradation of polysaccharides such as cellulose and hemicellulose were detected in the 477 CAZyme genes of this genome. This work on sequencing and bioinformatics analysis of the complete sequence of Bacillus sp. TK-2 promote the application and in-depth research of low-temperature biotechnology.

Biogeographic and Evolutionary Patterns of Trace Element Utilization in Marine Microbial World

Trace elements are required by all organisms, which are key components of many enzymes catalyzing important biological reactions. Many trace element-dependent proteins have been characterized; however, little is known about their occurrence in microbial communities in diverse environments, especially the global marine ecosystem. Moreover, the relationships between trace element utilization and different types of environmental stressors are unclear. In this study, we used metagenomic data from the Global Ocean Sampling expedition project to identify the biogeographic distribution of genes encoding trace element-dependent proteins (for copper, molybdenum, cobalt, nickel, and selenium) in a variety of marine and non-marine aquatic samples. More than 56,000 metalloprotein and selenoprotein genes corresponding to nearly 100 families were predicted, becoming the largest dataset of marine metalloprotein and selenoprotein genes reported to date. In addition, samples with enriched or depleted metalloprotein/selenoprotein genes were identified, suggesting an active or inactive usage of these micronutrients in various sites. Further analysis of interactions among the elements showed significant correlations between some of them, especially those between nickel and selenium/copper. Finally, investigation of the relationships between environmental conditions and metalloprotein/selenoprotein families revealed that many environmental factors might contribute to the evolution of different metalloprotein and/or selenoprotein genes in the marine microbial world. Our data provide new insights into the utilization and biological roles of these trace elements in extant marine microbes, and might also be helpful for the understanding of how these organisms have adapted to their local environments.

An Algorithm for Gene Fragment Reconstruction

Gene sequencing technology has been playing an important role in many aspects, such as life science, disease medicine and health medicine, particularly in the extremely tough process of fighting against 2019-novel coronavirus. Drawing DNA restriction map is a particularly important technology in genetic biology. The simplified partial digestion method (SPDP), a biological method, has been widely used to cut DNA molecules into DNA fragments and obtain the biological information of each fragment. In this work, we propose an algorithm based on 0–1 planning for the location of restriction sites on a DNA molecule, which is able to solve the problem of DNA fragment reconstruction just based on data of fragments’ length. Two specific examples are presented in detail. Furthermore, based on 1000 groups of original DNA sequences randomly generated, we define the coincidence rate and unique coincidence rate between the reconstructed DNA sequence and the original DNA sequence, and then analyze separately the effect of the number of fragments and the maximum length of DNA fragments on the coincidence rate and unique coincidence rate as defined. The effectiveness of the algorithm is proved. Besides, based on the existing optimization solution obtained, we simulate and discuss the influence of the error by computation method. It turns out that the error of position of one restriction site does not affect other restriction sites and errors of most restriction sites may lead to the failure of sequence reconstruction. Matlab 7.1 program is used to solve feasible solutions of the location of restriction sites, derive DNA fragment sequence and carry out the statistical analysis and error analysis. This paper focuses on basic computer algorithm implementation of rearrangement and sequencing rather than biochemical technology. The innovative application of the mathematical idea of 0–1 planning to DNA sequence mapping construction, to a certain extent, greatly simplifies the difficulty and complexity of calculation and accelerates the process of 'jigsaw' of DNA fragments.

Critical evaluation of short, long, and hybrid assembly for contextual analysis of antibiotic resistance genes in complex environmental metagenomes

In the fight to limit the global spread of antibiotic resistance, the assembly of environmental metagenomes has the potential to provide rich contextual information (e.g., taxonomic hosts, carriage on mobile genetic elements) about antibiotic resistance genes (ARG) in the environment. However, computational challenges associated with assembly can impact the accuracy of downstream analyses. This work critically evaluates the impact of assembly leveraging short reads, nanopore MinION long-reads, and a combination of the two (hybrid) on ARG contextualization for ten environmental metagenomes using seven prominent assemblers (IDBA-UD, MEGAHIT, Canu, Flye, Opera-MS, metaSpades and HybridSpades). While short-read and hybrid assemblies produced similar patterns of ARG contextualization, raw or assembled long nanopore reads produced distinct patterns. Based on an in-silico spike-in experiment using real and simulated reads, we show that low to intermediate coverage species are more likely to be incorporated into chimeric contigs across all assemblers and sequencing technologies, while more abundant species produce assemblies with a greater frequency of inversions and insertion/deletions (indels). In sum, our analyses support hybrid assembly as a valuable technique for boosting the reliability and accuracy of assembly-based analyses of ARGs and neighboring genes at environmentally-relevant coverages, provided that sufficient short-read sequencing depth is achieved.

Gene Sequence Assembly Algorithm Model Based on the DBG Strategy and Its Application

With the continuous development of sequencing technology, the amount of bioinformatics data has increased geometrically, and the massive amount of bioinformatics data puts forward more stringent requirements for sequence assembly problems. The sequence assembly algorithm based on DBG (De Bruijn graph) strategy is a key algorithm in bioinformatics, which is widely used in the domain of gene sequence assembly. Current research on the domain of sequence assembly always focuses on optimization of specific steps to a specific algorithm and lack of research on domain-level high-abstract algorithm frameworks. To some extent, it leads to the redundancy of the sequence assembly algorithm, and some problems may be caused by the artificial selection algorithm. This paper analyzes the domain of DBGSA and establishes a feature model of this domain. Based on the production programming method, the DBGSA algorithm component is interactively designed. With the support of the PAR platform, the DBGSA algorithm component library is formally implemented, and furthermore, the DBGSA component library is used to assemble the specific algorithm. This research adds domain-level research to the domain of sequence assembly and implements the DBGSA component library, which can assemble specific sequence assembly algorithms, ensuring the efficiency of algorithm development and the reliability of assembly generation algorithms. At the same time, it also provides a valuable reference for solving problems in the domain of sequence assembly.

Genome-scale de novo assembly using ALGA

Motivation

There are very few methods for de novo genome assembly based on the overlap graph approach. It is considered as giving more exact results than the so-called de Bruijn graph approach but in much greater time and of much higher memory usage. It is not uncommon that assembly methods involving the overlap graph model are not able to successfully compute greater datasets, mainly due to memory limitation of a computer. This was the reason for developing in last decades mainly de Bruijn-based assembly methods, fast and fairly accurate. However, the latter methods can fail for longer or more repetitive genomes, as they decompose reads to shorter fragments and lose a part of information. An efficient assembler for processing big datasets and using the overlap graph model is still looked out.

Results

We propose a new genome-scale de novo assembler based on the overlap graph approach, designed for short-read sequencing data. The method, ALGA, incorporates several new ideas resulting in more exact contigs produced in short time. Among these ideas, we have creation of a sparse but quite informative graph, reduction of the graph including a procedure referring to the problem of minimum spanning tree of a local subgraph, and graph traversal connected with simultaneous analysis of contigs stored so far. What is rare in genome assembly, the algorithm is almost parameter-free, with only one optional parameter to be set by a user. ALGA was compared with nine state-of-the-art assemblers in tests on genome-scale sequencing data obtained from real experiments on six organisms, differing in size, coverage, GC content and repetition rate. ALGA produced best results in the sense of overall quality of genome reconstruction, understood as a good balance between genome coverage, accuracy and length of resulting sequences. The algorithm is one of tools involved in processing data in currently realized national project Genomic Map of Poland.

Availability and implementation

ALGA is available at http://alga.put.poznan.pl.

Supplementary information

Supplementary data are available at Bioinformatics online.

Amyotrophic Lateral Sclerosis Genes in Drosophila melanogaster

Amyotrophic lateral sclerosis (ALS) is a devastating adult-onset neurodegenerative disease characterized by the progressive degeneration of upper and lower motoneurons. Most ALS cases are sporadic but approximately 10% of ALS cases are due to inherited mutations in identified genes. ALS-causing mutations were identified in over 30 genes with superoxide dismutase-1 (SOD1), chromosome 9 open reading frame 72 (C9orf72), fused in sarcoma (FUS), and TAR DNA-binding protein (TARDBP, encoding TDP-43) being the most frequent. In the last few decades, Drosophila melanogaster emerged as a versatile model for studying neurodegenerative diseases, including ALS. In this review, we describe the different Drosophila ALS models that have been successfully used to decipher the cellular and molecular pathways associated with SOD1, C9orf72, FUS, and TDP-43. The study of the known fruit fly orthologs of these ALS-related genes yielded significant insights into cellular mechanisms and physiological functions. Moreover, genetic screening in tissue-specific gain-of-function mutants that mimic ALS-associated phenotypes identified disease-modifying genes. Here, we propose a comprehensive review on the Drosophila research focused on four ALS-linked genes that has revealed novel pathogenic mechanisms and identified potential therapeutic targets for future therapy.

Yes, we can use it: a formal test on the accuracy of low-pass nanopore long-read sequencing for mitophylogenomics and barcoding research using the Caribbean spiny lobster Panulirus argus

Background

Whole mitogenomes or short fragments (i.e., 300–700 bp of the cox1 gene) are the markers of choice for revealing within- and among-species genealogies. Protocols for sequencing and assembling mitogenomes include ‘primer walking’ or ‘long PCR’ followed by Sanger sequencing or Illumina short-read low-coverage whole genome (LC-WGS) sequencing with or without prior enrichment of mitochondrial DNA. The aforementioned strategies assemble complete and accurate mitochondrial genomes but are time consuming and/or expensive. In this study, I first tested whether mitogenomes can be sequenced from long-read nanopore sequencing data exclusively. Second, I explored the accuracy of the long-read assembled genomes by comparing them to a ‘gold’ standard reference mitogenome retrieved from the same individual using Illumina sequencing. Third and lastly, I tested if the long-read assemblies are useful for mitophylogenomics and barcoding research. To accomplish these goals, I used the Caribbean spiny lobster Panulirus argus, an ecologically relevant species in shallow water coral reefs and target of the most lucrative fishery in the greater Caribbean region.

Results

LC-WGS using a MinION ONT device and various de-novo and reference-based assembly pipelines retrieved a complete and highly accurate mitogenome for the Caribbean spiny lobster Panulirus argus. Discordance between each of the long-read assemblies and the reference mitogenome was mostly due to indels at the flanks of homopolymer regions. Although not ‘perfect’, phylogenetic analyses using entire mitogenomes or a fragment of the cox1 gene demonstrated that mitogenomes assembled using long reads reliably identify the sequenced specimen as belonging to P. argus and distinguish it from other related species in the same genus, family, and superorder.

Conclusions

This study serves as a proof-of-concept for the future implementation of in-situ surveillance protocols using the MinION to detect mislabeling in P. argus across its supply chain. Mislabeling detection will improve fishery management in this overexploited lobster. This study will additionally aid in decreasing costs for exploring meta-population connectivity in the Caribbean spiny lobster and will aid with the transfer of genomics technology to low-income countries.

Accurate determination of node and arc multiplicities in de bruijn graphs using conditional random fields

Background

De Bruijn graphs are key data structures for the analysis of next-generation sequencing data. They efficiently represent the overlap between reads and hence, also the underlying genome sequence. However, sequencing errors and repeated subsequences render the identification of the true underlying sequence difficult. A key step in this process is the inference of the multiplicities of nodes and arcs in the graph. These multiplicities correspond to the number of times each k-mer (resp. k+1-mer) implied by a node (resp. arc) is present in the genomic sequence. Determining multiplicities thus reveals the repeat structure and presence of sequencing errors. Multiplicities of nodes/arcs in the de Bruijn graph are reflected in their coverage, however, coverage variability and coverage biases render their determination ambiguous. Current methods to determine node/arc multiplicities base their decisions solely on the information in nodes and arcs individually, under-utilising the information present in the sequencing data.

Results

To improve the accuracy with which node and arc multiplicities in a de Bruijn graph are inferred, we developed a conditional random field (CRF) model to efficiently combine the coverage information within each node/arc individually with the information of surrounding nodes and arcs. Multiplicities are thus collectively assigned in a more consistent manner.

Conclusions

We demonstrate that the CRF model yields significant improvements in accuracy and a more robust expectation-maximisation parameter estimation. True k-mers can be distinguished from erroneous k-mers with a higher F₁ score than existing methods. A C++11 implementation is available at https://github.com/biointec/detoxunder the GNU AGPL v3.0 license.

Tissue specific diversification, virulence and immune response to Mycobacterium bovis BCG in a patient with an IFN-γ R1 deficiency

Summary: We characterized Mycobacterium bovis BCG isolates found in lung and brain samples from a previously vaccinated patient with IFNγR1 deficiency. The isolates collected displayed distinct genomic and phenotypic features consistent with host adaptation and associated changes in antibiotic susceptibility and virulence traits. Background: We report a case of a patient with partial recessive IFNγR1 deficiency who developed disseminated BCG infection after neonatal vaccination (BCG-vaccine). Distinct M. bovis BCG-vaccine derived clinical strains were recovered from the patient's lungs and brain. Methods: BCG strains were phenotypically (growth, antibiotic susceptibility, lipid) and genetically (whole genome sequencing) characterized. Mycobacteria cell infection models were used to assess apoptosis, necrosis, cytokine release, autophagy, and JAK-STAT signaling. Results: Clinical isolates BCG-brain and BCG-lung showed distinct Rv0667 rpoB mutations conferring high- and low-level rifampin resistance; the latter displayed clofazimine resistance through Rv0678 gene (MarR-like transcriptional regulator) mutations. BCG-brain and BCG-lung showed mutations in fadA2, fadE5, and mymA operon genes, respectively. Lipid profiles revealed reduced levels of PDIM in BCG-brain and BCG-lung and increased TAGs and Mycolic acid components in BCG-lung, compared to parent BCG-vaccine. In vitro infected cells showed that the BCG-lung induced a higher cytokine release, necrosis, and cell-associated bacterial load effect when compared to BCG-brain; conversely, both strains inhibited apoptosis and altered JAK-STAT signaling. Conclusions: During a chronic-disseminated BCG infection, BCG strains can evolve independently at different sites likely due to particular microenvironment features leading to differential antibiotic resistance, virulence traits resulting in dissimilar responses in different host tissues.

Genomics and lipidomics analysis of the biotechnologically important oleaginous red yeast Rhodotorula glutinis ZHK provides new insights into its lipid and carotenoid metabolism

Background

Rhodotorula glutinis is recognized as a biotechnologically important oleaginous red yeast, which synthesizes numerous meritorious compounds with wide industrial usages. One of the most notable properties of R. glutinis is the formation of intracellular lipid droplets full of carotenoids. However, the basic genomic features that underlie the biosynthesis of these valuable compounds in R. glutinis have not been fully documented. To reveal the biotechnological potential of R. glutinis, the genomics and lipidomics analysis was performed through the Next-Generation Sequencing and HPLC-MS-based metabolomics technologies.

Results

Here, we firstly assemble the genome of R. glutinis ZHK into 21.8 Mb, containing 30 scaffolds and 6774 predicted genes with a N50 length of 14, 66,672 bp and GC content of 67.8%. Genome completeness assessment (BUSCO alignment: 95.3%) indicated the genome assembly with a high-quality features. According to the functional annotation of the genome, we predicted several key genes involved in lipids and carotenoids metabolism as well as certain industrial enzymes biosynthesis. Comparative genomics results suggested that most of orthologous genes have underwent the strong purifying selection within the five Rhodotorula species, especially genes responsible for carotenoids biosynthesis. Furthermore, a total of 982 lipids were identified using the lipidomics approaches, mainly including triacylglycerols, diacylglyceryltrimethylhomo-ser and phosphatidylethanolamine.

Conclusion

Using whole genome shotgun sequencing, we comprehensively analyzed the genome of R. glutinis and predicted several key genes involved in lipids and carotenoids metabolism. By performing comparative genomic analysis, we show that most of the ortholog genes have undergone strong purifying selection within the five Rhodotorula species. Furthermore, we identified 982 lipid species using lipidomic approaches. These results provided valuable resources to further advance biotechnological applications of R .glutinis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07244-z.

The interplay of colour and bioacoustic traits in the differentiation of a Southeast Asian songbird complex

Morphological traits have served generations of biologists as a taxonomic indicator, and have been the main basis for defining and classifying species diversity for centuries. A quantitative integration of behavioural characters, such as vocalizations, in studies on biotic differentiation has arisen more recently, and the relative importance of these different traits in the diversification process remains poorly understood. To provide a framework within which to interpret the evolutionary interplay between morphological and behavioural traits, we generated a draft genome of a cryptic Southeast Asian songbird, the limestone wren-babbler Napothera crispifrons. We resequenced whole genomes of multiple individuals of all three traditional subspecies and of a distinct leucistic population. We demonstrate strong genomic and mitochondrial divergence among all three taxa, pointing to the existence of three species-level lineages. Despite its great phenotypic distinctness, the leucistic population was characterized by shallow genomic differentiation from its neighbour, with only a few localized regions emerging as highly diverged. Quantitative bioacoustic analysis across multiple traits revealed deep differences especially between the two taxa characterized by limited plumage differentiation. Our study demonstrates that differentiation in these furtive songbirds has resulted in a complex mosaic of colour-based and bioacoustic differences among populations. Extreme colour differences can be anchored in few genomic loci and may therefore arise and subside rapidly.

The genomic timeline of cichlid fish diversification across continents

Cichlid fishes are celebrated for their vast taxonomic, phenotypic, and ecological diversity; however, a central aspect of their evolution - the timeline of their diversification - remains contentious. Here, we generate draft genome assemblies of 14 species representing the global cichlid diversity and integrate these into a new phylogenomic hypothesis of cichlid and teleost evolution that we time-calibrate with 58 re-evaluated fossil constraints and a new Bayesian model accounting for fossil-assignment uncertainty. Our results support cichlid diversification long after the breakup of the supercontinent Gondwana and lay the foundation for precise temporal reconstructions of the exceptional continental cichlid adaptive radiations.

Comparative Genomics Analysis Provides New Strategies for Bacteriostatic Ability of Bacillus velezensis HAB-2

Biocontrol formulations prepared from biocontrol bacteria are increasingly applied in sustainable agriculture. Notably, inoculants prepared from Bacillus strains have been proven efficient and environmentally friendly alternatives to chemical bactericides. The bacterium Bacillus velezensis HAB-2 (formerly classified as B. amyloliquefaciens HAB-2) is used as a biological control agent in agricultural fields. In this study, we reported a high-quality genome sequence of HAB-2 using third-generation sequencing technology (PacBio RS II). The 3.89 Mb genome encoded 3,820 predicted genes. Comparative analysis among the genome sequences of reference strains B. velezensis FZB42, B. amyloliquefaciens DSM7 and B. subtilis 168 with the HAB-2 genome revealed obvious differences in the variable part of the genomes, while the core genome shared by FZB42 and HAB-2 was similar (96.14%). However, there were differences in the prophage region among the four strains. The numbers of prophage regions and coding genes in HAB-2 and FZB42 were smaller than the other two strains. The HAB-2 genome showed superior ability to produce secondary metabolites and harbored 13 gene clusters involved in synthesis of antifungal and antibacterial acting secondary metabolites. Furthermore, there were two unique clusters: one cluster which encoded lanthipeptide was involved in mersacidin synthesis and another cluster which encoded ladderane was shown to direct an unknown compound. Multidomain enzymes, such as non-ribosomal peptide synthetase and polyketide synthase, control the biosynthesis of secondary metabolites and rely on 4'-phosphopantetheinyl transferases (PPTases). Key genes lpaH2 and a encoded PPTases in HAB-2 encoded 224 and 120 amino acids, respectively. The genomic features revealed that HAB-2 possesses immense potential to synthesize antimicrobial acting secondary metabolites by regulating and controlling gene clusters. The prophage regions and genes encoding PPTases may provide novel insight for the bacteriostatic mechanism of Bacillus in the biological control of plant diseases.