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

Impact factors and genetic characteristics of head lice infestation in schoolchildren: a cross-sectional study from 2018 to 2023 in central China

BACKGROUND

Head lice (Pediculus humanus capitis) infestation is a worldwide public health concern, especially in school-aged children. However, its main impact factors and genetic characteristics remain poorly understood in China. Hence, the purpose of the study was to explore the precise association between multiple impact factors and head lice infestation, characterize the genetic variation of the head lice, and develop a sensitive and specific mitochondrial (mt) DNA duplex polymerase chain reaction (PCR) for accurately distinguishing clades A and B.

METHODS

A cross-sectional study was conducted in Hunan Province, central China from January 2018 to July 2023. A total of 9254 schoolchildren from 48 primary schools in each administrative region were examined for head lice. Impact factors for infestation were analyzed using the data collected by a questionnaire. The mt cytb gene sequences of head lice collected in the current study were used for sequence analysis, then were added to the global pool to conduct the phylogenetic analyses. Primers designed on the basis of this gene sequence were used in duplex PCR to diagnose head lice clades A and B by amplicon size.

RESULTS

Head lice infestation was found in 93.8% (45/48) of the primary schools included in the study. Overall, 6.8% (630/9254) of the examined schoolchildren harbored head lice, with 94.6% (596/630) being girls. A total of 2132 adult head lice were collected from 630 infested cases. The impact factors for head lice infestation included gender, school location, family situation, per capita income, study mode, and hair washing per week (p < 0.01). However, season and age were not considered as impact factors for head lice prevalence (p > 0.05). Phylogenetic analysis based on mt cytb gene sequences showed that head lice are classified into two clades (A and B), with clade B being more dominant in Hunan Province, central China. The newly developed duplex PCR was able to differentiate clades A from B in China with 100% sensitivity and specificity.

CONCLUSIONS

Our findings revealed that head lice infestation is mostly associated with poverty and poor hygiene in Hunan Province, central China. It is crucial to consider the simultaneous surveillance of head lice infestation in schoolchildren in regions with low level of socioeconomic status; however, datasets from other provinces are warranted to confirm the findings. It further showed that clades A and B are common in central China and that the latter has emerged and become the dominant one.

A Genomic Surveillance Circuit for Emerging Viral Pathogens

Genomic surveillance has been crucial in monitoring the evolution and spread of SARS-CoV-2. In Andalusia (Spain), a coordinated genomic surveillance circuit was established to systematically sequence and analyze viral genomes across the region. This initiative organizes sample collection through 27 hospitals, which act as regional hubs within their respective health districts. Sequencing is performed at three reference laboratories, with downstream data analysis and reporting centralized at a bioinformatics platform. From 2021 to 2025, over 42,500 SARS-CoV-2 genomes were sequenced, enabling the identification of major variants and their evolutionary dynamics. The circuit tracked the transition from Alpha and Delta to successive Omicron waves, including both recombinant and non-recombinant clades. The integration of genomic and epidemiological data facilitated rapid variant detection, outbreak investigation, and public health decision making. This surveillance framework at a regional granularity demonstrates the feasibility of large-scale sequencing within a decentralized healthcare system and has expanded to monitor other pathogens, reinforcing its value for epidemic preparedness. Continued investment in genomic surveillance is critical for tracking viral evolution, guiding interventions, and mitigating future public health threats.

The vertebrate small leucine-rich proteoglycans: amplification of a clustered gene family and evolution of their transcriptional profile in jawed vertebrates

Small Leucine-Rich Proteoglycans (SLRPs) are a major family of vertebrate proteoglycans. In bony vertebrates, SLRPs have a variety of functions from structural to signaling and are found in extracellular matrices, notably in skeletal tissues. However, there is little or no data on the diversity, function and expression patterns of SLRPs in cartilaginous fishes, which hinders our understanding of how these genes evolved with the diversification of vertebrates, in particular regarding the early events of whole-genome duplications that shaped gnathostome and cyclostome genomes. We used a selection of chromosome-level assemblies of cartilaginous fish and other vertebrate genomes for phylogeny and synteny reconstructions, allowing better resolution and understanding of the evolution of this gene family in vertebrates. Novel SLRP members were uncovered together with specific loss events in different lineages. Our reconstructions support that the canonical SLRPs have originated from different series of tandem duplications that preceded the extant vertebrate last common ancestor, one of them even preceding the extant chordate last common ancestor. They then further expanded with additional tandem and whole-genome duplications during the diversification of extant vertebrates. Finally, we characterized the expression of several SLRP members in the small-spotted catshark Scyliorhinus canicula and from this, inferred conserved and derived SLRP expression in several skeletal and connective tissues in jawed vertebrates.

Dissecting Metabolic Functions and Sugar Transporters Using Genome and Transportome of Probiotic Limosilactobacillus fermentum KUB-D18

Background/Objectives:Limosilactobacillus fermentum KUB-D18, a heterofermentative lactic acid bacterium with promising probiotic properties, is known for promoting gut health and nutrient absorption. Originally isolated from chicken intestines, this strain demonstrates versatile metabolic capabilities in diverse gastrointestinal environments. However, the metabolic functions and sugar transport-related genes remain largely unexplored. This study thus aimed to dissect metabolic functions and sugar transports of L. fermentum KUB-D18. Methods: Next-generation and third-generation sequencing techniques using integrative genomic platform towards transportome analysis were performed. Results: The complete genome, sized at 2.12 Mbps with a GC content of 51.36%, revealed 2079 protein-encoding genes, of which 1876 protein functions were annotated and identified in top categories involved in amino acids, nucleotide, energy, and carbohydrate transports and metabolisms. Comparative genes analysis identified 50 core and 12 strain-specific genes linked to probiotic properties, e.g., acid resistances and bile tolerances, antioxidant functions, or anti-inflammatory properties. Further, sugar transportome analysis uncovered 57 transporter genes, demonstrating diverse carbon utilization and phosphotransferase (PTS) systems, corroborated by API 50 CHL test results for carbohydrate metabolism profile. Conclusions: These findings enhance the comprehensive metabolic understanding of L. fermentum KUB-D18, supporting its industrial potential and applications in engineered probiotics.

Morphological and multi-locus phylogenetic analyses reveal three new branched species of Clavaria (Clavariaceae, Agaricales) from China

Based on morphological and molecular evidence, 12 specimens have been identified as belonging to three previously unrecognized species of Clavaria, which are here described as C.divergens, C.orientalis, and C.tongdaoensis. Clavariadivergens is characterized by its branched, white basidiomata. Clavariaorientalis and C.tongdaoensis are very similar to C.zollingeri in the field. However, C.orientalis is distinguished by its more robust branches, while C.tongdaoensis differs by its varied or paler color of basidiomata. A concatenated sequence dataset (ITS-nrLSU-RPB2) was used for multi-locus phylogenetic analysis. The phylogenetic tree of Clavaria showed that the three branched species each formed a distinct lineage with strong support. A key to the known branched species of Clavaria in China is provided.

A phoenix in the greenhouse: characterization and phylogenomics of complete chloroplast genomes sheds light on the putatively extinct-in-the-wild Solanum ensifolium (Solanaceae)

BACKGROUND

The genus Solanum is a diverse group of flowering plants with significant economic importance. Within this genus, the subgenus Leptostemonum, comprising spiny solanums, is particularly noteworthy due to its high species diversity and endemism. Solanum ensifolium, a member of this subgenus, is a critically endangered species endemic to Puerto Rico and known locally as erubia. The species survives in greenhouses and botanical gardens and is thought to be extinct in the wild, but with reintroduction efforts in progress. Despite its conservation status, genomic data for S. ensifolium remains scarce, limiting our understanding of its evolutionary history and potential adaptations.

RESULTS

