Velvet: algorithms for de novo short read assembly using de Bruijn graphs

Mycolyltransferase is important for biofilm formation and pathogenesis of Tsukamurella keratitis

Tsukamurella, a group of multi-drug resistant, Gram-positive, aerobic, and partially acid-fast bacteria, are emerging causes of bacterial conjunctivitis and keratitis. However, the pathogenesis of Tsukamurella keratitis is largely unknown. To address this, we used New Zealand White rabbits to develop the first eye infection model and conducted in vitro tests to study the pathogenesis mechanisms of Tsukamurella. There is increasing evidence that biofilms play a significant role in ocular infections, leading us to hypothesize that biofilm formation is crucial for effective Tsukamurella infection. In order to look for potential candidate genes which are important in biofilm formation and Tsukamurella keratitis. We performed genome sequencing of two ocular isolates, T. pulmonis-PW1004 and T. tyrosinosolvens-PW899, to identify potential virulence factors. Through in vitro and in vivo studies, we characterized their biological roles in mediating Tsukamurella keratitis. Our findings confirmed that Tsukamurella is an ocular pathogen by fulfilling Koch's postulates, and using genome sequence data, we identified tmytC, encoding a mycolyltransferase, as a crucial gene in biofilm formation and causing Tsukamurella keratitis in the rabbit model. This is the first report demonstrating the novel role of mycolyltransferase in causing ocular infections. Overall, our findings contribute to a better understanding of Tsukamurella pathogenesis and provide a potential target for treatment. Specific inhibitors targeting TmytC could serve as an effective treatment option for Tsukamurella infections.

Oral microbiome insights: Tracing acidic culprits in dental caries with functional metagenomics

OBJECTIVE

This study aimed to investigate the presence and abundance of acid-producing bacteria in dental caries samples using functional gene prediction techniques.

DESIGN

A total of 24 dental caries samples were collected for analysis. DNA isolation was performed followed by shotgun metagenomic sequencing. Functional gene prediction techniques were used to identify enzymes responsible for acid production from primary metabolites. Enzymes responsible for converting primary metabolites into acids were identified from the KEGG database. Subsequently, 840 contigs were examined, and their genus and species were characterized.

RESULTS

Analysis of the obtained data revealed 31 KEGG IDs corresponding to enzymes involved in the conversion of primary metabolites into acids. All 117 identified genera from the contig analysis were found to be part of the oral microbiome. In addition, A higher prevalence of acid-producing bacteria was noted in dental caries samples compared to earlier reports.

CONCLUSION

The study indicates the significant role of acid-producing bacteria in the initiation and progression of dental caries. The findings highlight the importance of microbial activity in the demineralization process of tooth enamel. Methods for preventing dental decay may be promising if specific measures are implemented to reduce the amount of acid produced by oral bacteria.

Unitig-centered pan-genome machine learning approach for predicting antibiotic resistance and discovering novel resistance genes in bacterial strains

In current genomic research, the widely used methods for predicting antimicrobial resistance (AMR) often rely on prior knowledge of known AMR genes or reference genomes. However, these methods have limitations, potentially resulting in imprecise predictions owing to incomplete coverage of AMR mechanisms and genetic variations. To overcome these limitations, we propose a pan-genome-based machine learning approach to advance our understanding of AMR gene repertoires and uncover possible feature sets for precise AMR classification. By building compacted de Brujin graphs (cDBGs) from thousands of genomes and collecting the presence/absence patterns of unique sequences (unitigs) for Pseudomonas aeruginosa, we determined that using machine learning models on unitig-centered pan-genomes showed significant promise for accurately predicting the antibiotic resistance or susceptibility of microbial strains. Applying a feature-selection-based machine learning algorithm led to satisfactory predictive performance for the training dataset (with an area under the receiver operating characteristic curve (AUC) of > 0.929) and an independent validation dataset (AUC, approximately 0.77). Furthermore, the selected unitigs revealed previously unidentified resistance genes, allowing for the expansion of the resistance gene repertoire to those that have not previously been described in the literature on antibiotic resistance. These results demonstrate that our proposed unitig-based pan-genome feature set was effective in constructing machine learning predictors that could accurately identify AMR pathogens. Gene sets extracted using this approach may offer valuable insights into expanding known AMR genes and forming new hypotheses to uncover the underlying mechanisms of bacterial AMR.

Context-dependent genomic locus effects on antibody production in recombinant Chinese hamster ovary cells generated through random integration

High-yield production of therapeutic protein using Chinese hamster ovary (CHO) cells requires stable cell line development (CLD). CLD typically uses random integration of transgenes; however, this results in clonal variation and subsequent laborious clone screening. Therefore, site-specific integration of a protein expression cassette into a desired chromosomal locus showing high transcriptional activity and stability, referred to as a hot spot, is emerging. Although positional effects are important for therapeutic protein expression, the sequence-specific mechanisms by which hotspots work are not well understood. In this study, we performed whole-genome sequencing (WGS) to locate randomly inserted vectors in the genome of recombinant CHO cells expressing high levels of monoclonal antibodies (mAbs) and experimentally validated these locations and vector compositions. The integration site was characterized by active histone marks and potential enhancer activities, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) mediated indel mutations in the region upstream of the integration site led to a significant reduction in specific antibody productivity by up to 30%. Notably, the integration site and its core region did not function equivalently outside the native genomic context, showing a minimal effect on the increase in exogenous protein expression in the host cell line. We also observed a superior production capacity of the mAb expressing cell line compared to that of the host cell line. Collectively, this study demonstrates that developing recombinant CHO cell lines to produce therapeutic proteins at high levels requires a balance of factors including transgene configuration, genomic locus landscape, and host cell properties.

High-quality genome of a modern soybean cultivar and resequencing of 547 accessions provide insights into the role of structural variation

Soybean provides protein, oil and multiple health-related compounds. Understanding the effects of structural variations (SVs) on economic traits in modern breeding is important for soybean improvement. Here we assembled the high-quality genome of modern cultivar Nongdadou2 (NDD2) and identified 25,814 SV-gene pairs compared to 29 reported genomes, with 13 NDD2-private SVs validated in 547 deep-resequencing (average = 18.05-fold) accessions, which advances our understanding of genomic variation biology. We found some insertions/deletions involved in seed protein and weight formation, an inversion related to adaptation to drought and a large intertranslocation implicated in a key divergence event in soybean. Of 749,714 SVs from 547 accessions, 6,013 were significantly associated with 22 yield-related and seed-quality-related traits determined in ten location × year environments. We uncovered 1,761 associated SVs that hit genes or regulatory regions, with 12 in GmMQT influencing oil and isoflavone contents. Our work provides resources and insights into SV roles in soybean improvement.

S. Barkham T. Multilocus sequence typing database for Streptococcus agalactiae contains a spurious allele of the transketolase gene

The tkt (transketolase) gene is one of the seven gene fragments used in the multilocus sequence typing (MLST) system for Streptococcus agalactiae. We discovered that the tkt_134 allele is derived from a homologous gene (which we designate tktX) that is not present in all S. agalactiae; all known strains that contain a match to the tkt_134 allele also contain a gene sequence that is much closer in sequence identity to the other non-tkt_134 alleles (i.e., the canonical tkt gene) in the database. Based on these data, the tkt_134 allele has been removed from the MLST database as of September 2021, and all sequence types containing tkt_134 have also been removed.IMPORTANCEMultilocus sequence typing (MLST) databases are a common good and remain important for research, medical, and epidemiological purposes. This remains true even in the context of widespread whole-genome sequencing. We discovered a contaminating allele of the tkt gene in the S. agalactiae MLST database that led to unstable, ambiguous, or erroneous MLST assignment. The allele has since been removed from the public database based on the results presented in this manuscript.

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

