Parallelization of MAFFT for large-scale multiple sequence alignments

The tetraploid Camellia oleifera genome provides insights into evolution, agronomic traits, and genetic architecture of oil Camellia plants

Camellia oleifera is an economically important woody oil plant. Complex ploidy and lack of genomic information have seriously hindered the molecular breeding of C. oleifera. Here, we present an 11.43-Gb haplotype-resolved, chromosome-level genome assembly of tetraploid C. oleifera (COL-tetra). Methods employed in this study support the conclusion that COL-tetra is an autotetraploid and probably originates from genome doubling of the diploid C. brevistyla. In addition, DNA methylation plays a significant role in imbalanced allelic expression and seed development. Genetic divergence analyses reveal significant differentiation signals for flowering time between spring-flowering and autumn-flowering oil Camellia species, which probably account for reproductive isolation between species with distinct flowering times. Strong introgression signals are detected between COL-tetra and C. sasanqua and between C. vietnamensis and COL-hexa, which might affect the development of agronomic traits and environmental adaptability. This study provides valuable insights into the evolution, agronomic trait development, and genetic architecture of oil Camellia plants.

Membrane lipid homeostasis dually regulates conformational transition of phosphoethanolamine transferase EptA

The phosphoethanolamine transferase EptA utilizes phosphatidylethanolamine (PE) in the bacterial cell membrane to modify the structure of lipopolysaccharide, thereby conferring antimicrobial resistance on Gram-negative pathogens. Previous studies have indicated that excessive consumption of PE can disrupt the cell membrane, leading to cell death. This implies the presence of a regulatory mechanism for EptA catalysis to maintain a balance between antimicrobial resistance and bacterial growth. Through microsecond-scale all-atom molecular dynamics simulations, we demonstrate that membrane lipid homeostasis modulates the conformational transition and catalytic activation of EptA. The conformation of EptA oscillates between closed and open states, ensuring the precise spatiotemporal sequence of substrates binding. Interestingly, the conformation of EptA is significantly influenced by its surrounding lipid microenvironment, particularly the PE proportion in the membrane. PE-rich membrane conditions initiate and stabilize the open conformation of EptA through both orthosteric and allosteric effects. Importantly, the reaction mediated by EptA gradually depletes PE in the membrane, ultimately hindering its conformational transition and catalytic activation. These findings collectively establish a self-promoted model, illustrating the regulatory mechanism of EptA during the development of antibiotic resistance.

Phased chromosome-level genome provides insights into the molecular adaptation for migratory lifestyle and population diversity for Pacific saury, Cololabis saira

The Pacific saury (Cololabis saira) is a pelagic fish commonly found in the North Pacific Ocean. Its population diversity and migratory lifestyle have long captured global attention. Despite the inherent complexity of the C. saira genome, characterized by extremely high heterozygosity, we successfully assembled a phased chromosome-level genome. The genome analysis revealed the expansion and natural selection of numerous functional genes, likely contributing to its enduring and extensive migratory lifestyle. Notably, gpr35 and igh genes showed significant expansion in the C. saira genome, potentially associated with regulating the immune response against environmental parasites and pathogens. Moreover, genes involved in DNA repair/replication and peroxisome function, including atm, ercc6, pex14, and pex16, displayed evidence of positive selection. Based on genome-sequencing of 80 individuals from eight sampling sites, we demonstrated that the genomic divergence among C. saira populations is relatively low. However, the sampling sites could be grouped into two distinct clusters, roughly corresponding to the migratory route of C. saira. This suggests a possible genome-wide divergence for C. saira within the open ocean region. Furthermore, the trmu gene, responsible for controlling otolith development and sharpness, exhibited differentiation between the two groups, consistent with previously reported differences in otolith morphology. This study has provided a reference genome and insights into the evolution, ecology, and conservation of Pacific saury and closely-related species.

Viroid-like colonists of human microbiomes

Here, we describe "obelisks," a class of heritable RNA elements sharing several properties: (1) apparently circular RNA ∼1 kb genome assemblies, (2) predicted rod-like genome-wide secondary structures, and (3) open reading frames encoding a novel "Oblin" protein superfamily. A subset of obelisks includes a variant hammerhead self-cleaving ribozyme. Obelisks form their own phylogenetic group without detectable similarity to known biological agents. Surveying globally, we identified 29,959 distinct obelisks (clustered at 90% sequence identity) from diverse ecological niches. Obelisks are prevalent in human microbiomes, with detection in ∼7% (29/440) and ∼50% (17/32) of queried stool and oral metatranscriptomes, respectively. We establish Streptococcus sanguinis as a cellular host of a specific obelisk and find that this obelisk's maintenance is not essential for bacterial growth. Our observations identify obelisks as a class of diverse RNAs of yet-to-be-determined impact that have colonized and gone unnoticed in human and global microbiomes.

What are microsatellites and how to choose the best tool: a user-friendly review of SSR and 74 SSR mining tools

Microsatellites, also known as SSR or STR, are essential molecular markers in genomic research, playing crucial roles in genetic mapping, population genetics, and evolutionary studies. Their applications range from plant breeding to forensics, highlighting their diverse utility across disciplines. Despite their widespread use, traditional methods for SSR analysis are often laborious and time-consuming, requiring significant resources and expertise. To address these challenges, a variety of computational tools for SSR analysis have been developed, offering faster and more efficient alternatives to traditional methods. However, selecting the most appropriate tool can be daunting due to rapid technological advancements and the sheer number of options available. This study presents a comprehensive review and analysis of 74 SSR tools, aiming to provide researchers with a valuable resource for SSR analysis tool selection. The methodology employed includes thorough literature reviews, detailed tool comparisons, and in-depth analyses of tool functionality. By compiling and analyzing these tools, this study not only advances the field of genomic research but also contributes to the broader scientific community by facilitating informed decision-making in the selection of SSR analysis tools. Researchers seeking to understand SSRs and select the most appropriate tools for their projects will benefit from this comprehensive guide. Overall, this study enhances our understanding of SSR analysis tools, paving the way for more efficient and effective SSR research in various fields of study.

Description of three new species of Gyrodactylus von Nordmann, 1832 (Monopisthocotylea: Gyrodactylidae) on Triplophysa species (Nemacheilidae) from Qinghai-Tibet Plateau

