The COG database: a tool for genome-scale analysis of protein functions and evolution

Rational construction of genome-minimized Streptomyces host for the expression of secondary metabolite gene clusters

Streptomyces offer a wealth of naturally occurring compounds with diverse structures, many of which possess significant pharmaceutical values. However, new product exploration and increased yield of specific compounds in Streptomyces have been technically challenging due to their slow growth rate, complex culture conditions and intricate genetic backgrounds. In this study, we screened dozens of Streptomyces strains inhabiting in a plant rhizosphere for fast-growing candidates, and further employed CRISPR/Cas-based engineering techniques for stepwise refinement of a particular strain, Streptomyces sp. A-14 that harbors a 7.47 Mb genome. After strategic removal of nonessential genomic regions and most gene clusters, we reduced its genome size to 6.13 Mb, while preserving its growth rate to the greatest extent. We further demonstrated that cleaner metabolic background of this engineered strain was well suited for the expression and characterization of heterologous gene clusters, including the biosynthetic pathways of actinorhodin and polycyclic tetramate macrolactams. Moreover, this streamlined genome is anticipated to facilitate directing the metabolic flux towards the production of desired compounds and increasing their yields.

An upstream signaling gene calmodulin regulates the synthesis of insect wax via activating fatty acid biosynthesis pathway

Insect wax accumulates on the surface of insect cuticle, which acts as an important protective barrier against rain, ultraviolet light radiation, pathogens, etc. The waxing behavior, wax composition and molecular mechanism underling wax biosynthesis are unclear in dustywings. Herein, the current study determined the vital developmental stage for waxing behavior in dustywings, examined the components of waxy secretions, and identified key regulatory genes for wax biosynthesis. The wax glands were mainly located on the thorax and abdomen of dustywing adults. The adults spread the waxy secretions over their entire body surface. The metabolomics analysis identified 32 lipids and lipid-like molecules, 15 organic acids and derivatives, 7 benzenoids, etc. as the main components of waxy secretions. The fatty acids represented the largest proportion of the category of lipid and lipid-like molecules. The conjoint analysis of metabolomics and transcriptomics identified two crucial genes fatty acyl-CoA reductase (CsFAR) and calmodulin (CsCaM) for wax biosynthesis. The down-regulation of these genes via nanocarrier-mediated RNA interference technology significantly reduced the amount of wax particles. Notably, the RNAi of CsCaM apparently suppressed the expression of most genes in fatty acid biosynthesis pathway, indicating the CsCaM might act as a main upstream regulator of fatty acid biosynthesis pathway.

Genome analysis of Rossellomorea sp. y25, a deep sea bacterium isolated from the sediments of South China Sea

Rossellomorea sp. y25, a putative new species of yellow pigment-producing, aerobic and chemoheterotrophic bacterium belonging to the family Bacillaceae, was isolated from the sediments at the depth of 1829 m in the South China Sea. In this study, we present the complete genome sequences of strain y25, which consisted of only one circular chromosome with 4,633,006 bp and the content of G + C was 41.76%. A total of 4466 CDSs, 106 tRNA, 33 rRNA, and 101 sRNA genes were obtained. Genomic analysis of strain y25 showed that it has the ability to produce antioxidant carotenoids and a large number of heavy metal resistance genes, such as arsenic, cadmium and zinc. In addition, strain y25 contains a prophage that may contribute to host protection against lysis by related Bacillus-like phages. This is the first report of genome-wide information on a bacterium of the genus Rossellomorea isolated from the deep sea, providing insights into how microorganisms of this genus adapt to deep-sea environments.

Environmental magnesium ion affects global gene expression, motility, biofilm formation and virulence of Vibrio parahaemolyticus

Vibrio parahaemolyticus is widely distributed in marine ecosystems. Magnesium ion (Mg2+) is the second most abundant metal cation in seawater, and plays important roles in the growth and gene expression of V. parahaemolyticus, but lacks the detailed mechanisms. In this study, the RNA sequencing data demonstrated that a total of 1494 genes was significantly regulated by Mg2+. The majority of the genes associated with lateral flagella, exopolysaccharide, type III secretion system 2, type VI secretion system (T6SS) 1, T6SS2, and thermostable direct hemolysin were downregulated. A total of 18 genes that may be involved in c-di-GMP metabolism and more than 80 genes encoding putative regulators were also significantly and differentially expressed in response to Mg2+, indicating that the adaptation process to Mg2+ stress may be strictly regulated by complex regulatory networks. In addition, Mg2+ promoted the proliferative speed, swimming motility and cell adhesion of V. parahaemolyticus, but inhibited the swarming motility, biofilm formation, and c-di-GMP production. However, Mg2+ had no effect on the production of capsular polysaccharide and cytoxicity against HeLa cells. Therefore, Mg2+ had a comprehensive impact on the physiology and gene expression of V. parahaemolyticus.

Exploring the impact of breast cancer on colonization resistance of mouse microbiota using network node manipulation

Breast cancer, a global health concern affecting women, has been linked to alterations in the gut microbiota, impacting various aspects of human health. This study investigates the interplay between breast cancer and the gut microbiome, particularly focusing on colonization resistance-an essential feature of the microbiota's ability to prevent pathogenic overgrowth. Using a mouse model of breast cancer, we employ diversity analysis, co-occurrence network analysis, and robustness tests to elucidate the impact of breast cancer on microbiome dynamics. Our results reveal that breast cancer exposure affects the bacterial community's composition and structure, with temporal dynamics playing a role. Network analysis demonstrates that breast cancer disrupts microbial interactions and decreases network complexity, potentially compromising colonization resistance. Moreover, network robustness analysis shows the susceptibility of the microbiota to node removal, indicating potential vulnerability to pathogenic colonization. Additionally, predicted metabolic profiling of the microbiome highlights the significance of the enzyme EC 6.2.1.2 - Butyrate--CoA ligase, potentially increasing butyrate, and balancing the reduction of colonization resistance. The identification of Rubrobacter as a key contributor to this enzyme suggests its role in shaping the microbiota's response to breast cancer. This study uncovers the intricate relationship between breast cancer, the gut microbiome, and colonization resistance, providing insights into potential therapeutic strategies and diagnostic approaches for breast cancer patients.

Sugarcane ScCAX4 is a Negative Regulator of Resistance to Pathogen Infection

Calcium (Ca2+) is a second messenger in various physiological processes within plants. The significance of the Ca2+/H+ exchanger (CAX) has been established in facilitating Ca2+ transport in plants; however, disease resistance functions of the CAX gene remain elusive. In this study, we conducted sequence characterization and expression analysis for a sugarcane CAX gene, ScCAX4 (GenBank Accession Number: MW206380). In order to further investigate the disease resistance functions, this gene was then transiently overexpressed in Nicotiana benthamiana leaves, which were subsequently inoculated with Fusarium solani var. coeruleum. Results showed that ScCAX4 overexpression increased the susceptibility of N. benthamiana to pathogen infection by regulating the expression of genes related to salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) pathways, suggesting its negative role in disease resistance. Furthermore, we genetically transformed the ScCAX4 gene into N. benthamiana and obtained three positive T2 generation lines. Interestingly, the symptomatology of transgenic plants was consistent with that of transient overexpression after pathogen inoculation. Notably, the JA content in transgenic overexpression lines was significantly higher than that in the wild-type. RNA-seq revealed that ScCAX4 could mediate multiple signaling pathways, and the JA signaling pathway played a key role in modulating disease resistance. Finally, a regulatory model was depicted for the increased susceptibility to pathogen infection conferred by the ScCAX4 gene. This study provides genetic resources for sugarcane molecular breeding and the research direction for plant CAX genes.

Core proteome mediated subtractive approach for the identification of potential therapeutic drug target against the honeybee pathogen Paenibacillus larvae

Background & Objectives

American foulbrood (AFB), caused by the highly virulent, spore-forming bacterium Paenibacillus larvae, poses a significant threat to honey bee brood. The widespread use of antibiotics not only fails to effectively combat the disease but also raises concerns regarding honey safety. The current computational study was attempted to identify a novel therapeutic drug target against P. larvae, a causative agent of American foulbrood disease in honey bee.

Methods

We investigated effective novel drug targets through a comprehensive in silico pan-proteome and hierarchal subtractive sequence analysis. In total, 14 strains of P. larvae genomes were used to identify core genes. Subsequently, the core proteome was systematically narrowed down to a single protein predicted as the potential drug target. Alphafold software was then employed to predict the 3D structure of the potential drug target. Structural docking was carried out between a library of phytochemicals derived from traditional Chinese flora (n > 36,000) and the potential receptor using Autodock tool 1.5.6. Finally, molecular dynamics (MD) simulation study was conducted using GROMACS to assess the stability of the best-docked ligand.

Results

