The EMBL-EBI bioinformatics web and programmatic tools framework

The Comparative Characterization of a Hypervirulent Acinetobacter baumannii Bacteremia Clinical Isolate Reveals a Novel Mechanism of Pathogenesis

Acinetobacter baumannii is an opportunistic Gram-negative pathogen with exquisite survival capabilities under various environmental conditions and displays widespread resistance to common antibiotics. A. baumannii is a leading cause of nosocomial infections that result in high morbidity and mortality rates. Accordingly, when multidrug resistance rates surpass threshold levels, the percentage of A. baumannii clinical isolates surges. Research into A. baumannii has increased in the past decade, and multiple mechanisms of pathogenesis have been identified, including mechanisms underlying biofilm development, quorum sensing, exotoxin production, secretion system utilization, and more. To date, the two gold-standard strains used to investigate different aspects of A. baumannii pathogenesis include ATCC 17978 and ATCC 19606. Here, we report a comparative characterization study of three additional A. baumannii clinical isolates obtained from different infection types and derived from different anatomical regions of infected patients. The comparison of three clinical isolates in addition to the ATCC strains revealed that the hypervirulent bacteremia clinical isolate, known as HUMC1, employs a completely different mechanism of pathogenesis when compared to all its counterparts. In stark contrast to the other genetic variants, the hypervirulent HUMC1 isolate does not form biofilms, is antibiotic-susceptible, and has the capacity to reach higher levels of quorum compared to the other clinically relevant strains. Our data also reveal that HUMC1 does not shed endotoxin into the extracellular milieu, rather secretes the evolutionarily conserved, host-mimicking, Zonula occludens toxin (Zot). Taken together, our hypothesis that HUMC1 cells have the ability to reach higher levels of quorum and lack biofilm production and endotoxin shedding, accompanied by the substantial elaboration of Zot, suggests a novel mechanism of pathogenesis that appears to afford the hypervirulent pathogen with stealth-like capabilities when disseminating through the circulatory system in a state of bacteremia.

A novel step towards the heterologous biosynthesis of paclitaxel: Characterization of T1βOH taxane hydroxylase

In the quest for innovative cancer therapeutics, paclitaxel remains a cornerstone in clinical oncology. However, its complex biosynthetic pathway, particularly the intricate oxygenation steps, has remained a puzzle in the decades following the characterization of the last taxane hydroxylase. The high divergence and promiscuity of enzymes involved have posed significant challenges. In this study, we adopted an innovative approach, combining in silico methods and functional gene analysis, to shed light on this elusive pathway. Our molecular docking investigations using a library of potential ligands uncovered TB574 as a potential missing enzyme in the paclitaxel biosynthetic pathway, demonstrating auspicious interactions. Complementary in vivo assays utilizing engineered S. cerevisiae strains as novel microbial cell factory consortia not only validated TB574's critical role in forging the elusive paclitaxel intermediate, T5αAc-1β,10β-diol, but also achieved the biosynthesis of paclitaxel precursors at an unprecedented yield including T5αAc-1β,10β-diol with approximately 40 mg/L. This achievement is highly promising, offering a new direction for further exploration of a novel metabolic engineering approaches using microbial consortia. In conclusion, our study not only furthers study the roles of previously uncharacterized enzymes in paclitaxel biosynthesis but also forges a path for pioneering advancements in the complete understanding of paclitaxel biosynthesis and its heterologous production. The characterization of T1βOH underscores a significant leap forward for future advancements in paclitaxel production using heterologous systems to improve cancer treatment and pharmaceutical production, thereby holding immense promise for enhancing the efficacy of cancer therapies and the efficiency of pharmaceutical manufacturing.

DEAD box RNA helicases are pervasive protein kinase interactors and activators

DEAD box (DDX) RNA helicases are a large family of ATPases, many of which have unknown functions. There is emerging evidence that besides their role in RNA biology, DDX proteins may stimulate protein kinases. To investigate if protein kinase-DDX interaction is a more widespread phenomenon, we conducted three orthogonal large-scale screens, including proteomics analysis with 32 RNA helicases, protein array profiling, and kinome-wide in vitro kinase assays. We retrieved Ser/Thr protein kinases as prominent interactors of RNA helicases and report hundreds of binary interactions. We identified members of ten protein kinase families, which bind to, and are stimulated by, DDX proteins, including CDK, CK1, CK2, DYRK, MARK, NEK, PRKC, SRPK, STE7/MAP2K, and STE20/PAK family members. We identified MARK1 in all screens and validated that DDX proteins accelerate the MARK1 catalytic rate. These findings indicate pervasive interactions between protein kinases and DEAD box RNA helicases, and provide a rich resource to explore their regulatory relationships.

Phylogenetic Tracing of Evolutionarily Conserved Zonula Occludens Toxin Reveals a "High Value" Vaccine Candidate Specific for Treating Multi-Strain Pseudomonas aeruginosa Infections

Extensively drug-resistant Pseudomonas aeruginosa infections are emerging as a significant threat associated with adverse patient outcomes. Due to this organism's inherent properties of developing antibiotic resistance, we sought to investigate alternative strategies such as identifying "high value" antigens for immunotherapy-based purposes. Through extensive database mining, we discovered that numerous Gram-negative bacterial (GNB) genomes, many of which are known multidrug-resistant (MDR) pathogens, including P. aeruginosa, horizontally acquired the evolutionarily conserved gene encoding Zonula occludens toxin (Zot) with a substantial degree of homology. The toxin's genomic footprint among so many different GNB stresses its evolutionary importance. By employing in silico techniques such as proteomic-based phylogenetic tracing, in conjunction with comparative structural modeling, we discovered a highly conserved intermembrane associated stretch of 70 amino acids shared among all the GNB strains analyzed. The characterization of our newly identified antigen reveals it to be a "high value" vaccine candidate specific for P. aeruginosa. This newly identified antigen harbors multiple non-overlapping B- and T-cell epitopes exhibiting very high binding affinities and can adopt identical tertiary structures among the least genetically homologous P. aeruginosa strains. Taken together, using proteomic-driven reverse vaccinology techniques, we identified multiple "high value" vaccine candidates capable of eliciting a polarized immune response against all the P. aeruginosa genetic variants tested.

Fungal Trunk Diseases Causing Decline of Apricot and Plum Trees in the Czech Republic

Fungal trunk diseases (FTDs) have been a significant threat to the global stone fruit industry. FTDs are caused by a consortium of wood-decaying fungi. These fungi colonize woody tissues, causing cankers, dieback, and other decline-related symptoms in host plants. In this study, a detailed screening of the fungal microbiota associated with the decline of stone fruit trees in the Czech Republic was performed. The wood fragments of plum and apricot trees showing symptoms of FTDs were subjected to fungal isolation. The partial internal transcribed spacer region, partial beta-tubulin, and translation elongation factor 1-α genes were amplified from genomic DNA extracted from fungal cultures. All isolates were classified, and the taxonomic placement of pathogenic strains was illustrated in phylogenetic trees. The most abundant pathogenic genus was Dactylonectria (31%), followed by Biscogniauxia (13%), Thelonectria (10%), Eutypa (9%), Dothiorella (7%), Diplodia (6%), and Diaporthe (6%). The most frequent endophytic genus was Aposphaeria (17%). The pathogenicity of six fungal species (Cadophora daguensis, Collophorina africana, Cytospora sorbicola, Dothiorella sarmentorum, Eutypa lata, and E. petrakii var. petrakii) to four Prunus spp. was evaluated, and Koch's postulates were fulfilled. All tested isolates caused lesions on at least one Prunus sp. The most aggressive species was E. lata, which caused the largest lesions on all four tested Prunus spp., followed by E. petrakii var. petrakii and D. sarmentorum. Japanese plum (Prunus salicina) and almond (P. amygdalus) were the most susceptible hosts, while apricot (P. armeniaca) was the least susceptible host in the pathogenicity trial.

Using EMBL-EBI Services via Web Interface and Programmatically via Web Services

