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.

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.

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.

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.

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.

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.

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.

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.

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).

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.

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.

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.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Using short-read 16S rRNA sequencing of multiple variable regions to generate high-quality results to a species level

Introduction

Short-read amplicon sequencing studies have typically focused on 1-2 variable regions of the 16S rRNA gene. Species-level resolution is limited in these studies, as each variable region enables the characterisation of a different subsection of the microbiome. Although long-read sequencing techniques take advantage of all 9 variable regions by sequencing the entire 16S rRNA gene, they are substantially more expensive. This work assessed the feasibility of accurate species-level resolution and reproducibility using a relatively new sequencing kit and bioinformatics pipeline developed for short-read sequencing of multiple variable regions of the 16S rRNA gene. In addition, we evaluated the potential impact of different sample collection methods on our outcomes.

Methods

Using xGen™ 16S Amplicon Panel v2 kits, sequencing of all 9 variable regions of the 16S rRNA gene was carried out on an Illumina MiSeq platform. Mock cells and mock DNA for 8 bacterial species were included as extraction and sequencing controls respectively. Within-run and between-run replicate samples, and pairs of stool and rectal swabs collected at 0-5 weeks from the same participants, were incorporated. Observed relative abundances of each species were compared to theoretical abundances provided by ZymoBIOMICS. Paired Wilcoxon rank sum tests and distance-based intraclass correlation coefficients were used to statistically compare alpha and beta diversity measures, respectively, for pairs of replicates and stool/rectal swab sample pairs.

Results

Using multiple variable regions of the 16S ribosomal Ribonucleic Acid (rRNA) gene, we found that we could accurately identify taxa to a species level and obtain highly reproducible results at a species level. Yet, the microbial profiles of stool and rectal swab sample pairs differed substantially despite being collected concurrently from the same infants.

Conclusion

This protocol provides an effective means for studying infant gut microbial samples at a species level. However, sample collection approaches need to be accounted for in any downstream analysis.

Efficient and Robust Search of Microbial Genomes via Phylogenetic Compression

Comprehensive collections approaching millions of sequenced genomes have become central information sources in the life sciences. However, the rapid growth of these collections has made it effectively impossible to search these data using tools such as BLAST and its successors. Here, we present a technique called phylogenetic compression, which uses evolutionary history to guide compression and efficiently search large collections of microbial genomes using existing algorithms and data structures. We show that, when applied to modern diverse collections approaching millions of genomes, lossless phylogenetic compression improves the compression ratios of assemblies, de Bruijn graphs, and k -mer indexes by one to two orders of magnitude. Additionally, we develop a pipeline for a BLAST-like search over these phylogeny-compressed reference data, and demonstrate it can align genes, plasmids, or entire sequencing experiments against all sequenced bacteria until 2019 on ordinary desktop computers within a few hours. Phylogenetic compression has broad applications in computational biology and may provide a fundamental design principle for future genomics infrastructure.

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.

Proteogenomic Gene Structure Validation in the Pineapple Genome

MD2 pineapple (Ananas comosus) is the second most important tropical crop that preserves crassulacean acid metabolism (CAM), which has high water-use efficiency and is fast becoming the most consumed fresh fruit worldwide. Despite the significance of environmental efficiency and popularity, until very recently, its genome sequence has not been determined and a high-quality annotated proteome has not been available. Here, we have undertaken a pilot proteogenomic study, analyzing the proteome of MD2 pineapple leaves using liquid chromatography-mass spectrometry (LC-MS/MS), which validates 1781 predicted proteins in the annotated F153 (V3) genome. In addition, a further 603 peptide identifications are found that map exclusively to an independent MD2 transcriptome-derived database but are not found in the standard F153 (V3) annotated proteome. Peptide identifications derived from these MD2 transcripts are also cross-referenced to a more recent and complete MD2 genome annotation, resulting in 402 nonoverlapping peptides, which in turn support 30 high-quality gene candidates novel to both pineapple genomes. Many of the validated F153 (V3) genes are also supported by an independent proteomics data set collected for an ornamental pineapple variety. The contigs and peptides have been mapped to the current F153 genome build and are available as bed files to display a custom gene track on the Ensembl Plants region viewer. These analyses add to the knowledge of experimentally validated pineapple genes and demonstrate the utility of transcript-derived proteomics to discover both novel genes and genetic structure in a plant genome, adding value to its annotation.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BACKGROUND

Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.

