InterPro in 2019: improving coverage, classification and access to protein sequence annotations

Comprehensive physiological and genomic characterization of a potential probiotic strain, Lactiplantibacillus plantarum ILSF15, isolated from the gut of tribes of Odisha, India

Characterizing probiotic features of organisms isolated from diverse environments can lead to the discovery of novel strains with promising functional features and health attributes. The present study attempts to characterize a novel probiotic strain isolated from the gut of the tribal population of Odisha, India. Based on 16S rRNA-based phylogeny, the strain was identified as a species of the Lactiplantibacillus genus and was named Lactiplantibacillus plantarum strain ILSF15. The current investigation focuses on elucidating this strain's genetic and physiological properties associated with probiotic attributes such as biosafety risk, host adaptation/survival traits, and beneficial functional features. The novel strain was observed, in vitro, exhibiting features such as acid/bile tolerance, adhesion to the host enteric epithelial cells, cholesterol assimilation, and pathogen exclusion, indicating its ability to survive the harsh environment of the human GIT and resist the growth of harmful microorganisms. Additionally, the L. plantarum ILSF15 strain was found to harbor genes associated with the metabolism and synthesis of various bioactive molecules, including amino acids, carbohydrates, lipids, and vitamins, highlighting the organism's ability to efficiently utilize diverse resources and contribute to the host's nutrition and health. Several genes involved in host adaptation/survival strategies and host-microbe interactions were also identified from the ILSF15 genome. Moreover, L. plantarum strains, in general, were found to have an open pangenome characterized by high genetic diversity and the absence of specific lineages associated with particular habitats, signifying its versatile nature and potential applications in probiotic and functional food industries.

Nature-based bioreactors: Tackling antibiotic resistance in urban wastewater treatment

The overuse and misuse of antibiotics have accelerated the selection of antibiotic-resistant bacteria, significantly impacting human, animal, and environmental health. As aquatic environments are vulnerable to antibiotic resistance, suitable management practices should be adopted to tackle this phenomenon. Here we show an effective, nature-based solution for reducing antibiotic resistance from actual wastewater. We utilize a bioreactor that relies on benthic (biofilms) and planktonic microbial communities to treat secondary effluent from a small urban wastewater treatment plant (<10,000 population equivalent). This treated effluent is eventually released into the local aquatic ecosystem. We observe high removal efficiency for genes that provide resistance to commonly used antibiotic families, as well as for mobile genetic elements that could potentially aid in their spread. Importantly, we notice a buildup of sulfonamide (sul1 and sul2) and tetracycline (tet(C), tet(G), and tetR) resistance genes specifically in biofilms. This advancement marks the initial step in considering this bioreactor as a nature-based, cost-effective tertiary treatment option for small UWWTPs facing antibiotic resistance challenges.

Applications of bacteriophage in combination with nisin for controlling multidrug-resistant Bacillus cereus in broth and various food matrices

This study focused on the isolation and characterization of bacteriophages with specific activity against toxin-producing and multidrug-resistant strains of Bacillus cereus sensu stricto (B. cereus s. s.). Ten different samples yielded six bacteriophages by utilizing the double-layer agar technique. The most promising phage, vB_BceS-M2, was selected based on its broad host range and robust lytic activity against various B. cereus s. s. strains. The phage vB_BceS-M2 had a circular double-stranded DNA genome of 56,482 bp. This phage exhibited stability over a wide range of temperatures and pH values, which is crucial for its potential application in food matrices. The combined effect of phage vB_BceS-M2 and nisin, a widely used antimicrobial peptide, was investigated to enhance antimicrobial efficacy against B. cereus in food. The results suggested that nisin showed synergy and combined effect with the phage, potentially overcoming the growth of phage-resistant bacteria in the broth. Furthermore, practical applications were conducted in various liquid and solid food matrices, including whole and skimmed milk, boiled rice, cheese, and frozen meatballs, both at 4 and 25 °C. Phage vB_BceS-M2, either alone or in combination with nisin, reduced the growth rate of B. cereus in foods other than whole milk. The combination of bacteriophage and nisin showed promise for the development of effective antimicrobial interventions to counteract toxigenic and antibiotic-resistant B. cereus in food.

Advancing microbial production through artificial intelligence-aided biology

Microbial cell factories (MCFs) have been leveraged to construct sustainable platforms for value-added compound production. To optimize metabolism and reach optimal productivity, synthetic biology has developed various genetic devices to engineer microbial systems by gene editing, high-throughput protein engineering, and dynamic regulation. However, current synthetic biology methodologies still rely heavily on manual design, laborious testing, and exhaustive analysis. The emerging interdisciplinary field of artificial intelligence (AI) and biology has become pivotal in addressing the remaining challenges. AI-aided microbial production harnesses the power of processing, learning, and predicting vast amounts of biological data within seconds, providing outputs with high probability. With well-trained AI models, the conventional Design-Build-Test (DBT) cycle has been transformed into a multidimensional Design-Build-Test-Learn-Predict (DBTLP) workflow, leading to significantly improved operational efficiency and reduced labor consumption. Here, we comprehensively review the main components and recent advances in AI-aided microbial production, focusing on genome annotation, AI-aided protein engineering, artificial functional protein design, and AI-enabled pathway prediction. Finally, we discuss the challenges of integrating novel AI techniques into biology and propose the potential of large language models (LLMs) in advancing microbial production.

AcrR1, a novel TetR/AcrR family repressor, mediates acid and antibiotic resistance and nisin biosynthesis in Lactococcus lactis F44

Lactococcus lactis, widely used in the manufacture of dairy products, encounters various environmental stresses both in natural habitats and during industrial processes. It has evolved intricate machinery of stress sensing and defense to survive harsh stress conditions. Here, we identified a novel TetR/AcrR family transcription regulator, designated AcrR1, to be a repressor for acid and antibiotic tolerance that was derepressed in the presence of vancomycin or under acid stress. The survival rates of acrR1 deletion strain ΔAcrR1 under acid and vancomycin stresses were about 28.7-fold (pH 3.0, HCl), 8.57-fold (pH 4.0, lactic acid) and 2.73-fold (300 ng/mL vancomycin) greater than that of original strain F44. We also demonstrated that ΔAcrR1 was better able to maintain intracellular pH homeostasis and had a lower affinity to vancomycin. No evident effects of AcrR1 deletion on the growth and morphology of strain F44 were observed. Subsequently, we characterized that the transcription level of genes associated with amino acids biosynthesis, carbohydrate transport and metabolism, multidrug resistance, and DNA repair proteins significantly upregulated in ΔAcrR1 using transcriptome analysis and quantitative reverse transcription-PCR assays. Additionally, AcrR1 could repress the transcription of the nisin post-translational modification gene, nisC, leading to a 16.3% increase in nisin yield after AcrR1 deletion. Our results not only refined the knowledge of the regulatory mechanism of TetR/AcrR family regulator in L. lactis, but presented a potential strategy to enhance industrial production of nisin.

ADP-glucose pyrophosphorylase gene family in soybean and implications in drought stress tolerance

BACKGROUND

ADP-glucose pyrophosphorylase (AGPase) is the key rate-limiting enzyme in starch biosynthesis pathway, and has been identified as a potential target for manipulation strategies aimed at improving crop yield and quality.

OBJECTIVE

To identify the AGPase gene family members in soybean, and explore the potential implications of GmAGPS2 in drought stress tolerance.

METHODS