The European Bioinformatics Institute (EMBL-EBI)'s Job Dispatcher framework provides access to a wide range of core databases and analysis tools that are of key importance in bioinformatics. As well as providing web interfaces to these resources, web services are available using REST and SOAP protocols that enable programmatic access and allow their integration into other applications and analytical workflows and pipelines. This article describes the various options available to researchers and bioinformaticians who would like to use our resources via the web interface employing RESTful web services clients provided in Perl, Python, and Java or who would like to use Docker containers to integrate the resources into analysis pipelines and workflows. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Retrieving data from EMBL-EBI using Dbfetch via the web interface Alternate Protocol 1: Retrieving data from EMBL-EBI using WSDbfetch via the REST interface Alternate Protocol 2: Retrieving data from EMBL-EBI using Dbfetch via RESTful web services with Python client Support Protocol 1: Installing Python REST web services clients Basic Protocol 2: Sequence similarity search using FASTA search via the web interface Alternate Protocol 3: Sequence similarity search using FASTA via RESTful web services with Perl client Support Protocol 2: Installing Perl REST web services clients Basic Protocol 3: Sequence similarity search using NCBI BLAST+ RESTful web services with Python client Basic Protocol 4: Sequence similarity search using HMMER3 phmmer REST web services with Perl client and Docker Support Protocol 3: Installing Docker and running the EMBL-EBI client container Basic Protocol 5: Protein functional analysis using InterProScan 5 RESTful web services with the Python client and Docker Alternate Protocol 4: Protein functional analysis using InterProScan 5 RESTful web services with the Java client Support Protocol 4: Installing Java web services clients Basic Protocol 6: Multiple sequence alignment using Clustal Omega via web interface Alternate Protocol 5: Multiple sequence alignment using Clustal Omega with Perl client and Docker Support Protocol 5: Exploring the RESTful API with OpenAPI User Inferface.

Odorant receptors for floral- and plant-derived volatiles in the yellow fever mosquito, Aedes aegypti (Diptera: Culicidae)

Adult mosquitoes require regular sugar meals, including nectar, to survive in natural habitats. Both males and females locate potential sugar sources using sensory proteins called odorant receptors (ORs) activated by plant volatiles to orient toward flowers or honeydew. The yellow fever mosquito, Aedes aegypti (Linnaeus, 1762), possesses a large gene family of ORs, many of which are likely to detect floral odors. In this study, we have uncovered ligand-receptor pairings for a suite of Aedes aegypti ORs using a panel of environmentally relevant, plant-derived volatile chemicals and a heterologous expression system. Our results support the hypothesis that these odors mediate sensory responses to floral odors in the mosquito's central nervous system, thereby influencing appetitive or aversive behaviors. Further, these ORs are well conserved in other mosquitoes, suggesting they function similarly in diverse species. This information can be used to assess mosquito foraging behavior and develop novel control strategies, especially those that incorporate mosquito bait-and-kill technologies.

Marine introgressions and Andean uplift have driven diversification in neotropical Monkey tree frogs (Anura, Phyllomedusinae)

The species richness in the Neotropics has been linked to environmental heterogeneity and a complex geological history. We evaluated which biogeographic processes were associated with the diversification of Monkey tree frogs, an endemic clade from the Neotropics. We tested two competing hypotheses: the diversification of Phyllomedusinae occurred either in a "south-north" or a "north-south" direction in the Neotropics. We also hypothesized that marine introgressions and Andean uplift had a crucial role in promoting their diversification. We used 13 molecular markers in a Bayesian analysis to infer phylogenetic relationships among 57 species of Phyllomedusinae and to estimate their divergence times. We estimated ancestral ranges based on 12 biogeographic units considering the landscape modifications of the Neotropical region. We found that the Phyllomedusinae hypothetical ancestor range was probably widespread throughout South America, from Western Amazon to Southern Atlantic Forest, at 29.5 Mya. The Phyllomedusines' ancestor must have initially diverged through vicariance, generally followed by jump-dispersals and sympatric speciation. Dispersal among areas occurred mostly from Western Amazonia towards Northern Andes and the South American diagonal of dry landscapes, a divergent pattern from both "south-north" and "north-south" diversification hypotheses. Our results revealed a complex diversification process of Monkey tree frogs, occurring simultaneously with the orogeny of Northern Andes and the South American marine introgressions in the last 30 million years.

Highly Resolved Genomes of Two Closely Related Lineages of the Rodent Louse Polyplax serrata with Different Host Specificities

Sucking lice of the parvorder Anoplura are permanent ectoparasites with specific lifestyle and highly derived features. Currently, genomic data are only available for a single species, the human louse Pediculus humanus. Here, we present genomes of two distinct lineages, with different host spectra, of a rodent louse Polyplax serrata. Genomes of these ecologically different lineages are closely similar in gene content and display a conserved order of genes, with the exception of a single translocation. Compared with P. humanus, the P. serrata genomes are noticeably larger (139 vs. 111 Mbp) and encode a higher number of genes. Similar to P. humanus, they are reduced in sensory-related categories such as vision and olfaction. Utilizing genome-wide data, we perform phylogenetic reconstruction and evolutionary dating of the P. serrata lineages. Obtained estimates reveal their relatively deep divergence (∼6.5 Mya), comparable with the split between the human and chimpanzee lice P. humanus and Pediculus schaeffi. This supports the view that the P. serrata lineages are likely to represent two cryptic species with different host spectra. Historical demographies show glaciation-related population size (Ne) reduction, but recent restoration of Ne was seen only in the less host-specific lineage. Together with the louse genomes, we analyze genomes of their bacterial symbiont Legionella polyplacis and evaluate their potential complementarity in synthesis of amino acids and B vitamins. We show that both systems, Polyplax/Legionella and Pediculus/Riesia, display almost identical patterns, with symbionts involved in synthesis of B vitamins but not amino acids.

Mechanistic and structural studies reveal NRAP-1-dependent coincident activation of NMDARs

N-methyl-D-aspartate (NMDA)-type ionotropic glutamate receptors have essential roles in neurotransmission and synaptic plasticity. Previously, we identified an evolutionarily conserved protein, NRAP-1, that is required for NMDA receptor (NMDAR) function in C. elegans. Here, we demonstrate that NRAP-1 was sufficient to gate NMDARs and greatly enhanced glutamate-mediated NMDAR gating, thus conferring coincident activation properties to the NMDAR. Intriguingly, vertebrate NMDARs-and chimeric NMDARs where the amino-terminal domain (ATD) of C. elegans NMDARs was replaced by the ATD from vertebrate receptors-were spontaneously active when ectopically expressed in C. elegans neurons. Thus, the ATD is a primary determinant of NRAP-1- and glutamate-mediated gating of NMDARs. We determined the crystal structure of NRAP-1 at 1.9-Å resolution, which revealed two distinct domains positioned around a central low-density lipoprotein receptor class A domain. The NRAP-1 structure, combined with chimeric and mutational analyses, suggests a model where the three NRAP-1 domains work cooperatively to modify the gating of NMDARs.

Whole-genome sequence and characterization of a marine red yeast, Rhodosporidium sphaerocarpum GDMCC 60679, featuring the assimilation of ammonia nitrogen

A marine red yeast, Rhodosporidium sphaerocarpum, is generally used for the production of lipids and carotenoids. In a previous study, we demonstrated that a marine-derived R. sphaerocarpum GDMCC 60679 can efficiently remove ammonia nitrogen and exhibit multiple probiotic functions for shrimp, Litopenaeus vannamei. Here, we performed a genome assembly of the strain GDMCC 60679 using a combination of the data from Illumina PE and PacBio CLR reads. The genome has a size of 18.03 Mb and consists of 32 contigs with an N50 length of 1,074,774 bp and GC content of 63 %. The genome was predicted to contain 6092 protein-coding genes, 5962 of which were functionally annotated. Metabolic pathways responsible for the ammonia assimilation and the synthesis of lipids and carotenoids were particularly examined to explore and characterize genes contributing to these functions. Whole-genome sequence and annotation of the strain lays a foundation to reveal the molecular mechanism of its prominent biological functions and will facilitate us to further expand new applications of yeasts in Rhodosporidium.

scInterpreter: a knowledge-regularized generative model for interpretably integrating scRNA-seq data

BACKGROUND

The rapid emergence of single-cell RNA-seq (scRNA-seq) data presents remarkable opportunities for broad investigations through integration analyses. However, most integration models are black boxes that lack interpretability or are hard to train.

RESULTS

