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.

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.

Characterization of two virulent Salmonella phages and transient application in egg, meat and lettuce safety

Salmonella, a prominent foodborne pathogen, has posed enduring challenges to the advancement of food safety and global public health. The escalating concern over antibiotic misuse, resulting in the excessive presence of drug residues in animal-derived food products, necessitates urgent exploration of alternative strategies for Salmonella control. Bacteriophages emerge as promising green biocontrol agents against pathogenic bacteria. This study delineates the identification of two novel virulent Salmonella phages, namely phage vB_SalS_ABTNLsp11241 (referred to as sp11241) and phage 8-19 (referred to as 8-19). Both phages exhibited efficient infectivity against Salmonella enterica serotype Enteritidis (SE). Furthermore, this study evaluated the effectiveness of two phages to control SE in three different foods (whole chicken eggs, raw chicken meat, and lettuce) at different MOIs (1, 100, and 10000) at 4°C. It's worth noting that sp11241 and 8-19 achieved complete elimination of SE on eggs after 3 h and 6 h at MOI = 100, and after 2 h and 5 h at MOI = 10000, respectively. After 12 h of treatment with sp11241, a maximum reduction of 3.17 log10 CFU/mL in SE was achieved on raw chicken meat, and a maximum reduction of 3.00 log10 CFU/mL was achieved on lettuce. Phage 8-19 has the same effect on lettuce as sp11241, but is slightly less effective than sp11241 on chicken meat (a maximum 2.69 log10 CFU/mL reduction). In conclusion, sp11241 and 8-19 exhibit considerable potential for controlling Salmonella contamination in food at a low temperature and represent viable candidates as green antibacterial agents for food applications.

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.

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 of Plagiognathops microlepis based on PacBio HiFi and Hi-C sequencing

Plagiognathops microlepis is an economic freshwater fish in the subfamily Xenocyprinae of Cyprinidae. It is widely distributed in the freshwater ecosystem of China, with moderate economic value and broad development prospects. However, the lack of genomic resources has limited our understanding on the genetic basis, phylogenetic status and adaptive evolution strategies of this fish. Here, we assembled a chromosome-level reference genome of P. microlepis by integrating Pacbio HiFi long-reads, Illumina short-reads and Hi-C sequencing data. The size of this genome is 1004.34 Mb with a contig N50 of 38.80 Mb. Using Hi-C sequencing data, 99.59% of the assembled sequences were further anchored to 24 chromosomes. A total of 578.91 Mb repeat sequences and 28,337 protein-coding genes were predicted in the current genome, of which, 26,929 genes were functionally annotated. This genome provides valuable information for investigating the phylogeny and evolutionary history of cyprinid fishes, as well as the genetic basis of adaptive strategies and special traits in P. microlepis.

Characterization and efficacy of Salmonella phage cocktail PHA46 in the control of Salmonella Newport and Typhimurium internalized into cherry tomatoes

Non-typhoid Salmonella enterica causes salmonellosis illness, and this bacterium can contaminate food throughout the production chain, including those that are consumed as raw products. Salmonella enterica can adhere to and internalize into fresh produce such as cherry tomatoes. It has been reported that lytic bacteriophages (phages) can be used as a biocontrol agent in the agricultural field, being an alternative for the control of Salmonella in red meat, fish, lettuce, and cabbage. The aim of this study was to characterize the two phages present in the PHA46 cocktail to determine their morphology, genome, host range, and resistance to different temperatures and pHs values; and later evaluate their lytic activity to reduce the adherence to and internalization of Salmonella enterica serovars Newport and Typhimurium into cherry tomatoes. In addition, in this work, we also explored the effect of the PHA46 cocktail on the virulence of S. Newport-45 and S. Typhimurium SL1344, recovered from the interior of cherry tomatoes, on the lifespan of the animal model Caenorhabditis elegans. The nematode C. elegans, recently has been used to test the virulence of Salmonella and it is easy to maintain and work with in the laboratory. The results revealed that the morphology obtained by Transmission Electron Microscopy of two phages from the PHA46 cocktail correspond to a myovirus, the analyses of their genomes sequences did not report virulence or antimicrobial resistance genes. The PHA46 sample is specific for 33 different serovars from different Salmonella strains and shows stability at 7 °C and pH 6. Also, the PHA46 cocktail was effective in reducing the adherence of S. Newport-45 and S. Typhimurium SL1344 to cherry tomatoes, at an average of 0.9 log10, respectively. Regarding internalized bacteria, the reduction was at an average of 1.2 log10, of the serovars mentioned above. The lifespan experiments in C. elegans showed by itself, that the PHA46 cocktail was harmless to the nematode, and the virulence from both Salmonella strains grown in vitro is diminished in the presence of the PHA46 cocktail. In conclusion, these results showed that the PHA46 cocktail could be a good candidate to be used as a biocontrol agent against Salmonella enterica.

Human T-cell leukemia virus type 1 uses a specific tRNAPro isodecoder to prime reverse transcription

Human T-cell leukemia virus type 1 (HTLV-1) is the only oncogenic human retrovirus discovered to date. All retroviruses are believed to use a host cell tRNA to prime reverse transcription (RT). In HTLV-1, the primer-binding site (PBS) in the genomic RNA is complementary to the 3' 18 nucleotides (nt) of human tRNAPro The human genome encodes 20 cytoplasmic tRNAPro genes representing seven isodecoders, all of which share the same 3' 18 nt sequence but vary elsewhere. Whether all tRNAPro isodecoders are used to prime RT in cells is unknown. A previous study showed that a 3' 18 nt tRNAPro-derived fragment (tRFPro) is packaged into HTLV-1 particles and can serve as an RT primer in vitro. The role of this tRNA fragment in the viral life cycle is unclear. In retroviruses, N1-methylation of the tRNA primer at position A58 (m1A) is essential for successful plus-strand transfer. Using primer-extension assays performed in chronically HTLV-1-infected cells, we found that A58 of tRNAPro is m1A-modified, implying that full-length tRNAPro is capable of facilitating successful plus-strand transfer. Analysis of HTLV-1 RT primer extension products indicated that full-length tRNAPro is likely to be the primer. To determine which tRNAPro isodecoder is used as the RT primer, we sequenced the minus-strand strong-stop RT product containing the intact tRNA primer and established that HTLV-1 primes RT using a specific tRNAPro UGG isodecoder. Further studies are required to understand how this primer is annealed to the highly structured HTLV-1 PBS and to investigate the role of tRFPro in the viral life cycle.

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.

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.

Complete mitochondrial genome of Melia azedarach L., reveals two conformations generated by the repeat sequence mediated recombination