RESULTS

We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.

CONCLUSIONS

The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

A chromosome-level genome assembly for the amphibious plant Rorippa aquatica reveals its allotetraploid origin and mechanisms of heterophylly upon submergence

The ability to respond to varying environments is crucial for sessile organisms such as plants. The amphibious plant Rorippa aquatica exhibits a striking type of phenotypic plasticity known as heterophylly, a phenomenon in which leaf form is altered in response to environmental factors. However, the underlying molecular mechanisms of heterophylly are yet to be fully understood. To uncover the genetic basis and analyze the evolutionary processes driving heterophylly in R. aquatica, we assembled the chromosome-level genome of the species. Comparative chromosome painting and chromosomal genomics revealed that allopolyploidization and subsequent post-polyploid descending dysploidy occurred during the speciation of R. aquatica. Based on the obtained genomic data, the transcriptome analyses revealed that ethylene signaling plays a central role in regulating heterophylly under submerged conditions, with blue light signaling acting as an attenuator of ethylene signal. The assembled R. aquatica reference genome provides insights into the molecular mechanisms and evolution of heterophylly.

Serratia montpellierensis sp. nov., Isolated from Laboratory-Reared Parasitic Wasps Psyttalia lounsburyii Silvestri and Psyttalia ponerophaga Silvestri (Hymenoptera: Braconidae)

Two strains of bacteria, PsyLou2AT and PsyPon4B, were isolated from adult braconid wasps Psyttalia lounsburyii and Psyttalia ponerophaga, respectively. These laboratory-reared wasps were investigated as agents for biological control of the olive fruit fly, Bactrocera oleae. Analysis of 16S rRNA genes of the two isolates demonstrated that they were highly related and belonged to the genus Serratia. Genomic sequencing of these isolates revealed genomes of 5,152,551 bp and 5,154,385 bp for PsyLou2AT and PsyPon4B, respectively, and both genomes had a mol% G+C content of 59.6%. Phylogenetic analyses using BLAST-based average nucleotide identity (ANIb), and digital DNA-DNA hybridization methods indicated that PsyLou2AT was most closely related to Serratia nevei S15T, producing ANIb and dDDH values of 96.11% and 70.2%, respectively. Since these values were literally on the species cutoff threshold, additional S. nevei genome assemblies were analyzed using ANIb and dDDH calculations. This revealed that among assemblies that were clearly identifiable as S. nevei, S. nevei S15T was the most closely related to PsyLou2AT, and that a majority of assemblies produced dDDH values of 68.3-68.7% relative to PsyLou2AT. Additionally, PsyLou2AT differed biochemically from S. nevei S15T in that it produced positive Voges Proskauer tests, produced protease, lacked arginine dihydrolase, and did not utilize D-lactose. Hence, PsyLou2AT represents a novel taxon within the Serratia, for which we propose the name Serratia montpellierensis sp. nov. The type strain is PsyLou2AT (=LMG 32817T =NRRL B-65689T).

A Combined Approach for the Characterization of Small Ruminant Lentivirus Strains Circulating in the Islands and Mainland of Greece

Small ruminant lentiviruses are a group of viruses infecting goat and sheep worldwide. These viruses exhibit an extraordinary degree of genetic and antigenic variability that severely influence in vivo and in vitro features, as well as diagnostic test results. Small ruminant farming is the most important animal farming business in Greece, with a high impact on the Greek primary economy. Although SRLV infection and its impact on animal production are well established in the country, little is known about the circulating SRLV strains and their prevalence. The aim of this study was to characterize SRLVs circulating in Greece with a combined serological and molecular approach, using the bulk milk matrix collected from 60 farms in different municipalities. This study allowed us to estimate a seroprevalence of around 52% at the herd level. The B1, B2 and A3 subtypes and a novel A viral cluster were identified. Moreover, the amplicon sequencing method allowed us to identify more than one viral subtype in a sample. These results again confirm the high variability of these viruses and highlight the importance of the constant monitoring of viral evolution, in particular in antigens of diagnostic interest.

