tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

Comparative analyses of Pleurotus pulmonarius mitochondrial genomes reveal two major lineages of mini oyster mushroom cultivars

Pleurotus pulmonarius, commonly known as the mini oyster mushroom, is highly esteemed for its crisp texture and umami flavor. Limited genetic diversity among P. pulmonarius cultivars raises concerns regarding its sustainable industrial production. To delve into the maternal genetic diversity of the principal P. pulmonarius cultivars, 36 cultivars and five wild isolates were subjected to de novo sequencing and assembly to generate high-quality mitogenome sequences. The P. pulmonarius mitogenomes had lengths ranging from 69,096 to 72,905 base pairs. The mitogenome sizes of P. pulmonarius and those of other mushroom species in the Pleurotus genus showed a significant positive correlation with the counts of LAGLIDAG and GIY-YIG homing endonucleases encoded by intronic open reading frames. A comparison of gene arrangements revealed an inversion of a fragment containing atp9-nad3-nad2 between P. pulmonarius and P. ostreatus. The mitogenomes of P. pulmonarius were clustered into three distinct clades, two of which were crowded with commercial cultivars. Clade I, all of which possess an inserted dpo gene, shared a maternal origin linked to an ancestral cultivar from Taiwan. Primers were designed to target the dpo gene, potentially safeguarding intellectual property rights. The wild isolates in Clade III exhibited more divergent mitogenomes, rendering them valuable for breeding.

Comparative genomic and metabolomic analysis reveals the potential of a newly isolated Enterococcus faecium B6 involved in lipogenic effects

Evidence has indicated that Enterococcus plays a vital role in non-alcoholic fatty liver disease (NAFLD) development. However, the microbial genetic basis and metabolic potential in the disease are yet unknown. We previously isolated a bacteria Enterococcus faecium B6 (E. faecium B6) from children with NAFLD for the first time. Here, we aim to systematically investigate the potential of strain B6 in lipogenic effects. The lipogenic effects of strain B6 were explored in vitro and in vivo. The genomic and functional characterizations were investigated by whole-genome sequencing and comparative genomic analysis. Moreover, the metabolite profiles were unraveled by an untargeted metabolomic analysis. We demonstrated that strain B6 could effectively induce lipogenic effects in the liver of mice. Strain B6 contained a circular chromosome and two circular plasmids and posed various functions. Compared to the other two probiotic strains of E. faecium, strain B6 exhibited unique functions in pathways of ABC transporters, phosphotransferase system, and amino sugar and nucleotide sugar metabolism. Moreover, strain B6 produced several metabolites, mainly enriched in the protein digestion and absorption pathway. The unique potential of strain B6 in lipogenic effects was probably associated with glycolysis, fatty acid synthesis, and glutamine and choline transport. This study pioneeringly revealed the metabolic characteristics and specific detrimental traits of strain B6. The findings provided new insights into the underlying mechanisms of E. faecium in lipogenic effects, and laid essential foundations for further understanding of E. faecium-related disease.

Integrative genome-centric metagenomics for surface water surveillance: Elucidating microbiomes, antimicrobial resistance, and their associations

Surface water ecosystems are intimately intertwined with anthropogenic activities and have significant public health implications as primary sources of irrigation water in agricultural production. Our extensive metagenomic analysis examined 404 surface water samples from four different geological regions in Chile and Brazil, spanning irrigation canals (n = 135), rivers (n = 121), creeks (n = 74), reservoirs (n = 66), and ponds (n = 8). Overall, 50.25 % of the surface water samples contained at least one of the pathogenic or contaminant bacterial genera (Salmonella: 29.21 %; Listeria: 6.19 %; Escherichia: 35.64 %). Furthermore, a total of 1,582 antimicrobial resistance (AMR) gene clusters encoding resistance to 25 antimicrobial classes were identified, with samples from Brazil exhibiting an elevated AMR burden. Samples from stagnant water sources were characterized by dominant Cyanobacteriota populations, resulting in significantly reduced biodiversity and more uniform community compositions. A significant association between taxonomic composition and the resistome was supported by a Procrustes analysis (p < 0.001). Notably, regional signatures were observed regarding the taxonomic and resistome profiles, as samples from the same region clustered together on both ordinates. Additionally, network analysis illuminated the intricate links between taxonomy and AMR at the contig level. Our deep sequencing efforts not only mapped the microbial landscape but also expanded the genomic catalog with newly characterized metagenome-assembled genomes (MAGs), boosting the classification of reads by 12.85 %. In conclusion, this study underscores the value of metagenomic approaches in surveillance of surface waters, enhancing our understanding of microbial and AMR dynamics with far-reaching public health and ecological ramifications.

Differential adaptation of the yeast Candida anglica to fermented food

A single strain of Candida anglica, isolated from cider, is available in international yeast collections. We present here seven new strains isolated from French PDO cheeses. For one of the cheese strains, we achieved a high-quality genome assembly of 13.7 Mb with eight near-complete telomere-to-telomere chromosomes. The genomes of two additional cheese strains and of the cider strain were also assembled and annotated, resulting in a core genome of 5966 coding sequences. Phylogenetic analysis showed that the seven cheese strains clustered together, away from the cider strain. Mating-type locus analysis revealed the presence of a MATa locus in the cider strain but a MATalpha locus in all cheese strains. The presence of LINE retrotransposons at identical genome position in the cheese strains, and two different karyotypic profiles resulting from chromosomal rearrangements were observed. Together, these findings are consistent with clonal propagation of the cheese strains. Phenotypic trait variations were observed within the cheese population under stress conditions whereas the cider strain was found to have a much greater capacity for growth in all conditions tested.

The whole genome dataset of Ichthyscopus pollicaris

The classification of the Uranoscopidae species is controversial and the Ichthyscopus pollicaris belonging to Uranoscopidae was first reported in 2019. In the present study, the whole genome sequence of I. pollicaris were generated by PacBio and Illumina platforms for the first time. After de novo assembly and correction of the high-quality PacBio data, a 527.25 Mb I. pollicaris genome with an N50 length of 11.25 Mb was finally generated. Meanwhile, 170.41 Mb repeating sequence, 21,263 genes, 784 miRNAs, 2,225 tRNAs, 3004 rRNAs, and 1422 snRNAs were annotated in I. pollicaris genome. Furthermore, 3,168 single-copy orthologous genes were applied to reconstructed the phylogenetic relationship between I. pollicaris and other 11 species. The draft genome sequences have been deposited in NCBI database with the accession number of PRJNA1071810.

Combining Nanopore direct RNA sequencing with genetics and mass spectrometry for analysis of T-loop base modifications across 42 yeast tRNA isoacceptors