Applications of bacteriophage in combination with nisin for controlling multidrug-resistant Bacillus cereus in broth and various food matrices

This study focused on the isolation and characterization of bacteriophages with specific activity against toxin-producing and multidrug-resistant strains of Bacillus cereus sensu stricto (B. cereus s. s.). Ten different samples yielded six bacteriophages by utilizing the double-layer agar technique. The most promising phage, vB_BceS-M2, was selected based on its broad host range and robust lytic activity against various B. cereus s. s. strains. The phage vB_BceS-M2 had a circular double-stranded DNA genome of 56,482 bp. This phage exhibited stability over a wide range of temperatures and pH values, which is crucial for its potential application in food matrices. The combined effect of phage vB_BceS-M2 and nisin, a widely used antimicrobial peptide, was investigated to enhance antimicrobial efficacy against B. cereus in food. The results suggested that nisin showed synergy and combined effect with the phage, potentially overcoming the growth of phage-resistant bacteria in the broth. Furthermore, practical applications were conducted in various liquid and solid food matrices, including whole and skimmed milk, boiled rice, cheese, and frozen meatballs, both at 4 and 25 °C. Phage vB_BceS-M2, either alone or in combination with nisin, reduced the growth rate of B. cereus in foods other than whole milk. The combination of bacteriophage and nisin showed promise for the development of effective antimicrobial interventions to counteract toxigenic and antibiotic-resistant B. cereus in food.

Functional Identification and Genetic Analysis of O-Antigen Gene Clusters of Food-Borne Pathogen Yersinia enterocolitica O:10 and Other Uncommon Serotypes, Further Revealing Their Virulence Profiles

Yersinia enterocolitica is a globally distributed food-borne gastrointestinal pathogen. The O-antigen variation-determined serotype is an important characteristic of Y. enterocolitica, allowing intraspecies classification for diagnosis and epidemiology purposes. Among the 11 serotypes associated with human yersiniosis, O:3, O:5,27, O:8, and O:9 are the most prevalent, and their O-antigen gene clusters have been well defined. In addition to the O-antigen, several virulence factors are involved in infection and pathogenesis of Y. enterocolitica strains, and these are closely related to their biotypes, reflecting pathogenic properties. In this study, we identified the O-AGC of a Y. enterocolitica strain WL-21 of serotype O:10, and confirmed its functionality in O-antigen synthesis. Furthermore, we analyzed in silico the putative O-AGCs of uncommon serotypes, and found that the O-AGCs of Y. enterocolitica were divided into two genetic patterns: (1) O-AGC within the hemH-gsk locus, possibly synthesizing the O-antigen via the Wzx/Wzy dependent pathway, and (2) O-AGC within the dcuC-galU-galF locus, very likely assembling the O-antigen via the ABC transporter dependent pathway. By screening the virulence genes against genomes from GenBank, we discovered that strains representing different serotypes were grouped according to different virulence gene profiles, indicating strong links between serotypes and virulence markers and implying an interaction between them and the synergistic effect in pathogenicity. Our study provides a framework for further research on the origin and evolution of O-AGCs from Y. enterocolitica, as well as on differences in virulent mechanisms among distinct serotypes.

Draft genome sequences of Butyrivibrio hungatei DSM 14810 (JK 615T) and Butyrivibrio fibrisolvens DSM 3071 (D1T)

Here, we report the draft genome sequences of two Butyrivibrio-type strains isolated from rumen fluid. The genome sequence of Butyrivibrio hungatei DSM 14810 was 3.3 Mb with 3,093 predicted genes, while the Butyrivibrio fibrisolvens DSM 3071 genome sequence was 4.8 Mb with 4,132 predicted genes.

Isolation and Characterization of a Novel Phage against Vibrio alginolyticus Belonging to a New Genus

Vibrio alginolyticus causes substantial economic losses in the aquaculture industry. With the rise of multidrug-resistant Vibrio strains, phages present a promising solution. Here, a novel lytic Vibrio phage, vB_ValC_RH2G (RH2G), that efficiently infects the pathogenic strain V. alginolyticus ATCC 17749T, was isolated from mixed wastewater from an aquatic market in Xiamen, China. Transmission electron microscopy revealed that RH2G has the morphology of Siphoviruses, featuring an icosahedral head (73 ± 2 nm diameter) and long noncontractile tail (142 ± 4 nm). A one-step growth experiment showed that RH2G had a short latent period (10 min) and a burst size of 48 phage particles per infected cell. Additionally, RH2G was highly species-specific and was relatively stable at 4-55 °C and pH 4-10. A genomic analysis showed that RH2G has a 116,749 bp double-stranded DNA genome with 43.76% GC content. The intergenomic similarity between the genome sequence of RH2G and other phages recorded in the GenBank database was below 38.8%, suggesting that RH2G represents a new genus. RH2G did not exhibit any virulence or resistance genes. Its rapid lysis capacity, lytic activity, environmental resilience, and genetic safety suggested that RH2G may be a safe candidate for phage therapy in combatting vibriosis in aquaculture settings.

Development and implementation of a core genome multilocus sequence typing scheme for Yersinia enterocolitica: a tool for surveillance and outbreak detection

Yersinia enterocolitica (Y. enterocolitica) is the most frequent etiological agent of yersiniosis and has been responsible for several national outbreaks in Norway and elsewhere. A standardized high-resolution method, such as core genome Multilocus Sequence Typing (cgMLST), is needed for pathogen traceability at the national and international levels. In this study, we developed and implemented a cgMLST scheme for Y. enterocolitica. We designed a cgMLST scheme in SeqSphere + using high-quality genomes from different Y. enterocolitica biotype sublineages. The scheme was validated if more than 95% of targets were found across all tested Y. enterocolitica: 563 Norwegian genomes collected between 2012 and 2022 and 327 genomes from public data sets. We applied the scheme to known outbreaks to establish a threshold for identifying major complex types (CTs) based on the number of allelic differences. The final cgMLST scheme included 2,582 genes with a median of 97.9% (interquartile range 97.6%-98.8%) targets found across all tested genomes. Analysis of outbreaks identified all outbreak strains using single linkage clustering at four allelic differences. This threshold identified 311 unique CTs in Norway, of which CT18, CT12, and CT5 were identified as the most frequently associated with outbreaks. The cgMLST scheme showed a very good performance in typing Y. enterocolitica using diverse data sources and was able to identify outbreak clusters. We recommend the implementation of this scheme nationally and internationally to facilitate Y. enterocolitica surveillance and improve outbreak response in national and cross-border outbreaks.

Maptcha: an efficient parallel workflow for hybrid genome scaffolding

BACKGROUND

Genome assembly, which involves reconstructing a target genome, relies on scaffolding methods to organize and link partially assembled fragments. The rapid evolution of long read sequencing technologies toward more accurate long reads, coupled with the continued use of short read technologies, has created a unique need for hybrid assembly workflows. The construction of accurate genomic scaffolds in hybrid workflows is complicated due to scale, sequencing technology diversity (e.g., short vs. long reads, contigs or partial assemblies), and repetitive regions within a target genome.

RESULTS

