Comparative Genomics Research Articles

Discovery of Candidatus Nitrosomaritimum as a New Genus of Ammonia-Oxidizing Archaea Widespread in Anoxic Saltmarsh Intertidal Aquifers.

Ammonia-oxidizing archaea (AOA) are widely distributed in marine and terrestrial habitats, contributing significantly to global nitrogen and carbon cycles. However, their genomic diversity, ecological niches, and metabolic potentials in the anoxic intertidal aquifers remain poorly understood. Here, we discovered and named a novel AOA genus, Candidatus Nitrosomaritimum, from the intertidal aquifers of Yancheng Wetland, showing close metagenomic abundance to the previously acknowledged dominant Nitrosopumilus AOA. Further construction of ammonia monooxygenase-based phylogeny demonstrated the widespread distribution of Nitrosomaritimum AOA in global estuarine-coastal niches and marine sediment. Niche differentiation among sublineages of this new genus in anoxic intertidal aquifers is driven by salinity and dissolved oxygen gradients. Comparative genomics revealed that Candidatus Nitrosomaritimum has the genetic capacity to utilize urea and possesses high-affinity phosphate transporter systems (phnCDE) for surviving phosphorus-limited conditions. Additionally, it contains putative nosZ genes encoding nitrous-oxide (N2O) reductase for reducing N2O to nitrogen gas. Furthermore, we gained first genomic insights into the archaeal phylum Hydrothermarchaeota populations residing in intertidal aquifers and revealed their potential hydroxylamine-detoxification mutualism with AOA through utilizing the AOA-released extracellular hydroxylamine using hydroxylamine oxidoreductase. Together, this study unravels the overlooked role of priorly unknown but abundant AOA lineages of the newly discovered genus Candidatus Nitrosomaritimum in biological nitrogen transformation and their potential for nitrogen pollution mitigation in coastal environments.

Comparative Genomics Screens Identify a Novel Small Secretory Peptide, SlSolP12, which Activates Both Local and Systemic Immune Response in Tomatoes and Exhibits Broad-Spectrum Activity.

Small secreted peptides (SSPs) are essential for defense mechanisms in plant-microbe interactions, acting as danger-associated molecular patterns (DAMPs). Despite the first discovery of SSPs over three decades ago, only a limited number of SSP families, particularly within Solanaceae plants, have been identified due to inefficient approaches. This study employed comparative genomics screens with Solanaceae proteomes (tomato, tobacco, and pepper) to discover a novel SSP family, SolP. Bioinformatics analysis suggests that SolP may serve as an endogenous signal initiating the plant PTI response. Interestingly, SolP family members from tomato, tobacco, and pepper share an identical sequence (VTSNALALVNRFAD), named SlSolP12 (also referred to as NtSolP15 or CaSolP1). Biochemical and phenotypic analyses revealed that synthetic SlSolP12 peptide triggers multiple defense responses: ROS burst, MAPK activation, callose deposition, stomatal closure, and expression of immune defense genes. Furthermore, SlSolP12 enhances systemic resistance against Botrytis cinerea infection in tomato plants and interferes with classical peptides, flg22 and Systemin, which modulate the immune response. Remarkably, SolP12 activates ROS in diverse plant species, such as Arabidopsis thaliana, soybean, and rice, showing a broad spectrum of biological activities. This study provides valuable approaches for identifying endogenous SSPs and highlights SlSolP12 as a novel DAMP that could serve as a useful target for crop protection.

Genomic Analysis of Cronobacter condimenti s37: Identification of Resistance and Virulence Genes and Comparison with Other Cronobacter and Closely Related Species.

