A simple method to control over-alignment in the MAFFT multiple sequence alignment program

The complete plastome and phylogenetic analysis of Commelina benghalensis L.1753 (Commelinaceae)

Commelina benghalensis L. 1753, a member of the Commelinaceae family, holds significant medicinal and culinary value. This study represents the first documentation of the sequencing and assembly of the entire plastome of C. benghalensis. The genome spans a total length of 160,663 bp, exhibiting a conventional quadripartite architecture that comprises a large single-copy (LSC) region (87,750 bp), a small single-copy (SSC) region (18,417 bp), and two inverted repeats (IR) regions (both 27,248 bp). In its entirety, the C. benghalensis plastome encompasses 129 genes (with 108 being unique), incorporating 77 individual protein-coding genes, 37 unique tRNA genes, and four unique rRNA genes. Phylogenetic analysis revealed a close resemblance between C. benghalensis and C. communis. The sequencing of this plastome stands to expedite the development of molecular markers and significantly contribute to genetic assays involving this distinctive plant.

Phylogenomics and genetic analysis of solvent-producing Clostridium species

The genus Clostridium is a large and diverse group within the Bacillota (formerly Firmicutes), whose members can encode useful complex traits such as solvent production, gas-fermentation, and lignocellulose breakdown. We describe 270 genome sequences of solventogenic clostridia from a comprehensive industrial strain collection assembled by Professor David Jones that includes 194 C. beijerinckii, 57 C. saccharobutylicum, 4 C. saccharoperbutylacetonicum, 5 C. butyricum, 7 C. acetobutylicum, and 3 C. tetanomorphum genomes. We report methods, analyses and characterization for phylogeny, key attributes, core biosynthetic genes, secondary metabolites, plasmids, prophage/CRISPR diversity, cellulosomes and quorum sensing for the 6 species. The expanded genomic data described here will facilitate engineering of solvent-producing clostridia as well as non-model microorganisms with innately desirable traits. Sequences could be applied in conventional platform biocatalysts such as yeast or Escherichia coli for enhanced chemical production. Recently, gene sequences from this collection were used to engineer Clostridium autoethanogenum, a gas-fermenting autotrophic acetogen, for continuous acetone or isopropanol production, as well as butanol, butanoic acid, hexanol and hexanoic acid production.

Do birds play a role in the transmission of Toscana virus? Initial isolation results from birds in northernmost Türkiye

AIMS

Recent research has prioritized emerging and re-emerging diseases that affect human and animal health, particularly to describe how these diseases enter countries and determine their transmission cycles. Given that migratory birds play a significant role in spreading infections, the present study analysed their migration paths and specimens to investigate Orthoflavivirus, Orthonairovirus, Alphavirus and Phlebovirus in birds in Samsun province, Türkiye.

METHODS AND RESULTS

For these viruses, 312 samples from 56 birds were analysed using RT-PCR and qRT-PCR. Toscana virus (TOSV) was identified in 14 birds (four mallards, five partridges, four quails and one pigeon), representing 25% of the birds sampled. Genotype B was reported in all 14 birds. After inoculating the positive tissues in cell cultures, TOSV was isolated from the organs of pigeons, mallards and partridges.

CONCLUSIONS

This is the first time TOSV has been isolated in cell culture from birds and indicates that they may play a role in spreading TOSV in Türkiye. The results also suggest that TOSV might be carried between countries by migratory birds.

Complexity of avian evolution revealed by family-level genomes

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

A vast repertoire of secondary metabolites potentially influences community dynamics and biogeochemical processes in cold seeps

In deep-sea cold seeps, microbial communities thrive on the geological seepage of hydrocarbons and inorganic compounds, differing from photosynthetically driven ecosystems. However, their biosynthetic capabilities remain largely unexplored. Here, we analyzed 81 metagenomes, 33 metatranscriptomes, and 7 metabolomes derived from nine different cold seep areas to investigate their secondary metabolites. Cold seep microbiomes encode diverse and abundant biosynthetic gene clusters (BGCs). Most BGCs are affiliated with understudied bacteria and archaea, including key mediators of methane and sulfur cycling. The BGCs encode diverse antimicrobial compounds that potentially shape community dynamics and various metabolites predicted to influence biogeochemical cycling. BGCs from key players are widely distributed and highly expressed, with their abundance and expression levels varying with sediment depth. Sediment metabolomics reveals unique natural products, highlighting uncharted chemical potential and confirming BGC activity in these sediments. Overall, these results demonstrate that cold seep sediments serve as a reservoir of hidden natural products and sheds light on microbial adaptation in chemosynthetically driven ecosystems.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Species clarification of fairy inkcap ("Coprinellus disseminatus") in China

Coprinellus disseminatus and other morphologically similar species are widely dispersed worldwide and are commonly referred to as "fairy inkcap". Based on the molecular phylogenetic study and morphological observation, a thorough investigation was carried out utilising 74 collections of related species that were gathered from seventeen provinces and five Chinese fungaria between 1998 and 2023 and revealed 11 lineages of "fairy inkcap", nine of which were found in China, and which belonged to the two genera Coprinellus and Tulosesus. In sect. Disseminati, genetic diversities (π), and fixation index (Fst) amongst lineages were computed, and a haplotype-based network was established to ascertain the relationships amongst each clade. A new section of Coprinellus, sect. Aureodisseminati, were discovered. In addition, four new species (C. aureodisseminatus, C. austrodisseminatus, C. parcus, and C. velutipes), a new subspecies of C. disseminatus, a new combination (Tulosesus pseudodisseminatus), the first discovery of epigamous type of C. magnoliae and a new record to China (T. subdisseminatus) were also identified and thoroughly described with accompanying illustrations. Their differences in macro- and micro-features, as well as their character sequence, were discussed.

Characterization of the complete mitochondrial genome and phylogenetic analyses of Haemaphysalis tibetensis Hoogstraal, 1965 (Acari: Ixodidae)

Ticks are specialized ectoparasites that feed on blood, causing physical harm to the host and facilitating pathogen transmission. The genus Haemaphysalis contains vectors for numerous infectious agents. These agents cause various diseases in humans and animals. Mitochondrial genome sequences serve as reliable molecular markers, forming a crucial basis for evolutionary analyses, studying species origins, and exploring molecular phylogeny. We extracted mitochondrial genome from the enriched mitochondria of Haemaphysalis tibetensis and obtained a 14,714-bp sequence. The mitochondrial genome consists of 13 protein-coding genes (PCGs), two ribosomal RNA, 22 transfer RNAs (tRNAs), and two control regions. The nucleotide composition of H. tibetensis mitochondrial genome was 38.38 % for A, 9.61 % for G, 39.32 % for T, and 12.69 % for C. The A + T content of H. tibetensis mitochondrial genome was 77.7 %, significantly higher than the G + C content. The repeat units of H. tibetensis exhibited two identical repeat units of 33 bp in length, positioned downstream of nad1 and rrnL genes. Furthermore, phylogenetic analyses based on the 13 PCGs indicated that Haemaphysalis tibetensis (subgenus Allophysalis) formed a monophyletic clade with Haemaphysalis nepalensis (subgenus Herpetobia) and Haemaphysalis danieli (subgenus Allophysalis). Although the species Haemaphysalis inermis, Haemaphysalis kitaokai, Haemaphysalis kolonini, and Haemaphysalis colasbelcouri belong to the subgenus Alloceraea, which were morphologically primitive hemaphysalines just like H. tibetensis, these four tick species cannot form a single clade with H. tibetensis. In this study, the whole mitochondrial genome sequence of H. tibetensis from Tibet was obtained, which enriched the mitochondrial genome data of ticks and provided genetic markers to study the population heredity and molecular evolution of the genus Haemaphysalis.

