Parallelization of MAFFT for large-scale multiple sequence alignments

Description of three new species of Gyrodactylus von Nordmann, 1832 (Monogenea: Gyrodactylidae) on bitterling fishes (Acheilognathinae) from China

Three new species of Gyrodactylus are described from three species of bitterling in Donghu Lake, China: Gyrodactylus ocellorhodei n. sp. from Rhodeus ocellatus; G. sinenorhodei n. sp. from Rhodeus sinensis; and G. acheilorhodei n. sp. from Acheilognathus macropterus. All the three new species showed similar opisthaptor morphology, especially the marginal hooks: all had a slender and perpendicular sickle shaft, and flat sickle base with distinct heel and inner arch which was different from the G. rhodei-group species parasitic on bitterling. Multivariate analyses based on hamulus and marginal hooks suggested that these three new species cannot be completely distinguished, despite some morphology divergence observed in certain less reliable morphometric features, such as hamulus root length, ventral bar total length and process shape. These three new species shared an identical 18S ribosomal RNA gene sequence, while the variation in the Internal Transcribed Spacers (ITS1-ITS2) sequence among them (8.4-11.2%, K2P) far exceeded the 1% ITS sequence difference that had been suggested as a threshold for species delimitation of Gyrodactylus. Phylogenetic analysis based on ITS1-ITS2 showed that all these sequenced Gyrodactylus spp. parasitic on the subfamily Acheilognathinae host formed a monophyletic group. However, a clear differentiation (18.9-20.9%, K2P of ITS1-ITS2) could be found between the subgroup from China (G. ocellorhodei n. sp., G. sinenorhodei n. sp. and G. acheilorhodei n. sp.) and that from Europe (G. rhodei).

The virome of bubaline (Bubalus bubalis) tonsils reveals an unreported bubaline polyomavirus

Water buffalo (Bubalus bubalis) farming is increasing in many regions of the world due to the species' ability to thrive in environments where bovine cattle would struggle. Despite water buffaloes being known for their resistance to diseases, there is a lack of data about the diversity of the microbiome of the species. In this study, we examined the virome diversity in palatine tonsils collected from animals from the island of Marajó, northern Pará state, Brazil, which harbors the largest bubaline flock in the country. Tonsil fragments from 60 clinically healthy bubalines were randomly selected from a sample of 293 animals. The samples were purified, extracted, and randomly amplified with phi29 DNA polymerase. After amplification, the products were purified and sequenced. Circular DNA viruses were predominant in the tonsils' virome. Sequences of genome segments representative of members of the genera Alphapolyomavirus (including a previously unreported bubaline polyomavirus genome) and Gemycircularvirus were identified, along with other not yet classified circular virus genomes. As the animals were clinically healthy at the time of sampling, such viruses likely constitute part of the normal tonsillar virome of water buffaloes inhabiting the Ilha do Marajó biome.

Identification and characterization of a new species of Taxus - Taxus qinlingensis by multiple taxonomic methods

BACKGROUND

The taxonomy of Taxus Linn. remains controversial due to its continuous phenotypic variation and unstable topology, thus adversely affecting the formulation of scientific conservation strategies for this genus. Recently, a new ecotype, known as Qinling type, is mainly distributed in the Qinling Mountains and belongs to a monophyletic group. Here, we employed multiple methods including leaf phenotype comparison (leaf shapes and microstructure), DNA barcoding identification (ITS + trnL-trnF + rbcL), and niche analysis to ascertain the taxonomic status of the Qinling type.

RESULTS

Multiple comparisons revealed significant differences in the morphological characters (length, width, and length/width ratio) among the Qinling type and other Taxus species. Leaf anatomical analysis indicated that only the Qinling type and T. cuspidata had no papilla under the midvein or tannins in the epicuticle. Phylogenetic analysis of Taxus indicated that the Qinling type belonged to a monophyletic group. Moreover, the Qinling type had formed a relatively independent niche, it was mainly distributed around the Qinling Mountains, Ta-pa Mountains, and Taihang Mountains, situated at an elevation below 1500 m.

CONCLUSIONS

Four characters, namely leaf curvature, margin taper, papillation on midvein, and edges were put forward as primary indexes for distinguishing Taxus species. The ecotype Qingling type represented an independent evolutionary lineage and formed a unique ecological niche. Therefore, we suggested that the Qingling type should be treated as a novel species and named it Taxus qinlingensis Y. F. Wen & X. T. Wu, sp. nov.

Diversification of the Rho transcription termination factor in bacteria

Correct termination of transcription is essential for gene expression. In bacteria, factor-dependent termination relies on the Rho factor, that classically has three conserved domains. Some bacteria also have a functional insertion region. However, the variation in Rho structure among bacteria has not been analyzed in detail. This study determines the distribution, sequence conservation, and predicted features of Rho factors with diverse domain architectures by analyzing 2730 bacterial genomes. About half (49.8%) of the species analyzed have the typical Escherichia coli like Rho while most of the other species (39.8%) have diverse, atypical forms of Rho. Besides conservation of the main domains, we describe a duplicated RNA-binding domain present in specific species and novel variations in the bicyclomycin binding pocket. The additional regions observed in Rho proteins exhibit remarkable diversity. Commonly, however, they have exceptional amino acid compositions and are predicted to be intrinsically disordered, to undergo phase separation, or have prion-like behavior. Phase separation has recently been shown to play roles in Rho function and bacterial fitness during harsh conditions in one species and this study suggests a more widespread role. In conclusion, diverse atypical Rho factors are broadly distributed among bacteria, suggesting additional cellular roles.

Viperin immunity evolved across the tree of life through serial innovations on a conserved scaffold

Evolutionary arms races between cells and viruses drive the rapid diversification of antiviral genes in diverse life forms. Recent discoveries have revealed the existence of immune genes that are shared between prokaryotes and eukaryotes and show molecular and mechanistic similarities in their response to viruses. However, the evolutionary dynamics underlying the conservation and adaptation of these antiviral genes remain mostly unexplored. Here, we show that viperins constitute a highly conserved family of immune genes across diverse prokaryotes and eukaryotes and identify mechanisms by which they diversified in eukaryotes. Our findings indicate that viperins are enriched in Asgard archaea and widely distributed in all major eukaryotic clades, suggesting their presence in the last eukaryotic common ancestor and their acquisition in eukaryotes from an archaeal lineage. We show that viperins maintain their immune function by producing antiviral nucleotide analogues and demonstrate that eukaryotic viperins diversified through serial innovations on the viperin gene, such as the emergence and selection of substrate specificity towards pyrimidine nucleotides, and through partnerships with genes maintained through genetic linkage, notably with nucleotide kinases. These findings unveil biochemical and genomic transitions underlying the adaptation of immune genes shared by prokaryotes and eukaryotes. Our study paves the way for further understanding of the conservation of immunity across domains of life.

RNA-guided RNA silencing by an Asgard archaeal Argonaute