Cronobacter condimenti are environmental commensals that have not been associated with any clinical infections. To date, they are the least understood and described Cronobacter species within the genus. The objective of this study was to use a draft genome sequence (DGS) of the Cronobacter condimenti strain s37 to screen for genes encoding for antibiotic resistance, virulence, response to environmental stress, and biofilm formation. The strain was isolated in Poland from commercial small radish sprouts. This is the second genome of this species available in the GenBank database. The comparative genome analysis (cgMLST) of C. condimenti s37 with other Cronobacter spp. including the pathogenic species C. sakazakii and the plant-associated closely related genera Franconibacter and Siccibacter was also performed. The assembled and annotated genome of the C. condimenti s37 genome was 4,590,991 bp in length, with a total gene number of 4384, and a GC content of 55.7%. The s 37 genome encoded for genes associated with resistance to stressful environmental conditions (metal resistance genes: zinc, copper, osmotic regulation, and desiccation stress), 17 antimicrobial resistance genes encoding resistance to various classes of antibiotics and 50 genes encoding for the virulence factors. The latter were mainly genes associated with adhesion, chemotaxis, hemolysis, and biofilm formation. Cg-MLST analysis (3991 genes) revealed a greater similarity of C. condimenti s37 to S. turicensis, F. pulveris, and C. dublinensis than to other species of the genus Cronobacter. Studies on the diversity, pathogenicity, and virulence of Cronobacter species isolated from different sources are still insufficient and should certainly be continued. Especially the analysis of rare strains such as s37 is very important because it provides new information on the evolution of these bacteria. Comparative cgMLST analysis of s37 with other Cronobacter species, as well as closely related genera Franconibacter and Siccibacter, complements the knowledge on their adaptability to specific environments such as desiccation.

Comparative genomic sequencing to characterize Mycoplasma pneumoniae genome, typing, and drug resistance.

To analyze the characteristics of Mycoplasma pneumoniae as well as macrolide antibiotic resistance through whole-genome sequencing and comparative genomics. Thirteen clinical strains isolated from 2003 to 2019 were selected, 10 of which were resistant to erythromycin (MIC >64 µg/mL), including 8 P1-type I and 2 P1-type II. Three were sensitive (<1 µg/mL) and P1-type II. One resistant strain had an A→G point mutation at position 2064 in region V of the 23S rRNA, the others had it at position 2063, while the three sensitive strains had no mutation here. Genome assembly and comparative genome analysis revealed a high level of genome consistency within the P1 type, and the primary differences in genome sequences concentrated in the region encoding the P1 protein. In P1-type II strains, three specific gene mutations were identified: C162A and A430G in L4 gene and T1112G mutation in the CARDS gene. Clinical information showed seven cases were diagnosed with severe pneumonia, all of which were infected with drug-resistant strains. Notably, BS610A4 and CYM219A1 exhibited a gene multi-copy phenomenon and shared a conserved functional domain with the DUF31 protein family. Clinically, the patients had severe refractory pneumonia, with pleural effusion, necessitating treatment with glucocorticoids and bronchoalveolar lavage. The primary variations between strains occur among different P1-types, while there is a high level of genomic consistency within P1-types. Three mutation loci associated with specific types were identified, and no specific genetic alterations directly related to clinical presentation were observed.IMPORTANCEMycoplasma pneumoniae is an important pathogen of community-acquired pneumonia, and macrolide resistance brings difficulties to clinical treatment. We analyzed the characteristics of M. pneumoniae as well as macrolide antibiotic resistance through whole-genome sequencing and comparative genomics. The work addressed primary variations between strains that occur among different P1-types, while there is a high level of genomic consistency within P1-types. In P1-type II strains, three specific gene mutations were identified: C162A and A430G in L4 gene and T1112G mutation in the CARDS gene. All the strains isolated from severe pneumonia cases were drug-resistant, two of which exhibited a gene multi-copy phenomenon, sharing a conserved functional domain with the DUF31 protein family. Three mutation loci associated with specific types were identified, and no specific genetic alterations directly related to clinical presentation were observed.

Organellar Genomes of Sargassum hemiphyllum var. chinense Provide Insight into the Characteristics of Phaeophyceae.