Transfer RNAs (tRNAs) contain dozens of chemical modifications. These modifications are critical for maintaining tRNA tertiary structure and optimizing protein synthesis. Here we advance the use of Nanopore direct RNA-sequencing (DRS) to investigate the synergy between modifications that are known to stabilize tRNA structure. We sequenced the 42 cytosolic tRNA isoacceptors from wild-type yeast and five tRNA-modifying enzyme knockout mutants. These data permitted comprehensive analysis of three neighboring and conserved modifications in T-loops: 5-methyluridine (m5U54), pseudouridine (Ψ55), and 1-methyladenosine (m1A58). Our results were validated using direct measurements of chemical modifications by mass spectrometry. We observed concerted T-loop modification circuits-the potent influence of Ψ55 for subsequent m1A58 modification on more tRNA isoacceptors than previously observed. Growing cells under nutrient depleted conditions also revealed a novel condition-specific increase in m1A58 modification on some tRNAs. A global and isoacceptor-specific classification strategy was developed to predict the status of T-loop modifications from a user-input tRNA DRS dataset, applicable to other conditions and tRNAs in other organisms. These advancements demonstrate how orthogonal technologies combined with genetics enable precise detection of modification landscapes of individual, full-length tRNAs, at transcriptome-scale.

Complete mitochondrial genome assembly of Juglans regia unveiled its molecular characteristics, genome evolution, and phylogenetic implications

BACKGROUND

The Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology.

RESULTS

The 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress.

CONCLUSIONS

Our results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species.

Exploring microbial diversity and biosynthetic potential in zoo and wildlife animal microbiomes

Understanding human, animal, and environmental microbiota is essential for advancing global health and combating antimicrobial resistance (AMR). We investigate the oral and gut microbiota of 48 animal species in captivity, comparing them to those of wildlife animals. Specifically, we characterize the microbiota composition, metabolic pathways, AMR genes, and biosynthetic gene clusters (BGCs) encoding the production of specialized metabolites. Our results reveal a high diversity of microbiota, with 585 novel species-level genome bins (SGBs) and 484 complete BGCs identified. Functional gene analysis of microbiomes shows diet-dependent variations. Furthermore, by comparing our findings to wildlife-derived microbiomes, we observe the impact of captivity on the animal microbiome, including examples of converging microbiome compositions. Importantly, our study identifies AMR genes against commonly used veterinary antibiotics, as well as resistance to vancomycin, a critical antibiotic in human medicine. These findings underscore the importance of the 'One Health' approach and the potential for zoonotic transmission of pathogenic bacteria and AMR. Overall, our study contributes to a better understanding of the complexity of the animal microbiome and highlights its BGC diversity relevant to the discovery of novel antimicrobial compounds.

3t-seq: automatic gene expression analysis of single-copy genes, transposable elements, and tRNAs from RNA-seq data

RNA sequencing is the gold-standard method to quantify transcriptomic changes between two conditions. The overwhelming majority of data analysis methods available are focused on polyadenylated RNA transcribed from single-copy genes and overlook transcripts from repeated sequences such as transposable elements (TEs). These self-autonomous genetic elements are increasingly studied, and specialized tools designed to handle multimapping sequencing reads are available. Transfer RNAs are transcribed by RNA polymerase III and are essential for protein translation. There is a need for integrated software that is able to analyze multiple types of RNA. Here, we present 3t-seq, a Snakemake pipeline for integrated differential expression analysis of transcripts from single-copy genes, TEs, and tRNA. 3t-seq produces an accessible report and easy-to-use results for downstream analysis starting from raw sequencing data and performing quality control, genome mapping, gene expression quantification, and statistical testing. It implements three methods to quantify TEs expression and one for tRNA genes. It provides an easy-to-configure method to manage software dependencies that lets the user focus on results. 3t-seq is released under MIT license and is available at https://github.com/boulardlab/3t-seq.

Long-read de novo genome assembly of Gulf toadfish (Opsanus beta)

BACKGROUND

The family Batrachoididae are a group of ecologically important teleost fishes with unique life histories, behavior, and physiology that has made them popular model organisms. Batrachoididae remain understudied in the realm of genomics, with only four reference genome assemblies available for the family, with three being highly fragmented and not up to current assembly standards. Among these is the Gulf toadfish, Opsanus beta, a model organism for serotonin physiology which has recently been bred in captivity.

RESULTS

Here we present a new, de novo genome and transcriptome assemblies for the Gulf toadfish using PacBio long read technology. The genome size of the final assembly is 2.1 gigabases, which is among the largest teleost genomes. This new assembly improves significantly upon the currently available reference for Opsanus beta with a final scaffold count of 62, of which 23 are chromosome scale, an N50 of 98,402,768, and a BUSCO completeness score of 97.3%. Annotation with ab initio and transcriptome-based methods generated 41,076 gene models. The genome is highly repetitive, with ~ 70% of the genome composed of simple repeats and transposable elements. Satellite DNA analysis identified potential telomeric and centromeric regions.

CONCLUSIONS

This improved assembly represents a valuable resource for future research using this important model organism and to teleost genomics more broadly.

Haplotype-resolved genome assembly of the upas tree (Antiaris toxicaria)

The upas tree (Antiaris toxicaria Lesch.) is a medically important plant that contains various specialized metabolites with significant bioactivity. The lack of a reference genome hinders the in-depth study as well as rational exploitation and conservation of this plant. Here, we present the first holotype-resolved chromosome-scale genome of the upas tree. The assembled genome consisted of 26 chromosomes that contain 1.34 Gb of sequencing data with a contig N50 length of 60 Mb. Genome annotation identified 43,500 protein-coding genes in the upas tree genome, of which 98.75% were functionally annotated. This high-quality reference genome will lay the foundation for further studies on the evolution and functional genomics of the upas tree.

From raw milk cheese to the gut: investigating the colonization strategies of Bifidobacterium mongoliense

The microbial ecology of raw milk cheeses is determined by bacteria originating from milk and milk-producing animals. Recently, it has been shown that members of the Bifidobacterium mongoliense species may become transmitted along the Parmigiano Reggiano cheese production chain and ultimately may colonize the consumer intestine. However, there is a lack of knowledge regarding the molecular mechanisms that mediate the interaction between B. mongoliense and the human gut. Based on 128 raw milk cheeses collected from different Italian regions, we isolated and characterized 10 B. mongoliense strains. Comparative genomics allowed us to unveil the presence of enzymes required for the degradation of sialylated host-glycans in B. mongoliense, corroborating the appreciable growth on de Man-Rogosa-Sharpe (MRS) medium supplemented with 3'-sialyllactose (3'-SL) or 6'-sialyllactose (6'-SL). The B. mongoliense BMONG18 was chosen, due to its superior ability to utilize 3'-SL and mucin as representative strain, to investigate its behavior when co-inoculated with other bifidobacterial species. Conversely, members of other bifidobacterial species did not appear to benefit from the presence of BMONG18, highlighting a competitive scenario for nutrient acquisition. Transcriptomic data of BMONG18 reveal no significant differences in gene expression when cultivated in a gut simulating medium (GSM), regardless of whether cheese was included or not. Furthermore, BMONG18 was shown to exhibit high adhesion capabilities to HT29-MTX human cells, in line with its colonization ability of a human host.IMPORTANCEFermented foods are nourishments produced through controlled microbial growth that play an essential role in worldwide human nutrition. Research interest in fermented foods has increased since the 80s, driven by growing awareness of their potential health benefits beyond mere nutritional content. Bifidobacterium mongoliense, previously identified throughout the production process of Parmigiano Reggiano cheese, was found to be capable of establishing itself in the intestines of its consumers. Our study underscores molecular mechanisms through which this bifidobacterial species, derived from food, interacts with the host and other gut microbiota members.