Three new species of Gyrodactylus were identified from the body surface of the Triplophysa species from the Qinghai-Tibet Plateau, Gyrodactylus triplorienchili n. sp. on Triplophysa orientalis in northern Tibet, G. yellochili n. sp. on T. sellaefer and T. scleroptera and G. triplsellachili n. sp. on T. sellaefer and T. robusta in Lanzhou Reach of the Yellow River. The three newly identified species share the nemachili group species' characteristic of having inturning hamulus roots. Gyrodactylus triplorienchili n. sp. shared a quadrate sickle heel and a thin marginal hook sickle, two morphological traits that set them apart from G. yellochili n. sp. However, they may be identified by the distinct shapes of the sickle base and marginal hook sickle point. Gyrodactylus triplsellachili n. sp. had much larger opisthaptoral hard part size than the other two species. The three new species show relatively low interspecific differences of 2.9-5.3% p-distance for ITS1-5.85-ITS2 rDNA sequences. Phylogenetic analysis indicated that the three new species formed a well-supported monophyletic group (bp = 99) with the other nemachili group species.

Genome sequencing and transcriptome analysis provide an insight into the hypoxia resistance of Channa asiatica

Channa asiatica is an economically valuable fish species and excellent model for studying hypoxic tolerance. However, the underlying genetic and molecular mechanisms are poorly understood. In this study, we assembled a high-quality C. asiatica genome (23 chromosomes, totaling 722 Mb) using a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing. Repetitive elements accounted for 28.39%of the C. asiatica genome, and 23,949 protein-coding genes were predicted, with 96.63 % of these functionally annotated. Moreover, a comparative genomic analysis of 12 fish genomes showed that gene families associated with oxygen binding and transport were expanded in C. asiatica. In addition, transcriptome analysis revealed that multiple oxidative stress pathways were activated when C. asiatica was exposed to air. In conclusion, this study provided high-quality genome assembly and transcriptome data, both serving as critical resources for researching the genetic basis of hypoxic tolerance in C. asiatica.

Leaf pigmentation in Cannabis sativa: Characterization of anthocyanin biosynthesis in colorful Cannabis varieties

Cannabis plants produce a spectrum of secondary metabolites, encompassing cannabinoids and more than 300 non-cannabinoid compounds. Among these, anthocyanins have important functions in plants and also have well documented health benefits. Anthocyanins are largely responsible for the red/purple color phenotypes in plants. Although some well-known Cannabis varieties display a wide range of red/purple pigmentation, the genetic underpinnings of anthocyanin biosynthesis have not been well characterized in Cannabis. This study unveils the genetic diversity of anthocyanin biosynthesis genes found in Cannabis, and we characterize the diversity of anthocyanins and related phenolics found in four differently pigmented Cannabis varieties. Our investigation revealed that the genes 4CL, CHS, F3H, F3'H, FLS, DFR, ANS, and OMT exhibited the strongest correlation with anthocyanin accumulation in Cannabis leaves. The results of this study enhance our understanding of the anthocyanin biosynthetic pathway and shed light on the molecular mechanisms governing Cannabis leaf pigmentation.

Protein semisynthesis reveals plasticity in HECT E3 ubiquitin ligase mechanisms

Lys ubiquitination is catalysed by E3 ubiquitin ligases and is central to the regulation of protein stability and cell signalling in normal and disease states. There are gaps in our understanding of E3 mechanisms, and here we use protein semisynthesis, chemical rescue, microscale thermophoresis and other biochemical approaches to dissect the role of catalytic base/acid function and conformational interconversion in HECT-domain E3 catalysis. We demonstrate that there is plasticity in the use of the terminal side chain or backbone carboxylate for proton transfer in HECT E3 ubiquitin ligase reactions, with yeast Rsp5 orthologues appearing to be possible evolutionary intermediates. We also show that the HECT-domain ubiquitin covalent intermediate appears to eject the E2 conjugating enzyme, promoting catalytic turnover. These findings provide key mechanistic insights into how protein ubiquitination occurs and provide a framework for understanding E3 functions and regulation.

The key role of iroBCDN-lacking pLVPK-like plasmid in the evolution of the most prevalent hypervirulent carbapenem-resistant ST11-KL64 Klebsiella pneumoniae in China

AIMS

Hypervirulent carbapenem-resistant Klebsiella pneumoniae (hv-CRKP), coharboring hypervirulence and carbapenem-resistance genes mediated by plasmids, causes infections with extremely high mortality and seriously impacts public health. Exploring the transfer mechanisms of virulence/carbapenem-resistance plasmids, as well as the formation and evolution pathway of hv-CRKP is of great significance to the control of hv-CRKP infections.

METHODS

In this study, we identified the predominant clone of hv-CRKP in China and elucidated its genomic characteristics and formation route based on 239 multicenter clinical K. pneumoniae isolates and 1014 GenBank genomes by using comparative genomic analysis. Further, we revealed the factors affecting the transfer of virulence plasmids, and explained the genetic foundation for the prevalence of Chinese predominant hv-CRKP clone.

RESULTS

ST11-KL64 is the predominant clone of hv-CRKP in China and primarily evolved from ST11-KL64 CRKP by acquiring the pLVPK-like virulence plasmid from hvKP. Significantly, the virulence gene cluster iroBCDN was lost in the virulence plasmid of ST11-KL64 hv-CRKP but existed in that of hvKP. Moreover, the absence of iroBCDN didn't decrease the virulence of hv-CRKP, which was proved by bacterial test, cell-interaction test and mice infection model. On the contrary, loss of iroBCDN was observed to regulate virulence/carbapenem-resistance plasmid transfer and oxidative stress-related genes in strains and thus promoted the mobilization of nonconjugative virulence plasmid from hvKP into ST11-KL64 CRKP, forming hv-CRKP which finally had elevated antioxidant capacity and enhanced survival capacity in macrophages. The loss of iroBCDN increased the survival ability of hv-CRKP without decreasing its virulence, endowing it with an evolutionary advantage.

CONCLUSIONS

Our work provides new insights into the key role of iroBCDN loss in convergence of CRKP and hvKP, and the genetic and biological foundation for the widespread prevalence of ST11-KL64 hv-CRKP in China.