Sargassum hemiphyllum var. chinense, a prevalent seaweed along the Chinese coast, has economic and ecological significance. However, systematic positions within Sargassum and among the three orders of Phaeophyceae, Fucales, Ectocarpales, and Laminariales are in debate. Here, we reported the organellar genomes of S. hemiphyllum var. chinense (34,686-bp mitogenome with 65 genes and 124,323 bp plastome with 173 genes) and the investigation of comparative genomics and systematics of 37 mitogenomes and 22 plastomes of Fucales (including S. hemiphyllum var. chinense), Ectocarpales, and Laminariales in Phaeophyceae. Whole genome collinearity analysis showed gene number, type, and arrangement were consistent in organellar genomes of Sargassum with 360 SNP loci identified as S. hemiphyllum var. chinense and two genes (rps7 and cox2) identified as intrageneric classifications of Sargassum. Comparative genomics of the three orders of Phaeophyceae exhibited the same content and different types (petL was only found in plastomes of the order Fucales and Ectocarpales) and arrangements (most plastomes were rearranged, but trnA and trnD in the mitogenome represented different orders) in genes. We quantified the frequency of RNA-editing (canonical C-to-U) in both organellar genomes; the proportion of edited sites corresponded to 0.02% of the plastome and 0.23% of the mitogenome (in reference to the total genome) of S. hemiphyllum var. chinense. The repetition ratio of Fucales was relatively low, with scattered and tandem repeats (nine tandem repeats of 14-24 bp) dominating, while most protein-coding genes underwent negative selection (Ka/Ks < 1). Collectively, these findings provide valuable insights to guide future species identification and evolutionary status of three important Phaeophyceae order species.

Comparative genomic analysis of an emerging Pseudomonadaceae member, Thiopseudomonas alkaliphila.

Thiopseudomonas alkaliphila, an organism recently classified within the Pseudomonadaceae family, has been detected in diverse sources such as human tissues, animal guts, industrial fermenters, and decomposition environments, suggesting a diverse ecological role. However, a large knowledge gap exists in how T. alkaliphila functions. In this comparative genomic analysis, adaptations indicative of habitat specificity among strains and genomic similarity to known opportunistic pathogens are revealed. Genomic investigation reveals a core metabolic utilization of multiple oxidative and non-oxidative catabolic pathways, suggesting adaptability to varied environments and carbon sources. The genomic repertoire of T. alkaliphila includes secondary metabolites, such as antimicrobials and siderophores, indicative of its involvement in microbial competition and resource acquisition. Additionally, the presence of transposases, prophages, plasmids, and Clustered Regularly Interspaced Short Palindromic Repeats-Cas systems in T. alkaliphila genomes suggests mechanisms for horizontal gene transfer and defense against viral predation. This comprehensive genomic analysis expands our understanding on the ecological functions, community interactions, and potential virulence of T. alkaliphila, while emphasizing its adaptability and diverse capabilities across environmental and host-associated ecosystems.IMPORTANCEAs the microbial world continues to be explored, new organisms will emerge with beneficial and/or pathogenetic impact. Thiopseudomonas alkaliphila is a species originally isolated from clinical human tissue and fluid samples but has not been attributed to disease. Since its classification, T. alkaliphila has been found in animal guts, animal waste, decomposing remains, and biogas fermentation reactors. This is the first study to provide an in-depth view of the metabolic potential of publicly available genomes belonging to this species through a comparative genomics and draft pangenome calculation approach. It was found that T. alkaliphila is metabolically versatile and likely adapts to diverse energy sources and environments, which may make it useful for bioremediation and in industrial settings. A range of virulence factors and antibiotic resistances were also detected, suggesting T. alkaliphila may operate as an undescribed opportunistic pathogen.

Assessing a respiratory toxic infectious bronchitis virus (IBV) strain: isolation, identification, pathogenicity, and immunological failure insights.