Melia azedarach is a species of enormous value of pharmaceutical industries. Although the chloroplast genome of M. azedarach has been explored, the information of mitochondrial genome (Mt genome) remains surprisingly limited. In this study, we used a hybrid assembly strategy of BGI short-reads and Nanopore long-reads to assemble the Mt genome of M. azedarach. The Mt genome of M. azedarach is characterized by two circular chromosomes with 350,142 bp and 290,387 bp in length, respectively, which encodes 35 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes. A pair of direct repeats (R1 and R2) were associated with genome recombination, resulting in two conformations based on the Sanger sequencing and Oxford Nanopore sequencing. Comparative analysis identified 19 homologous fragments between Mt and chloroplast genome, with the longest fragment of 12,142 bp. The phylogenetic analysis based on PCGs were consist with the latest classification of the Angiosperm Phylogeny Group. Notably, a total of 356 potential RNA editing sites were predicted based on 35 PCGs, and the editing events lead to the formation of the stop codon in the rps10 gene and the start codons in the nad4L and atp9 genes, which were verified by PCR amplification and Sanger sequencing. Taken together, the exploration of M. azedarach gap-free Mt genome provides a new insight into the evolution research and complex mitogenome architecture.

Telomere-to-telomere Citrullus super-pangenome provides direction for watermelon breeding

To decipher the genetic diversity within the cucurbit genus Citrullus, we generated telomere-to-telomere (T2T) assemblies of 27 distinct genotypes, encompassing all seven Citrullus species. This T2T super-pangenome has expanded the previously published reference genome, T2T-G42, by adding 399.2 Mb and 11,225 genes. Comparative analysis has unveiled gene variants and structural variations (SVs), shedding light on watermelon evolution and domestication processes that enhanced attributes such as bitterness and sugar content while compromising disease resistance. Multidisease-resistant loci from Citrullus amarus and Citrullus mucosospermus were successfully introduced into cultivated Citrullus lanatus. The SVs identified in C. lanatus have not only been inherited from cordophanus but also from C. mucosospermus, suggesting additional ancestors beyond cordophanus in the lineage of cultivated watermelon. Our investigation substantially improves the comprehension of watermelon genome diversity, furnishing comprehensive reference genomes for all Citrullus species. This advancement aids in the exploration and genetic enhancement of watermelon using its wild relatives.

Unraveling genomic features and phylogenomics through the analysis of three Mexican endemic Myotis genomes

Background

Genomic resource development for non-model organisms is rapidly progressing, seeking to uncover molecular mechanisms and evolutionary adaptations enabling thriving in diverse environments. Limited genomic data for bat species hinder insights into their evolutionary processes, particularly within the diverse Myotis genus of the Vespertilionidae family. In Mexico, 15 Myotis species exist, with three-M. vivesi, M. findleyi, and M. planiceps-being endemic and of conservation concern.

Methods

We obtained samples of Myotis vivesi, M. findleyi, and M. planiceps for genomic analysis. Each of three genomic DNA was extracted, sequenced, and assembled. The scaffolding was carried out utilizing the M. yumanensis genome via a genome-referenced approach within the ntJoin program. GapCloser was employed to fill gaps. Repeat elements were characterized, and gene prediction was done via ab initio and homology methods with MAKER pipeline. Functional annotation involved InterproScan, BLASTp, and KEGG. Non-coding RNAs were annotated with INFERNAL, and tRNAscan-SE. Orthologous genes were clustered using Orthofinder, and a phylogenomic tree was reconstructed using IQ-TREE.

Results

We present genome assemblies of these endemic species using Illumina NovaSeq 6000, each exceeding 2.0 Gb, with over 90% representing single-copy genes according to BUSCO analyses. Transposable elements, including LINEs and SINEs, constitute over 30% of each genome. Helitrons, consistent with Vespertilionids, were identified. Values around 20,000 genes from each of the three assemblies were derived from gene annotation and their correlation with specific functions. Comparative analysis of orthologs among eight Myotis species revealed 20,820 groups, with 4,789 being single copy orthogroups. Non-coding RNA elements were annotated. Phylogenomic tree analysis supported evolutionary chiropterans' relationships. These resources contribute significantly to understanding gene evolution, diversification patterns, and aiding conservation efforts for these endangered bat species.

Exploring the regulatory role of tsRNAs in the TNF signaling pathway: Implications for cancer and non-cancer diseases

Transfer RNA-derived small RNAs (tsRNAs), a recently identified subclass of small non-coding RNAs (sncRNAs), emerge through the cleavage of mature transfer RNA (tRNA) or tRNA precursors mediated by specific enzymes. The tumor necrosis factor (TNF) protein, a signaling molecule produced by activated macrophages, plays a pivotal role in systemic inflammation. Its multifaceted functions include the capacity to eliminate or hinder tumor cells, enhance the phagocytic capabilities of neutrophils, confer resistance against infections, induce fever, and prompt the production of acute phase proteins. Notably, four TNF-related tsRNAs have been conclusively linked to distinct diseases. Examples include 5'tiRNA-Gly in skeletal muscle injury, tsRNA-21109 in systemic lupus erythematosus (SLE), tRF-Leu-AAG-001 in endometriosis (EMs), and tsRNA-04002 in intervertebral disk degeneration (IDD). These tsRNAs exhibit the ability to suppress the expression of TNF-α. Additionally, KEGG analysis has identified seven tsRNAs potentially involved in modulating the TNF pathway, exerting their influence across a spectrum of non-cancerous diseases. Noteworthy instances include aberrant tiRNA-Ser-TGA-001 and tRF-Val-AAC-034 in intrauterine growth restriction (IUGR), irregular tRF-Ala-AGC-052 and tRF-Ala-TGC-027 in obesity, and deviant tiRNA-His-GTG-001, tRF-Ser-GCT-113, and tRF-Gln-TTG-035 in irritable bowel syndrome with diarrhea (IBS-D). This comprehensive review explores the biological functions and mechanisms of tsRNAs associated with the TNF signaling pathway in both cancer and other diseases, offering novel insights for future translational medical research.

Phenotypic characterization and genomic analysis of a Salmonella phage L223 for biocontrol of Salmonella spp. in poultry

The escalating incidence of foodborne salmonellosis poses a significant global threat to food safety and public health. As antibiotic resistance in Salmonella continues to rise, there is growing interest in bacteriophages as potential alternatives. In this study, we isolated, characterized, and evaluated the biocontrol efficacy of lytic phage L223 in chicken meat. Phage L223 demonstrated robust stability across a broad range of temperatures (20-70 °C) and pH levels (2-11) and exhibited a restricted host range targeting Salmonella spp., notably Salmonella Typhimurium and Salmonella Enteritidis. Characterization of L223 revealed a short latent period of 30 min and a substantial burst size of 515 PFU/cell. Genomic analysis classified L223 within the Caudoviricetes class, Guernseyvirinae subfamily and Jerseyvirus genus, with a dsDNA genome size of 44,321 bp and 47.9% GC content, featuring 72 coding sequences devoid of antimicrobial resistance, virulence factors, toxins, and tRNA genes. Application of L223 significantly (p < 0.005) reduced Salmonella Typhimurium ATCC 14,028 counts by 1.24, 2.17, and 1.55 log CFU/piece after 2, 4, and 6 h of incubation, respectively, in experimentally contaminated chicken breast samples. These findings highlight the potential of Salmonella phage L223 as a promising biocontrol agent for mitigating Salmonella contamination in food products, emphasizing its relevance for enhancing food safety protocols.

Discovery of a novel symbiotic lineage associated with a hematophagous leech from the genus Haementeria