The genome of Lupinus mutabilis: Evolution and genetics of an emerging bio-based crop

Lupinus mutabilis is an under-domesticated legume species from the Andean region of South America. It belongs to the New World lupins clade, which groups several lupin species displaying large genetic variation and adaptability to highly different environments. L. mutabilis is attracting interest as a potential multipurpose crop to diversify the European supply of plant proteins, increase agricultural biodiversity, and fulfill bio-based applications. This study reports the first high-quality L. mutabilis genome assembly, which is also the first sequenced assembly of a New World lupin species. Through comparative genomics and phylogenetics, the evolution of L. mutabilis within legumes and lupins is described, highlighting both genomic similarities and patterns specific to L. mutabilis, potentially linked to environmental adaptations. Furthermore, the assembly was used to study the genetics underlying important traits for the establishment of L. mutabilis as a novel crop, including protein and quinolizidine alkaloids contents in seeds, genomic patterns of classic resistance genes, and genomic properties of L. mutabilis mycorrhiza-related genes. These analyses pointed out copy number variation, differential genomic gene contexts, and gene family expansion through tandem duplications as likely important drivers of the genomic diversity observed for these traits between L. mutabilis and other lupins and legumes. Overall, the L. mutabilis genome assembly will be a valuable resource to conduct genetic research and enable genomic-based breeding approaches to turn L. mutabilis into a multipurpose legume crop.

Isolation, characterization, and potential application of Acinetobacter baumannii phages against extensively drug-resistant strains

One of the significant issues in treating bacterial infections is the increasing prevalence of extensively drug-resistant (XDR) strains of Acinetobacter baumannii. In the face of limited or no viable treatment options for extensively drug-resistant (XDR) bacteria, there is a renewed interest in utilizing bacteriophages as a treatment option. Three Acinetobacter phages (vB_AbaS_Ftm, vB_AbaS_Eva, and vB_AbaS_Gln) were identified from hospital sewage and analyzed for their morphology, host ranges, and their genome sequences were determined and annotated. These phages and vB_AbaS_SA1 were combined to form a phage cocktail. The antibacterial effects of this cocktail and its combinations with selected antimicrobial agents were evaluated against the XDR A. baumannii strains. The phages exhibited siphovirus morphology. Out of a total of 30 XDR A. baumannii isolates, 33% were sensitive to vB_AbaS_Ftm, 30% to vB_AbaS_Gln, and 16.66% to vB_AbaS_Eva. When these phages were combined with antibiotics, they demonstrated a synergistic effect. The genome sizes of vB_AbaS_Ftm, vB_AbaS_Eva, and vB_AbaS_Gln were 48487, 50174, and 50043 base pairs (bp), respectively, and showed high similarity. Phage cocktail, when combined with antibiotics, showed synergistic effects on extensively drug-resistant (XDR) strains of A. baumannii. However, the need for further study to fully understand the mechanisms of action and potential limitations of using these phages is highlighted.

Complete genome sequence of Pseudomonas aeruginosa Y010, a taro rhizosphere strain producing potent antimicrobial agents

Pseudomonas aeruginosa Y010, isolated from the taro rhizosphere, exhibits great antagonistic abilities against Dickeya strains that cause soft-rot and blackleg diseases of plants by producing potent antimicrobial agents. The complete genome of Y010 was sequenced and annotated, which is 6,415,628 bp in length with 66.39% GC content.

Characterisation of the complete chloroplast genome of Solanum tuberosum cv. White Lady

Potato (Solanum tuberosum) is considered worldwide as one of the most important non-cereal food crops. As a result of its adaptability and worldwide production area, potato displays a vast phenotypical variability as well as genomic diversity. Chloroplast genomes have long been a core issue in plant molecular evolution and phylogenetic studies, and have an important role in revealing photosynthetic mechanisms, metabolic regulations and the adaptive evolution of plants. We sequenced the complete chloroplast genome of the Hungarian cultivar White Lady, which is 155 549 base pairs (bp) in length and is characterised by the typical quadripartite structure composed of a large- and small single-copy region (85 991 bp and 18 374 bp, respectively) interspersed by two identical inverted repeats (25 592 bp). The genome consists of 127 genes of which 82 are protein-coding, eight are ribosomal RNAs and 37 are transfer RNAs. The overall gene content and distribution of the genes on the White Lady chloroplast was the same as found in other potato chloroplasts. The alignment of S. tuberosum chloroplast genome sequences resulted in a highly resolved tree, with 10 out of the 13 nodes recovered having bootstrap values over 90%. By comparing the White Lady chloroplast genome with available S. tuberosum sequences we found that gene content and synteny are highly conserved. The new chloroplast sequence can support further studies of genetic diversity, resource conservation, evolution and applied agricultural research. The new sequence can support further potato genetic diversity and evolutionary studies, resource conservation, and also applied agricultural research.

The complete mitochondrial genome of the banana pathogen Fusarium oxysporum f. sp. cubense M5

We report the complete mitochondrial genome of a causal agent of banana fusarium wilt isolated in Mexico. The whole set of genes encoding proteins related to respiration and ATP synthesis, rRNAs, tRNAs are enlisted. Two open reading frames of unknown function conserved in Fusarium oxysporum were also identified.

Integrated analysis of facial microbiome and skin physio-optical properties unveils cutotype-dependent aging effects

BACKGROUND

Our facial skin hosts millions of microorganisms, primarily bacteria, crucial for skin health by maintaining the physical barrier, modulating immune response, and metabolizing bioactive materials. Aging significantly influences the composition and function of the facial microbiome, impacting skin immunity, hydration, and inflammation, highlighting potential avenues for interventions targeting aging-related facial microbes amidst changes in skin physiological properties.

RESULTS

We conducted a multi-center and deep sequencing survey to investigate the intricate interplay of aging, skin physio-optical conditions, and facial microbiome. Leveraging a newly-generated dataset of 2737 species-level metagenome-assembled genomes (MAGs), our integrative analysis highlighted aging as the primary driver, influencing both facial microbiome composition and key skin characteristics, including moisture, sebum production, gloss, pH, elasticity, and sensitivity. Further mediation analysis revealed that skin characteristics significantly impacted the microbiome, mostly as a mediator of aging. Utilizing this dataset, we uncovered two consistent cutotypes across sampling cities and identified aging-related microbial MAGs. Additionally, a Facial Aging Index (FAI) was formulated based on the microbiome, uncovering the cutotype-dependent effects of unhealthy lifestyles on skin aging. Finally, we distinguished aging related microbial pathways influenced by lifestyles with cutotype-dependent effect.

CONCLUSIONS

Together, our findings emphasize aging's central role in facial microbiome dynamics, and support personalized skin microbiome interventions by targeting lifestyle, skin properties, and aging-related microbial factors. Video Abstract.

Comparative Genomics of Fungi in Nectriaceae Reveals Their Environmental Adaptation and Conservation Strategies