To address the above issues, we propose scInterpreter, a deep learning-based interpretable model. scInterpreter substantially outperforms other state-of-the-art (SOTA) models in multiple benchmark datasets. In addition, scInterpreter is extensible and can integrate and annotate atlas scRNA-seq data. We evaluated the robustness of scInterpreter in a variety of situations. Through comparison experiments, we found that with a knowledge prior, the training process can be significantly accelerated. Finally, we conducted interpretability analysis for each dimension (pathway) of cell representation in the embedding space.

CONCLUSIONS

The results showed that the cell representations obtained by scInterpreter are full of biological significance. Through weight sorting, we found several new genes related to pathways in PBMC dataset. In general, scInterpreter is an effective and interpretable integration tool. It is expected that scInterpreter will bring great convenience to the study of single-cell transcriptomics.

Daisy: An integrated repeat protein curation service

Tandem repeats in proteins identification, classification and curation is a complex process that requires manual processing from experts, processing power and time. There are recent and relevant advances applying machine learning for protein structure prediction and repeat classification that are useful for this process. However, no service contemplates required databases and software to supplement researching on repeat proteins. In this publication we present Daisy, an integrated repeat protein curation web service. This service can process Protein Data Bank (PDB) and the AlphaFold Database entries for tandem repeats identification. In addition, it uses an algorithm to search a sequence against a library of Pfam hidden Markov model (HMM). Repeat classifications are associated with the identified families through RepeatsDB. This prediction is considered for enhancing the ReUPred algorithm execution and hastening the repeat units identification process. The service can also operate every associated PDB and AlphaFold structure with a UniProt proteome registry. Availability: The Daisy web service is freely accessible at daisy.bioinformatica.org.

Comparative analysis of cysteine proteases reveals gene family evolution of the group 1 allergens in astigmatic mites

BACKGROUND

Astigmatic mites contain potent allergens that can trigger IgE-mediated immune responses, leading to allergic diseases such as asthma, allergic rhinitis and atopic dermatitis. In house dust mites Dermatophagoides pteronyssinus and Dermatophagoides farinae, group 1 allergens (Der p 1 and Der f 1), characterized as papain-like cysteine proteases, have been defined as the major allergens that have high prevalence and potency. Previous studies of mite group 1 allergens mainly focused on identification, comparison of sequence and structure, as well as the investigation of cross-reactivity. To achieve a comprehensive view of mite group 1 allergens, we performed a comparative genomic analysis of all the cysteine proteases in six astigmatic mite species to elucidate the evolutionary relationships of group 1 allergens.

METHODS

Based on the high-quality and annotated genomes, all the cysteine proteases in six astigmatic mite species were identified by sequence homology search. The phylogenetic relationships, gene synteny and expression levels were revealed by bioinformatic tools. The allergenicity of recombinant cysteine proteases was evaluated by enzyme-linked immunosorbent assay.

RESULTS

Tandem duplication was revealed as the major feature of cysteine protease gene evolution in astigmatic mites. The high IgE-binding capacity and the significant expression level of the cysteine protease DP_007902.01 suggested its potential as a novel group 1 allergen of D. pteronyssinus. In addition, gene decay events were identified in the skin-burrowing parasitic mite Sarcoptes scabiei.

CONCLUSION

This comprehensive analysis provided insights into the evolution of cysteine proteases, as well as the component-resolved diagnosis of mite allergies.

Biological factors and statistical limitations prevent detection of most noncanonical proteins by mass spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry (MS) experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here, we leveraged recent advances in ribosome profiling and MS to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly expressed to be detected by shotgun MS at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for 4 noncanonical proteins in MS data, which were also supported by evolution and translation data. These results illustrate the power of MS to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly expressed proteins.

Expression of interferon regulatory factor family and its prognostic value in acute myeloid leukemia

Aim: To elucidate the clinicopathological and prognostic values of interferon regulatory factor (IRF) family genes in acute myeloid leukemia (AML). Patients & methods: Differential expression analysis and survival analysis from several reliable databases were conducted and further validated using patients with AML. Results: The expression level of IRF1/2/4/5/7/8/9 in patients with AML was upregulated, while IRF3/6 expression was downregulated. High IRF1/7/9 expression indicated a worse overall survival rate. Conclusion: Overexpression of IRF1/7/9 may be associated with poor survival in patients with AML, suggesting that the IRF family may be a promising therapeutic target.

An in silico analysis on the photoproteins Mnemiopsin 1 and Mnemiopsin 2 to explain the experimental results

Mnemiopsin 1 (Mn1) and Mnemiopsin 2 (Mn2) are photoproteins found in Mnemiopsis leidyi. We have tried to answer the question of whether the structural features of photoproteins can explain the observed activity data. According to the activity measurements data, they have the same characteristic wavelength. However, the initial intensity of Mn2 is significantly higher than that of Mn1, and decay time of Mn1 (0.92 s-1 ) is lower than that of Mn2 (1.46 s-1 ). The phylogenetic analysis demonstrates that, compared with Obelin and Aequorin from Obelia longissima and Aequorea victoria, respectively, a gene modification event may have caused the expansion of the N-terminal side of all photoproteins from M. leidyi. An in silico study has shown that the stability of the photoprotein-substrate complex of Mn2 is higher than that of Mn1, indicating a higher affinity of the substrate for Mn2 compared with Mn1. It was revealed that the active EF-hand loops 1 and III in Mn2 is locally more rigid compared with those in Mn1. We concluded that different stability of the photoprotein complexes leads to different initial intensity. While different patterns of the local dynamics of loops I and III may influence the decay rate.

Thimet oligopeptidase (THOP 1) distribution in cane toad (Bufo Marinus, Linnaeus, 1758) brain

Thimet oligopeptides (THOP 1) is a metal-dependent peptidase involved in the metabolism of neuropeptides and the presentation of peptides via MHC-1. It has been shown to play a role in the regulation of protein-protein interactions and the metabolism of intracellular peptides. THOP 1 is associated with important biological processes such as metabolism and neurodegenerative diseases, among others. The objective of this study is to elucidate the distribution of THOP 1 in the Bufo marinus brain. The analysis of THOP 1 amino acid sequences indicates that they have been conserved throughout evolution, with significant homology observed across various phyla. When comparing amphibians with other species, more than 70% identity can be identified. Immunohistochemistry analysis of the toad's brain has demonstrated that the enzyme has a ubiquitous distribution, consistent with previous findings in mammals. THOP 1 can be found in important areas of the brain, such as bulb, thalamic nuclei, striatum, hypothalamus, and among others. Nonetheless, THOP 1 is consistently localized within the nucleus, a pattern also observed in the rat brain. Therefore, based on these results, the toad appears to be an excellent model for studying the general biology of THOP 1, given the substantial homology of this enzyme with mammals and its similarity in distribution within the brain.

Comparative genomics of the medicinal plants Lonicera macranthoides and L. japonica provides insight into genus genome evolution and hederagenin-based saponin biosynthesis

Lonicera macranthoides (LM) and L. japonica (LJ) are medicinal plants widely used in treating viral diseases, such as COVID-19. Although the two species are morphologically similar, their secondary metabolite profiles are significantly different. Here, metabolomics analysis showed that LM contained ~86.01 mg/g hederagenin-based saponins, 2000-fold higher than LJ. To gain molecular insights into its secondary metabolite production, a chromosome-level genome of LM was constructed, comprising 9 pseudo-chromosomes with 40 097 protein-encoding genes. Genome evolution analysis showed that LM and LJ were diverged 1.30-2.27 million years ago (MYA). The two plant species experienced a common whole-genome duplication event that occurred ∼53.9-55.2 MYA before speciation. Genes involved in hederagenin-based saponin biosynthesis were arranged in clusters on the chromosomes of LM and they were more highly expressed in LM than in LJ. Among them, oleanolic acid synthase (OAS) and UDP-glycosyltransferase 73 (UGT73) families were much more highly expressed in LM than in LJ. Specifically, LmOAS1 was identified to effectively catalyse the C-28 oxidation of β-Amyrin to form oleanolic acid, the precursor of hederagenin-based saponin. LmUGT73P1 was identified to catalyse cauloside A to produce α-hederin. We further identified the key amino acid residues of LmOAS1 and LmUGT73P1 for their enzymatic activities. Additionally, comparing with collinear genes in LJ, LmOAS1 and LmUGT73P1 had an interesting phenomenon of 'neighbourhood replication' in LM genome. Collectively, the genomic resource and candidate genes reported here set the foundation to fully reveal the genome evolution of the Lonicera genus and hederagenin-based saponin biosynthetic pathway.