Infectious bronchitis virus (IBV) is caused by avian coronavirus and poses a global economic threat to the poultry industry. In 2023, a highly pathogenic IBV strain, IBV/CN/GD20230501, was isolated and identified from chickens vaccinated with IBV-M41 in Guangdong, China. This study comprehensively investigated the biological characteristics of the isolated IBV strain, including its genotype, whole genome sequence analysis of its S1 gene, pathogenicity, host immune response, and serum non-targeted metabolomics. Through the analysis of the S1 gene sequence, serum neutralization tests, and comparative genomics, it was proven that IBV/CN/GD20230501 belongs to the GI-I type of strain and is serotype II. One alanine residue in the S1 subunit of the isolated strain was mutated into serine, and some mutations were observed in the ORF1ab gene and the terminal region of the genome. Animal challenge experiments using the EID50 and TCID50 calculations showed that IBV/CN/GD20230501 possesses strong respiratory pathogenicity, with early and long-term shedding of viruses and rapid viral spread. Antibody detection indicated that chickens infected with IBV/CN/GD20230501 exhibited delayed expression of early innate immune genes, while those infected with M41 showed rapid gene induction and effective viral control. Metabolomics analysis demonstrated that this virus infection led to differential expression of 291 ions in chicken serum, mainly affecting the citric acid cycle (tricarboxylic acid cycle).IMPORTANCEThis study identified an infectious bronchitis virus (IBV) strain isolated from vaccinated chickens in an immunized population that had certain sequence differences compared to IBV-M41, resulting in significantly enhanced pathogenicity and host defense. This strain has the potential to replace M41 as a more suitable challenge model for drug research. The non-targeted metabolomics analysis highlighting the citric acid cycle provides a new avenue for studying this highly virulent strain.

Leveraging coevolutionary insights and AI-based structural modeling to unravel receptor–peptide ligand-binding mechanisms

Secreted signaling peptides are central regulators of growth, development, and stress responses, but specific steps in the evolution of these peptides and their receptors are not well understood. Also, the molecular mechanisms of peptide-receptor binding are only known for a few examples, primarily owing to the limited availability of protein structural determination capabilities to few laboratories worldwide. Plants have evolved a multitude of secreted signaling peptides and corresponding transmembrane receptors. Stress-responsive SERINE RICH ENDOGENOUS PEPTIDES (SCOOPs) were recently identified. Bioactive SCOOPs are proteolytically processed by subtilases and are perceived by the leucine-rich repeat receptor kinase MALE DISCOVERER 1-INTERACTING RECEPTOR-LIKE KINASE 2 (MIK2) in the model plant Arabidopsis thaliana. How SCOOPs and MIK2 have (co)evolved, and how SCOOPs bind to MIK2 are unknown. Using in silico analysis of 350 plant genomes and subsequent functional testing, we revealed the conservation of MIK2 as SCOOP receptor within the plant order Brassicales. We then leveraged AI-based structural modeling and comparative genomics to identify two conserved putative SCOOP-MIK2 binding pockets across Brassicales MIK2 homologues predicted to interact with the "SxS" motif of otherwise sequence-divergent SCOOPs. Mutagenesis of both predicted binding pockets compromised SCOOP binding to MIK2, SCOOP-induced complex formation between MIK2 and its coreceptor BRASSINOSTEROID INSENSITIVE 1-ASSOCIATED KINASE 1, and SCOOP-induced reactive oxygen species production, thus, confirming our in silico predictions. Collectively, in addition to revealing the elusive SCOOP-MIK2 binding mechanism, our analytic pipeline combining phylogenomics, AI-based structural predictions, and experimental biochemical and physiological validation provides a blueprint for the elucidation of peptide ligand-receptor perception mechanisms.

In Silico Evaluation of Lawsonia intracellularis Genes Orthologous to Genes Associated with Pathogenesis in Other Intracellular Bacteria.