Wasabi Gone Wild? Origin and Characterization of the Complete Plastomes of Ulleung Island Wasabi (Eutrema japonicum; Brassicaceae) and Other Cultivars in Korea

Korean wasabi occurs naturally on the young oceanic, volcanic Ulleung Island off the east coast of the Korean Peninsula. Although the Ulleung Island wasabi is reported as Eutrema japonicum and has been suggested to be morphologically identical to cultivars in Korea, very little is known about its taxonomic identity and relationship with other cultivars. In this study, we sequenced the complete chloroplast DNA sequences of three naturally occurring Ulleung Island wasabi plants and six cultivars ('Daewang', 'Daruma', 'Micado', 'Orochi', 'Green Thumb', and 'Shogun') from continental Korea and determined the taxonomic identity of Korean wasabi on Ulleung Island. The size and organization of the complete chloroplast genomes of the nine accessions were nearly identical to those of previously reported wasabi cultivars. In addition, phylogenetic analysis based on the complete plastomes suggested that Ulleung Island wasabi most likely comprises various wasabi cultivars with three chlorotypes ('Shogun', 'Green Thumb', and a unique Chusan type). Based on the complete plastomes, we identified eight chlorotypes for the major wasabi cultivars and the Ulleung Island wasabi. Two major groups (1-'Mazuma' and 'Daruma', and 2-'Fujidaruma'/'Shimane No. 3'/Ulleung Island wasabi/five cultivars in Korea) were also identified based on mother line genealogical history. Furthermore, different types of variations (mutations, insertions/deletions (indels), mononucleotide repeats, and inversions) in plastomes were identified to distinguish different cultivar lines and five highly divergent hotspots. The nine newly obtained complete plastomes are valuable organelle genomic resources for species identification and infraspecific phylogeographic studies on wild and cultivated wasabi.

GeneMiner: A tool for extracting phylogenetic markers from next-generation sequencing data