In this paper, we present a new parallel workflow for hybrid genome scaffolding that would allow combining pre-constructed partial assemblies with newly sequenced long reads toward an improved assembly. More specifically, the workflow, called Maptcha, is aimed at generating long scaffolds of a target genome, from two sets of input sequences-an already constructed partial assembly of contigs, and a set of newly sequenced long reads. Our scaffolding approach internally uses an alignment-free mapping step to build a ⟨ contig,contig ⟩ graph using long reads as linking information. Subsequently, this graph is used to generate scaffolds. We present and evaluate a graph-theoretic "wiring" heuristic to perform this scaffolding step. To enable efficient workload management in a parallel setting, we use a batching technique that partitions the scaffolding tasks so that the more expensive alignment-based assembly step at the end can be efficiently parallelized. This step also allows the use of any standalone assembler for generating the final scaffolds.

CONCLUSIONS

Our experiments with Maptcha on a variety of input genomes, and comparison against two state-of-the-art hybrid scaffolders demonstrate that Maptcha is able to generate longer and more accurate scaffolds substantially faster. In almost all cases, the scaffolds produced by Maptcha are at least an order of magnitude longer (in some cases two orders) than the scaffolds produced by state-of-the-art tools. Maptcha runs significantly faster too, reducing time-to-solution from hours to minutes for most input cases. We also performed a coverage experiment by varying the sequencing coverage depth for long reads, which demonstrated the potential of Maptcha to generate significantly longer scaffolds in low coverage settings ( 1 × - 10 × ).

Genomic Analysis of Cronobacter condimenti s37: Identification of Resistance and Virulence Genes and Comparison with Other Cronobacter and Closely Related Species

Cronobacter condimenti are environmental commensals that have not been associated with any clinical infections. To date, they are the least understood and described Cronobacter species within the genus. The objective of this study was to use a draft genome sequence (DGS) of the Cronobacter condimenti strain s37 to screen for genes encoding for antibiotic resistance, virulence, response to environmental stress, and biofilm formation. The strain was isolated in Poland from commercial small radish sprouts. This is the second genome of this species available in the GenBank database. The comparative genome analysis (cgMLST) of C. condimenti s37 with other Cronobacter spp. including the pathogenic species C. sakazakii and the plant-associated closely related genera Franconibacter and Siccibacter was also performed. The assembled and annotated genome of the C. condimenti s37 genome was 4,590,991 bp in length, with a total gene number of 4384, and a GC content of 55.7%. The s 37 genome encoded for genes associated with resistance to stressful environmental conditions (metal resistance genes: zinc, copper, osmotic regulation, and desiccation stress), 17 antimicrobial resistance genes encoding resistance to various classes of antibiotics and 50 genes encoding for the virulence factors. The latter were mainly genes associated with adhesion, chemotaxis, hemolysis, and biofilm formation. Cg-MLST analysis (3991 genes) revealed a greater similarity of C. condimenti s37 to S. turicensis, F. pulveris, and C. dublinensis than to other species of the genus Cronobacter. Studies on the diversity, pathogenicity, and virulence of Cronobacter species isolated from different sources are still insufficient and should certainly be continued. Especially the analysis of rare strains such as s37 is very important because it provides new information on the evolution of these bacteria. Comparative cgMLST analysis of s37 with other Cronobacter species, as well as closely related genera Franconibacter and Siccibacter, complements the knowledge on their adaptability to specific environments such as desiccation.

A robust genome assembly with transcriptomic data from the striped bark scorpion, Centruroides vittatus

Scorpions, a seemingly primitive, stinging arthropod taxa, are known to exhibit marked diversity in their venom components. These venoms are known for their human pathology, but they are also important as models for therapeutic and drug development applications. In this study, we report a high-quality genome assembly and annotation of the striped bark scorpion, Centruroides vittatus, created with several shotgun libraries. The final assembly is 760 Mb in size, with a BUSCO score of 97.8%, a 30.85% GC, and an N50 of 2.35 Mb. We estimated 36,189 proteins with 37.32% assigned to Gene Ontology (GO) terms in our GO annotation analysis. We mapped venom toxin genes to 18 contigs and 2 scaffolds. We were also able to identify expression differences between venom gland (telson) and body tissue (carapace) with 19 sodium toxin and 14 potassium toxin genes to 18 contigs and 2 scaffolds. This assembly, along with our transcriptomic data, provides further data to investigate scorpion venom genomics.

Characterization of selected phages for biocontrol of food-spoilage pseudomonads

Pseudomonas spp., such as P. fluorescens group, P. fragi, and P. putida, are the major psychrophilic spoilage bacteria in the food industry. Bacteriophages (phages) are a promising tool for controlling food-spoilage and food-poisoning bacteria; however, there are few reports on phages effective on food-spoilage bacteria such as Pseudomonas spp. In this study, 12 Pseudomonas phages were isolated from chicken and soil samples. Based on the host range and lytic activity at 30 °C and 4 °C and various combinations of phages, phages vB_PflP-PCS4 and vB_PflP-PCW2 were selected to prepare phage cocktails to control Pseudomonas spp. The phage cocktail consisting of vB_PflP-PCS4 and vB_PflP-PCW2 showed the strongest lytic activity and retarded regrowth of P. fluorescens and P. putida at 30 °C, 8 °C, and 4 °C at a multiplicity of infection of 100. Nucleotide sequence analysis of the genomic DNA indicated that vB_PflP-PCS4 and vB_PflP-PCW2 phages were lytic phages of the Podoviridae family and lacked tRNA, toxin, or virulence genes. A novel endolysin gene was found in the genomic DNA of phage vB_PflP-PCS4. The results of this study suggest that the phage cocktail consisting of vB_PflP-PCS4 and vB_PflP-PCW2 is a promising tool for the biocontrol of psychrophilic food-spoilage pseudomonads during cold storage and distribution.

Virulence phenotypes differ between toxigenic Vibrio parahaemolyticus isolated from western coasts of Europe

Vibrio parahaemolyticus is the leading bacterial cause of gastroenteritis associated with seafood consumption worldwide. Not all members of the species are thought to be pathogenic, thus identification of virulent organisms is essential to protect public health and the seafood industry. Correlations of human disease and known genetic markers (e.g. thermostable direct hemolysin (TDH), TDH-related hemolysin (TRH)) appear complex. Some isolates recovered from patients lack these factors, while their presence has become increasingly noted in isolates recovered from the environment. Here, we used whole-genome sequencing in combination with mammalian and insect models of infection to assess the pathogenic potential of V. parahaemolyticus isolated from European Atlantic shellfish production areas. We found environmental V. parahaemolyticus isolates harboured multiple virulence-associated genes, including TDH and/or TRH. However, carriage of these factors did not necessarily reflect virulence in the mammalian intestine, as an isolate containing TDH and the genes coding for a type 3 secretion system (T3SS) 2α virulence determinant, appeared avirulent. Moreover, environmental V. parahaemolyticus lacking TDH or TRH could be assigned to groups causing low and high levels of mortality in insect larvae, with experiments using defined bacterial mutants showing that a functional T3SS1 contributed to larval death. When taken together, our findings highlight the genetic diversity of V. parahaemolyticus isolates found in the environment, their potential to cause disease and the need for a more systematic evaluation of virulence in diverse V. parahaemolyticus to allow better genetic markers.

Closing the genome of Teredinibacter turnerae T7902 by long-read nanopore sequencing

We present the complete closed circular genome sequence derived from Oxford Nanopore sequencing of the shipworm endosymbiont Teredinibacter turnerae T7902 (DSM 15152, ATCC 39867), originally isolated from the shipworm Lyrodus pedicellatus (1). This sequence will aid in the comparative genomics of shipworm endosymbionts and the understanding of host-symbiont evolution.

Local Genomic Surveillance of Invasive Streptococcus pyogenes in Eastern North Carolina (ENC) in 2022-2023