The genome-wide identification and sequence analysis of soybean AGPase gene family was carried out by bioinformatics methods. The GmAGP gene expression was analyzed using transcriptome data and quantitative real-time PCR (qRT-PCR). Furthermore, transgenic yeast strains overexpressing GmAGPS2 were generated, and their growth was observed under drought stress.

RESULTS

In this study, we searched for AGPase genes (GmAGP) in the soybean genome and identified a total of 14 GmAGP genes. The GmAGP proteins had a unique conserved NTP_transferase domain and were mainly located in the chloroplast and cytosol. Evolutionarily, the GmAGP proteins can be clustered into two distinct subgroups; within the same subgroup, they displayed a similar distribution pattern of conserved motifs. The GmAGP genes exhibited an uneven distribution on 10 chromosomes, and segmental duplication contributed to AGPase gene family expansion in soybean. The GmAGP genes presented different tissue expression pattern, in which GmAGPL6, GmAGPL9, and GmAGPL10 mainly exhibited tissue-specific expression pattern. The promoter of GmAGP genes had multiple cis-acting elements related to phytohormones and stress responses, and 8 GmAGP genes contained drought-responsive cis-acting elements. qRT‒PCR analysis demonstrated a significant upregulation expression of GmAGPL6, GmAGPL10, and GmAGPS2 in response to drought stress. Further functional analysis indicated that GmAGPS2 gene could improve yeast growth under drought stress conditions and enhance the drought tolerance of yeast.

CONCLUSION

These results will contribute to further elucidation of the function of GmAGP genes, and offer important candidate genes for the genetic improvement of starch and yield-related traits and the breeding of high drought stress tolerance varieties in soybean.

Simulated microgravity impairs human NK cell cytotoxic activity against space radiation-relevant leukemic cells

Natural killer (NK) cells are an important first-line of defense against malignant cells. Because of the potential for increased cancer risk from astronaut exposure to space radiation, we determined whether microgravity present during spaceflight affects the body's defenses against leukemogenesis. Human NK cells were cultured for 48 h under normal gravity and simulated microgravity (sμG), and cytotoxicity against K-562 (CML) and MOLT-4 (T-ALL) cells was measured using standard methodology or under continuous sμG. This brief exposure to sμG markedly reduced NK cytotoxicity against both leukemias, and these deleterious effects were more pronounced in continuous sμG. RNA-seq performed on NK cells from two additional healthy donors provided insight into the mechanism(s) by which sμG reduced cytotoxicity. Given our prior report of space radiation-induced human T-ALL in vivo, the reduced cytotoxicity against MOLT-4 is striking and raises the possibility that μG may increase astronaut risk of leukemogenesis during prolonged missions beyond LEO.

Administration of intratumoral GD2-directed interleukin-2 immunocytokine and local radiation therapy to activate immune rejection of spontaneous canine melanoma

Canine malignant melanoma provides a clinically relevant, large animal parallel patient population to study the GD2-reactive hu14.18-IL-2 immunocytokine as it is similar to human melanoma and expresses GD2. The objectives of this study were to evaluate safety, radiation fractionation, and identify informative biomarkers of an in-situ tumor vaccine involving local radiation therapy plus intratumoral-immunocytokine in melanoma tumor-bearing dogs. Twelve dogs (six dogs/arm) with locally advanced or metastatic melanoma were randomized to receive a single 8 Gy fraction (arm A) or three 8 Gy fractions over 1 week (arm B) to the primary site and regional lymph nodes (when clinically involved) with the single or last fraction 5 days before intratumoral-immunocytokine at 12 mg/m 2 on 3 consecutive days. Serial tumor biopsies were obtained. All 12 dogs completed protocol treatment, and none experienced significant or unexpected adverse events. Evidence of antitumor activity includes one dog with a complete response at day 60, one dog with a partial response at day 60, and four dogs with mixed responses. Histology of serial biopsies shows a variably timed increase in intratumoral lymphocytic inflammation in some dogs. Canine NanoString analyses of serial biopsies identified changes in gene signatures of innate and adaptive cell types versus baseline. There were no significant differences in NanoString results between arm A and arm B. We conclude that intratumoral-immunocytokine in combination with local radiation therapy in canine melanoma is well tolerated and has antitumor activity with the potential to inform clinical development in melanoma patients.

A chromosome-level genome assembly provides insights into the local adaptation of Tamarix austromongolica in the Yellow River Basin, China

Tamarix austromongolica is endemic to the Yellow River Basin and has adapted to diverse ecological settings in the region, including the arid areas of northwestern China and the saline soil regions of the Yellow River Delta. However, the genetic basis of its local adaptation remains unclear. We report a chromosome-level assembly of the T. austromongolica genome based on PacBio high-fidelity sequencing and Hi-C technology. The 12 pseudochromosomes cover 98.44% of the 1.32 Gb assembly, with a contig N50 of 52.57 Mb and a BUSCO score of 98.2%. The genome comprises 913.6 Mb (68.83%) of repetitive sequences and 22,374 protein-coding genes. Genome evolution analyses suggest that genes under positive selection and significantly expanded gene families have facilitated T. austromongolica's adaptability to diverse environmental factors and high resistance to diseases. Using genotyping-by-sequencing, we conducted population structure and selection analyses of 114 samples from 15 sites. Two genetic groups were identified, and 114 and 289 candidate genes were assigned to the populations of the northwestern and eastern parts of the Yellow River, respectively. Furthermore, we discovered numerous candidate genes associated with high-altitude adaptability and salt tolerance. This research provides valuable genomic resources for the evolutionary study and genetic breeding of tamarisk.

"Tear down that wall"-a critical evaluation of bioinformatic resources available for lysin researchers

Lytic enzymes, or lysins for short, break down peptidoglycan and interrupt the continuity of the cell wall, which, in turn, causes osmotic lysis of the bacterium. Their ability to destroy bacteria from within makes them promising antimicrobial agents that can be used as alternatives or supplements to antibiotics. In this paper, we briefly summarize basic terms and concepts used to describe lysin sequences and delineate major lysin groups. More importantly, we describe the domain repertoire found in lysins and critically review bioinformatic tools or databases which are used in studies of these enzymes (with particular emphasis on the repositories of Hidden Markov models). Finally, we present a novel comprehensive, meticulously curated set of lysin-related family and domain models, sort them into clusters that reflect major families, and demonstrate that the selected models can be used to efficiently search for new lysins.

Chemical mutagenesis of Listeria monocytogenes for increased tolerance to benzalkonium chloride shows independent genetic underpinnings and off-target antibiotic resistance

Listeria monocytogenes, a potentially fatal foodborne pathogen commonly found in food processing facilities, creates a significant economic burden that totals more than $2 billion annually in the United States due to outbreaks. Quaternary ammonium compounds (QACs), including benzalkonium chloride (BAC), are among the most widely used sanitizers to inhibit the growth and spread of L. monocytogenes from food processing facilities. However, resistance to QACs has been increasing in L. monocytogenes and different genetic mechanisms conferring resistance have been discovered. Here, we used ethyl methanesulfonate (EMS) to chemically mutagenize the BAC-susceptible strain, L. monocytogenes FSL-N1-304. We isolated two mutants with increased tolerance to BAC compared to the parental strain. Next, we assessed the off-target effect of increased tolerance to BAC by measuring the minimum inhibitory concentrations (MICs) of a diverse set of antibiotics, revealing that mut-1 and mut-2 displayed significantly increased resistance to fluoroquinolone antibiotics compared to the parental strain. A hemolysis assay was then used to investigate a potential correlation between BAC tolerance and virulence. Interestingly, mut-1 and mut-2 both exhibited significantly higher hemolysis percentage than the parental strain. We then sequenced the genomes of the parental strain and both mutants to identify mutations that may be involved in the increased resistance to BAC. We identified 3 and 29 mutations in mut-1 and mut-2, respectively. mut-1 contained nonsynonymous mutations in dagK (a diacylglycerol kinase), lmo2768 (a permease-encoding gene), and lmo0186 (resuscitation promoting factor). mut-2 contained a nonsense mutation in the nucleotide excision repair enzyme UvrABC system protein B encoding gene, uvrB, which likely accounts for the higher number of mutations observed. Transcriptome analysis in the presence of BAC revealed that genes related to the phosphotransferase system and internalins were up-regulated in both mutants, suggesting their significance in the BAC stress response. These two mutants provide insights into alternative mechanisms for increased BAC tolerance and could further our understanding of how L. monocytogenes persists in the food processing environment.