Feature selection for effective prediction of SARS-COV-2 using machine learning

BACKGROUND

With rise in variants of SARS-CoV-2, it is necessary to classify the emerging SARS-CoV-2 for early detection and thereby reduce human transmission. Genomic and proteomic information have less frequently been used for classifying in a machine learning (ML) approach for detection of SARS-CoV-2.

OBJECTIVE

With this aim we used nucleoprotein and viral proteomic evolutionary information of SARS-CoV-2 along with the charge and basicity distribution of amino acids from various strains of SARS-CoV-2 to generate a disease severity model based on ML.

METHODS

All sequence and clinical data were obtained from GISAID. Proteomic level calculations were added to comprise the dataset. The training set was used for feature selection. Select K- Best feature selection method was employed which was cross validated with testing set and performance evaluated. Delong's test was also done. We also employed BIRCH clustering on SARS-CoV-2 for clustering the strains.

RESULTS

Out of six ML models four were successful in training and testing. Extra Trees algorithm generated a micro-averaged F1-score of 74.2% and a weighted averaged area under the receiver operating characteristic curve (AUC-ROC) score of 73.7% with multi-class option. The feature selection set to 5, enhanced the ROC AUC from 73.7 to 76.4%. Accuracy of the selected model of 86.9% was achieved.

CONCLUSION

The unique features identified in the ML approach was able to classify disease severity into classes and had potential for predicting risk in newer variants.

Molecular and phylogenetic characterization of the monkeypox outbreak in the South of Spain

Background and Aim

Until the May 2022 Monkeypox (MPXV) outbreak, which spread rapidly to many non-endemic countries, the virus was considered a viral zoonosis limited to some African countries. The Andalusian circuit of genomic surveillance was rapidly applied to characterize the MPXV outbreak in the South of Spain.

Methods

Whole genome sequencing was used to obtain the genomic profiles of samples collected across the south of Spain, representative of all the provinces of Andalusia. Phylogenetic analysis was used to study the relationship of the isolates and the available sequences of the 2022 outbreak.

Results

Whole genome sequencing of a total of 160 MPXV viruses from the different provinces that reported cases were obtained. Interestingly, we report the sequences of MPXV viruses obtained from two patients who died. While one of the isolates bore no noteworthy mutations that explain a potential heightened virulence, in another patient the second consecutive genome sequence, performed after the administration of tecovirimat, uncovered a mutation within the A0A7H0DN30 gene, known to be a prime target for tecovirimat in its Vaccinia counterpart. In general, a low number of mutations were observed in the sequences reported, which were very similar to the reference of the 2022 outbreak (OX044336), as expected from a DNA virus. The samples likely correspond to several introductions of the circulating MPXV viruses from the last outbreak. The virus sequenced from one of the two patients that died presented a mutation in a gene that bears potential connections to drug resistance. This mutation was absent in the initial sequencing before treatment.

Genetics of H5N1 and H5N8 High-Pathogenicity Avian Influenza Viruses Isolated in Japan in Winter 2021-2022

In winter 2021-2022, H5N1 and H5N8 high-pathogenicity avian influenza (HPAI) viruses (HPAIVs) caused serious outbreaks in Japan: 25 outbreaks of HPAI at poultry farms and 107 cases in wild birds or in the environment. Phylogenetic analyses divided H5 HPAIVs isolated in Japan in the winter of 2021-2022 into three groups-G2a, G2b, and G2d-which were disseminated at different locations and times. Full-genome sequencing analyses of these HPAIVs revealed a strong relationship of multiple genes between Japan and Siberia, suggesting that they arose from reassortment events with avian influenza viruses (AIVs) in Siberia. The results emphasize the complex of dissemination and reassortment events with the movement of migratory birds, and the importance of continual monitoring of AIVs in Japan and Siberia for early alerts to the intrusion of HPAIVs.

Three novel species of Aquapteridospora (Distoseptisporales, Aquapteridosporaceae) from freshwater habitats in Tibetan Plateau, China

During an investigation of lignicolous freshwater fungi in the Tibetan Plateau, three Aquapteridospora taxa were collected from freshwater habitats in Xizang, China. The new species possess polyblastic, sympodial, denticles conidiogenous cells and fusiform, septate, with or without sheath conidial, that fit within the generic concept of Aquapteridospora, and multi-gene phylogeny placed these species within Aquapteridospora. Detailed morphological observations clearly demarcate three of these from extant species and are hence described as new taxa. The multi-gene phylogeny of the combined LSU, TEF1-α, and ITS sequence data to infer phylogenetic relationships and discuss phylogenetic affinities with morphologically similar species. Based on morphological characteristics and phylogenetic analyses, three new species viz. A.linzhiensis, A.yadongensis, and A.submersa are introduced. Details of asexual morphs are described, and justifications for establishing these new species are also provided in this study.

The complete mitochondrial genome of Sarcophaga angarosinica (Diptera: Sarcophagidae)

Sarcophaga (Liosarcophaga) angarosinica (Rohdendorf, 1937) (Diptera: Sarcophagidae) is a species of both medical and ecological significance. In this study, the complete mitochondrial genome (mitogenome) of S. angarosinica was sequenced and characterized. The mitogenome has a total length of 15,215 bp, including 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and an adenine and thymine-rich region. This mitogenome comprises 39.5% adenine, 9.4% guanine, 14.4% cytosine, and 36.8% thymine. Phylogenetic analysis revealed that S. angarosinica is closely related to Sarcophaga similis. This study enriches the genetic data on S. angarosinica and will contribute to establishing the phylogenetic relationships among flesh flies.

Four new complete mitochondrial genomes of Gobioninae fishes (Teleostei: Cyprinidae) and their phylogenetic implications

The subfamily Gobioninae is one of the most diverse fish groups within Cyprinidae. Their taxonomy and phylogenetic relationships are not completely resolved. In this study, the complete mitochondrial genomes (mitogenome) of four Gobioninae species (Microphysogobio elongatus, Microphysogobio chinssuensis, Gobio rivuloides and Rhinogobio nasutus) were sequenced and compared. The mitogenomes of four species ranges from 16603 bp to 16609 bp in length, consisting of 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region. Most PCGs had significant codon usage bias. Except for the tRNASer (GCT), all the nucleotide substitutions of tRNA loops higher than the stems could fold into a stable secondary structure. The nucleotide compositions of Gobioninae mitogenome were biased toward A/T, and NAD4 was subjected to low purification selection and had a faster evolution rate among 13 PCGs. Bayesian inference and maximum likelihood phylogenetic analyses showed the consistent results. The four sequenced species clustered together with their congener species. However, more samples and mitogenome data are needed to untangle the phylogenetic relationships among genera Microphysogobio, Romanogobio, Hugobio, Biwia and Platysmacheilus.