The advancement of next-generation sequencing (NGS) technologies has been revolutionary for the field of evolutionary biology. This technology has led to an abundance of available genomes and transcriptomes for researchers to mine. Specifically, researchers can mine for various types of molecular markers that are vital for phylogenetic, evolutionary and ecological studies. Numerous tools have been developed to extract these molecular markers from NGS data. However, due to an insufficient number of well-annotated reference genomes for non-model organisms, it remains challenging to obtain these markers accurately and efficiently. Here, we present GeneMiner, an improved and expanded version of our previous tool, Easy353. GeneMiner combines the reference-guided de Bruijn graph assembly with seed self-discovery and greedy extension. Additionally, it includes a verification step using a parameter-bootstrap method to reduce the pitfalls associated with using a relatively distant reference. Our results, using both experimental and simulation data, showed GeneMiner can accurately acquire phylogenetic molecular markers for plants using transcriptomic, genomic and other NGS data. GeneMiner is designed to be user-friendly, fast and memory-efficient. Further, it is compatible with Linux, Windows and macOS. All source codes are publicly available on GitHub (https://github.com/sculab/GeneMiner) and Gitee (https://gitee.com/sculab/GeneMiner) for easy accessibility and transparency.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae)

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

Complete organelle genomes of Korean fir, Abies koreana and phylogenomics of the gymnosperm genus Abies using nuclear and cytoplasmic DNA sequence data

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus.

Correlation between the gut microbiome and neurodegenerative diseases: a review of metagenomics evidence

A growing body of evidence suggests that the gut microbiota contributes to the development of neurodegenerative diseases via the microbiota-gut-brain axis. As a contributing factor, microbiota dysbiosis always occurs in pathological changes of neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. High-throughput sequencing technology has helped to reveal that the bidirectional communication between the central nervous system and the enteric nervous system is facilitated by the microbiota's diverse microorganisms, and for both neuroimmune and neuroendocrine systems. Here, we summarize the bioinformatics analysis and wet-biology validation for the gut metagenomics in neurodegenerative diseases, with an emphasis on multi-omics studies and the gut virome. The pathogen-associated signaling biomarkers for identifying brain disorders and potential therapeutic targets are also elucidated. Finally, we discuss the role of diet, prebiotics, probiotics, postbiotics and exercise interventions in remodeling the microbiome and reducing the symptoms of neurodegenerative diseases.

In Vitro Profiling of Potential Fish Probiotics, Enterococcus hirae Strains, Isolated from Jade Perch, and Safety Properties Assessed Using Whole Genome Sequencing

The demands of intensified aquaculture production and escalating disease prevalence underscore the need for efficacious probiotic strategies to enhance fish health. This study focused on isolating and characterising potential probiotics from the gut microbiota of the emerging aquaculture species jade perch (Scortum barcoo). Eighty-seven lactic acid bacteria and 149 other bacteria were isolated from the digestive tract of five adult jade perch. The screening revealed that 24 Enterococcus hirae isolates inhibited the freshwater pathogens Aeromonas sobria and Streptococcus iniae. Co-incubating E. hirae with the host gut suspensions demonstrated a two- to five-fold increase in the size of growth inhibition zones compared to the results when using gut suspensions from tilapia (a non-host), indicating host-specificity. Genome analysis of the lead isolate, E. hirae R44, predicted the presence of antimicrobial compounds like enterolysin A, class II lanthipeptide, and terpenes, which underlay its antibacterial attributes. Isolate R44 exhibited desirable probiotic characteristics, including survival at pH values within the range of 3 to 12, bile tolerance, antioxidant activity, ampicillin sensitivity, and absence of transferable antimicrobial resistance genes and virulence factors commonly associated with hospital Enterococcus strains (IS16, hylEfm, and esp). This study offers a foundation for sourcing host-adapted probiotics from underexplored aquaculture species. Characterisation of novel probiotics like E. hirae R44 can expedite the development of disease mitigation strategies to support aquaculture intensification.

Identification of multiple transcription factor genes potentially involved in the development of electrosensory versus mechanosensory lateral line organs

In electroreceptive jawed vertebrates, embryonic lateral line placodes give rise to electrosensory ampullary organs as well as mechanosensory neuromasts. Previous reports of shared gene expression suggest that conserved mechanisms underlie electroreceptor and mechanosensory hair cell development and that electroreceptors evolved as a transcriptionally related "sister cell type" to hair cells. We previously identified only one transcription factor gene, Neurod4, as ampullary organ-restricted in the developing lateral line system of a chondrostean ray-finned fish, the Mississippi paddlefish (Polyodon spathula). The other 16 transcription factor genes we previously validated in paddlefish were expressed in both ampullary organs and neuromasts. Here, we used our published lateral line organ-enriched gene-set (arising from differential bulk RNA-seq in late-larval paddlefish), together with a candidate gene approach, to identify 25 transcription factor genes expressed in the developing lateral line system of a more experimentally tractable chondrostean, the sterlet (Acipenser ruthenus, a small sturgeon), and/or that of paddlefish. Thirteen are expressed in both ampullary organs and neuromasts, consistent with conservation of molecular mechanisms. Seven are electrosensory-restricted on the head (Irx5, Irx3, Insm1, Sp5, Satb2, Mafa and Rorc), and five are the first-reported mechanosensory-restricted transcription factor genes (Foxg1, Sox8, Isl1, Hmx2 and Rorb). However, as previously reported, Sox8 is expressed in ampullary organs as well as neuromasts in a catshark (Scyliorhinus canicula), suggesting the existence of lineage-specific differences between cartilaginous and ray-finned fishes. Overall, our results support the hypothesis that ampullary organs and neuromasts develop via largely conserved transcriptional mechanisms, and identify multiple transcription factors potentially involved in the formation of electrosensory versus mechanosensory lateral line organs.

Genomic diversity of the pathogenic fungus Aspergillus fumigatus in Japan reveals the complex genomic basis of azole resistance

Aspergillus fumigatus is a pathogenic fungus with a global distribution. The emergence of azole-resistant A. fumigatus (ARAf) other than the TR-mutants is a problem in Japan. Additionally, the genetic diversity of A. fumigatus strains in Japan remains relatively unknown. Here we show the diversity in the A. fumigatus strains isolated in Japan as well as the complexity in the global distribution of the pathogenic strains. First, we analyzed the genome sequences of 171 strains from Japan as well as the antifungal susceptibility of these strains. Next, we conducted a population analysis of 876 strains by combining the available genomic data for strains isolated worldwide, which were grouped in six clusters. Finally, a genome-wide association study identified the genomic loci associated with ARAf strains, but not the TR-mutants. These results highlight the complexity of the genomic mechanism underlying the emergence of ARAf strains other than the TR-mutants.

Inflammation and bacteriophages affect DNA inversion states and functionality of the gut microbiota

Reversible genomic DNA inversions control the expression of numerous gut bacterial molecules, but how this impacts disease remains uncertain. By analyzing metagenomic samples from inflammatory bowel disease (IBD) cohorts, we identified multiple invertible regions where a particular orientation correlated with disease. These include the promoter of polysaccharide A (PSA) of Bacteroides fragilis, which induces regulatory T cells (Tregs) and ameliorates experimental colitis. The PSA promoter was mostly oriented "OFF" in IBD patients, which correlated with increased B. fragilis-associated bacteriophages. Similarly, in mice colonized with a healthy human microbiota and B. fragilis, induction of colitis caused a decline of PSA in the "ON" orientation that reversed as inflammation resolved. Monocolonization of mice with B. fragilis revealed that bacteriophage infection increased the frequency of PSA in the "OFF" orientation, causing reduced PSA expression and decreased Treg cells. Altogether, we reveal dynamic bacterial phase variations driven by bacteriophages and host inflammation, signifying bacterial functional plasticity during disease.

Identification and Characterisation of a Myxoma Virus Detected in the Italian Hare (Lepus corsicanus)

Myxoma virus (MYXV) is a Leporipoxvirus (genus) belonging to the family Poxviridae; it is characterised by a genome of approximately 161 kb dsDNA encoding for several proteins that play an essential role in both host spectrum determination and immunomodulation. The healthy reservoir of the virus is Sylvilagus spp. At the same time, in wild and domestic European rabbits (Oryctolagus cuniculus), MYXV is the etiologic agent of myxomatosis, a disease with an extremely high mortality rate. In 2014, an interspecies jump of MYXV was reported in Lepus europaeus in the UK. In 2018, myxomatosis induced by a new recombinant strain called MYXV-To was identified during a large outbreak in Iberian hares (Lepus granatensis) in Spain. Here, we describe the case of myxomatosis in another hare species: an adult male Italian hare (Lepus corsicanus) found dead in 2018 in Sicily with lesions suggestive of myxomatosis and treponema infection. Laboratory tests, e.g., end-point PCR and negative staining electron microscopy, confirmed the presence of both pathogens. MYXV was then isolated from tissue samples in permissive cells and sequenced using NGS technology. Main genomic differences concerning known MYXV strains are discussed.

MetaCompass: Reference-guided Assembly of Metagenomes

Metagenomic studies have primarily relied on de novo assembly for reconstructing genes and genomes from microbial mixtures. While reference-guided approaches have been employed in the assembly of single organisms, they have not been used in a metagenomic context. Here we describe the first effective approach for reference-guided metagenomic assembly that can complement and improve upon de novo metagenomic assembly methods for certain organisms. Such approaches will be increasingly useful as more genomes are sequenced and made publicly available.

Rapid adaptive substitution of L226Q in HA protein increases the pathogenicity of H9N2 viruses in mice

Background

Since the first human infection with H9N2 virus was reported in 1998, the number of cases of H9N2 infection has exceeded one hundred by 2021. However, there is no systematic description of the biological characteristics of H9N2 viruses isolated from humans.

Methods

Therefore, this study analyzed the pathogenicity in mice of all available H9N2 viruses isolated from human cases in China from 2013 to 2021.

Results

Although most of the H9N2 viruses analyzed showed low or no pathogenicity in mice, the leucine to glutamine substitution at residue 226 (L226Q) in the hemagglutinin (HA) protein rapidly emerged during the adaptation of H9N2 viruses, and was responsible for severe infections and even fatalities. HA amino acid 226Q conferred a remarkable competitive advantage on H9N2 viruses in mice relative to viruses containing 226L, increasing their virulence, infectivity, and replication.

Conclusion

Thus, our study demonstrates that the adaptive substitution HA L226Q rapidly acquired by H9N2 viruses during the course of infection in mice contributed to their high pathogenicity.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.