Small molecule FICD inhibitors suppress endogenous and pathologic FICD-mediated protein AMPylation

The AMP transferase, FICD, is an emerging drug target finetuning stress signaling in the endoplasmic reticulum (ER). FICD is a bi-functional enzyme, catalyzing both AMP addition (AMPylation) and removal (deAMPylation) from the ER resident chaperone BiP/GRP78. Despite increasing evidence linking excessive BiP/GRP78 AMPylation to human diseases, small molecules to inhibit pathogenic FICD variants are lacking. Using an in-vitro high-throughput screen, we identify two small-molecule FICD inhibitors, C22 and C73. Both molecules significantly inhibit FICD-mediated BiP/GRP78 AMPylation in intact cells while only weakly inhibiting BiP/GRP78 deAMPylation. C22 and C73 also efficiently inhibit pathogenic FICD variants and improve proinsulin processing in β cells. Our study identifies and validates FICD inhibitors, highlighting a novel therapeutic avenue against pathologic protein AMPylation.

Novel splicing variant and gonadal mosaicism in DYRK1A gene identified by whole-genome sequencing in multiplex autism spectrum disorder families

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with considerable genetic heterogeneity. The disorder is clinically diagnosed based on DSM-5 criteria, featuring deficits in social communication and interaction, along with restricted and repetitive behaviours. Here, we performed whole-genome sequencing (WGS) on four individuals with ASD from two multiplex families (MPX), where more than one individual is affected, to identify potential single nucleotide variants (SNVs) and structural variants (SVs) in coding and non-coding regions. A rigorous bioinformatics pipeline was employed for variant detection, followed by segregation analysis. Our investigation revealed an unreported splicing variant in the DYRK1A gene (c.-77 + 2T > C; IVS1 + 2T > C; NM_001396.5), in heterozygote form in two affected children in one of the families (family B), which was absent in the healthy parents and siblings. This finding suggests the presence of gonadal mosaicism in one of the parents, representing the first documented instance of such inheritance for a variant in the DYRK1A gene associated with ASD. Furthermore, we identified a 50 bp deletion in intron 9 of the DLG2 gene in two affected patients from the same family, confirmed by PCR and Sanger sequencing. In Family A, we identified potential candidate variants associated with ASD shared by the two patients. These findings enhance our understanding of the genetic landscape of ASD, particularly in MPX families, and highlight the utility of WGS in uncovering novel genetic contributions to neurodevelopmental disorders.

Evolutionary trends indicate a coherent organization of sap operons

Human hosts possess a complex network of immune responses against microbial pathogens. The production of antimicrobial peptides (AMPs), which target the pathogen cell membranes and inhibit them from inhabiting the hosts, is one such mechanism. However, pathogens have evolved systems that encounter these host-produced AMPs. The Sap (sensitivity to antimicrobial peptides) transporter uptakes AMPs inside the microbial cell and proteolytically degrades them. The Sap transporters comprise five subunits encoded by genes in an operon. Despite its ubiquitous nature, its subunits are not found to be in tandem with many organisms. In this study, a total of 421 Sap transporters were analyzed for their operonic arrangement. Out of 421, a total of 352 operons were found to be in consensus arrangement, while the remaining 69 show a varying arrangement of genes. The analysis of the intergenic distance between the subunits of the sap operon suggests a signature pattern with sapAB (-4), sapBC (-14), sapCD (-1), and sapDF (-4 to 1). An evolutionary analysis of these operons favors the consensus arrangement of the Sap transporter systems, substantiating its prevalence in most of the Gram-negative pathogens. Overall, this study provides insight into bacterial evolution, favoring the maintenance of the genetic organization of essential pathogenicity factors.

Enhanced Genomic and Transcriptomic Resources for Trichinella pseudospiralis and T. spiralis to Underpin the Discovery of Molecular Differences between Stages and Species

Nematodes of the genus Trichinella are important pathogens of humans and animals. This study aimed to enhance the genomic and transcriptomic resources for T. pseudospiralis (non-encapsulated phenotype) and T. spiralis (encapsulated phenotype) and to explore transcriptional profiles. First, we improved the assemblies of the genomes of T. pseudospiralis (code ISS13) and T. spiralis (code ISS534), achieving genome sizes of 56.6 Mb (320 scaffolds, and an N50 of 1.02 Mb) and 63.5 Mb (568 scaffolds, and an N50 value of 0.44 Mb), respectively. Then, for each species, we produced RNA sequence data for three key developmental stages (first-stage muscle larvae [L1s], adults, and newborn larvae [NBLs]; three replicates for each stage), analysed differential transcription between stages, and explored enriched pathways and processes between species. Stage-specific upregulation was linked to cellular processes, metabolism, and host-parasite interactions, and pathway enrichment analysis showed distinctive biological processes and cellular localisations between species. Indeed, the secreted molecules calmodulin, calreticulin, and calsyntenin-with possible roles in modulating host immune responses and facilitating parasite survival-were unique to T. pseudospiralis and not detected in T. spiralis. These insights into the molecular mechanisms of Trichinella-host interactions might offer possible avenues for developing new interventions against trichinellosis.

Sirenian genomes illuminate the evolution of fully aquatic species within the mammalian superorder afrotheria

Sirenians of the superorder Afrotheria were the first mammals to transition from land to water and are the only herbivorous marine mammals. Here, we generated a chromosome-level dugong (Dugong dugon) genome. A comparison of our assembly with other afrotherian genomes reveals possible molecular adaptations to aquatic life by sirenians, including a shift in daily activity patterns (circadian clock) and tolerance to a high-iodine plant diet mediated through changes in the iodide transporter NIS (SLC5A5) and its co-transporters. Functional in vitro assays confirm that sirenian amino acid substitutions alter the properties of the circadian clock protein PER2 and NIS. Sirenians show evidence of convergent regression of integumentary system (skin and its appendages) genes with cetaceans. Our analysis also uncovers gene losses that may be maladaptive in a modern environment, including a candidate gene (KCNK18) for sirenian cold stress syndrome likely lost during their evolutionary shift in daily activity patterns. Genomes from nine Australian locations and the functionally extinct Okinawan population confirm and date a genetic break ~10.7 thousand years ago on the Australian east coast and provide evidence of an associated ecotype, and highlight the need for whole-genome resequencing data from dugong populations worldwide for conservation and genetic management.

Drug-target interaction prediction using knowledge graph embedding

The prediction of drug-target interactions (DTIs) is a critical phase in the sustainable drug development process, especially when the research focus is to capitalize on the repositioning of existing drugs. Computational approaches to predicting DTIs can provide important insights into drug mechanisms of action. However, current methods for predicting DTIs based on the structural information of the knowledge graph may suffer from the sparseness and incompleteness of the knowledge graph and neglect the latent type information of the knowledge graph. In this paper, we propose TTModel, a knowledge graph embedding model for DTI prediction. By exploiting biomedical text and type information, TTModel can learn latent text semantics and type information to improve the performance of representation learning. Comprehensive experiments on two public datasets demonstrate that our model outperforms the state-of-the-art methods significantly on the task of DTI prediction.