Biological Factors and Statistical Limitations Prevent Detection of Most Noncanonical Proteins by Mass Spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here we leveraged recent advances in ribosome profiling and mass spectrometry to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly-expressed to be detected by shotgun mass spectrometry at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for four noncanonical proteins in mass spectrometry data, which were also supported by evolution and translation data. These results illustrate the power of mass spectrometry to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly-expressed proteins.

Genome-wide analysis of the U-box E3 ubiquitin ligase family role in drought tolerance in sesame (Sesamum indicum L.)

Plant U-box (PUB) proteins belong to a class of ubiquitin ligases essential in various biological processes. Sesame (Sesamum indicum L.) is an important and worldwide cultivated oilseed crop. However few studies have been conducted to explore the role of PUBs in drought tolerance in sesame. This study identified a total of 56 members of the sesame PUB family (SiPUB) genes distributed unevenly across all 13 chromosomes. Based on phylogenetic analysis, all 56 SiPUB genes were classified into six groups with various structures and motifs. Cis-acting element analysis suggested that the SiPUB genes are involved in response to various stresses including drought. Based on RNA-seq analysis and quantitative real-time PCR, we identified nine SiPUB genes with significantly different expression profiles under drought stress. The expression patterns of six SiPUB genes in root, leaf and stem tissues corroborated the reliability of the RNA-seq datasets. These findings underscore the importance of SiPUB genes in enhancing drought tolerance in sesame plants. Our study provides novel insights into the evolutionary patterns and variations of PUB genes in sesame and lays the foundation for comprehending the functional characteristics of SiPUB genes under drought-induced stress conditions.

Two new species of Trichocomaceae (Eurotiales), accommodated in Rasamsonia and Talaromyces section Bacillispori, from the Czech Republic

During a previous study on microfungi associated with clematis roots, Penicillium-like fungi were isolated and identified based on morphology. In this study, we subjected those strains to a detailed examination which led to the proposal of two taxonomic novelties, named Rasamsonia chlamydospora and Talaromyces clematidis. The first taxon is characterized by rough-walled mycelium, acerose to flask shaped phialides, cylindrical conidia and by production of chlamydospore-like structures. The four-loci-based phylogeny analysis delineated the taxon as a taxonomic novelty in Rasamsonia. Talaromyces clematidis is characterized by restricted growth on Czapek yeast extract agar, dichloran 18% glycerol agar and yeast extract sucrose agar, and production of yellow ascomata on oatmeal agar. Phylogenetic analyses placed this taxon as a taxonomic novelty in Talaromyces sect. Bacillispori. Both taxa are introduced here with detailed descriptions, photoplates and information on their phylogenetic relationship with related species.

The Arf-GEF GBF1 undergoes multi-domain structural shifts to activate Arf at the Golgi

Golgi homeostasis require the activation of Arf GTPases by the guanine-nucleotide exchange factor requires GBF1, whose recruitment to the Golgi represents a rate limiting step in the process. GBF1 contains a conserved, catalytic, Sec7 domain (Sec7d) and five additional (DCB, HUS, HDS1-3) domains. Herein, we identify the HDS3 domain as essential for GBF1 membrane association in mammalian cells and document the critical role of HDS3 during the development of Drosophila melanogaster. We show that upon binding to Golgi membranes, GBF1 undergoes conformational changes in regions bracketing the catalytic Sec7d. We illuminate GBF1 interdomain arrangements by negative staining electron microscopy of full-length human GBF1 to show that GBF1 forms an anti-parallel dimer held together by the paired central DCB-HUS core, with two sets of HDS1-3 arms extending outward in opposite directions. The catalytic Sec7d protrudes from the central core as a largely independent domain, but is closely opposed to a previously unassigned α-helix from the HDS1 domain. Based on our data, we propose models of GBF1 engagement on the membrane to provide a paradigm for understanding GBF1-mediated Arf activation required for cellular and organismal function.

Exploring Proteus mirabilis Methionine tRNA Synthetase Active Site: Homology Model Construction, Molecular Dynamics, Pharmacophore and Docking Validation

Currently, the treatment of Proteus mirabilis infections is considered to be complicated as the organism has become resistant to numerous antibiotic classes. Therefore, new inhibitors should be developed, targeting bacterial molecular functions. Methionine tRNA synthetase (MetRS), a member of the aminoacyl-tRNA synthetase family, is essential for protein biosynthesis offering a promising target for novel antibiotics discovery. In the context of computer-aided drug design (CADD), the current research presents the construction and analysis of a comparative homology model for P. mirabilis MetRS, enabling development of novel inhibitors with greater selectivity. Molecular Operating Environment (MOE) software was used to build a homology model for P. mirabilis MetRS using Escherichia coli MetRS as a template. The model was evaluated, and the active site of the target protein predicted from its sequence using conservation analysis. Molecular dynamic simulations were performed to evaluate the stability of the modeled protein structure. In order to evaluate the predicted active site interactions, methionine (the natural substrate of MetRS) and several inhibitors of bacterial MetRS were docked into the constructed model using MOE. After validation of the model, pharmacophore-based virtual screening for a systemically prepared dataset of compounds was performed to prove the feasibility of the proposed model, identifying possible parent compounds for further development of MetRS inhibitors against P. mirabilis.

The first multi-tissue genome-scale metabolic model of a woody plant highlights suberin biosynthesis pathways in Quercus suber

Over the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behaviour at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871). The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyse the pathways associated with the synthesis of suberin monomers, namely the acyl-lipids, phenylpropanoids, isoprenoids, and flavonoids production. The models developed in this work provide a systematic overview of the metabolism of Q. suber, including its secondary metabolism pathways and cork formation.

Chromosomal level genome assembly of medicinal plant Sophora flavescens

Sophora flavescens is a medicinal plant in the genus Sophora of the Fabaceae family. The root of S. flavescens is known in China as Kushen and has a long history of wide use in multiple formulations of Traditional Chinese Medicine (TCM). In this study, we used third-generation Nanopore long-read sequencing technology combined with Hi-C scaffolding technology to de novo assemble the S. flavescens genome. We obtained a chromosomal level high-quality S. flavescens draft genome. The draft genome size is approximately 2.08 Gb, with more than 80% annotated as Transposable Elements (TEs), which have recently and rapidly proliferated. This genome size is ~5x larger than its closest sequenced relative Lupinus albus L. . We annotated 60,485 genes and examined their expression profiles in leaf, stem and root tissues, and also characterised the genes and pathways involved in the biosynthesis of major bioactive compounds, including alkaloids, flavonoids and isoflavonoids. The assembled genome highlights the very different evolutionary trajectories that have occurred in recently diverged Fabaceae, leading to smaller duplicated genomes.

Detection and Quantification of Botanical Impurities in Commercial Oregano (Origanum vulgare) Using Metabarcoding and Digital PCR

DNA technology for food authentication is already well established, and with the advent of Next Generation Sequencing (NGS) and, more specifically, metabarcoding, compositional analysis of food at the molecular level has rapidly gained popularity. This has led to several reports in the media about the presence of foreign, non-declared species in several food commodities. As herbs and spices are attractive targets for fraudulent manipulation, a combination of digital PCR and metabarcoding by NGS was employed to check the purity of 285 oregano samples taken from the European market. By using novel primers and analytical approaches, it was possible to detect and quantify both adulterants and contaminants in these samples. The results highlight the high potential of NGS for compositional analysis, although its quantitative information (read count percentages) is unreliable, and other techniques are therefore needed to complement the sequencing information for assessing authenticity ('true to the name') of food ingredients.

Analysis of macerated ticks using Boolean logic gating colorimetric isothermal nucleic acid assays for Lyme Borrelia and Ixodes scapularis ticks