Proteome mining led to the identification of Ketoacyl-ACP synthase III as a highly promising therapeutic target, making it a prime candidate for inhibitor screening. The subsequent virtual screening and MD simulation analyses further affirmed the selection of ZINC95910054 as a potent inhibitor, with the lowest binding energy. This finding presents significant promise in the battle against P. larvae.

Conclusions

Computer aided drug design provides a novel approach for managing American foulbrood in honey bee populations, potentially mitigating its detrimental effects on both bee colonies and the honey industry.

Study on the causes of changes in colour during Hibiscus syriacus flowering based on transcriptome and metabolome analyses

BACKGROUND

The flower colour of H. syriacus 'Qiansiban' transitions from fuchsia to pink-purple and finally to pale purple, thereby enhancing the ornamental value of the cultivars. However, the molecular mechanism underlying this change in flower colour in H. syriacus has not been elucidated. In this study, the transcriptomic data of H. syriacus 'Qiansiban' at five developmental stages were analysed to investigate the impact of flavonoid components on flower colour variation. Additionally, five cDNA libraries were constructed from H. syriacus 'Qiansiban' during critical blooming stages, and the transcriptomes were sequenced to investigate the molecular mechanisms underlying changes in flower colouration.

RESULTS

High-performance liquid chromatography‒mass spectrometry detected five anthocyanins in H. syriacus 'Qiansiban', with malvaccin-3-O-glucoside being the predominant compound in the flowers of H. syriacus at different stages, followed by petunigenin-3-O-glucoside. The levels of these five anthocyanins exhibited gradual declines throughout the flowering process. In terms of the composition and profile of flavonoids and flavonols, a total of seven flavonoids were identified: quercetin-3-glucoside, luteolin-7-O-glucoside, Santianol-7-O-glucoside, kaempferol-O-hexosyl-C-hexarbonoside, apigenin-C-diglucoside, luteolin-3,7-diglucoside, and apigenin-7-O-rutinoside. A total of 2,702 DEGs were identified based on the selected reference genome. Based on the enrichment analysis of differentially expressed genes, we identified 9 structural genes (PAL, CHS, FLS, DRF, ANS, CHI, F3H, F3'5'H, and UFGT) and 7 transcription factors (3 MYB, 4 bHLH) associated with flavonoid biosynthesis. The qRT‒PCR results were in good agreement with the high-throughput sequencing data.

CONCLUSION

This study will establish a fundamental basis for elucidating the mechanisms underlying alterations in the flower pigmentation of H. syriacus.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. In this manuscript, we present a genome for a laboratory-tractable Entomophthora muscae isolate that infects fruit flies. Our E. muscae assembly is 1.03 Gb, consists of 7810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with recent datasets from entomophthoralean fungi, we show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity. In addition, we find that E. muscae and its closest allies possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Finally, we offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex that support the need for taxonomic revision within this group. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Function of MYB8 in larch under PEG simulated drought stress

Larch, a prominent afforestation, and timber species in northeastern China, faces growth limitations due to drought. To further investigate the mechanism of larch's drought resistance, we conducted full-length sequencing on embryonic callus subjected to PEG-simulated drought stress. The sequencing results revealed that the differentially expressed genes (DEGs) primarily played roles in cellular activities and cell components, with molecular functions such as binding, catalytic activity, and transport activity. Furthermore, the DEGs showed significant enrichment in pathways related to protein processing, starch and sucrose metabolism, benzose-glucuronic acid interconversion, phenylpropyl biology, flavonoid biosynthesis, as well as nitrogen metabolism and alanine, aspartic acid, and glutamic acid metabolism. Consequently, the transcription factor T_transcript_77027, which is involved in multiple pathways, was selected as a candidate gene for subsequent drought stress resistance tests. Under PEG-simulated drought stress, the LoMYB8 gene was induced and showed significantly upregulated expression compared to the control. Physiological indices demonstrated an improved drought resistance in the transgenic plants. After 48 h of PEG stress, the transcriptome sequencing results of the transiently transformed LoMYB8 plants and control plants exhibited that genes were significantly enriched in biological process, cellular component and molecular function. Function analyses indicated for the enrichment of multiple KEGG pathways, including energy synthesis, metabolic pathways, antioxidant pathways, and other relevant processes. The pathways annotated by the differential metabolites mainly encompassed signal transduction, carbohydrate metabolism, amino acid metabolism, and flavonoid metabolism.

Chromosome-level Alstonia scholaris genome unveils evolutionary insights into biosynthesis of monoterpenoid indole alkaloids

Alstonia scholaris of the Apocynaceae family is a medicinal plant with a rich source of bioactive monoterpenoid indole alkaloids (MIAs), which possess anti-cancer activity like vinca alkaloids. To gain genomic insights into MIA biosynthesis, we assembled a high-quality chromosome-level genome for A. scholaris using nanopore and Hi-C data. The 444.95 Mb genome contained 35,488 protein-coding genes. A total of 20 chromosomes were assembled with a scaffold N50 of 21.75 Mb. The genome contained a cluster of strictosidine synthases and tryptophan decarboxylases with synteny to other species and a saccharide-terpene cluster involved in the monoterpenoid biosynthesis pathway of the MIA upstream pathway. The multi-omics data of A. scholaris provide a valuable resource for understanding the evolutionary origins of MIAs and for discovering biosynthetic pathways and synthetic biology efforts for producing pharmaceutically useful alkaloids.

Extracellular Vesicles From Mycoplasma gallisepticum: Modulators of Macrophage Activation and Virulence

Extracellular vesicles (EVs) mediate intercellular communication by transporting proteins. To investigate the pathogenesis of Mycoplasma gallisepticum, a major threat to the poultry industry, we isolated and characterized M. gallisepticum-produced EVs. Our study highlights the significant impact of M. gallisepticum-derived EVs on immune function and macrophage apoptosis, setting them apart from other M. gallisepticum metabolites. These EVs dose-dependently enhance M. gallisepticum adhesion and proliferation, simultaneously modulating Toll-like receptor 2 and interferon γ pathways and thereby inhibiting macrophage activation. A comprehensive protein analysis revealed 117 proteins in M. gallisepticum-derived EVs, including established virulence factors, such as GapA, CrmA, VlhA, and CrmB. Crucially, these EV-associated proteins significantly contribute to M. gallisepticum infection. Our findings advance our comprehension of M. gallisepticum pathogenesis, offering insights for preventive strategies and emphasizing the pivotal role of M. gallisepticum-derived EVs and their associated proteins. This research sheds light on the composition and crucial role of M. gallisepticum-derived EVs in M. gallisepticum pathogenesis, aiding our fight against M. gallisepticum infections.

Toward Characterizing Environmental Sources of Non-tuberculous Mycobacteria (NTM) at the Species Level: A Tutorial Review of NTM Phylogeny and Phylogenetic Classification

Nontuberculous mycobacteria (NTM) are any mycobacteria that do not cause tuberculosis or leprosy. While the majority of NTM are harmless and some of them are considered probiotic, a growing number of people are being diagnosed with NTM infections. Therefore, their detection in the environment is of interest to clinicians, environmental microbiologists, and water quality researchers alike. This review provides a tutorial on the foundational approaches for taxonomic classifications, with a focus on the phylogenetic relationships among NTM revealed by the 16S rRNA gene, rpoB gene, and hsp65 gene, and by genome-based approaches. Recent updates on the Mycobacterium genus taxonomy are also provided. A synthesis on the habitats of 189 mycobacterial species in a genome-based taxonomy framework was performed, with attention paid to environmental sources (e.g., drinking water, aquatic environments, and soil). The 16S rRNA gene-based classification accuracy for various regions was evaluated (V3, V3-V4, V3-V5, V4, V4-V5, and V1-V9), revealing overall excellent genus-level classification (up to 100% accuracy) yet only modest performance (up to 63.5% accuracy) at the species level. Future research quantifying NTM species in water systems, determining the effects of water treatment and plumbing conditions on their variations, developing high throughput species-level characterization tools for use in the environment, and incorporating the characterization of functions in a phylogenetic framework will likely fill critical knowledge gaps. We believe this tutorial will be useful for researchers new to the field of molecular or genome-based taxonomic profiling of environmental microbiomes. Experts may also find this review useful in terms of the selected key findings of the past 30 years, recent updates on phylogenomic analyses, as well as a synthesis of the ecology of NTM in a phylogenetic framework.

Synergy of lytic phage pB23 and meropenem combination against carbapenem-resistant Acinetobacter baumannii

