DNA methylation plays important roles in lifestyle transition of Arthrobotrys oligospora.

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.

MultiClassMetabo: A Superior Classification Model Constructed Using Metabolic Markers in Multiclass Metabolomics.

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

A Deep Learning Approach for Predicting Antigenic Variation of Influenza A H3N2.

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.

Mechanistic insights into SARS-CoV-2 epidemic via revealing the features of SARS-CoV-2 coding proteins and host responses upon its infection.

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.

Expansion of Adhesion Genes Drives Pathogenic Adaptation of Nematode-Trapping Fungi.

Understanding how fungi interact with other organisms has significant medical, environmental, and agricultural implications. Nematode-trapping fungi (NTF) can switch to pathogens by producing various trapping devices to capture nematodes. Here we perform comparative genomic analysis of the NTF with four representative trapping devices. Phylogenomic reconstruction of these NTF suggested an evolutionary trend of trapping device simplification in morphology. Interestingly, trapping device simplification was accompanied by expansion of gene families encoding adhesion proteins and their increasing adhesiveness on trap surfaces. Gene expression analysis revealed a consistent up-regulation of the adhesion genes during their lifestyle transition from saprophytic to nematophagous stages. Our results suggest that the expansion of adhesion genes in NTF genomes and consequential increase in trap surface adhesiveness are likely the key drivers of fungal adaptation in trapping nematodes, providing new insights into understanding mechanisms underlying infection and adaptation of pathogenic fungi.

The lifestyle transition of Arthrobotrys oligospora is mediated by microRNA-like RNAs.

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.

Protein profiling identified key chemokines that regulate the maintenance of human pluripotent stem cells.

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.

The cAMP-PKA pathway-mediated fat mobilization is required for cold tolerance in C. elegans.

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.

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations.

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.

Digital gene expression profiling analysis of DNA repair pathways in colon cancer stem population of HT29 cells.

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.

Effect of the Solvent Temperatures on Dynamics of Serine Protease Proteinase K.

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.

Insights into Protein-Ligand Interactions: Mechanisms, Models, and Methods.

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.

Molecular mechanisms of nematode-nematophagous microbe interactions: basis for biological control of plant-parasitic nematodes.

Plant-parasitic nematodes cause significant damage to a broad range of vegetables and agricultural crops throughout the world. As the natural enemies of nematodes, nematophagous microorganisms offer a promising approach to control the nematode pests. Some of these microorganisms produce traps to capture and kill the worms from the outside. Others act as internal parasites to produce toxins and virulence factors to kill the nematodes from within. Understanding the molecular basis of microbe-nematode interactions provides crucial insights for developing effective biological control agents against plant-parasitic nematodes. Here, we review recent advances in our understanding of the interactions between nematodes and nematophagous microorganisms, with a focus on the molecular mechanisms by which nematophagous microorganisms infect nematodes and on the nematode defense against pathogenic attacks. We conclude by discussing several key areas for future research and development, including potential approaches to apply our recent understandings to develop effective biocontrol strategies.

TWO INNEXINS OF Spodoptera litura INFLUENCES HEMICHANNEL AND GAP JUNCTION FUNCTIONS IN CELLULAR IMMUNE RESPONSES.

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.

Associated bacteria of different life stages of Meloidogyne incognita using pyrosequencing-based analysis.

The root knot nematode (RKN), Meloidogyne incognita, belongs to the most damaging plant pathogens worldwide, and is able to infect almost all cultivated plants, like tomato. Recent research supports the hypothesis that bacteria often associated with plant-parasitic nematodes, function as nematode parasites, symbionts, or commensal organisms etc. In this study, we explored the bacterial consortia associated with M. incognita at different developmental stages, including egg mass, adult female and second-stage juvenile using the pyrosequencing approach. The results showed that Proteobacteria, with a proportion of 71-84%, is the most abundant phylum associated with M. incognita in infected tomato roots, followed by Actinobacteria, Bacteroidetes, Firmicutes etc. Egg mass, female and second-stage juvenile of M. incognita harbored a core microbiome with minor difference in communities and diversities. Several bacteria genera identified in M. incognita are recognized cellulosic microorganisms, pathogenic bacteria, nitrogen-fixing bacteria and antagonists to M. incognita. Some genera previously identified in other plant-parasitic nematodes were also found in tomato RKNs. The potential biological control microorganisms, including the known bacterial pathogens and nematode antagonists, such as Actinomycetes and Pseudomonas, showed the largest diversity and proportion in egg mass, and dramatically decreased in second-stage juvenile and female of M. incognita. This is the first comprehensive report of bacterial flora associated with the RKN identified by pyrosequencing-based analysis. The results provide valuable information for understanding nematode-microbiota interactions and may be helpful in the development of novel nematode-control strategies.

A proposed adhesin AoMad1 helps nematode-trapping fungus Arthrobotrys oligospora recognizing host signals for life-style switching.

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.

Bacteria can mobilize nematode-trapping fungi to kill nematodes.

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.

Insight derived from molecular dynamics simulations into molecular motions, thermodynamics and kinetics of HIV-1 gp120.

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 DAF-16/FOXO transcription factor functions as a regulator of epidermal innate immunity.

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 Ca(2+) release. Ca(2+) 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-Ca(2+)-BLI-3-CST-1-DAF-16 signaling represents a previously unknown pathway to regulate epidermal damage response.