Global characterization of biosynthetic gene clusters in non-model eukaryotes using domain architectures

The majority of pharmaceuticals are derived from natural products, bioactive compounds naturally synthesized by organisms to provide evolutionary advantages. Although the rich evolutionary history of eukaryotic algal species implicates a high potential for natural product-based drug discovery, it remains largely untouched. This study investigates 2762 putative biosynthetic gene clusters (BGCs) from 212 eukaryotic algal genomes. To analyze a vast set of structurally diverse BGCs, we employed comparative analysis based on the vectorization of biosynthetic domains, referred to as biosynthetic domain architecture (BDA). By characterizing core biosynthetic machineries through BDA, we identified key BDAs of modular BGCs in diverse eukaryotes and introduced 16 candidate modular BGCs with similar BDAs to previously validated BGCs. This study provides a global characterization of eukaryotic algal BGCs, offering an alternative to laborious manual curation for BGC prioritization.

Chemolithoautotrophic diazotrophs dominate dark nitrogen fixation in mangrove sediments

Diazotrophic microorganisms regulate marine productivity by alleviating nitrogen limitation. So far chemolithoautotrophic bacteria are widely recognized as the principal diazotrophs in oligotrophic marine and terrestrial ecosystems. However, the contribution of chemolithoautotrophs to nitrogen fixation in organic-rich habitats remains unclear. Here, we utilized metagenomic and metatranscriptomic approaches integrated with cultivation assays to investigate the diversity, distribution, and activity of diazotrophs residing in Zhangzhou mangrove sediments. Physicochemical assays show that the studied mangrove sediments are typical carbon-rich, sulfur-rich, nitrogen-limited, and low-redox marine ecosystems. These sediments host a wide phylogenetic variety of nitrogenase genes, including groups I-III and VII-VIII. Unexpectedly diverse chemolithoautotrophic taxa including Campylobacteria, Gammaproteobacteria, Zetaproteobacteria, and Thermodesulfovibrionia are the predominant and active nitrogen fixers in the 0-18 cm sediment layer. In contrast, the 18-20 cm layer is dominated by active diazotrophs from the chemolithoautotrophic taxa Desulfobacterota and Halobacteriota. Further analysis of MAGs shows that the main chemolithoautotrophs can fix nitrogen by coupling the oxidation of hydrogen, reduced sulfur, and iron, with the reduction of oxygen, nitrate, and sulfur. Culture experiments further demonstrate that members of chemolithoautotrophic Campylobacteria have the nitrogen-fixing capacity driven by hydrogen and sulfur oxidation. Activity measurements confirm that the diazotrophs inhabiting mangrove sediments preferentially drain energy from diverse reduced inorganic compounds other than from organics. Overall, our results suggest that chemolithoautotrophs rather than heterotrophs are dominant nitrogen fixers in mangrove sediments. This study underscores the significance of chemolithoautotrophs in carbon-dominant ecosystems.

An insight on the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina, Asclepiadoideae (Apocynaceae)

Abstract We studied the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina from Asclepiadoideae subfamily; due to their medicinal importance and distribution worldwide their interest became high. In this study we analyzed the complete chloroplast genomes of G. siniacus and D. velutina using Illumina sequencing technology. The sequences were compared with the other species from Apocynaceae family. The complete genome of G. siniacus is 162,570 bp while D. velutina has154, 478 bp in length. Both genomes consist of 119 genes; encode 31 tRNA genes, and eight rRNA genes. Comparative studies of the two genomes showed variations in SSR markers in which G. siniacus possesses 223 while D. velutina has 186. This could be used for barcoding in order to aid in easy identification of the species. Phylogenetic analysis on the other hand reaffirms the tribal position of G. siniacus in Asclepiadeae and D. velutina in Ceropegieae. These findings could be used in subsequent research studies of angiosperms identification, genetic engineering, herb genomics and phylogenomic studies of Apocynaceae family.

Phylogenomics and topological conflicts in the tribe Anthospermeae (Rubiaceae)

Genome skimming (shallow whole-genome sequencing) offers time- and cost-efficient production of large amounts of DNA data that can be used to address unsolved evolutionary questions. Here we address phylogenetic relationships and topological incongruence in the tribe Anthospermeae (Rubiaceae), using phylogenomic data from the mitochondrion, the nuclear ribosomal cistron, and the plastome. All three genomic compartments resolve relationships in the Anthospermeae; the tribe is monophyletic and consists of three major subclades. Carpacoce Sond. is sister to the remaining clade, which comprises an African subclade and a Pacific subclade. Most results, from all three genomic compartments, are statistically well supported; however, not fully consistent. Intergenomic topological incongruence is most notable in the Pacific subclade but present also in the African subclade. Hybridization and introgression followed by organelle capture may explain these conflicts but other processes, such as incomplete lineage sorting (ILS), can yield similar patterns and cannot be ruled out based on the results. Whereas the null hypothesis of congruence among all sequenced loci in the individual genomes could not be rejected for nuclear and mitochondrial data, it was rejected for plastid data. Phylogenetic analyses of three subsets of plastid loci identified using the hierarchical likelihood ratio test demonstrated statistically supported intragenomic topological incongruence. Given that plastid genes are thought to be fully linked, this result is surprising and may suggest modeling or sampling error. However, biological processes such as biparental inheritance and inter-plastome recombination have been reported and may be responsible for the observed intragenomic incongruence. Mitochondrial insertions into the plastome are rarely documented in angiosperms. Our results indicate that a mitochondrial insertion event in the plastid trnS GGA - rps4 IGS region occurred in the common ancestor of the Pacific clade of Anthospermeae. Exclusion/inclusion of this locus in phylogenetic analyses had a strong impact on topological results in the Pacific clade.

Dietary- and host-derived metabolites are used by diverse gut bacteria for anaerobic respiration

Respiratory reductases enable microorganisms to use molecules present in anaerobic ecosystems as energy-generating respiratory electron acceptors. Here we identify three taxonomically distinct families of human gut bacteria (Burkholderiaceae, Eggerthellaceae and Erysipelotrichaceae) that encode large arsenals of tens to hundreds of respiratory-like reductases per genome. Screening species from each family (Sutterella wadsworthensis, Eggerthella lenta and Holdemania filiformis), we discover 22 metabolites used as respiratory electron acceptors in a species-specific manner. Identified reactions transform multiple classes of dietary- and host-derived metabolites, including bioactive molecules resveratrol and itaconate. Products of identified respiratory metabolisms highlight poorly characterized compounds, such as the itaconate-derived 2-methylsuccinate. Reductase substrate profiling defines enzyme-substrate pairs and reveals a complex picture of reductase evolution, providing evidence that reductases with specificities for related cinnamate substrates independently emerged at least four times. These studies thus establish an exceptionally versatile form of anaerobic respiration that directly links microbial energy metabolism to the gut metabolome.

Complete mitochondrial genome and phylogenetic analysis of the marine microalga Symbiochlorum hainanensis (Ulvophyceae, Chlorophyta)