Lyme disease, one of the most common tickborne diseases, has been rapidly spreading in parallel with the expansion of the range of its tick vector. Better tick surveillance efforts are needed to accurately estimate disease risk and to guide public health and clinical management. We have developed two multiplex loop-mediated isothermal amplification (LAMP) reactions coupled with oligonucleotide strand displacement (OSD) probes to identify the tick host, Ixodes scapularis, and the Lyme disease pathogen, Borrelia burgdorferi, they carry. In each multiplex LAMP-OSD assay the co-presence of two target sequences is computed at the DNA level by linking the two corresponding amplicons and detecting the co-product on colorimetric lateral flow dipsticks. In tests with synthetic DNA, the co-presence of as few as four copies of input DNA could be detected, without producing spurious signals. Most importantly, though, the LAMP-OSD assay is amenable to being carried out directly with macerated tick samples, without any sample preparation. In such field conditions, assays performed robustly and demonstrated 97-100% sensitivity and 100% specificity with both field-collected and lab-raised artificially infected ticks. Such easy-to-use, arthropod and pathogen-specific assays would be well suited to field and near patient use without relying on complex instrumentation or infrastructure.

In Silico Analysis of Honey Bee Peptides as Potential Inhibitors of Capripoxvirus DNA-Directed RNA Polymerase

The genus Capripoxvirus belongs to the Poxviridae family. The sheeppox, goatpox, and lumpy skin disease viruses are three species of this genus with 96% identity in their genomes. These are financially devastating viral infections among cattle, which cause a reduction in animal products and lead to a loss in livestock industries. In the current study, the phylogenetic analysis was carried out to reveal the evolutionary relationships of Capripoxvirus species (i.e., sheeppox virus (SPPV), goatpox virus (GTPV), and lumpy skin disease virus (LSDV)) with other viruses from the Poxviridae family with >96% query coverage to find the similarity index among all members. The three viruses (i.e., SPPV, GTPV, and LSDV) joined the clade of Capripoxvirus of the Poxviridae family in the phylogenetic tree and exhibited close evolutionary relationships. The multiple sequence alignment using ClustalOmega revealed significant variations in the protein sequences of the DNA-dependent RNA polymerase of SPPV, GTPV, and LSDV. The three-dimensional structures of five selected bee peptides and DNA-directed RNA polymerase of SPPV, GTPV, and LSDV were predicted using trRosetta and I-TASSER and used for molecular docking and simulation studies. The protein-protein docking was carried out using HADDOCK server to explore the antiviral activity of peptides as honey bee proteins against SPPV, GTPV, and LSDV. In total, five peptides were docked to DNA-directed RNA polymerase of these viruses. The peptides mellitin and secapin-1 displayed the lowest binding scores (-106.9 +/- 7.2 kcal/mol and -101.4 +/- 11.3 kcal/mol, respectively) and the best patterns with stable complexes. The molecular dynamics simulation indicated that the complex of protein DNA-dependent RNA polymerase and the peptide melittin stayed firmly connected and the peptide binding to the receptor protein was stable. The findings of this study provide the evidence of bee peptides as potent antimicrobial agents against sheeppox, goatpox, and lumpy skin disease viruses with no complexity.

An integrative bioinformatics investigation and experimental validation of chromobox family in diffuse large B-cell lymphoma

BACKGROUND

Diffuse large B-cell lymphoma (DLBCL) is one of the most aggressive malignant tumors. Chromobox (CBX) family plays the role of oncogenes in various malignancies.

METHODS

The transcriptional and protein levels of CBX family were confirmed by GEPIA, Oncomine, CCLE, and HPA database. Screening of co-expressed genes and gene function enrichment analysis were performed by GeneMANIA and DAVID 6.8. The prognostic value, immune cell infiltration and drug sensitivity analysis of CBX family in DLBCL were performed by Genomicscape, TIMER2.0, and GSCALite database. Confirmatory Tests of CBX family protein expression in DLBCL were performed by immunohistochemistry.

RESULTS

The mRNA and protein expressions of CBX1/2/3/5/6 were higher in DLBCL tissues than control groups. Enrichment analysis showed that the functions of CBX family were mainly related to chromatin remodeling, methylation-dependent protein binding, and VEGF signaling pathway. The high mRNA expressions of CBX2/3/5/6 were identified to be associated with short overall survival (OS) in DLBCL patients. Multivariate COX regression indicated that CBX3 was independent prognostic marker. Immune infiltration analysis revealed that the mRNA expressions of CBX family (especially CBX1, CBX5, and CBX6) in DLBCL were significantly correlated with the infiltration of most immune cells (including B cells, CD8 + T cells, CD4 + T cells, neutrophils, monocytes, macrophages, and Treg cells). Meanwhile, there was a strong correlation between the expression levels of CBX1/5/6 and surface markers of immune cells, such as the widely studied PVR-like protein receptor/ligand and PDL-1 immune checkpoint. Notably, our study found that DLBCL cells with CBX1 over-expression were resistant to the common anti-tumor drugs, but CBX2/5 had two polarities. Finally, we confirmed the higher expressions of CBX1/2/3/5/6 in DLBCL tissues compared with control groups by immunohistochemistry.

CONCLUSION

We provided a detailed analysis of the relationship between the CBX family and the prognosis of DLBCL. Distinguished from other studies, We found that high mRNA expressions of CBX2/3/5/6 were associated with poor prognosis in DLBCL patients, and Multivariate COX regression indicated that CBX3 was independent prognostic marker. Besides, our study also found an association between the CBX family and anti-tumour drug resistance, and provided a relationship between CBX family expression and immune cell infiltration.

Structural and functional analysis of the cerato-platanin-like protein Cpl1 suggests diverging functions in smut fungi

Plant-pathogenic fungi are causative agents of the majority of plant diseases and can lead to severe crop loss in infected populations. Fungal colonization is achieved by combining different strategies, such as avoiding and counteracting the plant immune system and manipulating the host metabolome. Of major importance are virulence factors secreted by fungi, which fulfil diverse functions to support the infection process. Most of these proteins are highly specialized, with structural and biochemical information often absent. Here, we present the atomic structures of the cerato-platanin-like protein Cpl1 from Ustilago maydis and its homologue Uvi2 from Ustilago hordei. Both proteins adopt a double-Ψβ-barrel architecture reminiscent of cerato-platanin proteins, a class so far not described in smut fungi. Our structure-function analysis shows that Cpl1 binds to soluble chitin fragments via two extended grooves at the dimer interface of the two monomer molecules. This carbohydrate-binding mode has not been observed previously and expands the repertoire of chitin-binding proteins. Cpl1 localizes to the cell wall of U. maydis and might synergize with cell wall-degrading and decorating proteins during maize infection. The architecture of Cpl1 harbouring four surface-exposed loop regions supports the idea that it might play a role in the spatial coordination of these proteins. While deletion of cpl1 has only mild effects on the virulence of U. maydis, a recent study showed that deletion of uvi2 strongly impairs U. hordei virulence. Our structural comparison between Cpl1 and Uvi2 reveals sequence variations in the loop regions that might explain a diverging function.

During heat stress in Myxococcus xanthus, the CdbS PilZ domain protein, in concert with two PilZ-DnaK chaperones, perturbs chromosome organization and accelerates cell death

C-di-GMP is a bacterial second messenger that regulates diverse processes in response to environmental or cellular cues. The nucleoid-associated protein (NAP) CdbA in Myxococcus xanthus binds c-di-GMP and DNA in a mutually exclusive manner in vitro. CdbA is essential for viability, and CdbA depletion causes defects in chromosome organization, leading to a block in cell division and, ultimately, cell death. Most NAPs are not essential; therefore, to explore the paradoxical cdbA essentiality, we isolated suppressor mutations that restored cell viability without CdbA. Most mutations mapped to cdbS, which encodes a stand-alone c-di-GMP binding PilZ domain protein, and caused loss-of-function of cdbS. Cells lacking CdbA and CdbS or only CdbS were fully viable and had no defects in chromosome organization. CdbA depletion caused post-transcriptional upregulation of CdbS accumulation, and this CdbS over-accumulation was sufficient to disrupt chromosome organization and cause cell death. CdbA depletion also caused increased accumulation of CsdK1 and CsdK2, two unusual PilZ-DnaK chaperones. During CdbA depletion, CsdK1 and CsdK2, in turn, enabled the increased accumulation and toxicity of CdbS, likely by stabilizing CdbS. Moreover, we demonstrate that heat stress, possibly involving an increased cellular c-di-GMP concentration, induced the CdbA/CsdK1/CsdK2/CdbS system, causing a CsdK1- and CsdK2-dependent increase in CdbS accumulation. Thereby this system accelerates heat stress-induced chromosome mis-organization and cell death. Collectively, this work describes a unique system that contributes to regulated cell death in M. xanthus and suggests a link between c-di-GMP signaling and regulated cell death in bacteria.