The recent increase in Group A Streptococcus (GAS) incidences in several countries across Europe and some areas of the Unites States (U.S.) has raised concerns. To understand GAS diversity and prevalence, we conducted a local genomic surveillance in Eastern North Carolina (ENC) in 2022-2023 with 95 isolates and compared its results to those of the existing national genomic surveillance in the U.S. in 2015-2021 with 13,064 isolates. We observed their epidemiological changes before and during the COVID-19 pandemic and detected a unique sub-lineage in ENC among the most common invasive GAS strain, ST28/emm1. We further discovered a multiple-copy insertion sequence, ISLgar5, in ST399/emm77 and its single-copy variants in some other GAS strains. We discovered ISLgar5 was linked to a Tn5801-like tetM-carrying integrative and conjugative element, and its copy number was associated with an ermT-carrying pRW35-like plasmid. The dynamic insertions of ISLgar5 may play a vital role in genome fitness and adaptation, driving GAS evolution relevant to antimicrobial resistance and potentially GAS virulence.

Exploring gene content with pangene graphs

MOTIVATION

The gene content regulates the biology of an organism. It varies between species and between individuals of the same species. Although tools have been developed to identify gene content changes in bacterial genomes, none is applicable to collections of large eukaryotic genomes such as the human pangenome.

RESULTS

We developed pangene, a computational tool to identify gene orientation, gene order and gene copy-number changes in a collection of genomes. Pangene aligns a set of input protein sequences to the genomes, resolves redundancies between protein sequences and constructs a gene graph with each genome represented as a walk in the graph. It additionally finds subgraphs, which we call bibubbles, that capture gene content changes. Applied to the human pangenome, pangene identifies known gene-level variations and reveals complex haplotypes that are not well studied before. Pangene also works with high-quality bacterial pangenome and reports similar numbers of core and accessory genes in comparison to existing tools.

AVAILABILITY AND IMPLEMENTATION

Source code at https://github.com/lh3/pangene; pre-built pangene graphs can be downloaded from https://zenodo.org/records/8118576 and visualized at https://pangene.bioinweb.org.

Genomic Analysis of Aspergillus Section Terrei Reveals a High Potential in Secondary Metabolite Production and Plant Biomass Degradation

Aspergillus terreus has attracted interest due to its application in industrial biotechnology, particularly for the production of itaconic acid and bioactive secondary metabolites. As related species also seem to possess a prosperous secondary metabolism, they are of high interest for genome mining and exploitation. Here, we present draft genome sequences for six species from Aspergillus section Terrei and one species from Aspergillus section Nidulantes. Whole-genome phylogeny confirmed that section Terrei is monophyletic. Genome analyses identified between 70 and 108 key secondary metabolism genes in each of the genomes of section Terrei, the highest rate found in the genus Aspergillus so far. The respective enzymes fall into 167 distinct families with most of them corresponding to potentially unique compounds or compound families. Moreover, 53% of the families were only found in a single species, which supports the suitability of species from section Terrei for further genome mining. Intriguingly, this analysis, combined with heterologous gene expression and metabolite identification, suggested that species from section Terrei use a strategy for UV protection different to other species from the genus Aspergillus. Section Terrei contains a complete plant polysaccharide degrading potential and an even higher cellulolytic potential than other Aspergilli, possibly facilitating additional applications for these species in biotechnology.

A comprehensive review of deep learning-based variant calling methods

Genome sequencing data have become increasingly important in the field of personalized medicine and diagnosis. However, accurately detecting genomic variations remains a challenging task. Traditional variation detection methods rely on manual inspection or predefined rules, which can be time-consuming and prone to errors. Consequently, deep learning-based approaches for variation detection have gained attention due to their ability to automatically learn genomic features that distinguish between variants. In our review, we discuss the recent advancements in deep learning-based algorithms for detecting small variations and structural variations in genomic data, as well as their advantages and limitations.

Genetic diversity, determinants, and dissemination of Burkholderia pseudomallei lineages implicated in melioidosis in Northeast Thailand

Melioidosis is an often-fatal neglected tropical disease caused by an environmental bacterium Burkholderia pseudomallei. However, our understanding of the disease-causing bacterial lineages, their dissemination, and adaptive mechanisms remains limited. To address this, we conduct a comprehensive genomic analysis of 1,391 B. pseudomallei isolates collected from nine hospitals in northeast Thailand between 2015 and 2018, and contemporaneous isolates from neighbouring countries, representing the most densely sampled collection to date. Our study identifies three dominant lineages, each with unique gene sets potentially enhancing bacterial fitness in the environment. We find that recombination drives lineage-specific gene flow. Transcriptome analyses of representative clinical isolates from each dominant lineage reveal increased expression of lineage-specific genes under environmental conditions in two out of three lineages. This underscores the potential importance of environmental persistence for these dominant lineages. The study also highlights the influence of environmental factors such as terrain slope, altitude, and river direction on the geographical dispersal of B. pseudomallei. Collectively, our findings suggest that environmental persistence may play a role in facilitating the spread of B. pseudomallei, and as a prerequisite for exposure and infection, thereby providing useful insights for informing melioidosis prevention and control strategies.

Evolution and host-specific adaptation of Pseudomonas aeruginosa

The major human bacterial pathogen Pseudomonas aeruginosa causes multidrug-resistant infections in people with underlying immunodeficiencies or structural lung diseases such as cystic fibrosis (CF). We show that a few environmental isolates, driven by horizontal gene acquisition, have become dominant epidemic clones that have sequentially emerged and spread through global transmission networks over the past 200 years. These clones demonstrate varying intrinsic propensities for infecting CF or non-CF individuals (linked to specific transcriptional changes enabling survival within macrophages); have undergone multiple rounds of convergent, host-specific adaptation; and have eventually lost their ability to transmit between different patient groups. Our findings thus explain the pathogenic evolution of P. aeruginosa and highlight the importance of global surveillance and cross-infection prevention in averting the emergence of future epidemic clones.

Acute Promyelocytic Leukemia With Torque Teno Mini Virus::RARA Fusion: An Approach to Screening and Diagnosis

Acute promyelocytic leukemia (APL) with variant RARA translocation is linked to over 15 partner genes. Recent publications encompassing 6 cases have expanded the spectrum of RARA partners to torque teno mini virus (TTMV). This entity is likely underrecognized due to the lack of clinician and pathologist familiarity, inability to detect the fusion using routine testing modalities, and informatic challenges in its recognition within next-generation sequencing (NGS) data. We describe a clinicopathologic approach and provide the necessary tools to screen and diagnose APL with TTMV::RARA using existing clinical DNA- or RNA-based NGS assays, which led to the identification of 4 cases, all without other known cytogenetic/molecular drivers. One was identified prospectively and 3 retrospectively, including 2 from custom automated screening of multiple data sets (50,257 cases of hematopoietic malignancy, including 4809 acute myeloid leukemia/myeloid sarcoma/APL cases). Two cases presented as myeloid sarcoma, including 1 with multiple relapses after acute myeloid leukemia-type chemotherapy and hematopoietic stem cell transplant. Two cases presented as leukemia, had a poor response to induction chemotherapy, but achieved remission upon reinduction (including all-trans retinoic acid in 1 case) and subsequent hematopoietic stem cell transplant. Neoplastic cells demonstrated features of APL including frequent azurophilic granules and dim/absent CD34 and HLA-DR expression. RARA rearrangement was not detected by karyotype or fluorescent in situ hybridization. Custom analysis of NGS fusion panel data identified TTMV::RARA rearrangements and, in the prospectively identified case, facilitated monitoring in sequential bone marrow samples. APL with TTMV::RARA is a rare leukemia with a high rate of treatment failure in described cases. The diagnosis should be considered in leukemias with features of APL that lack detectable RARA fusions and other drivers, and may be confirmed by appropriate NGS tests with custom informatics. Incorporation of all-trans retinoic acid may have a role in treatment but requires accurate recognition of the fusion for appropriate classification as APL.