Phage-antibiotic combination treatment is a novel noteworthy drug delivery method in anti-infection. In the current study, we have isolated a new phage, pB23, against carbapenem-resistant Acinetobacter baumannii 2023. Synergistic antibacterial effect between phage pB23 and meropenem combination could be more stable, using moderate doses of phage (multiplicity of infection ranging from 0.1 to 1,000) based on results of in vitro antibacterial activity. Phage pB23 and meropenem combination could effectively clear mature biofilms and prevent biofilm formation of carbapenem-resistant Acinetobacter baumannii in vitro. Phage pB23 and meropenem combination also has good synergistic antibacterial effects against carbapenem-resistant Acinetobacter baumannii in different growth phases under static culture conditions. The pig skin explant model shows that phage pB23 and meropenem combination has a synergistic effect to remove bacteria from wounds ex vivo. Phage pB23 and meropenem combination also exhibited a synergistic antibacterial effect in vivo using a zebrafish infection mode. The potential promotion of phage proliferation by meropenem and the sensitivity recovery of phage-resistant bacteria to meropenem might elucidate the mechanism of the synergistic antimicrobial activity. In summary, our study illustrates that phage pB23 and meropenem combination could produce synergistic antibacterial effects against carbapenem-resistant Acinetobacter baumannii under static growth conditions. This study also demonstrates that phage-antibiotic combination will become an effective strategy to enhance antibacterial activity of individual drug and provide a new idea of the drug development for the treatment of infections due to carbapenem-resistant Acinetobacter baumannii and other multidrug-resistant bacteria.

Genome-Wide Identification and Characterization of Homeobox Transcription Factors in Phoma sorghina var. saccharum Causing Sugarcane Twisted Leaf Disease

A homeobox transcription factor is a conserved transcription factor, ubiquitous in eukaryotes, that regulates the tissue formation of structure, cell differentiation, proliferation, and cancer. This study identified the homeobox transcription factor family and its distribution in Phoma sorghina var. saccharum at the whole genome level. It elucidated the gene structures and evolutionary characteristics of this family. Additionally, knockout experiments were carried out and the preliminary function of these transcription factors was studied. Through bioinformatics approaches, nine homeobox transcription factors (PsHOX1-PsHOX9) were identified in P. sorghina var. saccharum, and these contained HOX-conserved domains and helix-turn-helix secondary structures. Nine homeobox gene deletion mutants were obtained using the homologous recombinant gene knockout technique. Protoplast transformation was mediated by polyethylene glycol (PEG) and the transformants were identified using PCR. The knockouts of PsHOX1, PsHOX2, PsHOX3, PsHOX4, PsHOX6, PsHOX8, and PsHOX9 genes resulted in a smaller growth diameter in P. sorghina var. saccharum. In contrast, the knockouts of the PsHOX3, PsHOX6, and PsHOX9 genes inhibited the formation of conidia and led to a significant decrease in the pathogenicity. This study's results will provide insights for understanding the growth and development of P. sorghina var. saccharum. The pathogenic mechanism of the affected sugarcane will provide an essential theoretical basis for preventing and controlling sugarcane twisted leaf disease.

The gills and skin microbiota of five pelagic fish species from the Atlantic Ocean

The gills and skin microbiota and microbiome of wild fish remain far more under-investigated compared to that of farmed fish species, despite that these animal-microbe interactions hold the same ecophysiological roles in both cases. In this study, the gills and skin bacterial microbiota profiles and their presumptive bacterial metabolisms were investigated in five open-sea fishes: bullet tuna (Auxis sp.), common dolphinfish (Coryphaena hippurus), Atlantic little tunny (Euthynnus alletteratus), Atlantic bonito (Sarda sarda) and Atlantic white marlin (Kajikia albida). Gills and skin tissues were collected from two to three individuals per species, from specimens caught by recreational trolling during summer of 2019, and their bacterial 16S rRNA gene diversity was analysed by high-throughput sequencing. The gills bacterial communities among the five species were clearly different but not the skin bacterial microbiota. The dominant operational taxonomic units belonged to the Moraxellaceae, Pseudomonadaceae, Rhodobacteraceae, Staphylococcaceae and Vibrionaceae families. Despite the differences in taxonomic composition, the presumptive bacterial metabolisms between the gills and skin of the five fishes investigated here were ≥ 94% similar and were dominated by basic metabolism, most likely reflecting the continuous exposure of these tissues in the surrounding seawater.

Transcriptome dynamics of Gossypium purpurascens in response to abiotic stresses by Iso-seq and RNA-seq data

Gossypium purpurascens is a member of the Malvaceae family, holds immense economic significance as a fiber crop worldwide. Abiotic stresses harm cotton crops, reduce yields, and cause economic losses. Generating high-quality reference genomes and large-scale transcriptomic datasets across diverse conditions can offer valuable insights into identifying preferred agronomic traits for crop breeding. The present research used leaf tissues to conduct PacBio Iso-seq and RNA-seq analysis. We carried out an in-depth analysis of DEGs using both correlations with cluster analysis and principal component analysis. Additionally, the study also involved the identification of both lncRNAs and CDS. We have prepared RNA-seq libraries from 75 RNA samples to study the effects of drought, salinity, alkali, and saline-alkali stress, as well as control conditions. A total of 454.06 Gigabytes of transcriptome data were effectively validated through the identification of differentially expressed genes and KEGG and GO analysis. Overwhelmingly, gene expression profiles and full-length transcripts from cotton tissues will aid in understanding the genetic mechanism of abiotic stress tolerance in G. purpurascens.

A metagenomics pipeline reveals insertion sequence-driven evolution of the microbiota

Insertion sequence (IS) elements are mobile genetic elements in bacterial genomes that support adaptation. We developed a database of IS elements coupled to a computational pipeline that identifies IS element insertions in the microbiota. We discovered that diverse IS elements insert into the genomes of intestinal bacteria regardless of human host lifestyle. These insertions target bacterial accessory genes that aid in their adaptation to unique environmental conditions. Using IS expansion in Bacteroides, we show that IS activity leads to the insertion of "hot spots" in accessory genes. We show that IS insertions are stable and can be transferred between humans. Extreme environmental perturbations force IS elements to fall out of the microbiota, and many fail to rebound following homeostasis. Our work shows that IS elements drive bacterial genome diversification within the microbiota and establishes a framework for understanding how strain-level variation within the microbiota impacts human health.

Characterization of clumpy adhesion of Escherichia coli to human cells and associated factors influencing antibiotic sensitivity

Escherichia coli intestinal infection pathotypes are characterized by distinct adhesion patterns, including the recently described clumpy adhesion phenotype. Here, we identify and characterize the genetic factors contributing to the clumpy adhesion of E. coli strain 4972. In this strain, the transcriptome and proteome of adhered bacteria were found to be distinct from planktonic bacteria in the supernatant. A total of 622 genes in the transcriptome were differentially expressed in bacteria present in clumps relative to the planktonic bacteria. Seven genes targeted for disruption had variable distribution in different pathotypes and nonpathogenic E. coli, with the pilV and spnT genes being the least frequent or absent from most groups. Deletion (Δ) of five differentially expressed genes, flgH, ffp, pilV, spnT, and yggT, affected motility, adhesion, or antibiotic stress. ΔflgH exhibited 80% decrease and ΔyggT depicted 184% increase in adhesion, and upon complementation, adhesion was significantly reduced to 13%. ΔflgH lost motility and was regenerated when complemented, whereas Δffp had significantly increased motility, and reintroduction of the same gene reduced it to the wild-type level. The clumps produced by Δffp and ΔspnT were more resistant and protected the bacteria, with ΔspnT showing the best clump formation in terms of ampicillin stress protection. ΔyggT had the lowest tolerance to gentamicin, where the antibiotic stress completely eliminated the bacteria. Overall, we were able to investigate the influence of clump formation on cell surface adhesion and antimicrobial tolerance, with the contribution of several factors crucial to clump formation on susceptibility to the selected antibiotics.

IMPORTANCE

The study explores a biofilm-like clumpy adhesion phenotype in Escherichia coli, along with various factors and implications for antibiotic susceptibility. The phenotype permitted the bacteria to survive the onslaught of high antibiotic concentrations. Profiles of the transcriptome and proteome allowed the differentiation between adhered bacteria in clumps and planktonic bacteria in the supernatant. The deletion mutants of genes differentially expressed between adhered and planktonic bacteria, i.e., flgH, ffp, pilV, spnT, and yggT, and respective complementations in trans cemented their roles in multiple capacities. ffp, an uncharacterized gene, is involved in motility and resistance to ampicillin in a clumpy state. The work also affirms for the first time the role of the yggT gene in adhesion and its involvement in susceptibility against another aminoglycoside antibiotic, i.e., gentamicin. Overall, the study contributes to the mechanisms of biofilm-like adhesion phenotype and understanding of the antimicrobial therapy failures and infections of E. coli.

Critical steps in an environmental metaproteomics workflow