Attenuation of Dopaminergic Neurodegeneration in a C. elegans Parkinson's Model through Regulation of Xanthine Dehydrogenase (XDH-1) Expression by the RNA Editase, ADR-2

Differential RNA editing by adenosine deaminases that act on RNA (ADARs) has been implicated in several neurological disorders, including Parkinson's disease (PD). Here, we report results of a RNAi screen of genes differentially regulated in adr-2 mutants, normally encoding the only catalytically active ADAR in Caenorhabditis elegans, ADR-2. Subsequent analysis of candidate genes that alter the misfolding of human α-synuclein (α-syn) and dopaminergic neurodegeneration, two PD pathologies, reveal that reduced expression of xdh-1, the ortholog of human xanthine dehydrogenase (XDH), is protective against α-synuclein-induced dopaminergic neurodegeneration. Further, RNAi experiments show that WHT-2, the worm ortholog of the human ABCG2 transporter and a predicted interactor of XDH-1, is the rate-limiting factor in the ADR-2, XDH-1, WHT-2 system for dopaminergic neuroprotection. In silico structural modeling of WHT-2 indicates that the editing of one nucleotide in the wht-2 mRNA leads to the substitution of threonine with alanine at residue 124 in the WHT-2 protein, changing hydrogen bonds in this region. Thus, we propose a model where wht-2 is edited by ADR-2, which promotes optimal export of uric acid, a known substrate of WHT-2 and a product of XDH-1 activity. In the absence of editing, uric acid export is limited, provoking a reduction in xdh-1 transcription to limit uric acid production and maintain cellular homeostasis. As a result, elevation of uric acid is protective against dopaminergic neuronal cell death. In turn, increased levels of uric acid are associated with a decrease in ROS production. Further, downregulation of xdh-1 is protective against PD pathologies because decreased levels of XDH-1 correlate to a concomitant reduction in xanthine oxidase (XO), the form of the protein whose by-product is superoxide anion. These data indicate that modifying specific targets of RNA editing may represent a promising therapeutic strategy for PD.

The predictive value of prognosis and therapeutic response for STAT family in pancreatic cancer

Background

Signal transducers and activators of transcription (STAT) proteins, well-known cytoplasmic transcription factors, were found to be abnormally expressed in various cancers and play essential parts in the initiation, progression and therapy resistance of cancer. Nevertheless, the functions of different STATs in pancreatic cancer (PC) and their relationship to the prognosis and immune infiltration as well as drug efficacy in PC patients have not been systematically elucidated.

Methods

Expression, prognosis, genetic alterations and pathway enrichment analyses of the STAT family were investigated via Oncomine, GEPIA, Kaplan Meier-plotter, cBioPortal, Metascape and GSEA. Analysis of tumor immune microenvironment was conducted by ESTIMATE and TIMER. "pRRophetic" packages were used for analysis of chemotherapeutic response. Finally, the diagnostic and prognostic value of key STATs were further validated through public datasets and immunohistochemistry.

Results

In this study, only STAT1 mRNA level was significantly increased in tumor tissues and highly expressed in PC cell lines via multiple datasets. PC patients with higher STAT1/4/6 expression had a worse overall survival (OS) and progression-free survival (PFS), while higher STAT5B expression was correlated with better prognosis in the TCGA cohort. The STATs-associated genes were enriched in pathways about the remodeling of tumor immune microenvironment. The STATs levels were significantly correlated with immune infiltration, except STAT6. The STAT1 was identified as a potential biomarker and its diagnostic and prognostic value were further validated at mRNA and protein levels. GSEA showed that STAT1 may be involved in the progression and immune regulations of PC. Moreover, STAT1 expression was significantly related to the level of immune checkpoint, and predicted immunotherapy and chemotherapy responses.

Conclusion

STAT family members were comprehensively analyzed and STAT1 was identified as an effective biomarker for predicting the survival and therapeutic response, which could be beneficial to develop better treatment strategies.

Members of the tomato NRC4 h-NLR family augment each other in promoting basal immunity

Plants possess an efficient, two-tiered immune system to combat pathogens and pests. Several decades of research have characterized different features of these two well-known tiers, PTI and ETI (Pattern/ Effector-triggered Immunity). NLR (Nucleotide-binding domain Leucine-rich Repeat) receptors have been found to link PTI to ETI, and be required for full potentiation of plant immune responses in several systems. Intra-cellular helper-NLRs (h-NLRs) mediate ETI and have been focused on extensively in recent research. Previously, we investigated the roles of the h-NLR SlNRC4a in tomato immunity, finding that a specific mutation in this gene results in gain of function constitutive defense activation and broad disease resistance. Deletion of the entire NRC4 clade, which contains 3 genes, can compromise tomato immunity. Here, we decided to investigate the role of an additional clade member, SlNRC4b, in basal immunity. We generated a gain of function mutant in SlNRC4b using CRISPR-Cas9, as well as a double gain of function mutant in both genes. Similarly to the slnrc4a mutant, a slnrc4b mutant also possessed increased basal immunity and broad spectrum disease resistance. The double mutant displayed additive effects in some cases, with significant increases in resistance to fungal phytopathogens as compared with each of the single mutants. Our work confirms that the NRC4 family h-NLRs are important in the plant immune system, suggesting that this gene family has the potential to be promising in targeted agricultural adaptation in the Solanaceae family, promoting disease resistance and prevention of yield loss to pathogens.

Comparative genomics analyses reveal sequence determinants underlying interspecies variations in injury-responsive enhancers

BACKGROUND

Injury induces profound transcriptional remodeling events, which could lead to only wound healing, partial tissue repair, or perfect regeneration in different species. Injury-responsive enhancers (IREs) are cis-regulatory elements activated in response to injury signals, and have been demonstrated to promote tissue regeneration in some organisms such as zebrafish and flies. However, the functional significances of IREs in mammals remain elusive. Moreover, whether the transcriptional responses elicited by IREs upon injury are conserved or specialized in different species, and what sequence features may underlie the functional variations of IREs have not been elucidated.

RESULTS

We identified a set of IREs that are activated in both regenerative and non-regenerative neonatal mouse hearts upon myocardial ischemia-induced damage by integrative epigenomic and transcriptomic analyses. Motif enrichment analysis showed that AP-1 and ETS transcription factor binding motifs are significantly enriched in both zebrafish and mouse IREs. However, the IRE-associated genes vary considerably between the two species. We further found that the IRE-related sequences in zebrafish and mice diverge greatly, with the loss of IRE inducibility accompanied by a reduction in AP-1 and ETS motif frequencies. The functional turnover of IREs between zebrafish and mice is correlated with changes in transcriptional responses of the IRE-associated genes upon injury. Using mouse cardiomyocytes as a model, we demonstrated that the reduction in AP-1 or ETS motif frequency attenuates the activation of IREs in response to hypoxia-induced damage.

CONCLUSIONS

By performing comparative genomics analyses on IREs, we demonstrated that inter-species variations in AP-1 and ETS motifs may play an important role in defining the functions of enhancers during injury response. Our findings provide important insights for understanding the molecular mechanisms of transcriptional remodeling in response to injury across species.

CRISPR-Cas9-mediated knockout of CYP79D1 and CYP79D2 in cassava attenuates toxic cyanogen production

Cassava (Manihot esculenta) is a starchy root crop that supports over a billion people in tropical and subtropical regions of the world. This staple, however, produces the neurotoxin cyanide and requires processing for safe consumption. Excessive consumption of insufficiently processed cassava, in combination with protein-poor diets, can have neurodegenerative impacts. This problem is further exacerbated by drought conditions which increase this toxin in the plant. To reduce cyanide levels in cassava, we used CRISPR-mediated mutagenesis to disrupt the cytochrome P450 genes CYP79D1 and CYP79D2 whose protein products catalyze the first step in cyanogenic glucoside biosynthesis. Knockout of both genes eliminated cyanide in leaves and storage roots of cassava accession 60444; the West African, farmer-preferred cultivar TME 419; and the improved variety TMS 91/02324. Although knockout of CYP79D2 alone resulted in significant reduction of cyanide, mutagenesis of CYP79D1 did not, indicating these paralogs have diverged in their function. The congruence of results across accessions indicates that our approach could readily be extended to other preferred or improved cultivars. This work demonstrates cassava genome editing for enhanced food safety and reduced processing burden, against the backdrop of a changing climate.