Genetic variation of TLR3 gene is associated with the outcome of hepatitis b infection in mauritanian patients: case control study

BACKGROUND

Toll-Like receptors (TLRs) play an important role in the immune response during hepatitis B virus (HBV) infection. In this study, we evaluated the association between two SNP variants (TLR3 rs3775290 and TLR4 rs4986790) and susceptibility to chronic HBV infection in Mauritania.

SUBJECTS AND METHODS

A total of 188 subjects were recruited for this study: 102 chronically infected patients and 86 individuals with spontaneously resolved HBV infection who were considered controls. Targeted PCR products were sequenced using Sanger sequencing.

RESULTS

We found that TLR3 rs3775290 was significantly more frequent in patients with chronic HBV than in the control population (p = 0.03). However, no association was found between the TLR4 rs3775290 polymorphism and chronic infection.

CONCLUSION

Our results suggest that the TLR3 rs3775290 polymorphism may be a risk factor for susceptibility to chronic HBV infection in the Mauritanian population.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

Summary

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based translated sequence annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long-read sequencing data and in the context of pseudogenes.

Availability and implementation

The software is available at https://github.com/TravisWheelerLab/BATH.

A genomic compendium of cultivated human gut fungi characterizes the gut mycobiome and its relevance to common diseases

The gut fungal community represents an essential element of human health, yet its functional and metabolic potential remains insufficiently elucidated, largely due to the limited availability of reference genomes. To address this gap, we presented the cultivated gut fungi (CGF) catalog, encompassing 760 fungal genomes derived from the feces of healthy individuals. This catalog comprises 206 species spanning 48 families, including 69 species previously unidentified. We explored the functional and metabolic attributes of the CGF species and utilized this catalog to construct a phylogenetic representation of the gut mycobiome by analyzing over 11,000 fecal metagenomes from Chinese and non-Chinese populations. Moreover, we identified significant common disease-related variations in gut mycobiome composition and corroborated the associations between fungal signatures and inflammatory bowel disease (IBD) through animal experimentation. These resources and findings substantially enrich our understanding of the biological diversity and disease relevance of the human gut mycobiome.

Chromosome-level genome assembly of Solanum pimpinellifolium

Solanum pimpinellifolium, the closest wild relative of the domesticated tomato, has high potential for use in breeding programs aimed at developing multi-pathogen resistance and quality improvement. We generated a chromosome-level genome assembly of S. pimpinellifolium LA1589, with a size of 833 Mb and a contig N50 of 31 Mb. We anchored 98.80% of the contigs into 12 pseudo-chromosomes, and identified 74.47% of the sequences as repetitive sequences. The genome evaluation revealed BUSCO and LAI score of 98.3% and 14.49, respectively, indicating high quality of this assembly. A total of 41,449 protein-coding genes were predicted in the genome, of which 89.17% were functionally annotated. This high-quality genome assembly serves as a valuable resource for accelerating the biological discovery and molecular breeding of this important horticultural crop.

HortGenome Search Engine, a universal genomic search engine for horticultural crops

Horticultural crops comprising fruit, vegetable, ornamental, beverage, medicinal and aromatic plants play essential roles in food security and human health, as well as landscaping. With the advances of sequencing technologies, genomes for hundreds of horticultural crops have been deciphered in recent years, providing a basis for understanding gene functions and regulatory networks and for the improvement of horticultural crops. However, these valuable genomic data are scattered in warehouses with various complex searching and displaying strategies, which increases learning and usage costs and makes comparative and functional genomic analyses across different horticultural crops very challenging. To this end, we have developed a lightweight universal search engine, HortGenome Search Engine (HSE; http://hort.moilab.net), which allows for the querying of genes, functional annotations, protein domains, homologs, and other gene-related functional information of more than 500 horticultural crops. In addition, four commonly used tools, including 'BLAST', 'Batch Query', 'Enrichment analysis', and 'Synteny Viewer' have been developed for efficient mining and analysis of these genomic data.

Checkpoint inhibitor-induced lichen planus differs from spontaneous lichen planus on the clinical, histological, and gene expression level

Background

Although highly efficacious, immune checkpoint inhibitors induce a multitude of immune-related adverse events including lichenoid skin reactions (irLP) that are often therapy-resistant.

Objectives

To compare the clinical, histological, and transcriptional features of irLP with spontaneous lichen planus (LP).

Methods

Clinical and histological presentations of irLP and LP, as well as the gene expression profiles of irLP and LP lesional and healthy skin were assessed.

Results

irLP differed considerably from LP with regard to the distribution pattern of skin lesions with irLP appearing mostly in an exanthematous form, whereas lesions were more localized in the LP group. Histologically, dermal lymphocyte infiltration was significantly lower in irLP compared with LP, whereas lymphocyte exocytosis and apoptotic keratinocytes were significantly higher in irLP. Gene expression analysis revealed irLP to have a more inflammatory profile with elevated IFNG levels and a possible role of phagosome signaling compared with LP.

Limitations

The study is descriptive and necessitates further investigation with larger cohorts and broader analyses.

Conclusion

irLP differs from spontaneous LP on the clinical, histopathological, and gene expression level. The inflammatory gene signature in irLP suggests that topical JAK inhibitors could be an effective treatment, targeting local skin inflammation without systemic immunosuppression.

POLAT: Protein function prediction based on soft mask graph network and residue-Label ATtention

MOTIVATION

Elucidating protein function is a central problem in biochemistry, genetics, and molecular biology. Developing computational methods for protein function prediction is critical due to the significant gap between sequence and functional data. Recent advances in protein structure prediction, which strongly correlates with function, make it feasible to use structure to predict function. However, current structure-based methods overlook the fact that individual residues may contribute differently to the protein's function and do not take into account the correlation between protein residues and their functions. The challenge of effectively utilizing the relationship between protein residues and function-level information to predict protein function remains unsolved.

RESULT

We proposed a protein function prediction method based on Soft Mask Graph Networks and Residue-Label Attention (POLAT), which could combine sequence features, predicted structure features, and function-level information to get an accurate prediction. We use soft mask graph networks to adaptively extract the residues relevant to functions. A residue-label attention mechanism is adopted to obtain the protein-level encoded features of a protein, which are then concatenated with a protein-level embedding and fed into a dense classifier to determine the probabilities of each function. POLAT achieves 0.670, 0.515, 0.578 Fmax and 0.677, 0.409, 0.507 AUPR on the PDB cdhit test set for the MFO, BPO, and CCO domains, respectively, outperforming the existing structure-based SOTA method GAT-GO (Fmax 0.633, 0.492, 0.547; AUPR 0.660, 0.381, 0.479). POLAT is also competitive in extensive experiments among sequence-based and multimodal methods and achieves the SOTA performance in three out of six metrics.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape-all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs-questioning the differences between them and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the consensus protein derived from our full Ribonuclease H sequence alignment is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted set of sequences is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order sequence correlations using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Identification and Characterization of Three Novel Solemo-like Viruses in the White-Backed Planthopper, Sogatella furcifera

Agricultural insects play a crucial role in transmitting plant viruses and host a considerable number of insect-specific viruses (ISVs). Among these insects, the white-backed planthoppers (WBPH; Sogatella furcifera, Hemiptera: Delphacidae) are noteworthy rice pests and are responsible for disseminating the southern rice black-streaked dwarf virus (SRBSDV), a significant rice virus. In this study, we analyzed WBPH transcriptome data from public sources and identified three novel viruses. These newly discovered viruses belong to the plant-associated viral family Solemoviridae and were tentatively named Sogatella furcifera solemo-like virus 1-3 (SFSolV1-3). Among them, SFSolV1 exhibited a prevalent existence in different laboratory populations, and its complete genome sequence was obtained using rapid amplification of cDNA ends (RACE) approaches. To investigate the antiviral RNA interference (RNAi) response in WBPH, we conducted an analysis of virus-derived small interfering RNAs (vsiRNAs). The vsiRNAs of SFSolV1 and -2 exhibited typical patterns associated with the host's siRNA-mediated antiviral immunity, with a preference for 21- and 22-nt vsiRNAs derived equally from both the sense and antisense genomic strands. Furthermore, we examined SFSolV1 infection and distribution in WBPH, revealing a significantly higher viral load of SFSolV1 in nymphs' hemolymph compared to other tissues. Additionally, in adult insects, SFSolV1 exhibited higher abundance in male adults than in female adults.

Enhanced amphiregulin exposure promotes modulation of the high grade serous ovarian cancer tumor immune microenvironment

High grade serous ovarian cancer (HGSOC) is a lethal gynecologic malignancy in which chemoresistant recurrence rates remain high. Furthermore, HGSOC patients have demonstrated overall low response rates to clinically available immunotherapies. Amphiregulin (AREG), a low affinity epidermal growth factor receptor ligand is known to be significantly upregulated in HGSOC patient tumors following neoadjuvant chemotherapy exposure. While much is known about AREG's role in oncogenesis and classical immunity, it is function in tumor immunology has been comparatively understudied. Therefore, the objective of this present study was to elucidate how increased AREG exposure impacts the ovarian tumor immune microenvironment (OTIME). Using NanoString IO 360 and protein analysis, it was revealed that treatment with recombinant AREG led to prominent upregulation of genes associated with ovarian pathogenesis and immune evasion (CXCL8, CXCL1, CXCL2) along with increased STAT3 activation in HGSOC cells. In vitro co-culture assays consisting of HGSOC cells and peripheral blood mononuclear cells (PBMCs) stimulated with recombinant AREG (rAREG) led to significantly enhanced tumor cell viability. Moreover, PBMCs stimulated with rAREG exhibited significantly lower levels of IFNy and IL-2. In vivo rAREG treatment promoted significant reductions in circulating levels of IL-2 and IL-5. Intratumoral analysis of rAREG treated mice revealed a significant reduction in CD8+ T cells coupled with an upregulation of PD-L1. Finally, combinatorial treatment with an AREG neutralizing antibody and carboplatin led to a synergistic reduction of cell viability in HGSOC cell lines OVCAR8 and PEA2. Overall, this study demonstrates AREG's ability to modulate cytotoxic responses within the OTIME and highlights its role as a novel HGSOC immune target.

Genomic identification and evolutionary analysis of chemosensory receptor gene families in two Phthorimaea pest species: insights into chemical ecology and host adaptation

BACKGROUND

Insects rely on sophisticated sensitive chemosensory systems to sense their complex chemical environment. This sensory process involves a combination of odorant receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs) in the chemosensory system. This study focused on the identification and characterization of these three types of chemosensory receptor genes in two closely related Phthorimaea pest species, Phthorimaea operculella (potato tuber moth) and Phthorimaea absoluta (tomato leaf miner).