Characterization of a Bifidobacterium animalis subsp. lactis reference strain based on ecology and transcriptomics

Bifidobacteria are recognized as health-promoting bacteria that reside in the human gut, helping in the digestion of fiber, preventing infections, and producing essential compounds like vitamins. To date, Bifidobacterium animalis subsp. lactis, together with Bifidobacterium adolescentis, Bifidobacterium bifidum, Bifidobacterium breve, and Bifidobacterium longum, represents one of the species that are used as probiotic bacteria. Despite the extensive and detailed scientific research conducted on this microbial taxon, the molecular mechanisms by which B. animalis subsp. lactis exerts health benefits to its host are still largely unknown. Thus, we dissected the genetic repertoire and phylogenetic relationship of 162 strains of B. animalis subsp. lactis to select a representative reference strain of this taxon suitable for investigating its interaction with the host. The B. animalis subsp. lactis PRL2013 strain, which was isolated by a mucosal sample of a healthy adult, was chosen as the reference of the monophyletic cluster of human origin and revealed a greater adhesion index than that observed for another B. animalis subsp. lactis strain used in the industry as a probiotic supplement. Transcriptomics analyses of PRL2013 strain, when exposed to human cell monolayers, revealed 291 significantly upregulated genes, among which were found genes predicted to encode extracellular structures that may directly interact with human cells, such as extracellular polymeric substances, wall teichoic acids, and pili.

IMPORTANCE

To date, many Bifidobacterium animalis subsp. lactis strains have been isolated from human fecal samples. However, their presence in these samples does not necessarily suggest an ability to colonize the human gut. Furthermore, probiotics of non-human origin may not effectively interact with the gut epithelium, resulting in transient bacteria of the gut microbiota. In vitro experiments with human cells revealed that B. animalis subsp. lactis PRL2013, an autochthonous member of the human gut, shows colonization capability, leading to future applications in functional foods.

Sarbecovirus RBD indels and specific residues dictating multi-species ACE2 adaptiveness

Our comprehensive understanding of the multi-species ACE2 adaptiveness of sarbecoviruses remains elusive, particularly for those with various receptor binding motif (RBM) insertions/deletions (indels). Here, we analyzed RBM sequences from 268 sarbecoviruses categorized into four RBM indel types. We examined the ability of 20 representative sarbecovirus Spike glycoproteins (S) and derivatives in utilizing ACE2 from various bats and several other mammalian species. We reveal that sarbecoviruses with long RBMs (type-I) can achieve broad ACE2 tropism, whereas viruses with single deletions in Region 1 (type-II) or Region 2 (type-III) exhibit narrower ACE2 tropism. Sarbecoviruses with double region deletions (type-IV) completely lost ACE2 usage, which is restricted by clade-specific residues within and outside RBM. Lastly, we propose the evolution of sarbecovirus RBM indels and illustrate how loop lengths, disulfide, and residue determinants shape multi-species ACE2 adaptiveness. This study provides profound insights into the mechanisms governing ACE2 usage and spillover risks of sarbecoviruses.

From light into shadow: comparative plastomes in Petrocosmea and implications for low light adaptation

BACKGROUND

Plastids originated from an ancient endosymbiotic event and evolved into the photosynthetic organelles in plant cells. They absorb light energy and carbon dioxide, converting them into chemical energy and oxygen, which are crucial for plant development and adaptation. However, little is known about the plastid genome to light adaptation. Petrocosmea, a member of the Gesneriaceae family, comprises approximately 70 species with diverse light environment, serve as an ideal subject for studying plastomes adapt to light.

RESULTS

In this study, we selected ten representative species of Petrocosmea from diverse light environments, assembled their plastid genomes, and conducted a comparative genomic analysis. We found that the plastid genome of Petrocosmea is highly conserved in both structure and gene content. The phylogenetic relationships reconstructed based on the plastid genes were divided into five clades, which is consistent with the results of previous studies. The vast majority of plastid protein-coding genes were under purifying selection, with only the rps8 and rps16 genes identified under positive selection in different light environments. Notably, significant differences of evolutionary rate were observed in NADH dehydrogenase, ATPase ribosome, and RNA polymerase between Clade A and the other clades. Additionally, we identified ycf1 and several intergenic regions (trnH-psbA, trnK-rps16, rpoB-trnC, petA-psbJ, ccsA-trnL, rps16-trnQ, and trnS-trnG) as candidate barcodes for this emerging ornamental horticulture.

CONCLUSION

We newly assembled ten plastid genomes of Petrocosmea and identified several hypervariable regions, providing genetic resources and candidate markers for this promising emerging ornamental horticulture. Furthermore, our study suggested that rps8 and rps16 were under positive selection and that the evolutionary patterns of NADH dehydrogenase, ATPase ribosome, and RNA polymerase were related to the diversity light environment in Petrocosmea. This revealed an evolutionary scenario for light adaptation of the plastid genome in plants.

Molecular identification and studies on genetic diversity and structure-related GC heterogeneity of Spatholobus Suberectus based on ITS2

To determine the role of internal transcribed spacer 2 (ITS2) in the identification of Spatholobus suberectus and explore the genetic diversity of S. suberectus. A total of 292 ITS2s from S. suberectus and 17 other plant species were analysed. S. suberectus was clustered separately in the phylogenetic tree. The genetic distance between species was greater than that within S. suberectus. Synonymous substitution rate (Ks) analysis revealed that ITS2 diverged the most recently within S. suberectus (Ks = 0.0022). These findings suggested that ITS2 is suitable for the identification of S. suberectus. The ITS2s were divided into 8 haplotypes and 4 evolutionary branches on the basis of secondary structure, indicating that there was variation within S. suberectus. Evolutionary analysis revealed that the GC content of paired regions (pGC) was greater than that of unpaired regions (upGC), and the pGC showed a decreasing trend, whereas the upGC remained unchanged. Single-base mutation was the main cause of base pair substitution. In both the initial state and the equilibrium state, the substitution rate of GC was higher than that of AU. The increase in the GC content was partly attributed to GC-biased gene conversion (gBGC). High GC content reflected the high recombination and mutation rates of ITS2, which is the basis for species identification and genetic diversity. We characterized the sequence and structural characteristics of S. suberectus ITS2 in detail, providing a reference and basis for the identification of S. suberectus and its products, as well as the protection and utilization of wild resources.