This study presents the first genome assembly of the freshwater saprobe fungus Neonectria lugdunensis and a comprehensive phylogenomics analysis of the Nectriaceae family, examining genomic traits according to fungal lifestyles. The Nectriaceae family, one of the largest in Hypocreales, includes fungi with significant ecological roles and economic importance as plant pathogens, endophytes, and saprobes. The phylogenomics analysis identified 2684 single-copy orthologs, providing a robust evolutionary framework for the Nectriaceae family. We analyzed the genomic characteristics of 17 Nectriaceae genomes, focusing on their carbohydrate-active enzymes (CAZymes), biosynthetic gene clusters (BGCs), and adaptations to environmental temperatures. Our results highlight the adaptation mechanisms of N. lugdunensis, emphasizing its capabilities for plant litter degradation and enzyme activity in varying temperatures. The comparative genomics of different Nectriaceae lifestyles revealed significant differences in genome size, gene content, repetitive elements, and secondary metabolite production. Endophytes exhibited larger genomes, more effector proteins, and BGCs, while plant pathogens had higher thermo-adapted protein counts, suggesting greater resilience to global warming. In contrast, the freshwater saprobe shows less adaptation to warmer temperatures and is important for conservation goals. This study underscores the importance of understanding fungal genomic adaptations to predict ecosystem impacts and conservation targets in the face of climate change.

The assembly and comparative analysis of the first complete mitogenome of Lindera aggregata

Lindera aggregata, a member belongs to the genus Lindera of Lauraceae family. Its roots and leaves have been used as traditional Chinese medicine or functional food for thousands of years. However, its mitochondrial genome has not been explored. Our aim is to sequence and assemble the mitogenome of L. aggregata to elucidate the genetic mechanism and evolutionary pathway. The results had shown that the mitogenome was extremely complex and had a unique multi-branched conformation with total size of 912,473 bp. Comprehensive analysis of protein coding genes of 7 related species showed that there were 40 common genes in their mitogenome. Interestingly, positive selection had become an important factor in the evolution of ccmB, ccmFC, rps10, rps11 and rps7 genes. Furthermore, our data highlighted the repeated trend of homologous fragment migrations between chloroplast and mitochondrial organelles, and 38 homologous fragments were identified. Phylogenetic analysis identified a tree that was basically consistent with the phylogeny of Laurales species described in the APG IV system. To sum up, this study will be helpful to the study of population genetics and evolution of Lindera species.

Single-cell RNA sequencing facilitates the elucidation of the complete biosynthesis of the antidepressant hyperforin in St. John's wort

Hyperforin is the compound responsible for the effectiveness of St. John's wort (Hypericum perforatum) as an antidepressant, but its complete biosynthetic pathway remains unknown. Gene discovery based on co-expression analysis of bulk RNA-sequencing data or genome mining failed to discover the missing steps in hyperforin biosynthesis. In this study, we sequenced the 1.54-Gb tetraploid H. perforatum genome assembled into 32 chromosomes with the scaffold N50 value of 42.44 Mb. By single-cell RNA sequencing, we identified a type of cell, "Hyper cells", wherein hyperforin biosynthesis de novo takes place in both the leaves and flowers. Through pathway reconstitution in yeast and tobacco, we identified and characterized four transmembrane prenyltransferases (HpPT1-4) that are localized at the plastid envelope and complete the hyperforin biosynthetic pathway. The hyperforin polycyclic scaffold is created by a reaction cascade involving an irregular isoprenoid coupling and a tandem cyclization. Our findings reveal how and where hyperforin is biosynthesized, enabling synthetic-biology reconstitution of the complete pathway. Thus, this study not only deepens our comprehension of specialized metabolism at the cellular level but also provides strategic guidance for elucidation of the biosynthetic pathways of other specializied metabolites in plants.

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.

The Spartina alterniflora genome sequence provides insights into the salt-tolerance mechanisms of exo-recretohalophytes

Spartina alterniflora is an exo-recretohalophyte Poaceae species that is able to grow well in seashore, but the genomic basis underlying its adaptation to salt tolerance remains unknown. Here, we report a high-quality, chromosome-level genome assembly of S. alterniflora constructed through PacBio HiFi sequencing, combined with high-throughput chromosome conformation capture (Hi-C) technology and Illumina-based transcriptomic analyses. The final 1.58 Gb genome assembly has a contig N50 size of 46.74 Mb. Phylogenetic analysis suggests that S. alterniflora diverged from Zoysia japonica approximately 21.72 million years ago (MYA). Moreover, whole-genome duplication (WGD) events in S. alterniflora appear to have expanded gene families and transcription factors relevant to salt tolerance and adaptation to saline environments. Comparative genomics analyses identified numerous species-specific genes, significantly expanded genes and positively selected genes that are enriched for 'ion transport' and 'response to salt stress'. RNA-seq analysis identified several ion transporter genes including the high-affinity K+ transporters (HKTs), SaHKT1;2, SaHKT1;3 and SaHKT1;8, and high copy number of Salt Overly Sensitive (SOS) up-regulated under high salt conditions, and the overexpression of SaHKT2;4 in Arabidopsis thaliana conferred salt tolerance to the plant, suggesting specialized roles for S. alterniflora to adapt to saline environments. Integrated metabolomics and transcriptomics analyses revealed that salt stress activate glutathione metabolism, with differential expressions of several genes such as γ-ECS, GSH-S, GPX, GST and PCS in the glutathione metabolism. This study suggests several adaptive mechanisms that could contribute our understanding of evolutional basis of the halophyte.

A novel genus of Pectobacterium bacteriophages display broad host range by targeting several species of Danish soft rot isolates

The bacterial diseases black leg and soft rot in potatoes cause heavy losses of potatoes worldwide. Bacteria within the genus Pectobacteriaceae are the causative agents of black leg and soft rot. The use of antibiotics in agriculture is heavily regulated and no other effective treatment currently exists, but bacteriophages (phages) have shown promise as potential biocontrol agents. In this study we isolated soft rot bacteria from potato tubers and plant tissue displaying soft rot or black leg symptoms collected in Danish fields. We then used the isolated bacterial strains as hosts for phage isolation. Using organic waste, we isolated phages targeting different species within Pectobacterium. Here we focus on seven of these phages representing a new genus primarily targeting P. brasiliense; phage Ymer, Amona, Sabo, Abuela, Koroua, Taid and Pappous. TEM image of phage Ymer showed siphovirus morphotype, and the proposed Ymer genus belongs to the class Caudoviricetes, with double-stranded DNA genomes varying from 39 kb to 43 kb. In silico host range prediction using a CRISPR-Cas spacer database suggested both P. brasiliense, P. polaris and P. versatile as natural hosts for phages within the proposed Ymer genus. A following host range experiment, using 47 bacterial isolates from Danish tubers and plants symptomatic with soft rot or black leg disease verified the in silico host range prediction, as the genus as a group were able to infect all three Pectobacterium species. Phages did, however, primarily target P. brasiliense isolates and displayed differences in host range even within the species level. Two of the phages were able to infect two or more Pectobacterium species. Despite no nucleotide similarity with any phages in the NCBI database, the proposed Ymer genus did share some similarity at the protein level, as well as gene synteny, with currently known phages. None of the phages encoded integrases or other genes typically associated with lysogeny. Similarly, no virulence factors nor antimicrobial resistance genes were found, and combined with their ability to infect several soft rot-causing Pectobacterium species from Danish fields, demonstrates their potential as biocontrol agents against soft rot and black leg diseases in potatoes.