The S. ensifolium chloroplast genome (155,295 bp) exhibits a typical quadripartite structure and encodes 151 genes, including 95 protein-coding genes involved in photosynthesis, transcription, translation, and other essential cellular functions. Gene content and genome organization are similar to those observed in closely related Solanum species. Comparative genomic analysis of the annotated genome with that of closely related Solanum species revealed differences in nucleotide diversity between the large single-copy (LSC) and small single-copy regions (SSC), and the inverted repeat (IR) regions. Additionally, phylogenetic analyses confirmed placement of S. ensifolium within the Leptostemonum subgenus, affirming its suspected close relationship with S. crotonoides and S. aturense. Furthermore, of the three individuals of S. ensifolium for which chloroplast genomes were obtained, no genetic variation was observed.

CONCLUSIONS

The availability of the S. ensifolium chloroplast genome provides insights into its evolutionary history and conservation needs. Comparative genomics uncovered evolutionary differences in Solanum chloroplast genomes, including nucleotide diversity and structural variations. Phylogenetic analyses confirmed the close relationship between S. ensifolium and other Leptostemonum species. These findings enhance our understanding of this critically endangered species' evolution, guiding effective conservation strategies like using chloroplast variation to assess genetic diversity for ex situ conservation and reintroduction programs. The uniformity of the chloroplast genome in S. ensifolium may reveal that this species has undergone a genetic bottleneck. To prevent inbreeding depression and maintain evolutionary adaptability, efforts should be made to generate and preserve as much genetic diversity as possible.

Complete Mitochondrial Genome of Platygyra daedalea and Characteristics Analysis of the Mitochondrial Genome in Merulinidae

BACKGROUND

The Merulinidae family belonging to the order Scleractinia is mainly distributed in the Indo-Pacific and Caribbean regions and often constitute the most dominant species of coral reefs. Mitochondrial genome is a key tool for studying the phylogeny and adaptation. Only a few studies have conducted the characteristics analyses of mitochondrial genome in the Merulinidae family.

METHODS

Therefore, we used high-throughput sequencing technology to describe the mitochondrial genome of Platygyra daedalea, a member of this family. Bioinformatics was used to analyze the composition characteristics of the mitochondrial genome of 10 Merulinidae species.

RESULTS

The mitochondrial genome of P. daedalea had a total length of 16,462 bp and a GC content of 33.0%. Thirteen unique protein-coding genes (PCGs), two transfer RNA (tRNA) genes, and two ribosomal RNA (rRNA) genes were annotated. Each species of Merulinidae had 13 unique PCGs in the mitochondrial genome. In contrast, the number of tRNAs and rRNAs significantly varied in Merulinidae species. Collinearity and gene rearrangement analyses indicated that the mitochondrial evolution of species in the Merulinidae family was relatively conserved. Divergence time analysis indicated that Merulinidae originated in the Oligocene, whereas the Platygyra genus originated in the Miocene. The formation and intraspecific divergence of coral species were consistent with geological changes in the ocean.

CONCLUSIONS

The results of this study help better understand the characteristics of the mitochondrial genome in the Merulinidae family and provide insights into the utility of mitochondrial genes as molecular markers of phylogeny.

Backbone phylogeny of Salix based on genome skimming data

The genus Salix is a common component of the Northern Hemisphere dendroflora with important ecological and economic value. However, taxonomy and systematics of Salix is extremely difficult and relationships between main lineages, especially deep phylogenies, remain largely unresolved. In this study, we used genome-skimming, plastome assembly, and single-copy orthologs (SCOs) from 66 Salix accessions, along with publicly available plastome and sequence read archive (SRA) datasets to obtain a robust backbone phylogeny of Salix, clarify relationships between its main lineages, and gain a more precise understanding of the origin and diversification of this species-rich genus. The plastome and SCO datasets resolved Salix into two robust clades, with plastome-based phylogenies lacking inner resolution and SCO offering fully resolved phylogenies. Our results support the classification of Salix into five subgenera: Salix, Urbaniana, Triandrae, Longifoliae and Vetrix. We observed a significant acceleration in the diversification rate within the Chamaetia-Vetrix clade, while Salix exhibited increased rates of diversification spanning from the early Oligocene to the late Miocene. These changes coincided with contemporaneous tectonic and climate change events. Our results provide a foundation for future systematic and evolutionary studies of Salix. Additionally, we showed that genome skimming data is an efficient, rapid, and reliable approach for obtaining extensive genomic data for phylogenomic studies, enabling the comprehensive elucidation of Salix relationships.

JN.1 variants circulating in Italy from October 2023 to April 2024: genetic diversity and immune recognition

BACKGROUND

The continuous emergence of SARS-CoV-2 variants and subvariants poses significant public health challenges. The latest designated subvariant JN.1, with all its descendants, shows more than 30 mutations in the spike gene. JN.1 has raised concerns due to its genomic diversity and its potential to enhance transmissibility and immune evasion. This study aims to analyse the molecular characteristics of JN.1-related lineages (JN.1*) identified in Italy from October 2023 to April 2024 and to evaluate the neutralization activity against JN.1 of a subsample of sera from individuals vaccinated with XBB.1.5 mRNA.

METHODS

The genomic diversity of the spike gene of 794 JN.1* strain was evaluated and phylogenetic analysis was conducted to compare the distance to XBB.1.5. Moreover, serum neutralization assays were performed on a subsample of 19 healthcare workers (HCWs) vaccinated with the monovalent XBB.1.5 mRNA booster to assess neutralizing capacity against JN.1.

RESULTS

Sequence analysis displayed high spike variability between JN.1* and phylogenetic investigation confirmed a substantial differentiation between JN.1* and XBB.1.5 spike regions with 29 shared mutations, of which 17 were located within the RBD region. Pre-booster neutralization activity against JN.1 was observed in 42% of HCWs sera, increasing significantly post-booster, with all HCWs showing neutralization capacity three months after vaccination. A significant correlation was found between anti-trimeric Spike IgG levels and neutralizing titers against JN.1.

CONCLUSIONS

The study highlights the variability of JN.1* in Italy. Results on a subsample of sera from HCWs vaccinated with XBB.1.5 mRNA booster vaccine suggested enhanced neutralization activity against JN.1.

Whole plastome sequence of Hibiscus moscheutos L. 1753 (Malvaceae) and its phylogenetic analysis

Hibiscus moscheutos, a perennial herb in the Malvaceae family, has medicinal properties but lacks molecular data. This study sequenced and analyzed its first complete plastome using next-generation sequencing. The 160,208 bp circular plastome has a typical quadripartite structure with 130 genes and 36.93% GC content. Phylogenetic analysis showed that H. mutabilis is sister to H. taiwanensis. This study provides essential molecular data for future research on Hibiscus phylogeny, taxonomy, and evolution.

Symbionts of predatory protists are widespread in the oceans and related to animal pathogens

Protists are major predators of ocean microbial life, with an ancient history of entanglements with prokaryotes, but their delicate cell structures and recalcitrance to culturing hinder exploration of marine symbioses. We report that tiny oceanic protistan predators, specifically choanoflagellates-the closest living unicellular relatives of animals-and uncultivated MAST-3 form symbioses with four bacterial lineages related to animal symbionts. By targeting living phagotrophs on ship expeditions, we recovered genomes from physically associated uncultivated Legionellales and Rickettsiales. The evolutionary trajectories of Marinicoxiellaceae, Cosmosymbacterales, Simplirickettsiaceae, and previously named Gamibacteraceae vary, including host-engagement mechanisms unknown in marine bacteria, horizontally transferred genes that mediate pathogen-microbiome interactions, and nutritional pathways. These symbionts and hosts occur throughout subtropical and tropical oceans. Related bacteria were detected in public data from freshwater, fish, and human samples. Symbiont associations with animal-related protists, alongside relationships to animal pathogens, suggest an unexpectedly long history of shifting associations and possibilities for host expansion as environments change.