RESULTS

Based on manual annotation of the genome, we identified a total of 349 chemoreceptor genes from the genome of P. operculella, including 93 OR, 206 GR and 50 IR genes, while for P. absoluta, we identified 72 OR, 122 GR and 46 IR genes. Through phylogenetic analysis, we observed minimal differences in the number and types of ORs and IRs between the potato tuber moth and tomato leaf miner. In addition, we found that compared with those of tomato leaf miners, the gustatory receptor branch of P. operculella has undergone a large expansion, which may be related to P. absoluta having a narrower host range than P. operculella. Through analysis of differentially expressed genes (DEGs) of male and female antennae, we uncovered 45 DEGs (including 32ORs, 9 GRs, and 4 IRs).

CONCLUSIONS

Our research provides a foundation for exploring the chemical ecology of these two pests and offers new insights into the dietary differentiation of lepidopteran insects, while simultaneously providing molecular targets for developing environmentally friendly pest control methods based on insect chemoreception.

Assembly and analysis of the genome of Notholithocarpus densiflorus

Tanoak (Notholithocarpus densiflorus) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes, and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks (Lithocarpus spp.), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks (Quercus) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to "sudden oak death" (SOD), a plant pathogen (Phytophthora ramorum) that has caused widespread deaths of tanoaks. In this study, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using cobarcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore long reads to 30× coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published to date, with a contig N50 of ∼1.2 Mb, a scaffold N50 of ∼2.1 Mb, and a complete gene score of 95.5% through BUSCO analysis. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome, enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using cobarcoded data, we were able to generate a complete haplotype coverage of all 11 genomes.

Protein function prediction through multi-view multi-label latent tensor reconstruction

BACKGROUND

In last two decades, the use of high-throughput sequencing technologies has accelerated the pace of discovery of proteins. However, due to the time and resource limitations of rigorous experimental functional characterization, the functions of a vast majority of them remain unknown. As a result, computational methods offering accurate, fast and large-scale assignment of functions to new and previously unannotated proteins are sought after. Leveraging the underlying associations between the multiplicity of features that describe proteins could reveal functional insights into the diverse roles of proteins and improve performance on the automatic function prediction task.

RESULTS

We present GO-LTR, a multi-view multi-label prediction model that relies on a high-order tensor approximation of model weights combined with non-linear activation functions. The model is capable of learning high-order relationships between multiple input views representing the proteins and predicting high-dimensional multi-label output consisting of protein functional categories. We demonstrate the competitiveness of our method on various performance measures. Experiments show that GO-LTR learns polynomial combinations between different protein features, resulting in improved performance. Additional investigations establish GO-LTR's practical potential in assigning functions to proteins under diverse challenging scenarios: very low sequence similarity to previously observed sequences, rarely observed and highly specific terms in the gene ontology.

IMPLEMENTATION

The code and data used for training GO-LTR is available at https://github.com/aalto-ics-kepaco/GO-LTR-prediction .

Structural characterization of lytic transglycosylase MltD of Pseudomonas aeruginosa, a target for the natural product bulgecin A

Natural product bulgecin A potentiates the activity of β-lactam antibiotics by inhibition of three lytic transglycosylases in Pseudomonas aeruginosa, of which MltD is one. MltD exhibits both endolytic and exolytic reactions in the turnover of the cell-wall peptidoglycan and tolerates the presence or absence of stem peptides in its substrates. The present study reveals structural features of the multimodular MltD, presenting a catalytic module and four cell-wall-binding LysM modules that account for these attributes. Three X-ray structures are reported herein for MltD that disclose one unpredicted LysM module tightly attached to the catalytic domain, whereas the other LysM modules are mobile, and connected to the catalytic domain through long flexible linkers. The formation of crystals depended on the presence of bulgecin A. The expansive active-site cleft is highlighted by the insertion of a helical region, a hallmark of the family 1D of lytic transglycosylases, which was mapped out in a ternary complex of MltD:bulgecinA:chitotetraose, revealing at the minimum the presence of eight subsites (from -4 to +4, with the seat of reaction at subsites -1 and + 1) for binding of sugars of the substrate for the endolytic reaction. The mechanism of the exolytic reaction is revealed in one of the structures, showing how the substrate's terminal anhydro-NAM moiety could be sequestered at subsite +2. Our results provide the structural insight for both the endolytic and exolytic activities of MltD during cell-wall-turnover events.