Characterization of gut symbionts from wild-caught Drosophila and other Diptera: description of Utexia brackfieldae gen. nov., sp. nov., Orbus sturtevantii sp. nov., Orbus wheelerorum sp. nov, and Orbus mooreae sp. nov

Non-culture based surveys show that the bacterial family Orbaceae is widespread in guts of insects, including wild Drosophila. Relatively few isolates have been described, and none has been described from Drosophila. We present the isolation and characterization of five strains of Orbaceae from wild-caught flies of the genera Drosophila (Diptera: Drosophilidae) and Neogriphoneura (Diptera: Lauxaniidae). Cells are generally rod-shaped, mesophilic, and measure 0.8-2.0 µm long by 0.3-0.5 µm wide. Optimal growth was observed under ambient atmosphere. Reconstruction of phylogenies from the 16S rRNA gene and from single-copy orthologs verify placement of these strains within Orbaceae. Cells exhibited similar fatty acid profiles to those of other Orbaceae. Strain lpD01T shared 74% average nucleotide identity (ANI) with its closest relatives Ca. Schmidhempelia bombi Bimp and Zophobihabitans entericus IPMB12T. Results from multiple genome-wide similarity comparisons indicate lpD01T should be classified as a novel species within a novel genus. The major respiratory quinone for lpD01T is ubiquinone Q-8. lpD02T, lpD03, lpD04T, and BiBT are more closely related to Orbus hercynius CN3T (76, 77, 76, and 77% ANI, respectively) than to other described Orbaceae. Genomic and phylogenetic analyses suggest that lpD03 and lpD04T belong to the same species and that lpD02T, lpD03/lpD04T, and BiBT are each novel species of the genus Orbus. The proposed names of these strains are Utexia brackfieldae gen. nov., sp. nov. (type strain lpD01T =NCIMB 15517T =ATCC TSD-399T), Orbus sturtevantii sp. nov (type strain lpD02T =NCIMB 15518T =ATCC TSD-400T), Orbus wheelerorum sp. nov. (type strain lpD04T =NCIMB 15520T =ATCC TSD-401T), and Orbus mooreae sp. nov (type strain BiBT=NCIMB 15516T =ATCC TSD-402T). The isolation and characterization of these strains expands the repertoire of culturable bacteria naturally associated with insects, including the model organism D. melanogaster.

Molecular diversity of green-colored microbial mats from hot springs of northern Japan

We acquired and analyzed metagenome and 16S/18S rRNA gene amplicon data of green-colored microbial mats from two hot springs within the Onikobe geothermal region (Miyagi Prefecture, Japan). The two collection sites-Tamago and Warabi-were in proximity and had the same temperature (40 °C), but the Tamago site was connected to a nearby stream, whereas the Warabi site was isolated. Both the amplicon and metagenome data suggest the bacterial, especially cyanobacterial, dominance of the mats; other abundant groups include Chloroflexota, Pseudomonadota, Bacteroidota/Chlorobiota, and Deinococcota. At finer resolution, however, the taxonomic composition entirely differed between the mats. A total of 5 and 21 abundant bacterial 16S rRNA gene OTUs were identified for Tamago and Warabi, respectively; of these, 12 are putative chlorophyll- or rhodopsin-based phototrophs. The presence of phylogenetically diverse microbial eukaryotes was noted, with ciliates and amoebozoans being the most abundant eukaryote groups for Tamago and Warabi, respectively. Fifteen metagenome-assembled genomes (MAGs) were obtained, represented by 13 bacteria, one ciliate (mitochondrion), and one giant virus. A total of 15 novel taxa, including a new deeply branching Chlorobiota species, is noted from the amplicon and MAG data, highlighting the importance of environmental sequencing in uncovering hidden microorganisms.

Draft hybrid genome sequence of phytopathogen Marasmius tenuissimus strain MS-2, isolated from cacao in Ghana

Here, we report a hybrid genome of Marasmius tenuissimus strain MS-2, a cacao thread blight disease causing isolate that was collected from cacao leaves in Tafo, Eastern region, Ghana. The final assembly consists of 2,083 contigs spanning 69,843,039 bp, with 49.21% GC content, 92.6% BUSCO completeness, and scaffold N50 186,871.

A chromosome-level genome assembly of the gall maker pest inquiline, Diomorus aiolomorphi Kamijo (Hymenoptera: Torymidae)

Diomorus aiolomorphi Kamijo (Hymenoptera: Torymidae) is an inquiline of gall maker Aiolomorphus rhopaloides Walker (Hymenoptera: Eurytomidae). They are of significant economic significance and predominantly inhabit bamboo forest. So far, only four scaffold-level genomes have been published for the family Torymidae. In this study, we present a high-quality genome assembly of D. aiolomorphi at the chromosome level, achieved through the integration of Nanopore (ONT) long-read, Illumina pair-end DNA short-read, and High-through Chromosome Conformation Capture (Hi-C) sequencing methods. The final assembly was 1,084.56 Mb in genome size, with 1,083.41 Mb (99.89%) assigned to five pseudochromosomes. The scaffold N50 length reached 224.87 Mb, and the complete Benchmarking Universal Single-Copy Orthologs (BUSCO) score was 97.3%. The genome contained 762.12 Mb of repetitive elements, accounting for 70.27% of the total genome size. A total of 18,011 protein-coding genes were predicted, with 17,829 genes being functionally annotated. The high-quality genome assembly of D. aiolomorphi presented in this study will serve as a valuable genomic resource for future research on parasitoid wasps. The results of this study may also contribute to the development of biological control strategies for pest management in bamboo forests, enhancing ecological balance and economic sustainability.

A gap-free genome assembly of Fusarium oxysporum f. sp. conglutinans, a vascular wilt pathogen

Fusarium oxysporum is an asexual filamentous fungus that causes vascular wilt in hundreds of crop plants and poses a threat to public health through Fusariosis. F. oxysporum f. sp. conglutinans strain Fo5176, originally isolated from Brassica oleracea, is pathogenic to Arabidopsis, making it a model pathosystem for dissecting the molecular mechanisms underlying host-pathogen interactions. Assembling the F. oxysporum genome is notoriously challenging due to the presence of repeat-rich accessory chromosomes. Here, we report a gap-free genome assembly of Fo5176 using PacBio HiFi and Hi-C data. The 69.56 Mb assembly contained 18 complete chromosomes, including all centromeres and most telomeres (20/36), representing the first gap-free genome sequence of a pathogenic F. oxysporum strain. In total, 21,460 protein-coding genes were annotated, a 26.3% increase compared to the most recent assembly. This high-quality reference genome for F. oxysporum f. sp. conglutinans Fo5176 provides a valuable resource for further research into fungal pathobiology and evolution.