Similarly to other strict blood feeders, leeches from the Haementeria genus (Hirudinida: Glossiphoniidae) have established a symbiotic association with bacteria harbored intracellularly in esophageal bacteriomes. Previous genome sequence analyses of these endosymbionts revealed co-divergence with their hosts, a strong genome reduction, and a simplified metabolism largely dedicated to the production of B vitamins, which are nutrients lacking from a blood diet. 'Candidatus Providencia siddallii' has been identified as the obligate nutritional endosymbiont of a monophyletic clade of Mexican and South American Haementeria spp. However, the Haementeria genus includes a sister clade of congeners from Central and South America, where the presence or absence of the aforementioned symbiont taxon remains unknown. In this work, we report on a novel bacterial endosymbiont found in a representative from this Haementeria clade. We found that this symbiont lineage has evolved from within the Pluralibacter genus, known mainly from clinical but also environmental strains. Similarly to Ca. Providencia siddallii, the Haementeria-associated Pluralibacter symbiont displays clear signs of genome reduction, accompanied by an A+T-biased sequence composition. Genomic analysis of its metabolic potential revealed a retention of pathways related to B vitamin biosynthesis, supporting its role as a nutritional endosymbiont. Finally, comparative genomics of both Haementeria symbiont lineages suggests that an ancient Providencia symbiont was likely replaced by the novel Pluralibacter one, thus constituting the first reported case of nutritional symbiont replacement in a leech without morphological changes in the bacteriome.

IMPORTANCE

Obligate symbiotic associations with a nutritional base have likely evolved more than once in strict blood-feeding leeches. Unlike those symbioses found in hematophagous arthropods, the nature, identity, and evolutionary history of these remains poorly studied. In this work, we further explored obligate nutritional associations between Haementeria leeches and their microbial symbionts, which led to the unexpected discovery of a novel symbiosis with a member of the Pluralibacter genus. When compared to Providencia siddallii, an obligate nutritional symbiont of other Haementeria leeches, this novel bacterial symbiont shows convergent retention of the metabolic pathways involved in B vitamin biosynthesis. Moreover, the genomic characteristics of this Pluralibacter symbiont suggest a more recent association than that of Pr. siddallii and Haementeria. We conclude that the once-thought stable associations between blood-feeding Glossiphoniidae and their symbionts (i.e., one bacteriome structure, one symbiont lineage) can break down, mirroring symbiont turnover observed in various arthropod lineages.

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

fENko-Kae01 is a flagellum-specific jumbo phage infecting Klebsiella aerogenes

BACKGROUND

Klebsiella aerogenes is an opportunistic pathogen that causes a wide variety of infections. Due to the rising problem of antibiotic resistance, novel antibiotics and strategies to combat bacterial infections are needed. Host-specific bacteriophages are natural enemies of bacteria and can be used in phage therapy as an alternative form of treatment against bacterial infections. Jumbo phages are defined as phages with genomes larger than 200 kb. Relatively few studies have been done on jumbo phages compared to smaller phages.

RESULTS

A novel phage, fENko-Kae01, was isolated from a commercial phage cocktail. Genomic analysis revealed that fENko-Kae01 is a lytic jumbo phage with a 360 kb genome encoding 578 predicted genes. No highly similar phage genomes were identified and fENko-Kae01 may be a completely new genus representative. No known genes associated with lysogenic life cycle, bacterial virulence, or antibiotic resistance were identified. The phage had myovirus morphology and a narrow host range. Phage resistant bacterial mutants emerged under phage selection. Whole genome sequencing revealed that the biogenesis of the flagellum was affected in four mutants and the lack of functional flagellum was confirmed in motility assays. Furthermore, phage fENKo-Kae01 failed to adsorb on the non-motile mutants indicating that the bacterial flagellum is the phage-binding receptor.

CONCLUSIONS

fENko-Kae01 is a novel jumbo bacteriophage that is considered safe for phage therapy. fENko-Kae01 uses the flagellum as the phage-binding receptor and may represent a completely novel genus.

Facultatively intrabacterial localization of a planthopper endosymbiont as an adaptation to its vertical transmission

Transovarial transmission is the most reliable way of passing on essential nutrient-providing endosymbionts from mothers to offspring. However, not all endosymbiotic microbes follow the complex path through the female host tissues to oocytes on their own. Here, we demonstrate an unusual transmission strategy adopted by one of the endosymbionts of the planthopper Trypetimorpha occidentalis (Hemiptera: Tropiduchidae) from Bulgaria. In this species, an Acetobacteraceae endosymbiont is transmitted transovarially within deep invaginations of cellular membranes of an ancient endosymbiont Sulcia-strikingly resembling recently described plant virus transmission. However, in males, Acetobacteraceae colonizes the same bacteriocytes as Sulcia but remains unenveloped. Then, the unusual endobacterial localization of Acetobacteraceae observed in females appears to be a unique adaptation to maternal transmission. Further, the symbiont's genomic features, including encoding essential amino acid biosynthetic pathways and its similarity to a recently described psyllid symbiont, suggest a unique combination of the ability to horizontally transmit among species and confer nutritional benefits. The close association with Acetobacteraceae symbiont correlates with the so-far-unreported level of genomic erosion of ancient nutritional symbionts of this planthopper. In Sulcia, this is reflected in substantial changes in genomic organization, reported for the first time in the symbiont renowned for its genomic stability. In Vidania, substantial gene loss resulted in one of the smallest genomes known, at 108.6 kb. Thus, the symbionts of T. occidentalis display a combination of unusual adaptations and genomic features that expand our understanding of how insect-microbe symbioses may transmit and evolve.IMPORTANCEReliable transmission across host generations is a major challenge for bacteria that associate with insects, and independently established symbionts have addressed this challenge in different ways. The facultatively endobacterial localization of Acetobacteraceae symbiont, enveloped by cells of ancient nutritional endosymbiont Sulcia in females but not males of the planthopper Trypetimorpha occidentalis, appears to be a unique adaptation to maternal transmission. Acetobacteraceae's genomic features indicate its unusual evolutionary history, and the genomic erosion experienced by ancient nutritional symbionts demonstrates the apparent consequences of such close association. Combined, this multi-partite symbiosis expands our understanding of the diversity of strategies that insect symbioses form and some of their evolutionary consequences.

Novel phage infecting the Roseobacter CHUG lineage reveals a diverse and globally distributed phage family

Bacteriophages play an essential role in shaping the diversity and metabolism of bacterial communities. Marine Roseobacter group is an abundant heterotrophic bacterial group that is involved in many major element cycles, especially carbon and sulfur. Members of the Roseobacter CHUG (Clade Hidden and Underappreciated Globally) lineage are globally distributed and are activated in pelagic marine environments. In this study, we isolated and characterized a phage, CRP-810, that infects the CHUG strain FZCC0198. The genome of CRP-810 was dissimilar to those of other known phages. Additionally, 251 uncultured viral genomes (UViGs) closely related to CRP-810 were obtained from the uncultivated marine viral contig databases. Comparative genomic and phylogenetic analyses revealed that CRP-810 and these related UViGs exhibited conserved genome synteny, representing a new phage family with at least eight subgroups. Most of the CRP-810-type phages contain an integrase gene, and CRP-810 can be integrated into the host genome. Further analysis revealed that three CRP-810-type members were prophages found in the genomes of marine SAR11, Poseidonocella, and Sphingomonadaceae. Finally, viromic read-mapping analysis showed that CRP-810-type phages were globally distributed and displayed distinct biogeographic patterns related to temperature and latitude. Many members with a lower G + C content were mainly distributed in the trade station, whereas members with a higher G + C content were mainly distributed in polar and westerlies station, indicating that the niche differentiation of phages was subject to host adaptation. Collectively, these findings identify a novel phage family and expand our understanding of phylogenetic diversity, evolution, and biogeography of marine phages.