The Evolution of Extreme Genetic Variability in a Parasite-Resistance Complex

Genomic regions that play a role in parasite defense are often found to be highly variable, with the major histocompatibility complex serving as an iconic example. Single nucleotide polymorphisms may represent only a small portion of this variability, with Indel polymorphisms and copy number variation further contributing. In extreme cases, haplotypes may no longer be recognized as orthologous. Understanding the evolution of such highly divergent regions is challenging because the most extreme variation is not visible using reference-assisted genomic approaches. Here we analyze the case of the Pasteuria Resistance Complex in the crustacean Daphnia magna, a defense complex in the host against the common and virulent bacterium Pasteuria ramosa. Two haplotypes of this region have been previously described, with parts of it being nonhomologous, and the region has been shown to be under balancing selection. Using pan-genome analysis and tree reconciliation methods to explore the evolution of the Pasteuria Resistance Complex and its characteristics within and between species of Daphnia and other Cladoceran species, our analysis revealed a remarkable diversity in this region even among host species, with many nonhomologous hyper-divergent haplotypes. The Pasteuria Resistance Complex is characterized by extensive duplication and losses of Fucosyltransferase (FuT) and Galactosyltransferase (GalT) genes that are believed to play a role in parasite defense. The Pasteuria Resistance Complex region can be traced back to common ancestors over 250 million years. The unique combination of an ancient resistance complex and a dynamic, hyper-divergent genomic environment presents a fascinating opportunity to investigate the role of such regions in the evolution and long-term maintenance of resistance polymorphisms. Our findings offer valuable insights into the evolutionary forces shaping disease resistance and adaptation, not only in the genus Daphnia, but potentially across the entire Cladocera class.

Identification and Characterization of Four Novel Viruses in Balclutha incisa

Balclutha incisa (Cicadellidae: Deltocephalinae), a leafhopper prevalent in tropical and temperate regions, is notably abundant in grasses and rice. The virome of B. incisa was investigated using deep transcriptome sequencing, leading to the first identification of four viruses belonging to the families Aliusviridae, Iflaviridae, and Totiviridae in B. incisa. These viruses have been provisionally named B. incisa ollusvirus 1 (BiOV1), B. incisa ollusvirus 2 (BiOV2), B. incisa iflavirus 1 (BiIV1), and B. incisa totivirus 1 (BiTV1). The complete genome sequences of these viruses were obtained through rapid amplification of cDNA ends (RACE). BiOV1 has a linear genome of 15,125 nucleotides (nt), while BiOV2 possesses a circular genome of 14,853 nt. The BiIV1 genome, excluding the poly(A) tail, is 10,903 nt in length and encodes a single open reading frame (ORF) for a polyprotein consisting of 3194 amino acids (aa). The BiTV1 genome is 4357 nt long and contains two overlapping ORFs, with the viral RNA-dependent RNA polymerase (RdRp) translated via a -1 ribosomal frameshift. Phylogenetic and sequence identity analyses suggest that all these viruses are novel members of their respective families. This study significantly expands our understanding of the virome associated with B. incisa by reporting and characterizing these novel viruses.

Comparative genomics unravels a rich set of biosynthetic gene clusters with distinct evolutionary trajectories across fungal species (Termitomyces) farmed by termites

The use of compounds produced by hosts or symbionts for defence against antagonists has been identified in many organisms, including in fungus-farming termites (Macrotermitinae). The obligate mutualistic fungus Termitomyces plays a pivotal role in plant biomass decomposition and as the primary food source for these termites. Despite the isolation of various specialized metabolites from different Termitomyces species, our grasp of their natural product repertoire remains incomplete. To address this knowledge gap, we conducted a comprehensive analysis of 39 Termitomyces genomes, representing 21 species associated with members of five termite host genera. We identified 754 biosynthetic gene clusters (BGCs) coding for specialized metabolites and categorized 660 BGCs into 61 biosynthetic gene cluster families (GCFs) spanning five compound classes. Seven GCFs were shared by all 21 Termitomyces species and 21 GCFs were present in all genomes of subsets of species. Evolutionary constraint analyses on the 25 most abundant GCFs revealed distinctive evolutionary histories, signifying that millions of years of termite-fungus symbiosis have influenced diverse biosynthetic pathways. This study unveils a wealth of non-random and largely undiscovered chemical potential within Termitomyces and contributes to our understanding of the intricate evolutionary trajectories of biosynthetic gene clusters in the context of long-standing symbiosis.

Chromosome-scale genome assembly of Astragalus membranaceus using PacBio and Hi-C technologies

Astragalus membranaceus (Fisch.) Bge (AM) is a medicinal herb plant belonging to the Leguminosae family. In this study, we present a chromosome-scale genome assembly of AM, aiming to enhance the molecular biology and functional studies of Astragali Radix. The genome size of AM is about 1.43 Gb, with a contig N50 value of 1.67 Mb. A total of 98.16% of the assembly anchored to 9 pseudochromosomes using Hi-C technology. The assembly completeness was estimated to be 97.27% using BUSCO with the long terminal repeat assembly index (LAI) of 16.22 and quality value (QV) of 48.58. Additionally, the genome contained 67.98% repetitive sequences. Genome annotation predicted 29,914 protein-coding genes, including 73 genes involved in the flavonoid biosynthetic pathway and 2,048 transcription factors. The high-quality genome assembly and gene annotation resources will greatly facilitate future functional genomic studies in Leguminosae species.

Oreocharis scutifolia (Gesneriaceae), a Peltate-Leaved New Species From the Dry-Hot Valley of the Jinsha River Basin, Yunnan, China

A peltate-leaved new species, Oreocharis scutifolia Z. Xie & H. H. Kong, endemic to the Dry-Hot Valley of the Jinsha River Basin, Yunnan, China, is described and illustrated here. It is similar to O. cordatula (Craib) Pellegr. and O. aurantiaca Baill. in floral characters but differs in its peltate leaf blades, which are unique in the genus (and only occur in one population of O. henryana Oliv., but its flowers are smaller, campanulate, and deep purple). Molecular phylogenetic analysis based on transcriptome data confirmed its systematic position to be sister group with O. henryana, O. cordatula, O. minor Pellegr., and O. aurantiaca Baill. (LPP = 1), and well apart from the remaining members of the genus. The new species O. scutifolia is assessed as "Critically Endangered" following the IUCN categories and criteria, due to its small and single population, thus making it face serious threats from human disturbance, invasive plants, shrinking habitat, and decreasing habitat quality.