Environmental metaproteomics is a rapidly advancing field that provides insights into the structure, dynamics, and metabolic activity of microbial communities. As the field is still maturing, it lacks consistent workflows, making it challenging for non-expert researchers to navigate. This review aims to introduce the workflow of environmental metaproteomics. It outlines the standard practices for sample collection, processing, and analysis, and offers strategies to overcome the unique challenges presented by common environmental matrices such as soil, freshwater, marine environments, biofilms, sludge, and symbionts. The review also highlights the bottlenecks in data analysis that are specific to metaproteomics samples and provides suggestions for researchers to obtain high-quality datasets. It includes recent benchmarking studies and descriptions of software packages specifically built for metaproteomics analysis. The article is written without assuming the reader's familiarity with single-organism proteomic workflows, making it accessible to those new to proteomics or mass spectrometry in general. This primer for environmental metaproteomics aims to improve accessibility to this exciting technology and empower researchers to tackle challenging and ambitious research questions. While it is primarily a resource for those new to the field, it should also be useful for established researchers looking to streamline or troubleshoot their metaproteomics experiments.

Comparative genomic analysis and functional investigations for MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria in ecology

Microcystins (MCs) significantly threaten the ecosystem and public health. Biodegradation has emerged as a promising technology for removing MCs. Many MCs-degrading bacteria have been identified, including an indigenous bacterium Sphingopyxis sp. YF1 that could degrade MC-LR and Adda completely. Herein, we gained insight into the MCs biodegradation mechanisms and evolutionary dynamics of MCs-degrading bacteria, and revealed the toxic risks of the MCs degradation products. The biochemical characteristics and genetic repertoires of strain YF1 were explored. A comparative genomic analysis was performed on strain YF1 and six other MCs-degrading bacteria to investigate their functions. The degradation products were investigated, and the toxicity of the intermediates was analyzed through rigorous theoretical calculation. Strain YF1 might be a novel species that exhibited versatile substrate utilization capabilities. Many common genes and metabolic pathways were identified, shedding light on shared functions and catabolism in the MCs-degrading bacteria. The crucial genes involved in MCs catabolism mechanisms, including mlr and paa gene clusters, were identified successfully. These functional genes might experience horizontal gene transfer events, suggesting the evolutionary dynamics of these MCs-degrading bacteria in ecology. Moreover, the degradation products for MCs and Adda were summarized, and we found most of the intermediates exhibited lower toxicity to different organisms than the parent compound. These findings systematically revealed the MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria. Consequently, this research contributed to the advancement of green biodegradation technology in aquatic ecology, which might protect human health from MCs.

Metabolomic and transcriptomic reveal flavonoid biosynthesis and regulation mechanism in Phlomoides rotata from different habitats

Phlomoides rotata is a traditional medical plant at 3100-5200 m altitude in the Tibet Plateau. In this study, flavonoid metabolites were investigated in P. rotata from Henan County (HN), Guoluo County (GL), Yushu County (YS), and Chengduo County (CD) habitats in Qinghai. The level of kaempferol 3-neohesperidoside, sakuranetin, and biochanin A was high in HN. The content of limocitrin and isoquercetin was high in YS. The levels of ikarisoside A and chrysosplenol D in GL were high. Schaftoside, miquelianin, malvidin chloride, and glabrene in CD exhibited high levels. The results showed a significant correlation between 59 flavonoids and 29 DEGs. Eleven flavonoids increased with altitude. PAL2, UFGT6, COMT1, HCT2, 4CL4, and HCT3 genes were crucial in regulating flavonoid biosynthesis. Three enzymes CHS, 4CL, and UFGT, were crucial in regulating flavonoid biosynthesis. This study provided biological and chemical evidence for the different uses of various regional plants of P. rotata.

Screening, identification and evaluation of an acidophilic strain of Bacillus velezensis B4-7 for the biocontrol of tobacco bacterial wilt

Tobacco (Nicotiana tabacum L.) bacterial wilt, caused by Ralstonia solanacearum, is indeed a highly destructive plant disease, leading to substantial damage in tobacco production. While biological control is considered an effective measure for managing bacterial wilt, related research in this area has been relatively limited compared to other control methods. In order to discover new potential antagonistic bacteria with high biocontrol efficacy against tobacco bacterial wilt, we conducted an analysis of the microbial composition differences between disease-suppressive and disease-conducive soils using Illumina sequencing. As a result, we successfully isolated six strains from the disease-suppressive soil that exhibited antibacterial activity against Ralstonia solanacearum. Among these strains, B4-7 showed the strongest antibacterial activity, even at acidic conditions with a pH of 4.0. Based on genome analysis using Average Nucleotide Identity (ANI), B4-7 was identified as Bacillus velezensis. In greenhouse and field trials, strain B4-7 significantly reduced the disease index of tobacco bacterial wilt, with control efficiencies reaching 74.03% and 46.88% respectively. Additionally, B4-7 exhibited plant-promoting abilities that led to a 35.27% increase in tobacco production in field conditions. Quantitative real-time (qPCR) analysis demonstrated that strain B4-7 effectively reduced the abundance of R. solanacearum in the rhizosphere. Genome sequencing and Liquid Chromatography-Mass Spectrometry (LC-MS) analysis revealed that strain B4-7 potentially produces various lipopeptide metabolites, such as microlactin, bacillaene, difficidin, bacilysin, and surfactin. Furthermore, B4-7 influenced the structure of the rhizosphere soil microbial community, increasing bacterial abundance and fungal diversity, while also promoting the growth of different beneficial microorganisms. In addition, B4-7 enhanced tobacco's resistance to R. solanacearum by increasing the activities of defense enzymes, including superoxide dismutase (SOD), phenylalanine ammonia-lyase (PAL), peroxidase (POD), catalase (CAT), and polyphenol oxidase (PPO). Collectively, these findings suggest that B. velezensis B4-7 holds significant biocontrol potential and can be considered a promising candidate strain for eco-friendly management of tobacco bacterial wilt.

A comprehensive analysis of SARS-CoV-2 missense mutations indicates that all possible amino acid replacements in the viral proteins occurred within the first two-and-a-half years of the pandemic

The surveillance of COVID-19 pandemic has led to the determination of millions of genome sequences of the SARS-CoV-2 virus, with the accumulation of a wealth of information never collected before for an infectious disease. Exploring the information retrieved from the GISAID database reporting at that time >13 million genome sequences, we classified the 141,639 unique missense mutations detected in the first two-and-a-half years (up to October 2022) of the pandemic. Notably, our analysis indicates that 98.2 % of all possible conservative amino acid replacements occurred. Even non-conservative mutations were highly represented (73.9 %). For a significant number of residues (3 %), all possible replacements with the other nineteen amino acids have been observed. These observations strongly indicate that, in this time interval, the virus explored all possible alternatives in terms of missense mutations for all sites of its polypeptide chain and that those that are not observed severely affect SARS-CoV-2 integrity. The implications of the present findings go well beyond the structural biology of SARS-CoV-2 as the huge amount of information here collected and classified may be valuable for the elucidation of the sequence-structure-function relationships in proteins.

Transcriptome Analysis of mfs2-Defective Penicillium digitatum Mutant to Reveal Importance of Pdmfs2 in Developing Fungal Prochloraz Resistance

Demethylation inhibitors (DMIs), including prochloraz, are popular fungicides to control citrus postharvest pathogens such as Penicillium digitatum (green mold). However, many P. digitatum strains have developed prochloraz resistance, which decreases drug efficacy. Specific major facilitator superfamily (MFS) transporter gene mfs2, encoding drug-efflux pump protein MFS2, has been identified in P. digitatum strain F6 (PdF6) to confer fungal strain prochloraz resistance. However, except for the drug-efflux pump function of MFS2, other mechanisms relating to the Pdmfs2 are not fully clear. The present study reported a transcriptome investigation on the mfs2-defective P. digitatum strain. Comparing to the wild-type strain, the mfs2-defective strain showed 717 differentially expressed genes (DEGs) without prochloraz induction, and 1221 DEGs with prochloraz induction. The obtained DEGs included multiple isoforms of MFS transporter-encoding genes, ATP-binding cassette (ABC) transporter-encoding genes, and multidrug and toxic compound extrusion (MATE) family protein-encoding genes. Many of these putative drug-efflux pump protein-encoding genes had significantly lower transcript abundances in the mfs2-defective P. digitatum strain at prochloraz induction, as compared to the wild-type strain, including twenty-two MFS transporter-encoding genes (MFS1 to MFS22), two ABC transporter-encoding genes (ABC1 and ABC2), and three MATE protein-encoding genes (MATE1 to MATE3). The prochloraz induction on special drug-efflux pump protein genes in the wild-type strain was not observed in the mfs2-defective strain, including MFS21, MFS22, ABC2, MATE1, MATE2, and MATE3. On the other hand, the up-regulation of other drug-efflux pump protein genes in the mfs2-defective strain cannot recover the fungal prochloraz resistance, including MFS23, MFS26, MFS27, MFS31, MFS33, and ABC3 to ABC8. The functional enrichment of DEGs based on Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and euKaryotic Orthologous Groups (KOG) database resources suggested some essential contributors to the mfs2-relating prochloraz resistance, including ribosome biosynthesis-related genes, oxidative phosphorylation genes, steroid biosynthesis-related genes, fatty acid and lipid metabolism-related genes, and carbon- and nitrogen-metabolism-related genes. The results indicated that the MFS2 transporter might be involved in the regulation of multiple drug-efflux pump protein gene expressions and multiple metabolism-related gene expressions, thus playing an important role in developing P. digitatum prochloraz resistance.