Unravelling the genomic and environmental diversity of the ubiquitous Solirubrobacter

Solirubrobacter, though widespread in soils and rhizospheres, has been relatively unexplored despite its ubiquity. Previously acknowledged as a common soil bacterium, our research explores its phylogenomics, pangenomics, environmental diversity, and interactions within bacterial communities. By analysing seven genomic sequences, we have identified a pangenome consisting of 19,645 protein families, of which 2644 are shared across all studied genomes, forming the core genome. Interestingly, despite the non-motility of reported isolates, we discovered genes for flagellin and a partial flagellum assembly pathway. Examining the 16S ribosomal RNA genes of Solirubrobacter revealed substantial diversity, with 3166 operational taxonomic units identified in Mexican soils. Co-occurrence network analysis further demonstrated its significant integration within bacterial communities. Through phylogenomic scrutiny, we conclusively excluded the NCBI's GCA_009993245.1 genome from being classified as a Solirubrobacter. Our research into the metagenomic diversity of Solirubrobacter across various environments confirmed its presence in rhizospheres and certain soils, underscoring its adaptability. The geographical ubiquity of Solirubrobacter in rhizospheres raises intriguing questions regarding its potential interactions with plant hosts and the biotic and abiotic factors influencing its presence in soil. Given its ecological significance and genetic diversity, Solirubrobacter warrants further investigation as a potentially crucial yet underappreciated keystone species.

Geographical migration and fitness dynamics of Streptococcus pneumoniae

Streptococcus pneumoniae is a leading cause of pneumonia and meningitis worldwide. Many different serotypes co-circulate endemically in any one location1,2. The extent and mechanisms of spread and vaccine-driven changes in fitness and antimicrobial resistance remain largely unquantified. Here using geolocated genome sequences from South Africa (n = 6,910, collected from 2000 to 2014), we developed models to reconstruct spread, pairing detailed human mobility data and genomic data. Separately, we estimated the population-level changes in fitness of strains that are included (vaccine type (VT)) and not included (non-vaccine type (NVT)) in pneumococcal conjugate vaccines, first implemented in South Africa in 2009. Differences in strain fitness between those that are and are not resistant to penicillin were also evaluated. We found that pneumococci only become homogenously mixed across South Africa after 50 years of transmission, with the slow spread driven by the focal nature of human mobility. Furthermore, in the years following vaccine implementation, the relative fitness of NVT compared with VT strains increased (relative risk of 1.68; 95% confidence interval of 1.59-1.77), with an increasing proportion of these NVT strains becoming resistant to penicillin. Our findings point to highly entrenched, slow transmission and indicate that initial vaccine-linked decreases in antimicrobial resistance may be transient.

Transition of survival strategies under global climate shifts in the grape family

Faced with environmental changes, plants may either move to track their ancestral niches or evolve to adapt to new niches. Vitaceae, the grape family, has evolved diverse adaptive traits facilitating a global expansion in wide-ranging habitats, making it ideal for investigating transition between move and evolve strategies and exploring the underlying mechanisms. Here we inferred the patterns of biogeographic diversification and trait evolution in Vitaceae based on a robust phylogeny with dense sampling including 495 species (~52% of Vitaceae species). Vitaceae probably originated from Asia-the diversity centre of extant genera and the major source of dispersals. Boundaries of the Eocene, Oligocene and Miocene were identified as turning points in shifting strategies. A significant decrease in move strategy was identified during the Oligocene, followed by increases in move and evolve. After the Miocene, evolve began to dominate, during which increased niche opportunities and key trait innovations played important roles.

UnCoVar: a reproducible and scalable workflow for transparent and robust virus variant calling and lineage assignment using SARS-CoV-2 as an example

BACKGROUND

At a global scale, the SARS-CoV-2 virus did not remain in its initial genotype for a long period of time, with the first global reports of variants of concern (VOCs) in late 2020. Subsequently, genome sequencing has become an indispensable tool for characterizing the ongoing pandemic, particularly for typing SARS-CoV-2 samples obtained from patients or environmental surveillance. For such SARS-CoV-2 typing, various in vitro and in silico workflows exist, yet to date, no systematic cross-platform validation has been reported.

RESULTS

In this work, we present the first comprehensive cross-platform evaluation and validation of in silico SARS-CoV-2 typing workflows. The evaluation relies on a dataset of 54 patient-derived samples sequenced with several different in vitro approaches on all relevant state-of-the-art sequencing platforms. Moreover, we present UnCoVar, a robust, production-grade reproducible SARS-CoV-2 typing workflow that outperforms all other tested approaches in terms of precision and recall.

CONCLUSIONS

In many ways, the SARS-CoV-2 pandemic has accelerated the development of techniques and analytical approaches. We believe that this can serve as a blueprint for dealing with future pandemics. Accordingly, UnCoVar is easily generalizable towards other viral pathogens and future pandemics. The fully automated workflow assembles virus genomes from patient samples, identifies existing lineages, and provides high-resolution insights into individual mutations. UnCoVar includes extensive quality control and automatically generates interactive visual reports. UnCoVar is implemented as a Snakemake workflow. The open-source code is available under a BSD 2-clause license at github.com/IKIM-Essen/uncovar.

Acinetobacter pittii: the emergence of a hospital-acquired pathogen analyzed from the genomic perspective

Acinetobacter pittii has increasingly been associated with several types of hospital-acquired severe infections. Genes implicated in carbapenem resistance, tigecycline resistance, or genes encoding extended spectrum cephalosporinases, such as blaADC, are commonly found in isolates implicated in these infections. A. pittii strains that are pandrug resistant have occasionally been identified. Food for human consumption, animals and plants are environmental sources of this pathogen. An alarming situation is that A. pitti has been identified as responsible for outbreaks in different regions worldwide. In this study, 384 genomes of A. pittii were analyzed, comprising sequences from clinical and non-clinical origins from 32 countries. The objective was to investigate if clinical strains possess genetic traits facilitating hospital adaptation. Results indicate significant genomic variability in terms of size and gene content among A. pittii isolates. The core genome represents a small portion (25-36%) of each isolate's genome, while genes associated with antibiotic resistance and virulence predominantly belong to the accessory genome. Notably, antibiotic resistance genes are encoded by a diverse array of plasmids. As the core genome between environmental and hospital isolates is the same, we can assume that hospital isolates acquired ARGs due to a high selective pressure in these settings. The strain's phylogeographic distribution indicates that there is no geographical bias in the isolate distribution; isolates from different geographic regions are dispersed throughout a core genome phylogenetic tree. A single clade may include isolates from extremely distant geographical areas. Furthermore, strains isolated from the environment or animal, or plant sources frequently share the same clade as hospital isolates. Our analysis showed that the clinical isolates do not already possess specific genes, other than antibiotic-resistant genes, to thrive in the hospital setting.

Genomic and biological insights of bacteriophages JNUWH1 and JNUWD in the arms race against bacterial resistance