Ubiquitination and degradation of plant helper NLR by the Ralstonia solanacearum effector RipV2 overcome tomato bacterial wilt resistance

The Ralstonia solanacearum species complex causes bacterial wilt in a variety of crops. Tomato cultivar Hawaii 7996 is a widely used resistance resource; however, the resistance is evaded by virulent strains, with the underlying mechanisms still unknown. Here, we report that the phylotype Ⅱ strain ES5-1 can overcome Hawaii 7996 resistance. RipV2, a type Ⅲ effector specific to phylotype Ⅱ strains, is vital in overcoming tomato resistance. RipV2, which encodes an E3 ubiquitin ligase, suppresses immune responses and Toll/interleukin-1 receptor/resistance nucleotide-binding/leucine-rich repeat (NLR) (TNL)-mediated cell death. Tomato helper NLR N requirement gene 1 (NRG1), enhanced disease susceptibility 1 (EDS1), and senescence-associated gene 101b (SAG101b) are identified as RipV2 target proteins. RipV2 is essential for ES5-1 virulence in Hawaii 7996 but not in SlNRG1-silenced tomato, demonstrating SlNRG1 to be an RipV2 virulence target. Our results dissect the mechanisms of RipV2 in disrupting immunity and highlight the importance of converged immune components in conferring bacterial wilt resistance.

Association with TFIIIC limits MYCN localisation in hubs of active promoters and chromatin accumulation of non-phosphorylated RNA polymerase II

MYC family oncoproteins regulate the expression of a large number of genes and broadly stimulate elongation by RNA polymerase II (RNAPII). While the factors that control the chromatin association of MYC proteins are well understood, much less is known about how interacting proteins mediate MYC's effects on transcription. Here, we show that TFIIIC, an architectural protein complex that controls the three-dimensional chromatin organisation at its target sites, binds directly to the amino-terminal transcriptional regulatory domain of MYCN. Surprisingly, TFIIIC has no discernible role in MYCN-dependent gene expression and transcription elongation. Instead, MYCN and TFIIIC preferentially bind to promoters with paused RNAPII and globally limit the accumulation of non-phosphorylated RNAPII at promoters. Consistent with its ubiquitous role in transcription, MYCN broadly participates in hubs of active promoters. Depletion of TFIIIC further increases MYCN localisation to these hubs. This increase correlates with a failure of the nuclear exosome and BRCA1, both of which are involved in nascent RNA degradation, to localise to active promoters. Our data suggest that MYCN and TFIIIC exert an censoring function in early transcription that limits promoter accumulation of inactive RNAPII and facilitates promoter-proximal degradation of nascent RNA.

Hybrid assembly and comparative genomics unveil insights into the evolution and biology of the red-legged partridge

The red-legged partridge Alectoris rufa plays a crucial role in the ecosystem of southwestern Europe, and understanding its genetics is vital for conservation and management. Here we sequence, assemble, and annotate a highly contiguous and nearly complete version of its genome. This assembly encompasses 96.9% of the avian genes flagged as essential in the BUSCO aves_odb10 dataset. Moreover, we pinpointed RNA and protein-coding genes, 95% of which had functional annotations. Notably, we observed significant chromosome rearrangements in comparison to quail (Coturnix japonica) and chicken (Gallus gallus). In addition, a comparative phylogenetic analysis of these genomes suggests that A. rufa and C. japonica diverged roughly 20 million years ago and that their common ancestor diverged from G. gallus 35 million years ago. Our assembly represents a significant advancement towards a complete reference genome for A. rufa, facilitating comparative avian genomics, and providing a valuable resource for future research and conservation efforts for the red-legged partridge.

A high-quality chromosome-level genome assembly of the topmouth culter (Culter alburnus Basilewsky, 1855)

Culter alburnus is extensively distributed in various rivers and lakes across China. As a widely adaptive fish species, it has significant economic values and special ecological roles. To meet research demands and provide better genomic resources, in this research, a chromosome-level genome assembly was constructed using HiFi long-reads and Hi-C sequencing data. Compared with the published versions, our genome assembly is of higher quality with only 31 gaps and closer to its true structure and sequence. The genome size was 1.052 Gb, with a contig N50 of 32.92 Mb and a scaffold N50 of 43.09 Mb. 55 contigs were anchored to 24 chromosomes on the basis of Hi-C data. A total of 598.23 Mb of repetitive sequences were annotated and 28,228 protein-coding genes were predicted. Additionally, BUSCO assessment indicated assembly and annotation scores of 98.3% and 99.2%, respectively. This high-quality genome will provide scientific support for excavating the species characteristics of C. alburnus and exploring its molecular mechanisms in response to environmental changes and stress.

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.

Comprehensive analysis of the Lycopodium japonicum mitogenome reveals abundant tRNA genes and cis-spliced introns in Lycopodiaceae species

Lycophytes and ferns represent one of the earliest-diverging lineages of vascular plants, with the Lycopodiaceae family constituting the basal clade among lycophytes. In this research, we successfully assembled and annotated the complete Lycopodium japonicum Thunb. (L. japonicum) mitochondrial genome (mitogenome) utilizing PacBio HiFi sequencing data, resulting in a single circular molecule with a size of 454,458 bp. 64 unique genes were annotated altogether, including 34 protein-coding genes, 27 tRNAs and 3 rRNAs. It also contains 32 group II introns, all of which undergo cis-splicing. We identified 195 simple sequence repeats, 1,948 dispersed repeats, and 92 tandem repeats in the L. japonicum mitogenome. Collinear analysis indicated that the mitogenomes of Lycopodiaceae are remarkably conserved compared to those of other vascular plants. We totally identified 326 RNA editing sites in 31 unique protein-coding genes with 299 sites converting cytosine to uracil and 27 sites the reverse. Notably, the L. japonicum mitogenome has small amounts foreign DNA from plastid or nuclear origin, accounting for only 2.81% of the mitogenome. The maximum likelihood phylogenetic analysis based on 23 diverse land plant mitogenomes and plastid genomes supports the basal position of lycophytes within vascular plants and they form a sister clade to all other vascular lineages, which is consistent with the PPG I classification system. As the first reported mitogenome of Lycopodioideae subfamily, this study enriches our understanding of Lycopodium mitogenomes, and sets the stage for future research on mitochondrial diversity and evolution within the lycophytes and ferns.

Comparative analyses of the banded alder borer (Rosalia funebris) and Asian longhorned beetle (Anoplophora glabripennis) genomes reveal significant differences in genome architecture and gene content among these and other Cerambycidae