WGCNA Reveals Hub Genes and Key Gene Regulatory Pathways of the Response of Soybean to Infection by Soybean mosaic virus

Soybean mosaic virus (SMV) is one of the main pathogens that can negatively affect soybean production and quality. To study the gene regulatory network of soybeans in response to SMV SC15, the resistant line X149 and susceptible line X97 were subjected to transcriptome analysis at 0, 2, 8, 12, 24, and 48 h post-inoculation (hpi). Differential expression analysis revealed that 10,190 differentially expressed genes (DEGs) responded to SC15 infection. Weighted gene co-expression network analysis (WGCNA) was performed to identify highly related resistance gene modules; in total, eight modules, including 2256 DEGs, were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of 2256 DEGs revealed that the genes significantly clustered into resistance-related pathways, such as the plant-pathogen interaction pathway, mitogen-activated protein kinases (MAPK) signaling pathway, and plant hormone signal transduction pathway. Among these pathways, we found that the flg22, Ca2+, hydrogen peroxide (H2O2), and abscisic acid (ABA) regulatory pathways were fully covered by 36 DEGs. Among the 36 DEGs, the gene Glyma.01G225100 (protein phosphatase 2C, PP2C) in the ABA regulatory pathway, the gene Glyma.16G031900 (WRKY transcription factor 22, WRKY22) in Ca2+ and H2O2 regulatory pathways, and the gene Glyma.04G175300 (calcium-dependent protein kinase, CDPK) in Ca2+ regulatory pathways were highly connected hub genes. These results indicate that the resistance of X149 to SC15 may depend on the positive regulation of flg22, Ca2+, H2O2, and ABA regulatory pathways. Our study further showed that superoxide dismutase (SOD) activity, H2O2 content, and catalase (CAT) and peroxidase (POD) activities were significantly up-regulated in the resistant line X149 compared with those in 0 hpi. This finding indicates that the H2O2 regulatory pathway might be dependent on flg22- and Ca2+-pathway-induced ROS generation. In addition, two hub genes, Glyma.07G190100 (encoding F-box protein) and Glyma.12G185400 (encoding calmodulin-like proteins, CMLs), were also identified and they could positively regulate X149 resistance. This study provides pathways for further investigation of SMV resistance mechanisms in soybean.

Deciphering the microbial communities of alkaline hot spring in Panamik, Ladakh, India using a high-throughput sequencing approach

Due to their distinctive physicochemical characteristics, hot springs are extremely important. The whole genome metagenomic sequencing technology can be utilized to analyze the diverse microbial community that thrives in this habitat due to the particular selection pressure that prevails there. The current investigation emphasizes on culture-independent metagenomic study of the Panamik hot spring and its nearby areas from Ladakh, India. Based on different diversity indices, sequence analysis of the soil reservoir showed higher species richness and diversity in comparison to water and sediment samples. The mineral content and various physicochemical pameters like temperature, pH had an impact on the composition of the microbial community of the geothermal springs. The phyla Proteobacteria, Cyanobacteria, Bacteroidetes, Actinobacter, Firmicutes, and Verrucomicrobia in bacterial domain dominate the thermos-alkaline spring at Panamik in different concentrations. Economically significant microbes from the genera Actinobacter, Thermosynechoccus, Candidatus Solibacter, Chthoniobacter, Synechoccus, Pseudomonas and Sphingomonas, were prevalent in hot spring. In the archaeal domain, the most dominant phylum and genera were Euryarchaeota and Thermococcus in all the samples. Further, the most abundant species were Methanosarcina barkeri, Nitrospumilus maritimus and Methanosarcina acetivorans. The present study which only examined one of the several thermal springs present in the Himalayan geothermal area, should be regarded as a preliminary investigation of the microbiota that live in the hot springs on these remote areas. These findings suggest that further investigations should be undertaken to characterize the ecosystems of the Panamik hot spring, which serve as a repository for unidentified microbial lineages.

Full-length transcriptome analysis of a bloom-forming dinoflagellate Prorocentrum shikokuense (Dinophyceae)

Prorocentrum shikokuense (formerly P. donghaiense) is a pivotal dinoflagellate species associating with the HABs in the East China Sea. The complexity of its large nuclear genome hindered us from understanding its genomic characteristics. Full-length transcriptome sequencing offers a practical solution to decipher the physiological mechanisms of a species without the reference genome. In this study, we employed single-molecule real-time (SMRT) sequencing technology to sequence the full-length transcriptome of Prorocentrum shikokuense. We successfully generated 41.73 Gb of clean SMRT sequencing reads and isolated 105,249 non-redundant full-length non-chimeric reads. Our trial has led to the identification of 11,917 long non-coding RNA transcripts, 514 alternative splicing events, 437 putative transcription factor genes from 17 TF gene families, and 34,723 simple sequence repeats. Additionally, a total of 78,265 open reading frames were identified, of them 15,501 were the protein coding sequences. This dataset is valuable for annotating P. shikokuense genome, and will contribute significantly to the in-depth studies on the molecular mechanisms underlining the dinoflagellate bloom formation.

Genomic characterization and probiotic potential assessment of an exopolysaccharide-producing strain Pediococcus pentosaceus LL-07 isolated from fermented meat

BACKGROUND

The genomic information available for Pediococcus pentosaceus is primarily derived from fermented fruits and vegetables, with less information available from fermented meat. P. pentosaceus LL-07, a strain isolated from fermented meat, has the capability of producing exopolysaccharides (EPS). To assess the probiotic attributes of P. pentosaceus LL-07, we conducted whole-genome sequencing (WGS) using the PacBio SequelIIe and Illumina MiSeq platforms, followed by in vitro experiments to explore its probiotic potential.

RESULTS

The genome size of P. pentosaceus LL-07 is 1,782,685 bp, comprising a circular chromosome and a circular plasmid. Our investigation revealed the absence of a CRISPR/Cas system. Sugar fermentation experiments demonstrated the characteristics of carbohydrate metabolism. P. pentosaceus LL-07 contains an EPS synthesis gene cluster consisting of 13 genes, which is different from the currently known gene cluster structure. NO genes associated with hemolysis or toxin synthesis were detected. Additionally, eighty-six genes related to antibiotic resistance were identified but not present in the prophage, transposon or plasmid. In vitro experiments demonstrated that P. pentosaceus LL-07 was comparable to the reference strain P. pentosaceus ATCC25745 in terms of tolerance to artificial digestive juice and bile, autoaggregation and antioxidation, and provided corresponding genomic evidence.

CONCLUSION

This study confirmed the safety and probiotic properties of P. pentosaceus LL-07 via complete genome and phenotype analysis, supporting its characterization as a potential probiotic candidate.

Genomic characterization and related functional genes of γ- poly glutamic acid producing Bacillus subtilis

γ- poly glutamic acid (γ-PGA), a high molecular weight polymer, is synthesized by microorganisms and secreted into the extracellular space. Due to its excellent performance, γ-PGA has been widely used in various fields, including food, biomedical and environmental fields. In this study, we screened natto samples for two strains of Bacillus subtilis N3378-2at and N3378-3At that produce γ-PGA. We then identified the γ-PGA synthetase gene cluster (PgsB, PgsC, PgsA, YwtC and PgdS), glutamate racemase RacE, phage-derived γ-PGA hydrolase (PghB and PghC) and exo-γ-glutamyl peptidase (GGT) from the genome of these strains. Based on these γ-PGA-related protein sequences from isolated Bacillus subtilis and 181 B. subtilis obtained from GenBank, we carried out genotyping analysis and classified them into types 1-5. Since we found B. amyloliquefaciens LL3 can produce γ-PGA, we obtained the B. velezensis and B. amyloliquefaciens strains from GenBank and classified them into types 6 and 7 based on LL3. Finally, we constructed evolutionary trees for these protein sequences. This study analyzed the distribution of γ-PGA-related protein sequences in the genomes of B. subtilis, B. velezensis and B. amyloliquefaciens strains, then the evolutionary diversity of these protein sequences was analyzed, which provided novel information for the development and utilization of γ-PGA-producing strains.

Identification of positively selected genes in Mycobacterium tuberculosis from southern Xinjiang Uygur autonomous region of China

Background

Tuberculosis (TB), mainly caused by Mycobacterium tuberculosis (Mtb), remains a serious public health problem. Increasing evidence supports that selective evolution is an important force affecting genomic determinants of Mtb phenotypes. It is necessary to further understand the Mtb selective evolution and identify the positively selected genes that probably drive the phenotype of Mtb.

Methods

This study mainly focused on the positive selection of 807 Mtb strains from Southern Xinjiang of China using whole genome sequencing (WGS). PAML software was used for identifying the genes and sites under positive selection in 807 Mtb strains.