IMPORTANCE

The Roseobacter CHUG lineage, affiliated with the Pelagic Roseobacter Cluster (PRC), is widely distributed in the global oceans and is active in oligotrophic seawater. However, knowledge of the bacteriophages that infect CHUG members is limited. In this study, a CHUG phage, CRP-810, that infects the CHUG strain FZCC0198, was isolated and shown to have a novel genomic architecture. In addition, 251 uncultured viral genomes closely related to CRP-810 were recovered and included in the analyses. Phylogenomic analyses revealed that the CRP-810-type phages represent a new phage family containing at least eight genus-level subgroups. Members of this family were predicted to infect various marine bacteria. We also demonstrated that the CRP-810-type phages are widely distributed in global oceans and display distinct biogeographic patterns related to latitude. Collectively, this study provides important insights into the genomic organization, diversity, and ecology of a novel phage family that infect ecologically important bacteria in the global ocean.

Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome

KEY MESSAGE

We reported the mitochondrial genome of Cinnamomum camphora for the first time, revealing frequent rearrangement events in the non-coding regions of Magnoliids mitochondrial genomes. As one of the representative species in the Lauraceae family of Magnoliids, Cinnamomum camphora holds significant economic and ecological value. In this study, the mitochondrial genome (mitogenome) of C. camphora was complete assembled and annotated using PacBio HiFi sequencing. The C. camphora mitogenome is characterized by a branch structure, spans 900,894 bp, and contains 43 protein-coding genes (PCGs), 24 tRNAs, and 3 rRNAs. Most of these PCGs are under purifying selection, with only two (ccmFc and rps7) exhibiting signs of positive selection. The C. camphora mitogenome contains numerous repetitive sequences and intracellular gene transfers, with a total of 36 mitochondrial plastid DNAs, amounting to a combined length of 23,816 bp. Comparative analysis revealed that the non-coding regions of Magnoliids mitogenomes have undergone frequent rearrangements during evolution, but the coding sequences remain highly conserved (more than 98% similarity for protein-coding sequences). Furthermore, a maximum-likelihood phylogenetic tree was reconstructed based on 25 PCGs from 23 plant mitogenomes. The analysis supports the closest relationship between C. camphora and C. chekiangense, consistent with the APG IV classification system. This study elucidates the unique evolutionary features of the C. camphora mitogenome, which will provide valuable insights into the study of genetics and evolution of the family Lauraceae.

Nuclear and mitochondrial genomes of the plum fruit moth Grapholita funebrana

The plum fruit moth Grapholita funebrana (Tortricidae, Lepidoptera) is an important pest of many wild and cultivated stone fruits and other plants in the family Rosaceae. Here, we assembled its nuclear and mitochondrial genomes using Illumina, Nanopore, and Hi-C sequencing technologies. The nuclear genome size is 570.9 Mb, with a repeat rate of 51.28%, and a BUCSO completeness of 97.7%. The karyotype for males is 2n = 56. We identified 17,979 protein-coding genes, 5,643 tRNAs, and 94 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and 2 rRNA. These genomes provide resources to understand the genetics, ecology, and genome evolution of the tortricid moths.

Complete chloroplast genomes of eight Delphinium taxa (Ranunculaceae) endemic to Xinjiang, China: insights into genome structure, comparative analysis, and phylogenetic relationships

BACKGROUND

Delphinium L. represents a taxonomically intricate genus of significant phylogenetic and economic importance in Ranunculaceae. Despite the existence of few chloroplast genome datasets, a comprehensive understanding of genome structures and selective pressures within the genus remains unknown. Furthermore, several taxa in this genus are exclusively found in Xinjiang, China, a region renowned for its distribution and diversity of Chinese and Central Asian Delphinium species. Therefore, investigating the features of chloroplast genomes in this area will provide valuable insights into the evolutionary processes and phylogenetic relationships of the genus.

RESULTS

In this study, the eight newly completed chloroplast genomes are examined, ranging in length from 153,979 bp to 154,284 bp. Alongside these, analysing six previously reported taxa re-annotated in Delphinium, 111 unique genes are identified across all samples. Genome structure, distributions of simple sequence repeats and short dispersed repeats, as well as gene content are similar among these Delphinium taxa. Nine hypervariable intergenic spacers and protein coding regions, including ndhF-trnL(TAG), rpl16-intron, rpl33, rps15, rps18, trnK(TTT)-trnQ(TTG), trnP(TGG)-psaJ, trnT(GGT)-psbD and ycf1, are identified among 13 perennial Delphinium. Selective pressure and codon usage bias of all the plastid genes are performed within 14 Delphinium taxa. Phylogenetic analysis based on 14 Delphinium plastomes, alongside two Aconitum (Ranunculaceae) species serving as outgroup taxa, reveals the monophyletic nature of Delphinium. Our findings further discern Delphinium into two distinct clades: perennial species (clade I) and annual species (clade II). In addition, compared with the nrDNA ITS topology, cytological data and morphological characters, D. mollifolium and D. maackianum showed potential involvement in hybridization or polyploidization processes. Excluding these two species, the perennial Delphinium (clade I) exhibits a stronger consistency with the morphology-based system that utilized seed morphology.

CONCLUSION

This study represents the first comprehensive analysis of plastomic variations among Delphinium taxa, based on the examination of 14 complete plastomes. The chloroplast genome structure of Delphinium is similar to other angiosperms and possesses the typical quadripartite structure with the conserved genome arrangement and gene features. In addition, the variation of non-coding regions is larger than coding regions of the chloroplast genome. Through DNA sequence divergence across Delphinium plastomes and subsequent phylogenomic analyses ndhF-trnL(TAG) and ycf1 are identified as promising molecular markers. These highly variable loci held significant potential for future phylogenetic and phylogeographic studies on Delphinium. Our phylogenomic analyses based on the whole plastomes, concatenation of 132 unique intergenic spacer regions, concatenation of 77 unique protein-coding genes and nrDNA ITS, all support the monophyly of Delphinium and perennial taxa clusters together into one clade within this genus. These findings provide crucial data for systematic, phylogenomic and evolutionary research in the genus for future studies.

Integration of large and diverse angiosperm DNA fragments into Asian Gnetum mitogenomes

BACKGROUND

Horizontal gene transfer (HGT) events have rarely been reported in gymnosperms. Gnetum is a gymnosperm genus comprising 25‒35 species sympatric with angiosperms in West African, South American, and Southeast Asian rainforests. Only a single acquisition of an angiosperm mitochondrial intron has been documented to date in Asian Gnetum mitogenomes. We wanted to develop a more comprehensive understanding of frequency and fragment length distribution of such events as well as their evolutionary history in this genus.