Proliferative enteropathy is an enteric disease caused by the bacterium Lawsonia intracellularis, which affects several species of domestic and wild animals. The mechanisms underlying the mechanisms employed by L. intracellularis to cause host cell proliferation are poorly understood, mostly because this bacterium is extremely difficult to isolate and propagate in vitro. Comparative genomics methods for searching for genes orthologous to genes known to be associated with pathogenesis allow identification of genes potentially involved in pathogenesis by the pathogen of interest. The goal of this study was to carry out in silico research on L. intracellularis genes orthologous to genes required for intracellular invasion and survival present in other pathogenic bacteria, particularly Brucella abortus, B. melitensis, B. suis, Listeria monocytogenes, Mycobacterium tuberculosis, Mycobacterium avium subspecies paratuberculosis, Salmonella enterica, Yersinia pestis, Y. enterocolitica, and Y. pseudotuberculosis. A total of 127 genes associated with invasion and intracellular survival from five known intracellular bacteria were mapped against the predicted proteomes of all L. intracellularis strains publicly available on GenBank, using the OrthoFinder program. A total of 45 L. intracellularis genes were orthologous to genes associated with pathogenesis of other intracellular bacteria. Genes putatively associated with signal the transduction of chemotaxis and cell motility were identified. Genes related to DNA binding and repair were also identified, with some of them supporting a possible association of bacteria with macrophages or inducing pro-inflammatory responses. The homology-based identification of these genes suggests their potential involvement in the virulence and pathogenicity of L. intracellularis, opening avenues for future research and insights into the molecular mechanisms of Lawsonia-elicited proliferative enteropathy.

Taxonomic characterization and comparative genomic analysis of a novel Devosia species revealed that phenolic acid-degrading traits are ubiquitous in the Devosia genus

Phenolic acids (PAs) are widely distributed allelochemicals in various environments. To better understand the fate of PAs in environments, a halotolerant PAs-degrading bacterium (named strain RR2S18T) isolated from rhizosphere soil was identified as a novel species of Devosia, named Devosia rhizosphaerae sp. nov. The strain initially degraded PAs into central ring-fission intermediates (protocatechuic acid) using the CoA-dependent non-β-oxidation pathway. The produced ring-fission intermediates were then consecutively degraded by an ortho-cleavage reaction and the β-ketoadipic acid pathway. A comparative genomics analysis of 62 Devosia strains revealed that PAs-degrading genes were ubiquitous in their genomes, indicating that PAs degradation is universal among members of this genus. The analysis also suggested that the genes involved in CoA-dependent non-β-oxidation are inherent to Devosia strains, while those involved in ring-fission and β-ketoadipic acid pathways were obtained by horizontal gene transfer.

Chalcopyrite bioleaching efficacy by extremely thermoacidophilic archaea leverages balanced iron and sulfur biooxidation

Factors that contribute to optimal chalcopyrite bioleaching by extremely thermoacidophilic archaea were examined for ten species belonging to the order Sulfolobales from the genera Acidianus (A. brierleyi), Metallosphaera (M. hakonensis, M. sedula, M. prunae), Sulfuracidifex (S. metallicus, S. tepriarius), Sulfolobus (S. acidocaldarius), Saccharlobus (S. solfataricus) and Sulfurisphaera (S. ohwakuensis, S. tokodaii). Only A. brierleyi, M. sedula, S. metallicus, S. tepriarius, S. ohwakuensis, and S. tokodai exhibited significant amounts of bioleaching and were investigated further. At 70–75 °C, Chalcopyrite loadings of 10 g/l were leached for 21 days during which pH, redox potential, planktonic cell density, iron concentrations and sulfate levels were monitored, in addition to copper mobilization. S. ohwakuensis proved to be the most prolific bioleacher. This was attributed to balanced iron and sulfur oxidation, thereby reducing by-product (e.g., jarosites) formation and minimizing surface passivation. Comparative genomics suggest markers for bioleaching potential, but the results here point to the need for experimental verification.

De Novo Emerged Gene Search in Eukaryotes with DENSE.

The discovery of de novo emerged genes, originating from previously noncoding DNA regions, challenges traditional views of species evolution. Indeed, the hypothesis of neutrally evolving sequences giving rise to functional proteins is highly unlikely. This conundrum has sparked numerous studies to quantify and characterize these genes, aiming to understand their functional roles and contributions to genome evolution. Yet, no fully automated pipeline for their identification is available. Therefore, we introduce DENSE (DE Novo emerged gene SEarch), an automated Nextflow pipeline based on two distinct steps: detection of taxonomically restricted genes (TRGs) through phylostratigraphy, and filtering of TRGs for de novo emerged genes via genome comparisons and synteny search. DENSE is available as a user-friendly command-line tool, while the second step is accessible through a web server upon providing a list of TRGs. Highly flexible, DENSE provides various strategy and parameter combinations, enabling users to adapt to specific configurations or define their own strategy through a rational framework, facilitating protocol communication, and study interoperability. We apply DENSE to seven model organisms, exploring the impact of its strategies and parameters on de novo gene predictions. This thorough analysis across species with different evolutionary rates reveals useful metrics for users to define input datasets, identify favorable/unfavorable conditions for de novo gene detection, and control potential biases in genome annotations. Additionally, predictions made for the seven model organisms are compiled into a requestable database, which we hope will serve as a reference for de novo emerged gene lists generated with specific criteria combinations.