Results

Lineage 2 (62.70%) strains were the dominant strains in this area, followed by lineage 3 (19.45%) and lineage 4 (17.84%) strains. There were 239 codons in 47 genes under positive selection, and the genes were majorly associated with the functions of transcription, defense mechanisms, and cell wall/membrane/envelope biogenesis. There were 28 codons (43 mutations) in eight genes (gyrA, rpoB, rpoC, katG, pncA, embB, gid, and cut1) under positive selection in multi-drug resistance (MDR) strains but not in drug-susceptible (DS) strains, in which 27 mutations were drug-resistant loci, 9 mutations were non-drug-resistant loci but were in drug-resistant genes, 2 mutations were compensatory mutations, and 5 mutations were in unknown drug-resistant gene of cut1. There was a codon in Rv0336 under positive selection in L3 strains but not in L2 and L4 strains. The epitopes of T and B cells were both hyper-conserved, particularly in the T-cell epitopes.

Conclusion

This study revealed the ongoing selective evolution of Mtb. We found some special genes and sites under positive selection which may contribute to the advantage of MDR and L3 strains. It is necessary to further study these mutations to understand their impact on phenotypes for providing more useful information to develop new TB interventions.

A chromosomal-scale genome assembly of modern cultivated hybrid sugarcane provides insights into origination and evolution

Sugarcane is a vital crop with significant economic and industrial value. However, the cultivated sugarcane's ultra-complex genome still needs to be resolved due to its high ploidy and extensive recombination between the two subgenomes. Here, we generate a chromosomal-scale, haplotype-resolved genome assembly for a hybrid sugarcane cultivar ZZ1. This assembly contains 10.4 Gb genomic sequences and 68,509 annotated genes with defined alleles in two sub-genomes distributed in 99 original and 15 recombined chromosomes. RNA-seq data analysis shows that sugar accumulation-associated gene families have been primarily expanded from the ZZSO subgenome. However, genes responding to pokkah boeng disease susceptibility have been derived dominantly from the ZZSS subgenome. The region harboring the possible smut resistance genes has expanded significantly. Among them, the expansion of WAK and FLS2 families is proposed to have occurred during the breeding of ZZ1. Our findings provide insights into the complex genome of hybrid sugarcane cultivars and pave the way for future genomics and molecular breeding studies in sugarcane.

Inter-species Transcriptomic Analysis Reveals a Constitutive Adaptation Against Oxidative Stress for the Highly Virulent Leptospira Species

Transcriptomic analyses across large scales of evolutionary distance have great potential to shed light on regulatory evolution but are complicated by difficulties in establishing orthology and limited availability of accessible software. We introduce here a method and a graphical user interface wrapper, called Annotator-RNAtor, for performing interspecies transcriptomic analysis and studying intragenus evolution. The pipeline uses third-party software to infer homologous genes in various species and highlight differences in the expression of the core-genes. To illustrate the methodology and demonstrate its usefulness, we focus on the emergence of the highly virulent Leptospira subclade known as P1+, which includes the causative agents of leptospirosis. Here, we expand on the genomic study through the comparison of transcriptomes between species from P1+ and their related P1- counterparts (low-virulent pathogens). In doing so, we shed light on differentially expressed pathways and focused on describing a specific example of adaptation based on a differential expression of PerRA-controlled genes. We showed that P1+ species exhibit higher expression of the katE gene, a well-known virulence determinant in pathogenic Leptospira species correlated with greater tolerance to peroxide. Switching PerRA alleles between P1+ and P1- species demonstrated that the lower repression of katE and greater tolerance to peroxide in P1+ species was solely controlled by PerRA and partly caused by a PerRA amino-acid permutation. Overall, these results demonstrate the strategic fit of the methodology and its ability to decipher adaptive transcriptomic changes, not observable by comparative genome analysis, that may have been implicated in the emergence of these pathogens.

Whole transcriptome sequencing analyses of islets reveal ncRNA regulatory networks underlying impaired insulin secretion and increased β-cell mass in high fat diet-induced diabetes mellitus

AIM

Our study aims to identify novel non-coding RNA-mRNA regulatory networks associated with β-cell dysfunction and compensatory responses in obesity-related diabetes.

METHODS

Glucose metabolism, islet architecture and secretion, and insulin sensitivity were characterized in C57BL/6J mice fed on a 60% high-fat diet (HFD) or control for 24 weeks. Islets were isolated for whole transcriptome sequencing to identify differentially expressed (DE) mRNAs, miRNAs, IncRNAs, and circRNAs. Regulatory networks involving miRNA-mRNA, lncRNA-mRNA, and lncRNA-miRNA-mRNA were constructed and functions were assessed through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses.

RESULTS

Despite compensatory hyperinsulinemia and a significant increase in β-cell mass with a slow rate of proliferation, HFD mice exhibited impaired glucose tolerance. In isolated islets, insulin secretion in response to glucose and palmitic acid deteriorated after 24 weeks of HFD. Whole transcriptomic sequencing identified a total of 1324 DE mRNAs, 14 DE miRNAs, 179 DE lncRNAs, and 680 DE circRNAs. Our transcriptomic dataset unveiled several core regulatory axes involved in the impaired insulin secretion in HFD mice, such as miR-6948-5p/Cacna1c, miR-6964-3p/Cacna1b, miR-3572-5p/Hk2, miR-3572-5p/Cckar and miR-677-5p/Camk2d. Additionally, proliferative and apoptotic targets, including miR-216a-3p/FKBP5, miR-670-3p/Foxo3, miR-677-5p/RIPK1, miR-802-3p/Smad2 and ENSMUST00000176781/Caspase9 possibly contribute to the increased β-cell mass in HFD islets. Furthermore, competing endogenous RNAs (ceRNA) regulatory network involving 7 DE miRNAs, 15 DE lncRNAs and 38 DE mRNAs might also participate in the development of HFD-induced diabetes.

CONCLUSIONS

The comprehensive whole transcriptomic sequencing revealed novel non-coding RNA-mRNA regulatory networks associated with impaired insulin secretion and increased β-cell mass in obesity-related diabetes.

Expanded phylogeny of extremely halophilic archaea shows multiple independent adaptations to hypersaline environments

Extremely halophilic archaea (Haloarchaea, Nanohaloarchaeota, Methanonatronarchaeia and Halarchaeoplasmatales) thrive in saturating salt concentrations where they must maintain osmotic equilibrium with their environment. The evolutionary history of adaptations enabling salt tolerance remains poorly understood, in particular because the phylogeny of several lineages is conflicting. Here we present a resolved phylogeny of extremely halophilic archaea obtained using improved taxon sampling and state-of-the-art phylogenetic approaches designed to cope with the strong compositional biases of their proteomes. We describe two uncultured lineages, Afararchaeaceae and Asbonarchaeaceae, which break the long branches at the base of Haloarchaea and Nanohaloarchaeota, respectively. We obtained 13 metagenome-assembled genomes (MAGs) of these archaea from metagenomes of hypersaline aquatic systems of the Danakil Depression (Ethiopia). Our phylogenomic analyses including these taxa show that at least four independent adaptations to extreme halophily occurred during archaeal evolution. Gene-tree/species-tree reconciliation suggests that gene duplication and horizontal gene transfer played an important role in this process, for example, by spreading key genes (such as those encoding potassium transporters) across extremely halophilic lineages.

Transcriptomic analysis and oxidative stress induced by sodium dichloroisocyanurate in the intestine of Phascolosoma esculenta

Sodium dichloroisocyanurate (NaDCC, C3Cl2N3NaO3) is a solid chlorine-containing product that is widely used as a disinfectant in living environments, which has potential toxic effects on human and rats. Phascolosoma esculenta is a species native to the southeast coast of China and can be used as an indicator organism. In the present study, 150 P. esculenta were used to determine the LC50 of NaDCC for P. esculenta, then 100 P. esculenta were used to analysis the change of histopathology, oxidative stress and transcriptome after NaDCC exposure. The results showed that the LC50 of NaDCC for 48 h was 50 mg/L. NaDCC stress induced pathological events in P. esculenta, including blisters, intestinal structural damage and epithelial cell ruptured or even loss. The highest and lowest intestinal activity of superoxide dismutase in individual survivors was detected at 12 h and 72 h, respectively. Malondialdehyde levels in the intestine declined gradually from 3 h and increased at 9 h, and peaked at 12 h. Total antioxidant capacity declined at 3 h and dropped below the levels of control group after 9 h. Transcriptome sequencing analysis yielded a total of 48.65 Gb of clean data. A total of 34,759 new genes were found including 957 differentially expressed genes (DEGs). The DEGs were significantly enriched in ferroptosis, response to chemicals, response to stress, immune system, ion transport, cell death, oxidation-reduction, cellular homeostasis, protein ubiquitination, and protein neddylation. Additionally, the levels of detoxification enzymes, such as glutathione-S-transferase, cytochrome P450, ABC, UDP-glycosyltransferase and SLC transporters of endogenous and exogenous solutes were significantly changed. Overall, the results provide reference for reasonable use of disinfectants during farming, and also provide insight into the mechanisms related to NaDCC toxicity in P. esculenta.