The survival of Amblyomma sculptum ticks upon blood-feeding depends on the expression of an inhibitor of apoptosis protein

BACKGROUND

The tick Amblyomma sculptum is the major vector of Rickettsia rickettsii, the causative agent of the highly lethal Brazilian spotted fever. It has been shown that R. rickettsii inhibits apoptosis in both human endothelial cells and tick cells. Apoptosis is regulated by different factors, among which inhibitors of apoptosis proteins (IAPs) play a central role. In the study reported here, we selected an IAP of A. sculptum that has not yet been characterized to assess its role in cell death and to determine the effects of its gene silencing on tick fitness and R. rickettsii infection.

METHODS

An A. sculptum cell line (IBU/ASE-16) was treated with specific double-stranded RNA (dsRNA) for either IAP (dsIAP) or green fluorescent protein (dsGFP; as a control). The activity of caspase-3 and the exposure of phosphatidylserine were determined in both groups. In addition, unfed adult ticks, infected or not infected with R. rickettsii, were treated with either dsIAP or dsGFP and allowed to feed on noninfected rabbits. In parallel, noninfected ticks were allowed to feed on an R. rickettsii-infected rabbit. Ticks (infected or not with R. rickettsii) that remained unfed were used as a control.

RESULTS

Caspase-3 activity and the externalization of phosphatidylserine were significantly higher in IBU/ASE-16 cells treated with dsIAP than in those treated with dsGFP. The mortality rates of ticks in the dsIAP group were much higher than those in the dsGFP group when they were allowed to feed on rabbits, independent of the presence of R. rickettsii. Conversely, lower mortality rates were recorded in unfed ticks.

CONCLUSIONS

Our results show that IAP negatively regulates apoptosis in A. sculptum cells. Moreover, IAP-silenced ticks experienced higher mortality rates following the acquisition of a blood meal, suggesting that feeding may trigger the activation of apoptosis in the absence of this physiological regulator. These findings indicate that IAP is a potential antigen for an anti-tick vaccine.

Interaction between the α-glucosidases, sucrase-isomaltase and maltase-glucoamylase, in human intestinal brush border membranes and its potential impact on disaccharide digestion

The two major intestinal α-glycosidases, sucrase-isomaltase (SI) and maltase-glucoamylase (MGAM), are active towards α-1,4 glycosidic linkages that prevail in starch. These enzymes share striking structural similarities and follow similar biosynthetic pathways. It has been hypothesized that starch digestion can be modulated via “toggling” of activities of these mucosal α-glycosidases, suggesting a possible interaction between these two enzyme complexes in the intestinal brush border membrane (BBM). Here, the potential interaction between SI and MGAM was investigated in solubilized BBMs utilizing reciprocal pull down assays, i.e., immunoprecipitation with anti-SI antibody followed by Western blotting with anti-MGAM antibody and vice versa. Our results demonstrate that SI interacts avidly with MGAM concomitant with a hetero-complex assembly in the BBMs. This interaction is resistant to detergents, such as Triton X-100 or Triton X-100 in combination with sodium deoxycholate. By contrast, inclusion of sodium deoxycholate into the solubilization buffer reduces the enzymatic activities towards sucrose and maltose substantially, most likely due to alterations in the quaternary structure of either enzyme. In view of their interaction, SI and MGAM regulate the final steps in starch digestion in the intestine, whereby SI assumes the major role by virtue of its predominant expression in the intestinal BBMs, while MGAM acts in auxiliary supportive fashion. These findings will help understand the pathophysiology of carbohydrate malabsorption in functional gastrointestinal disorders, particularly in irritable bowel syndrome, in which gene variants of SI are implicated.

Primate protein-ligand interfaces exhibit significant conservation and unveil human-specific evolutionary drivers

Despite the vast phenotypic differences observed across primates, their protein products are largely similar to each other at the sequence level. We hypothesized that, since proteins accomplish all their functions via interactions with other molecules, alterations in the sites that participate in these interactions may be of critical importance. To uncover the extent to which these sites evolve across primates, we built a structurally-derived dataset of ~4,200 one-to-one orthologous sequence groups across 18 primate species, consisting of ~68,000 ligand-binding sites that interact with DNA, RNA, small molecules, ions, or peptides. Using this dataset, we identify functionally important patterns of conservation and variation within the amino acid residues that facilitate protein-ligand interactions across the primate phylogeny. We uncover that interaction sites are significantly more conserved than other sites, and that sites binding DNA and RNA further exhibit the lowest levels of variation. We also show that the subset of ligand-binding sites that do vary are enriched in components of gene regulatory pathways and uncover several instances of human-specific ligand-binding site changes within transcription factors. Altogether, our results suggest that ligand-binding sites have experienced selective pressure in primates and propose that variation in these sites may have an outsized effect on phenotypic variation in primates through pleiotropic effects on gene regulation.

On the design of linked datasets mapping networks of collaboration in the genomic sequencing of Saccharomyces cerevisiae, Homo sapiens, and Sus scrofa

This data note describes a unique two-step methodology to construct six linked datasets covering the sequencing of Saccharomyces cerevisiae, Homo sapiens, and Sus scrofa genomes. The datasets were used as evidence in a project that investigated the history of genomic science. To design the datasets, we first retrieved all sequence submission data from the European Nucleotide Archive (ENA), including accession numbers associated with each of our three species. Second, we used these accession numbers to construct queries to retrieve peer-reviewed scientific publications that first described these sequence submissions in the scientific literature. For each species, this resulted in two associated datasets: 1) A .csv file documenting the PMID of each article describing new sequences, all paper authors, all institutional affiliations of each author, countries of institution, year of first submission to the ENA (when available), and the year of article publication, and 2) A .csv file documenting all institutions submitting to the ENA, number of nucleotides sequenced and years of submission to the database. We utilised these datasets to understand how institutional collaboration shaped sequencing efforts, and to systematically identify important institutions and changes in the structure of research communities throughout the history of genomics and across our three target species. This data note, therefore, should aid researchers who would like to use these data for future analyses by making the methodology that underpins it transparent. Further, by detailing our methodology, researchers may be able to utilise our approach to construct similar datasets in the future.

Characteristics of Microbiota in Different Segments of the Digestive Tract of Lycodon rufozonatus

The gastrointestinal tract of animals contains microbiota, forming a complex microecosystem. Gut microbes and their metabolites can regulate the development of host innate and adaptive immune systems. Animal immune systems maintain intestinal symbiotic microbiota homeostasis. However, relatively few studies have been published on reptiles, particularly snakes, and even fewer studies on different parts of the digestive tracts of these animals. Herein, we used 16S rRNA gene sequencing to investigate the microbial community composition and adaptability in the stomach and small and large intestines of Lycodon rufozonatus. Proteobacteria, Bacteroidetes, and Firmicutes were most abundant in the stomach; Fusobacteria in the small intestine; and Proteobacteria, Bacteroidetes, Fusobacteria, and Firmicutes in the large intestine. No dominant genus could be identified in the stomach; however, dominant genera were evident in the small and large intestines. The microbial diversity index was significantly higher in the stomach than in the small and large intestines. Moreover, the influence of the microbial community structure on function was clarified through function prediction. Collectively, the gut microbes in the different segments of the digestive tract revealed the unique features of the L. rufozonatus gut microbiome. Our results provide insights into the co-evolutionary relationship between reptile gut microbiota and their hosts.

An immunoinformatics analysis: design of a multi-epitope vaccine against Cryptosporidium hominis by employing heat shock protein triggers the innate and adaptive immune responses