Age-Related Differences in the Mouse Corneal Epithelial Transcriptome and Their Impact on Corneal Wound Healing

Purpose

Aging is a risk factor for dry eye. We sought to identify changes in the aged mouse corneal epithelial transcriptome and determine how age affects corneal sensitivity, re-epithelialization, and barrier reformation after corneal debridement.

Methods

Corneal epithelium of female C57BL/6J (B6) mice of different ages (2, 12, 18, and 24 months) was collected, RNA extracted, and bulk RNA sequencing performed. Cornea sensitivity was measured with an esthesiometer in 2- to 3-month-old, 12- to 13-month-old, 18- to 19-month-old, and 22- to 25-month-old female and male mice. The 2-month-old and 18-month-old female and male mice underwent unilateral corneal debridement using a blunt blade. Wound size and fluorescein staining were visualized and photographed at different time points, and a re-epithelialization rate curve was calculated.

Results

There were 157 differentially expressed genes in aged mice compared with young mice. Several pathways downregulated with age control cell migration, proteoglycan synthesis, and collagen trimerization, assembly, biosynthesis, and degradation. Male mice had decreased corneal sensitivity compared with female mice at 12 and 24 months of age. Aged mice, irrespective of sex, had delayed corneal re-epithelialization in the first 48 hours and worse corneal fluorescein staining intensity at day 14 than young mice.

Conclusions

Aged corneal epithelium has an altered transcriptome. Aged mice regardless of sex heal more slowly and displayed more signs of corneal epithelial defects after wounding than young mice. These results indicate that aging significantly alters the corneal epithelium and its ability to coordinate healing.

Blue-light irradiation induced partial nitrification

The role of ray radiation from the sunlight acting on organisms has long-term been investigated. However, how the light with different wavelengths affects nitrification and the involved nitrifiers are still elusive. Here, we found more than 60 % of differentially expressed genes (DEGs) in nitrifiers were observed under irradiation of blue light with wavelengths of 440-480 nm, which were 13.4 % and 20.3 % under red light and white light irradiation respectively. Blue light was more helpful to achieve partial nitrification rather than white light or red light, where ammonium oxidization by ammonia-oxidizing archaea (AOA) with the increased relative abundance from 8.6 % to 14.2 % played a vital role. This was further evidenced by the enhanced TCA cycle, reactive oxygen species (ROS) scavenge and DNA repair capacity in AOA under blue-light irradiation. In contrast, nitrite-oxidizing bacteria (NOB) was inhibited severely to achieve partial nitrification, and the newly discovered encoded blue light photoreceptor proteins made them more sensitive to blue light and hindered cell activity. Ammonia-oxidizing bacteria (AOB) expressed genes for DNA repair capacity under blue-light irradiation, which ensured their tiny impact by light irradiation. This study provided valuable insights into the photosensitivity mechanism of nitrifiers and shed light on the diverse regulatory by light with different radiation wavelengths in artificial systems, broadening our comprehension of the nitrogen cycle on earth.

Simulated Microgravity Impairs Human NK Cell Cytotoxic Activity Against Space Radiation-Relevant Leukemic Cells

Natural killer (NK) cells are important effectors of the innate immune system. Unlike T cells, NK cells do not require antigen-priming, making them an important first-line of defense against malignant cells. Because of the potential for increased cancer risk as a result of astronaut exposure to space radiation, we performed studies to determine whether conditions of microgravity present during spaceflight affects the body's natural defenses against leukemogenesis. Human NK cells were cultured for 48 hours under normal gravity and simulated microgravity (sµG), and cytotoxicity against K-562 (CML) and MOLT-4 (T-ALL) cell lines was measured using standard methodology or under continuous conditions of sµG. Even this brief exposure to sµG markedly reduced NK cytotoxicity against both leukemic cells using standard assay procedures, and these deleterious effects were even more pronounced in continuous sµG. RNA-seq performed on NK cells from two healthy donors provided insight into the mechanism(s) by which sµG reduced cytotoxicity. Given our prior report that human HSC exposed to simulated space radiation gave rise to T-ALL in vivo , the reduced cytotoxicity against MOLT-4 is striking and raises the possibility that µG may add to astronaut risk of leukemogenesis during prolonged missions beyond LEO.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Freshwater genome-reduced bacteria exhibit pervasive episodes of adaptive stasis

The emergence of bacterial species is rooted in their inherent potential for continuous evolution and adaptation to an ever-changing ecological landscape. The adaptive capacity of most species frequently resides within the repertoire of genes encoding the secreted proteome (SP), as it serves as a primary interface used to regulate survival/reproduction strategies. Here, by applying evolutionary genomics approaches to metagenomics data, we show that abundant freshwater bacteria exhibit biphasic adaptation states linked to the eco-evolutionary processes governing their genome sizes. While species with average to large genomes adhere to the dominant paradigm of evolution through niche adaptation by reducing the evolutionary pressure on their SPs (via the augmentation of functionally redundant genes that buffer mutational fitness loss) and increasing the phylogenetic distance of recombination events, most of the genome-reduced species exhibit a nonconforming state. In contrast, their SPs reflect a combination of low functional redundancy and high selection pressure, resulting in significantly higher levels of conservation and invariance. Our findings indicate that although niche adaptation is the principal mechanism driving speciation, freshwater genome-reduced bacteria often experience extended periods of adaptive stasis. Understanding the adaptive state of microbial species will lead to a better comprehension of their spatiotemporal dynamics, biogeography, and resilience to global change.

A coordinated switch in sucrose and callose metabolism enables enhanced symplastic unloading in potato tubers

One of the early changes upon tuber induction is the switch from apoplastic to symplastic unloading. Whether and how this change in unloading mode contributes to sink strength has remained unclear. In addition, developing tubers also change from energy to storage-based sucrose metabolism. Here, we investigated the coordination between changes in unloading mode and sucrose metabolism and their relative role in tuber sink strength by looking into callose and sucrose metabolism gene expression combined with a model of apoplastic and symplastic unloading. Gene expression analysis suggests that callose deposition in tubers is decreased by lower callose synthase expression. Furthermore, changes in callose and sucrose metabolism are strongly correlated, indicating a well-coordinated developmental switch. Modelling indicates that symplastic unloading is not the most efficient unloading mode per se. Instead, it is the concurrent metabolic switch that provides the physiological conditions necessary to potentiate symplastic transport and thereby enhance tuber sink strength .

Homology modeling and molecular docking studies to decrease glutamine affinity of Yarrowia lipolytica L-asparaginase

L-Asparaginase is a key component in the treatment of leukemias and lymphomas. However, the glutamine affinity of this therapeutic enzyme is an off-target activity that causes several side effects. The modeling and molecular docking study of Yarrowia lipolytica L-asparaginase (YL-ASNase) to reduce its l-glutamine affinity and increase its stability was the aim of this study. Protein-ligand interactions of wild-type and different mutants of YL-ASNase against L-asparagine compared to l-glutamine were assessed using AutoDock Vina tools because the crystal structure of YL-ASNase does not exist in the protein data banks. The results showed that three mutants, T171S, T171S-N60A, and T171A-T223A, caused a considerable increase in L-asparagine affinity and a decrease in l-glutamine affinity as compared to the wild-type and other mutants. Then, molecular dynamics simulation and MM/GBSA free energy were applied to assess the stability of protein structure and its interaction with ligands. The three mutated proteins, especially T171S-N60A, had higher stability and interactions with L-asparagine than l-glutamine in comparison with the wild-type. The YL-ASNase mutants could be introduced as appropriate therapeutic candidates that might cause lower side effects. However, the functional properties of these mutated enzymes need to be confirmed by genetic manipulation and in vitro and in vivo studies.