The coevolution of bacteria and bacteriophages has created a great diversity of mechanisms by which bacteria fight phage infection, and an equivalent diversity of mechanisms by which phages subvert bacterial immunity. Effective and continuous evolution by phages is necessary to deal with coevolving bacteria. In this study, to better understand the connection between phage genes and host range, we examine the isolation and genomic characterization of two bacteriophages, JNUWH1 and JNUWD, capable of infecting Escherichia coli. Sourced from factory fermentation pollutants, these phages were classified within the Siphoviridae family through TEM and comparative genomic analysis. Notably, the phages exhibited a viral burst size of 500 and 1,000 PFU/cell, with latent periods of 15 and 20 min, respectively. They displayed stability over a pH range of 5 to 10, with optimal activity at 37°C. The complete genomes of JNUWH1 and JNUWD were 44,785 bp and 43,818 bp, respectively. Phylogenetic analysis revealed their close genetic relationship to each other. Antibacterial assays demonstrated the phages' ability to inhibit E. coli growth for up to 24 h. Finally, through laboratory-driven adaptive evolution, we successfully identified strains for both JNUWH1 and JNUWD with mutations in receptors specifically targeting lipopolysaccharides (LPS) and the lptD gene. Overall, these phages hold promise as additives in fermentation products to counter E. coli, offering potential solutions in the context of evolving bacterial resistance.

Molecular and Phylogenomic Analysis of a Vancomycin Intermediate Resistance USA300LV Strain in Chile

Antimicrobial resistance is a major global health problem, and, among Gram-positive bacteria, methicillin-resistant Staphylococcus aureus (MRSA) represents a serious threat. MRSA causes a wide range of infections, including bacteremia, which, due to the limited use of β-lactams, is difficult to treat. This study aimed to analyze 51 MRSA isolates collected in 2018 from samples of patients with bacteremia from two hospitals of the Metropolitan Health Service of Santiago, Chile, both in their resistance profile and in the identification of virulence factors. In addition, genomic characterization was carried out by the WGS of an isolate that was shown to be the one of greatest concern (N°. 42) due to its intermediate resistance to vancomycin, multiple virulence factors and being classified as ST8 PVL-positive. In our study, most of the isolates turned out to be multidrug-resistant, but there are still therapeutic options, such as tetracycline, rifampicin, chloramphenicol and vancomycin, which are currently used for MRSA infections; however, 18% were PVL positive, which suggests greater virulence of these isolates. It was determined that isolate N°42 is grouped within the USA300-LV strains (ST8, PVL+, COMER+); however, it has been suggested that, in Chile, a complete displacement of the PVL-negative ST5 clone has not occurred.

CT135 mediates the resistance of Chlamydia trachomatis to primate interferon gamma stimulated immune defenses

Evading host innate immune defenses is a critical feature of Chlamydia trachomatis infections, and the mechanisms used by C. trachomatis to subvert these pathways are incompletely understood. We screened a library of chimeric C. trachomatis mutants for genetic factors important for interference with cell-autonomous immune defenses. Mutant strains with predicted truncations of the inclusion membrane protein CT135 were susceptible to interferon gamma-activated immunity in human cells. CT135 functions to prevent host-driven recruitment of ubiquitin and p62/SQSTM to the inclusion membrane. In a nonhuman primate model of C. trachomatis infection, a CT135-deficient strain was rapidly cleared, highlighting the importance of this virulence factor for C. trachomatis pathogenesis. Analysis of CT135 phenotypes in primary macaque cells revealed that cell-autonomous immune defenses against C. trachomatis are conserved between humans and nonhuman primates and connects mechanistic findings with in vivo infection outcomes.

New mitogenomes of Runcinidae and Facelinidae: two understudied heterobranch families (Mollusca: Gastropoda)

Here, we present the mitochondrial sequences of two sea slugs (Heterobranchia): Runcina aurata and Facelina auriculata, the latter being the type species of the family. The mitochondrial genomes are 14,282 and 14,171bp in length, respectively, with a complete set of 13 PCGs, 2 rRNAs, and 22 tRNAs. None of the mitogenomes show gene reorganization, keeping the standard mitogenomic structure of Heterobranchia. Nucleotide composition differs significantly between them, with R. aurata showing the most AT-rich mitogenome (25.7% GC content) reported to date in Heterobranchia, and F. auriculata showing a rich GC content (35%) compared with other heterobranch mitochondrial genomes.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

Isolation and identification of specific Enterococcus faecalis phage C-3 and G21-7 against Avian pathogenic Escherichia coli and its application to one-day-old geese

Colibacillosis caused by Avian pathogenic Escherichia coli (APEC), including peritonitis, respiratory tract inflammation and ovaritis, is recognized as one of the most common and economically destructive bacterial diseases in poultry worldwide. In this study, the characteristics and inhibitory potential of phages were investigated by double-layer plate method, transmission electron microscopy, whole genome sequencing, bioinformatics analysis and animal experiments. The results showed that phages C-3 and G21-7 isolated from sewage around goose farms infected multiple O serogroups (O1, O2, O18, O78, O157, O26, O145, O178, O103 and O104) Escherichia coli (E.coli) with a multiplicity of infection (MOI) of 10 and 1, respectively. According to the one-step growth curve, the incubation time of both bacteriophage C-3 and G21-7 was 10 min. Sensitivity tests confirmed that C-3 and G21-6 are stable at 4 to 50 °C and pH in the range of 4 to 11. Based on morphological and phylogenetic analysis, phages C-3 and G21-7 belong to Enterococcus faecalis (E. faecalis) phage species of the genus Saphexavirus of Herelleviridae family. According to genomic analysis, phage C-3 and G21-7 were 58,097 bp and 57,339 bp in size, respectively, with G+C content of 39.91% and 39.99%, encoding proteins of 97 CDS (105 to 3,993 bp) and 96 CDS (105 to 3,993 bp), and both contained 2 tRNAs. Both phages contained two tail proteins and holin-endolysin system coding genes, and neither carried resistance genes nor virulence factors. Phage mixture has a good safety profile and has shown good survival probability and feed efficiency in both treatment and prophylaxis experiments with one-day-old goslings. These results suggest that phage C-3 and G21-7 can be used as potential antimicrobials for the prevention and treatment of APEC.

Combining machine learning with high-content imaging to infer ciprofloxacin susceptibility in isolates of Salmonella Typhimurium

Antimicrobial resistance (AMR) is a growing public health crisis that requires innovative solutions. Current susceptibility testing approaches limit our ability to rapidly distinguish between antimicrobial-susceptible and -resistant organisms. Salmonella Typhimurium (S. Typhimurium) is an enteric pathogen responsible for severe gastrointestinal illness and invasive disease. Despite widespread resistance, ciprofloxacin remains a common treatment for Salmonella infections, particularly in lower-resource settings, where the drug is given empirically. Here, we exploit high-content imaging to generate deep phenotyping of S. Typhimurium isolates longitudinally exposed to increasing concentrations of ciprofloxacin. We apply machine learning algorithms to the imaging data and demonstrate that individual isolates display distinct growth and morphological characteristics that cluster by time point and susceptibility to ciprofloxacin, which occur independently of ciprofloxacin exposure. Using a further set of S. Typhimurium clinical isolates, we find that machine learning classifiers can accurately predict ciprofloxacin susceptibility without exposure to it or any prior knowledge of resistance phenotype. These results demonstrate the principle of using high-content imaging with machine learning algorithms to predict drug susceptibility of clinical bacterial isolates. This technique may be an important tool in understanding the morphological impact of antimicrobials on the bacterial cell to identify drugs with new modes of action.

Beneficial effects of Bacillus mojavensis strain MTC-8 on plant growth, immunity and disease resistance against Magnaporthe oryzae