Insight on the heterocyte patterning and the proheterocyte division in the toxic cyanobacterium Kaarinaea lacus gen. nov., sp. nov., and its genomic potential for natural products

Nostoc sp. 152 (= PCC 9237T), a toxic cyanobacterium, isolated from mixed cyanobacterial bloom in a Finnish freshwater lake was reassessed using phylogenetic, morphological and genomic analyses. Multilocus and 16S rRNA gene phylogenetic analyses confirmed that this strain represents a novel Nostocalean genus for which we propose the name Kaarinaea lacus gen. nov., sp. nov., in accordance with the International Code of Nomenclature of Prokaryotes. The most intriguing morphological feature exhibited by PCC 9237T is the occasional division of proheterocytes. Based on our observations, we propose a hypothesis for the sequence of events taking place leading to this phenomenon. The other interesting feature includes the unusual heterocyte patterning resulting in the development of heterocytes in series in the old culture followed by its fragmentation at this site. Among the genes involved in heterocyte differentiation and patterning, patC and three homologs of hetP were not found from the genome of PCC 9237T. Furthermore, genomic investigations revealed a variant of heterocyte glycolipid gene cluster with a reduced hglB, but having additional gene coding for a protein with TubCN terminal docking domain. Finally, the presence of a pks2-like gene cluster, whose product may interfere with the cellular differentiation, and an anachelin-like gene cluster demands further investigations.

Insight into the codon usage patterns and adaptation of Tembusu Virus

Since its emergence in 2010, Tembusu virus (TMUV) has rapidly spread across poultry populations in Asia, leading to substantial economic losses in these areas. Here, we investigate the codon usage patterns (CUPs) underlying TMUV's adaptation and evolutionary dynamics within host environments. Phylogenetic and compositional analyses consistently classify TMUV into four evolutionary lineages-Clusters 1, 2, 3, and ancestral TMUV-with Cluster 2 emerging as the dominant lineage. Codon adaptation index (CAI) analysis reveals that this lineage of TMUV show best adapted to the CUPs of duck than other lineages, underscoring the role of natural selection in shaping viral evolution, a finding in line with evidence that CUPs in the TMUV genome is predominantly shaped by natural selection. Furthermore, TMUV exhibits markedly higher adaptation to the CUPs of poultry hosts (duck, goose, and chicken) compared to potential host humans or vector mosquito. Thus, species-specific adaptability to the host environment may be a reason account for the distinct infectivity and clinic outcome of TMUV acted on hosts. Analysis of dinucleotide distribution reveals significant suppression of CpG and UpA dinucleotides in the TMUV genome, reflecting adaptive pressures to evade vertebrate immune responses. During transmission, TMUV shows increasing alignment with host CUPs and a continuous reduction in CpG dinucleotides, potentially enhancing its fitness within host microenvironments. This work advances our understanding of the basic biology underlying TMUV epidemiology, pathogenicity, and species-specific adaptation.

Protists as mediators of complex microbial and viral associations

Microbial eukaryotes (aka protists) are known for their important roles in nutrient cycling across different ecosystems. However, the composition and function of protist-associated microbiomes remains largely elusive. Here, we employ cultivation-independent single-cell isolation and genome-resolved metagenomics to provide detailed insights into underexplored microbiomes and viromes of over 100 currently uncultivable ciliates and amoebae isolated from diverse environments. Our findings reveal unique microbiome compositions and hint at an intricate network of complex interactions and associations with bacterial symbionts and viruses. We observed stark differences between ciliates and amoebae in terms of microbiome and virome compositions, highlighting the specificity of protist-microbe interactions. Over 115 of the recovered microbial genomes were affiliated with known endosymbionts of eukaryotes, including diverse members of the Holosporales, Rickettsiales, Legionellales, Chlamydiae, Dependentiae , and more than 250 were affiliated with possible host-associated bacteria of the phylum Patescibacteria. We also identified more than 80 giant viruses belonging to diverse viral lineages, of which some were actively expressing genes in single cell transcriptomes, suggesting a possible association with the sampled protists. We also revealed a wide range of other viruses that were predicted to infect eukaryotes or host-associated bacteria. Our results provide further evidence that protists serve as mediators of complex microbial and viral associations, playing a critical role in ecological networks. The frequent co-occurrence of giant viruses and diverse microbial symbionts in our samples suggests multipartite associations, particularly among amoebae. Our study provides a preliminary assessment of the microbial diversity associated with lesser-known protist lineages and paves the way for a deeper understanding of protist ecology and their roles in environmental and human health.

First Report and Phylogenetic Analysis of Mitochondrial Genomes of Chrysomya villeneuvi and Sarcophaga genuforceps

The mitochondrial genome, highly conserved across species, is crucial for species identification, phylogenetic analysis, and evolutionary research. Chrysomya villeneuvi and Sarcophaga genuforceps, two species with significant forensic value, have been understudied in terms of genetic data. In this study, the complete mitochondrial genomes of C. villeneuvi (15,623 bp) and S. genuforceps (15,729 bp) were sequenced and analyzed. All thirteen protein-coding genes (PCGs) exhibited Ka/Ks ratios below one, indicating purifying selection and supporting their utility as barcoding markers. Phylogenetic analysis and genetic distance calculations based on PCGs showed that C. villeneuvi is closely related to Chrysomya rufifacies and Chrysomya albiceps, and S. genuforceps aligns more closely with Sarcophaga kentejana and Sarcophaga schuetzei. This research is the first to provide mitochondrial genome data for C. villeneuvi and S. genuforceps, expanding the genetic resources available for Calliphoridae and Sarcophagidae and offering a foundation for further forensic and evolutionary studies.

Description of Cohnella rhizoplanae sp. nov., isolated from the root surface of soybean (Glycine max)

A Gram-staining-positive, aerobic bacterium, designated strain JJ-181 T, was isolated from the root surface of soybean. Based on the 16S rRNA gene sequence similarities, strain JJ-181 T was grouped into the genus Cohnella, most closely related to Cohnella hashimotonis F6_2S_P_1T (98.85%) and C. ginsengisoli GR21-5 T (98.3%). The pairwise average nucleotide identity and digital DNA-DNA hybridisation values of the JJ-181 T genome assembly against publicly available Cohnella type strain genomes were below 84% and 28%, respectively. The fatty acid profile from whole cell hydrolysates, the cell wall diaminoacid, the quinone system, the polar lipid profile, and the polyamine pattern supported the allocation of strain JJ-181 T to the genus Cohnella. In addition, the results of physiological and biochemical tests also allowed phenotypic differentiation of strain JJ-181 T from its closely-related Cohnella species. Functional analysis revealed that strain JJ-181 T has different gene clusters related to swarming motility, chemotaxis ability, and endospore and biofilm formations. The gene content also suggests the ability of strain JJ-181 T to synthesise biotin and riboflavin, as well as indole-3-acetic acid, an important phytohormone for plant growth. Based on polyphasic analyses, strain JJ-181 T can be classified as a new species of the genus Cohnella, for which we propose the name Cohnella rhizoplanae sp. nov., with strain JJ-181 T (= LMG 31678 T = CIP 112018 T = CCM 9031 T = DSM 110650 T) as the type strain.

Phylogenetic and Pathogenic Analysis of H5N1 and H5N6 High Pathogenicity Avian Influenza Virus Isolated from Poultry Farms (Layer and Broiler Chickens) in Japan in the 2023/2024 Season