Symbiochlorum hainanensis Gong et al. (2018). is a unicellular green alga belonging to Ulvophyceae, Chlorophyta, and considered as an important species in coral-algae symbiont areas. In this study, we revealed firstly the mitochondrial genome sequence of the S. hainanensis. This mitochondrial genome was a circular DNA molecule of 59,508 bp, including 24 transfer RNA genes, 3 ribosomal RNA genes, and 31 protein-coding genes. The GC content of the genome was 35.4%. The phylogenetic tree suggested that S. hainandiae was a sister to the OUU clade within the class Ulvophyceae. The mitochondrial genome structure and gene content of S. hainanensis supported that S. hainanensis was a new unidentified green alga in Ulvophyceae.

Safety assessment, whole genome sequence, and metabolome analysis of Streptococcus thermophilus CICC 20372 for bone cement fermentation

Bone is a kind of meat processing by-product with high nutritional value but low in calorie, which is a typical food in China and parts of East Asian countries. Microbial fermentation by lactic acid bacteria showed remarkable advantages to increase the absorption of nutrients from bone cement by human body. Streptococcus thermophilus CICC 20372 is proven to be a good starter for bone cement fermentation. No genes encoding virulence traits or virulence factors were found in the genome of S. thermophilus CICC 20372 by a thorough genomic analysis. Its notable absence of antibiotic resistance further solidifies the safety. Furthermore, the genomic analysis identified four types of gene clusters responsible for the synthesis of antimicrobial metabolites. A comparative metabolomic analysis was performed by cultivating the strain in bone cement at 37 °C for 72 h, with the culture in de Man, Rogosa, and Sharpe (MRS) medium as control. Metabolome analysis results highlighted the upregulation of pathways involved in 2-oxocarboxylic acid metabolism, ATP-binding cassette (ABC) transporters, amino acid synthesis, and nucleotide metabolism during bone cement fermentation. S. thermophilus CICC 20372 produces several metabolites with health-promoting function during bone cement fermentation, including indole-3-lactic acid, which is demonstrated ameliorative effects on intestinal inflammation, tumor growth, and gut dysbiosis. In addition, lots of nucleotide and organic acids were accumulated at higher levels, which enriched the fermented bone cement with a variety of nutrients. Collectively, these features endow S. thermophilus CICC 20372 a great potential strain for bone food processing.

Characterization of the Complete Mitochondrial Genome and Phylogenetic Analyses of Eurytrema coelomaticum (Trematoda: Dicrocoeliidae)

Eurytrema coelomaticum, a pancreatic fluke, is recognized as a causative agent of substantial economic losses in ruminants. This infection, commonly referred to as eurytrematosis, is a significant concern due to its detrimental impact on livestock production. However, there is a paucity of knowledge regarding the mitochondrial genome of E. coelomaticum. In this study, we performed the initial sequencing of the complete mitochondrial genome of E. coelomaticum. Our findings unveiled that the mitochondrial genome of E. coelomaticum spans a length of 15,831 bp and consists of 12 protein-coding genes, 22 tRNA genes, two rRNA genes, and two noncoding regions. The A+T content constituted 62.49% of the genome. Moreover, all 12 protein-coding genes of E. coelomaticum exhibit the same arrangement as those of E. pancreaticum and other published species belonging to the family Dicrocoeliidae. The presence of a short string of additional amino acids (approximately 20~23 aa) at the N-terminal of the cox1 protein in both E. coelomaticum and E. pancreaticum mitochondrial genomes has contributed to the elongation of the cox1 gene in genus Eurytrema, surpassing that of all previously sequenced Dicrocoeliidae. The phylogenetic analysis displayed a close relationship between E. coelomaticum and E. pancreaticum, along with a genus-level association between Eurytrema and Lyperosomum. These findings underscore the importance of mitochondrial genomic data for comparative studies of Dicrocoeliidae and even Digenea, offering valuable DNA markers for future investigations in the systematic, epidemiological, and population genetic studies of this parasite and other digenean trematodes.

Integrative taxonomy of Serrasentis gibsoni n. sp. (Acanthocephala: Isthmosacanthidae) from flatfishes in the Gulf of Mexico

The Isthmosacanthidae acanthocephalan species of the genus Serrasentis are parasites of marine teleosts and an elasmobranch. In this study, Serrasentis gibsoni n. sp. is described from the intestines of four flatfish species (Paralichthyidae), namely Ancyclopsetta quadrocellata, Cyclopsetta chittendeni, Syacium gunteri, and S. papillosum from 10 oceanic sites in the Gulf of Mexico (GoM). Twenty sequences of the 'barcoding' region of cytochrome C oxidase subunit I gene were obtained from 20 adults of Serrasentis gibsoni n. sp. Additionally, five sequences of the barcoding region were obtained from five adults of rhadinorhynchid Gorgorhynchus lepidus from C. chittendeni, S. papillosum and one species of Haemulidae, Haemulom aurolineatum, from five oceanic sites from the GoM. Two phylogenetic approaches were followed: Bayesian inference and maximum likelihood. In both phylogenetic reconstructions, the sequences of Serrasentis gibsoni n. sp. were recovered as a monophyletic group within the genus Serrasentis and placed as a sister group to G. lepidus. However, due to the lack of molecular data for species of the Isthmosacanthidae and Rhadinorhynchidea, these phylogenetic inferences must be taken with caution. Serrasentis gibsoni n. sp. is the first species of Serrasentis described from Paralichthyidae flatfish species from marine waters of the Americas and from the GoM. Based on the barcoding data set analyzed, Serrasentis gibsoni n. sp. appears to have high intraspecific genetic variation; thus, it is necessary to continue exploring the genetic diversity of this species to infer its intraspecific evolutionary patterns.

Characterization of the chloroplast genome of the marine microalga Tetraselmis marina (Cienkowski) R.E.Norris, Hori & Chihara 1980

Tetraselmis marina (Cienkowski) R.E.Norris, Hori & Chihara 1980, a costal green microalga, is considered as a promising animal feed in aquaculture due to the high content of fatty acids and carotenoid. Furthermore, T. marina plays important roles in bioremediation. In this study, we assembled the complete chloroplast genome of T. marina. Results showed that the full length of the complete chloroplast genome was 96,151 bp, containing a large single-copy region of 62,574 bp, a small single-copy region of 1261 bp, and a pair of inverted repeat regions of 16,158 bp. The GC content of the genome was 36.6%. A total of 125 genes were annotated, including 81 protein coding genes, 38 tRNA genes, and six rRNA genes. Phylogenetic analysis based on 22 chloroplast genomes suggested that T. marina was closely related to Tetraselmis sp. CCMP 881.

Leveraging genomic redundancy to improve inference and alignment of orthologous proteins

Identifying protein sequences with common ancestry is a core task in bioinformatics and evolutionary biology. However, methods for inferring and aligning such sequences in annotated genomes have not kept pace with the increasing scale and complexity of the available data. Thus, in this work, we implemented several improvements to the traditional methodology that more fully leverage the redundancy of closely related genomes and the organization of their annotations. Two highlights include the application of the more flexible k-clique percolation algorithm for identifying clusters of orthologous proteins and the development of a novel technique for removing poorly supported regions of alignments with a phylogenetic hidden Markov model (phylo-HMM). In making the latter, we wrote a fully documented Python package Homomorph that implements standard HMM algorithms and created a set of tutorials to promote its use by a wide audience. We applied the resulting pipeline to a set of 33 annotated Drosophila genomes, generating 22,813 orthologous groups and 8,566 high-quality alignments.