The perception and evolution of flagellin, cold shock protein and elongation factor Tu from vector-borne bacterial plant pathogens

Vector-borne bacterial pathogens cause devastating plant diseases that cost billions of dollars in crop losses worldwide. These pathogens have evolved to be host- and vector-dependent, resulting in a reduced genome size compared to their free-living relatives. All known vector-borne bacterial plant pathogens belong to four different genera: 'Candidatus Liberibacter', 'Candidatus Phytoplasma', Spiroplasma and Xylella. To protect themselves against pathogens, plants have evolved pattern recognition receptors that can detect conserved pathogen features as non-self and mount an immune response. To gain an understanding of how vector-borne pathogen features are perceived in plants, we investigated three proteinaceous features derived from cold shock protein (csp22), flagellin (flg22) and elongation factor Tu (elf18) from vector-borne bacterial pathogens as well as their closest free-living relatives. In general, vector-borne pathogens have fewer copies of genes encoding flagellin and cold shock protein compared to their closest free-living relatives. Furthermore, epitopes from vector-borne pathogens were less likely to be immunogenic compared to their free-living counterparts. Most Liberibacter csp22 and elf18 epitopes do not trigger plant immune responses in tomato or Arabidopsis. Interestingly, csp22 from the citrus pathogen 'Candidatus Liberibacter asiaticus' triggers immune responses in solanaceous plants, while csp22 from the solanaceous pathogen 'Candidatus Liberibacter solanacearum' does not. Our findings suggest that vector-borne plant pathogenic bacteria evolved to evade host recognition.

Diverse RNA viruses of parasitic nematodes can elicit antibody responses in vertebrate hosts

Parasitic nematodes have an intimate, chronic and lifelong exposure to vertebrate tissues. Here we mined 41 published parasitic nematode transcriptomes from vertebrate hosts and identified 91 RNA viruses across 13 virus orders from 24 families in ~70% (28 out of 41) of parasitic nematode species, which include only 5 previously reported viruses. We observe widespread distribution of virus-nematode associations across multiple continents, suggesting an ancestral acquisition event and host-virus co-evolution. Characterization of viruses of Brugia malayi (BMRV1) and Onchocerca volvulus (OVRV1) shows that these viruses are abundant in reproductive tissues of adult parasites. Importantly, the presence of BMRV1 RNA in B. malayi parasites mounts an RNA interference response against BMRV1 suggesting active viral replication. Finally, BMRV1 and OVRV1 were found to elicit antibody responses in serum samples from infected jirds and infected or exposed humans, indicating direct exposure to the immune system.

PanKA: Leveraging population pangenome to predict antibiotic resistance

Machine learning has the potential to be a powerful tool in the fight against antimicrobial resistance (AMR), a critical global health issue. Machine learning can identify resistance mechanisms from DNA sequence data without prior knowledge. The first step in building a machine learning model is a feature extraction from sequencing data. Traditional methods like single nucleotide polymorphism (SNP) calling and k-mer counting yield numerous, often redundant features, complicating prediction and analysis. In this paper, we propose PanKA, a method using the pangenome to extract a concise set of relevant features for predicting AMR. PanKA not only enables fast model training and prediction but also improves accuracy. Applied to the Escherichia coli and Klebsiella pneumoniae bacterial species, our model is more accurate than conventional and state-of-the-art methods in predicting AMR.

Approaching Optimal pH Enzyme Prediction with Large Language Models

Enzymes are widely used in biotechnology due to their ability to catalyze chemical reactions: food making, laundry, pharmaceutics, textile, brewing─all these areas benefit from utilizing various enzymes. Proton concentration (pH) is one of the key factors that define the enzyme functioning and efficiency. Usually there is only a narrow range of pH values where the enzyme is active. This is a common problem in biotechnology to design an enzyme with optimal activity in a given pH range. A large part of this task can be completed in silico, by predicting the optimal pH of designed candidates. The success of such computational methods critically depends on the available data. In this study, we developed a language-model-based approach to predict the optimal pH range from the enzyme sequence. We used different splitting strategies based on sequence similarity, protein family annotation, and enzyme classification to validate the robustness of the proposed approach. The derived machine-learning models demonstrated high accuracy across proteins from different protein families and proteins with lower sequence similarities compared with the training set. The proposed method is fast enough for the high-throughput virtual exploration of protein space for the search for sequences with desired optimal pH levels.

Impact of weight variation on the microbiome of yak dams and calves

Introduction

Limited information exists regarding the microbiome composition of yak calves of varying weights. Therefore, this study aimed to investigate the microbiomes of mother-calf pairs with different weight profiles.

Methods

Fecal and blood samples were collected from both lower-weight (CB) and higher-weight (HB) yak calves, along with their corresponding female yaks (CA, HA).

Results

The results revealed significantly higher levels of T-AOC (total antioxidant capacity) and GSH-Px (glutathione peroxidase) in HB animals (p < 0.001). Sequencing yielded 652,181 and 643,369 filtered reads in female and calf yaks, respectively. Alpha diversity analysis indicated that Chao1, Faith_pd, and Observed species were significantly higher in CA compared to HA (p < 0.01). Furthermore, nine genera were notably different between HA and CA yaks, including Avispirillum, Fimenecus, CAG-1031, Odoribacter 865974, and Jeotgalicoccus A 310962. Compared to CB yaks, CA animals exhibited significant differences in one phylum and six genera, including CAG-485 (p < 0.05), CAG-83 (p < 0.01), Copromorpha (p < 0.01), Phocaeicola A 858004 (p < 0.05), and UBA2253 (p < 0.05).

Conclusion

In summary, higher-weight yak calves demonstrated increased oxidative resistance, and weight profiles were linked to the microbiomes of both female yaks and their calves. These findings offer valuable insights for optimizing yak breeding practices in high-altitude regions.

Engineering TadA ortholog-derived cytosine base editor without motif preference and adenosine activity limitation