During the 2023-2024 winter, 11 high pathogenicity avian influenza (HPAI) outbreaks caused by clade 2.3.4.4b H5N1 and H5N6 HPAI viruses were confirmed in Japanese domestic poultry among 10 prefectures (n = 10 and 1, respectively). In this study, we aimed to genetically and pathologically characterize these viruses. Phylogenetic analysis revealed that H5N1 viruses were classified into the G2d-0 genotype, whereas the H5N6 virus was a novel genotype in Japan, designated as G2c-12. The G2c-12 virus shared PB2, PB1, PA, HA, and M genes with previous G2c viruses, but had NP and NS genes originating from avian influenza viruses in wild birds abroad. The N6 NA gene was derived from an H5N6 HPAI virus that was different from the viruses responsible for the outbreaks in Japan in 2016-2017 and 2017-2018. Experimental infections in chickens infected with H5N1(G2d-0) and H5N6(G2c-12) HPAI viruses showed no significant differences in the 50% chicken lethal dose, mean death time, or virus shedding from the trachea and cloaca, or in the histopathological findings. Different genotypes of the viruses worldwide, their introduction into the country, and their stable lethality in chickens may have triggered the four consecutive seasons of HPAI outbreaks in Japan.

Evolution of copper resistance in Xanthomonas euvesicatoria pv. perforans population

The widespread use of antimicrobials that target bacterial pathogens has driven evolution of resistance, compromising the efficacy of these bactericides. Understanding the emergence and spread of resistance genes via mobile genetic elements is crucial for combating antimicrobial resistance. Copper resistance (CuR) in Xanthomonas euvesicatoria pv. perforans has severely affected the efficacy of copper-based bactericides for controlling bacterial leaf spot disease of tomato and pepper. Here, we investigated the evolutionary pathways of CuR acquisition and dissemination in X. euvesicatoria pv. perforans using an extensive collection of strains. We determined that chromosomally encoded CuR predominates over plasmid-borne CuR in multiple distinct phylogenetic groups of X. euvesicatoria pv. perforans. Our analysis revealed a single site of chromosomal integration by a CuR genomic island, although the genomic island showed sequence variation among phylogenetic groups. While chromosomal CuR was more prevalent, strains with plasmid-borne resistance conferred greater copper tolerance. Additionally, we identified strains carrying two copies of CuR genes, on plasmid and chromosome, that exhibited increased copper tolerance. Strains of X. euvesicatoria pv. perforans from the USA shared identical CuR gene sequences whether on plasmids or chromosome while different alleles were found in strains from other countries. In contrast to X. euvesicatoria pv. perforans, plasmid-borne CuR predominated in closely related pathovar, X. euvesicatoria pv. euvesicatoria. Overall, these findings contribute to a better understanding of the evolution and persistence of CuR in X. euvesicatoria pv. perforans and its closest relatives.IMPORTANCEThe emergence of antimicrobial resistance is a significant threat to agricultural production as it reduces the efficacy of various antimicrobials including copper-based bactericides that are widely used to control plant diseases. The challenge of increasing antimicrobial resistance entering a production system necessitates a deeper understanding of the dynamics and mechanisms by which pathogens acquire resistance. As a result of this research, we have identified different mechanisms of copper resistance acquisition as well as levels of copper resistance in a devastating plant pathogen, X. euvesicatoria pv. perforans. The evolution and dissemination of copper resistance in strains through plasmid or chromosomally integrated genomic island or both presents barriers to current management approaches, where growers rely heavily on copper-based bactericides to manage disease outbreaks. This knowledge is crucial when considering the continued use of existing antimicrobials or adopting alternative antimicrobials in efforts to implement enhanced antimicrobial stewardship strategies in agriculture.

Biodiversity of Herbivores Triggers Species Differentiation of Coprophilous Fungi: A Case Study of Snow Inkcap (Coprinopsis sect. Niveae)

Coprophilous species of Coprinopsis sect. Niveae, commonly known as "snow inkcap", are widespread in pastoral areas; however, wide sampling approaches are needed to discover new taxa and to clarify the taxonomic status of the so-called "snow inkcap". Nationwide field work was conducted in China with a detailed record collected of the distribution and the animal origin of the dung. A four-loci phylogenetic study of Coprinopsis sect. Niveae was conducted based on the internal transcribed spacer regions (ITS), the ribosomal large subunit (LSU), and translation elongation factor 1-α (tef1-α)], and the mitochondria small ribosomal RNA subunit (mtSSU). Fourteen phylogenetic species were assigned to this section, including six novel species, namely Coprinopsis furfuracea, C. iliensis, C. khorqinensis, C. sericivia, C. subigarashii, and C. tenuipes. Macro-, micro-, and ultramicro-morphological observations of species collected from China were also conducted and the detailed descriptions and illustrations of the novel species are provided. Our studies revealed that the different origin of herbivore dung, the distribution, the color and thickness of the pileus, the shape of stipes, the shape and size of basidiospores, and the presence or absence of pleurocystidia can be used as characteristics for distinguishing species in sect. Niveae. The key to species belonging to this section is also provided.

New Contributions on Species Diversity of Genus Hydnum and Lentaria s.l. in China

Southwest China is extremely rich in fungal resources, and a large number of new taxa have been discovered in recent years. In the present study, we examined 26 specimens of the genera Hydnum and Lentaria sensu lato, most of which were obtained in Yunnan Province. Through ITS-nrLSU-tef1 phylogenetic analysis, combined with morphological studies and geographic analyses, five new species were described, viz. H. cremeum (), H. flavoquamosum (), H. roseoalbum (), H. roseotangerinum (), and L. subalpina (). Furthermore, we also supplied new information on some known species, including host plants and new distribution records. We re-examined the holotype sequences of two known taxa, H. flabellatum and H. pallidomarginatum, treating H. flabellatum as a synonym of H. pallidomarginatum. Additionally, a key to Lentaria s.l. in China was provided.

Arsenic-induced enhancement of diazotrophic recruitment and nitrogen fixation in Pteris vittata rhizosphere

Heavy metal contamination poses an escalating global challenge to soil ecosystems, with hyperaccumulators playing a crucial role in environmental remediation and resource recovery. The enrichment of diazotrophs and resulting nitrogen accumulation promoted hyperaccumulator growth and facilitated phytoremediation. Nonetheless, the regulatory mechanism of hyperaccumulator biological nitrogen fixation has remained elusive. Here, we report the mechanism by which arsenic regulates biological nitrogen fixation in the arsenic-hyperaccumulator Pteris vittata. Field investigations and greenhouse experiments, based on multi-omics approaches, reveal that elevated arsenic stress induces an enrichment of key diazotrophs, enhances plant nitrogen acquisition, and thus improves plant growth. Metabolomic analysis and microfluidic experiments further demonstrate that the upregulation of specific root metabolites plays a crucial role in recruiting key diazotrophic bacteria. These findings highlight the pivotal role of nitrogen-acquisition mechanisms in the arsenic hyperaccumulation of Pteris vittata, and provide valuable insights into the plant stress resistance.

Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

New Helminthosporium (Massarinaceae, Dothideomycetes) and Nigrospora (Incertae sedis, Sordariomycetes) species associated with walnut (Juglansregia L.) in China

Six collections of ascomycetes were obtained from samples collected from dead branches and leaves of Juglansregia in Guizhou and Yunnan provinces, China. By incorporating multigene phylogenetic analysis (ITS, LSU, rpb2, SSU, tef1-α, tub2) supplemented by morphological data, we establish two novel species, namely Helminthosporiumguizhouense and Nigrosporayunnanensis. In morphology, H.guizhouense can be distinguished from H.caespitosum by its narrower conidia (13-16 µm vs. 27.3-35.5 µm), and N.yunnanensis is characterized by black, globose conidia (16.2 × 14.4 µm). The phylogenetic results further substantiated them as novel taxa. The present study contributes to our comprehension of the range of fungi found in Juglansregia, thereby expanding our knowledge of the diversity of fungi within this host.

Phylogenetic and Morphological Perspectives on Crepidotus subg. Dochmiopus: Exploratively Unveiling Hidden Diversity in China