Phosphate Limitation Responses in Marine Green Algae Are Linked to Reprogramming of the tRNA Epitranscriptome and Codon Usage Bias

Marine algae are central to global carbon fixation, and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change; however, there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA modification enzymes and biased codon usage of transcripts more abundant during phosphate-limiting versus phosphate-replete conditions, consistent with the role of transfer RNA modifications in regulating codon recognition. To associate the observed shift in the expression of the transfer RNA modification enzyme complement with the transfer RNAs encoded by M. commoda, we also determined the transfer RNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing phosphate limitation and photosystem-associated proteins that have ill-characterized putative functions in "light stress." The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in transfer RNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation that extends beyond the catalog of up- or downregulated protein-encoding genes to the cell biological controls that underpin acclimation to changing environmental conditions.

Single-cell genomics reveals new rozellid lineages and supports their sister relationship to Microsporidia

The phylum Rozellomycota has been proposed for a group of early-branching holomycotan lineages representing obligate parasites and hyperparasites of zoosporic fungi, oomycotes or phytoplankton. Given their predominantly intracellular lifestyle, rozellids are typically known from environmental ribosomal DNA data, except for the well-studied Rozella species. To date, the phylogenetic relationship between rozellids and microsporidians (Microsporidia) is not fully understood and most reliable hypotheses are based on phylogenomic analyses that incorporate the only publicly available rozellid genome of Rozella allomycis. Here, we provide genomic data of three new rozellid lineages obtained by single-cell sequencing from environmental samples and show with a phylogenomic approach that rozellids form a monophyletic group that is sister to microsporidians, corroborating the previously proposed phylum Rozellomycota. Whereas no mitochondrial genes coding for the respiratory Complex I could be found, we discovered a gene coding for a nucleotide phosphate transporter in one of the three draft genomes. The scattered absence of Complex I genes and scattered presence of nucleotide transporter genes across diverse microsporidian and rozellid lineages suggest that these adaptations to a parasitic lifestyle, which reduce the parasite's capability to synthesize ATP but enables it to steal ATP from its host, evolved independently in microsporidians and rozellids.

Geographical Distribution and Multimethod Species Identification of Forensically Important Necrophagous Flies on Hainan Island

Forensic entomology offers unique advantages for the minimum postmortem interval (PMImin) estimation of decomposed corpses in forensic investigations. Accurate species identification and up-to-date locality information are essential. Hainan Island has a tropical rainforest climate and a vast territory. In this study, the community structure of necrophagous flies on Hainan Island was investigated in detail according to geographical environment. The results showed that the dominant species included C. megacephala, S. peregrina, C. rufifacies, S. misera, H. ligurriens, S. sericea, S. cinerea, S. dux, C. pinguis, and M. domestica. Furthermore, C. rufifacies and C. villeneuvi were found only in the high-altitude areas of Wuzhi Mountain, while S. cinerea was distributed only in coastal areas; the latter is a representative species of Hainan Island and has not been reported before. Furthermore, a GenBank database of forensically important flies was established, whilst a high-resolution melt (HRM) curve analysis was applied to identify the common species of Hainan Island for the first time. This study enriches the database of forensically important flies in tropical rainforest regions.

Systematic assessment of nymphs of Diesingiinae (Pentastomida: Sebekidae) parasitizing Palloceros harpagos (Cyprinodontiformes: Poeciliidae) from Parque Nacional Iguazú, Argentina: filling gaps in an incomplete phylogenetic framework

During a search for parasites in fishes from Iguazu River, Argentina, a nymph of pentastomid was found encysted on mesenteries of Phalloceros harpagos. The 28S rDNA and COI mt-DNA were used to compare with the sequences deposited in the GenBank. Pentastomid nymphs belong to the subfamily Diesingiinae (Sebekidae) for having chloride cell pores distributed in a single row per annulus; also, the hooks and rows of spines of our material differ to other genera. Present specimens are most likely Diesingia sp., having 63-74 annuli, a U shaped oral cadre with fibers closing anteriorly and a peg like extension of the oral cadre. The 28S rDNA analysis places our samples into a sister group of Alofia sp., but the COI mt-DNA situate them on the base of the clade. In conclusion, our pentastomid positively belongs to Diesingia sp., however, indisputable assignation requires a reliable description of the nymph, or the availability of sequences linking nymphs and adults, which even could provide evidence of a new hitherto undescribed genus. Aditionally, the systematic position of Sebekia oxycephala previously described by P. harpagos is discussed on the basis of our results, allowing us to suggest a re-assignation of such specimens to the genus Diesingia.

The complete plastome of Campsis radicans (L.) Bureau 1864 and its phylogenetic analysis

Campsis radicans (L.) Bureau 1864, a species of Bignoniaceae, has a widespread paleotropical distribution and is utilized for horticultural and traditional Chinese medicinal purposes. Despite the plant's significance, its genetic diversity must be better understood. In this study, we have successfully assembled and characterized the complete plastome of C. radicans, marking a significant advancement toward comprehending its genetic composition. The plastome is 153,630 bp long and harbors 130 genes, including 86 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Our phylogenomic analysis of the representative species of Bignoniaceae indicated that C. radicans formed a monophyletic sister clade of Campsis with C. grandiflora. These findings are crucial for conserving and utilizing this important plant species. They also highlight the potential for future research into the evolution and preservation of C. radicans, which could be advantageous in pharmaceutical applications.

Comparative Population Transcriptomics Provide New Insight into the Evolutionary History and Adaptive Potential of World Ocean Krill

Genetic variation is instrumental for adaptation to changing environments but it is unclear how it is structured and contributes to adaptation in pelagic species lacking clear barriers to gene flow. Here, we applied comparative genomics to extensive transcriptome datasets from 20 krill species collected across the Atlantic, Indian, Pacific, and Southern Oceans. We compared genetic variation both within and between species to elucidate their evolutionary history and genomic bases of adaptation. We resolved phylogenetic interrelationships and uncovered genomic evidence to elevate the cryptic Euphausia similis var. armata into species. Levels of genetic variation and rates of adaptive protein evolution vary widely. Species endemic to the cold Southern Ocean, such as the Antarctic krill Euphausia superba, showed less genetic variation and lower evolutionary rates than other species. This could suggest a low adaptive potential to rapid climate change. We uncovered hundreds of candidate genes with signatures of adaptive evolution among Antarctic Euphausia but did not observe strong evidence of adaptive convergence with the predominantly Arctic Thysanoessa. We instead identified candidates for cold-adaptation that have also been detected in Antarctic fish, including genes that govern thermal reception such as TrpA1. Our results suggest parallel genetic responses to similar selection pressures across Antarctic taxa and provide new insights into the adaptive potential of important zooplankton already affected by climate change.

A new cryptic species of Imparfinis (Siluriformes: Heptapteridae) from the Orinoco River basin, revealed by an iterative approach