RESULTS

We sequenced and assembled mitogenomes from five Asian Gnetum species. These genomes vary remarkably in size and foreign DNA content. We identified 15 mitochondrion-derived and five plastid-derived (MTPT) foreign genes. Our phylogenetic analyses strongly indicate that these foreign genes were transferred from diverse eudicots-mostly from the Rubiaceae genus Coptosapelta and ten genera of Malpighiales. This indicates that Asian Gnetum has experienced multiple independent HGT events. Patterns of sequence evolution strongly suggest DNA-mediated transfer between mitochondria as the primary mechanism giving rise to these HGT events. Most Asian Gnetum species are lianas and often entwined with sympatric angiosperms. We therefore propose that close apposition of Gnetum and angiosperm stems presents opportunities for interspecific cell-to-cell contact through friction and wounding, leading to HGT.

CONCLUSIONS

Our study reveals that multiple HGT events have resulted in massive amounts of angiosperm mitochondrial DNA integrated into Asian Gnetum mitogenomes. Gnetum and its neighboring angiosperms are often entwined with each other, possibly accounting for frequent HGT between these two phylogenetically remote lineages.

The complete mitochondrial genome of a jujube geometrid, Sucra jujuba (Lepidoptera: Geometridae) and its phylogenetic analysis

Sucra jujuba Chu, 1979 (Lepidoptera: Geometridae) is a major insect pest in jujube plantation. In this study, we have sequenced the complete mitochondrial genome of S. jujuba. The circular genome was 15,557 bp in length and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and one AT-rich region (GenBank accession no. MZ507574). The nucleotide composition was significantly biased (A, T, C, and G were 41.85%, 39.65%, 10.97%, and 7.53%, respectively) with A + T contents of 81.50%. The Bayesian phylogenetic analysis of the concatenated nucleotide sequences of 13 PCGs from 30 species in the subfamily Ennominae and two outgroup species was performed. The results indicated that S. jujuba was closely related to Amraica recursaria in the subfamily Ennominae.

A multi-kingdom collection of 33,804 reference genomes for the human vaginal microbiome

The human vagina harbours diverse microorganisms-bacteria, viruses and fungi-with profound implications for women's health. Genome-level analysis of the vaginal microbiome across multiple kingdoms remains limited. Here we utilize metagenomic sequencing data and fungal cultivation to establish the Vaginal Microbial Genome Collection (VMGC), comprising 33,804 microbial genomes spanning 786 prokaryotic species, 11 fungal species and 4,263 viral operational taxonomic units. Notably, over 25% of prokaryotic species and 85% of viral operational taxonomic units remain uncultured. This collection significantly enriches genomic diversity, especially for prevalent vaginal pathogens such as BVAB1 (an uncultured bacterial vaginosis-associated bacterium) and Amygdalobacter spp. (BVAB2 and related species). Leveraging VMGC, we characterize functional traits of prokaryotes, notably Saccharofermentanales (an underexplored yet prevalent order), along with prokaryotic and eukaryotic viruses, offering insights into their niche adaptation and potential roles in the vagina. VMGC serves as a valuable resource for studying vaginal microbiota and its impact on vaginal health.

A Genomics-Based Discovery of Secondary Metabolite Biosynthetic Gene Clusters in the Potential Novel Strain Streptomyces sp. 21So2-11 Isolated from Antarctic Soil

Streptomyces species are attractive sources of secondary metabolites that serve as major sources of antibiotics and other drugs. In this study, genome mining was used to determine the biosynthetic potential of Streptomyces sp. 21So2-11 isolated from Antarctic soil. 16S rRNA gene sequencing revealed that this strain is most closely related to Streptomyces drozdowiczii NBRC 101007T, with a similarity of 98.02%. Genome comparisons based on average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) showed that strain 21So2-11 represents a novel species of the genus Streptomyces. In addition to a large number of genes related to environmental adaptation and ecological function, a total of 28 putative biosynthetic gene clusters (BGCs) responsible for the biosynthesis of known and/or novel secondary metabolites, including terpenes, lantipeptides, polyketides, nonribosomal peptides, RiPPs and siderophores, were detected in the genome of strain 21So2-11. In addition, a total of 1456 BGCs were predicted to contribute to the biosynthesis of more than 300 secondary metabolites based on the genomes of 47 Streptomyces strains originating from polar regions. The results indicate the potential of Streptomyces sp. 21So2-11 for bioactive secondary metabolite production and are helpful for understanding bacterial adaptability and ecological function in cold terrestrial environments.

Temporal turnover of viral biodiversity and functional potential in intertidal wetlands

As the central members of the microbiome networks, viruses regulate the composition of microbial communities and drive the nutrient cycles of ecosystems by lysing host cells. Therefore, uncovering the dynamic patterns and the underlying ecological mechanisms mediating the tiniest viral communities across space and through time in natural ecosystems is of crucial importance for better understanding the complex microbial world. Here, the temporal dynamics of intertidal viral communities were investigated via a time-series sampling effort. A total of 1911 viral operational taxonomic units were recovered from 36 bimonthly collected shotgun metagenomes. Functionally important auxiliary metabolic genes involved in carbohydrate, sulfur, and phosphorus metabolism were detected, some of which (e.g., cysH gene) were stably present within viral genomes over time. Over the sampling period, strong and comparable temporal turnovers were observed for intertidal viromes and their host microbes. Winter was determined as the pivotal point for the shifts in viral diversity patterns. Notably, the viral micro-diversity covaried with the macro-diversity, following similar temporal patterns. The relative abundances of viral taxa also covaried with their host prokaryotes. Meanwhile, the virus-host relationships at the whole community level were relatively stable. Further statistical analyses demonstrated that the dynamic patterns of viral communities were highly deterministic, for which temperature was the major driver. This study provided valuable mechanistic insights into the temporal turnover of viral communities in complex ecosystems such as intertidal wetlands.

Characterization and Genome Analysis of the Delftia lacustris Strain LzhVag01 Isolated from Vaginal Discharge

Delftia has been separated from freshwater, sludge, and soil and has emerged as a novel opportunistic pathogen in the female vagina. However, the genomic characteristics, pathogenicity, and biotechnological properties still need to be comprehensively investigated. In this study, a Delftia strain was isolated from the vaginal discharge of a 43-year-old female with histologically confirmed cervical intraepithelial neoplasm (CIN III), followed by whole-genome sequencing. Phylogenetic analysis and average nucleotide identity (ANI) analysis demonstrated that it belongs to Delftia lacustris, named D. lacustris strain LzhVag01. LzhVag01 was sensitive to β-lactams, macrolides, and tetracyclines but exhibited resistance to lincoamines, nitroimidazoles, aminoglycosides, and fluoroquinolones. Its genome is a single, circular chromosome of 6,740,460 bp with an average GC content of 66.59%. Whole-genome analysis identified 16 antibiotic resistance-related genes, which match the antimicrobial susceptibility profile of this strain, and 11 potential virulence genes. These pathogenic factors may contribute to its colonization in the vaginal environment and its adaptation and accelerate the progression of cervical cancer. This study sequenced and characterized the whole-genome of Delftia lacustris isolated from vaginal discharge, which provides investigators and clinicians with valuable insights into this uncommon species.