Argonaute proteins are the central effectors of RNA-guided RNA silencing pathways in eukaryotes, playing crucial roles in gene repression and defense against viruses and transposons. Eukaryotic Argonautes are subdivided into two clades: AGOs generally facilitate miRNA- or siRNA-mediated silencing, while PIWIs generally facilitate piRNA-mediated silencing. It is currently unclear when and how Argonaute-based RNA silencing mechanisms arose and diverged during the emergence and early evolution of eukaryotes. Here, we show that in Asgard archaea, the closest prokaryotic relatives of eukaryotes, an evolutionary expansion of Argonaute proteins took place. In particular, a deep-branching PIWI protein (HrAgo1) encoded by the genome of the Lokiarchaeon 'Candidatus Harpocratesius repetitus' shares a common origin with eukaryotic PIWI proteins. Contrasting known prokaryotic Argonautes that use single-stranded DNA as guides and/or targets, HrAgo1 mediates RNA-guided RNA cleavage, and facilitates gene silencing when expressed in human cells and supplied with miRNA precursors. A cryo-EM structure of HrAgo1, combined with quantitative single-molecule experiments, reveals that the protein displays structural features and target-binding modes that are a mix of those of eukaryotic AGO and PIWI proteins. Thus, this deep-branching archaeal PIWI may have retained an ancestral molecular architecture that preceded the functional and mechanistic divergence of eukaryotic AGOs and PIWIs.

Novel Qnr Families as Conserved and Intrinsic Quinolone Resistance Determinants in Aeromonas spp

The environment has been identified as an origin, reservoir, and transmission route of antibiotic resistance genes (ARGs). Among diverse environments, freshwater environments have been recognized as pivotal in the transmission of ARGs between opportunistic pathogens and autochthonous bacteria such as Aeromonas spp. In this study, five environmental strains of Aeromonas spp. exhibiting multidrug resistance (MDR) were selected for whole-genome sequencing to ascertain their taxonomic assignment at the species-level and to delineate their ARG repertoires. Analyses of their genomes revealed the presence of one protein almost identical to AhQnr (A. hydrophila Qnr protein) and four novel proteins similar to AhQnr. To scrutinize the classification and taxonomic distribution of these proteins, all Aeromonas genomes deposited in the NCBI RefSeq genome database (1,222 genomes) were investigated. This revealed that these Aeromonas Qnr (AQnr) proteins are conserved intrinsic resistance determinants of the genus, exhibiting species-specific diversity. Additionally, structure prediction and analysis of contribution to quinolone resistance by AQnr proteins of the isolates, confirmed their functionality as quinolone resistance determinants. Given the origin of mobile qnr genes from aquatic bacteria and the crucial role of Aeromonas spp. in ARG dissemination in aquatic environments, a thorough understanding and strict surveillance of AQnr families prior to the clinical emergence are imperative. In this study, using comparative genome analyses and functional characterization of AQnr proteins in the genus Aeromonas, novel Aeromonas ARGs requiring surveillance has suggested.

Data-Driven Engineering of Phages with Tunable Capsule Tropism for Klebsiella pneumoniae

Klebsiella pneumoniae, a major clinical pathogen known for causing severe infections, is attracting heightened attention due to its escalating antibiotic resistance. Phages are emerging as a promising alternative to antibiotics; however, their specificity to particular hosts often restricts their use. In this study, a collection of 114 phages is obtained and subjected to analysis against 238 clinical K. pneumoniae strains, revealing a spectrum of lytic behaviors. A correlation between putative tail protein clusters and lysis patterns leads to the discovery of six receptor-binding protein (RBP) clusters that determine host capsule tropism. Significantly, RBPs with cross-capsular lysis capabilities are identified. The newly-identified RBPs provide a toolbox for customizing phages to target diverse capsular types. Building on the toolbox, the engineered phages with altered RBPs successfully shifted and broadened their host capsule tropism, setting the stage for tunable phage that offer a precise and flexible solution to combat K. pneumoniae infections.

Impact of whole-genome duplications on structural variant evolution in Cochlearia

Polyploidy, the result of whole-genome duplication (WGD), is a major driver of eukaryote evolution. Yet WGDs are hugely disruptive mutations, and we still lack a clear understanding of their fitness consequences. Here, we study whether WGDs result in greater diversity of genomic structural variants (SVs) and how they influence evolutionary dynamics in a plant genus, Cochlearia (Brassicaceae). By using long-read sequencing and a graph-based pangenome, we find both negative and positive interactions between WGDs and SVs. Masking of recessive mutations due to WGDs leads to a progressive accumulation of deleterious SVs across four ploidal levels (from diploids to octoploids), likely reducing the adaptive potential of polyploid populations. However, we also discover putative benefits arising from SV accumulation, as more ploidy-specific SVs harbor signals of local adaptation in polyploids than in diploids. Together, our results suggest that SVs play diverse and contrasting roles in the evolutionary trajectories of young polyploids.

The complete chloroplast genome of the first registered Paeonia Itoh hybrid cv. Hexie in China

The first registered Paeonia Itoh hybrid cv. Hexie in China is a naturally occurring intersectional hybrid of Sect. Paeonia and Sect. Moutan. In this study, we sequenced, assembled, and analyzed the complete chloroplast genome of Paeonia Itoh hybrid cv. Hexie. The result showed that the chloroplast genome of Hexie, with a typical circular tetrad structure, is 152,958 bp in length, comprising a large single copy (LSC) region of 84,613 bp, a small single copy (SSC) region of 17,051 bp, and two reverse complementary sequences (IRs) of 25,647 bp. The chloroplast genome encoded 116 genes, including 80 protein-coding genes, 32 tRNA genes, and 4 rRNA genes. Phylogenetic analysis inferred from the shared protein-coding genes showed that the Paeonia Itoh hybrid cv. Hexie had the closest phylogenetic relationship with P. suffruticosa, followed by P. ostii, indicating that P. suffruticosa was its maternal parent. This study provides a molecular resource for phylogenetic and maternal parent studies of Paeonia Itoh hybrid, contributing to a basis for Paeonia Itoh hybrid breeding strategies in the future.

Characterization of the family-level Borreliaceae pan-genome and development of an episomal typing protocol

Background

The Borreliaceae family includes many obligate parasitic bacterial species which are etiologically associated with a myriad of zoonotic borrelioses including Lyme disease and vector-borne relapsing fevers. Infections by the Borreliaceae are difficult to detect by both direct and indirect methods, often leading to delayed and missed diagnoses. Efforts to improve diagnoses center around the development of molecular diagnostics (MDx), but due to deep tissue sequestration of the causative spirochaetes and the lack of persistent bacteremias, even MDx assays suffer from a lack of sensitivity. Additionally, the highly extensive genomic heterogeneity among isolates, even within the same species, contributes to the lack of assay sensitivity as single target assays cannot provide universal coverage. This within-species heterogeneity is partly due to differences in replicon repertoires and genomic structures that have likely arisen to support the complex Borreliaceae lifecycle in which these parasites have to survive in multiple hosts each with unique immune responses.

Results

We constructed a Borreliaceae family-level pangenome and characterized the phylogenetic relationships among the constituent taxa which supports the recent taxonomy of splitting the family into at least two genera. Gene content pro les were created for the majority of the Borreliaceae replicons, providing for the first time their unambiguous molecular typing.

Conclusion

Our characterization of the Borreliaceae pan-genome supports the splitting of the former Borrelia genus into two genera and provides for the phylogenetic placement of several non-species designated isolates. Mining this family-level pangenome will enable precision diagnostics corresponding to gene content-driven clinical outcomes while also providing targets for interventions.