The engineered TadA variants used in cytosine base editors (CBEs) present distinctive advantages, including a smaller size and fewer off-target effects compared to cytosine base editors that rely on natural deaminases. However, the current TadA variants demonstrate a preference for base editing in DNA with specific motif sequences and possess dual deaminase activity, acting on both cytosine and adenosine in adjacent positions, limiting their application scope. To address these issues, we employ TadA orthologs screening and multi sequence alignment (MSA)-guided protein engineering techniques to create a highly effective cytosine base editor (aTdCBE) without motif and adenosine deaminase activity limitations. Notably, the delivery of aTdCBE to a humanized mouse model of Duchenne muscular dystrophy (DMD) mice achieves robust exon 55 skipping and restoration of dystrophin expression. Our advancement in engineering TadA ortholog for cytosine editing enriches the base editing toolkits for gene-editing therapy and other potential applications.

PPI-hotspotID for detecting protein-protein interaction hot spots from the free protein structure

Experimental detection of residues critical for protein-protein interactions (PPI) is a time-consuming, costly, and labor-intensive process. Hence, high-throughput PPI-hot spot prediction methods have been developed, but they have been validated using relatively small datasets, which may compromise their predictive reliability. Here, we introduce PPI-hotspotID, a novel method for identifying PPI-hot spots using the free protein structure, and validated it on the largest collection of experimentally confirmed PPI-hot spots to date. We explored the possibility of detecting PPI-hot spots using (i) FTMap in the PPI mode, which identifies hot spots on protein-protein interfaces from the free protein structure, and (ii) the interface residues predicted by AlphaFold-Multimer. PPI-hotspotID yielded better performance than FTMap and SPOTONE, a webserver for predicting PPI-hot spots given the protein sequence. When combined with the AlphaFold-Multimer-predicted interface residues, PPI-hotspotID yielded better performance than either method alone. Furthermore, we experimentally verified several PPI-hotspotID-predicted PPI-hot spots of eukaryotic elongation factor 2. Notably, PPI-hotspotID can reveal PPI-hot spots not obvious from complex structures, including those in indirect contact with binding partners. PPI-hotspotID serves as a valuable tool for understanding PPI mechanisms and aiding drug design. It is available as a web server (https://ppihotspotid.limlab.dnsalias.org/) and open-source code (https://github.com/wrigjz/ppihotspotid/).

The Comparative Full-Length Genome Characterization of African Swine Fever Virus Detected in Thailand

African swine fever virus (ASFV) has been responsible for the globally devastating epidemics in wild and domesticated pigs. Of the 24 identified ASFV genotypes, genotype II is the primary cause for the pandemic occurring in Europe and Asia since its emergence in Georgia in 2007. The current study aimed to characterize the full-length genomic pattern of the ASFV strain from Thailand, TH1_22/CR (Accession No. PP915735), which was then compared with genomic diversity across other Asian isolates using Georgia 2007/1 (Accession No. FR682468) as the reference. Viral DNA was isolated from the pig spleen sample following library preparation and paired-end sequencing using the MiSeq Illumina platform. The sequenced TH1_22/CR isolate spanned 189,395 nucleotides encoding 193 open reading frames (ORFs), exhibiting maximum nucleotide similarity (99.99%) with Georgian (Georgia 2007/1) and Chinese (Wuhan 2019-1 and China HLJ) isolates. Based on phylogenetic analysis, the TH1_22/CR isolate (Accession No. PP915735) was characterized as genotype II, serogroup 8, and IGR-II due to the presence of three tandem repeat sequences (TRSs). Genetic variations including SNPs and single and polynucleotide indels were identified in TH1_22/CR in agreement with other Asian isolates. For comprehensive analysis, the genome was divided into four regions (I-IV) based on gene location. Overall, the TH1_22/CR isolate demonstrated eight SNPs and indels in its genome. Two unique SNPs were reported in the coding regions of the TH1_22/CR isolate, out of which, a C-591-T substitution was seen in MGF 360-4L and a C-297-T was found in A238L, and four unique SNPs were reported in non-coding regions (NCRs). Furthermore, a 29 bp deletion was observed in the IGR between MGF 110-13La and MGF 110-13Lb, as well as 52 bp deletion in the ASFV G ACD 00350 gene. This comparative analysis establishes the foundational information for future studies on the diversity and phylogeography of this regionally significant genetic sub-group of ASFV.

Isolation and Characterization of Plant-Growth-Promoting Bacteria Associated with Salvinia auriculata Aublet

Salvinia auriculata Aublet is a floating aquatic plant, capable of absorbing the excess of nutrients and water contaminants and can be used in effluent treatment plants. The ability to survive in degraded areas may be related to the association with beneficial bacteria capable of promoting plant growth. However, little is known about the microbiota associated with this aquatic plant and its potential application to the aquatic environment. In this sense, this work aims to identify bacteria associated with S. auriculata that could be able to promote plant growth. Eighteen bacterial strains were identified by sequencing of the 16S rRNA gene, belonging to the genera Agrobacterium, Bacillus, Curtobacterium, Enterobacter, Pseudomonas, Siccibacter, and Stenotrophomonas. All isolates produced indole compounds, 12 fixed N2, and 16 solubilized phosphate. A new strain of Enterobacter (sp 3.1.3.0.X.18) was selected for inoculation into S. auriculata. For this purpose, 500 mL of nutrient solution and 1 g of the plant were used in the control and inoculated conditions. Enterobacter inoculation promoted a significant increase (p ≤ 0.05) in fresh plant biomass (17%) after 4 days of cultivation. In summary, the present study characterized 18 plant-growth-promoting bacteria isolated from S. auriculata with potential for biotechnological application, such as the production of bioinoculants or biomass resources, to protect or improve plant growth under conditions of stress.

Medicago truncatula SOBIR1 controls pathogen immunity and specificity in the Rhizobium-legume symbiosis

Medicago truncatula Nod Factor Perception (MtNFP) plays a role in both the Rhizobium-Legume (RL) symbiosis and plant immunity, and evidence suggests that the immune-related function of MtNFP is relevant for symbiosis. To better understand these roles of MtNFP, we sought to identify new interacting partners. We screened a yeast-2-hybrid cDNA library from Aphanomyces euteiches infected and noninfected M. truncatula roots. The M. truncatula leucine-rich repeat (LRR) receptor-like kinase SUPPRESSOR OF BIR1 (MtSOBIR1) was identified as an interactor of MtNFP and was characterised for kinase activity, and potential roles in symbiosis and plant immunity. We showed that the kinase domain of MtSOBIR1 is active and can transphosphorylate the pseudo-kinase domain of MtNFP. MtSOBIR1 could functionally complement Atsobir1 and Nbsobir1/sobir1-like mutants for defence activation, and Mtsobir1 mutants were defective in immune responses to A. euteiches. For symbiosis, we showed that Mtsobir1 mutant plants had both a strong, early infection defect and defects in the defence suppression in nodules, and both effects were plant genotype- and rhizobial strain-specific. This work highlights a conserved function for MtSOBIR1 in activating defence responses to pathogen attack, and potentially novel symbiotic functions of downregulating defence in association with the control of symbiotic specificity.

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.

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.

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.

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.

Efficient inference of large prokaryotic pangenomes with PanTA

Pangenome inference is an indispensable step in bacterial genomics, yet its scalability poses a challenge due to the rapid growth of genomic collections. This paper presents PanTA, a software package designed for constructing pangenomes of large bacterial datasets, showing unprecedented efficiency levels multiple times higher than existing tools. PanTA introduces a novel mechanism to construct the pangenome progressively without rebuilding the accumulated collection from scratch. The progressive mode is shown to consume orders of magnitude less computational resources than existing solutions in managing growing datasets. The software is open source and is publicly available at https://github.com/amromics/panta and at 10.6084/m9.figshare.23724705 .

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.

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).