In-depth exploration of the genomic diversity in tea varieties based on a newly constructed pangenome of Camellia sinensis

Tea, one of the most widely consumed beverages globally, exhibits remarkable genomic diversity in its underlying flavour and health-related compounds. In this study, we present the construction and analysis of a tea pangenome comprising a total of 11 genomes, with a focus on three newly sequenced genomes comprising the purple-leaved assamica cultivar "Zijuan", the temperature-sensitive sinensis cultivar "Anjibaicha" and the wild accession "L618" whose assemblies exhibited excellent quality scores as they profited from latest sequencing technologies. Our analysis incorporates a detailed investigation of transposon complement across the tea pangenome, revealing shared patterns of transposon distribution among the studied genomes and improved transposon resolution with long read technologies, as shown by long terminal repeat (LTR) Assembly Index analysis. Furthermore, our study encompasses a gene-centric exploration of the pangenome, exploring the genomic landscape of the catechin pathway with our study, providing insights on copy number alterations and gene-centric variants, especially for Anthocyanidin synthases. We constructed a gene-centric pangenome by structurally and functionally annotating all available genomes using an identical pipeline, which both increased gene completeness and allowed for a high functional annotation rate. This improved and consistently annotated gene set will allow for a better comparison between tea genomes. We used this improved pangenome to capture the core and dispensable gene repertoire, elucidating the functional diversity present within the tea species. This pangenome resource might serve as a valuable resource for understanding the fundamental genetic basis of traits such as flavour, stress tolerance, and disease resistance, with implications for tea breeding programmes.

Chromosome-level genome assembly and annotation of a potential model organism Gossypium arboreum ZB-1

Recent advancements in plant regeneration and synthetic polyploid creation have been documented in Gossypium arboreum ZB-1. These developments make ZB-1 a potential model within the Gossypium genus for investigating gene function and polyploidy. This work generated the sequence and annotation of the ZB-1 genome. The contig-level genome was constructed using the PacBio high-fidelity reads, encompassing 81 contigs with an N50 length of 112.12 Mb. The Hi-C data assisted the construction of the chromosome-level genome, which consists of 13 pseudo-chromosomes and 39 un-anchored contigs, with a total length of about 1.67 Gb. Repetitive sequences accounted for about 69.7% of the genome in length. Based on ab initio and evidence-based prediction, we have identified 48,021 protein-coding genes in the ZB-1 genome. Comparative genomics analysis revealed conserved gene content and arrangement between ZB-1 and G. arboreum SXY1. The single nucleotide polymorphism occurrence rate between ZB-1 and SXY1 was about 0.54 per 1,000 nucleotides. This study enriched the genomic resources for further exploration into cotton regeneration and polyploidy mechanisms.

In vivo efficacy of phage cocktails against carbapenem resistance Acinetobacter baumannii in the rat pneumonia model

Acinetobacter baumannii, an opportunistic pathogen, poses a significant threat in intensive care units, leading to severe nosocomial infections. The rise of multi-drug-resistant strains, particularly carbapenem-resistant A. baumannii, has created formidable challenges for effective treatment. Given the prolonged development cycle and high costs associated with antibiotics, phages have garnered clinical attention as an alternative for combating infections caused by drug-resistant bacteria. However, the utilization of phage therapy encounters notable challenges, including the narrow host spectrum, where each phage targets a limited subset of bacteria, increasing the risk of phage resistance development. Additionally, uncertainties in immune system dynamics during treatment hinder tailoring symptomatic interventions based on patient-specific states. In this study, we isolated two A. baumannii phages from wastewater and conducted a comprehensive assessment of their potential applications. This evaluation included sequencing analysis, genome classification, pH and temperature stability assessments, and in vitro bacterial inhibition assays. Further investigations involved analyzing histological and cytokine alterations in rats undergoing phage cocktail treatment for pneumonia. The therapeutic efficacy of the phages was validated, and transcriptomic studies of rat lung tissue during phage treatment revealed crucial changes in the immune system. The findings from our study underscore the potential of phages for future development as a treatment strategy and offer compelling evidence regarding immune system dynamics throughout the treatment process.IMPORTANCEDue to the growing problem of multi-drug-resistant bacteria, the use of phages is being considered as an alternative to antibiotics, and the genetic safety and application stability of phages determine the potential of phage application. The absence of drug resistance genes and virulence genes in the phage genome can ensure the safety of phage application, and the fact that phage can remain active in a wide range of temperatures and pH is also necessary for application. In addition, the effect evaluation of preclinical studies is especially important for clinical application. By simulating the immune response situation during the treatment process through mammalian models, the changes in animal immunity can be observed, and the effect of phage therapy can be further evaluated. Our study provides compelling evidence that phages hold promise for further development as therapeutic agents for Acinetobacter baumannii infections.

Genome assembly of autotetraploid Actinidia arguta highlights adaptive evolution and enables dissection of important economic traits

Actinidia arguta, the most widely distributed Actinidia species and the second cultivated species in the genus, can be distinguished from the currently cultivated Actinidia chinensis on the basis of its small and smooth fruit, rapid softening, and excellent cold tolerance. Adaptive evolution of tetraploid Actinidia species and the genetic basis of their important agronomic traits are still unclear. Here, we generated a chromosome-scale genome assembly of an autotetraploid male A. arguta accession. The genome assembly was 2.77 Gb in length with a contig N50 of 9.97 Mb and was anchored onto 116 pseudo-chromosomes. Resequencing and clustering of 101 geographically representative accessions showed that they could be divided into two geographic groups, Southern and Northern, which first diverged 12.9 million years ago. A. arguta underwent two prominent expansions and one demographic bottleneck from the mid-Pleistocene climate transition to the late Pleistocene. Population genomics studies using paleoclimate data enabled us to discern the evolution of the species' adaptation to different historical environments. Three genes (AaCEL1, AaPME1, and AaDOF1) related to flesh softening were identified by multi-omics analysis, and their ability to accelerate flesh softening was verified through transient expression assays. A set of genes that characteristically regulate sexual dimorphism located on the sex chromosome (Chr3) or autosomal chromosomes showed biased expression during stamen or carpel development. This chromosome-level assembly of the autotetraploid A. arguta genome and the genes related to important agronomic traits will facilitate future functional genomics research and improvement of A. arguta.

Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica

MAIN CONCLUSION