Transcriptomics and metabolomics analyses provide insights into resistance genes of tree ferns

As ancient organisms, tree ferns play a crucial role as an evolutionary bridge between lower and higher plant species, providing various utilitarian benefits. However, they face challenges such as overexploitation, climate change, adverse environmental conditions, and insect pests, resulting in conservation concerns. In this study, we provide an overview of metabolic and transcriptomic resources of leaves in two typical tree ferns, A. spinulosa and A. metteniana, and explore the resistance genes for the first time. The landscape of metabolome showed that the compound skimmin may hold medicinal significance. A total of 111 differentially accumulated metabolites (DAMs) were detected, with pathway enrichment analysis highlighting 14 significantly enriched pathways, including 2-oxocarboxylic acid metabolism possibly associated with environmental adaptations. A total of 14,639 differentially expressed genes (DEGs) were found, among which 606 were resistance (R) genes. We identified BAM1 as a significantly differentially expressed R gene, which is one of the core genes within the R gene interaction network. Both the maximum-likelihood phylogenetic tree and the PPI network revealed a close relationship between BAM1, FLS2, and TMK. Moreover, BAM1 showed a significant positive correlation with neochlorogenic acid and kaempferol-7-O-glucoside. These metabolites, known for their antioxidant and anti-inflammatory properties, likely play a crucial role in the defense response of tree ferns. This research provides valuable insights into the metabolic and transcriptomic differences between A. spinulosa and A. metteniana, enhancing our understanding of resistance genes in tree ferns.

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.

PHDtools: A platform for pathogen detection and multi-dimensional genetic signatures decoding to realize pathogen genomics data analyses online

OBJECTIVES

Facing the emerging diseases, rapid identification of the pathogen and multi-dimensional characterization of the genomic features at both isolate-level and population-level through high-throughput sequencing data can provide invaluable information to guide the development of antiviral agents and strategies. However, a user-friendly program is in urgent need for clinical laboratories without bioinformatics background to decode the complex big genomics data.

METHODS

In this study, we developed an interactive online platform named PHDtools with a total of 15 functions to analyze metagenomics data to identify the potential pathogen and decode multi-dimensional genetic signatures including intra-/inter-host variations and lineage-level variations. The platform was applied to analyze the meta-genomic data of the samples collected from the 172 imported COVID-19 cases.

RESULTS

According to the analytical results of mNGS, 27 patients were found to have the co-infections of SARS-CoV-2 with either influenza virus (n = 9) or human picobirnavirus (n = 19). Enough coverages of all the assembled SARS-CoV-2 genomes provided the sub-lineages of Omicron variant, and the number of mutations in the non-structural genes and M gene was increased, as well as the intra-host variations occurred in E and M gene were under positive selection (Ka/Ks > 1). These findings of increased or changed mutations in the SARS-CoV-2 genome characterized the current adaptive evolution patterns of Omicron sub-lineages, and revealed the evolution speed of these sub-lineages might increase.

CONCLUSIONS

Consequently, the application of PHDtools has proved that this platform is accurate, user-friendly and convenient for clinical users who are deficient in bioinformatics, and the clinical monitor of SARS-CoV-2 genomes by PHDtools also highlighted the potential evolution features of current SARS-CoV-2 and indicated that the development of anti-SARS-CoV-2 agents and new-designed vaccines should incorporate the gene variations other than S gene.

AMRViz enables seamless genomics analysis and visualization of antimicrobial resistance

We have developed AMRViz, a toolkit for analyzing, visualizing, and managing bacterial genomics samples. The toolkit is bundled with the current best practice analysis pipeline allowing researchers to perform comprehensive analysis of a collection of samples directly from raw sequencing data with a single command line. The analysis results in a report showing the genome structure, genome annotations, antibiotic resistance and virulence profile for each sample. The pan-genome of all samples of the collection is analyzed to identify core- and accessory-genes. Phylogenies of the whole genome as well as all gene clusters are also generated. The toolkit provides a web-based visualization dashboard allowing researchers to interactively examine various aspects of the analysis results. Availability: AMRViz is implemented in Python and NodeJS, and is publicly available under open source MIT license at https://github.com/amromics/amrviz .

Strain-specific evolution and host-specific regulation of transposable elements in the model plant symbiont Rhizophagus irregularis