Phototrophs as the central components of the plastisphere microbiome in coastal environments

Upon entering the marine environment, plastics are colonized by a plethora of microorganisms to form a plastisphere, influencing the fate and transport of the plastic debris and the health of marine ecosystems. The assembly of marine plastisphere is generally believed to be dominated by stochastic processes. However, it remains elusive whether microbial interaction in the assembly of plastisphere microbial communities is conserved or not. We analyzed the plastisphere microbiomes of 137 plastic debris samples from intertidal zones at different geographical locations and habitats (seagrass, coral, mangrove, beach, and open ocean) and compared them with the surrounding sediment and seawater microbiomes. Microbial community structures of the plastisphere from different locations were more similar to each other but differed substantially from the surrounding sediment and water microbiomes, implying a common mechanism of plastisphere assembly. We used different machine learning algorithms (Multinomial Logistic Regression, Support Vector Machine, Decision Trees, Random Forest, and Artificial Neural Networks) to classify plastic debris samples with high sensitivity based on the microbiome composition. Eukaryotic and prokaryotic phototrophic organisms such as green algae, diatoms, and cyanobacteria, were found to be enriched on the plastic surfaces. Network analysis revealed the central role of the phototrophic organisms in the formation and sustenance of the plastispheres. We found that phototrophs served as core members interacting strongly with heterotrophic organisms in marine plastisphere, irrespective of the sampling location, habitats, and polymer types. This would explain the stochastic assembly of the plastisphere along with conserved properties driven by the phototrophs in the surrounding environment. Our results highlight the importance of phototrophic organisms in shaping the marine plastisphere microbial communities.

The first complete mitochondrial genome of Grossulariaceae: Molecular features, structure recombination, and genetic evolution

BACKGROUND

Mitochondria play crucial roles in the growth, development, and adaptation of plants. Blackcurrant (Ribes nigrum L.) stands out as a significant berry species due to its rich nutritional profile, medicinal properties, and health benefits. Despite its importance, the mitochondrial genome of blackcurrant remains unassembled.

RESULTS

This study presents the first assembly of the mitochondrial genome of R. nigrum in the Grossulariaceae family. The genome spans 450,227 base pairs (bp) and encompasses 39 protein-coding genes (PCGs), 19 transfer RNAs (tRNAs), and three ribosomal RNAs (rRNAs). Protein-coding regions constitute 8.88% of the entire genome. Additionally, we identified 180 simple sequence repeats, 12 tandem repeats, and 432 pairs of dispersed repeats. Notably, the dispersed sequence R1 (cotig3, 1,129 bp) mediated genome recombination, resulting in the formation of two major conformations, namely master and double circles. Furthermore, we identified 731 C-to-U RNA editing sites within the PCGs. Among these, cox1-2, nad1-2, and nad4L-2 were associated with the creation of start codons, whereas atp6-718 and rps10-391 were linked to termination codons. We also detected fourteen plastome fragments within the mitogenome, constituting 1.11% of the total length. Phylogenetic analysis suggests that R. nigrum might have undergone multiple genomic reorganization and/or gene transfer events, resulting in the loss of two PCGs (rps2 and rps11) during its evolutionary history.

CONCLUSIONS

This investigation unveils the molecular characteristics of the R. nigrum mitogenome, shedding light on its evolutionary trajectory and phylogenetic implications. Furthermore, it serves as a valuable reference for evolutionary research and germplasm identification within the genus.

Phylogenetic Conflict Between Species Tree and Maternally Inherited Gene Trees in a Clade of Emberiza Buntings (Aves: Emberizidae)