An iterative analysis of Imparfinis, combining phylogenetic analysis based on cytochrome oxidase gene and multivariate morphometrics, revealed a new cryptic species from the Andean tributaries of the Orinoco River basin, which is described here. The new species is sister to a clade constituted by Imparfinis hasemani and Imparfinis pijpersi, both from the river basins of the Guiana Shield, being also the most geographically proximate species. Nonetheless, the new species is most similar in general appearance to Imparfinis guttatus from the Madeira and Paraguay River drainages, being almost undistinguishable by conventional characters of external morphology, differing only by morphometric attributes overall. The new species can be distinguished from the remaining congeners by a unique combination of characters, including lower lobe of caudal fin darker than upper lobe, maxillary barbel reaching or surpassing pelvic-fin insertion, 12-15 gill rakers on first gill arch, 40-42 total vertebrae and 9-10 ribs. The new species constitutes the only representative from the Orinoco River basin belonging to Imparfinis sensu stricto.

Anelloviridae taxonomy update 2023

The family Anelloviridae comprises negative single-stranded circular DNA viruses. Within this family, there are 30 established genera. Anelloviruses in the genus Gyrovirus have been identified infecting various avian species, whereas those in the remaining 29 genera have been found primarily infecting various mammal species. We renamed the 146 anellovirus species with binomial species names, as required by the International Committee on Taxonomy of Viruses (ICTV) using a "genus + freeform epithet" format.

Identification and characterization of novel alphacoronaviruses in Tadarida brasiliensis (Chiroptera, Molossidae) from Argentina: insights into recombination as a mechanism favoring bat coronavirus cross-species transmission

Bats are reservoirs of various coronaviruses that can jump between bat species or other mammalian hosts, including humans. This article explores coronavirus infection in three bat species (Tadarida brasiliensis, Eumops bonariensis, and Molossus molossus) of the family Molossidae from Argentina using whole viral metagenome analysis. Fecal samples of 47 bats from three semiurban or highly urbanized areas of the province of Santa Fe were investigated. After viral particle enrichment, total RNA was sequenced using the Illumina NextSeq 550 instrument; the reads were assembled into contigs and taxonomically and phylogenetically analyzed. Three novel complete Alphacoronavirus (AlphaCoV) genomes (Tb1-3) and two partial sequences were identified in T. brasiliensis (Tb4-5), and an additional four partial sequences were identified in M. molossus (Mm1-4). Phylogenomic analysis showed that the novel AlphaCoV clustered in two different lineages distinct from the 15 officially recognized AlphaCoV subgenera. Tb2 and Tb3 isolates appeared to be variants of the same virus, probably involved in a persistent infectious cycle within the T. brasiliensis colony. Using recombination analysis, we detected a statistically significant event in Spike gene, which was reinforced by phylogenetic tree incongruence analysis, involving novel Tb1 and AlphaCoVs identified in Eptesicus fuscus (family Vespertilionidae) from the U.S. The putative recombinant region is in the S1 subdomain of the Spike gene, encompassing the potential receptor-binding domain of AlphaCoVs. This study reports the first AlphaCoV genomes in molossids from the Americas and provides new insights into recombination as an important mode of evolution of coronaviruses involved in cross-species transmission. IMPORTANCE This study generated three novel complete AlphaCoV genomes (Tb1, Tb2, and Tb3 isolates) identified in individuals of Tadarida brasiliensis from Argentina, which showed two different evolutionary patterns and are the first to be reported in the family Molossidae in the Americas. The novel Tb1 isolate was found to be involved in a putative recombination event with alphacoronaviruses identified in bats of the genus Eptesicus from the U.S., whereas isolates Tb2 and Tb3 were found in different collection seasons and might be involved in persistent viral infections in the bat colony. These findings contribute to our knowledge of the global diversity of bat coronaviruses in poorly studied species and highlight the different evolutionary aspects of AlphaCoVs circulating in bat populations in Argentina.

Three coralloid species of the genus Trechispora (Trechisporales, Basidiomycota) in China: two newly discovered taxa and one reported for the first time

Two new species of Trechispora indigenous to southern China, T.laxa and T.tongdaoensis, are described and illustrated, and the first record of T.khokpasiensis in China is reported. Molecular phylogenetic analyses of the concatenated nuclear rDNA ITS1-5.8S-ITS2 and nuclear large subunit sequences supported the inclusion of the three species within the Trechispora clade, together with species formerly classified in Scytinopogon. The new species are similar in micromorphology to species of Trechispora (as traditionally circumscribed) but are distinguished by having coralloid basidiomata. A key to the known coralloid Trechispora species in China is provided.

FREEDA: An automated computational pipeline guides experimental testing of protein innovation

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that lead to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA that provides a simple graphical user interface requiring only a gene name; integrates widely used molecular evolution tools to detect positive selection in rodents, primates, carnivores, birds, and flies; and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 centromere proteins, we find statistical evidence of positive selection within loops and turns of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of mouse CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

Isolation, characterization, and genome analysis of a broad host range Salmonella phage vB_SenS_TUMS_E4: a candidate bacteriophage for biocontrol

Salmonella enteritidis is one of the most important foodborne pathogens that cause numerous outbreaks worldwide. Some strains of Salmonella have become progressively resistant to antibiotics, so they could represent a critical threat to public health and have led to the use of alternative therapeutic approaches like phage therapy. In this study, a lytic phage, vB_SenS_TUMS_E4 (E4), was isolated from poultry effluent and characterized to evaluate its potential and efficacy for bio-controlling S. enteritidis in foods. Transmission electron microscopy revealed that E4 has a siphovirus morphotype, with an isometric head and non-contractile tail. Determining the host range showed that this phage can effectively infect different motile as well as non-motile Salmonella enterica serovars. The biological characteristics of E4 showed that it has a short latent period of about 15 min and a large burst size of 287 PFU/cell, and is also significantly stable in a broad range of pHs and temperatures. The E4 whole genome contains 43,018 bp and encodes 60 coding sequences (CDSs) but no tRNA genes. Bioinformatics analysis revealed that the genome of E4 lacks any genes related to lysogeny behavior, antibiotic resistance, toxins, or virulence factors. The efficacy of phage E4 as a bio-control agent was assessed in various foodstuffs inoculated with S. enteritidis at 4°C and 25°C, and the resulting data indicated that it could eradicate S. enteritidis after a very short time of 15 min. The findings of the present study showed that E4 is a hopeful candidate as a bio-control agent against S. enteritidis and has the potential to be used in various foodstuffs.

A contribution to the knowledge of the genus Infundibulicybe (Tricholomatineae, Agaricales) in China: Two new species and five redescribed taxa

In China, species of the genus Infundibulicybe are often confused with other taxa and misdetermined. Here we describe two newly discovered species, I. bispora and I. ellipsospora, and redescribe five known taxa of this genus present in China. These identifications are supported by both morphological and DNA-based phylogenetic evidence. Additionally, a key to all known species of Infundibulicybe is provided.

Long noncoding RNAs contribute to DNA damage resistance in Arabidopsis thaliana

Efficient repair of DNA lesions is essential for the faithful transmission of genetic information between somatic cells and for genome integrity across generations. Plants have multiple, partially redundant, and overlapping DNA repair pathways, probably due to the less constricted germline and the inevitable exposure to light including higher energy wavelengths. Many proteins involved in DNA repair and their mode of actions are well described. In contrast, a role for DNA damage-associated RNA components, evident from many other organisms, is less well understood. Here, we have challenged young Arabidopsis thaliana plants with two different types of genotoxic stress and performed de novo assembly and transcriptome analysis. We identified three long noncoding RNAs (lncRNAs) that are lowly or not expressed under regular conditions but up-regulated or induced by DNA damage. We generated CRISPR/Cas deletion mutants and found that the absence of the lncRNAs impairs the recovery capacity of the plants from genotoxic stress. The genetic loci are highly conserved among world-wide distributed Arabidopsis accessions and within related species in the Brassicaceae group. Together, these results suggest that the lncRNAs have a conserved function in connection with DNA damage and provide a basis for mechanistic analysis of their role.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