Crepidotus subg. Dochmiopus contributes to more than half of Crepidotus species and exhibits highly hidden diversity. However, C. subg. Dochmiopus is challenging to study because the basidiomata of C. subg. Dochmiopus species are usually small and white, inconspicuous interspecific distinctions, and possess a familiar complex. In this study, we utilized a variety of characteristics for species identification, including habitat, presence or absence of a stipe in mature specimens, pileipellis and cheilocystidia patterns, whether the lamellae edges are fimbriated, and other characteristics. Above all, cheilocystidia and pileipellis patterns will be important in C. subg. Dochmiopus research. Based on the present specimens, we constructed a multigene phylogenetic tree (ITS + LSU) and recognized four new species: C. lamellomaculatus sp. nov., C. capitatocystidiatus sp. nov., C. succineus sp. nov., C. clavocystidiatustustus sp. nov. Detailed morphological descriptions, photographs, line drawings and comparisons with closely related taxa for the new species are provided. The current phylogenetic analysis does not support the previously classifications, indicating that the classification of Crepidotus requires re-evaluation. But the existing molecular datasets and species' descriptions are insufficient to fully resolve the classification. Further integration of new gene segments and a comprehensive review of morphological characteristics will reveal a natural classification for Crepidotus.

Multi-gene analysis of the Russula crown clade (Russulales, Basidiomycota) revealed six new species and Alboflavinae subsect. nov. from Fagaceae forests in China

Introduction

The crown clade is one of two major groups in the Russula subg. Russula.

Methods/material

An analysis of Chinese samples was performed based on the morphology, internal transcribed spacer (ITS) sequences, and multi-gene phylogenies of 28S nrLSU, 16S mtSSU, rpb1, rpb2, and tef1-α.

Results

The results supported the independence of six new species: Russula alboflava (sect. Amethystinae), R. chrysantha (subsect. Chamaeleontinae), R. liyui (subsect. Laricinae), R. lutescens (subsect. Olivaceinae), R. paraxerampelina, and R. prunicolor (subsect. Xerampelinae) from Fagaceae forest habitats. Subsect. Alboflavinae was newly proposed in sect. Amethystinae. Members of the new subsection include R. alboflava, R. burlinghamiae, and possibly R. ballouii.

Discussion

Our analyses also supported the claim that two species of R. fulvograminea (subsect. Laricinae) and R. subrubens (subsect. Xerampelinae) have a Eurasian distribution. The habitat and primary hosts of the main phylogenetic clades within related subsections were summarized and discussed.

The complete chloroplast genome of Camellia huulungensis Rosmann et Ninh, a golden Camellia species endemic to Vietnam

Camellia huulungensis Rosmann & Ninh 1997, belonging to the sect. Chrysantha, holds important ornamental value and medicinal value. In this study, the complete chloroplast genome sequence of C. huulungensis was assembled using high-throughput sequencing technology. The entire length of chloroplast genome is 156,546 bp and contains a small single-copy region (18,257 bp), a large single-copy region (86,219 bp), and a pair of inverted repeat regions (26,035 bp). A total of 133 genes were annotated, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The overall GC content is 37.33%. The phylogenetic analysis showed that C. huulungensis is sister to C. aurea. The results can provide genetic data for further phylogenetic studies of Camellia.

The complete chloroplast genome of Tabebuia rosea (Bignoniaceae)

Tabebuia rosea is a world-renowned woody plant with colorful flowers in full bloom. In addition to its high ornamental value, it also has ecological and medicinal value. In this study, the complete circular chloroplast genome of T. rosea was reconstructed and annotated using Illumina sequencing. The chloroplast genome was 158,919 bp in size with GC content of 38.21%, including a large single-copy region of 85,823 bp, a small single-copy region of 12,816 bp, and a pair of inverted repeats of 30,140 bp. It encoded 132 genes, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Based on current available chloroplast genome sequences, the phylogenetic analysis indicated that T. rosea was clustered with T. nodosa and H. chrysanthus. This study provided insights into the evolutionary relationships among different species of Bignoniaceae.

New insights into honey bee viral and bacterial seasonal infection patterns using third-generation nanopore sequencing on honey bee haemolymph

Honey bees are rapidly declining, which poses a significant threat to our environment and agriculture industry. These vital insects face a disease complex believed to be caused by a combination of parasites, viruses, pesticides, and nutritional deficiencies. However, the real aetiology is still enigmatic. Due to the conventional analysis methods, we still lack complete insights into the honey bee virome and the presence of pathogenic bacteria. To fill this knowledge gap, we employed third-generation nanopore metagenomic sequencing on honey bee haemolymph to monitor the presence of pathogens over almost a year. This study provides valuable insights into the changes in bacterial and viral loads within honey bee colonies. We identified different pathogens in the honey bee haemolymph, which are not included in honey bee screenings. These pathogens comprise the Apis mellifera filamentous virus, Apis rhabdoviruses, and various bacteria such as Frischella sp. and Arsenophonus sp. Furthermore, a sharp contrast was observed between young and old bees. Our research proposes that transgenerational immune priming may play a role in shaping infection patterns in honey bees. We observed a significant increase in pathogen loads in the spring, followed by a notable decrease in pathogen presence during the summer and autumn months. However, certain pathogens seem to be able to evade this priming effect, making them particularly intriguing as potential factors contributing to mortality. In the future, we aim to expand our research on honey bee transgenerational immune priming and investigate its potential in natural settings. This knowledge will ultimately enhance honey bee health and decrease colony mortality.

The complete plastome of Rorippa palustris Besser 1821 and its phylogenetic analysis

Rorippa palustris Besser 1821, a species of Brassicaceae, is widely distributed around the world and used for both food 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 a complete plastome of R. palustris, marking a significant advancement toward comprehending its genetic composition. The plastome is 154,674 bp long and harbors 128 genes, including 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Our phylogenomic analysis indicated that R. palustris is closely related to R. curvipes. 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 R. palustris, which could be advantageous in pharmaceutical applications.

Comparison of organelle genomes between endangered mangrove plant Dolichandrone spathacea to terrestrial relative provides insights into its origin and adaptative evolution

Dolichandrone spathacea is a mangrove associate with high medicinal and ecological values. However, due to the dual-pressure of climate change and human activities, D. spathacea has become endangered in China. Moreover, misidentification between D. spathacea and its terrestrial relative D. cauda-felina poses further challenges to field protection and proper medicinal usage of D. spathacea. Thus, to address these problems, we sequenced and assembled mitochondrial (mt) and chloroplast (cp) genomes for both D. spathacea and D. cauda-felina. Comparative analysis revealed apparently different size and scaffold number between the two mt genomes, but a high similarity between the cp genomes. Eight regions with high sequence divergence were identified between the two cp genomes, which might be used for developing candidate DNA markers for distinguishing the two species. The splitting between D. spathacea and D. cauda-felina was inferred to occur at ~6.8 - 7.7 million years ago (Mya), which may be driven by the environment fluctuations in late Miocene. In the cp genome, 12 genes related to the expression of photosynthesis-associated proteins were detected with signatures of positive selection, which may contribute to the origin and evolutionary adaptation of Dolichandrone mangrove species. These new findings do not only enrich organelle genomic resources of Dolichandrone species, but also provide important genetic clues for improving the conservation and proper usage of endangered mangrove associate D. spathacea.

Epidemiology and genetic characterization of influenza viruses circulating in Bhutan in 2022

INTRODUCTION

Influenza (Flu) causes considerable morbidity and mortality globally, and in Bhutan, Flu viruses are a leading cause of acute respiratory infection and cause outbreaks during Flu seasons. In this study, we aim to analyze the epidemiology and the genetic characterization of Flu viruses circulated in Bhutan in 2022.

METHOD

Respiratory specimens were collected from patients who meet the case definition for influenza-like illness (ILI) and severe acute respiratory infection (SARI) from sentinel sites. Specimens were tested for Flu and SARS-CoV-2 viruses by RT-PCR using the Multiplex Assay. Selected positive specimens were utilized for Flu viral genome sequencing by next-generation sequencing. Descriptive analysis was performed on patient demographics to see the proportion of Flu-associated ILI and SARI. All data were analyzed using Epi Info7 and QGIS 3.16 software.