Rosalia funebris (RFUNE; Cerambycidae), the banded alder borer, is a longhorn beetle whose larvae feed on the wood of various economically and ecologically significant trees in western North America. Adults are short-lived and not known to consume plant material substantially. We sequenced, assembled, and annotated the RFUNE genome using HiFi and RNASeq data. We documented genome architecture and gene content, focusing on genes putatively involved in plant feeding (phytophagy). Comparisons were made to the well-studied genome of the Asian longhorned beetle (AGLAB; Anoplophora glabripennis) and other Cerambycidae. The 814 Mb RFUNE genome assembly was distributed across 42 contigs, with an N50 of 30.18 Mb. Repetitive sequences comprised 60.27% of the genome, and 99.0% of expected single-copy orthologous genes were fully assembled. We identified 12,657 genes, fewer than in the four other species studied, and 46.4% fewer than for Aromia moschata (same subfamily as RFUNE). Of the 7,258 orthogroups shared between RFUNE and AGLAB, 1,461 had more copies in AGLAB and 1,023 had more copies in RFUNE. We identified 240 genes in RFUNE that putatively arose via horizontal transfer events. The RFUNE genome encoded substantially fewer putative plant cell wall degrading enzymes than AGLAB, which may relate to the longer-lived plant-feeding adults of the latter species. The RFUNE genome provides new insights into cerambycid genome architecture and gene content and provides a new vantage point from which to study the evolution and genomic basis of phytophagy in beetles.

Characterization of the novel phage vB_BceP_LY3 and its potential role in controlling Bacillus cereus in milk and rice

Bacillus cereus is a foodborne pathogen that induces vomiting and diarrhea in affected individuals. It exhibits resistance to traditional sterilization methods and has a high contamination rate in dairy products and rice. Therefore, the development of a new food safety controlling strategy is necessary. In this research, we isolated and identified a novel phage named vB_BceP_LY3, which belongs to a new genus of the subfamily Northropvirinae. This phage demonstrates a short latency period and remains stable over a wide range of temperatures (4-60 °C) and pH levels (4-11). The 28,124 bp genome of LY3 does not contain any antibiotic-resistance genes or virulence factors. With regards to its antibacterial properties, LY3 not only effectively inhibits the growth of B. cereus in TSB (tryptic soy broth), but also demonstrates significant inhibitory effects in various food matrices. Specifically, LY3 treatment at 4 °C with a high MOI (MOI = 10,000) can maintain B. cereus levels below the detection limit for up to 24 h in milk. LY3 represents a safe and promising biocontrol agent against B. cereus, possessing long-term antibacterial capabilities and stability.

Climate-driven succession in marine microbiome biodiversity and biogeochemical function

Seasonal and El Niño-Southern Oscillation (ENSO) warming result in similar ocean changes as predicted with climate change. Climate-driven environmental cycles have strong impacts on microbiome diversity, but impacts on microbiome function are poorly understood. We quantified changes in microbial genomic diversity and functioning over 11 years covering seasonal and ENSO cycles at a coastal site in the southern California Current. We observed seasonal oscillations between large genome lineages during cold, nutrient rich conditions in winter and spring versus small genome lineages, including Prochlorococcus and Pelagibacter , in summer and fall. Parallel interannual changes separated communities depending on ENSO condition. Biodiversity shifts translated into clear oscillations in microbiome functional potential. Ocean warming induced an ecosystem with less iron but more macronutrient stress genes, depressed organic carbon degradation potential and biomass, and elevated carbon-to-nutrient biomass ratios. The consistent microbial response observed across time-scales points towards large climate-driven changes in marine ecosystems and biogeochemical cycles.

Comparative genomics and stable isotope analysis reveal the saprotrophic-pathogenic lifestyle of a neotropical fungus

In terrestrial forested ecosystems, fungi may interact with trees in at least three distinct ways: (i) associated with roots as symbionts; (ii) as pathogens in roots, trunks, leaves, flowers, and fruits; or (iii) decomposing dead tree tissues on soil or even on dead tissues in living trees. Distinguishing the latter two nutrition modes is rather difficult in Hymenochaetaceae (Basidiomycota) species. Herein, we have used an integrative approach of comparative genomics, stable isotopes, host tree association, and bioclimatic data to investigate the lifestyle ecology of the scarcely known neotropical genus Phellinotus, focusing on the unique species Phellinotus piptadeniae. This species is strongly associated with living Piptadenia gonoacantha (Fabaceae) trees in the Atlantic Forest domain on a relatively high precipitation gradient. Phylogenomics resolved P. piptadeniae in a clade that also includes both plant pathogens and typical wood saprotrophs. Furthermore, both genome-predicted Carbohydrate-Active Enzymes (CAZy) and stable isotopes (δ13C and δ15N) revealed a rather flexible lifestyle for the species. Altogether, our findings suggest that P. piptadeniae has been undergoing a pathotrophic specialization in a particular tree species while maintaining all the metabolic repertoire of a wood saprothroph.

IMPORTANCE

This is the first genomic description for Phellinotus piptadeniae. This basidiomycete is found across a broad range of climates and ecosystems in South America, including regions threatened by extensive agriculture. This fungus is also relevant considering its pathotrophic-saprotrophic association with Piptadenia goanocantha, which we began to understand with these new results that locate this species among biotrophic and necrotrophic fungi.

Chromosome-level genome assembly and annotation of Flueggea virosa (Phyllanthaceae)

Flueggea virosa (Roxb. ex Willd.) Royle, an evergreen shrub and small tree in the Phyllanthaceae family, holds significant potential in garden landscaping and pharmacological applications. However, the lack of genomic data has hindered further scientific understanding of its horticultural and medicinal values. In this study, we have assembled a haplotype-resolved genome of F. virosa for the first time. The two haploid genomes, named haplotype A genome and haplotype B genome, are 487.33 Mb and 477.53 Mb in size, respectively, with contig N50 lengths of 31.45 Mb and 32.81 Mb. More than 99% of the assembled sequences were anchored to 13 pairs of pseudo-chromosomes. Furthermore, 21,587 and 21,533 protein-coding genes were predicted in haplotype A and haplotype B genomes, respectively. The availability of this chromosome-level genome fills the gap in genomic data for F. virosa and provides valuable resources for molecular studies of this species, supporting future research on speciation, functional genomics, and comparative genomics within the Phyllanthaceae family.

Complete genome sequence of Escherichia coli O157:H7 phage Φ241

We report the genome sequence of phage Φ241 infecting Escherichia coli O157:H7. Phage Φ241 was isolated from an industrial cucumber fermentation at high acidity (pH 3.7) and high salinity (5% NaCl). The phage genome consists of a 157,291 bp circular double-stranded DNA with 203 coding regions and 44.96% GC content.

Chromosome-level and haplotype-resolved genome assembly of Dracaena cambodiana (Asparagaceae)

Dracaena cambodiana Pierre ex Gagn. (Asparagaceae) is the source plant of Dragon's blood and has high ornamental values in gardening. Currently, this species is classified as the second-class state-protected species in the National Key Protected Wild Plants (NKPWP) of China. However, limited genomic data has hindered a more comprehensive scientific understanding of the processes involved in the production of Dragon's blood and the related conservation genomics research. In this study, we assembled a haplotype-resolved genome of D. cambodiana. The haploid genomes, haplotype A and haplotype B, are 1,015.22 Mb and 1,003.13 Mb in size, respectively. The completeness of haplotype A and haplotype B genomes was 98.60% and 98.20%, respectively, using the "embryophyta_10" dataset. Haplotype A and haplotype B genomes contained 27,361 and 27,066 protein-coding genes, respectively, with nearly all being functionally annotated. These findings provide new insights into the genomic characteristics of D. cambodiana and will offer additional genomic resources for studying the biosynthesis mechanism of Dragon's blood and the horticultural application of Dragon trees.