The fragmented mitochondrial genomes of two Linognathus lice reveal active minichromosomal recombination and recombination hotspots

Evidence for recombination between mitochondrial (mt) minichromosomes has been reported in sucking lice, but it is still not clear how frequent mt minichromosomal recombination occurs. We sequenced the mt genomes of the cattle louse Linognathus vituli and the goat louse L. africanus. Both Linognathus species have 10 mt minichromosomes, and seven of them have the same gene content and gene arrangement. Comparison of mt karyotypes revealed numerous inter-minichromosomal recombination events in the evolution of Linognathus species. Minichromosome merger, gene duplication and gene translocation occurred in the lineage leading to Linognathus lice. After the divergence of L. vituli and L. africanus, duplication, degeneration, deletion and translocation of genes also occurred independently in each species. Most of the recombination events in the Linognathus species occurred upstream of either cox3 or nad2, indicating these two locations were hotspots for inter-minichromosomal recombination. Our results provide an important perspective on mt genome evolution in metazoans.

Genomic insights into biased allele loss and increased gene numbers after genome duplication in autotetraploid Cyclocarya paliurus

BACKGROUND

Autopolyploidy is a valuable model for studying whole-genome duplication (WGD) without hybridization, yet little is known about the genomic structural and functional changes that occur in autopolyploids after WGD. Cyclocarya paliurus (Juglandaceae) is a natural diploid-autotetraploid species. We generated an allele-aware autotetraploid genome, a chimeric chromosome-level diploid genome, and whole-genome resequencing data for 106 autotetraploid individuals at an average depth of 60 × per individual, along with 12 diploid individuals at an average depth of 90 × per individual.

RESULTS

Autotetraploid C. paliurus had 64 chromosomes clustered into 16 homologous groups, and the majority of homologous chromosomes demonstrated similar chromosome length, gene numbers, and expression. The regions of synteny, structural variation and nonalignment to the diploid genome accounted for 81.3%, 8.8% and 9.9% of the autotetraploid genome, respectively. Our analyses identified 20,626 genes (69.18%) with four alleles and 9191 genes (30.82%) with one, two, or three alleles, suggesting post-polyploid allelic loss. Genes with allelic loss were found to occur more often in proximity to or within structural variations and exhibited a marked overlap with transposable elements. Additionally, such genes showed a reduced tendency to interact with other genes. We also found 102 genes with more than four copies in the autotetraploid genome, and their expression levels were significantly higher than their diploid counterparts. These genes were enriched in enzymes involved in stress response and plant defense, potentially contributing to the evolutionary success of autotetraploids. Our population genomic analyses suggested a single origin of autotetraploids and recent divergence (~ 0.57 Mya) from diploids, with minimal interploidy admixture.

CONCLUSIONS

Our results indicate the potential for genomic and functional reorganization, which may contribute to evolutionary success in autotetraploid C. paliurus.

Phage tailspike modularity and horizontal gene transfer reveals specificity towards E. coli O-antigen serogroups

BACKGROUND

The interaction between bacteriophages and their hosts is intricate and highly specific. Receptor-binding proteins (RBPs) of phages such as tail fibers and tailspikes initiate the infection process. These RBPs bind to diverse outer membrane structures, including the O-antigen, which is a serogroup-specific sugar-based component of the outer lipopolysaccharide layer of Gram-negative bacteria. Among the most virulent Escherichia coli strains is the Shiga toxin-producing E. coli (STEC) pathotype dominated by a subset of O-antigen serogroups.

METHODS

Extensive phylogenetic and structural analyses were used to identify and validate specificity correlations between phage RBP subtypes and STEC O-antigen serogroups, relying on the principle of horizontal gene transfer as main driver for RBP evolution.

RESULTS

We identified O-antigen specific RBP subtypes for seven out of nine most prevalent STEC serogroups (O26, O45, O103, O104, O111, O145 and O157) and seven additional E. coli serogroups (O2, O8, O16, O18, 4s/O22, O77 and O78). Eight phage genera (Gamaleya-, Justusliebig-, Kaguna-, Kayfuna-, Kutter-, Lederberg-, Nouzilly- and Uetakeviruses) emerged for their high proportion of serogroup-specific RBPs. Additionally, we reveal sequence motifs in the RBP region, potentially serving as recombination hotspots between lytic phages.

CONCLUSION

The results contribute to a better understanding of mosaicism of phage RBPs, but also demonstrate a method to identify and validate new RBP subtypes for current and future emerging serogroups.

Trichinella T9 in wild bears in Japan: Prevalence, species/genotype identification, and public health implications

In Japan, the recent series of sporadic outbreaks of human trichinellosis caused by Trichinella (Nematoda: Trichocephalida) has occurred owing to the consumption of raw or insufficiently cooked meat from wild bears. However, the infection status and molecular characteristics of Trichinella larvae in Japanese wild bears remain poorly understood. This study investigated the prevalence of Trichinella spp. in brown bears (Ursus arctos) from Hokkaido, and Japanese black bears (Ursus thibetanus japonicus) from three prefectures (Aomori, Akita, and Iwate) in northern Japan, between April 2019 and August 2022. Trichinella larvae were detected in 2.5% (6/236) of the brown bears and 0.9% (1/117) of the Japanese black bears. Sequence analysis using two genetic loci, the internal transcribed spacer region of nuclear ribosomal DNA and the mitochondrial cytochrome c oxidase subunit I gene, revealed that the larvae collected from the seven infected bears were identical to one of the two haplotypes of Trichinella T9. The prevalence of Trichinella T9 is low but is maintained in bears in the Hokkaido and Iwate prefectures suggesting that undercooked meat from these animals could cause human infection. Thus, continued health education campaigns are needed to raise awareness of the potential risk of trichinellosis among hunters, meat suppliers, consumers, and local governmental health agencies.

Expansion of Betatorquevirus and/or Gammatorquevirus in Patients with Severe Clinical Outcomes of the Liver Diseases

Anellovirus (AV) is a ubiquitous virus in the human population. Individuals can be infected with multiple AV genera and species to form a heterogeneous repertoire, termed the anellome. Using advanced methods, we examined the anellomes from 12 paired serum and liver samples, as well as 2701 subjects with different clinical diagnoses. Overall, anellomes are remarkably individualized, with significant among-group differences (Kruskal-Wallis test p = 6.6 × 10-162 for richness and p = 7.48 × 10-162 for Shannon entropy). High dissimilarity scores (beta diversity) were observed between patient groups, except for paired serum and liver samples. At the population level, the relative abundance of combinational AV genus Betatorquevirus (torque teno mini viruses, TTMV), and Gammatorquevirus (torque teno midi viruses, TTMDV) exhibited an exponential distribution with a low bound point at 32%. Defined by this value, the AV TTMV/TTMDV-expanded anellome was significantly enriched among patients with acute liver failure (31.7%) and liver transplantation (40.7%), compared with other patient groups (χ2 test: p = 4.1 × 10-8-3.2 × 10-3). Therefore, anellome heterogeneity may be predictive of clinical outcomes in certain diseases, such as liver disease. The consistency of anellome between paired serum and liver samples indicates that a liquid biopsy approach would be suitable for longitudinal studies to clarify the causality of the AV TTMV/TTMDV-expanded anellome in the outcomes of liver disease.

