Trap formation is the key indicator of carnivorous lifestyle transition of nematode-trapping fungi (NTF). Here, the DNA methylation profile was explored during trap induction of Arthrobotrys oligospora, a typical NTF that captures nematodes by developing adhesive networks. Whole-genome bisulfite sequencing identified 871 methylation sites and 1979 differentially methylated regions (DMRs). This first-of-its-kind investigation unveiled the widespread presence of methylation systems in NTF, and suggested potential regulation of ribosomal RNAs through DNA methylation. Functional analysis indicated DNA methylation's involvement in complex gene regulations during trap induction, impacting multiple biological processes like response to stimulus, transporter activity, cell reproduction and molecular function regulator. These findings provide a glimpse into the important roles of DNA methylation in trap induction and offer new insights for understanding the molecular mechanisms driving carnivorous lifestyle transition of NTF.

Multiclass metabolomics has become a popular technique for revealing the mechanisms underlying certain physiological processes, different tumor types, or different therapeutic responses. In multiclass metabolomics, it is highly important to uncover the underlying biological information on biosamples by identifying the metabolic markers with the most associations and classifying the different sample classes. The classification problem of multiclass metabolomics is more difficult than that of the binary problem. To date, various methods exist for constructing classification models and identifying metabolic markers consisting of well-established techniques and newly emerging machine learning algorithms. However, how to construct a superior classification model using these methods remains unclear for a given multiclass metabolomic data set. Herein, MultiClassMetabo has been developed for constructing a superior classification model using metabolic markers identified in multiclass metabolomics. MultiClassMetabo can enable online services, including (a) identifying metabolic markers by marker identification methods, (b) constructing classification models by classification methods, and (c) performing a comprehensive assessment from multiple perspectives to construct a superior classification model for multiclass metabolomics. In summary, MultiClassMetabo is distinguished for its capability to construct a superior classification model using the most appropriate method through a comprehensive assessment, which makes it an important complement to other available tools in multiclass metabolomics. MultiClassMetabo can be accessed at http://idrblab.cn/multiclassmetabo/.

The LeaderRank algorithm identifies influential nodes using the idea of the random walk and performs well in the identification of hub genes in biological networks. The rapid development of biological technologies promotes the evolution of networks to be highly complex and increasing largely in size. In this study, we introduce LRCu, for computing LeaderRank using CUDA, to identify hub genes. The results show that the parallel implementation is highly effective and achieves a speedup of more than 42 times for a 13,000-gene network, thereby enhancing large-scale network computation.

Modeling antigenic variation in influenza (flu) virus A H3N2 using amino acid sequences is a promising approach for improving the prediction accuracy of immune efficacy of vaccines and increasing the efficiency of vaccine screening. Antigenic drift and antigenic jump/shift, which arise from the accumulation of mutations with small or moderate effects and from a major, abrupt change with large effects on the surface antigen hemagglutinin (HA), respectively, are two types of antigenic variation that facilitate immune evasion of flu virus A and make it challenging to predict the antigenic properties of new viral strains. Despite considerable progress in modeling antigenic variation based on the amino acid sequences, few studies focus on the deep learning framework which could be most suitable to be applied to this task. Here, we propose a novel deep learning approach that incorporates a convolutional neural network (CNN) and bidirectional long-short-term memory (BLSTM) neural network to predict antigenic variation. In this approach, CNN extracts the complex local contexts of amino acids while the BLSTM neural network captures the long-distance sequence information. When compared to the existing methods, our deep learning approach achieves the overall highest prediction performance on the validation dataset, and more encouragingly, it achieves prediction agreements of 99.20% and 96.46% for the strains in the forthcoming year and in the next two years included in an existing set of chronological amino acid sequences, respectively. These results indicate that our deep learning approach is promising to be applied to antigenic variation prediction of flu virus A H3N2.