Transcriptomic analysis provides an insight into the function of CmGH9B3, a key gene of β-1, 4-glucanase, during the graft union healing of oriental melon scion grafted onto squash rootstock

The melon (Cucumis melo L.) is a globally cherished and economically significant crop. The grafting technique has been widely used in the vegetative propagation of melon to promote environmental tolerance and disease resistance. However, mechanisms governing graft healing and potential incompatibilities in melons following the grafting process remain unknown. To uncover the molecular mechanism of healing of grafted melon seedlings, melon wild type (Control) and TRV-CmGH9B3 lines were obtained and grafted onto the squash rootstocks (C. moschata). Anatomical differences indicated that the healing process of the TRV-CmGH9B3 plants was slower than that of the control. A total of 335 significantly differentially expressed genes (DEGs) were detected between two transcriptomes. Most of these DEGs were down-regulated in TRV-CmGH9B3 grafted seedlings. GO and KEGG analysis showed that many metabolic, physiological, and hormonal responses were involved in graft healing, including metabolic processes, plant hormone signaling, plant MAPK pathway, and sucrose starch pathway. During the healing process of TRV-CmGH9B3 grafted seedlings, gene synthesis related to hormone signal transduction (auxin, cytokinin, gibberellin, brassinolide) was delayed. At the same time, it was found that most of the DEGs related to the sucrose pathway were down-regulated in TRV-CmGH9B3 grafted seedlings. The results showed that sugar was also involved in the healing process of melon grafted onto squash. These results deepened our understanding of the molecular mechanism of GH9B3, a key gene of β-1, 4-glucanase. It also provided a reference for elucidating the gene mechanism and function analysis of CmGH9B3 in the process of graft union healing.

Similarities and differences: species and diet impact gut microbiota of captive pheasants

The fecal microbiota plays an important role in maintaining animal health and is closely related to host life activities. In recent years, there have been an increasing number of studies on the fecal microbiota from birds. An exploration of the effects of species and living environments on the composition of gut microbiota will provide better protection for wildlife. In this study, non-injury sampling and 16S rDNA high-throughput sequencing were used to investigate the bacterial composition and diversity of the fecal microbiota in silver pheasants (Lophura nycthemera) and golden pheasants (Chrysolophus pictus) from Tianjin Zoo and Beijing Wildlife Park. The results showed that the abundance of Firmicutes was the highest in all fecal samples. At the genus level, Bacteroides was the common dominant bacteria, while there were some differences in other dominant bacteria genera. There were significant differences in fecal microbial composition between the golden pheasants from Tianjin Zoo and Beijing Wildlife Park. The metabolic analysis and functional prediction suggested that the gut microbiota composition and host metabolism were influenced by dietary interventions and living conditions. The results of this study provide the basis for further research of intestinal microbial of L. nycthemera and C. pictus, and valuable insights for conservation of related species.

Metagenomics untangles potential adaptations of Antarctic endolithic bacteria at the fringe of habitability

Survival and growth strategies of Antarctic endolithic microbes residing in Earth's driest and coldest desert remain virtually unknown. From 109 endolithic microbiomes, 4539 metagenome-assembled genomes were generated, 49.3 % of which were novel candidate bacterial species. We present evidence that trace gas oxidation and atmospheric chemosynthesis may be the prevalent strategies supporting metabolic activity and persistence of these ecosystems at the fringe of life and the limits of habitability.

The Home Environment Is a Reservoir for Methicillin-Resistant Coagulase-Negative Staphylococci and Mammaliicocci

Coagulase-negative staphylococci (CoNS) and mammaliicocci are opportunistic human and animal pathogens, often resistant to multiple antimicrobials, including methicillin. Methicillin-resistant CoNS (MRCoNS) have traditionally been linked to hospitals and healthcare facilities, where they are significant contributors to nosocomial infections. However, screenings of non-hospital environments have linked MRCoNS and methicillin-resistant mammaliicocci (MRM) to other ecological niches. The aim of this study was to explore the home environment as a reservoir for MRCoNS and MRM. A total of 33 households, including households with a dog with a methicillin-resistant staphylococcal infection, households with healthy dogs or cats and households without pets, were screened for MRCoNS and MRM by sampling one human, one pet (if present) and the environment. Samples were analyzed by a selective culture-based method, and bacterial species were identified by MALDI-TOF MS and tested for antibiotic susceptibility by the agar disk diffusion method. Following whole-genome sequencing, a large diversity of SCCmec elements and sequence types was revealed, which did not indicate any clonal dissemination of specific strains. Virulome and mobilome analyses indicated a high degree of species specificity. Altogether, this study documents that the home environment is a reservoir for a variety of MRCoNS and MRM regardless of the type of household.

Towards modelling tick-virus interactions using the weakly pathogenic Sindbis virus: Evidence that ticks are competent vectors

Most tick-borne viruses (TBVs) are highly pathogenic and require high biosecurity, which severely limits their study. We found that Sindbis virus (SINV), predominantly transmitted by mosquitoes, can replicate in ticks and be subsequently transmitted, with the potential to serve as a model for studying tick-virus interactions. We found that both larval and nymphal stages of Rhipicephalus haemaphysaloides can be infected with SINV-wild-type (WT) when feeding on infected mice. SINV replicated in two species of ticks (R. haemaphysaloides and Hyalomma asiaticum) after infecting them by microinjection. Injection of ticks with SINV expressing enhanced Green Fluorescent Protein (eGFP) revealed that SINV-eGFP specifically aggregated in the tick midguts for replication. During blood-feeding, SINV-eGFP migrated from the midguts to the salivary glands and was transmitted to a new host. SINV infection caused changes in expression levels of tick genes related to immune responses, substance transport and metabolism, cell growth and death. SINV mainly induced autophagy during the early stage of infection; with increasing time of infection, the level of autophagy decreased, while the level of apoptosis increased. During the early stages of infection, the transcript levels of immune-related genes were significantly upregulated, and then decreased. In addition, SINV induced changes in the transcription levels of some functional genes that play important roles in the interactions between ticks and tick-borne pathogens. These results confirm that the SINV-based transmission model between ticks, viruses, and mammals can be widely used to unravel the interactions between ticks and viruses.

Annotating microbial functions with ProkFunFind

Analyzing microbial genomes has become an essential part of microbiology research, giving valuable insights into the functions and evolution of microbial species. Identifying genes of interest and assigning putative annotations to those genes is a central task in genome analysis, and a plethora of tools and approaches have been developed for this task. The ProkFunFind tool was developed to bridge the gap between these various annotation approaches, providing a flexible and customizable search approach to annotate microbial functions. ProkFunFind is designed around hierarchical definitions of biological functions, where individual genes can be identified using heterogeneous search terms consisting of sequences, profile hidden Markov models, protein domains, and orthology groups. This flexible and customizable search approach allows for searches to be tailored to specific biological functions, and the search results are output in multiple formats to facilitate downstream analyses. The utility of the ProkFunFind search tool was demonstrated through its application in searching for bacterial flagella, which are complex organelles composed of multiple genes. Overall, ProkFunFind provides an accessible and flexible way to integrate multiple types of annotation and sequence data while annotating biological functions in microbial genomes.IMPORTANCEGenome sequencing and analysis are increasingly important parts of microbiology, providing a way to predict metabolic functions, identify virulence factors, and understand the evolution of microbes. The expanded use of genome sequencing has also brought an abundance of search and annotation methods, but integrating the information from these different methods can be challenging and is often done through ad hoc approaches. To bridge the gap between different types of annotations, we developed ProkFunFind, a flexible and customizable search tool incorporating multiple search approaches and annotation types to annotate microbial functions. We demonstrated the utility of ProkFunFind by searching for gene clusters encoding flagellar genes using a combination of different annotation types and searches. Overall, ProkFunFind provides a reproducible and flexible way to identify gene clusters of interest, facilitating the meaningful analysis of new and existing microbial genomes.

Cyberlindnera jadinii yeast as a functional protein source: Modulation of immunoregulatory pathways in the intestinal proteome of zebrafish (Danio rerio)

Yeasts contain bioactive components that can enhance fish immune robustness and disease resistance. Our study focused on analyzing intestinal immunoregulatory pathways in zebrafish (Danio rerio) using iTRAQ and 2D LC-MS/MS to quantify intestinal proteins. Zebrafish were fed either control diet (C) or diet supplemented with autolyzed Cyberlindnera jadinii (ACJ). KEGG analysis revealed that ACJ yeast diet induced increased abundance of proteins related to arginine and proline metabolism, phagosome, C-lectin receptor signaling, ribosome and PPAR signaling pathways, which can modulate and enhance innate immune responses. ACJ yeast diet also showed decreased abundance of proteins associated with inflammatory pathways, including apoptosis, necroptosis and ferroptosis. These findings indicate boosted innate immune response and control of inflammation-related pathways in zebrafish intestine. Our findings in the well annotated proteome of zebrafish enabled a detailed investigation of intestinal responses and provide insight into health-beneficial effects of yeast species C. jadinii, which is relevant for aquaculture species.