Diversity, Relationship, and Distribution of Virophages and Large Algal Viruses in Global Ocean Viromes

Virophages are a group of small double-stranded DNA viruses that replicate and proliferate with the help of the viral factory of large host viruses. They are widely distributed in aquatic environments but are more abundant in freshwater ecosystems. Here, we mined the Global Ocean Viromes 2.0 (GOV 2.0) dataset for the diversity, distribution, and association of virophages and their potential host large viruses in marine environments. We identified 94 virophage sequences (>5 kbp in length), of which eight were complete genomes. The MCP phylogenetic tree showed that the GOV virophages were widely distributed on the global virophage tree but relatively clustered on three major branches. The gene-sharing network divided GOV virophages into 21 outliers, 2 overlaps, and 14 viral clusters, of which 4 consisted of only the GOV virophages. We also identified 45 large virus sequences, 8 of which were >100 kbp in length and possibly involved in cell-virus-virophage (C-V-v) trisome relationships. The potential eukaryotic hosts of these eight large viruses and the eight virophages with their complete genomes identified are likely to be algae, based on comparative genomic analysis. Both homologous gene and codon usage analyses support a possible interaction between a virophage (GOVv18) and a large algal virus (GOVLV1). These results indicate that diverse and novel virophages and large viruses are widespread in global marine environments, suggesting their important roles and the presence of complicated unknown C-V-v relationships in marine ecosystems.

Insights into the Genomics of Clownfish Adaptive Radiation: The Genomic Substrate of the Diversification

Clownfishes are an iconic group of coral reef fishes that evolved a mutualistic interaction with sea anemones, which triggered the rapid diversification of the group. Following the emergence of this mutualism, clownfishes diversified into different ecological niches and developed convergent phenotypes associated with their host use. The genetic basis of the initial acquisition of the mutualism with host anemones has been described, but the genomic architecture underlying clownfish diversification once the mutualism was established and the extent to which clownfish phenotypic convergence originated through shared genetic mechanisms are still unknown. Here, we investigated these questions by performing comparative genomic analyses on the available genomic data of five pairs of closely related but ecologically divergent clownfish species. We found that clownfish diversification was characterized by bursts of transposable elements, an overall accelerated coding evolution, incomplete lineage sorting, and ancestral hybridization events. Additionally, we detected a signature of positive selection in 5.4% of the clownfish genes. Among them, five presented functions associated with social behavior and ecology, and they represent candidate genes involved in the evolution of the size-based hierarchical social structure so particular to clownfishes. Finally, we found genes with patterns of either relaxation or intensification of purifying selection and signals of positive selection linked with clownfish ecological divergence, suggesting some level of parallel evolution during the diversification of the group. Altogether, this work provides the first insights into the genomic substrate of clownfish adaptive radiation and integrates the growing collection of studies investigating the genomic mechanisms governing species diversification.

Paenibacillus plantiphilus sp. nov. from the plant environment of Zea mays

A Gram-strain positive, aerobic, endospore-forming bacterial strain (JJ-246^T) was isolated from the rhizosphere of Zea mays. The 16S rRNA gene sequence similarity comparisons showed a most closely relationship to Paenibacillus oenotherae DT7-4^T (98.4%) and Paenibacillus xanthinolyticus 11N27^T (98.0%). The pairwise average nucleotide identity and digital DNA-DNA hybridisation values of the JJ-246^T genome assembly against publicly available Paenibacillus type strain genomes were below 82% and 33%, respectively. The draft genome of JJ-246^T shared many putative plant-beneficial functions contributing (PBFC) genes, related to plant root colonisation, oxidative stress protection, degradation of aromatic compounds, plant growth-promoting traits, disease resistance, drug and heavy metal resistance, and nutrient acquisition. The quinone system of strain JJ-246^T, the polar lipid profile and the major fatty acids were congruent with those reported for members of the genus Paenibacillus. JJ-246^T was shown to represent a novel species of the genus Paenibacillus, for which the name Paenibacillus plantiphilus sp. nov. is proposed, with JJ-246^T (= LMG 32093^T = CCM 9089^T = CIP 111893^T) as the type strain.

Novel Gene Rearrangement Pattern in Pachycrepoideus vindemmiae Mitochondrial Genome: New Gene Order in Pteromalidae (Hymenoptera: Chalcidoidea)

The mitochondrial genomes of Muscidifurax similadanacus, M. sinesensilla, Nasonia vitripennis, and Pachycrepoideus vindemmiae were sequenced to better understand the structural evolution of Pteromalidae mitogenomes. These newly sequenced mitogenomes all contained 37 genes. Nucleotide composition was AT-biased and the majority of the protein-coding genes exhibited a negative AT skew. All 13 protein-coding genes (PCGs) initiated with the standard start codon of ATN, excepted for nad1 of N. vitripennis, which started with TTG, and terminated with a typical stop codon TAA/TAG or an incomplete stop codon T. All transfer RNA (tRNA) genes were predicted to fold into the typical clover-leaf secondary structures, except for trnS1, which lacks the DHU arm in all species. In P. vindemmiae, trnR and trnQ lack the DHU arm and TΨC arm, respectively. Although most genes evolved under a strong purifying selection, the Ka/Ks value of the atp8 gene of P. vindemmiae was greater than 1, indicating putative positive selection. A novel transposition of trnR in P. vindemmiae was revealed, which was the first of this kind to be reported in Pteromalidae. Two kinds of datasets (PCG12 and AA) and two inference methods (maximum likelihood and Bayesian inference) were used to reconstruct a phylogenetic hypothesis for the newly sequenced mitogenomes of Pteromalidae and those deposited in GenBank. The topologies obtained recovered the monophyly of the three subfamilies included. Pachyneurinae and Pteromalinae were recovered as sister families, and both appeared sister to Sycophaginae. The pairwise breakpoint distances of mitogenome rearrangements were estimated to infer phylogeny among pteromalid species. The topology obtained was not totally congruent with those reconstructed using the ML and BI methods.

Taxonomy and Phylogenetic Relationships of Clavulinopsis (Clavariaceae, Agaricales): Description of Six New Species and One Newly Recorded Species from China

Specimens of Clavulinopsis (Clavariaceae, Agaricales) collected in China were studied using morphological and molecular methods. Six species-C. aspersa, C. bicolor, C. bispora, C. erubescens, C. incarnata, and C. tropicalis-are described as new to science, and C. trigonospora is a newly recorded species in China. Phylogenetic analysis was conducted based on a combined dataset of internal transcribed spacer and nuclear ribosomal RNA large subunit sequences. The phylogenetic reconstruction revealed that the six new species each formed an independent lineage, and the samples of C. trigonospora from China were nested with accessions of C. trigonospora collected from Italy. The morphology of the seven Chinese species is described in detail, and is illustrated with line drawings and photographs. A key to the known Clavulinopsis species in China is provided.