In this study, we assembled the first complete mitochondrial genome of Setaria italica and confirmed the multi-branched architecture. The foxtail millet (Setaria italica) holds significant agricultural importance, particularly in arid and semi-arid regions. It plays a pivotal role in diversifying dietary patterns and shaping planting strategies. Although the chloroplast genome of S. italica has been elucidated in recent studies, the complete mitochondrial genome remains largely unexplored. In this study, we employed PacBio HiFi sequencing platforms to sequence and assemble the complete mitochondrial genome. The mitochondrial genome spans a total length of 446,614 base pairs and harbors a comprehensive set of genetic elements, including 33 unique protein-coding genes (PCGs), encompassing 24 unique mitochondrial core genes and 9 variable genes, along with 20 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. Our analysis of mitochondrial PCGs revealed a pronounced codon usage preference. For instance, the termination codon exhibits a marked preference for UAA, while alanine (Ala) exhibits a preference for GCU, and glutamine (Gln) favors CAA. Notably, the maximum Relative Synonymous Codon Usage (RSCU) values for cysteine (Cys) and phenylalanine (Phe) are both below 1.2, indicating a lack of strong codon usage preference for these amino acids. Phylogenetic analyses consistently place S. italica in close evolutionary proximity to Chrysopogon zizanioides, relative to other Panicoideae plants. Collinearity analysis showed that a total of 39 fragments were identified to display homology with both the mitochondrial and chloroplast genomes. A total of 417 potential RNA-editing sites were discovered across the 33 mitochondrial PCGs. Notably, all these editing events involved the conversion of cytosine (C) to uracil (U). Through the employment of PCR validation coupled with Sanger sequencing for the anticipated editing sites of these codons, RNA-editing events were conclusively identified at two specific loci: nad4L-2 and atp6-1030. The results of this study provide a pivotal foundation for advanced genomic breeding research in foxtail millet. Furthermore, they impart essential insights that will be instrumental for forthcoming investigations into the evolutionary and molecular dynamics of Panicoideae species.

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.

Haplotype-resolved chromosome-level genome assembly of Ehretia macrophylla

Ehretia macrophylla Wall, known as wild loquat, is an ecologically, economically, and medicinally significant tree species widely grown in China, Japan, Vietnam, and Nepal. In this study, we have successfully generated a haplotype-resolved chromosome-scale genome assembly of E. macrophylla by integrating PacBio HiFi long-reads, Illumina short-reads, and Hi-C data. The genome assembly consists of two haplotypes, with sizes of 1.82 Gb and 1.58 Gb respectively, and contig N50 lengths of 28.11 Mb and 21.57 Mb correspondingly. Additionally, 99.41% of the assembly was successfully anchored into 40 pseudo-chromosomes. We predicted 58,886 protein-coding genes, of which 99.60% were functionally annotated from databases. We furthermore detected 2.65 Gb repeat sequences, 659,290 rRNAs, 4,931 tRNAs and 4,688 other ncRNAs. The high-quality assembly of the genome offers a solid basis for furthering the fields of molecular breeding and functional genomics of E. macrophylla.

Mitogenomic Insights into the Evolution, Divergence Time, and Ancestral Ranges of Coturnix Quails

The Old-World quails, Coturnix coturnix (common quail) and Coturnix japonica (Japanese quail), are morphologically similar yet occupy distinct geographic ranges. This study aimed to elucidate their evolutionary trajectory and ancestral distribution patterns through a thorough analysis of their mitochondrial genomes. Mitogenomic analysis revealed high structural conservation, identical translational mechanisms, and similar evolutionary pressures in both species. Selection analysis revealed significant evidence of positive selection across the Coturnix lineage for the nad4 gene tree owing to environmental changes and acclimatization requirements during its evolutionary history. Divergence time estimations imply that diversification among Coturnix species occurred in the mid-Miocene (13.89 Ma), and their current distributions were primarily shaped by dispersal rather than global vicariance events. Phylogenetic analysis indicates a close relationship between C. coturnix and C. japonica, with divergence estimated at 2.25 Ma during the Pleistocene epoch. Ancestral range reconstructions indicate that the ancestors of the Coturnix clade were distributed over the Oriental region. C. coturnix subsequently dispersed to Eurasia and Africa, and C. japonica to eastern Asia. We hypothesize that the current geographic distributions of C. coturnix and C. japonica result from their unique dispersal strategies, developed to evade interspecific territoriality and influenced by the Tibetan Plateau's geographic constraints. This study advances our understanding of the biogeographic and evolutionary processes leading to the diversification of C. coturnix and C. japonica, laying important groundwork for further research on this genus.

Mitochondrial genome complexity in Stemona sessilifolia: nanopore sequencing reveals chloroplast gene transfer and DNA rearrangements

Mitochondria are semi-autonomous organelles in eukaryotic cells with their own genome. Plant mitogenomes differ from animal mitogenomes in size, structure, and repetitive DNA sequences. Despite larger sizes, plant mitogenomes do not have significantly more genes. They exhibit diverse structures due to variations in size, repetitive DNA, recombination frequencies, low gene densities, and reduced nucleotide substitution rates. In this study, we analyzed the mitochondrial genome of Stemona sessilifolia using Nanopore and Illumina sequencing. De-novo assembly and annotation were conducted using Unicycler, Geseq, tRNAscan-SE and BLASTN, followed by codon usage, repeat sequence, RNA-editing, synteny, and phylogenetic analyses. S. sessilifolia's mitogenome consisted of one linear contig and six circular contigs totaling 724,751 bp. It had 39 protein-coding genes, 27 tRNA genes, and 3 rRNA genes. Transfer of chloroplast sequences accounted for 13.14% of the mitogenome. Various analyses provided insights into genetic characteristics, evolutionary dynamics, and phylogenetic placement. Further investigations can explore transferred genes' functions and RNA-editing's role in mitochondrial gene expression in S. sessilifolia.

Chromosome-level genome assembly of Solanum pimpinellifolium

Solanum pimpinellifolium, the closest wild relative of the domesticated tomato, has high potential for use in breeding programs aimed at developing multi-pathogen resistance and quality improvement. We generated a chromosome-level genome assembly of S. pimpinellifolium LA1589, with a size of 833 Mb and a contig N50 of 31 Mb. We anchored 98.80% of the contigs into 12 pseudo-chromosomes, and identified 74.47% of the sequences as repetitive sequences. The genome evaluation revealed BUSCO and LAI score of 98.3% and 14.49, respectively, indicating high quality of this assembly. A total of 41,449 protein-coding genes were predicted in the genome, of which 89.17% were functionally annotated. This high-quality genome assembly serves as a valuable resource for accelerating the biological discovery and molecular breeding of this important horticultural crop.

Wolbachia symbionts control sex in a parasitoid wasp using a horizontally acquired gene

Host reproduction can be manipulated by bacterial symbionts in various ways. Parthenogenesis induction is the most effective type of reproduction manipulation by symbionts for their transmission. Insect sex is determined by regulation of doublesex (dsx) splicing through transformer2 (tra2) and transformer (tra) interaction. Although parthenogenesis induction by symbionts has been studied since the 1970s, its underlying molecular mechanism is unknown. Here we identify a Wolbachia parthenogenesis-induction feminization factor gene (piff) that targets sex-determining genes and causes female-producing parthenogenesis in the haplodiploid parasitoid Encarsia formosa. We found that Wolbachia elimination repressed expression of female-specific dsx and enhanced expression of male-specific dsx, which led to the production of wasp haploid male offspring. Furthermore, we found that E. formosa tra is truncated and non-functional, and Wolbachia has a functional tra homolog, termed piff, with an insect origin. Wolbachia PIFF can colocalize and interact with wasp TRA2. Moreover, Wolbachia piff has coordinated expression with tra2 and dsx of E. formosa. Our results demonstrate the bacterial symbiont Wolbachia has acquired an insect gene to manipulate the host sex determination cascade and induce parthenogenesis in wasps. This study reveals insect-to-bacteria horizontal gene transfer drives the evolution of animal sex determination systems, elucidating a striking mechanism of insect-microbe symbiosis.