Different genomic regions may reflect conflicting phylogenetic topologies primarily due to incomplete lineage sorting and/or gene flow. Genomic data are necessary to reconstruct the true species tree and explore potential causes of phylogenetic conflict. Here, we investigate the phylogenetic relationships of 4 Emberiza species (Aves: Emberizidae) and discuss the potential causes of the observed mitochondrial non-monophyly of Emberiza godlewskii (Godlewski's bunting) using phylogenomic analyses based on whole genome resequencing data from 41 birds. Analyses based on both the whole mitochondrial genome and ~39 kilobases from the non-recombining W chromosome reveal sister relationships between each the northern and southern populations of E. godlewskii with E. cioides and E. cia, respectively. In contrast, the monophyly of E. godlewskii is reflected by the phylogenetic signal of autosomal and Z chromosomal sequence data as well as demographic inference analyses, which-in combination-support the following tree topology: ([{E. godlewskii, E. cia}, E. cioides], E. jankowskii). Using D-statistics, we detected multiple gene flow events among different lineages, indicating pervasive introgressive hybridization within this clade. Introgression from an unsampled lineage that is sister to E. cioides or introgression from an unsampled mitochondrial + W chromosomal lineage of E. cioides into northern E. godlewskii may explain the phylogenetic conflict between the species tree estimated from genome-wide data versus mtDNA/W tree topologies. These results underscore the importance of using genomic data for phylogenetic reconstruction and species delimitation.

AMRomics: a scalable workflow to analyze large microbial genome collections

Whole genome analysis for microbial genomics is critical to studying and monitoring antimicrobial resistance strains. The exponential growth of microbial sequencing data necessitates a fast and scalable computational pipeline to generate the desired outputs in a timely and cost-effective manner. Recent methods have been implemented to integrate individual genomes into large collections of specific bacterial populations and are widely employed for systematic genomic surveillance. However, they do not scale well when the population expands and turnaround time remains the main issue for this type of analysis. Here, we introduce AMRomics, an optimized microbial genomics pipeline that can work efficiently with big datasets. We use different bacterial data collections to compare AMRomics against competitive tools and show that our pipeline can generate similar results of interest but with better performance. The software is open source and is publicly available at https://github.com/amromics/amromics under an MIT license.

Jamaican fruit bat (Artibeus jamaicensis) insusceptibility to mucosal inoculation with SARS-CoV-2 Delta variant is not caused by receptor compatibility

The ancestral sarbecovirus giving rise to SARS-CoV-2 is posited to have originated in bats. While SARS-CoV-2 causes asymptomatic to severe respiratory disease in humans, little is known about the biology, virus tropism, and immunity of SARS-CoV-2-like sarbecoviruses in bats. SARS-CoV-2 has been shown to infect multiple mammalian species, including various rodent species, non-human primates, and Egyptian fruit bats. We show that SARS-CoV-2 can utilize Jamaican fruit bat (Artibeus jamaicensis) ACE2 spike for entry in vitro. Therefore, we investigate the Jamaican fruit bat as a possible in vivo model to study reservoir responses. We find that SARS-CoV-2 Delta does not efficiently replicate in Jamaican fruit bats in vivo. We observe infectious viruses in the lungs of only one animal on day 1 post-inoculation and find no evidence of shedding or seroconversion. This is possibly due to host factors restricting virus egress after aborted replication. Furthermore, we observe no significant immune gene expression changes in the respiratory tract but do observe changes in the intestinal metabolome after inoculation. This suggests that, despite its broad host range, SARS-CoV-2 is unable to infect all bat species, and Jamaican fruit bats are not an appropriate model to study SARS-CoV-2 reservoir infection.

Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae

The mitochondrial genomes of three species of Yangiella were sequenced, annotated, and analyzed. The genome length of the three species of the genus is 15,070-15,202 bp, with a typical gene number, including a control region, 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs). It was found that the mitochondrial genome of Yangiella had AT bias. Except for the lack of a DHU arm of the trnS1 gene, the other tRNAs had a typical cloverleaf structure, and the codon usage preferences of the three species exhibited high similarity. In addition, tRNA gene rearrangements were observed among the three subfamilies of Aradidae (Mezirinae, Calisiinae, Aradinae), and it was found that codon usage preferences appeared to be less affected by base mutation and more by natural selection. The Pi and Ka/Ks values indicated that cox1 was the most conserved gene in the mitochondrial genome of Aradidae, while atp8 and nad6 were rapidly evolved genes. Substitution saturation level analysis showed that the nucleic acid sequence of mitochondrial protein-coding genes in Aradidae did not reach saturation, suggesting the rationality of the phylogenetic analysis data. Bayesian and maximum likelihood methods were used to analyze the phylogeny of 16 species of Hemiptera insects, which supported the monophyly of Aneurinae, Carventinae, and Mezirinae, as well as the monophyly of Yangiella. Based on fossils and previous studies, the differentiation time was inferred, indicating that Yangiella diverged about 57 million years ago.

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.

Aurantiacibacter hainanensis sp. nov. and Qipengyuania zhejiangensis sp. nov., two novel Erythrobacteraceae species isolated from tidal flat sediments

Two Gram-stain-negative, rod-shaped, non-motile, aerobic and carotenoid-producing strains, belonging to the family Erythrobacteraceae, designated as H149T and Z2T, were isolated from tidal flat sediment samples collected in Hainan and Zhejiang, PR China, respectively. Growth of strain H149T occurred at 15-42 °C, 0-10.0 % (w/v) NaCl, and pH 6.0-8.5, with the optima at 35-37 °C, 3.0-3.5 % (w/v) NaCl and pH 7.0. Strain Z2T grew at 15-37 °C, 0-6.0 % (w/v) NaCl, and pH 6.0-9.5, with the optima at 25-30 °C, 0.5-1.0 % (w/v) NaCl and pH 6.0-6.5. Ubiquinone-10 was the sole ubiquinone in two strains. The predominant cellular fatty acids of strain H149T were C16 : 0, summed feature 3 and summed feature 8, while those of strain Z2T were C17 : 1  ω6c, summed feature 3 and summed feature 8. Strains H149T and Z2T shared diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, phosphatidylcholine and sphingoglycolipid as major polar lipids. The 16S rRNA gene sequence identity analysis indicated that strain H149T had the highest sequence identity of 98.4 % with Aurantiacibacter odishensis KCTC 23981T, and strain Z2T had that of 98.2 % with Qipengyuania pacifica NZ-96T. Phylogenetic trees based on 16S rRNA gene and core-genome sequences revealed that strains H149T and Z2T formed two independent clades in the genera Aurantiacibacter and Qipengyuania, respectively. Strain H149T had average nucleotide identity values of 74.0-81.3 % and in silico DNA-DNA hybridization values of 18.5-23.1 % with Aurantiacibacter type strains, while strain Z2T had values of 73.3-78.7 % and 14.5-33.3 % with Qipengyuania type strains. The genomic DNA G+C contents of strains H149T and Z2T were 64.3 and 61.8 %, respectively. Based on the genetic, genomic, phylogenetic, physiological and chemotaxonomic results, strains H149T (=KCTC 8397T=MCCC 1K08920T) and Z2T (=KCTC 8396T=MCCC 1K08946T) are concluded to represent two novel Erythrobacteraceae species for which the names Aurantiacibacter hainanensis sp. nov. and Qipengyuania zhejiangensis sp. nov. are proposed, respectively.

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.

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.