The complete chloroplast genome of Rhododendronambiguum and comparative genomics of related species.

Rhododendron Linnaeus, 1753, the largest genus of woody plants in the Northern Hemisphere, includes some of the most significant species in horticulture. Rhododendronambiguum Hemsl, 1911, a member of subsection Triflora Sleumer 1947, exemplifies typical alpine Rhododendron species. The analysis of the complete chloroplast genome of R.ambiguum offers new insights into the evolution of Rhododendron species and enhances the resolution of phylogenetic relationships. This genome is composed of 207,478 base pairs, including a pair of inverted repeats (IRs) of 47,249 bp each, separated by a large single-copy (LSC) region of 110,367 bp and a small single-copy (SSC) region of 2,613 bp. It contains 110 genes: 77 protein-coding genes, 29 tRNAs, four unique rRNAs (4.5S, 5S, 16S, and 23S), with 16 genes duplicated in the IRs. Comparative analyses reveal substantial diversity in the Rhododendron chloroplast genome structures, identifying a fourth variant pattern. Specifically, four highly divergent regions (trnI-rpoB, ndhE-psaC, rpl32-ndhF, rrn16S-trnI) were noted in the intergenic spacers. Additionally, 76 simple sequence repeats were identified. Positive selection signals were detected in four genes (cemA, rps4, rpl16, and rpl14), evidenced by high Ka/Ks ratios. Phylogenetic reconstruction based on two datasets (shared protein-coding genes and complete chloroplast genomes) suggests that R.ambiguum is closely related to R.concinnum Hemsley, 1889. However, the phylogenetic positions of subsection Triflora Pojarkova, 1952 species remain unresolved, indicating that the use of complete chloroplast genomes for phylogenetic research in Rhododendron requires careful consideration. Overall, our findings provide valuable genetic information that will enhance understanding of the evolution, molecular biology, and genetic improvement of Rhododendron spieces.

Genomic Underpinnings of Cytoplasmic Incompatibility: CIF Gene-Neighborhood Diversification Through Extensive Lateral Transfers and Recombination in Wolbachia.

Cytoplasmic incompatibility (CI), a non-Mendelian genetic phenomenon, involves the manipulation of host reproduction by Wolbachia, a maternally transmitted alphaproteobacterium. The underlying mechanism is centered around the CI Factor (CIF) system governed by two genes, cifA and cifB, where cifB induces embryonic lethality, and cifA counteracts it. Recent investigations have unveiled intriguing facets of this system, including diverse cifB variants, prophage association in specific strains, copy number variation, and rapid component divergence, hinting at a complex evolutionary history. We utilized comparative genomics to systematically classify CIF systems, analyze their locus structure and domain architectures, and reconstruct their diversification and evolutionary trajectories. Our new classification identifies ten distinct CIF types, featuring not just versions present in Wolbachia, but also other intracellular bacteria, and eukaryotic hosts. Significantly, our analysis of CIF loci reveals remarkable variability in gene composition and organization, encompassing an array of diverse endonucleases, variable toxin domains, deubiquitinating peptidases (DUBs), prophages, and transposons. We present compelling evidence that the components within the loci have been diversifying their sequences and domain architectures through extensive, independent lateral transfers and interlocus recombination involving gene conversion. The association with diverse transposons and prophages, coupled with selective pressures from host immunity, likely underpins the emergence of CIF loci as recombination hotspots. Our investigation also posits the origin of CifB-REase domains from mobile elements akin to CR (Crinkler-RHS-type) effectors and Tribolium Medea1 factor, which is linked to another non-Mendelian genetic phenomenon. This comprehensive genomic analysis offers novel insights into the molecular evolution and genomic foundations of Wolbachia-mediated host reproductive control.