Cryptosporidium hominis, an anthropologically transferred species in the Cryptosporidium genus, represents many clinical studies in several countries. Its growth in the recent decade is primarily owing to epidemiologic studies. This parasite has complicated life cycles that require differentiation through a variety of phases of development and passage across two or more hosts throughout their lifetimes. As they move from host to host and environment to environment, pathogenic organisms are continually exposed to unexpected changes in the circumstances under which they develop. Heat shock proteins (HSPs) are targets of the host immune response; they are involved in the progression of diseases and play a significant part in this process. It has been discovered that the immunodominant immunogenic antigens in parasite infections HSPs. In this study, we have generated a multi-epitope vaccine against Cryptosporidium hominis (C. hominis) by using heat shock proteins. The epitopes that were selected had a substantial binding affinity for the B- and T-cell reference set of alleles, a high antigenicity score, a nature that was not allergic, a high solubility, non-toxicity and good binders. The epitopes were incorporated into a chimeric vaccine by using appropriate linkers. In order to increase the immunogenicity of the connected epitopes and effectively activate both innate and adaptive immunity, an adjuvant was attached to the epitopes. We have also analyzed the physiochemical characteristics of the vaccine which were satisfactory and then lead to the development of a 3D model. In addition, the binding confirmation of the vaccine to the TLR-4 innate immune receptor was also determined using molecular docking and molecular dynamics (MD) simulation. The results of this simulation show that the vaccine has a strong binding affinity for TLR4, which indicates that the vaccine is highly effective. In general, the vaccine that has been described here has a good potential for inducing protective and targeted immunogenicity, however, this hypothesis is contingent upon more experimental testing.Communicated by Ramaswamy H. Sarma.

psnGPCRdb: The Structure-network Database of G Protein Coupled Receptors

G protein coupled receptors (GPCRs) are critical eukaryotic signal transduction gatekeepers and represent the largest protein superfamily in the human proteome, with more than 800 members. They share seven transmembrane helices organized in an up-down bundle architecture. GPCR-mediated signaling pathways have been linked to numerous human diseases, and GPCRs are the targets of approximately 35% of all drugs currently on the market. Structure network analysis, a graph theory-based approach, represents a cutting-edge tool to deeply understand GPCR function, which strongly relies on communication between the extracellular and intracellular poles of their structure. psnGPCRdb stores the structure networks (i.e., linked nodes, hubs, communities and communication pathways) computed on all updated GPCR structures in the Protein Data Bank, in their isolated states or in complex with extracellular and/or intracellular molecules. The structure communication signatures of a sub-family or family of GPCRs as well as of their small-molecule activators or inhibitors are stored as consensus networks. The database stores also all meaningful structure network-based comparisons (i.e., difference networks) of functionally different states (i.e., inactive or active) of a given receptor sub-type, or of consensus networks representative of a receptor sub-type, type, sub-family or family. Single or consensus GPCR networks hold also information on amino acid conservation. The database allows to graphically analyze 3D structure networks together with interactive data-tables. Ligand-centric networks can be analyzed as well. psnGPCRdb is unique and represents a powerful resource to unravel GPCR function with important implications in cell signaling and drug design. psnGPCRdb is freely available at: http://webpsn.hpc.unimo.it/psngpcr.php.

Diversity of the Hydroxylamine Oxidoreductase (HAO) Gene and Its Enzyme Active Site in Agricultural Field Soils

Nitrification is a key process in the biogeochemical nitrogen cycle and a major emission source of the greenhouse gas nitrous oxide (N2O). The periplasmic enzyme hydroxylamine oxidoreductase (HAO) is involved in the oxidation of hydroxylamine to nitric oxide in the second step of nitrification, producing N2O as a byproduct. Its three-dimensional structure demonstrates that slight differences in HAO active site residues have inhibitor effects. Therefore, a more detailed understanding of the diversity of HAO active site residues in soil microorganisms is important for the development of novel nitrification inhibitors using structure-guided drug design. However, this has not yet been examined. In the present study, we investigated hao gene diversity in beta-proteobacterial ammonia-oxidizing bacteria (β-AOB) and complete ammonia-oxidizing (comammox; Nitrospira spp.) bacteria in agricultural fields using a clone library ana-lysis. A total of 1,949 hao gene sequences revealed that hao gene diversity in β-AOB and comammox bacteria was affected by the fertilizer treatment and field type, respectively. Moreover, hao sequences showed the almost complete conservation of the six HAO active site residues in both β-AOB and comammox bacteria. The diversity of nitrifying bacteria showed similarity between hao and amoA genes. The nxrB amplicon sequence revealed the dominance of Nitrospira cluster II in tea field soils. The present study is the first to reveal hao gene diversity in agricultural soils, which will accelerate the efficient screening of HAO inhibitors and evaluations of their suppressive effects on nitrification in agricultural soils.

Complex of Proline-Specific Peptidases in the Genome and Gut Transcriptomes of Tenebrionidae Insects and Their Role in Gliadin Hydrolysis

A detailed analysis of the complexes of proline-specific peptidases (PSPs) in the midgut transcriptomes of the larvae of agricultural pests Tenebrio molitor and Tribolium castaneum and in the genome of T. castaneum is presented. Analysis of the T. castaneum genome revealed 13 PSP sequences from the clans of serine and metal-dependent peptidases, of which 11 sequences were also found in the gut transcriptomes of both tenebrionid species' larvae. Studies of the localization of PSPs, evaluation of the expression level of their genes in gut transcriptomes, and prediction of the presence of signal peptides determining secretory pathways made it possible to propose a set of peptidases that can directly participate in the hydrolysis of food proteins in the larvae guts. The discovered digestive PSPs of tenebrionids in combination with the post-glutamine cleaving cysteine cathepsins of these insects effectively hydrolyzed gliadins, which are the natural food substrates of the studied pests. Based on the data obtained, a hypothetical scheme for the complete hydrolysis of immunogenic gliadin peptides by T. molitor and T. castaneum digestive peptidases was proposed. These results show promise regarding the development of a drug based on tenebrionid digestive enzymes for the enzymatic therapy of celiac disease and gluten intolerance.

Large-scale genotyping and phenotyping of a worldwide winter wheat genebank for its use in pre-breeding

Plant genetic resources (PGR) stored at genebanks are humanity's crop diversity savings for the future. Information on PGR contrasted with modern cultivars is key to select PGR parents for pre-breeding. Genotyping-by-sequencing was performed for 7,745 winter wheat PGR samples from the German Federal ex situ genebank at IPK Gatersleben and for 325 modern cultivars. Whole-genome shotgun sequencing was carried out for 446 diverse PGR samples and 322 modern cultivars and lines. In 19 field trials, 7,683 PGR and 232 elite cultivars were characterized for resistance to yellow rust - one of the major threats to wheat worldwide. Yield breeding values of 707 PGR were estimated using hybrid crosses with 36 cultivars - an approach that reduces the lack of agronomic adaptation of PGR and provides better estimates of their contribution to yield breeding. Cross-validations support the interoperability between genomic and phenotypic data. The here presented data are a stepping stone to unlock the functional variation of PGR for European pre-breeding and are the basis for future breeding and research activities.

Analyzing the postulated inhibitory effect of Manumycin A on farnesyltransferase

Manumycin A is postulated to be a specific inhibitor against the farnesyltransferase (FTase) since this effect has been shown in 1993 for yeast FTase. Since then, plenty of studies investigated Manumycin A in human cells as well as in model organisms like Caenorhabditis elegans. Some studies pointed to additional targets and pathways involved in Manumycin A effects like apoptosis. Therefore, these studies created doubt whether the main mechanism of action of Manumycin A is FTase inhibition. For some of these alternative targets half maximal inhibitory concentrations (IC₅₀) of Manumycin A are available, but not for human and C. elegans FTase. So, we aimed to 1) characterize missing C. elegans FTase kinetics, 2) elucidate the IC₅₀ and K_i values of Manumycin A on purified human and C. elegans FTase 3) investigate Manumycin A dependent expression of FTase and apoptosis genes in C. elegans. C. elegans FTase has its temperature optimum at 40°C with K_M of 1.3 µM (farnesylpyrophosphate) and 1.7 µM (protein derivate). Whilst other targets are inhibitable by Manumycin A at the nanomolar level, we found that Manumycin A inhibits cell-free FTase in micromolar concentrations (K_{i human} 4.15 μM; K_{i C. elegans} 3.16 μM). Furthermore, our gene expression results correlate with other studies indicating that thioredoxin reductase 1 is the main target of Manumycin A. According to our results, the ability of Manumycin A to inhibit the FTase at the micromolar level is rather neglectable for its cellular effects, so we postulate that the classification as a specific FTase inhibitor is no longer valid.