ABSTRACT Summary There are seven known coronaviruses that infect humans: four mild coronaviruses, including HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1, only cause mild respiratory diseases, and three severe coronaviruses, including SARS-CoV, MERS-CoV and SARS-CoV-2, can cause severe respiratory diseases even death of infected patients. Both infection and death caused by SARS-CoV-2 are still rapidly increasing worldwide. In this study, we demonstrate that viral coding proteins of SARS-CoV-2 have distinct features and are most, medium and least conserved with SARS-CoV, MERS-CoV and the rest four mild coronaviruses (HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1), respectively. Moreover, expression of host responsive genes (HRG), HRG-enriched biological processes and HRG-enriched KEGG pathways upon infection of SARS-CoV-2 shows slightly overlapping with SARS-CoV and MERS-CoV but distinctive to the four mild coronaviruses. Interestingly, enrichment of overactivation of neutrophil by HRGs is only and commonly found in infections of severe SARS-CoV-2, SARS-CoV and MERS-CoV but not in the other four mild coronaviruses, and the related gene networks show different patterns. Clinical data support that overactivation of neutrophil for severe patients can be one major factor for the similar clinical symptoms observed in SARS-CoV-2 infection compared to infections of the other two severe coronavirus (SARS-CoV and MERS-CoV). Taken together, our study provides a mechanistic insight into SARS-CoV-2 epidemic via revealing the conserved and distinct features of SARS-CoV-2, raising the critical role of dysregulation of neutrophil for SARS-CoV-2 infection. Availability and Implementation All data sources and analysis methods related to this manuscript are available in the methods, supplementary materials and GEO database (https://www.ncbi.nlm.nih.gov/geo/). Supplementary information Supplementary data are available at Bioinformatics online.

Arthrobotrys oligospora, a model nematophagous fungus that produces specific adhesive networks to capture nematodes, has been proposed as a potentially effective biological agent to control harmful plant-parasitic nematodes. Although thioredoxin has been characterized as playing important roles in many cellular processes in other species, its function in nematophagous fungi has not been studied. Here, the function of a thioredoxin homolog, Aotrx1, was investigated in A. oligospora. The encoding gene of Aotrx1 in the nematophagous fungus A. oligospora was knocked out by homologous recombination; strain growth was assessed. The ΔAotrx1 strain of A. oligospora showed a significant decrease in growth rate on different media (PDA, CMY, and TG), a 70% decrease of conidia production, and a lower germination rate compared with the wild type. The mutant strain was unable to form traps to capture nematodes and was more sensitive to SDS and H2O2. Thioredoxin is involved in conidia development, trap formation, normal mycelial growth, and resistance to environmental stresses in the nematode-trapping fungus A. oligospora.

The lifestyle transition of fungi, defined as switching from taking organic material as nutrients to pathogens, is a fundamental phenomenon in nature. However, the mechanisms of such transition remain largely unknown. Here we show microRNA-like RNAs(milRNAs) play a key role in fungal lifestyle transition for the first time. We identified milRNAs by small RNA sequencing in Arthrobotrys oligospora, a known nematode-trapping fungus. Among them, 7 highly expressed milRNAs were confirmed by northern-blot analysis. Knocking out two milRNAs significantly decreased A. oligospora’s ability to switch lifestyles. We further identified that two of these milRNAs were associated with argonaute protein QDE-2 by RNA-immunoprecipitation(RIP) analysis. Three of the predicted target genes of milRNAs were found in immunoprecipitation(IP)products of QDE-2. Disruption of argonaute gene qde-2 also led to serious defects in lifestyle transition. Interestingly, knocking out individual milRNAs or qde-2 lead to diverse responses under different conditions, and qde-2 itself may be targeted by the milRNAs. Collectively, it indicates the lifestyle transition of fungi is mediated by milRNAs through RNA interference(RNAi)machinery, revealing the wide existence of miRNAs in fungi kingdom and providing new insights into understanding the adaptation of fungi from scavengers to predators and the mechanisms underlying fungal infections.

Low temperature has a great impact on animal life. Homoiotherms such as mammals increase their energy expenditure to produce heat by activating the cAMP-protein kinase A (PKA)-hormone-sensitive lipase (HSL) pathway under cold stress. Although poikilothermic animals do not have the ability to regulate body temperature, whether this pathway is required for cold tolerance remains unknown. We have now achieved this using the genetically tractable model animal Caenorhabditis elegans. We demonstrate that cold stress activates PKA signaling, which in turn up-regulates the expression of a hormone-sensitive lipase hosl-1. The lipase induces fat mobilization, leading to glycerol accumulation, thereby protecting worms against cold stress. Our findings provide an example of an evolutionarily conserved mechanism for cold tolerance that has persisted in both poikilothermic and homoeothermic animals.

Microenvironment (or niche)-providing chemokines regulate many important biological functions of tissue-specific stem cells. However, to what extent chemokines influence human pluripotent stem cells (hPSCs) is not yet completely understood. In this study, we applied protein array to screen chemokines found within the cytokine pool in the culture supernatant of hPSCs. Our results showed that chemokines were the predominant supernatant components, and came from three sources: hPSCs, feeder cells, and culture media. Chemotaxis analysis of IL-8, SDF-1α, and IP-10 suggested that chemokines function as uniform chemoattractants to mediate in vitro migration of the hPSCs. Chemokines mediate both differentiated and undifferentiated states of hPSCs. However, balanced chemokine signaling tends to enhance their stemness in vitro . These results indicate that chemokines secreted from both stem cells and feeder cells are essential to mobilize hPSCs and maintain their stemness.

Cancer stem cells (CSCs) contribute to the relapse and development of new neoplasm lesions. While most available clinical approaches, such as chemical and radiation therapies, will kill the majority of cancer cells, they do not kill them all. Some resisting cells, like CSCs, are able to survive due to their excellent self-maintaining capabilities, even in challenging environments. In the present study, we investigated the mRNA level of DNA repair genes of colon CSCs from the HT29 cell line in response to single-strand damage and double-strand breaks, as well as the evident upregulation of key genes in base excision repair, mismatch repair, non-homologous end-joining, and homologous recombination pathways in these cells. Digital gene expression analysis identified upregulated genes in CD44+ HT29 cells that may play important roles in DNA repair. Our results reveal that colon CSCs bear efficient DNA repair abilities, which might explain the survival of colon CSCs after repeated chemical and radiation therapy.

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.

Molecular recognition, which is the process of biological macromolecules interacting with each other or various small molecules with a high specificity and affinity to form a specific complex, constitutes the basis of all processes in living organisms. Proteins, an important class of biological macromolecules, realize their functions through binding to themselves or other molecules. A detailed understanding of the protein–ligand interactions is therefore central to understanding biology at the molecular level. Moreover, knowledge of the mechanisms responsible for the protein-ligand recognition and binding will also facilitate the discovery, design, and development of drugs. In the present review, first, the physicochemical mechanisms underlying protein–ligand binding, including the binding kinetics, thermodynamic concepts and relationships, and binding driving forces, are introduced and rationalized. Next, three currently existing protein-ligand binding models—the “lock-and-key”, “induced fit”, and “conformational selection”—are described and their underlying thermodynamic mechanisms are discussed. Finally, the methods available for investigating protein–ligand binding affinity, including experimental and theoretical/computational approaches, are introduced, and their advantages, disadvantages, and challenges are discussed.

To obtain detailed information about the effect of the solvent temperatures on protein dynamics, multiple long molecular dynamics (MD) simulations of serine protease proteinase K with the solute and solvent coupled to different temperatures (either 300 or 180 K) have been performed. Comparative analyses demonstrate that the internal flexibility and mobility of proteinase K are strongly dependent on the solvent temperatures but weakly on the protein temperatures. The constructed free energy landscapes (FELs) at the high solvent temperatures exhibit a more rugged surface, broader spanning range, and higher minimum free energy level than do those at the low solvent temperatures. Comparison between the dynamic hydrogen bond (HB) numbers reveals that the high solvent temperatures intensify the competitive HB interactions between water molecules and protein surface atoms, and this in turn exacerbates the competitive HB interactions between protein internal atoms, thus enhancing the conformational flexibility and facilitating the collective motions of the protein. A refined FEL model was proposed to explain the role of the solvent mobility in facilitating the cascade amplification of microscopic motions of atoms and atomic groups into the global collective motions of the protein.

The nematode-trapping fungus Arthrobotrys oligospora is an important natural enemy of nematodes. It can capture nematodes by producing a special mycelial structure called adhesive network or trap. The trap is also a signature of the fungus switching from the saprophytic lifestyle to the predacious lifestyle. At present, little is known about the mechanism of lifestyle switch in nematode-trapping fungi. Here we describe the effect of a cell wall protein called AoMad1 on lifestyle switch. The disruption of the AoMad1-encoding gene in A. oligospora resulted in the formation of more traps in the presence of nematodes. Interestingly, the mutant strain was more sensitive to certain nitrogen sources as trap inducers than the wild type strain. The microscopic examinations revealed that the AoMad1-deletion mutant lacked cell surface adhesive materials and the cell wall structures were more porous than wild-type strains. A great of genes were differentially expressed by transcriptomic analysis when trap formation was induced by sodium nitrate compared to the wild type strain, many of them were related to nitrogen metabolism, host-pathogen interaction and mycelia development. The results suggest that AoMad1 plays an important role in life style switching in A. oligospora.

Insect cellular immune responses include encapsulation, nodule formation, and phagocytosis. Hemichannels and gap junctions are involved in these cellular actions. Innexins (Inxs: analogous to the vertebrate connexins) form hemichannels and gap junctions, but the molecular mechanisms underlying their biology is still unclear. In this article, we reported a steady-state level of Inxs (SpliInxs) in hemocytes of Spodoptera litura, which formed nonfunctional hemichannels on the cell surface to maintain normal metabolism. We also reported that two innnexins (SpliInx2 and SpliInx3) were expressed significantly higher in hemocytes compared to other tissues, suggesting that they play important roles in hemocytes. Amino acid analysis found that two cysteine residues in two extracellular loops provided the capability for SpliInx2 and SpliInx3 hemichannels to dock into gap junctions. Western blotting demonstrated that both extracellular and intracellular loops of SpliInx3 and the extracellular loops of SpliInx2 might undergo posttranslational modification during the formation of a steady-state hemichannel. During hemichannel formation, SpliInx2 presented as one isoform, while SpliInx3 presented as three isoforms. These results provide fundamental knowledge for further study of how steady-state levels of SpliInxs are dynamically adjusted to perform cellular immune responses under immune challenge.

In their natural habitat, bacteria are consumed by bacterivorous nematodes; however, they are not simply passive preys. Here we report a defensive mechanism used by certain bacteria to mobilize nematode-trapping fungi to kill nematodes. These bacteria release urea, which triggers a lifestyle switch in the fungus Arthrobotrys oligospora from saprophytic to nematode–predatory form; this predacious form is characterized by formation of specialized cellular structures or ‘traps’. The bacteria significantly promote the elimination of nematodes by A. oligospora . Disruption of genes involved in urea transport and metabolism in A. oligospora abolishes the urea-induced trap formation. Furthermore, the urea metabolite ammonia functions as a signal molecule in the fungus to initiate the lifestyle switch to form trap structures. Our findings highlight the importance of multiple predator–prey interactions in prey defense mechanisms.

Although the crystal structures of the HIV-1 gp120 core bound and pre-bound by CD4 are known, the details of dynamics involved in conformational equilibrium and transition in relation to gp120 function have remained elusive. The homology models of gp120 comprising the N- and C-termini and loops V3 and V4 in the CD4-bound and CD4-unbound states were built and subjected to molecular dynamics (MD) simulations to investigate the differences in dynamic properties and molecular motions between them. The results indicate that the CD4-bound gp120 adopted a more compact and stable conformation than the unbound form during simulations. For both the unbound and bound gp120, the large concerted motions derived from essential dynamics (ED) analyses can influence the size/shape of the ligand-binding channel/cavity of gp120 and, therefore, were related to its functional properties. The differences in motion direction between certain structural components of these two forms of gp120 were related to the conformational interconversion between them. The free energy calculations based on the metadynamics simulations reveal a more rugged and complex free energy landscape (FEL) for the unbound than for the bound gp120, implying that gp120 has a richer conformational diversity in the unbound form. The estimated free energy difference of ∼-6.0 kJ/mol between the global minimum free energy states of the unbound and bound gp120 indicates that gp120 can transform spontaneously from the unbound to bound states, revealing that the bound state represents a high-probability "ground state" for gp120 and explaining why the unbound state resists crystallization. Our results provide insight into the dynamics-and-function relationship of gp120, and facilitate understandings of the thermodynamics, kinetics and conformational control mechanism of HIV-1 gp120.

The Caenorhabditis elegans DAF-16 transcription factor is critical for diverse biological processes, particularly longevity and stress resistance. Disruption of the DAF-2 signaling cascade promotes DAF-16 activation, and confers resistance to killing by pathogenic bacteria, such as Pseudomonas aeruginosa, Staphylococcus aureus, and Enterococcus faecalis. However, daf-16 mutants exhibit similar sensitivity to these bacteria as wild-type animals, suggesting that DAF-16 is not normally activated by these bacterial pathogens. In this report, we demonstrate that DAF-16 can be directly activated by fungal infection and wounding in wild-type animals, which is independent of the DAF-2 pathway. Fungal infection and wounding initiate the Gαq signaling cascade, leading to Ca2+ release. Ca2+ mediates the activation of BLI-3, a dual-oxidase, resulting in the production of reactive oxygen species (ROS). ROS then activate DAF-16 through a Ste20-like kinase-1/CST-1. Our results indicate that DAF-16 in the epidermis is required for survival after fungal infection and wounding. Thus, the EGL-30-Ca2+-BLI-3-CST-1-DAF-16 signaling represents a previously unknown pathway to regulate epidermal damage response.

Protein folding is a complex problem and involves a series of complicated processes. This problem was further divided into two main problems: (i) the thermodynamic question of how a native structur...

Sporadic canine colorectal cancers (CRCs) should make excellent models for studying the corresponding human cancers. To molecularly characterize canine CRC, we investigated exonic sequence mutations of adenomatous polyposis coli (APC), the best known tumor suppressor gene of human CRC, in 23 sporadic canine colorectal tumors, including 8 adenomas and 15 adenocarcinomas, via exon-resequencing analysis. As a comparison, we also performed the same sequencing analysis on 10 other genes, either located at human 5q22 (the same locus as APC) or 18q21 (also frequently altered in human CRC), or known to play a role in human carcinogenesis. We noted that APC was the most significantly mutated gene in both canine adenomas and adenocarcinomas among the 11 genes examined. Significantly, we detected large deletions of ≥ 10 bases, many clustered near the mutation cluster region, as well as single or two base deletions in ~70% canine tumors of both subtypes. These observations indicate that like in the human, APC is also frequently altered in sporadic colorectal tumors in the dog and its alteration is an early event in canine colorectal tumorigenesis. Our study provides further evidence demonstrating the molecular similarity in pathogenesis between sporadic human and canine CRCs. This work, along with our previous copy number abnormality study, supports that sporadic canine CRCs are valid models of human CRCs at the molecular level.