Seed coating with phages for sustainable plant biocontrol of plant pathogens and influence of the seed coat mucilage

Pathogens resistant to classical control strategies pose a significant threat to crop yield, with seeds being a major transmission route. Bacteriophages, viruses targeting bacteria, offer an environmentally sustainable biocontrol solution. In this study, we isolated and characterized two novel phages, Athelas and Alfirin, which infect Pseudomonas syringae and Agrobacterium fabrum, respectively, and included the recently published Pfeifenkraut phage infecting Xanthomonas translucens. Using a simple immersion method, phages coated onto seeds successfully lysed bacteria post air-drying. The seed coat mucilage (SCM), a polysaccharide-polymer matrix exuded by seeds, plays a critical role in phage binding. Seeds with removed mucilage formed five to 10 times less lysis zones compared to those with mucilage. The podovirus Athelas showed the highest mucilage dependency. Phages from the Autographiviridae family also depended on mucilage for seed adhesion. Comparative analysis of Arabidopsis SCM mutants suggested the diffusible cellulose as a key component for phage binding. Long-term activity tests demonstrated high phage stability on seed surfaces and significantly increasing seedling survival rates in the presence of pathogens. Using non-virulent host strains enhanced phage presence on seeds but also has potential limitations. These findings highlight phage-based interventions as promising, sustainable strategies for combating pathogen resistance and improving crop yield.

Telomere-to-telomere genome assembly of Oldenlandia diffusa

We report the complete telomere-to-telomere genome assembly of Oldenlandia diffusa which renowned in traditional Chinese medicine, comprising 16 chromosomes and spanning 499.7 Mb. The assembly showcases 28 telomeres and minimal gaps, with a total of only five. Repeat sequences constitute 46.41% of the genome, and 49,701 potential protein-coding genes have been predicted. Compared with O. corymbosa, O. diffusa exhibits chromosome duplication and fusion events, diverging 20.34 million years ago. Additionally, a total of 11 clusters of terpene synthase have been identified. The comprehensive genome sequence, gene catalog, and terpene synthase clusters of O. diffusa detailed in this study will significantly contribute to advancing research in this species' genetic, genomic, and pharmacological aspects.

Engineering natural microbiomes toward enhanced bioremediation by microbiome modeling

Engineering natural microbiomes for biotechnological applications remains challenging, as metabolic interactions within microbiomes are largely unknown, and practical principles and tools for microbiome engineering are still lacking. Here, we present a combinatory top-down and bottom-up framework to engineer natural microbiomes for the construction of function-enhanced synthetic microbiomes. We show that application of herbicide and herbicide-degrader inoculation drives a convergent succession of different natural microbiomes toward functional microbiomes (e.g., enhanced bioremediation of herbicide-contaminated soils). We develop a metabolic modeling pipeline, SuperCC, that can be used to document metabolic interactions within microbiomes and to simulate the performances of different microbiomes. Using SuperCC, we construct bioremediation-enhanced synthetic microbiomes based on 18 keystone species identified from natural microbiomes. Our results highlight the importance of metabolic interactions in shaping microbiome functions and provide practical guidance for engineering natural microbiomes.

Complete genome assembly provides a high-quality skeleton for pan-NLRome construction in melon

Melon (Cucumis melo L.), being under intensive domestication and selective breeding, displays an abundant phenotypic diversity. Wild germplasm with tolerance to stress represents an untapped genetic resource for discovery of disease-resistance genes. To comprehensively characterize resistance genes in melon, we generate a telomere-to-telomere (T2T) and gap-free genome of wild melon accession PI511890 (C. melo var. chito) with a total length of 375.0 Mb and a contig N50 of 31.24 Mb. The complete genome allows us to dissect genome architecture and identify resistance gene analogs. We construct a pan-NLRome using seven melon genomes, which include 208 variable and 18 core nucleotide-binding leucine-rich repeat receptors (NLRs). Multiple disease-related transcriptome analyses indicate that most up-regulated NLRs induced by pathogens are shell or cloud NLRs. The T2T gap-free assembly and the pan-NLRome not only serve as essential resources for genomic studies and molecular breeding of melon but also provide insights into the genome architecture and NLR diversity.

Genomic and taxonomic evaluation of 38 Treponema prophage sequences

BACKGROUND

Despite Spirochetales being a ubiquitous and medically important order of bacteria infecting both humans and animals, there is extremely limited information regarding their bacteriophages. Of the genus Treponema, there is just a single reported characterised prophage.

RESULTS

We applied a bioinformatic approach on 24 previously published Treponema genomes to identify and characterise putative treponemal prophages. Thirteen of the genomes did not contain any detectable prophage regions. The remaining eleven contained 38 prophage sequences, with between one and eight putative prophages in each bacterial genome. The prophage regions ranged from 12.4 to 75.1 kb, with between 27 and 171 protein coding sequences. Phylogenetic analysis revealed that 24 of the prophages formed three distinct sequence clusters, identifying putative myoviral and siphoviral morphology. ViPTree analysis demonstrated that the identified sequences were novel when compared to known double stranded DNA bacteriophage genomes.

CONCLUSIONS

In this study, we have started to address the knowledge gap on treponeme bacteriophages by characterising 38 prophage sequences in 24 treponeme genomes. Using bioinformatic approaches, we have been able to identify and compare the prophage-like elements with respect to other bacteriophages, their gene content, and their potential to be a functional and inducible bacteriophage, which in turn can help focus our attention on specific prophages to investigate further.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.

Genomic characteristics of nine Nitrospirota metagenome-assembled genomes in deep-sea sediments from East Pacific polymetallic nodules zone

Previously studies have reported that MAGs (Metagenome-assembled genomes) belong to "Candidatus Manganitrophaceae" of phylum Nitrospirota with chemolithoautotrophic manganese oxidation potential exist in freshwater and hydrothermal environments. However, Nitrospirota members with chemolithoautotrophic manganese oxidation potential have not been reported in other marine environments. Through metagenomic sequencing, assembly and binning, nine metagenome-assembled genomes belonging to Nitrospirota are recovered from sediment of different depths in the polymetallic nodule area. Through the key functional genes annotation results, we find that these Nitrospirota have limited potential to oxidize organic carbon because of incomplete tricarboxylic acid cycle and most of them (6/9) have carbon dioxide fixation potential through different pathway (rTCA, WL or CBB). One MAG belongs to order Nitrospirales has the potential to use manganese oxidation to obtain energy for carbon fixation. In addition to manganese ions, the oxidation of inorganic nitrogen, sulfur, hydrogen and carbon monoxide may also provide energy for the growth of these Nitrospirota. In addition, different metal ion transport systems can help those Nitrospirota to resist heavy metal in sediment. Our work expands the understanding of the metabolic potential of Nitrospirota in sediment of polymetallic nodule region and may contributes to promoting the study of chemolithoautotrophic manganese oxidation.