Transcriptomic profiling of Poa pratensis L. under treatment of various phytohormones

Poa pratensis L. (Poaceae) is a valuable grass across the north hemisphere, inhabiting diverse environments with wide altitudinal span, where ubiquitous various kinds of stresses. Phytohormones would be helpful to improve tolerance to abiotic and biotic stresses, but the responses of transcriptome regulation of P. pratensis to exogenous phytohormones application remain unclear. In this study, we explored the alteration of plant physiological responses by the application of phytohormones. Aiming to achieve this knowledge, we got full-length transcriptome data 42.76 Gb, of which 74.9% of transcripts were completed. Then used 27 samples representing four treatments conducted at two time points (1 h and 6 h after application) to generate RNA-seq data. 371 and 907 common DEGs were identified in response to four phytohormones application, respectively, these DEGs were involved in "plant hormone signal transduction", "carbon metabolism" and "plant-pathogen interaction". Finally, P. pratensis basic research can gain valuable information regarding the responses to exogenous application of phytohormones in physiological indicators and transcriptional regulations in order to facilitate the development of new cultivars.

Genome Sequence Analysis of Calcifying Bacteria Bacillus paranthracis CT5 and Its Biomineralization Efficacy to Improve the Strength and Durability Properties of Civil Structures

Bacteria producing urea amidohydrolases (UA) and carbonic anhydrases (CA) are of great importance in civil engineering as these enzymes are responsible for microbially induced calcium carbonate precipitation (MICCP). In this investigation, genomic insights of Bacillus paranthracis CT5 and the expression of genes underlying in MICCP were studied. B. paranthracis produced a maximum level of UA (669.3 U/ml) and CA (125 U/ml) on 5th day of incubation and precipitated 197 mg/100 ml CaCO3 after 7 days of incubation. After 28 days of curing, compressive strength of bacterial admixed and bacterial cured (B-B) specimens was 13.7% higher compared to water-mixed and water-cured (W-W) specimens. A significant decrease in water absorption was observed in bacterial-cured specimens compared to water-cured specimens after 28 days of curing. For genome analysis, reads were assembled de novo producing 5,402,771 bp assembly with N50 of 273,050 bp. RAST annotation detected six amidohydrolase and three carbonic anhydrase genes. Among 5700 coding sequences found in genome, COG gene annotation grouped 4360 genes into COG categories with highest number of genes to transcription (435 genes), amino acid transport and metabolism (362 genes) along with cell wall/membrane/envelope biogenesis and ion transport and metabolism. KEGG functional classification predicted 223 pathways consisting of 1,960 genes and the highest number of genes belongs to two-component system (101 genes) and ABC transporter pathways (98 genes) enabling bacteria to sense and respond to environmental signals and actively transport various minerals and organic molecules, which facilitate the active transport of molecules required for MICCP.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. This is in part due to their large genome sizes and difficulty in obtaining sequenceable material. In this manuscript, we leveraged a recently-isolated, laboratory-tractable Entomophthora muscae isolate and improved long-read sequencing to obtain a largely-complete entomophthoralean genome. Our E. muscae assembly is 1.03 Gb, consists of 7,810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with other available (transcriptomic and genomic) datasets from entomophthoralean fungi, we provide new insight into the biology of these understudied pathogens. We offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex. Our findings suggest that substantial taxonomic revision is needed to define species within this group and we provide recommendations for differentiating strains and species in the context of the existing body of E. muscae scientific literature. We show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity, despite the presence of machinery to defend against transposable elements(RNAi). In addition, we find that E. muscae and its closest allies are enriched for M16A peptidases and possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We find that E. muscae has an expanded group of acid-trehalases, consistent with trehalose being the primary sugar component of fly (and insect) hemolymph. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

A metagenomics pipeline reveals insertion sequence-driven evolution of the microbiota

Insertion sequence (IS) elements are mobile genetic elements in bacterial genomes that support adaptation. We developed a database of IS elements coupled to a computational pipeline that identifies IS element insertions in the microbiota. We discovered that diverse IS elements insert into the genomes of intestinal bacteria regardless of human host lifestyle. These insertions target bacterial accessory genes that aid in their adaptation to unique environmental conditions. Using IS expansion in Bacteroides, we show that IS activity leads to insertion "hot spots" in accessory genes. We show that IS insertions are stable and can be transferred between humans. Extreme environmental perturbations force IS elements to fall out of the microbiota and many fail to rebound following homeostasis. Our work shows that IS elements drive bacterial genome diversification within the microbiota and establishes a framework for understanding how strain level variation within the microbiota impacts human health.

Characterization of PYL gene family and identification of HaPYL genes response to drought and salt stress in sunflower

In the context of global climate change, drought and soil salinity are some of the most devastating abiotic stresses affecting agriculture today. PYL proteins are essential components of abscisic acid (ABA) signaling and play critical roles in responding to abiotic stressors, including drought and salt stress. Although PYL genes have been studied in many species, their roles in responding to abiotic stress are still unclear in the sunflower. In this study, 19 HaPYL genes, distributed on 15 of 17 chromosomes, were identified in the sunflower. Fragment duplication is the main cause of the expansion of PYL genes in the sunflower genome. Based on phylogenetic analysis, HaPYL genes were divided into three subfamilies. Members in the same subfamily share similar protein motifs and gene exon-intron structures, except for the second subfamily. Tissue expression patterns suggested that HaPYLs serve different functions when responding to developmental and environmental signals in the sunflower. Exogenous ABA treatment showed that most HaPYLs respond to an increase in the ABA level. Among these HaPYLs, HaPYL2a, HaPYL4d, HaPYL4g, HaPYL8a, HaPYL8b, HaPYL8c, HaPYL9b, and HaPYL9c were up-regulated with PEG6000 treatment and NaCl treatment. This indicates that they may play a role in resisting drought and salt stress in the sunflower by mediating ABA signaling. Our findings provide some clues to further explore the functions of PYL genes in the sunflower, especially with regards to drought and salt stress resistance.

bacLIFE: a user-friendly computational workflow for genome analysis and prediction of lifestyle-associated genes in bacteria

Bacteria have an extensive adaptive ability to live in close association with eukaryotic hosts, exhibiting detrimental, neutral or beneficial effects on host growth and health. However, the genes involved in niche adaptation are mostly unknown and their functions poorly characterized. Here, we present bacLIFE ( https://github.com/Carrion-lab/bacLIFE ) a streamlined computational workflow for genome annotation, large-scale comparative genomics, and prediction of lifestyle-associated genes (LAGs). As a proof of concept, we analyzed 16,846 genomes from the Burkholderia/Paraburkholderia and Pseudomonas genera, which led to the identification of hundreds of genes potentially associated with a plant pathogenic lifestyle. Site-directed mutagenesis of 14 of these predicted LAGs of unknown function, followed by plant bioassays, showed that 6 predicted LAGs are indeed involved in the phytopathogenic lifestyle of Burkholderia plantarii and Pseudomonas syringae pv. phaseolicola. These 6 LAGs encompassed a glycosyltransferase, extracellular binding proteins, homoserine dehydrogenases and hypothetical proteins. Collectively, our results highlight bacLIFE as an effective computational tool for prediction of LAGs and the generation of hypotheses for a better understanding of bacteria-host interactions.

Identifying the suite of genes central to swimming in the biocontrol bacterium Pseudomonas protegens Pf-5

Swimming motility is a key bacterial trait, important to success in many niches. Biocontrol bacteria, such as Pseudomonas protegens Pf-5, are increasingly used in agriculture to control crop diseases, where motility is important for colonization of the plant rhizosphere. Swimming motility typically involves a suite of flagella and chemotaxis genes, but the specific gene set employed for both regulation and biogenesis can differ substantially between organisms. Here we used transposon-directed insertion site sequencing (TraDIS), a genome-wide approach, to identify 249 genes involved in P. protegens Pf-5 swimming motility. In addition to the expected flagella and chemotaxis, we also identified a suite of additional genes important for swimming, including genes related to peptidoglycan turnover, O-antigen biosynthesis, cell division, signal transduction, c-di-GMP turnover and phosphate transport, and 27 conserved hypothetical proteins. Gene knockout mutants and TraDIS data suggest that defects in the Pst phosphate transport system lead to enhanced swimming motility. Overall, this study expands our knowledge of pseudomonad motility and highlights the utility of a TraDIS-based approach for analysing the functions of thousands of genes. This work sets a foundation for understanding how swimming motility may be related to the inconsistency in biocontrol bacteria performance in the field.