A high-quality chromosome-level genome assembly of the Chinese medaka Oryzias sinensis

Oryzias sinensis, also known as Chinese medaka or Chinese ricefish, is a commonly used animal model for aquatic environmental assessment in the wild as well as gene function validation or toxicology research in the lab. Here, a high-quality chromosome-level genome assembly of O. sinensis was generated using single-tube long fragment read (stLFR) reads, Nanopore long-reads, and Hi-C sequencing data. The genome is 796.58 Mb, and a total of 712.17 Mb of the assembled sequences were anchored to 23 pseudo-chromosomes. A final set of 22,461 genes were annotated, with 98.67% being functionally annotated. The Benchmarking Universal Single-Copy Orthologs (BUSCO) benchmark of genome assembly and gene annotation reached 95.1% (93.3% single-copy) and 94.6% (91.7% single-copy), respectively. Furthermore, we also use ATAC-seq to uncover chromosome transposase-accessibility as well as related genome area function enrichment for Oryzias sinensis. This study offers a new improved foundation for future genomics research in Chinese medaka.

Phosphoregulation in the N-terminus of NRT2.1 affects nitrate uptake by controlling the interaction of NRT2.1 with NAR2.1 and kinase HPCAL1 in Arabidopsis

NRT2.1, the major high affinity nitrate transporter in roots, can be phosphorylated at five different sites within the N- and the C-terminus. Here, we characterized the functional relationship of two N-terminal phosphorylation sites, S21 and S28, in Arabidopsis. Based on a site-specific correlation network, we identified a receptor kinase (HPCAL1, AT5G49770), phosphorylating NRT2.1 at S21 and resulting in active nitrate uptake. HPCAL1 itself was regulated by phosphorylation at S839 and S870 within its kinase domain. In the active state, when S839 was dephosphorylated and S870 was phosphorylated, HPCAL1 was found to interact with the N-terminus of NRT2.1, mainly when S28 was dephosphorylated. Phosphorylation of NRT2.1 at S21 resulted in a reduced interaction of NRT2.1 with its activator NAR2.1, but nitrate transport activity remained. By contrast, phosphorylated NRT2.1 at S28 enhanced the interaction with NAR2.1, but reduced the interaction with HPCAL1. Here we identified HPCAL1 as the kinase affecting this phospho-switch through phosphorylation of NRT2.1 at S21.

A novel laccase from Trametes polyzona with high performance in the decolorization of textile dyes

Laccases are multicopper oxidases able to oxidize several phenolic compounds and find application in numerous industrial applications. Among laccase producers, white-rot fungi represent a valuable source of multiple isoforms and isoenzymes of these multicopper oxidases. Here we describe the identification, biochemical characterization, and application of laccase 2 from Trametes polyzona (TP-Lac2), a basidiomycete fungus emerged among others that have been screened by plate assay. This enzyme has an optimal temperature of 50 °C and in acidic conditions it is able to oxidize both phenolic and non-phenolic compounds. The ability of TP-Lac2 to decolorize textile dyes was tested in the presence of natural and synthetic mediators at 30 °C and 50 °C. Our results indicate that TP-Lac2 most efficiently decolorizes (decolorization rate > 75%) malachite green oxalate, orange G, amido black10B and bromocresol purple in the presence of acetosyringone and 2,2'-azinobis (3-ethylbenzthiazoline-6-sulfonate)-ABTS. Overall, the laccase mediator system consisting of TP-Lac2 and the natural mediator acetosyringone has potential as an environmentally friendly alternative for wastewater treatment in the textile industry.

Distinct characteristics of putative archaeal 5-methylcytosine RNA methyltransferases unveil their substrate specificities and evolutionary ancestries

5-Methylcytosine methyltransferases (m5C MTases) are known to be involved in the modification of RNA. Although these enzymes have been relatively well characterized in bacteria and eukarya, a complete understanding of the archaeal counterparts is lacking. In this study, the identification and characterization of archaeal RNA m5C MTases were performed. As a case study, a hyperthermophilic archaeon, Pyrococcus horikoshii OT3, which possesses five putative RNA m5C MTases, was chosen. Among the five putative RNA m5C MTases, two proteins (PH0851 and PH1991) have been characterized as homologs of a bacterial rRNA MTase (RsmB) and eukaryal tRNA MTase (NSUN6), respectively. The in-depth characterization of the remaining three putative RNA m5C MTases (PH1078, PH1374, and PH1537) in this study suggests the presence of the signature architecture and catalytic residues plausibly involved in the binding of their cognate RNA substrates. Additionally, the results also suggest the existence of two RsmB-like proteins (PH0851 and PH1078) belonging to the same subfamily IV of m5C RNA MTase. However, the proteins PH1374 and PH1537 belong to the same subfamily V but bind to different substrates, rRNA and tRNA, respectively. The findings further indicate that archaeal RNA m5C MTases link those from bacteria and eukarya.Communicated by Ramaswamy H. Sarma.

Cyclo-diphenylalanine production in Aspergillus nidulans through stepwise metabolic engineering

Cyclo-diphenylalanine (cFF) is a symmetrical aromatic diketopiperazine (DKP) found wide-spread in microbes, plants, and resulting food products. As different bioactivities continue being discovered and relevant food and pharmaceutical applications gradually emerge for cFF, there is a growing need for establishing convenient and efficient methods to access this type of compound. Here, we present a robust cFF production system which entailed stepwise engineering of the filamentous fungal strain Aspergillus nidulans A1145 as a heterologous expression host. We first established a preliminary cFF producing strain by introducing the heterologous nonribosomal peptide synthetase (NRPS) gene penP1 to A. nidulans A1145. Key metabolic pathways involving shikimate and aromatic amino acid biosynthetic support were then engineered through a combination of gene deletions of competitive pathway steps, over-expressing feedback-insensitive enzymes in phenylalanine biosynthesis, and introducing a phosphoketolase-based pathway, which diverted glycolytic flux toward the formation of erythrose 4-phosphate (E4P). Through the stepwise engineering of A. nidulans A1145 outlined above, involving both heterologous pathway addition and native pathway metabolic engineering, we were able to produce cFF with titers reaching 611 mg/L in shake flask culture and 2.5 g/L in bench-scale fed-batch bioreactor culture. Our study establishes a production platform for cFF biosynthesis and successfully demonstrates engineering of phenylalanine derived diketopiperazines in a filamentous fungal host.

Functional annotation of Candida albicans hypothetical proteins: a bioinformatics approach

Candida albicans is a member of the ascomycetes class of fungi and it is an opportunistic pathogen species responsible for a wide range of fungal infections in humans. Bioinformatics and sequencing analysis of Candida proteomics has disclosed that around 69% proteome is still uncharacterized which needs to be annotated with functions. The NCBI-Genome has termed them as hypothetical proteins (HPs) in the whole proteome of Candida. Interpretation of this substantial portion of the proteome can reveal novel pharmacological targets for markers, drug development, and other therapeutics and so on. In this article, we have assigned functional annotation to these hypothetical proteins using bioinformatics methodologies. The advanced and robust computational models have been used to assign the preliminary functions to these putative HPs with high level of confidence. The findings of this study unveil some novel pharmacological targets for drug therapy and vaccines and it would help to identify novel molecular mechanisms underlying the fungal pathogenesis.