RESULT

A weekly average of 16.2 ILI cases per 1000 outpatient visits and 18 SARI cases per 1000 admitted cases were reported in 2022. The median age among ILI was 12 years (IQR: 5-28) and SARI was 6.2 (IQR: 2.5-15) years. Flu A(H3N2) (70.2%) subtype was the most predominant circulating strain. Flu A(H1N1)pdm09 and Flu B viruses belonged to subclades that were mismatched to the vaccine strains recommended for the 2021-2022 season but matched the vaccine strain for the 2022-2023 season with vaccine efficacy 85.14% and 88.07% respectively. Flu A(H3N2) virus belonged to two subclades which differed from the vaccine strains recommended in both the 2021-2022 and 2022-2023 seasons with vaccine efficacy 68.28%.

CONCLUSION

Flu virus positivity rates were substantially elevated during the Flu season in 2022 compared to 2021. Flu A(H3N2) subtype was the most predominant circulating strain in the country and globally. Genetic characterization of the Flu viruses in Bhutan showed a close relatedness of high vaccine efficacy with the vaccine strain that WHO recommended for the 2022-23 season.

SARS-CoV-2 variant replacement constrains vaccine-specific viral diversification

Coronavirus disease 2019 (COVID-19) vaccine breakthrough infections have been important for all circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant periods, but the contribution of vaccine-specific SARS-CoV-2 viral diversification to vaccine failure remains unclear. This study analyzed 595 SARS-CoV-2 sequences collected from the Military Health System beneficiaries between December 2020 and April 2022 to investigate the impact of vaccination on viral diversity. By comparing sequences based on the vaccination status of the participant, we found limited evidence indicating that vaccination was associated with increased viral diversity in the SARS-CoV-2 spike, and we show little to no evidence of a substantial sieve effect within major variants; rather, we show that rapid variant replacement constrained intragenotype COVID-19 vaccine strain immune escape. These data suggest that, during past and perhaps future periods of rapid SARS-CoV-2 variant replacement, vaccine-mediated effects were subsumed with other drivers of viral diversity due to the massive scale of infections and vaccinations that occurred in a short time frame. However, our results also highlight some limitations of using sieve analysis methods outside of placebo-controlled clinical trials.

Genetic diversity of H5N1 and H5N2 high pathogenicity avian influenza viruses isolated from poultry in Japan during the winter of 2022-2023

High pathogenicity avian influenza viruses (HPAIVs) of the H5N1 and H5N2 subtypes were responsible for 84 HPAI outbreaks on poultry premises in Japan during October 2022-April 2023. The number of outbreaks during the winter of 2022-2023 is the largest ever reported in Japan. In this study, we performed phylogenetic analyses using the full genetic sequences of HPAIVs isolated in Japan during 2022-2023 and those obtained from a public database to identify their genetic origin. Based on the hemagglutinin genes, these HPAIVs were classified into the G2 group of clade 2.3.4.4b, whose ancestors were H5 HPAIVs that circulated in Europe in late 2020, and were then further divided into three subgroups (G2b, G2d, and G2c). Approximately one-third of these viruses were classified into the G2b and G2d groups, which also included H5N1 HPAIVs detected in Japan during 2021-2022. In contrast, the remaining two-thirds were classified into the G2c group, which originated from H5N1 HPAIVs isolated in Asian countries and Russia during the winter of 2021-2022. Unlike the G2b and G2d viruses, the G2c viruses were first detected in Japan in the fall of 2022. Importantly, G2c viruses caused the largest number of outbreaks throughout Japan over the longest period during the season. Phylogenetic analyses using eight segment genes revealed that G2b, G2d, and G2c viruses were divided into 2, 4, and 11 genotypes, respectively, because they have various internal genes closely related to those of avian influenza viruses detected in wild birds in recent years in Asia, Russia, and North America, respectively. These results suggest that HPAIVs were disseminated among migratory birds, which may have generated numerous reassortant viruses with various gene constellations, resulting in a considerable number of outbreaks during the winter of 2022-2023.

Three new ramarioid species of Phaeoclavulina (Gomphaceae, Gomphales) from China

Three new species of Phaeoclavulina from China are described: Phaeoclavulinabicolor, P.echinoflava, and P.jilinensis. Recognition of the new species is supported by morphological and molecular evidence. Phylogenetic analyses of concatenated ITS1-5.8S-ITS2 and nuclear large subunit sequences support the establishment of the new species and their placement within the Phaeoclavulina clade. A key to the known Phaeoclavulina species in China is provided.

The complete chloroplast genome of Camellia flava (Pitard) Sealy, a golden camellia of Vietnam

Camellia flava (Pit.) Sealy 1949 is a rare and precious species with golden flowers, which hold important ornamental and breeding values. In this study, the complete chloroplast genome of C. flava is reported for the first time. The chloroplast genome exhibits a typical quadripartite structure with a total length of 156,670 bp and a GC content of 37.32%, including a large single-copy region (86,250 bp), a small single-copy region (18,292 bp), and a pair of inverted repeat regions (26,064 bp). A total of 133 genes, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes were annotated. The phylogenetic analysis revealed a close relationship between C. flava and C. tamdaoensis. The chloroplast genome sequence of C. flava serves as a valuable resource for further breeding research and genetic phylogenetic studies.

The complete chloroplast genome of Lindernia crustacea (L.) F. Muell 1882 (Linderniaceae) and its phylogenetic analysis

Lindernia crustacea (L.) F. Muell 1882, a species in the Linderniaceae family, holds traditional medicinal value in China. To investigate its genetic diversity, we assembled, annotated, and characterized the first complete chloroplast genome of L. crustacea using Illumina sequencing data and various bioinformatics tools. The genome is 153,647 bp in length, with a GC content of 37.6%. It exhibits a typical quadripartite structure, consisting of a large single-copy region (LSC) of 85,411 bp, a small single-copy region (SSC) of 18,724 bp, and two inverted repeat sequences (IRa and IRb) of 25,816 bp each. The genome was predicted to contain 131 genes, including 87 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenomic analysis indicated that L. crustacea is closely related to L. stricta. These findings provide a foundation for further research on the evolution and potential medicinal applications of the Linderniaceae family.

Phylogenomics including new sequence data of phytoplankton-infecting chytrids reveals multiple independent lifestyle transitions across the phylum

Parasitism is the most common lifestyle on Earth and has emerged many times independently across the eukaryotic tree of life. It is frequently found among chytrids (Chytridiomycota), which are early-branching unicellular fungi that feed osmotrophically via rhizoids as saprotrophs or parasites. Chytrids are abundant in most aquatic and terrestrial environments and fulfil important ecosystem functions. As parasites, they can have significant impacts on host populations. They cause global amphibian declines and influence the Earth's carbon cycle by terminating algal blooms. To date, the evolution of parasitism within the chytrid phylum remains unclear due to the low phylogenetic resolution of rRNA genes for the early diversification of fungi, and because few parasitic lineages have been cultured and genomic data for parasites is scarce. Here, we combine transcriptomics, culture-independent single-cell genomics and a phylogenomic approach to overcome these limitations. We newly sequenced 29 parasitic taxa and combined these with existing data to provide a robust backbone topology for the diversification of Chytridiomycota. Our analyses reveal multiple independent lifestyle transitions between parasitism and saprotrophy among chytrids and multiple host shifts by parasites. Based on these results and the parasitic lifestyle of other early-branching holomycotan lineages, we hypothesise that the chytrid last common ancestor was a parasite of phytoplankton.

Ecological genomics in the Northern krill uncovers loci for local adaptation across ocean basins