U, R. K, M.W. R, J. P, F. A, J. H. J. L. Draft genome sequence of Collimonas sp. strain H4R21, an effective mineral-weathering bacterial strain isolated from the beech rhizosphere

We present the draft genome sequence of Collimonas sp. strain H4R21, isolated from the rhizosphere of Fagus sylvatica in the forest experimental site of Montiers (France). This genome features coding capacity for plant growth promotion, such as the ability to solubilize minerals, to produce siderophores and antifungal secondary metabolites.

Chromosome-level genome assembly and annotation of the Rhabdophis nuchalis (Hubei keelback)

Rhabdophis nuchalis, a snake widely distributed in China, possesses a unique trait: glands beneath the skin on its neck and back, known as nucho-dorsal glands. These features make it a valuable subject for studying genetic diversity and the evolution of complex traits. In this study, we obtained a high-quality chromosome-level reference genome of R. nuchalis using MGI short-read sequencing, PacBio Revio long-read sequencing, and Hi-C sequencing techniques. The final assembly comprised 1.92 Gb of the R. nuchalis genome, anchored to 20 chromosomes (including 9 macrochromosomes and 11 microchromosomes), with a contig N50 of 104.79 Mb, a scaffold N50 of 204.96 Mb, and a BUSCO completeness of 97.50%. Additionally, we annotated a total of 1.09 Gb of repetitive sequences (which constitute 56.51% of the entire genome) and identified 22,057 protein-coding genes. This high-quality reference genome of R. nuchalis furnishes essential genomic data for comprehending the genetic diversity and evolutionary history of the species, as well as for facilitating species conservation efforts and comparative genomics studies.

Tyrosine transfer RNA levels and modifications during blood-feeding and vitellogenesis in the mosquito, Aedes aegypti

Mosquitoes such as Aedes aegypti must consume a blood meal for the nutrients necessary for egg production. Several transcriptome and proteome changes occur post-blood meal that likely corresponds with codon usage alterations. Transfer RNA (tRNA) is the adapter molecule that reads messenger RNA codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codon decoding, improving the accuracy and efficiency of protein synthesis. Here, we examined tRNA modifications and transcripts associated with the blood meal and subsequent periods of vitellogenesis in A. aegypti. More specifically, we assessed tRNA transcript abundance and modification levels in the fat body at critical times post blood-feeding. Based on a combination of alternative codon usage and identification of particular modifications, we discovered that increased transcription of tyrosine tRNAs is likely critical during the synthesis of egg yolk proteins in the fat body following a blood meal. Altogether, changes in both the abundance and modification of tRNA are essential factors in the process of vitellogenin production after blood-feeding in mosquitoes.

The first mitogenome of the subfamily Stenoponiinae (Siphonaptera: Ctenophthalmidae) and implications for its phylogenetic position

Fleas are the most important insect vectors that parasitize warm-blooded animals and are known vectors of zoonotic pathogens. A recent study showed that Stenoponia polyspina parasitizing Eospalax baileyi in Zoige County have carried Bartonella spp. and Spotted fever group Rickettsia (SFGR). Accurate identification and differentiation of fleas are essential for prevention and control of zoonotic pathogens. To understand phylogenetic relationship of the subfamily Stenoponiinae, we described morphological characteristics of S. polyspina and sequenced its mitogenome with 14,933 bp in size and high A + T content (~ 79%). The S. polyspina mitogenome retained the ancestral pattern of mitochondrial gene arrangement of arthropods without rearrangement. The start codons of 13 protein-coding genes (PCGs) are traditional ATN and the stop codons are TAA or TAG. Anticodon loops of all tRNA genes were 7 bp except for trnL2 and trnD had anticodon loops with 9 bp and the abnormal anticodon loops may be associated with frameshifting mutation. Genetic distance and Ka/Ks ratios indicated that all 13 PCGs of S. polyspina were subjected to purifying selection, with cox1 at the slowest rate and atp8 at the fastest rate. The mitogenomes of 24 species representing 7 families in the order Siphonaptera were selected to reconstruct phylogenetic tree based on concatenated nucleotide sequences of two datasets (PCGRNA matrix and PCG12RNA matrix) using Bayesian inference (BI) and Maximum likelihood (ML) methods. Phylogenetic tree supported that the superfamilies Ceratophylloidea, Vermipsylloidea, Pulicoidea were monophyletic and the superfamily Hystrichopsylloidea was paraphyletic. The family Ctenophthalmidae was monophyletic in PCGRNA-ML (codon partition) and paraphyletic in the remain trees. S. polyspina belongs to the subfamily Stenoponiinae was closely more related to the subfamily Rhadinopsyllinae. This paper explored phylogenetic position of diverse clades within the order Siphonaptera based on morphological and mitogenome data of S. polyspina. Our research enriched NCBI database of the order Siphonaptera.

Impact of the chemical modification of tRNAs anticodon loop on the variability and evolution of codon usage in proteobacteria

Despite the highly conserved nature of the genetic code, the frequency of usage of each codon can vary significantly. The evolution of codon usage is shaped by two main evolutionary forces: mutational bias and selection pressures. These pressures can be driven by environmental factors, but also by the need for efficient translation, which depends heavily on the concentration of transfer RNAs (tRNAs) within the cell. The data presented here supports the proposal that tRNA modifications play a key role in shaping the overall preference of codon usage in proteobacteria. Interestingly, some codons, such as CGA and AGG (encoding arginine), exhibit a surprisingly low level of variation in their frequency of usage, even across genomes with differing GC content. These findings suggest that the evolution of GC content in proteobacterial genomes might be primarily driven by changes in the usage of a specific subset of codons, whose usage is itself influenced by tRNA modifications.

Chromosome-level genome assembly and annotation of the black sea urchin Arbacia lixula (Linnaeus, 1758)

The black sea urchin (Arbacia lixula) is a keystone species inhabiting the coastal shallow waters of the Mediterranean Sea, which is a key driver of littoral communities' structure. Here, we present the first genome assembly and annotation of this species, standing as the first Arbacioida genome, including both nuclear and mitochondrial genomes. To obtain a chromosome-level assembly, we used a combination of PacBio high fidelity (HiFi) reads and chromatin capture reads (Omni-C). In addition, we generated a high-quality nuclear annotation of both coding and non-coding genes, by using published RNA-Seq data from several individuals of A. lixula and gene models from closely related species. The nuclear genome assembly has a total span of 607.91 Mb, being consistent with its experimentally estimated genome size. The assembly contains 22 chromosome-scale scaffolds (96.52% of the total length), which coincides with its known karyotype. A total of 72,767 transcripts were predicted from the nuclear genome, 24,171 coding, and 48,596 non-coding that included lncRNA, snoRNA, and tRNAs. The circularized mitochondrial genome had 15,740 bp comprising 13 protein-coding genes, 2 rRNA, and 22 tRNA. This reference genome will enhance ongoing A. lixula studies and benefit the wider sea urchin scientific community.

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.