Pontiella agarivorans sp. nov., a novel marine anaerobic bacterium capable of degrading macroalgal polysaccharides and fixing nitrogen

Marine macroalgae produce abundant and diverse polysaccharides, which contribute substantially to the organic matter exported to the deep ocean. Microbial degradation of these polysaccharides plays an important role in the turnover of macroalgal biomass. Various members of the Planctomycetes-Verrucomicrobia-Chlamydia (PVC) superphylum are degraders of polysaccharides in widespread anoxic environments. In this study, we isolated a novel anaerobic bacterial strain NLcol2T from microbial mats on the surface of marine sediments offshore Santa Barbara, CA, USA. Based on 16S ribosomal RNA (rRNA) gene and phylogenomic analyses, strain NLcol2T represents a novel species within the Pontiella genus in the Kiritimatiellota phylum (within the PVC superphylum). Strain NLcol2T is able to utilize various monosaccharides, disaccharides, and macroalgal polysaccharides such as agar and ɩ-carrageenan. A near-complete genome also revealed an extensive metabolic capacity for anaerobic degradation of sulfated polysaccharides, as evidenced by 202 carbohydrate-active enzymes (CAZymes) and 165 sulfatases. Additionally, its ability of nitrogen fixation was confirmed by nitrogenase activity detected during growth on nitrogen-free medium, and the presence of nitrogenases (nifDKH) encoded in the genome. Based on the physiological and genomic analyses, this strain represents a new species of bacteria that may play an important role in the degradation of macroalgal polysaccharides and with relevance to the biogeochemical cycling of carbon, sulfur, and nitrogen in marine environments. Strain NLcol2T (= DSM 113125T = MCCC 1K08672T) is proposed to be the type strain of a novel species in the Pontiella genus, and the name Pontiella agarivorans sp. nov. is proposed.IMPORTANCEGrowth and intentional burial of marine macroalgae is being considered as a carbon dioxide reduction strategy but elicits concerns as to the fate and impacts of this macroalgal carbon in the ocean. Diverse heterotrophic microbial communities in the ocean specialize in these complex polymers such as carrageenan and fucoidan, for example, members of the Kiritimatiellota phylum. However, only four type strains within the phylum have been cultivated and characterized to date, and there is limited knowledge about the metabolic capabilities and functional roles of related organisms in the environment. The new isolate strain NLcol2T expands the known substrate range of this phylum and further reveals the ability to fix nitrogen during anaerobic growth on macroalgal polysaccharides, thereby informing the issue of macroalgal carbon disposal.

A bacterial sialidase mediates early-life colonization by a pioneering gut commensal

The early microbial colonization of the gastrointestinal tract can have long-term impacts on development and health. Keystone species, including Bacteroides spp., are prominent in early life and play crucial roles in maintaining the structure of the intestinal ecosystem. However, the process by which a resilient community is curated during early life remains inadequately understood. Here, we show that a single sialidase, NanH, in Bacteroides fragilis mediates stable occupancy of the intestinal mucosa in early life and regulates a commensal colonization program. This program is triggered by sialylated glycans, including those found in human milk oligosaccharides and intestinal mucus. NanH is required for vertical transmission from dams to pups and promotes B. fragilis dominance during early life. Furthermore, NanH facilitates commensal resilience and recovery after antibiotic treatment in a defined microbial community. Collectively, our study reveals a co-evolutionary mechanism between the host and microbiota mediated through host-derived glycans to promote stable colonization.

Streptococcus pneumoniae secretion chaperones PrsA, SlrA, and HtrA are required for competence, antibiotic resistance, colonization, and invasive disease

Streptococcus pneumoniae is a Gram-positive bacterium and a significant health threat with the populations most at risk being children, the elderly, and the immuno-compromised. To colonize and transition into an invasive infectious organism, S. pneumoniae secretes virulence factors that are translocated across the bacterial membrane and destined for surface exposure, attachment to the cell wall, or secretion into the host. The surface exposed protein chaperones PrsA, SlrA, and HtrA facilitate S. pneumoniae protein secretion; however, the distinct roles contributed by each of these secretion chaperones have not been well defined. Tandem Mass-Tagged Mass Spectrometry and virulence, adhesion, competence, and cell wall integrity assays were used to interrogate the individual and collective contributions of PrsA, SlrA, and HtrA to multiple aspects of S. pneumoniae physiology and virulence. PrsA, SlrA, and HtrA were found to play critical roles in S. pneumoniae host cell infection and competence, and the absence of each of these secretion chaperones significantly altered the S. pneumoniae secretome in distinct ways. PrsA and SlrA were additionally found to contribute to cell wall assembly and resistance to cell wall-active antimicrobials and were important for enabling S. pneumoniae host cell adhesion during colonization and invasive infection. These findings serve to further illustrate the pivotal contributions of PrsA, SlrA, and HtrA to S. pneumoniae protein secretion and virulence.

Charge cluster occurrence in land plants' mitochondrial proteomes with functional and structural insights

The Charge Clusters (CCs) are involved in key functions and are distributed according to the organism, the protein's type, and the charge of amino acids. In the present study, we have explored the occurrence, position, and annotation as a first large-scale study of the CCs in land plants mitochondrial proteomes. A new python script was used for data curation. The Finding Clusters Charge in Protein Sequences Program was performed after adjusting the reading window size. A 44316 protein sequences belonging to 52 species of land plants were analysed. The occurrence of Negative Charge Clusters (NCCs) (1.2%) is two times more frequent than the Positive Charge Clusters (PCCs) (0.64%). Moreover, 39 and 30 NCCs were conserved in 88 and 41 proteins in intra and in inter proteomes respectively, while 14 and 21 PCCs were conserved in 53 and 85 protein sequences in intra and inter proteomes consecutively. Sequences carrying mixed CCs are rare (0.12%). Despite this low abundance, CCs play a crucial role in protein function. The CCs tend to be located mainly in the terminal regions of proteins which guarantees specific protein targeting and import into the mitochondria. In addition, the functional annotation of CCs according to Gene Ontology shows that CCs are involved in binding functions of either proteins or macromolecules which are deployed in different metabolic and cellular processes such as RNA editing and transcription. This study may provide valuable information while considering the CCs in understanding the environmental adaptation of plants.Communicated by Ramaswamy H. Sarma.

Passionfruit Genomic Database (PGD): a comprehensive resource for passionfruit genomics

Passionfruit (Passiflora edulis) is a significant fruit crop in the commercial sector, owing to its high nutritional and medicinal value. The advent of high-throughput genomics sequencing technology has led to the publication of a vast amount of passionfruit omics data, encompassing complete genome sequences and transcriptome data under diverse stress conditions. To facilitate the efficient integration, storage, and analysis of these large-scale datasets, and to enable researchers to effectively utilize these omics data, we developed the first passionfruit genome database (PGD). The PGD platform comprises a diverse range of functional modules, including a genome browser, search function, heatmap, gene expression patterns, various tools, sequence alignment, and batch download, thereby providing a user-friendly interface. Additionally, supplementary practical tools have been developed for the PGD, such as gene family analysis tools, gene ontology (GO) terms, a pathway enrichment analysis, and other data analysis and mining tools, which enhance the data's utilization value. By leveraging the database's robust scalability, the intention is to continue to collect and integrate passionfruit omics data in the PGD, providing comprehensive and in-depth support for passionfruit research. The PGD is freely accessible via http://passionfruit.com.cn .