Rice blast, a prevalent and highly destructive rice disease that significantly impacts rice yield, is caused by the rice blast fungus. In the present study, a strain named MTC-8, identified as Bacillus mojavensis, was demonstrated has strong antagonistic activity against the rice blast fungus, Rhizoctonia solani, Ustilaginoidea virens, and Bipolaria maydis. The potential biocontrol agents were identified using ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) analysis and chromatography. Further investigations elucidated the inhibitory mechanism of the isolated compound and demonstrated its ability to suppress spore germination, alter hyphal morphology, disrupt cell membrane integrity, and induce defense-related gene expression in rice. MTC-8 promoted plant growth and may lead to the development of a biocontrol agent that meets agricultural standards. Overall, the Bacillus mojavensis MTC-8 strain exerted beneficial effects on plant growth, immunity and disease resistance against rice blast fungus. In this study, we isolated and purified a bioactive substance from fermentation broth, and the results provide a foundation for the development and application of biopesticides. Elucidation of the inhibitory mechanism against rice blast fungus provides theoretical support for the identification of molecular targets. The successful development of a biocontrol agent lays the groundwork for its practical application in agriculture.

Genomic profiling and characterization of ocular Chlamydia trachomatis reference strain B/HAR36

The human pathogen Chlamydia trachomatis has multiple serovariants that have distinct organotropisms. We recently revised genomic sequence data linked to ocular reference strain, B/HAR36. Now linked to its correct genomic data in the European Nucleotide Archive, we describe its genomic features.

The GEA pipeline for characterizing Escherichia coli and Salmonella genomes

Salmonella enterica and Escherichia coli are major food-borne human pathogens, and their genomes are routinely sequenced for clinical surveillance. Computational pipelines designed for analyzing pathogen genomes should both utilize the most current information from annotation databases and increase the coverage of these databases over time. We report the development of the GEA pipeline to analyze large batches of E. coli and S. enterica genomes. The GEA pipeline takes as input paired Illumina raw reads files which are then assembled followed by annotation. Alternatively, assemblies can be provided as input and directly annotated. The pipeline provides predictive genome annotations for E. coli and S. enterica with a focus on the Center for Genomic Epidemiology tools. Annotation results are provided as a tab delimited text file. The GEA pipeline is designed for large-scale E. coli and S. enterica genome assembly and characterization using the Center for Genomic Epidemiology command-line tools and high-performance computing. Large scale annotation is demonstrated by an analysis of more than 14,000 Salmonella genome assemblies. Testing the GEA pipeline on E. coli raw reads demonstrates reproducibility across multiple compute environments and computational usage is optimized on high performance computers.

Readon: a novel algorithm to identify read-through transcripts with long-read sequencing data

MOTIVATION

There are many clustered transcriptionally active regions in the human genome, in which the transcription complex cannot immediately terminate transcription at the upstream gene termination site, but instead continues to transcribe intergenic regions and downstream genes, resulting in read-through transcripts. Several studies have demonstrated the regulatory roles of read-through transcripts in tumorigenesis and development. However, limited by the read length of next-generation sequencing, discovery of read-through transcripts has been slow. For long but also erroneous third-generation sequencing data, this study developed a novel minimizer sketch algorithm to accurately and quickly identify read-through transcripts.

RESULTS

Readon initially splits the reference sequence into distinct active regions. It employs a sliding window approach within each region, calculates minimizers, and constructs the specialized structured arrays for query indexing. Following initial alignment anchor screening of candidate read-through transcripts, further confirmation steps are executed. Comparative assessments against existing software reveal Readon's superior performance on both simulated and validated real data. Additionally, two downstream tools are provided: one for predicting whether a read-through transcript is likely to undergo nonsense-mediated decay or encodes a protein, and another for visualizing splicing patterns.

AVAILABILITY AND IMPLEMENTATION