Transposable elements (TEs) are repetitive DNA that can create genome structure and regulation variability. The genome of Rhizophagus irregularis, a widely studied arbuscular mycorrhizal fungus (AMF), comprises ∼50% repetitive sequences that include TEs. Despite their abundance, two-thirds of TEs remain unclassified, and their regulation among AMF life stages remains unknown. Here, we aimed to improve our understanding of TE diversity and regulation in this model species by curating repeat datasets obtained from chromosome-level assemblies and by investigating their expression across multiple conditions. Our analyses uncovered new TE superfamilies and families in this model symbiont and revealed significant differences in how these sequences evolve both within and between R. irregularis strains. With this curated TE annotation, we also found that the number of upregulated TE families in colonized roots is 4 times higher than in the extraradical mycelium, and their overall expression differs depending on the plant host. This work provides a fine-scale view of TE diversity and evolution in model plant symbionts and highlights their transcriptional dynamism and specificity during host-microbe interactions. We also provide Hidden Markov Model profiles of TE domains for future manual curation of uncharacterized sequences (https://github.com/jordana-olive/TE-manual-curation/tree/main).

Decoding huge phage diversity: a taxonomic classification of Lak megaphages

High-throughput sequencing for uncultivated viruses has accelerated the understanding of global viral diversity and uncovered viral genomes substantially larger than any that have so far been cultured. Notably, the Lak phages are an enigmatic group of viruses that present some of the largest known phage genomes identified in human and animal microbiomes, and are dissimilar to any cultivated viruses. Despite the wealth of viral diversity that exists within sequencing datasets, uncultivated viruses have rarely been used for taxonomic classification. We investigated the evolutionary relationships of 23 Lak phages and propose a taxonomy for their classification. Predicted protein analysis revealed the Lak phages formed a deeply branching monophyletic clade within the class Caudoviricetes which contained no other phage genomes. One of the interesting features of this clade is that all current members are characterised by an alternative genetic code. We propose the Lak phages belong to a new order, the 'Grandevirales'. Protein and nucleotide-based analyses support the creation of two families, three sub-families, and four genera within the order 'Grandevirales'. We anticipate that the proposed taxonomy of Lak megaphages will simplify the future classification of related viral genomes as they are uncovered. Continued efforts to classify divergent viruses are crucial to aid common analyses of viral genomes and metagenomes.

The variation and evolution of complete human centromeres

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

BS, Raghuraman V, Rajyaguru U, Llamera KE, Andrew L, Anderson AS, Hovius JW, Liberator PA, Simon R, Hao L. Development of a sequence-based in silico OspA typing method for Borrelia burgdorferi sensu lato

Lyme disease (LD), caused by spirochete bacteria of the genus Borrelia burgdorferi sensu lato, remains the most common vector-borne disease in the northern hemisphere. Borrelia outer surface protein A (OspA) is an integral surface protein expressed during the tick cycle, and a validated vaccine target. There are at least 20 recognized Borrelia genospecies, that vary in OspA serotype. This study presents a new in silico sequence-based method for OspA typing using next-generation sequence data. Using a compiled database of over 400 Borrelia genomes encompassing the 4 most common disease-causing genospecies, we characterized OspA diversity in a manner that can accommodate existing and new OspA types and then defined boundaries for classification and assignment of OspA types based on the sequence similarity. To accommodate potential novel OspA types, we have developed a new nomenclature: OspA in silico type (IST). Beyond the ISTs that corresponded to existing OspA serotypes 1-8, we identified nine additional ISTs that cover new OspA variants in B. bavariensis (IST9-10), B. garinii (IST11-12), and other Borrelia genospecies (IST13-17). The IST typing scheme and associated OspA variants are available as part of the PubMLST Borrelia spp. database. Compared to traditional OspA serotyping methods, this new computational pipeline provides a more comprehensive and broadly applicable approach for characterization of OspA type and Borrelia genospecies to support vaccine development.

Peloplasma aerotolerans gen. nov., sp. nov., a Novel Anaerobic Free-Living Mollicute Isolated from a Terrestrial Mud Volcano

A novel aerotolerant anaerobic bacterium (strain M4AhT) was isolated from a terrestrial mud volcano (Taman Peninsula, Russia). Cells were small, cell-wall-less, non-motile cocci, 0.32-0.65 μm in diameter. The isolate was a mesophilic, neutrophilic chemoorganoheterotroph, growing on carbohydrates (D-glucose, D-trehalose, D-ribose, D-mannose, D-xylose, D-maltose, D-lactose, D-cellobiose, D-galactose, D-fructose, and D-sucrose), proteinaceous compounds (yeast extract, tryptone), and pyruvate. Strain M4AhT tolerated 2% oxygen in the gas phase, was catalase-positive, and showed sustainable growth under microaerobic conditions. The dominant cellular fatty acids of strain M4AhT were C16:0 and C18:0. The G+C content of the genomic DNA was 32.42%. The closest phylogenetic relative of strain M4AhT was Mariniplasma anaerobium from the family Acholeplasmataceae (order Acholeplasmatales, class Mollicutes). Based on the polyphasic characterization of the isolate, strain M4AhT is considered to represent a novel species of a new genus, for which the name Peloplasma aerotolerans gen. nov., sp. nov. is proposed. The type strain of Peloplasma aerotolerans is M4AhT (=DSM 112561T = VKM B-3485T = UQM 41475T). This is the first representative of the order Acholeplasmatales, isolated from a mud volcano.

Systematic revision of species of Atractilina and Spiropes hyperparasitic on Meliolales (Ascomycota) in the tropics

Atractilina Dearn. & Barthol. and Spiropes Cif. are genera of asexual fungi that comprise species mainly hyperparasitic on black mildews (Meliolales, Ascomycota). Although a common group of anamorphic fungi, they have been described up to now only by morphology and their systematic position is unknown. The present study provides a morphological treatise of all known species of Atractilina and Spiropes hyperparasitic on Meliolales, including insights into their systematic position, based on DNA sequences generated here for the first time. The study was conducted, based on 33 herbarium specimens and 23 specimens recently collected in Benin and Panama. The obtained DNA sequence data (28S rDNA and ITS rDNA) of A.parasitica and of two species of Spiropes show systematic placements in the Dothideomycetes and Leotiomycetes, respectively. The sequence data of the two Spiropes spp. do not group together. Moreover, the anamorph-teleomorph connection between Atractilinaparasitica and Malacariameliolicola, a pseudothecioid fungus, is confirmed. Three species in the genus Spiropes are proposed as new to science, namely S.angylocalycis, S.carpolobiae and S.croissantiformis. Four species are reported for Benin for the first time, three species for Panama and one species for mainland America. Atractilina and Spiropes are currently two genera with highly heterogeneous species and they might have to be split in the future, once the taxonomic concepts are validated by morphology and molecular sequence data.

Annotation of 2,507 Saccharomyces cerevisiae genomes

Saccharomyces cerevisiae (baker's yeast, budding yeast) is one of the most important model organisms for biological research and is a crucial microorganism in industry. Currently, a huge number of Saccharomyces cerevisiae genome sequences are available at the public domain. However, these genomes are distributed at different websites and a large number of them are released without annotation information. To provide one complete annotated genome data resource, we collected 2,507 Saccharomyces cerevisiae genome assemblies and re-annotated 2,506 assemblies using a custom annotation pipeline, producing a total of 15,407,164 protein-coding gene models. With a custom pipeline, all these gene sequences were clustered into families. A total of 1,506 single-copy genes were selected as marker genes, which were then used to evaluate the genome completeness and base qualities of all assemblies. Pangenomic analyses were performed based on a selected subset of 847 medium-high-quality genomes. Statistical comparisons revealed a number of gene families showing copy number variations among different organism sources. To the authors' knowledge, this study represents the largest genome annotation project of S. cerevisiae so far, providing rich genomic resources for the future studies of the model organism S. cerevisiae and its relatives.IMPORTANCESaccharomyces cerevisiae (baker's yeast, budding yeast) is one of the most important model organisms for biological research and is a crucial microorganism in industry. Though a huge number of Saccharomyces cerevisiae genome sequences are available at the public domain, these genomes are distributed at different websites and most are released without annotation, hindering the efficient reuse of these genome resources. Here, we collected 2,507 genomes for Saccharomyces cerevisiae, performed genome annotation, and evaluated the genome qualities. All the obtained data have been deposited at public repositories and are freely accessible to the community. This study represents the largest genome annotation project of S. cerevisiae so far, providing one complete annotated genome data set for S. cerevisiae, an important workhorse for fundamental biology, biotechnology, and industry.

The first complete mitochondrial genome data of the Afghan pika Ochotona rufescens (Lagomorpha, Ochotonidae), near the type locality

The Afghan pika Ochotona rufescens (Gray, 1842) is widely distributed across the mountains of Afghanistan, Iran, Pakistan, and southwestern Turkmenistan, most often at elevations between 2,000 and 3,000 m. Here we present, for the first time, the complete mitochondrial genomes of two specimens of the nominotypical subspecies Ochotona rufescens rufescens, de novo assembled from Illumina short reads of fragmented probe-enriched DNA. The lengths of the circular mitogenomes are 16,408 bp and 16,407 bp, respectively. Both mitogenomes contain 13 protein-coding genes (PCGs), two ribosomal RNAs (16S rRNA and 12S rRNA), 22 transfer RNA genes, and a control region. The gene NAD6 and the tRNA (Gln), tRNA (Ala), tRNA (Asn), tRNA (Cys), tRNA (Tyr), tRNA (Ser), tRNA (Glu), and tRNA (Pro) are encoded on the light strand while the rest are encoded on the heavy strand. The overall nucleotide composition was ∼30% for A, 25% for T, 15% for G, and 29% for C. The mitogenome data are available in the GenBank under the accession numbers ON859136 and ON859137.

Molecular and phylogenetic evidence of parallel expansion of anion channels in plants

Aluminum-activated malate transporters (ALMTs) and slow anion channels (SLACs) are important in various physiological processes in plants, including stomatal regulation, nutrient uptake, and in response to abiotic stress such as aluminum toxicity. To understand their evolutionary history and functional divergence, we conducted phylogenetic and expression analyses of ALMTs and SLACs in green plants. Our findings from phylogenetic studies indicate that ALMTs and SLACs may have originated from green algae and red algae, respectively. The ALMTs of early land plants and charophytes formed a monophyletic clade consisting of three subgroups. A single duplication event of ALMTs was identified in vascular plants and subsequent duplications into six clades occurred in angiosperms, including an identified clade, 1-1. The ALMTs experienced gene number losses in clades 1-1 and 2-1 and expansions in clades 1-2 and 2-2b. Interestingly, the expansion of clade 1-2 was also associated with higher expression levels compared to genes in clades that experienced apparent loss. SLACs first diversified in bryophytes, followed by duplication in vascular plants, giving rise to three distinct clades (I, II, and III), and clade II potentially associated with stomatal control in seed plants. SLACs show losses in clades II and III without substantial expansion in clade I. Additionally, ALMT clade 2-2 and SLAC clade III contain genes specifically expressed in reproductive organs and roots in angiosperms, lycophytes, and mosses, indicating neofunctionalization. In summary, our study demonstrates the evolutionary complexity of ALMTs and SLACs, highlighting their crucial role in the adaptation and diversification of vascular plants.

The genome of Stephania japonica provides insights into the biosynthesis of cepharanthine

Stephania japonica is an early-diverging eudicotyledon plant with high levels of cepharanthine, proven to be effective in curing coronavirus infections. Here, we report a high-quality S. japonica genome. The genome size is 688.52 Mb, and 97.37% sequences anchor to 11 chromosomes. The genome comprises 67.46% repetitive sequences and 21,036 genes. It is closely related to two Ranunculaceae species, which diverged from their common ancestor 55.90-71.02 million years ago (Mya) with a whole-genome duplication 85.59-96.75 Mya. We further reconstruct ancestral karyotype of Ranunculales. Several cepharanthine biosynthesis genes are identified and verified by western blot. Two genes (Sja03G0243 and Sja03G0241) exhibit catalytic activity as shown by liquid chromatography-mass spectrometry. Then, cepharanthine biosynthesis genes, transcription factors, and CYP450 family genes are used to construct a comprehensive network. Finally, we construct an early-diverging eudicotyledonous genome resources (EEGR) database. As the first genome of the Menispermaceae family to be released, this study provides rich resources for genomic studies.

Spatiotemporal genotype replacement of H5N8 avian influenza viruses contributed to H5N1 emergence in 2021/2022 panzootic

Since 2020, clade 2.3.4.4b highly pathogenic avian influenza H5N8 and H5N1 viruses have swept through continents, posing serious threats to the world. Through comprehensive analyses of epidemiological, genetic, and bird migration data, we found that the dominant genotype replacement of the H5N8 viruses in 2020 contributed to the H5N1 outbreak in the 2021/2022 wave. The 2020 outbreak of the H5N8 G1 genotype instead of the G0 genotype produced reassortment opportunities and led to the emergence of a new H5N1 virus with G1's HA and MP genes. Despite extensive reassortments in the 2021/2022 wave, the H5N1 virus retained the HA and MP genes, causing a significant outbreak in Europe and North America. Furtherly, through the wild bird migration flyways investigation, we found that the temporal-spatial coincidence between the outbreak of the H5N8 G1 virus and the bird autumn migration may have expanded the H5 viral spread, which may be one of the main drivers of the emergence of the 2020-2022 H5 panzootic.IMPORTANCESince 2020, highly pathogenic avian influenza (HPAI) H5 subtype variants of clade 2.3.4.4b have spread across continents, posing unprecedented threats globally. However, the factors promoting the genesis and spread of H5 HPAI viruses remain unclear. Here, we found that the spatiotemporal genotype replacement of H5N8 HPAI viruses contributed to the emergence of the H5N1 variant that caused the 2021/2022 panzootic, and the viral evolution in poultry of Egypt and surrounding area and autumn bird migration from the Russia-Kazakhstan region to Europe are important drivers of the emergence of the 2020-2022 H5 panzootic. These findings provide important targets for early warning and could help control the current and future HPAI epidemics.

Endogenous virophages are active and mitigate giant virus infection in the marine protist Cafeteria burkhardae

Endogenous viral elements (EVEs) are common genetic passengers in various protists. Some EVEs represent viral fossils, whereas others are still active. The marine heterotrophic flagellate Cafeteria burkhardae contains several EVE types related to the virophage mavirus, a small DNA virus that parasitizes the lytic giant virus CroV. We hypothesized that endogenous virophages may act as an antiviral defense system in protists, but no protective effect of virophages in wild host populations has been shown so far. Here, we tested the activity of virophage EVEs and studied their impact on giant virus replication. We found that endogenous mavirus-like elements (EMALEs) from globally distributed Cafeteria populations produced infectious virus particles specifically in response to CroV infection. However, reactivation was stochastic, often inefficient, and poorly reproducible. Interestingly, only one of eight EMALE types responded to CroV infection, implying that other EMALEs may be linked to different giant viruses. We isolated and cloned several reactivated virophages and characterized their particles, genomes, and infection dynamics. All tested virophages inhibited the production of CroV during coinfection, thereby preventing lysis of the host cultures in a dose-dependent manner. Comparative genomics of different C. burkhardae strains revealed that inducible EMALEs are common and are not linked to specific geographic locations. We demonstrate that naturally occurring virophage EVEs reactivate upon giant virus infection, thus providing a striking example that eukaryotic EVEs can become active under specific conditions. Moreover, our results support the hypothesis that virophages can act as an adaptive antiviral defense system in protists.

A long-read sequencing strategy with overlapping linkers on adjacent fragments (OLAF-Seq) for targeted resequencing and enrichment

In this report, we present OLAF-Seq, a novel strategy to construct a long-read sequencing library such that adjacent fragments are linked with end-terminal duplications. We use the CRISPR-Cas9 nickase enzyme and a pool of multiple sgRNAs to perform non-random fragmentation of targeted long DNA molecules (> 300kb) into smaller library-sized fragments (about 20 kbp) in a manner so as to retain physical linkage information (up to 1000 bp) between adjacent fragments. DNA molecules targeted for fragmentation are preferentially ligated with adaptors for sequencing, so this method can enrich targeted regions while taking advantage of the long-read sequencing platforms. This enables the sequencing of target regions with significantly lower total coverage, and the genome sequence within linker regions provides information for assembly and phasing. We demonstrated the validity and efficacy of the method first using phage and then by sequencing a panel of 100 full-length cancer-related genes (including both exons and introns) in the human genome. When the designed linkers contained heterozygous genetic variants, long haplotypes could be established. This sequencing strategy can be readily applied in both PacBio and Oxford Nanopore platforms for both long and short genes with an easy protocol. This economically viable approach is useful for targeted enrichment of hundreds of target genomic regions and where long no-gap contigs need deep sequencing.

Comparative evolutionary analyses of peste des petits ruminants virus genetic lineages

Peste des petits ruminants virus (PPRV) causes a highly infectious disease affecting mainly goats and sheep in large parts of Africa, Asia, and the Middle East and has an important impact on the global economy and food security. Full genome sequencing of PPRV strains has proved to be critical to increasing our understanding of PPR epidemiology and to inform the ongoing global efforts for its eradication. However, the number of full PPRV genomes published is still limited and with a heavy bias towards recent samples and genetic Lineage IV (LIV), which is only one of the four existing PPRV lineages. Here, we generated genome sequences for twenty-five recent (2010-6) and seven historical (1972-99) PPRV samples, focusing mainly on Lineage II (LII) in West Africa. This provided the first opportunity to compare the evolutionary pressures and history between the globally dominant PPRV genetic LIV and LII, which is endemic in West Africa. Phylogenomic analysis showed that the relationship between PPRV LII strains was complex and supported the extensive transboundary circulation of the virus within West Africa. In contrast, LIV sequences were clearly separated per region, with strains from West and Central Africa branched as a sister clade to all other LIV sequences, suggesting that this lineage also has an African origin. Estimates of the time to the most recent common ancestor place the divergence of modern LII and LIV strains in the 1960s-80s, suggesting that this period was particularly important for the diversification and spread of PPRV globally. Phylogenetic relationships among historical samples from LI, LII, and LIII and with more recent samples point towards a high genetic diversity for all these lineages in Africa until the 1970s-80s and possible bottleneck events shaping PPRV's evolution during this period. Molecular evolution analyses show that strains belonging to LII and LIV have evolved under different selection pressures. Differences in codon usage and adaptative selection pressures were observed in all viral genes between the two lineages. Our results confirm that comparative genomic analyses can provide new insights into PPRV's evolutionary history and molecular epidemiology. However, PPRV genome sequencing efforts must be ramped up to increase the resolution of such studies for their use in the development of efficient PPR control and surveillance strategies.

Redefining the treponemal history through pre-Columbian genomes from Brazil

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Comparative genomics reveals the adaptation of ammonia-oxidising Thaumarchaeota to arid soils

Thaumarchaeota are predominant in oligotrophic habitats such as deserts and arid soils, but their adaptations to these arid conditions are not well understood. In this study, we assembled 23 Thaumarchaeota genomes from arid and semi-arid soils collected from the Inner Mongolia Steppe and the Qinghai-Tibet Plateau. Using a comparative genomics approach, integrated with 614 Thaumarchaeota genomes from public databases, we identified the traits and evolutionary forces that contribute to their adaptations to aridity. Our results showed that the newly assembled genomes represent an early diverging group within the lineage of ammonia-oxidising Thaumarchaeota. While the genomic functions previously identified in arid soil lineages were conserved across terrestrial, shallow-ocean and deep-ocean lineages, several traits likely contribute to Thaumarchaeota's adaptation to aridity. These include chlorite dismutase, arsenate reductase, V-type ATPase and genes dealing with oxidative stresses. The acquisition and loss of traits at the last common ancestor of arid soil lineages may have facilitated the specialisation of Thaumarchaeota in arid soils. Additionally, the acquisition of unique adaptive traits, such as a urea transporter, Ca2+ :H+ antiporter, mannosyl-3-phosphoglycerate synthase and phosphatase, DNA end-binding protein Ku and phage shock protein A, further distinguishes arid soil Thaumarchaeota. This study provides evidence for the adaptations of Thaumarchaeota to arid soil, enhancing our understanding of the nitrogen and carbon cycling driven by Thaumarchaeota in drylands.

CNCA aligns small annotated genomes

BACKGROUND

To explore the evolutionary history of sequences, a sequence alignment is a first and necessary step, and its quality is crucial. In the context of the study of the proximal origins of SARS-CoV-2 coronavirus, we wanted to construct an alignment of genomes closely related to SARS-CoV-2 using both coding and non-coding sequences. To our knowledge, there is no tool that can be used to construct this type of alignment, which motivated the creation of CNCA.

RESULTS

CNCA is a web tool that aligns annotated genomes from GenBank files. It generates a nucleotide alignment that is then updated based on the protein sequence alignment. The output final nucleotide alignment matches the protein alignment and guarantees no frameshift. CNCA was designed to align closely related small genome sequences up to 50 kb (typically viruses) for which the gene order is conserved.

CONCLUSIONS

CNCA constructs multiple alignments of small genomes by integrating both coding and non-coding sequences. This preserves regions traditionally ignored in conventional back-translation methods, such as non-coding regions.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

High-quality genome of the zoophytophagous stink bug, Nesidiocoris tenuis, informs their food habit adaptation

The zoophytophagous stink bug, Nesidiocoris tenuis, is a promising natural enemy of micro-pests such as whiteflies and thrips. This bug possesses both phytophagous and entomophagous food habits, enabling it to obtain nutrition from both plants and insects. This trait allows us to maintain its population density in agricultural fields by introducing insectary plants, even when the pest prey density is extremely low. However, if the bugs' population becomes too dense, they can sometimes damage crop plants. This dual character seems to arise from the food preferences and chemosensation of this predator. To understand the genomic landscape of N. tenuis, we examined the whole genome sequence of a commercially available Japanese strain. We used long-read sequencing and Hi-C analysis to assemble the genome at the chromosomal level. We then conducted a comparative analysis of the genome with previously reported genomes of phytophagous and hematophagous stink bugs to focus on the genetic factors contributing to this species' herbivorous and carnivorous tendencies. Our findings suggest that the gustatory gene set plays a pivotal role in adapting to food habits, making it a promising target for selective breeding. Furthermore, we identified the whole genomes of microorganisms symbiotic with this species through genomic analysis. We believe that our results shed light on the food habit adaptations of N. tenuis and will accelerate breeding efforts based on new breeding techniques for natural enemy insects, including genomics and genome editing.

Genetic characterization of parvoviruses identified in stray cats in Nigeria

Parvoviruses are a major cause of haemorrhagic gastroenteritis, leukopenia and high mortality in cats and dogs. In this study, the presence and genetic characteristics of parvoviruses circulating among cats in Nigeria are reported. Faecal samples of stray cats from live animal markets in southwestern (Oyo and Osun States) and north-central (Kwara State) Nigeria were screened for the presence of parvoviral DNA using a qPCR. Positive samples were further characterized using a qPCR based on minor groove binder probes. Overall, 85/102 (83.3 %) stray cats tested positive for feline panleukopenia virus (FPV) DNA and one cat was co-infected with canine parvovirus-2 type a. Sequence analysis of the complete capsid region of 15 Nigerian FPV strains revealed that they were up to 99.9 % similar to the American reference strain FPV-b at the nucleotide level, and three of them presented amino acid mutations in key capsid residues. This is the first report of identification and molecular characterization of FPV strains in cats in Nigeria. The high prevalence of the virus emphasizes the need for constant surveillance of the circulation of parvoviruses in Nigeria and underscores the need to deploy an effective vaccination strategy.

Whole-Genome Sequencing Analyses Reveal the Whip-like Tail Formation, Innate Immune Evolution, and DNA Repair Mechanisms of Eupleurogrammus muticus

Smallhead hairtail (Eupleurogrammus muticus) is an important marine economic fish distributed along the northern Indian Ocean and the northwest Pacific coast; however, little is known about the mechanism of its genetic evolution. This study generated the first genome assembly of E. muticus at the chromosomal level using a combination of PacBio SMRT, Illumina Nova-Seq, and Hi-C technologies. The final assembled genome size was 709.27 Mb, with a contig N50 of 25.07 Mb, GC content of 40.81%, heterozygosity rate of 1.18%, and repetitive sequence rate of 35.43%. E. muticus genome contained 21,949 protein-coding genes (97.92% of the genes were functionally annotated) and 24 chromosomes. There were 143 expansion gene families, 708 contraction gene families, and 4888 positively selected genes in the genome. Based on the comparative genomic analyses, we screened several candidate genes and pathways related to whip-like tail formation, innate immunity, and DNA repair in E. muticus. These findings preliminarily reveal some molecular evolutionary mechanisms of E. muticus at the genomic level and provide important reference genomic data for the genetic studies of other trichiurids.

Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene βglu-1 with large variations in white spruce but not Norway spruce

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and βglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of βglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of βglu-1 and Ugt5 genes. We observed very large copy numbers of βglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of βglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

The complete chloroplast genome sequence of Malus × adstringens Zabel 'Hopa' (Rosaceae)

Malus × adstringens Zabel 'Hopa' is an important crabapple cultivar with significant ornamental value. Here, we assembled its complete chloroplast (cp) genome using the next-generation sequencing technology to clarify the phylogenetic relationships in Malus. The total length of the complete chloroplast genome was 160,230 base pairs (bp) with a GC content of 36.50%, consisting of a large single-copy (LSC) region with a sequence length of 88,310 bp, a small single-copy (SSC) region with a sequence length of 19,196 bp, and a pair of inverted repeat (IR) regions of 26,362 bp. The complete chloroplast genome contained 128 genes, namely 84 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. In addition, 73 SSRs were found in the M. 'Hopa' cp genome. The phylogenetic relationship of M. 'Hopa' in Malus is closely related to M. spectabilis (Aiton) Borkh. and then to M. sieversii (Lebed.) M. Roem. Our results demonstrate that it is feasible to resolve the phylogenetic relationships of crabapple cultivars and identify their putative maternal lineages using cp genomic data.

Complete chloroplast genome of Petrocosmea qinlingensis (Gesneriaceae), a protected wild plant in the Qinling mountains

Petrocosmea qinlingensis is a protected wild plant endemic in China, inhabiting low-light limestone cliffs but the complete chloroplast genome has not been reported. In this study, we first sequenced and assembled the complete chloroplast genome of P. qinlingensis. The total size of this genome was 153,865 bp, including a large single-copy (LSC) region (84,737 bp), a small single-copy (SSC) region (18,244 bp), and two inverted repeats (IRs) regions (25,442 bp). This genome encoded 111 uniquegenes, consisted of 77 protein-coding genes, four ribosomal RNA genes, and 30 transfer RNA genes. Phylogenomic analysis based on the chloroplast protein-coding genes and showed that the genus Petrocosmea was the closest relative to Raphiocarpus. Our results will support further phylogeographic, population genetic studies of this species.

Phylogenetically evolutionary analysis provides insights into the genetic diversity and adaptive evolution of porcine deltacoronavirus

BACKGROUND

Porcine deltacoronavirus (PDCoV) is one of the emerging swine enteric coronaviruses (SECoVs), which has been widely prevalent in the North America and Asia. In addition to causing severe diarrhea in piglets, PDCoV also shows the potential to infect diverse host species, including calves, chickens, turkey poults, and humans. However, the clinical pathogenicity and genetic evolution of PDCoV is still not fully understood.

RESULTS

Here, we recorded an outbreak of a novel recombinant PDCoV strain (CHN-HeN06-2022) in a large nursery fattening pig farm. Genomic analysis showed that the CHN-HeN06-2022 strain shared 98.3-98.7% sequence identities with the Chinese and American reference strains. To clarify the evolutionary relationships, phylogenetic analysis was performed using the PDCoV genome sequences available in the GenBank database. Based on genetic distance and geographical distribution, the phylogenetic tree clearly showed that all the PDCoV sequences could be divided into lineage 1 and lineage 2, which were further classified into sublineage 1.1 (Chinese strains), 1.2 (the North American strains), 2.1 (the Southeast Asian strains), and 2.2 (Chinese strains). Corresponding to the evolutionary tree, we found that, compared to lineage 1, lineage 2 strains usually contain a continuous 6-nt deletion in Nsp2 and a 9-nt deletion in Nsp3, respectively. Furthermore, recombination analysis suggested that the CHN-HeN06-2022 occurred segments exchange crossed Nsp2 and Nsp3 region between sublineage 1.1 and sublineage 2.1. Combined with previously reported recombinant strains, the highest recombination frequency occurred in Nsp2, Nsp3, and S gene. Additionally, we identified a total of 14 amino acid sites under positive selection in spike protein, most of which are located in the regions related with the viral attachment, receptor binding, and membrane fusion.

CONCLUSIONS

Taken together, our studies provide novel insights into the genetic diversity and adaptive evolution of PDCoV. It would be helpful to the development of vaccine and potential antiviral agent.

Chromosomal-scale genome assembly and annotation of the land slug (Meghimatium bilineatum)

Meghimatium bilineatum is a notorious pest land slug used as a medicinal resource to treat ailments in China. Although this no-model species is unique in terms of their ecological security and medicinal value, the genome resource of this slug is lacking to date. Here, we used the Illumina, PacBio, and Hi-C sequencing techniques to construct a chromosomal-level genome of M. bilineatum. With the Hi-C correction, the sequencing data from PacBio system generated a 1.61 Gb assembly with a scaffold N50 of 68.08 Mb, and anchored to 25 chromosomes. The estimated assembly completeness at 91.70% was obtained using BUSCO methods. The repeat sequence content in the assembled genome was 72.51%, which mainly comprises 34.08% long interspersed elements. We further identified 18631 protein-coding genes in the assembled genome. A total of 15569 protein-coding genes were successfully annotated. This genome assembly becomes an important resource for studying the ecological adaptation and potential medicinal molecular basis of M. bilineatum.

RAF-like protein kinases mediate a deeply conserved, rapid auxin response

The plant-signaling molecule auxin triggers fast and slow cellular responses across land plants and algae. The nuclear auxin pathway mediates gene expression and controls growth and development in land plants, but this pathway is absent from algal sister groups. Several components of rapid responses have been identified in Arabidopsis, but it is unknown if these are part of a conserved mechanism. We recently identified a fast, proteome-wide phosphorylation response to auxin. Here, we show that this response occurs across 5 land plant and algal species and converges on a core group of shared targets. We found conserved rapid physiological responses to auxin in the same species and identified rapidly accelerated fibrosarcoma (RAF)-like protein kinases as central mediators of auxin-triggered phosphorylation across species. Genetic analysis connects this kinase to both auxin-triggered protein phosphorylation and rapid cellular response, thus identifying an ancient mechanism for fast auxin responses in the green lineage.

Effect of operational parameters on the performance of an anaerobic sequencing batch reactor (AnSBR) treating protein-rich wastewater

Treating protein-rich wastewater using cost-effective and simple-structured single-stage reactors presents several challenges. In this study, we applied an anaerobic sequencing batch reactor (AnSBR) to treat protein-rich wastewater from a slaughterhouse. We focused on identifying the key factors influencing the removal of chemical oxygen demand (COD) and the settling performance of the sludge. The AnSBR achieved a maximum total COD removal of 90%, a protein degradation efficiency exceeding 80%, and a COD to methane conversion efficiency of over 70% at organic loading rates of up to 6.2 g COD L-1 d-1. We found that the variations in both the organic loading rate within the reactor and the hydraulic retention time in the buffer tank had a significant effect on COD removal. The hydraulic retention time in the buffer tank and the reactor, which determined the ammonification efficiencies and the residual carbohydrate concentrations in the reactor liquid, affected the sludge settleability. Furthermore, the genus Clostridium sensu stricto 1, known as protein- and lipids-degraders, was predominant in the reactor. Statistical analysis showed a significant correlation between the core microbiome and ammonification efficiency, highlighting the importance of protein degradation as the governing process in the treatment. Our results will provide valuable insights to optimise the design and operation of AnSBR for efficient treatment of protein-rich wastewater.

Genome-wide identification of the glutamate receptor-like gene family in Vanilla planifolia and their response to Fusarium oxysporum infection

Glutamate receptor-like genes (GLRs) are essential for plant growth and development and for coping with environmental (biological and non-biological) stresses. In this study, 13 GLR members were identified in the Vanilla planifolia genome and attributed to two subgroups (Clade I and Clade III) based on their physical relationships. Cis-acting element analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations indicated the GLR gene regulation's complexity and their functional diversity. Expression analysis revealed a relatively higher and more general expression pattern of Clade III members compared to the Clade I subgroup in tissues. Most GLRs showed significant differences in expression during Fusarium oxysporum infection. This suggested that GLRs play a critical role in the response of V. planifolia to pathogenic infection. These results provide helpful information for further functional research and crop improvement of VpGLRs.

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

Genotypic Characterization of Uropathogenic Escherichia coli from Companion Animals: Predominance of ST372 in Dogs and Human-Related ST73 in Cats

Extraintestinal pathogenic Escherichia coli (ExPEC) account for over 80% and 60% of bacterial urinary tract infections (UTIs) in humans and animals, respectively. As shared uropathogenic E. coli (UPEC) strains have been previously reported among humans and pets, our study aimed to characterize E. coli lineages among UTI isolates from dogs and cats and to assess their overlaps with human UPEC lineages. We analysed 315 non-duplicate E. coli isolates from the UT of dogs (198) and cats (117) collected in central Germany in 2019 and 2020 utilizing whole genome sequencing and in silico methods. Phylogroup B2 (77.8%), dog-associated sequence type (ST) 372 (18.1%), and human-associated ST73 (16.6%), were predominant. Other STs included ST12 (8.6%), ST141 (5.1%), ST127 (4.8%), and ST131 (3.5%). Among these, 58.4% were assigned to the ExPEC group and 51.1% to the UPEC group based on their virulence associated gene (VAG) profile (ExPEC, presence of ≥VAGs: papAH and/or papC, sfa/focG, afaD/draBC, kpsMTII, and iutA; UPEC, additionally cnf1 or hlyD). Extended-spectrum cephalosporin (ESC) resistance mediated by extended-spectrum β-lactamases (ESBL) and AmpC-β-lactamase was identified in 1.9% of the isolates, along with one carbapenemase-producing isolate and one isolate carrying a mcr gene. Low occurrence of ESC-resistant or multidrug-resistant (MDR) isolates (2.9%) in the two most frequently detected STs implies that E. coli isolated from UTIs of companion animals are to a lesser extent associated with resistance, but possess virulence-associated genes enabling efficient UT colonization and carriage. Detection of human-related pandemic lineages suggests interspecies transmission and underscores the importance of monitoring companion animals.

Revealing viral diversity in the Napahai plateau wetland based on metagenomics

We focused on exploring the diversity of viruses in the Napahai plateau wetland, a unique ecosystem located in Yunnan, China. While viruses in marine environments have been extensively studied for their influence on microbial metabolism and biogeochemical cycles, little is known about their composition and function in plateau wetlands. Metagenomic analysis was employed to investigate the viral diversity and biogeochemical impacts in the Napahai wetland. It revealed that the Caudoviricetes and Malgrandaviricetes class level was the most abundant viral category based on phylogenetic analysis. Additionally, a gene-sharing network highlighted the presence of numerous unexplored viruses and demonstrated their unique characteristics and significant variation within the viral community of the Napahai wetland. Furthermore, the study identified the auxiliary metabolic genes (AMGs). AMGs provide phages with additional functions, such as protection against host degradation and involvement in metabolic pathways, such as the pentose phosphate pathway and DNA biosynthesis. The viruses in the Napahai wetland were found to influence carbon, nitrogen, sulfur, and amino acid metabolism, indirectly contributing to biogeochemical cycling through these AMGs. Overall, the research sheds light on the diverse and unique viral communities in the Napahai plateau wetland and emphasizes the significant roles of viruses in microbial ecology. The findings contribute to a deeper understanding of the characteristics and ecological functions of viral communities in plateau wetland ecosystems.

The genomic and epidemiological investigations of enteric viruses of domestic caprine (Capra hircus) revealed the presence of multiple novel viruses related to known strains of humans and ruminant livestock species

IMPORTANCE

Compared with other domestic animals, the virome and viral diversity of small ruminants especially in caprine are less studied even of its zoonotic potential. In this study, the enteric virome of caprine was investigated in detail using next-generation sequencing and reverse transcription PCR techniques. The complete or nearly complete genomes of seven novel viruses were determined which show a close phylogenetic relationship to known human and ruminant viruses. The high similarity between the identified caprine tusavirus (family Parvoviridae) and an unassigned CRESS DNA virus with closely related human strains could indicate the (reverse) zoonotic potential of these viruses. Others, like astroviruses (family Astroviridae), enteroviruses, or novel caripiviruses (named after the term caprine picornavirus) of family Picornaviridae found mostly in multiple co-infections in caprine and ovine, could indicate the cross-species transmission capabilities of these viruses between small ruminants.

Nightmare or delight: Taxonomic circumscription meets reticulate evolution in the phylogenomic era

Phylogenetic studies in the phylogenomics era have demonstrated that reticulate evolution greatly impedes the accuracy of phylogenetic inference, and consequently can obscure taxonomic treatments. However, the systematics community lacks a broadly applicable strategy for taxonomic delimitation in groups characterized by pervasive reticulate evolution. The red-fruit genus, Stranvaesia, provides an ideal model to examine the influence of reticulation on generic circumscription, particularly where hybridization and allopolyploidy dominate the evolutionary history. In this study, we conducted phylogenomic analyses integrating data from hundreds of single-copy nuclear (SCN) genes and plastomes, and interrogated nuclear paralogs to clarify the inter/intra-generic relationship of Stranvaesia and its allies in the framework of Maleae. Analyses of phylogenomic discord and phylogenetic networks showed that allopolyploidization and introgression promoted the origin and diversification of the Stranvaesia clade, a conclusion further bolstered by cytonuclear and gene tree discordance. With a well-inferred phylogenetic backbone, we propose an updated generic delimitation of Stranvaesia and introduce a new genus, Weniomeles. This new genus is distinguished by its purple-black fruits, thorns trunk and/or branches, and a distinctive fruit core anatomy characterized by multilocular separated by a layer of sclereids and a cluster of sclereids at the top of the locules. Through this study, we highlight a broadly-applicable workflow that underscores the significance of reticulate evolution analyses in shaping taxonomic revisions from phylogenomic data.