Deuterostome Ancestors and Chordate Origins.

The Deuterostomia are a monophyletic group, consisting of the Ambulacraria, with two phyla, Hemichordata and Echinodermata, and the phylum Chordata, containing the subphyla Cephalochordata (lancelets or Amphioxus), Tunicata (Urochordata) and Vertebrata. Hemichordates and echinoderms are sister groups and are critical for understanding the deuterostome ancestor and the origin and evolution of the chordates within the deuterostomes. Enteropneusta, worm-like hemichordates, share many chordate features as adults, including a post-anal tail, gill slits, and a Central Nervous System (CNS) that deploy similar developmental Genetic Regulatory Networks (GRNs). Genomic comparisons show that cephalochordates share synteny and a vermiform body plan similar to vertebrates, but phylogenomic analyses place tunicates as the sister group of vertebrates. Tunicates have a U-shaped gut and a very different adult body plan than the rest of the chordates, and all tunicates have small genomes and many gene losses, although the GRNs underlying specific tissues, such as notochord and muscle, are conserved. Echinoderms and vertebrates have extensive fossil records, with fewer specimens found for tunicates and enteropneusts, or worm-like hemichordates. The data is mounting that the deuterostome ancestor was a complex benthic worm, with gill slits, a cartilaginous skeleton, and a CNS. Two extant groups, echinoderms and tunicates, have evolved highly derived body plans, remarkably different than the deuterostome ancestor. We review the current genomic and GRN data on the different groups of deuterostomes' characters to re-evaluate different hypotheses of chordate origins. Notochord loss in echinoderms and hemichordates is as parsimonious as notochord gain in the chordates but has implications for the deuterostome ancestor. The chordate ancestor lost an ancestral nerve net, retained the central nervous system, and evolved neural crest cells.

Elucidating Scarab Divergence in an Evolutionary-Ecological Context through the Comprehensive Analysis of the Complete Mitogenome of Anomala.

Anomala Samouelle, 1819 is one of the specious genera of Coleoptera, with over 1000 known species, and includes some of the most destructive pests of crops or forests. Morphological convergence is a common phenomenon within this genus, making the identification of closely related species very difficult. To explore the phylogenetic placement of Anomalini and provide a basis for the classification and identification of Anomala, we comparatively analyzed the complete mitogenome of three Anomala species (A. exoleta, A. perplexa diana, and A. praecoxalis). Based on all accessible mitogenome data, we performed comparative mitochondrial genomics analysis of this genus and reconstructed the phylogenetic relationships of Scarabaeidae based on two datasets (protein-coding genes and amino acids) and two methods (Bayesian approach and maximum likelihood). The phylogenetic relationships found in this study highly support that the groups of Aphodiinae, Cetoniinae, Dynastinae, Rutelinae and Scarabaeinae are monophyletic. Interestingly, the phylogenetic clustering relationship was highly consistent with the Scarabaeidae diet, indicating that the herbivorous species and dung-eating species are clustered separately. The phylogenetic tree showed that the subfamily Melolonthinae and the genus Anomala are not monophyletic, suggesting that these two groups should be further revised with more data.

Characteristics and comparative genome analysis of Yersinia enterocolitica and related species associated with human infections in Switzerland 2019–2023