Readon is freely available on GitHub (https://github.com/Bulabula45/Readon).

Transcriptomic investigation and biomarker discovery for zinc response in oysters Crassostrea gasar

In an era of unprecedented industrial and agricultural growth, metal contamination in marine environments is a pressing concern. Sentinel organisms such as the mangrove oyster Crassostrea gasar provide valuable insights into these environments' health. However, a comprehensive understanding of the molecular mechanisms underlying their response to metal exposure remains elusive. To address this gap, we reanalyzed the 454-sequencing data of C. gasar, utilizing an array of bioinformatics workflow of CDTA (Combined De Novo Transcriptome Assembly) to generate a more representative assembly. In parallel, C. gasar individuals were exposed to two concentrations of zinc (850 and 4500 μg L-1 Zn) for 48 h to understand their molecular responses. We utilized Trinotate workflow for the 11,684-CDTA unigenes annotation, with most transcripts aligning with the genus Crassostrea. Our analysis indicated that 67.3% of transcript sequences showed homology with Pfam, while 51.4% and 54.5%, respectively had GO and KO terms annotated. We identified potential metal pollution biomarkers, focusing on metal-related genes, such as those related to the GSH biosynthesis (CHAC1 and GCLC-like), to zinc transporters (ZNT2-like), and metallothionein (MT-like). The evolutionary conservation of these genes within the Crassostrea genus was assessed through phylogenetic analysis. Further, these genes were evaluated by qPCR in the laboratory exposed oysters. All target genes exhibited significant upregulation upon exposure to Zn at both 850 and 4500 μg L-1, except for GCLC-like, which showed upregulation only at the higher concentration of 4500 μg L-1. This result suggests distinct activation thresholds and complex interactions among these genes in response to varying Zn concentrations. Our study provides insights into the molecular responses of C. gasar to Zn, adding valuable tools for monitoring metal pollution in marine ecosystems using the mangrove oyster as a sentinel organism.

Synthesizing environmental, epidemiological and vector and parasite genetic data to assist decision making for disease elimination

We present a framework for identifying when conditions are favourable for transmission of vector-borne diseases between communities by incorporating predicted disease prevalence mapping with landscape analysis of sociological, environmental and host/parasite genetic data. We explored the relationship between environmental features and gene flow of a filarial parasite of humans, Onchocerca volvulus, and its vector, blackflies in the genus Simulium. We generated a baseline microfilarial prevalence map from point estimates from 47 locations in the ecological transition separating the savannah and forest in Ghana, where transmission of O. volvulus persists despite onchocerciasis control efforts. We generated movement suitability maps based on environmental correlates with mitochondrial population structure of 164 parasites from 15 communities and 93 vectors from only four sampling sites, and compared these to the baseline prevalence map. Parasite genetic distance between sampling locations was significantly associated with elevation (r = .793, p = .005) and soil moisture (r = .507, p = .002), while vector genetic distance was associated with soil moisture (r = .788, p = .0417) and precipitation (r = .835, p = .0417). The correlation between baseline prevalence and parasite resistance surface maps was stronger than that between prevalence and vector resistance surface maps. The centre of the study area had high prevalence and suitability for parasite and vector gene flow, potentially contributing to persistent transmission and suggesting the importance of re-evaluating transmission zone boundaries. With suitably dense sampling, this framework can help delineate transmission zones for onchocerciasis and would be translatable to other vector-borne diseases.

Molecular characterization and pathogenicity of a novel monopartite geminivirus infecting tobacco in China

The occurrence of geminiviruses causes significant economic losses in many economically important crops. In this study, a novel geminivirus isolated from tobacco in Sichuan province of China, named tomato leaf curl Chuxiong virus (TLCCxV), was characterized by small RNA-based deep sequencing. The full-length of TLCCxV genome was determined to be 2744 nucleotides (nt) encoding six open reading frames. Phylogenetic and genome-wide pairwise identity analysis revealed that TLCCxV shared less than 91% identities with reported geminiviruses. A TLCCxV infectious clone was constructed and successfully infected Nicotiana benthamiana, N. tabacum, N. glutinosa, Solanum lycopersicum and Petunia hybrida plants. Furthermore, expression of the V2, C1 and C4 proteins through a potato virus X vector caused severe chlorosis or necrosis symptom in N. benthamiana. Taken together, we identified a new geminivirus in tobacco plants, and found that V2, C1 and C4 contribute to symptom development.

Comparative Genomics of the First Resistant Candida auris Strain Isolated in Mexico: Phylogenomic and Pan-Genomic Analysis and Mutations Associated with Antifungal Resistance

Candida auris is an emerging multidrug-resistant and opportunistic pathogenic yeast. Whole-genome sequencing analysis has defined five major clades, each from a distinct geographic region. The current study aimed to examine the genome of the C. auris 20-1498 strain, which is the first isolate of this fungus identified in Mexico. Based on whole-genome sequencing, the draft genome was found to contain 70 contigs. It had a total genome size of 12.86 Mbp, an N50 value of 1.6 Mbp, and an average guanine-cytosine (GC) content of 45.5%. Genome annotation revealed a total of 5432 genes encoding 5515 proteins. According to the genomic analysis, the C. auris 20-1498 strain belongs to clade IV (containing strains endemic to South America). Of the two genes (ERG11 and FKS1) associated with drug resistance in C. auris, a mutation was detected in K143R, a gene located in a mutation hotspot of ERG11 (lanosterol 14-α-demethylase), an antifungal drug target. The focus on whole-genome sequencing and the identification of mutations linked to the drug resistance of fungi could lead to the discovery of new therapeutic targets and new antifungal compounds.

Exploring gene content with pangene graphs

Motivation

The gene content regulates the biology of an organism. It varies between species and between individuals of the same species. Although tools have been developed to identify gene content changes in bacterial genomes, none is applicable to collections of large eukaryotic genomes such as the human pangenome.

Results

We developed pangene, a computational tool to identify gene orientation, gene order and gene copy-number changes in a collection of genomes. Pangene aligns a set of input protein sequences to the genomes, resolves redundancies between protein sequences and constructs a gene graph with each genome represented as a walk in the graph. It additionally finds subgraphs, which we call bibubbles, that capture gene content changes. Applied to the human pangenome, pangene identifies known gene-level variations and reveals complex haplotypes that are not well studied before. Pangene also works with high-quality bacterial pangenome and reports similar numbers of core and accessory genes in comparison to existing tools.

Availability and implementation

Source code at https://github.com/lh3/pangene; pre-built pangene graphs can be downloaded from https://zenodo.org/records/8118576 and visualized at https://pangene.bioinweb.org.

Comparative Binding Study of Gliptins to Bacterial DPP4-like Enzymes for the Treatment of Type 2 Diabetes Mellitus (T2DM)

The role of the gut microbiota and its interplay with host metabolic health, particularly in the context of type 2 diabetes mellitus (T2DM) management, is garnering increasing attention. Dipeptidyl peptidase 4 (DPP4) inhibitors, commonly known as gliptins, constitute a class of drugs extensively used in T2DM treatment. However, their potential interactions with gut microbiota remain poorly understood. In this study, we employed computational methodologies to investigate the binding affinities of various gliptins to DPP4-like homologs produced by intestinal bacteria. The 3D structures of DPP4 homologs from gut microbiota species, including Segatella copri, Phocaeicola vulgatus, Bacteroides uniformis, Parabacteroides merdae, and Alistipes sp., were predicted using computational modeling techniques. Subsequently, molecular dynamics simulations were conducted for 200 ns to ensure the stability of the predicted structures. Stable structures were then utilized to predict the binding interactions with known gliptins through molecular docking algorithms. Our results revealed binding similarities of gliptins toward bacterial DPP4 homologs compared to human DPP4. Specifically, certain gliptins exhibited similar binding scores to bacterial DPP4 homologs as they did with human DPP4, suggesting a potential interaction of these drugs with gut microbiota. These findings could help in understanding the interplay between gliptins and gut microbiota DPP4 homologs, considering the intricate relationship between the host metabolism and microbial communities in the gut.

Global soil metagenomics reveals distribution and predominance of Deltaproteobacteria in nitrogen-fixing microbiome

BACKGROUND

Biological nitrogen fixation is a fundamental process sustaining all life on earth. While distribution and diversity of N2-fixing soil microbes have been investigated by numerous PCR amplicon sequencing of nitrogenase genes, their comprehensive understanding has been hindered by lack of de facto standard protocols for amplicon surveys and possible PCR biases. Here, by fully leveraging the planetary collections of soil shotgun metagenomes along with recently expanded culture collections, we evaluated the global distribution and diversity of terrestrial diazotrophic microbiome.

RESULTS

After the extensive analysis of 1,451 soil metagenomic samples, we revealed that the Anaeromyxobacteraceae and Geobacteraceae within Deltaproteobacteria are ubiquitous groups of diazotrophic microbiome in the soils with different geographic origins and land usage types, with particular predominance in anaerobic soils (paddy soils and sediments).

CONCLUSION

Our results indicate that Deltaproteobacteria is a core bacterial taxon in the potential soil nitrogen fixation population, especially in anaerobic environments, which encourages a careful consideration on deltaproteobacterial diazotrophs in understanding terrestrial nitrogen cycling. Video Abstract.

Characterization of CRISPR-Cas Systems in Shewanella algae and Shewanella haliotis: Insights into the Adaptation and Survival of Marine Pathogens

CRISPR-Cas systems are adaptive immune mechanisms present in most prokaryotes that play an important role in the adaptation of bacteria and archaea to new environments. Shewanella algae is a marine zoonotic pathogen with worldwide distribution, which accounts for the majority of clinical cases of Shewanella infections. However, the characterization of Shewanella algae CRISPR-Cas systems has not been well investigated yet. Through whole genome sequence analysis, we characterized the CRISPR-Cas systems in S. algae. Our results indicate that CRISPR-Cas systems are prevalent in S. algae, with the majority of strains containing the Type I-F system. This study provides new insights into the diversity and function of CRISPR-Cas systems in S. algae and highlights their potential role in the adaptation and survival of these marine pathogens.

Exploring the genetic makeup of Xanthomonas species causing bacterial spot in Taiwan: evidence of population shift and local adaptation

The introduction of plant pathogens can quickly reshape disease dynamics in island agro-ecologies, representing a continuous challenge for local crop management strategies. Xanthomonas pathogens causing tomato bacterial spot were probably introduced in Taiwan several decades ago, creating a unique opportunity to study the genetic makeup and adaptive response of this alien population. We examined the phenotypic and genotypic identity of 669 pathogen entries collected across different regions of Taiwan in the last three decades. The analysis detected a major population shift, where X. euvesicatoria and X. vesicatoria races T1 and T2 were replaced by new races of X. perforans. After its introduction, race T4 quickly became dominant in all tomato-growing areas of the island. The genomic analysis of 317 global genomes indicates that the Xanthomonas population in Taiwan has a narrow genetic background, most likely resulting from a small number of colonization events. However, despite the apparent genetic uniformity, X. perforans race T4 shows multiple phenotypic responses in tomato lines. Additionally, an in-depth analysis of effector composition suggests diversification in response to local adaptation. These include unique mutations on avrXv3 which might allow the pathogen to overcome Xv3/Rx4 resistance gene. The findings underscore the dynamic evolution of a pathogen when introduced in a semi-isolated environment and provide insights into the potential management strategies for this important disease of tomato.