Krill are vital as food for many marine animals but also impacted by global warming. To learn how they and other zooplankton may adapt to a warmer world we studied local adaptation in the widespread Northern krill (Meganyctiphanes norvegica). We assemble and characterize its large genome and compare genome-scale variation among 74 specimens from the colder Atlantic Ocean and warmer Mediterranean Sea. The 19 Gb genome likely evolved through proliferation of retrotransposons, now targeted for inactivation by extensive DNA methylation, and contains many duplicated genes associated with molting and vision. Analysis of 760 million SNPs indicates extensive homogenizing gene-flow among populations. Nevertheless, we detect signatures of adaptive divergence across hundreds of genes, implicated in photoreception, circadian regulation, reproduction and thermal tolerance, indicating polygenic adaptation to light and temperature. The top gene candidate for ecological adaptation was nrf-6, a lipid transporter with a Mediterranean variant that may contribute to early spring reproduction. Such variation could become increasingly important for fitness in Atlantic stocks. Our study underscores the widespread but uneven distribution of adaptive variation, necessitating characterization of genetic variation among natural zooplankton populations to understand their adaptive potential, predict risks and support ocean conservation in the face of climate change.

Analysis of the complete mitochondrial genome sequence of Hipposideros pratti

In order to explore the characteristics of the mitochondrial genome sequence of Pratt's leaf-nosed bat (Hipposideros pratti Thomas 1891) and understand their phylogenetic status in Chiroptera, this study determined the mitochondrial genome sequences of H. pratti from five regions in China using high-throughput sequencing technology, sequence assembly, and genome annotation. The results showed that these sequences contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 non-coding region, all exhibiting a significant AT bias. Based on the phylogenetic tree constructed using 13 protein-coding genes from 15 Chiroptera species, the study found that H. pratti from the five regions clustered together, and then clustered with H. lylei into a single clade. Meanwhile, H. pratti from Jiangxi, Fujian, and Guangdong regions of China showed closer genetic relationships, while H. pratti from Yunnan and Henan regions of China exhibited closer genetic relationships. This study not only supplemented the mitochondrial genome database of H. pratti but also laid a foundation for genetic variation, molecular classification, and evolutionary studies of H. pratti.

The complete mitochondrial genome of the rodent flea Nosopsyllus laeviceps: genome description, comparative analysis, and phylogenetic implications

BACKGROUND

Fleas are one of the most common and pervasive ectoparasites worldwide, comprising at least 2500 valid species. They are vectors of several disease-causing agents, such as Yersinia pestis. Despite their significance, however, the molecular genetics, biology, and phylogenetics of fleas remain poorly understood.

METHODS

We sequenced, assembled, and annotated the complete mitochondrial (mt) genome of the rodent flea Nosopsyllus laeviceps using next-generation sequencing technology. Then we combined the new mitogenome generated here with mt genomic data available for 23 other flea species to perform comparative mitogenomics, nucleotide diversity, and evolutionary rate analysis. Subsequently, the phylogenetic relationship within the order Siphonaptera was explored using the Bayesian inference (BI) and maximum likelihood (ML) methods based on concentrated data for 13 mt protein-coding genes.

RESULTS

The complete mt genome of the rodent flea N. laeviceps was 16,533 base pairs (bp) in a circular DNA molecule, containing 37 typical genes (13 protein-coding genes, 22 transfer RNA [tRNA] genes, and two ribosomal RNA [rRNA] genes) with one large non-coding region (NCR). Comparative analysis among the order Siphonaptera showed a stable gene order with no gene arrangement, and high AT content (76.71-83.21%) with an apparent negative AT and GC skew except in three fleas Aviostivalius klossi bispiniformis, Leptopsylla segnis, and Neopsylla specialis. Moreover, we found robust evidence that the cytochrome c oxidase subunit 1 (cox1) gene was the most conserved protein-coding gene (Pi = 0.15, non-synonymous/synonymous [Ka/Ks] ratio = 0.13) of fleas. Phylogenomic analysis conducted using two methods revealed different topologies, but both results strongly indicated that (i) the families Ceratophyllidae and Leptopsyllidae were paraphyletic and were the closest to each other, and (ii) the family Ctenophthalmidae was paraphyletic.

CONCLUSIONS

In this study, we obtained a high-quality mt genome of the rodent flea N. laeviceps and performed comparative mitogenomics and phylogeny of the order Siphonaptera using the mt database. The results will enrich the mt genome data for fleas, lay a foundation for the phylogenetic analysis of fleas, and promote the evolutionary analysis of Siphonaptera.

Genetic characterization, structural analysis, and detection of positive selection in small heat shock proteins of Cypriniformes and Clupeiformes

In the current investigation, a total of 42 full-length, non-redundant small heat shock proteins (sHsp) were detected in Cyprinus carpio, Labeo rohita, Danio rerio, Salmo salar, Oncorhynchus mykiss, and Clupea harengus. The sHsp genes were classified into three groups based on phylogenetic analysis. All the sHsps were shown to have higher aliphatic index values, which is an indication that these proteins are more thermally stable. The hydrophilic nature of sHsps was deduced from the fact that all fish species had negative GRAVY scores. In all of the representative fish species, sHsp genes were assigned to distinct chromosomes in an inconsistent and unequal manner. Segmental duplications are the main events that have contributed to the expansion of the sHsp genes in all species. We were also able to determine the selective pressure that was placed on particular codons and discovered several significant coding sites within the coding region of sHsps. Eventually, diversifying positive selection was found to be connected with evolutionary changes in sHsp proteins, which showed that gene evolution controlled the fish adaption event in response to environmental conditions. Clarification of the links between sHsps and environmental stress in fish will be achieved through rigorous genomic comparison, which will also yield substantial new insights.

Comparative analysis of gut microbiome of mangrove brachyuran crabs revealed patterns of phylosymbiosis and codiversification

The acquisition of microbial symbionts enables animals to rapidly adapt to and exploit novel ecological niches, thus significantly enhancing the evolutionary fitness and success of their hosts. However, the dynamics of host-microbe interactions and their evolutionary implications remain largely underexplored in marine invertebrates. Crabs of the family Sesarmidae (Crustacea: Brachyura) are dominant inhabitants of mangrove forests and are considered keystone species there. Their rapid diversification, particularly after adopting a plant-feeding lifestyle, is believed to have been facilitated by symbiotic gut microbes, enabling successful colonization of intertidal and terrestrial environments. To investigate the patterns and mechanisms shaping the microbial communities and the role of microbes in the evolution of Sesarmidae, we characterized and compared the gut microbiome compositions across 43 crab species from Sesarmidae and other mangrove-associated families using 16S metabarcoding. We found that the gut microbiome assemblages in crabs are primarily determined by host identity, with a secondary influence from environmental factors such as microhabitat and sampling location, and to a lesser extent influenced by biological factors such as sex and gut region. While patterns of phylosymbiosis (i.e. when microbial community relationships recapitulate the phylogeny of their hosts) were consistently observed in all beta-diversity metrics analysed, the strength of phylosymbiosis varied across crab families. This suggests that the bacterial assemblages in each family were differentially shaped by different degrees of host filtering and/or other evolutionary processes. Notably, Sesarmidae displayed signals of cophylogeny with its core gut bacterial genera, which likely play crucial functional roles in their hosts by providing lignocellulolytic enzymes, essential amino acids, and fatty acids supplementation. Our results support the hypothesis of microbial contribution to herbivory and terrestrialization in mangrove crabs, highlighting the tight association and codiversification of the crab holobiont.

Cryptic diversity and speciation in an endemic copepod crustacean Harpacticella inopinata within Lake Baikal