PurposeWe aimed to characterise Yersinia enterocolitica from human clinical specimens in Switzerland using epidemiological, microbiological and whole-genome sequencing (WGS) data. MethodsIsolates (n = 149) were collected between January 2019 and December 2023. Epidemiological data was noted and strains were characterized by biochemical and serological typing, antimicrobial susceptibility testing (AST), and WGS-based analysis. ResultsMost of the isolates (86%) were from stool specimens and 52% were from male patients. The patients' median age was 28 years (range < 1–94 years). Typing assigned the isolates to bioserotype 4/O:3 (44%), biotype 1A (34%), bioserotype 2/O:9 (21%), and bioserotype 3/O:3 (1%). WGS identified Y. enterocolitica (n = 147), Y. alsatica (n = 1) and Y. proxima (n = 1). Seven isolates were multidrug resistant (MDR) and harboured plasmid pAB829 carrying aph(3″)-Ib, aph(6)-Id, and tet(Y) (n = 1), pAC120 carrying aph(6)-Id and tet(A) (n = 2), or a 12.6 kb Tn2670-like transposon containing catA1, aadA12, sul1, and qacEΔ1 (n = 4). Virulence factors (VFs) included ail (n = 99), invB, (n = 145), ystA (n = 99), ystB (n = 48) and pYV-associated VFs (n = 93). MLST and cgMLST analysis showed that BT 1A strains consisted of several STs and were highly diverse, whereas BT 2/O:9 strains were all ST12 and clustered closely, and BT 4/O:3 strains mostly belonged to ST18 but were more diverse. SNP analysis revealed two highly clonal BT 4/O:3 subpopulations with wide spatio-temporal distribution. ConclusionsY. enterocolitica BT 1A, BT 2/O:9 and BT 4/O:3 are frequently associated with human yersiniosis in Switzerland. WGS-based subtyping of Y. enterocolitica is a powerful tool to explore the genetic diversity and the pathogenic potential of human isolates.

Turbinal bones are still one of the last frontiers of the tetrapod skull: hypotheses, challenges and perspectives.

Turbinals are bony or cartilaginous structures that are present in the nasal cavity of most tetrapods. They are involved in key functions such as olfaction, heat, and moisture conservation, as well as protection of the respiratory tract. Despite recent studies that challenged long-standing hypotheses about their physiological and genomic correlation, turbinals remain largely unexplored, particularly for non-mammalian species. Herein, we review and synthesise the current knowledge of turbinals using an integrative approach that includes comparative anatomy, physiology, histology and genomics. In addition, we provide synonyms and correspondences of tetrapod turbinals from about 80 publications. This work represents a first step towards drawing hypotheses of homology for the whole clade, and provides a strong basis to develop new research avenues.

Unravelling the novel sex determination genotype with 'ZY' and a distinctive 2.15-2.95 Mb inversion among poplar species through haplotype-resolved genome assembly and comparative genomics analysis.

Populus tomentosa, an indigenous tree species, is widely distributed and cultivated over 1,000,000 km2 in China, contributing significantly to forest production, ecological conservation and urban-rural greening. Although a reference genome is available for P. tomentosa, the intricate interspecific hybrid origins, chromosome structural variations (SVs) and sex determination mechanisms remain confusion and unclear due to its broad and even overlapping geographical distribution, extensive morphological variations and cross infiltration among white poplar species. We conducted a haplotype-resolved de novo assembly of P. tomentosa elite individual GM107, which comprises subgenomes a and b with a total genome size of 714.9 Mb. We then analysed the formation of hybrid species and the phylogenetic evolution and sex differentiation across the entire genus. Phylogenomic analyses suggested that GM107 likely originated from a hybridisation event between P. alba (♀) and P. davidiana (♂) which diverged at approximately 3.8 Mya. A total of 1551 chromosome SVs were identified between the two subgenomes. More noteworthily, a distinctive inversion structure spanning 2.15-2.95 Mb was unveiled among Populus, Tacamahaca, Turaga, Aigeiros poplar species and Salix, highlighting a unique evolutionary feature. Intriguingly, a novel sex genotype of the ZY type, which represents a crossover between XY and ZW systems, was identified and confirmed through both natural and artificial hybrids populations. These novel insights offer significant theoretical value for the study of the species' evolutionary origins and serve as a valuable resource for ecological genetics and forest biotechnology.

Comparative Genomics of Three Hybrid-Pathogen Multidrug-Resistant Escherichia coli Strains Isolated from Healthy Donors’ Feces

Comparative Genomics of Three Hybrid-Pathogen Multidrug-Resistant Escherichia coli Strains Isolated from Healthy Donors’ Feces