Ancient lakes are hotspots of species diversity, posing challenges and opportunities for exploration of the dynamics of endemic diversification. Lake Baikal in Siberia, the oldest lake in the world, hosts a particularly rich crustacean fauna, including the largest known species flock of harpacticoid copepods with some 70 species. Here, we focused on exploring the diversity and evolution within a single nominal species, Harpacticella inopinata Sars, 1908, using molecular markers (mitochondrial COI, nuclear ITS1 and 28S rRNA) and a set of qualitative and quantitative morphological traits. Five major mitochondrial lineages were recognized, with model-corrected COI distances of 0.20-0.37. A concordant pattern was seen in the nuclear data set, and qualitative morphological traits also distinguish a part of the lineages. All this suggests the presence of several hitherto unrecognized cryptic taxa within the baikalian H. inopinata, with long independent histories. The abundances, distributions and inferred demographic histories were different among taxa. Two taxa, H. inopinata CE and H. inopinata CW, were widespread on the eastern and western coasts, respectively, and were largely allopatric. Patterns in mitochondrial variation, that is, shallow star-like haplotype networks, suggest these taxa have spread through the lake relatively recently. Three other taxa, H. inopinata RE, RW and RW2, instead were rare and had more localized distributions on either coast, but showed deeper intraspecies genealogies, suggesting older regional presence. The rare taxa were often found in sympatry with the others and occasionally introgressed by mtDNA from the common ones. The mitochondrial divergence between and within the H. inopinata lineages is still unexpectedly deep, suggesting an unusually high molecular rate. The recognition of true systematic diversity in the evaluation and management of ecosystems is important in hotspots, as it is everywhere else, while the translation of the diversity into a formal taxonomy remains a challenge.

Phylogenomics of the leafhopper genus Neoaliturus Distant, 1918 (Hemiptera: Cicadellidae: Deltocephalinae) reveals genetically divergent lineages in the invasive beet leafhopper

Phylogenomic analysis based on nucleotide sequences of 398 nuclear gene loci for 67 representatives of the leafhopper genus Neoaliturus yielded well-resolved estimates of relationships among species of the genus. Subgenus Neoaliturus (Neoaliturus) is consistently paraphyletic with respect to Neoaliturus (Circulifer). The analysis revealed the presence of at least ten genetically divergent clades among specimens consistent with the previous morphology-based definition of the leafhopper genus "Circulifer" which includes three previously recognized "species complexes." Specimens of the American beet leafhopper, N. tenellus (Baker), collected from the southwestern USA consistently group with one of these clades, comprising specimens from the eastern Mediterranean. Some of the remaining lineages are consistent with ecological differences previously observed among eastern Mediterranean populations and suggest that N. tenellus, as previously defined, comprises multiple monophyletic species, distinguishable by slight morphological differences.

Evolutionary history of arbuscular mycorrhizal fungi and genomic signatures of obligate symbiosis

BACKGROUND

The colonization of land and the diversification of terrestrial plants is intimately linked to the evolutionary history of their symbiotic fungal partners. Extant representatives of these fungal lineages include mutualistic plant symbionts, the arbuscular mycorrhizal (AM) fungi in Glomeromycota and fine root endophytes in Endogonales (Mucoromycota), as well as fungi with saprotrophic, pathogenic and endophytic lifestyles. These fungal groups separate into three monophyletic lineages but their evolutionary relationships remain enigmatic confounding ancestral reconstructions. Their taxonomic ranks are currently fluid.

RESULTS

In this study, we recognize these three monophyletic linages as phyla, and use a balanced taxon sampling and broad taxonomic representation for phylogenomic analysis that rejects a hard polytomy and resolves Glomeromycota as sister to a clade composed of Mucoromycota and Mortierellomycota. Low copy numbers of genes associated with plant cell wall degradation could not be assigned to the transition to a plant symbiotic lifestyle but appears to be an ancestral phylogenetic signal. Both plant symbiotic lineages, Glomeromycota and Endogonales, lack numerous thiamine metabolism genes but the lack of fatty acid synthesis genes is specific to AM fungi. Many genes previously thought to be missing specifically in Glomeromycota are either missing in all analyzed phyla, or in some cases, are actually present in some of the analyzed AM fungal lineages, e.g. the high affinity phosphorus transporter Pho89.

CONCLUSION

Based on a broad taxon sampling of fungal genomes we present a well-supported phylogeny for AM fungi and their sister lineages. We show that among these lineages, two independent evolutionary transitions to mutualistic plant symbiosis happened in a genomic background profoundly different from that known from the emergence of ectomycorrhizal fungi in Dikarya. These results call for further reevaluation of genomic signatures associated with plant symbiosis.

The complete plastome of Acalypha australis (Euphorbiaceae) and its phylogenetic analysis

Acalypha australis L. 1753 is a potherb popular among Asian populations and is a traditional herbal medicine. In the current study, the overall genetic diversity of A. australis still needs to be better. Here, we assembled and characterized the complete plastome of A. australis. The plastome is 168,885 bp in length with a large single-copy (LSC) of 94,576 bp, a small single-copy (SSC) of 19,715 bp, and two copies of inverted repeat region (IRs) of 27,297 bp each. The overall GC content is 34.9%. The plastome contains 127 genes, including 83 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenomic analysis of the representative species of Euphorbiaceae showed that A. australis and A. hispida formed a monophyletic sister clade. The results of this study will support further research on the evolution and conservation of the Euphorbiaceae species; they will benefit pharmaceutical applications and ornamentation of the medicinal plant A. australis.

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

Characterization of the Mitogenome of the Genus Dendrocerus Ratzeburg (Hymenoptera: Megaspilidae) with the Specific Designed Primers

In Hymenoptera, the monophyly of Evaniomorpha has been the focus of debate among different scholars. In this study, we sequenced two mitochondrial genomes of Dendrocerus (Hymenoptera: Megaspilidae) to analyze the mitochondrial genomic features of Dendrocerus and provide new molecular data for phylogenetic studies of Evaniomorpha. The mitogenome sizes of D. bellus and D. anisodontus were 15,445 bp and 15,373 bp, respectively, with the trnG of D. bellus missing. The nucleotide composition was significantly biased toward adenine and thymine, with A + T contents of 81.2% (D. bellus) and 82.4% (D. anisodontus). Using Ceraphron sp. (Ceraphronidae) as reference, the Ka/Ks values of NAD4L and NAD6 in D. anisodontus were both greater than one, indicating that non-synonymous mutations are favored by Darwinian selection, which is rare in other hymenopteran species. Compared with Ceraphon sp. gene order, nine operations were identified in D. anisodontus, including four reversals, four TDRLs (tandem duplication random losses) and one transposition, or four reversals and five TDRLs. Phylogenetic analysis of 40 mitochondrial genomes showed that Evaniomorpha was not a monophyletic group, which was also supported by the PBD values. Ceraphronoidea is a monophyletic group and is a sister to Aulacidae + Gasteruptiidae. Based on the conserved region of the newly sequenced mitochondrial genomes, a pair of specific primers MegaF/MegaR was designed for sequencing the COX1 genes in Megaspilidae and a 60% rate of success was achieved in the genus Dendrocerus.

The complete chloroplast genome and phylogenetic analysis of Cyperus brevifolius (Rottb.) Hassk. 1844 (Cyperaceae)

The medicinal plant Cyperus brevifolius (Rottb.) Hassk. 1844 has a long history of use in traditional Chinese medicine. In this study, we determined and systematically analyzed the complete chloroplast (cp) genome of C. brevifolius. The genome is 183,717 bp in length with a GC content of 33.24%. It comprises four distinct regions: a large-single copy (LSC) region of 101,190 bp, a small-single copy (SSC) region of 10,366 bp, and two inverted repeat (IR) regions of 36,079 bp each. A total of 137 genes are present in the genome including 89 protein-coding genes, 40 tRNA genes, and eight rRNA genes. Phylogenetic analysis reveals that C. brevifolius belongs to the Cyperus genus. This newly sequenced cp genome provides valuable insights for future genetic and genomic studies on Cyperus.